SEARCH · Engineering Papers

Results for “Transient Modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,441 records · Page 80

All sources NASA NTRS DOE OSTI

NASA NTRS2023

Hybrid-Electric Aero-Propulsion Controls Testbed Results

NASA is supporting the development of Electrified Aircraft Propulsion (EAP) technology due to its potential to reduce aircraft fuel burn, emissions, and noise as well as improving safety and performance. One focus of this research is the electrification of conventional turbomachinery propulsion systems, which offers ways to improve the performance and operability of turbine-engine powered aircraft through the addition of electro-mechanical systems. These hybrid-electric turbine engines provide additional actuation and energy management control opportunities for improving stability and transient response behavior. This paper summarizes the results of a Hardware-in-the-Loop (HIL) test performed at the NASA Electric Aircraft Testbed (NEAT) during the summer of 2022. The test demonstrates the feasibility and performance of an advanced energy management control strategy by integrating a simulated turbofan engine with scaled electro-mechanical hardware. A full-scale real-time reference model of a geared turbofan was run alongside a scaled electro-mechanical system representing the electrified turbofan components operating at a megawatt-scale power level. The model was interfaced with the hardware through a novel closed-loop control and scaling algorithm that emulated the dynamic speed and torque response of the turbofan shafts. The control strategy was implemented on the electrical machines connected to the emulated turbomachinery shafts. The results from the testbed are compared against simulations that predict the testbed and geared turbofan model operation. The energy management control strategy successfully changed the operating point of the engine model and improved its stability during throttle transients. These results also demonstrate the success of the novel closed loop control and scaling approach for emulating turbomachinery and elevate the Technology Readiness Level (TRL) of the energy management control strategy.

Aeronautics↗

DOE OSTI2021

SN 2018ijp: the explosion of a stripped-envelope star within a dense H-rich shell?

In this work, we discuss the outcomes of the follow-up campaign of SN 2018ijp, discovered as part of the Zwicky Transient Facility survey for optical transients. Its first spectrum shows similarities to broad-lined Type Ic supernovae around maximum light, whereas later spectra display strong signatures of interaction between rapidly expanding ejecta and a dense H-rich circumstellar medium, coinciding with a second peak in the photometric evolution of the transient. This evolution, along with the results of modeling of the first light-curve peak, suggests a scenario where a stripped star exploded within a dense circumstellar medium. The two main phases in the evolution of the transient could be interpreted as a first phase dominated by radioactive decays, and a later interaction-dominated phase where the ejecta collide with a pre-existing shell. We therefore discuss SN 2018jp within the context of a massive star depleted of its outer layers exploding within a dense H-rich circumstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2003

Magnetic Field Effects on Convection and Solidification in Normal and Microgravity

It has been well understood that convective flows induced by g-jitter forces associated with spacecraft are responsible for defects formation and irregularity in product quality during melt growth of single crystals in microgravity. This research is concerned about numerical simulations and experimental measurements for the purpose of developing a fundamental understanding of the g-jitter induced fluid flows and their effects on solidification in microgravity with and without the presence of additional damping forces that are derived from the applied DC magnetic fields. The numerical models include both 2-D and 3-D transient fluid flow, heat transfer, mass transfer and solidification under the combined action of g-jitter and magnetic fields. Numerical simulations using both the 2-D and 3-D models are conducted for both idealized, synthesized and real g-jitter forces, and 2-D simulations are tested against the experimental measurements taken on the thermal oscillator. 2-D solidification models have also been developed and simulations are conducted. Results show that the numerical model predictions compare well with the measurements. Analysis of these results illustrates that an applied magnetic field can have a drastic influence on the convective flows induced by g-jitter and can be particularly useful to suppress the effects resulting from the spiking of g-jitter signatures, which are considered the most detrimental effects on quality of crystals grown in space. Work in progress includes developing 3-D numerical models for solidification phenomena with the presence of both g-jitter and magnetic fields and measurements of flow fluid and its effect on solidification in both transparent fluid and low melting point melts to verify numerical predictions.

Li, Ben Q.↗

DOE OSTI2025

Progress Towards the Validation of a new RELAP5-3D model of the High Temperature Test Facility

Validation is a key step in the development of any type of systems model. As the next generation of reactors approaches, the need for codes that have been validated for these new types of systems continues to grow. An example of a prominent option is the Reactor Excursion Leak Analysis Program (RELAP5-3D), developed by Idaho National Laboratory. This code was developed for the purpose of systems level thermal-hydraulic modeling of light water reactors (LWRS) and postulated transients that can occur in LWRS.RELAP5-3D has been substantially validated against LWR data. Due to its long history as a reactor safety analysis tool, there has been an effort to adapt RELAP5-3D for the purposes of advanced reactor concepts such as prismatic high-temperature gas-cooled reactors (HTGRs). However, RELAP5-3D has not nearly been validated and verified for HTGRs to the degree of LWRs, warranting verification and validation opportunities with computational benchmarks and existing experimental facilities. Examples of such facilities include the modular high-temperature gas-cooled reactor (MHTGR) 350 and the high temperature engineering test reactor (HTTR) from Japan. The MHTGR 350 is a benchmark design concept for code-to-code verification purposes; therefore, it does not provide any experimental data for validation opportunities The HTTR provides useful multiphysics validation data but does not have the in-core instruments to generate thermal-hydraulic experimental data to help with RELAP5-3D validation. Consequently, a facility that could provide key in-core temperatures for thermal-hydraulic validation was still needed. The High Temperature Test Facility (HTTF) is an integral effects facility for HTGR thermal hydraulics developed and operated by Oregon State University. HTTF represents ¼ length scale of the General Atomics MHTGR and is rated for a total power of 2.2 MW. Axially, the core consists of an upper and lower reflector and 10 blocks, numbered from bottom to top (Block 1 is right above lower reflector). The core is heated via graphite resistive heater rods, with respective channels distributed throughout the core. The primary coolant is helium and heat can radiate out of the core to the reactor cavity cooling system (RCCS), which is cooled by water. The primary purpose of the facility is to investigate pressurized conduction cooldown (PCC) and depressurized conduction cooldown (DCC) transients, which are also referred to as the pressurized and depressurized loss of forced cooling respectively. Two experiments were chosen to perform the validation study with a RELAP5-3D model of HTTF. These experiments are PG-27 (PCC) and PG-29 (DCC). These were chosen based off of the quality of available experimental data before and during the experiment which led to their inclusion in the HTGR Thermal Hydraulics Benchmark.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Improvement of two-phase closure models in CTF using Bayesian inference

Under the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, extensive capabilities have been developed in CTF to analyze light-water reactors (LWRs) for normal operating conditions, departure from nucleate boiling (DNB), and system transients. However, further improvements are required in the modeling and simulation of boiling water reactors (BWRs), which is a focus of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, CTF validation results were used to optimize selected modeling coefficients by calibrating to experimental data using a Bayesian inference approach. Here, calibration studies were conducted to improve (vapor) void fraction prediction without worsening the two-phase pressure drop prediction, as well as to improve the two-phase pressure drop prediction. Calibration was performed for interfacial drag and wall shear models. Surrogates were developed to alleviate the computational expense required for sampling the parameter space using Markov chain Monte Carlo (MCMC). An assessment performed with calibrated models demonstrated an improvement of CTF in its prediction of key parameters such as void fraction and two-phase pressure drop.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2001

Simulation of the Effect of Realistic Space Vehicle Environments on Binary Metal Alloys

Simulations that assess the effect of space vehicle acceleration environments on the solidification of Pb-Sb alloys are reported. Space microgravity missions are designed to provide a near zero-g acceleration environment for various types of scientific experiments. Realistically. these space missions cannot provide a perfect environment. Vibrations caused by crew activity, on-board experiments, support systems stems (pumps, fans, etc.), periodic orbital maneuvers, and water dumps can all cause perturbations to the microgravity environment. In addition, the drag on the space vehicle is a source of acceleration. Therefore, it is necessary to predict the impact of these vibration-perturbations and the steady-state drag acceleration on the experiments. These predictions can be used to design mission timelines. so that the experiment is run during times that the impact of the acceleration environment is acceptable for the experiment of interest. The simulations reported herein were conducted using a finite element model that includes mass, species, momentum, and energy conservation. This model predicts the existence of "channels" within the processing mushy zone and subsequently "freckles" within the fully processed solid, which are the effects of thermosolutal convection. It is necessary to mitigate thermosolutal convection during space experiments of metal alloys, in order to study and characterize diffusion-controlled transport phenomena (microsegregation) that are normally coupled with macrosegregation. The model allows simulation of steady-state and transient acceleration values ranging from no acceleration (0 g). to microgravity conditions (10(exp -6) to 10(exp -3) g), to terrestrial gravity conditions (1 g). The transient acceleration environments simulated were from the STS-89 SpaceHAB mission and from the STS-94 SpaceLAB mission. with on-orbit accelerometer data during different mission periods used as inputs for the simulation model. Periods of crew exercise, quiet (no crew activity), and nominal conditions from STS-89 were used as simulation inputs as were periods of nominal. overboard water-dump, and free-drift (no orbit maneuvering operations) from STS-94. Steady-state acceleration environments of 0.0 and 10(exp -6) to 10(exp -1) g were also simulated, to serve as a comparison to the transient data and to assess an acceptable magnitude for the steady-state vehicle drag

Westra, Douglas G.↗

NASA NTRS1990

An experimental study of fault propagation in a jet-engine controller

An experimental analysis of the impact of transient faults on a microprocessor-based jet engine controller, used in the Boeing 747 and 757 aircrafts is described. A hierarchical simulation environment which allows the injection of transients during run-time and the tracing of their impact is described. Verification of the accuracy of this approach is also provided. A determination of the probability that a transient results in latch, pin or functional errors is made. Given a transient fault, there is approximately an 80 percent chance that there is no impact on the chip. An empirical model to depict the process of error exploration and degeneration in the target system is derived. The model shows that, if no latch errors occur within eight clock cycles, no significant damage is likely to happen. Thus, the overall impact of a transient is well contained. A state transition model is also derived from the measured data, to describe the error propagation characteristics within the chip, and to quantify the impact of transients on the external environment. The model is used to identify and isolate the critical fault propagation paths, the module most sensitive to fault propagation and the module with the highest potential of causing external pin errors.

Choi, Gwan Seung↗

NASA NTRS1986

Adaptive finite element methods for compressible flow problems

Some recent work on adaptive FEMs for solving transient Euler equations in two-dimensional domains is summarized. The formulation of an FEM model of the Euler equations is shown, and the application of the adaptive strategies to data management schemes is addressed. Sample numerical results from the application of the model and strategies to the flow over a step and to transient cases are given.

Oden, J. T.↗

NASA NTRS2003

Transient Simulation of the Integrated Powerhead Demonstrator (IPD) Rocket Engine

The IPD rocket engine is the product of a joint Air Force/NASA program to demonstrate the concept of a full-flow staged combustion power cycle. In this type of rocket engine cycle, both the entire fuel flow and the entire oxidizer flow are combusted in one of two pre-bumers and used to drive two turbopumps, thus utilizing the flow work of the total propellant flow. The basis of the IPD engine program relies on integrating newly developed and pre-existing hardware to demonstrate the component and material technologies to make this concept feasible, while simultaneously saving development time and costs. To provide insight to the project team and contractors during engine development and test phases, the engine system was modeled at MSFC using the Rocket Engine Transient Simulation (ROCETS) software to analyze system performance and determine component integration issues. The ROCETS software is used extensively a MSFC to perform steady-state power-balances and transient simulations of thermodynamic power and general fluid systems. The software is favored for its capability to solve large systems of non-linear equations, its librarie of fluid properties and flow devices, its flexibility to modify existing code to improve the physics-derived approximations of real fluid thermodynamic behavior, and the ability to add unique system constraints. The purpose of this paper is to present the methodology used to model the IPD engine system, detail the pitfalls encountered with the software, and explain the approximations made to more accurately represent engine component and fluid behavior. Engine system performance output from the model will be presented and explained in comparison with real fluid and component behavior.

Guidos, Mike↗

NASA NTRS2012

A Turbine Based Combined Cycle Engine Inlet Model and Mode Transition Simulation Based on HiTECC Tool

An inlet system is being tested to evaluate methodologies for a turbine based combined cycle propulsion system to perform a controlled inlet mode transition. Prior to wind tunnel based hardware testing of controlled mode transitions, simulation models are used to test, debug, and validate potential control algorithms. One candidate simulation package for this purpose is the High Mach Transient Engine Cycle Code (HiTECC). The HiTECC simulation package models the inlet system, propulsion systems, thermal energy, geometry, nozzle, and fuel systems. This paper discusses the modification and redesign of the simulation package and control system to represent the NASA large-scale inlet model for Combined Cycle Engine mode transition studies, mounted in NASA Glenn s 10- by 10-Foot Supersonic Wind Tunnel. This model will be used for designing and testing candidate control algorithms before implementation.

Csank, Jeffrey T.↗

NASA NTRS2012

A Turbine Based Combined Cycle Engine Inlet Model and Mode Transition Simulation Based on HiTECC Tool

An inlet system is being tested to evaluate methodologies for a turbine based combined cycle propulsion system to perform a controlled inlet mode transition. Prior to wind tunnel based hardware testing of controlled mode transitions, simulation models are used to test, debug, and validate potential control algorithms. One candidate simulation package for this purpose is the High Mach Transient Engine Cycle Code (HiTECC). The HiTECC simulation package models the inlet system, propulsion systems, thermal energy, geometry, nozzle, and fuel systems. This paper discusses the modification and redesign of the simulation package and control system to represent the NASA large-scale inlet model for Combined Cycle Engine mode transition studies, mounted in NASA Glenn s 10-foot by 10-foot Supersonic Wind Tunnel. This model will be used for designing and testing candidate control algorithms before implementation.

Csank, Jeffrey↗

DOE OSTI2021

BISON microstructure-based pulverization criterion in high burnup structure

To improve the economics of commercial nuclear power production, utilities are seeking to increase the allowable burnup limit of UO$_2$ fuel. One of the main factors that contributes to the current burnup limit of 62 GWd/MTU in commercial light water reactors (LWRs) is the risk of fine fragmentation or pulverization during a loss of coolant accident (LOCA). Pulverization primarily occurs at high burnups, especially when the high burnup structure (HBS) has formed. To allow the industry to pursue increased burnup and develop mitigation strategies, it is essential to have improved capability to predict the onset of pulverization. However, the mechanism of pulverization is not well understood, and the existing predictive capabilities implemented in the BISON fuel performance code are empirical in nature. In this report, mesoscale simulations are used to improve understanding of the formation mechanism of the HBS and how it responds during a LOCA transient, and inform development of a BISON pulverization criterion. A phase-field model was used to simulate the evolution of bubble pressure as a result of HBS formation. The simulations showed that gas atoms diffuse from grain interiors to the new grain boundaries created during HBS formation, and diffuse rapidly along these grain boundaries to reach existing bubbles. This causes an increase in bubble pressure in existing bubbles, leading to bubble growth during steady-state operation. To simulate the response of HBS bubbles to a LOCA transient, a newly developed phase-field model was used; in agreement with preliminary results from FY20, bubble size did not change significantly during the duration of the transient. A phase-field fracture model was used to study fragmentation patterns in the HBS, including using input from the phase-field model as initial conditions. Phase-field fracture simulations were used to determine a pulverization criterion for BISON. A function for the critical pressure for pulverization to occur was fit to data from the phase-field fracture simulations, and this function was implemented as a material property in BISON. For comparison, an analytical criterion for pulverization was developed and implemented within the same material property.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1981

Analytical performance prediction for large antennas

A linear mathematical modeling of the static and dynamic characteristics of a deployable offset wrap rib antenna was attempted using small motions in an effort to simulate the surface distortion resulting from thermal loading and transient oscillations from control. Both cyclic symmetry and model synthesis were used for parametric analysis. Although the ribs can be modeled as slender beams varying cross sections and equivalent linear material properties, the mesh poses a problem since its equivalent stress strain relations are highly nonlinear and sensitive to the biaxial prestressed state. Consequently, the macroscopic material properties were measured experimentally in the neighborhood of the design prestress. An equivalent tangent modulus technique was then adopted about this equilibrium state and assumed to be valid within a small range of incremental stresses.

El-Raheb, M.↗

DOE OSTI2024

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2024

Improved Modelling of SEP Event Onset Within the WSA–Enlil–SEPMOD Framework

Multi-spacecraft observations of solar energetic particle (SEP) events not only enable a deeper understanding and development of particle acceleration and transport theories but also provide important constraints for model validation efforts. However, because of computational limitations, a given physics-based SEP model is usually best suited to capture a particular phase of an SEP event, rather than its whole development from onset through decay. For example, magnetohydrodynamic (MHD) models of the heliosphere often incorporate solar transients only at the outer boundary of their so-called coronal domain – usually set at a heliocentric distance of 20–30 R⊙. This means that particle acceleration at coronal mass ejection (CME)-driven shocks is also computed from this boundary onwards, leading to simulated SEP event onsets that can be many hours later than observed, since shock waves can form much lower in the solar corona. In this work, we aim to improve the modelled onset of SEP events by inserting a “fixed source” of particle injection at the outer boundary of the coronal domain of the coupled WSA–Enlil 3D MHD model of the heliosphere. The SEP model that we employ for this effort is Solar Energetic Particle MODel (SEPMOD), a physics-based test-particle code based on a field line tracer and adiabatic invariant conservation. We apply our initial tests and results of SEPMOD’s fixed-source option to the 2021 October 9 SEP event, which was detected at five well-separated locations in the inner heliosphere – Parker Solar Probe, STEREO-A, Solar Orbiter, BepiColombo, and near-Earth spacecraft.

Solar eruptions↗

NASA NTRS2018

Improved Mars Helicopter Aerodynamic Rotor Model for Comprehensive Analyses

The Mars Helicopter is part of the NASA Mars 2020 rover mission scheduled to launch in July of 2020. Its goal is to demonstrate the viability and potential of heavier-than-air vehicles in the Martian atmosphere. Ultimately, it aims to bridge the resolution gap between orbiters and the rover as well as allow access to otherwise inaccessible regions. The low density of the Martian atmosphere and the relatively small-scale rotor result in very low Reynolds number flows. The low density and low Reynolds numbers reduce the lifting force and lifting efficiency, respectively. This paper describes the generation of the improved Mars Helicopter aerodynamic rotor model. The goal is to generate a performance model for the Mars Helicopter rotor using a free wake analysis, since this has a low computational cost for design. The improvements in the analysis are two-fold and are expanded on from two prior publications. First, the fidelity of the simulations is increased by performing higher-order two-dimensional time-accurate OVERFLOW simulations allowing for higher accuracy aerodynamic coefficients and a better understanding of the boundary layer behavior as well as its transient features. Second, a version of the model is generated to duplicate the exact testing conditions in the 25-ft. diameter Space Simulator at the Jet Propulsion Laboratory, which allows for better correlation of rotor performance figures. Previous work correlated performance with that test, but did not consider the higher temperatures in the experiment compared to those of the Martian atmosphere. The higher temperatures in the experiment are expected to give conservative performance estimates, as they give rise to an increase in speed of sound and decrease in observed Reynolds numbers.

Analyses↗

NASA NTRS2018

Improved Mars Helicopter Aerodynamic Rotor Model for Comprehensive Analyses

The Mars Helicopter is part of the NASA Mars 2020 rover mission scheduled to launch in July of 2020. Its goal is to demonstrate the viability and potential of heavier-than-air vehicles in the Martian atmosphere. Ultimately, it aims to bridge the resolution gap between orbiters and the rover as well as allow access to otherwise inaccessible regions. The low density of the Martian atmosphere and the relatively small-scale rotor result in very low Reynolds number flows. The low density and low Reynolds numbers reduce the lifting force and lifting efficiency, respectively. This paper describes the generation of the improved Mars Helicopter aerodynamic rotor model. The goal is to generate a performance model for the Mars Helicopter rotor using a free wake analysis, since this has a low computational cost for design. The improvements in the analysis are two-fold and are expanded on from two prior publications. First, the fidelity of the simulations is increased by performing higher-order two-dimensional time-accurate OVERFLOW simulations allowing for higher accuracy aerodynamic coefficients and a better understanding of the boundary layer behavior as well as its transient features. Second, a version of the model is generated to duplicate the exact testing conditions in the 25–ft. diameter Space Simulator at the Jet Propulsion Laboratory, which allows for better correlation of rotor performance figures. Previous work correlated performance with that test, but did not consider the higher temperatures in the experiment compared to those of the Martian atmosphere. The higher temperatures in the experiment are expected to give conservative performance estimates, as they give rise to an increase in speed of sound and decrease in observed Reynolds numbers.

Grip, Håvard F.↗

DOE OSTI2023

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↗