modelATES_confinedLayer
Information
1. Purpose of model
Top-level assembly model of a confined high-temperature aquifer thermal energy storage (HT-ATES) system. Couples the aquifer (AquiferL4) with upper and lower impermeable confining layers (ImpermeableLayer arrays) and the injection/extraction well (Well_L4). Provides a single fluid inlet port for connection to the district heating network.
2. Level of detail, physical effects considered, and physical insight
L4: Full 2D (radial-vertical) discretization of aquifer and aquitards. Boundary conditions at the outer radial and top/bottom thermal boundaries are represented by fixed temperature sources. Hydrostatic pressure gradient is applied at the outer radial fluid boundary. Physics toggles (buoyancy, dispersion) are propagated from the Setting record.
3. Limits of validity
- Outer radial boundary: fixed pressure (hydrostatic) and fixed temperature.
- Top and bottom thermal boundaries of the aquitard: fixed temperature.
- No capillary effects or multi-phase flow.
4. Interfaces
fluidPortIn: fluid inlet/outlet connection to external network (Fluid, bidirectional via well)
5. Nomenclature
N_r |
number of radial control volumes [-] | |
N_z |
number of vertical aquifer control volumes [-] | |
fluidPortIn.m_flow |
injection/extraction mass flow rate [kg/s] | |
T_Boundary |
far-field boundary temperature [K] | |
p_Boundary |
reference pressure at aquifer top [Pa] |
6. Governing Equations
(no equations — assembly model; governing equations are implemented in sub-models AquiferL4, ImpermeableLayer, Well_L4)
7. Remarks for Usage
Pass site-specific material properties via the Parameters record (e.g., Molz1983) and grid specification via the Setting record. An inner ClaRa.SimCenter must be present in the enclosing model.
8. Validation
Tested in check model "ATES.Validierung.Buoyancy.Val_homogenLayer"
9. References
(no remarks)
10. Version History
Model created by Markus Gillner (markus.gillner@tuhh.de), May 2026
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| Integer | N_r | geo.N_r | |
| Integer | N_z | geo.N_z | |
| Integer | N_z_c | geo.N_z_c | |
| TILMedia.VLEFluidTypes.BaseVLEFluid | water | simCenter.fluid1 | Heat carrier |
| SI.Pressure | p_Boundary | Parameters.p_initial | Pressure of the fluid boundary in bar |
| SI.Temperature | T_Boundary | Parameters.T_initial | Temperature of the temperatur boundaries in degree celsius |
| SI.Length[N_r,N_z_c] | r_new | {{geo.r[j] for i in 1:N_z_c} for j in 1:N_r} | |
| SI.Length[N_r,N_z] | r_newA | {{geo.r[j] for i in 1:N_z} for j in 1:N_r} | |
| SI.Area[N_r,N_z_c] | A_q_new | {{geo.A_q[j] for i in 1:N_z_c} for j in 1:N_r} | |
| SI.Area[N_r,N_z] | A_q_newA | {{geo.A_q[j] for i in 1:N_z} for j in 1:N_r} | |
| Boolean | with_buoyancy | setting.buoyancy | logical variable for switching buoyancy physics on (true) or off (false) |
| Real | f_buoyancy | if with_buoyancy then 1 else 0 | factor for switching buoyancy physics on (1) or off (0) |
| SI.Length[N_z] | H_cell | {geo.z_A[1, j] for j in 1:N_z} |
Connectors
| Type | Name | Default | Description |
|---|---|---|---|
| TransiEnt.Basics.Interfaces.Thermal.FluidPortIn | fluidPortIn |