.Buildings.Fluid.Geothermal.Aquifer.MultiWell .Buildings.Fluid.Geothermal.Aquifer.MultiWellThis model simulates aquifer thermal energy storage, using one or multiple pairs of cold and hot wells.
To calculate aquifer temperature at different locations over time, the model applies physical principles of water flow and heat transfer phenomena. The model is based on the partial differential equation (PDE) for 1D conductive-convective transient radial heat transport in porous media
ρ c (∂ T(r,t) ⁄ ∂t) = k (∂² T(r,t) ⁄ ∂r²) - ρw cw u(∂ T(r,t) ⁄ ∂t),
where ρ is the mass density, c is the specific heat capacity per unit mass, T is the temperature at location r and time t, u is water velocity and k is the heat conductivity. The subscript w indicates water. The first term on the right hand side of the equation describes the effect of conduction, while the second term describes the fluid flow.
The pressure losses in the aquifer are calculated using the Darcy's law
Δp = ṁ g ⁄ (2 π K h ln(rMax ⁄ rWB)),
where ṁ is the water mass flow rate, g is the gravitational acceleration, K is the hydraulic conductivity, h is the thickness of the aquifer, rMax is the domain radius and rWB is the well radius. The pressure losses in the wells are calculated using Modelica.Fluid.Pipes.BaseClasses.WallFriction.Detailed.pressureLoss_m_flow.
To discretize the conductive-convective equation, the domain is divided into a series of thermal capacitances and thermal resistances along the radial direction. The implementation uses an array of Modelica.Thermal.HeatTransfer.Components.HeatCapacitor and Modelica.Thermal.HeatTransfer.Components.ThermalResistor. Fluid flow was modelled by adding a series of fluid volumes, which are connected to the thermal capacitances via heat ports. The fluid stream was developed using the model Buildings.Fluid.MixingVolumes.MixingVolume. The geometric representation of the model is illustrated in the figure below.
By default, the component consists of a single pair of wells:
one cold well and one warm well. The number of paired wells can be
increased by modifing the parameters nPai. The effect
is a proportional increase of thermal capacity, and the mass flow
rate at port_a and port_b is equally
distributed to the pairs of well, thus all pairs have the same mass
flow rates and temperatures, and the quantities at the fluid ports
is for all well combined, as is the electricity consumption
PTot for the well pumps.
To ensure conservation of energy, the wells are connected via
fluid ports. To avoid thermal interferences, make sure that the
aquifer domain radius rMax is large enough for your
specific use case.
Circulation pumps are included in the model and they can be
controlled by acting on the input connector. The input must vary
between [1, -1]. A positive value will circulate water
clockwise (from port_Hot to port_Col,
thus extraction from the cold well and injection into the warm
well). A negative value will circulate water anticlockwise.
The temperature values in the warm and cold aquifers can be
accessed using TAquHot and TAquCol. These
temperatures correspond to the temperatures of each thermal
capacitance in the discretized domain. The location of the thermal
capacitance is expressed by rVol.
The model computes the flow resistance of the cold well and the cold acquifer. Because of symmetry, the warm side has the same flow resistances. To reduce the number of equations in the model, the flow resistance of the warm side is set to the same flow resistance as on the cold side, but with the sign reversed.
| Name | Description |
|---|---|
 Medium |
Medium in the component |