Using a geothermal gradient with Richards mode

[Orlando Silva Rojas]

Hi PFLOTRAN group

Is there a way that the Richards flow mode considers a geothermal gradient when coupled to GIRT? I'm using the initial and boundary conditions shown below, the input file runs without crashing but the temperature is kept uniform and constant at 25 ºC throughout the domain.

FLOW_CONDITION initial
  TYPE
    LIQUID_PRESSURE hydrostatic
  /
  GRADIENT
    TEMPERATURE 0.d0 0.d0 -0.0141d0
  /
  LIQUID_PRESSURE 636608.0
  TEMPERATURE 6.d0
  DATUM 0.d0 0.d0 0.0d0
END

FLOW_CONDITION surface
  TYPE
    LIQUID_PRESSURE Dirichlet
    TEMPERATURE Dirichlet
  /
  LIQUID_PRESSURE DATASET hydrostaticpressure
  TEMPERATURE 6.d0
/

Many thanks in advance

Orlando

[Hammond, Glenn E]

Orlando,

 

The requested capability does not exist at this point in time. RICHARDS assumes a uniform constant temperature, though the temperature can be set to another value using the REFERENCE_TEMPERATURE card.

 

Glenn

 

From: pflotra...@googlegroups.com <pflotra...@googlegroups.com> On Behalf Of Orlando Silva Rojas
Sent: Thursday, May 2, 2024 8:48 AM
To: pflotra...@googlegroups.com
Subject: [pflotran-users: 7705] Using a geothermal gradient with Richards mode

 

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 on the web visit https://groups.google.com/d/msgid/pflotran-users/CACSNO9Pi9cG-2WckJi53aibG48_KXMTmXAtd1nUDfaFTm6iJxQ%40mail.gmail.com.

[Orlando Silva Rojas]

Thank you for the response, Glenn. Just to completely understand this, I assume the REFERENCE_TEMPERATURE may be considered changing in time (at least, by running a series of input files with different REFERENCE_TEMPERATURE). Is that correct?

All the best

Orlando

To view this discussion on the web visit https://groups.google.com/d/msgid/pflotran-users/PH0PR09MB7836496CC18A66E3664C67B09A1F2%40PH0PR09MB7836.namprd09.prod.outlook.com.

[Hammond, Glenn E]

I believe that will work. Temperature is not checkpointed for RICHARDS. 

Glenn

On May 4, 2024, at 12:56 AM, Orlando Silva Rojas <orlando.s...@gmail.com> wrote:



To view this discussion on the web visit https://groups.google.com/d/msgid/pflotran-users/CACSNO9MzeLD9NQi49toSBROnUtFkFyT5f5xoyTQ%2BLd0vwFxF6w%40mail.gmail.com.

[Yiwei Cai]

Dear all,

I am trying to changing the temperature to 300 C, and I add the reference_temperature 300. to the input file. However, it was found that the output temperature did not change. The system temperature is still at the default 25. I don't know if something went wrong.

REFERENCE_TEMPERATURE 300.d0

SIMULATION                               #
  SIMULATION_TYPE SUBSURFACE             # set the simulation type
  PROCESS_MODELS                         # use the defined model
    SUBSURFACE_FLOW flow                 # flow process model
      MODE RICHARDS                      # single phase variably-saturation flow
    /
    SUBSURFACE_TRANSPORT transport       # use reaction-transport model
      MODE GIRT                          # golbal implicit reactive transport
    /
  /
END

SUBSURFACE

Best,

Yiwei

[Hammond, Glenn E]

Yiwei,

 

Was REFERENCE_TEMPERATURE placed within the SUBSURFACE block or at the top of the file as shown below? It needs to be placed within the SUBSURFACE block.

 

Glenn

 

To view this discussion on the web visit https://groups.google.com/d/msgid/pflotran-users/d046e3ef-d99b-4f4d-bf5c-32f5bb71a8ben%40googlegroups.com.

[Yiwei Cai]

Dear Glenn,

I tried to put the REFERENCE_TEMPERATURE within the SUBSURFACE, but the PFLOTRAN will complain: Keyword "REFERENCE_TEMPERATURE" not recognized in Richards Options.

SIMULATION                               #
  SIMULATION_TYPE SUBSURFACE             # set the simulation type
  PROCESS_MODELS                         # use the defined model
    SUBSURFACE_FLOW flow                 # flow process model
      MODE RICHARDS                      # single phase variably-saturation flow

        OPTIONS
          REFERENCE_TEMPERATURE 300.d0
        /
      /

    /
    SUBSURFACE_TRANSPORT transport       # use reaction-transport model
      MODE GIRT                          # golbal implicit reactive transport

#      OPTIONS
#       ISOTHERMAL
#      /
    /
  /
END

SUBSURFACE

Best,

Yiwei