Is it possible to put the temperature of 'Solution 0' as a function over time?

[Jeonghwan Hwang]

Hello, this is Jeonghwan Hwang.
I am studying about the thermal fluid transports in high-level waste disposal repository using 1-D transport model.

To make the accurate model, the 'solution 0' as a boundary solution should has a diverse temperature as increasing time.
Is it possible to put the temperature of 'Solution 0' as a function over time?
My model based on Exercise 12 in PHREEQC, and codes attached as below;

Thank you for your time
Sincerely,

Jeonghwan Hwang


============================================================
TITLE Example 12.--Advective and diffusive transport of heat and solutes.

PRINT       
-reset false

SELECTED_OUTPUT
   -file    Exe.txt
   -high_precision  true
   -reset   false
   -step true
   -solution true
   -time ture
   -distance    true
   -temp    true
   -pH   true
   -percent_error true   
   
USER_PUNCH
   -head Na_mmol K_mmol Cl_mmol
   10 PUNCH TOT("Na")*1000, TOT("K")*1000, TOT("Cl")*1000
END

PRINT       
-SELECTED_OUTPUT false
SOLUTION 1-25   
    units mol/kgw
    temp 25   
    pH   7.0
    pe  12.0   O2(g) -0.67
    K    1e-6   
    Cl   1e-6

SOLUTION 0   
    units mol/kgw
    temp  0                 
    pH    7.0
    pe   12.0   O2(g) -0.67
    K    24.e-3      
    N(5) 24.e-3   
END
   
PRINT
-SELECTED_OUTPUT true   
TRANSPORT                   
   -boundary_conditions   constant  constant
   -cells   25
   -lengths  0.02
   -shifts  25   
   -time_step 1.0e+1 year     
   -flow_direction  diffusion_only
   -thermal_diffusion  1.0
   -diffusion_coefficient  0.3e-9
   -dispersivities    0.0
END

[dlparkhurst]

No. All you can do is run TRANSPORT for a number of shifts, redefine the temperature of SOLUTION 0, and run TRANSPORT for another period.

The simplest way to change the temperature is SOLUTION_MODIFY. To run another transport period, you need only define -shifts and the other parameters that are different from the previous TRANSPORT calculation.

[Jeonghwan Hwang]

Thank you for your answer.

I decided to use the matlab based IPHREEQC to implement the temperature over time as a function.
The model will be implemented by calculating the temperature at each time intervals and performing transport for that time.
I hope to get meaningful results.

Thanks again.
Sincerely,

Jeonghwan Hwang

[djacques]

Dear,

for 1D solute transport coupled with heat transport (and phreeqc), HPx is an option as well (free - on https://www.pc-progress.com/en/Default.aspx?hydrus-1d )