Welcome,
Guest
. Please
login
or
register
.
Did you miss your
activation email
?
1 Hour
1 Day
1 Week
1 Month
Forever
Login with username, password and session length
Forum Home
Login
Register
PhreeqcUsers Discussion Forum
»
Processes
»
Reactive transport modelling
»
1-D heat transport in two different matrix
« previous
next »
Print
Pages: [
1
]
Go Down
Author
Topic: 1-D heat transport in two different matrix (Read 80 times)
Jeonghwan Hwang
Top Contributor
Posts: 37
1-D heat transport in two different matrix
«
on:
August 30, 2021, 08:23:47 AM »
Hello, this is Jeonghwan Hwang.
I want to make the
1-D transport for heat and solute at two different matrix, bentonite and granite
.
I made a total of 200 cells divided with 50 cells of bentonite and 150 cells of granite.
Then, I set up the 120 cel at inlet boundary and 35 cel at ourlet boundary.
Finally, the flow_direction was set up to "diffusion only" because the bentonite and granite were assumed to have very low hydraulic conductivity.
At this part, I have a question. The bentonite and granite have different thermal diffusion coefficients.
However, at the "TRANSPORT", I can define only one value,
Example code which I set up was written as below;
The code was made to apply for IPHREEQC at matlab.
----------------------------------------------------------------------------------------------
Transport = {' TRANSPORT', ...
' -boundary_conditions constant constant', ...
[' -cells ', num2str(Cell_No)], ...
[' -lengths ', num2str((Bentonite_thickness+Granite_thickness)/Cell_No)], ...
' -shifts 1', ...
[' -time_step ',num2str(Time_Step), ' year'], ...
' -flow_direction diffusion_only', ...
[' -thermal_diffusion ', num2str(Bentonite_Thermal_Retardation), ' ', num2str(Bentonite_Thermal_Diffusion_Coefficient)],
...
' -diffusion_coefficient 0.3e-9', ...
' -dispersivities 0.0', ...
' END'};
----------------------------------------------------------------------------------------------
Can I set the different two diffusion coefficients at one transport?
Thank you for reading.
Sincerely,
Jeonghwan Hwang.
Logged
dlparkhurst
Top Contributor
Posts: 2409
Re: 1-D heat transport in two different matrix
«
Reply #1 on:
August 31, 2021, 04:14:34 PM »
As far as I know, you only get one heat retardation factor.
My only thought is that if you can calculate the temperature distribution with time from another model or an analytic equation, you can use IPhreeqc to set the temperature in each cell with REACTION_TEMPERATURE n data blocks, or modify the temperature of the solutions with SOLUTION_MODIFY data blocks.
Logged
dlparkhurst
Top Contributor
Posts: 2409
Re: 1-D heat transport in two different matrix
«
Reply #2 on:
September 05, 2021, 09:09:14 PM »
Tony Appelo gives the following suggestion:
I'd advise to do transport in matlab using the formulas from Chpt. 11 in Appelo
and Postma, and to use iphreeqc to obtain the chemical reactions and concentrations.
Logged
Print
Pages: [
1
]
Go Up
« previous
next »
PhreeqcUsers Discussion Forum
»
Processes
»
Reactive transport modelling
»
1-D heat transport in two different matrix