HydroGeoSphere/Pore Water Freezing and Thawing
To calculate the freezing and thawing of pore water, HydroGeoSphere requires the following set of instructions:
Pore water freezing-thawing...End
editgrok reads instructions that define the pore water freezing and thawing parameters for the domain until it encounters an End instruction.
- • • •
- • • •
For example:
pore water freezing-thawing
surface temperature
0 -7.0
2592000 -6.1
5184000 -0.9
7776000 6.2
10368000 12.9
12960000 18.0
15552000 20.2
18144000 19.3
20736000 14.8
23328000 8.6
25920000 2.4
28512000 -4.0
31104000 -7.0
end
thermal diffusivity
0.2d-6
background temperature
6.0d0
melting temperature
0.5d0
freezing temperature
-0.5d0
maximum freezing depth
2.0d0
integration convergence criteria
1.0e-4
memory length for convolution integral
15552000.0d0 ! 6*30*86400 sec = 6 months
end
Surface temperature
edit- time(i), val(i)...end The time and surface temperature value list.
At time time(i) the value val(i) is applied and maintained until time time(i+1). The last value entered in the list will be applied until the end of the simulation.
- • • •
- • • •
Thermal diffusivity
edit- val Thermal diffusivity [L2 T−1].
- • • •
- • • •
Background temperature
edit- val Background temperature [Θ].
- • • •
- • • •
Melting temperature
edit- val Melting temperature [Θ].
- • • •
- • • •
Freezing temperature
edit- val Freezing temperature [Θ].
- • • •
- • • •
Maximum freezing depth
edit- val Maximum freezing depth [L].
- • • •
- • • •
Integration convergence criteria
edit- val Integration convergence criteria.
- • • •
- • • •
Memory length for convolution integral
edit- val Memory length for convolution integral.
- • • •
- • • •