Differences
This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
test_cases:seven [2014/12/19 14:23] grenier |
test_cases:seven [2015/01/08 15:42] (current) grenier |
||
---|---|---|---|
Line 1: | Line 1: | ||
- | Advection with constant velocity is here added to conduction and phase change. | + | Heat advection via constant |
- | The paper Kurylyk et al. 2014 describes at depth the analytical solutions | + | The paper Kurylyk et al. (2014, AWR) details several alternative |
- | A presentation of the TH1 Case by Barret Kurylyk | + | A presentation of the TH1 Case was made by Barret Kurylyk during the kick off meeting ({{: |
- | The recommanded benchmark | + | Benchmark cases 2 and 3 recommended by Kurylyk et al. (2014) are included in the InterFrost project as two TH1 cases differing by the flow velocity |
+ | {{ : | ||
+ | where X is the depth to the thawing front, α is the thermal diffusivity, | ||
+ | Benchmark 1 (Neumann solution) recommended by Kurylyk et al 2014 is an option to complement the Lunardini case provided in the InterFrost project as T1. For clarification purposes, it should be noted that the initial temperature (Ti) term in the Neumann solution is expressed as the number of degrees below 0°C (i.e. it is a positive number). | ||
**TH1 Test Case** | **TH1 Test Case** | ||
- | Initial | + | The initially uniform temperature is at 0°C. This condition simplifies the energy balance at the thawing front by ensuring that there is no thermal gradient (or conductive flux) below the thawing front. The specified surface temperature is above 0°C and hence induces thaw. The water is advected through the entire medium, but the divergence of the advective flux is zero below the thawing front due to the uniform thermal |
- | {{ : | + | //Initial and boundary conditions//: |
+ | {{ : | ||
- | Parameter set: | + | The //parameter table// is accessible as {{:test_cases: |
- | {{ : | + | The // |
- | The analytical solutions of Kurylyk et al 2014 are contained in this {{: | + | Code developers |
The approach used for SUTRA code is summed up in two figures from Kurylyk et al 2014 as a source of inspiration. They provide the simulated domain and the approach for the freezing curve function in order to approximate the step function with a linear curve: | The approach used for SUTRA code is summed up in two figures from Kurylyk et al 2014 as a source of inspiration. They provide the simulated domain and the approach for the freezing curve function in order to approximate the step function with a linear curve: | ||
- | {{ :test_cases:figkurylyk.jpg?300 |}} {{ :test_cases:freezingcurvesutra.jpg? | + | {{test_cases: |
+ | {{test_cases: | ||