Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient storage model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

Online Optimization of NSLS-II Dynamic Aperture and Injection Transient (Facility Improvement Project)

The Facility Improvement Project (FIP) “Methods of online optimization of NSLS-II storage ring concurrent with user operations” proposed and approved in 2018 is now complete. The first goal is the online optimization of nonlinear beam dynamics to increase the beam lifetime and injection efficiency. We have developed a model-independent optimization technique using advanced algorithms. Using this technique, we increased the NSLS-II dynamic aperture by 20% and reduced the amplitude-dependent tune shift by a factor of two. We applied sextupole optimization to the new high-chromaticity lattice, which has been developed to improve the beam stability and to increase the single-bunch beam intensity. We were able to double the injection efficiency and exceed 90%. The second goal of this project is to provide the top-off injection with minimized perturbations of the beamline user operations. To achieve the transparent injection, we optimize the matching of four storage ring injection kickers reducing the perturbation of the stored beam orbit. The first application of the online optimization resulted in a reduction of the amplitude of residual beam oscillation by a factor of six, from 1300 μm down to 200 μm. Further improvement was limited by the timing jitter of the trigger boards. Then, we replaced all trigger boards with a new design and reduced the timing jitter by a factor of five, from 5 ns down to 1 ns. This upgrade resulted in the injection transient reduction from 200 μm to 120 μm. To optimize the full set of kicker parameters, including the trigger timing, amplitude, and pulse width, we upgraded all kicker power supplies with the capability of tunable waveform width. As a result, we have reduced the injection transient to the limit of 60 μm.

43 PARTICLE ACCELERATORS↗

New insights on carbon black suspension rheology—Anisotropic thixotropy and antithixotropy

We report a detailed experimental study of peculiar thixotropic dynamics of carbon black (CB, Vulcan XC-72) suspensions in mineral oil, specifically the observation of sequential stress increase then decrease at a fixed shear rate in a step-down test. We verify that such dynamics, though peculiar, come from a true material response rather than experimental artifacts. We also reveal how this long-time stress decay is associated with antithixotropic structural change rather than viscoelastic stress relaxation by using orthogonal superposition (OSP) rheometry to probe viscoelastic moduli during the step-down tests. The orthogonal storage and loss moduli are present, showing this two-timescale recovery then decay response, which demonstrates that this response is antithixotropic, and it involves shear-induced structuring. We further show a mechanical anisotropy in the CB suspension under shear using OSP. Based on the rheological results, a microstructural schematic is proposed, considering qualitatively thixotropic structure build-up, antithixotropic densification, and anisotropic structure evolution. Our observation for these CB suspensions is outside the standard paradigm of thixotropic structure-parameter models, and the elastic response provides us with new insight into the transient dynamics of CB suspensions.

Mechanics↗

Targeted modulation of pennycress lipid droplet proteins impacts droplet morphology and seed oil content

SUMMARY Lipid droplets (LDs) are unusual organelles that have a phospholipid monolayer surface and a hydrophobic matrix. In oilseeds, this matrix is nearly always composed of triacylglycerols (TGs) for efficient storage of carbon and energy. Various proteins play a role in their assembly, stability and turnover, and even though the major structural oleosin proteins in seed LDs have been known for decades, the factors influencing LD formation and dynamics are still being uncovered mostly in the “model oilseed” Arabidopsis. Here we identified several key LD biogenesis proteins in the seeds of pennycress, a potential biofuel crop, that were correlated previously with seed oil content and characterized here for their participation in LD formation in transient expression assays and stable transgenics. One pennycress protein, the lipid droplet associated protein‐interacting protein (LDIP), was able to functionally complement the Arabidopsis ldip mutant, emphasizing the close conservation of lipid storage among these two Brassicas. Moreover, loss‐of‐function ldip mutants in pennycress exhibited increased seed oil content without compromising plant growth, raising the possibility that LDIP or other LD biogenesis factors may be suitable targets for improving yields in oilseed crops more broadly.

Guzha, Athanas↗

Thermal Modeling and Simulation of the Packed-bed Thermal Energy Storage Combined with INL Thermal Energy Distribution System

Dynamic Energy Transport and Integration Lab (DETAIL) at Idaho National Laboratory is to support experimental demonstration and validation research on Nuclear-Renewable Hybrid Energy System [1]. The Thermal Energy Distribution System (TEDS) is a thermal-hydraulic flow loop in DETAIL with its own dedicated control system to support the integration of co-located multiple experimental systems, where a packed-bed thermal energy storage (TES) is installed as a thermal buffer and storage unit. Among various TES options, the packed-bed TES is adopted in TEDS because it offers a low-cost single-tank thermal storage option compared to the traditional two-tank TES. However, since the TES tank is filled with granular fillers having different thermophysical properties from those of TES tank wall, there is a potential thermo-mechanical issue to be carefully addressed like thermal ratcheting which may pose a significant design concern for the packed-bed TES tank. Thermal ratcheting is a thermomechanical process caused by the repeated rearrangement of granular filler inside a TES tank during continuous thermal cycling operation of the packed-bed TES system. If the thermally induced stress exceeds the yield strength of a TES tank wall, catastrophic consequences may happen like rupture of the TES tank. Thus, it is important to understand the phenomenon to ensure the robust operation of the packed-bed TES tank. Given that thermal ratcheting is caused by complex interaction of thermal transport in the porous bed and solid mechanics, the accurate prediction of transient thermal behavior of the packed-bed TES tank, which is the focus of this paper, is critical to the reliable thermal ratcheting analysis. This paper discusses the numerical modeling, simulation, and validation studies that are ongoing at INL to investigate the transient thermal behavior of the packed-bed TES. Of particular concern is the transient thermal process occurring in the packed-bed TES unit that is operated in conjunction with the INL TEDS. The main goal of this research is three-fold: (i) provide preliminary insights into the transient thermal behavior of the TEDS TES tank, (ii) support the thermal measurement and validation plan for TEDS experiment, and (iii) provide transient thermal boundary conditions to support the reliable thermal ratcheting analysis of the TEDS TES tank. For the transient thermal modeling and analysis, a CFD model was developed, and the validity of the modeling approach was examined via comparing the numerical simulation results with the experimental data obtained from various design characteristics of packed-bed TES tanks. Then, the present modeling method was applied for the transient thermal analysis of the TEDS TES tank, and the results are discussed along with the potential improvement of data acquisition strategy for the future TEDS experiments for more precise validation study.

25 ENERGY STORAGE↗

Wind Power as a Virtual Synchronous Generator (WindVSG)

This project investigated the theory, implemented it in hardware, and validated the Wind as a Virtual Isochronous Generator (WindVSG) concept by combining the advantages of modern dynamic inverter technologies with static, dynamic, and transient electromechanical properties of synchronous machines. During this project we demonstrated how to control the inverters of wind turbine generators (wind alone or in parallel with other GFM sources, such battery energy storage) so that wind power behaves like a synchronous machine-based power plant with a conventional prime mover. For this purpose, testing was conducted at NLR ARIES facility with real 2.5 MW wind-turbine generator operating in GFM mode under dynamic and transient conditions. The team also developed models and conducted simulations for GFM wind power to evaluate stability impacts of GFM operation on power grid. This report describes efforts by the NLR team working in collaboration GE Vernova during 3-year project.

17 WIND ENERGY↗

Design and optimization of flexible decoupled high-temperature gas-cooled reactor plants with thermal energy storage

Advanced nuclear power plants are well-positioned for future zero-carbon grids, however, the need for flexible power generation will be required over the traditional emphasis on baseload generation for meeting historical demands. To achieve such flexibility, this work examines viable configurations for coupling nuclear energy production with thermal energy storage. Previous designs on nuclear-thermal energy storage configurations for advanced reactor designs, which utilized reactor steam as the heat source for charging the thermal energy storage, are restricted by the heat diversion ratio and efficiency losses, thus their impacts can be limited. In this context, this study proposes configurations for fully decoupling the nuclear reactor from the power cycle and positioning the storage as an intermediate loop, thereby achieving an unconstrained heat diversion ratio and improved efficiency. Compared with a standard high-temperature gas-cooled reactor’s power cycle, steady-state thermodynamic modeling and dispatch optimizations quantify the benefits of a steam reheat cycle within the fully decoupled thermal energy system to separate the plant cycle from the high-pressure primary side. These benefits are further detailed, compatible with required high-temperature and high-pressure conditions, through (1) open-source dynamic transient models that examine the impact of off-design operation on the systems, (2) the investigation of components design and costing and finally (3) sizing and dispatch optimization. The fully-decoupled design achieves a cycle efficiency of 43.1%, an enhancement over the vendor’s standard efficiency of 42.2% (Xe-100 design). Here, the proposed design offers strengthened physical barriers from the nuclear island as well as superior operational flexibility and power boosting. Dispatch optimization and market analysis reveal that thermal energy storage size is highly dependent on the peak patterns of electricity prices and the minimum generation level constraint imposed on the balance of plant. Evaluation of off-design operation demonstrates that the full decoupling design with the suggested fail-safe control mechanisms ensures a minimal impact on reactor parameters, even during rapid power ramping.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Tundra SGS – Critical Importance of Site-Specific Injectivity Testing for Carbon Dioxide Injection Well Simulation Calibration

Paper presented at 16th International Conference on Greenhouse Gas Control Technologies (GHGT-16), October 23–27, 2022, Lyon, France. Brine injectivity tests were conducted on both the Broom Creek (BC) and Deadwood (DW) storage formations in one of the cored stratigraphic test wells prior to finalized simulations. The injectivity tests entailed injection of brine into the formations, followed by shutting in the wells. The shut-in duration represents the fall-off period which was used to carry out pressure transient analyses of the tests. The analysis of injectivity tests indicated that core-indicated and modeled permeability distributions should be conservatively multiplied by 2.5× and 5× in the BC and DW geologic models, respectively, prior to numerical simulations of CO 2 injection using Computer Modelling Group’s (CMG’s) GEM equation-of-state reservoir simulator. The scaler permeability multipliers allowed the fluid flow simulations to fully represent well test-measured injection well rates in the respective formations.

20 FOSSIL-FUELED POWER PLANTS↗

Enhancing EnergyPlus capabilities to model dynamic building envelopes using python plugin

Nearly half of the energy consumption in the United States is related to buildings, resulting in an urgent need to develop innovative technologies to improve building energy efficiency. Dynamic building envelopes, comprising switchable insulation and thermal energy storage materials, have been proposed recently as a promising solution to reduce buildings' heating and cooling loads by thermally coupling the indoor environment with the ambient environment when beneficial while decoupling them when outdoor conditions are not favorable. Although various related technologies are still underway, the whole-building energy modeling tools, like EnergyPlus, do not have the capability to simulate the transient and dynamic nature of dynamic envelope materials and components to accurately predict their impact on building energy use. The objective of this study is to formulate a method in EnergyPlus simulation engine to model multilayer envelopes, comprising dynamic building materials with variable thermophysical properties, and discuss the changes made to the program using a Python plugin. Furthermore, the thermal performance of the dynamic envelopes using the proposed method is compared and verified with the results from a well-established commercial code, COMSOL Multiphysics. A parametric assessment is also conducted to evaluate the energy efficiency benefits of dynamic envelopes in a single-family residential building, demonstrating total annual energy savings up to 11.6 %, when a dynamic envelope operates alone, and up to 18.2 % when it is combined with a thin layer of phase change material as a thermal storage medium. Finally, a United States wide energy efficiency assessment is presented to showcase the geographical spread of the energy savings. The method designed and implemented in this study provides the researchers with the ability to implement their dynamic insulation methods in EnergyPlus and evaluate the whole building energy impact.

25 ENERGY STORAGE↗

Validation of the Fault Ride-Through Response of a Generic EMT Inverter Model by Laboratory Testing

This paper presents the validation of the fault ride-through response of a large-scale inverter model using laboratory testing of a commercial inverter. First, we present the generic inverter model specifications. The generic inverter control is developed based on the fault ride-through response performance requirements of IEEE Std 2800-2022 for inverter-based resources. Next, we present laboratory tests that emulate voltages at the inverter terminal resulting from balanced and unbalanced faults on the transmission system. The tests are performed on a commercia12.2 MVA battery storage inverter at the U.S. National Renewable Energy Laboratory. Comparisons of the model's response to the laboratory measurements show the generic model presented can predict the controlled response of the commercial equipment. The differences between the model's response and the laboratory measurements during the initial transient at fault inception are discussed.

fault ride-through↗

Thermal Conductivity of a Novel Salt-Hydrogel Complex as a Phase Change Material for Building Thermal Management

Use of electricity for space heating, ventilation and air conditioning in buildings may be better managed through temporary energy storage. Compared with electrochemical storage, the levelized cost of thermal energy storage systems can be lower depending on geographical location and thermophysical properties [1]. In any phase change material (PCM), the rate of absorption and release of thermal energy increases with higher thermal conductivity [2]. However, leading PCMs like paraffin waxes have thermal conductivities on the order of 0.1 W/mK. Glauber’s salt, a salt hydrate, possesses a relatively high thermal conductivity but suffers from supercooling and phase segregation. Research has shown that polymer-based composites with salt hydrates can prevent phase segregation and limit supercooling over hundreds of cycles [3] though the thermal conductivity of the composite may be reduced leading to low energy storage rates. Here, we report thermal conductivity measurements on polymer hydrogels which are stable in solutions of Glauber’s salt. To measure thermal conductivity of soft gels, we use a 3-omega method [4] where heat flows bidirectionally through a substrate below and through a gel placed on top. This suits measurements of both gels and liquids unlike transient plane heat source or transient hot wire methods. We report measurements of thermal conductivity of salt-hydrogel complexes as a function of temperature, cross-linking and salt concentration. In addition, we present modeling of thermal conductivity of sodium sulfate solutions with available theories. We finally employ polarized optical microscopy to observe inside the transparent hydrogel networks, showing that precipitating crystals may influence thermal conductivity measurements for high salt concentrations. Our thermal conductivity results combined with enthalpy of phase change, mass density and specific heat are essential to accurately design thermal storage systems for energy-efficient buildings.

composite, thermal conductivity↗

Thermal design analysis for SuperTruck II lithium-titanate battery pack

Here, this paper presents a systematic thermal management analysis for a new lithium-titanate-oxide battery pack to be installed in a SuperTruck II, Class 8 hybrid truck. The authors investigate the feasibility of mounting the battery pack inside the vehicle and air-cooling it with fans supplying conditioned air from the cabin. Moreover, the cells within each module are to be immersed in a heat-transfer fluid to improve temperature homogeneity. A multi-stage thermal analysis is performed to ensure adequate thermal regulation of the battery under the proposed design, considering installation and operational constraints. The authors perform thermal and electrical characterization tests on a single cell. Results inform the development of computational fluid dynamics models of the cells, modules, and pack. Module-level analysis shows that the cell temperatures can be maintained below their upper operational limit of 55 °C with module wall-to-air heat transfer coefficients between 10 and 20 Wm -2 K -1 . Pack-level analysis of distinct configurations determines effective airflow paths for adequate heat transfer and delivers a final battery pack design that achieves sufficient cooling and temperature uniformity. A reduced-order electrothermal model is developed to rapidly predict the transient battery performance and to develop a temperature control strategy.

25 ENERGY STORAGE↗

Design of a Combined Proportional Integral Derivative Controller to Regulate the Temperature Inside a High-Temperature Tubular Solar Reactor

Solar fuels are proven to be promising candidates for thermochemical energy storage. However, the transient nature of solar radiation is an obstacle to maintaining a stable operational temperature inside a solar reactor. To overcome this challenge, the temperature of a solar reactor can be regulated by controlling the incoming solar radiation or the feedstock flowrate inside the reactor. In this work, a combined proportional integral derivative (PID) controller is implemented to regulate the temperature inside a high-temperature tubular solar reactor with counter-current flowing gas/particles. The control model incorporates two control systems to regulate incoming solar radiation and gas flow simultaneously. The design of the controller is based on a reduced-order numerical model of a high-temperature tubular solar reactor that is vertically oriented with an upward gas flow and downward particle flow. The reactor receives heat circumferentially through its wall over a finite segment of its length. Formulation of the heat transfer model is presented by applying the energy balance for the reactor tube and considering heat and mass transfer inside. A set of governing differential equations are solved numerically by using the finite volume method to obtain reactor wall, particles, and gas temperatures along the reactor length with various boundary conditions. Simulation results are used to tune the PID controller parameters by utilizing the Ziegler–Nichols tuning method. Both the simulation results and the controller performance are visualized on the LABVIEW platform. The controller is challenged to track different temperature setpoints with different scenarios of transient solar radiation. Here, the performance of the PID controller was compared to experimental results obtained from an industrial PID controller embedded in a 7 kW electric furnace. Results show that the combined PID controller is successful in maintaining a stable temperature inside the reactor by regulating the incoming solar radiation and the flowrate via small steady-state error and reasonable settling time and overshoot.

14 SOLAR ENERGY↗

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage (CRADA Report)

NREL and Element 16 collaborated on sulfur thermal energy storage modeling using NREL’s high performance computing (HPC) resources to assist its application in industrial processes. Industrial process heat (IPH) accounts for ~70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately, 1500 TWht (~60% Terawatt hour thermal) of IPH demand is in the temperature range of 100-300°C. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. This HPC4EI project facilitated Element 16’s development of low-cost and high-impact molten sulfur TES for dispatchable IPH. The development of a high-fidelity model validated by experimental data and HPC simulations enabled the successful resolution of the complex interplay between fluid dynamics and heat transfer processes during transient operation of sulfur TES, overcoming the numerical challenges posed by the non-linear temperature-dependent physical properties of sulfur. The project helped accelerate Element 16’s molten sulfur TES product design and support its broad applications.

25 ENERGY STORAGE↗

The encasulation effect on thermal performance of micro-encapsulated phase change materials during energy absorption

Understanding the performance of the thermal storage materials is required to successfully meet the needs of next generation concentrating solar power projects. This study focuses on the thermal characteristics of the micro-encapsulated phase change material (MEPCM) with metallic core under high temperature. A numerical model was developed and verified to study the transient heat absorption process of MEPCMs made of tin, aluminum and copper core. Convective heat transfer rate, total energy absorption and volumetric energy density were presented. The effect of encapsulation on theses variables was analyzed under various boundary and initial conditions. Its impact on characteristic time variables was discussed. The method presented can be adopted to analyzing both manufacturing and application processes for MEPCMs of various compositions, during which melting and solidification happen. Given certain operating conditions, the results from the model can be referenced to obtain optimal thermal performance of these materials.

Benner, Jingru↗

A reduced-order modeling of a tubular solar reactor for long duration thermochemical energy storage

The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.

14 SOLAR ENERGY↗

Transient Pressure Interference During CO2 Injection in Saline Aquifers

Abstract CO2 injection in subsurface geological formations (e.g. deep saline aquifers) causes pressure perturbations over a large area surrounding the injection well. Observation wells are widely considered in geologic CO2 storage (GCS) projects where the pressure perturbation induced by CO2 injection is measured. In this work, we use analytical and numerical modeling tools along with field data to examine the pressure behavior in GCS projects before and after CO2 arrival at an observation well. Prior to CO2 arrival, a baseline pressure trend is established which corresponds to single-phase brine flow across the observation well (approximated by Theis solution). Therefore, analysis of early-time pressure data is straightforward, provides the single-phase flow characteristics (mobility and storativity), and helps establishing a baseline pressure change that can be extended beyond the single-phase flow period at the observation well. Upon CO2 arrival, a departure from this baseline trend is expected. For the pressure to detect the CO2 arrival at an observation well, the departure from baseline pressure behavior must be significant and well above the background noise levels. We use existing analytical models to determine the strength of the expected pressure departure signal from the baseline trend upon CO2 arrival. The strength of the expected pressure departure is found to be directly proportional to the mobility ratio. Accordingly, we establish a criterion to determine whether the pressure at an observation well can detect the CO2 arrival. We present an analysis approach through application to synthetic and field data and show the characteristic pressure behavior before and after CO2 arrival. We show that while generally the pressure can be either above or below the expected baseline pressure trend, it would be likely above the baseline upon CO2 arrival. This is because the mobility ratio becomes less than unity after CO2 arrival. We show that depending on the reservoir characteristics, changes in the pressure trend may or may not be sufficient to detect the CO2 arrival.

Engineering↗

Deep learning model for fast, science-based forecasting of fluid migration along faults in geologic carbon storage scenarios

Effective long-term geologic storage depends on robust site selection and credible, science-based forecasting of subsurface behavior to ensure storage integrity. For this work, we develop a deep learning–based reduced-order model (ROM) to quantify potential carbon dioxide (CO₂) and brine migration through geological faults. The ROM combines a Transformer model for binary classification and a Stacked Ensemble for regression, trained on a comprehensive dataset generated from 1400 physics-based reservoir simulations. Key geologic and operational parameters—including fault geometry, reservoir structure, and injection conditions—were systematically varied to capture a wide range of fluid migration scenarios. The ROM accurately predicts the onset of migration, cumulative migration volumes of both CO₂ and brine, and associated migration rates, as compared to an independent set of validation simulations, while significantly reducing computational cost compared to traditional simulation methods. Model performance was evaluated across diverse fault configurations, revealing that shallow reservoir geometry and fault angle are among the most influential factors governing migration behavior. Sensitivity analysis using SHapley Additive exPlanations (SHAP) provided interpretability, revealing distinct patterns in how geological and operational features drive transient versus cumulative migration outcomes. The ROM’s ability to rapidly simulate fault migration scenarios enables efficient sensitivity analyses, scenario evaluations, and decision support for site selection and monitoring design. This approach enhances the safety, scalability, and long-term operational performance of geologic carbon storage (GCS) systems by providing a robust, interpretable tool for predicting subsurface fluid migration and assessing fault-related migration potential.

42 ENGINEERING↗

Low-order aquifer thermal energy storage model for geothermal system simulation

This paper presents a low-order aquifer thermal energy storage (ATES) model for simulation of combined subsurface and above-surface energy systems. The model is included in the Modelica IBPSA Library, which is a free open-source library with basic models for building and district energy and control systems. The model uses a lumped-component method, in which the transient conductive-convective heat and mass transfer equation is radially discretized. To verify the accuracy of the model, we present an inter-model comparison from a simulation test suite. Results show that the Modelica ATES model is in good agreement, with a normalized mean bias error for yearly variation of aquifer temperatures of 1.6×10−2 and 9×10−5 at 1 m and 10 m distance from the well.

Maccarini, Alessandro↗