Residual vectors and Jacobian matrices for coupled systems

[Bex Abylkhani]

Hi All! 

I have a question on printing jacobian and residual vectors. My simulation uses Richard's mode where coupled flow and transport is being solved sequentially. Currently, Jacobian and residual vectors are only being printed for the flow equation. Is there any way to obtain Jacobian matrix and residual vectors for transport equation as well?

________

DEBUG
  PRINT_SOLUTION
  PRINT_RESIDUAL
  PRINT_JACOBIAN
  FORMAT MATLAB
  APPEND_COUNTS_TO_FILENAMES
/

________

Many thanks! 

Sincerely, 

Bex. 

[Hammond, Glenn E]

Assuming that you have reactive transport enabled, you should see the residual and Jacobian output as:

 

RTxx_*.mat

RTresidual_*.mat

RTjacobian_*.mat

RTjacobianLog_*.mat

 

Glenn

 

From: pflotra...@googlegroups.com <pflotra...@googlegroups.com> on behalf of Bex Abylkhani <bexu...@stanford.edu>
Date: Wednesday, May 7, 2025 at 7:42 AM
To: pflotran-users <pflotra...@googlegroups.com>
Subject: [pflotran-users: 8385] Residual vectors and Jacobian matrices for coupled systems

Check twice before you click! This email originated from outside PNNL.

 

--
You received this message because you are subscribed to the Google Groups "pflotran-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pflotran-user...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/pflotran-users/0bf65b01-4c28-44ff-ab30-9fa586f89258n%40googlegroups.com.

[Bex Abylkhani]

Hi Glenn! 

Thank you for your response! 

I do not have reactive transport but "NUMERICAL_METHODS TRANSPORT" enabled. Currently, I have only one jacobian and residual being printed for transport (without ts, tc, or ni information). The screenshot is attached:  

 

Sincerely, 
Bex. 

[Hammond, Glenn E]

Bex,

 

Please give the attached input file that includes the correct debugging block a try:

 

DEBUG

  PRINT_SOLUTION

  PRINT_RESIDUAL

  PRINT_JACOBIAN

  FORMAT MATLAB

  APPEND_COUNTS_TO_FILENAMES

/ 

 

You must have the reactive transport process model enabled:

 

SIMULATION

  SIMULATION_TYPE SUBSURFACE

  PROCESS_MODELS

    SUBSURFACE_FLOW flow

      MODE RICHARDS

    /

    SUBSURFACE_TRANSPORT transport

      MODE GIRT

    /

  /

END

 

Glenn

 

From: pflotra...@googlegroups.com <pflotra...@googlegroups.com> on behalf of Bex Abylkhani <bexu...@stanford.edu>
Date: Tuesday, May 13, 2025 at 9:31 AM
To: pflotran-users <pflotra...@googlegroups.com>
Subject: Re: [pflotran-users: 8393] Residual vectors and Jacobian matrices for coupled systems

Hi Glenn! 

 

Thank you for your response! 

 

I do not have reactive transport but "NUMERICAL_METHODS TRANSPORT" enabled. Currently, I have only one jacobian and residual being printed for transport (without ts, tc, or ni information). The screenshot is attached:  

 

To view this discussion visit https://groups.google.com/d/msgid/pflotran-users/e68b2aec-34de-4b31-8672-92df8e6d7c1fn%40googlegroups.com.