Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time transfer”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Systems Engineering Analysis of National Ignition Facility Industrial Controls Systems and Safety Interlock Systems Remote Input/Output Networking Migration from ControlNet to EtherNet/IP

The ControlNet industrial communications protocol and modules used in the Industrial Control System (ICS) and Safety Interlock System (SIS) at the National Ignition Facility (NIF) are no longer necessary and the ICS and SIS would be better served by migrating the communications structure to use EtherNet/Industrial Protocol (IP) and EtherNet bridge modules instead. By the admission of the vendor of ControlNet hardware, Rockwell Automation, in literature by Bill Petro [1], “Moving forward, customers will be able to optimize their asset utilization better using EtherNet/IP protocol than with ControlNet.” The NIF is one of the key elements of the Inertial Confinement Fusion (ICF) program at Lawrence Livermore National Laboratory (LLNL), a federally funded research and development center (FFRDC). The NIF contains the systems and provides the operational capacity to perform ICF, high energy density (HED), and discovery science experiments utilizing 192 individual beamlines, a host of diagnostics, and all the industrial systems required to facilitate these beamlines and diagnostics. The industrial systems are governed by the ICS and SIS, with the ICS providing control and the SIS providing monitoring and permissives. Construction on the NIF began in 1997 and was certified complete in 2009 and, as a result, the ICS and SIS were developed during this time using the tools that were available then. This includes the communications structure and protocols for these systems, much of which was, and still is, ControlNet. ControlNet, particularly during the time that the ICS and SIS were being built, has several attractive features. ControlNet hardware is exclusive to Rockwell Automation, which was the automation hardware chosen for the ICS and SIS. One feature that could be considered an advantage or a disadvantage depending on the communication needs of the system is that ControlNet also utilizes no active network components, excluding repeaters which are not always necessary. According to the architect of the ICS system at the NIF, Gordon Lau, one of the most attractive features of the ControlNet protocol during development of the ICS and SIS was that it is deterministic, providing timing of data transfer that is executed exactly as it is defined by the developer.

42 ENGINEERING↗
DOE OSTI2023

Spark Channel Dynamics of Electrostatic Discharges (Rev. 1) [Slides]

To model spark discharges, this work implemented, verified, and validated: 1-D radial conservative Lagrangian hydrodynamic scheme (with heat conduction); Eddington/P1 radiative transfer approximation, time-implicit scheme; RLC circuit solver, time-implicit scheme (adaptable to other circuit types). Novel analytics EOS fit for air up to temperatures of ~150,000 K. Novel "seed" electron concept to initialize finite conductivity in spark without pre-heating the air in the channel, as typically used in spark modeling literature. Reasonable, if conservative (from an engineering viewpoint), agreement with experimental data for indirect (circuit) and direct (interferometry) measurements of spark. Many potential avenues of future work: Two-temperature dryodynamics (requires two-temperature EOS development); 2-D axisymmetric geometry for axial (z) variation, modeling delectric electrodes; Kinetic equations for time-dependent charged species populations.

42 ENGINEERING↗
DOE OSTI2022

EPCAPE – Proton Transfer Reaction Mass Spectrometer Field Campaign Report

The Eastern Pacific Cloud Aerosol Precipitation Experiment – Proton Transfer Reaction Mass Spectrometer (EPCAPE-PTRMS) campaign was created in response to the urgent need expressed by EPCAPE co-investigators for measurements of volatile organic compounds (VOCs) during the EPCAPE intensive operational period. A proton transfer reaction time-of-flight mass spectrometer (PTR-ToF-MS, model 8000, Ionicon Analytik) was deployed to Mt. Soledad, San Diego, California, on 15 May 2023 for the measurement of VOCs. The instrument was installed in the same Atmospheric Radiation Measurement (ARM) user facility guest container that the Smith group used for the EPCAPE – Ultrafine Particle Properties (EPCAPE-UPP) campaign. The Smith group provided the sample inlet and operated and maintained the instrument. A complete VOC mass spectrum was recorded every 10 sec, continuously. Calibrations using a traceable VOC mixture were performed onsite at the end of the campaign.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Filament Formation and Melt Spinning of Coal-Based Mesophase Pitch for Carbon Fiber Production

High-performance carbon fibers excel as high specific strength and modulus materials, and are utilized in composites applications ranging across aerospace, automotive, energy, infrastructure, and sports equipment sectors. Compared to PAN-based carbon fiber, mesophase pitch-based carbon fiber has lower tensile strength but provides a higher modulus, higher thermal conductivity, and a lower cost potential. However, challenges in stable, continuous melt spinning processing are a serious limiting factor. Unlike typical melt spinning of long, linear chain polymers, mesophase pitch is a comprised of relatively shorter polycyclic aromatic hydrocarbons which form a liquid crystal. In its ‘green’ or ‘as-spun’ state, the fibers are very fragile. Moreover, its temperature of processing is quite high, often approaching 400ºC. Complex flow dynamics combined with short length and time scale heat transfer render its melt spinning a formidable processing challenge. Improved understanding of the root causes for nascent filament breakage and spinning instabilities are needed. This work aims to determine fundamental phenomena that govern mesophase pitch filament formation by comparison of required draw force and uninterrupted spinning minutes with applied draw down ratio. The effect of perturbations on spinning stability are analyzed using capillary rheometry, microscopy of elongating filaments, and filtration.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Atmospheric Gases as Quantitative Indicators of Urbanization Across Land-Use Gradients

The urban heat island (UHI) effect is a significant phenomenon characterized by higher surface temperatures in cities compared to surrounding rural areas. This change in the Earth's energy balance presents considerable challenges for urban populations, increasing heat exposure, straining energy infrastructure, and exacerbating health vulnerabilities. Quantifying UHI intensity is challenging due to blurred distinctions between urban and rural areas as land-use patterns evolve. To address this challenge, the present study builds on the framework established by Salvador et al. (2025) and investigates the use of atmospheric constituents, specifically volatile organic compounds (VOCs) and inorganic gases, as tracers of urbanization levels. A key improvement in this approach is the inclusion of a designated true rural site for comparison. Field measurements were conducted at two locations in the Southeastern United States during the summer of 2025. Inorganic gases were monitored continuously, while VOC samples were collected biweekly and analyzed using a proton transfer reaction time-of-flight mass spectrometer. The urban site exhibited elevated concentrations of anthropogenic inorganic gases, with carbon monoxide and nitrogen dioxide averaging 878 ppb and 17.99 ppb, respectively. In contrast, the rural site showed lower levels of these pollutants, 796 ppb for carbon monoxide and 12.46 ppb for nitrogen dioxide, while displaying higher concentrations of biogenic VOCs, including isoprene and monoterpenes at 1.24 ppb and 0.08 ppb, respectively. Extended mass spectral analyses revealed distinct chemical differences between the two sites, supporting the observed variability of major gases. Overall, these findings are consistent with our previous research and highlight the value of atmospheric constituents as effective indicators of urbanization, thereby improving our understanding of UHI dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Deposition Behavior of Supersaturated Silicic Acid on the Surfaces of Chlorite and Sericite around Geological Repository - 20138

Around the repository for HLW underground, cementitious materials used for the repository may complicate the migration of radionuclides in the host rock by the alteration of groundwater to a high alkaline condition. Specially, the alkaline groundwater dissolves silicate minerals of host rocks but the plenty of silicic acid dissolved from host rocks in groundwater will become supersaturated and deposit on the surface of flow-paths after decreasing the pH value to about 8 with dispersion and dilution. The deposition of silicic acid is presumed to recover and improve the barrier performance of host rocks for the radionuclide migration because of the narrowing of flow paths and the decrease of the permeability. In this study, the deposition behavior of supersaturated silicic acid was examined by the deposition experiments and the deposition rate constants, k [m/s], were evaluated from the experimental data. As the minerals of the surface of flow-paths, this study focused on chlorite and sericite, hydrothermally-altered minerals distributed to fill fractures of granite. In the experiments, deposition behavior of silicic acid on the powder of chlorite or sericite was observed after the pH of the solution was adjusted to 8 in order to make the silicic acid supersaturated. The solution of CaCl{sub 2} and NaCl were added into a silicic acid to simulate Ca{sup 2+} from the cementitious materials of the repository, Na{sup +} of the saline groundwater in Japan. From the results, the value of k was estimated to be about the order of 10{sup -10}∼10{sup -9} m/s for both chlorite and sericite and the value was almost the same as those for amorphous silica in the previous study. Therefore, the deposition of silicic acid may be assessed by similar apparent deposition rate constants even in actual underground containing various silicate minerals. Moreover, the Damkoehler number which is the dimensionless number expressed by the ratio of characteristic time of mass transfer and chemical reaction was calculated to be 10{sup 4} or more from the value of k. This value shows that the deposition reaction of silicic acid was estimated as much more predominant than the groundwater flow under the situation simulated in this study. These results suggest that the deposition of supersaturated silicic acid on the flow paths composed with various silicate minerals will be able to retard the migration of radionuclides because of narrowing micro flow-paths. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Influences of Mg Ions and Al Ions on Deposition Behavior of Supersaturated Silicic Acid - 20261

Silicic acid can change into various species such as deposit, colloidal, and soluble form under certain conditions: for example, gradual pH changes around geological disposal repositories. In particular, silicic acid deposits narrow groundwater flow paths, thereby restricting the flow of groundwater. Such effects may contribute to a delay in radionuclide migration. Previous studies have reported that silicic acid supersaturated in sodium and calcium ions can narrow the flow paths. This study examined the effect of magnesium (Mg) and aluminum (Al) ions on the deposition of supersaturated silicic acid on the solid phase. The deposition behavior of supersaturated silicic acid on the solid phase was investigated in a batch system while changing the initial supersaturation concentration, the coexisting ion concentration, and the amount of added solid phase. In the presence of Mg ions, silicic acid precipitation increased with increasing Mg ion concentration due to the influence of the electrolyte, and crosslinking reactions. The coexistence of Al ions significantly suppressed silicic acid deposition due to the formation of a covering on the solid surface during Al ion hydrolysis under conditions of high Al concentration (1.4 mM). On the basis of these experimental results, the apparent deposition rate constant of silicic acid was estimated to be in the range from 10{sup -11} to 10{sup -10} m/s. The Damkoehler number (the dimensionless number expressed by the ratio of the characteristic time of mass transfer and the chemical reaction constant) was sufficient to cause the clogging of the host rock by deposition. This suggests that supersaturated silicic acid deposition, regardless of the presence of Mg ions (and Al ions), contributes to improving barrier performance against radionuclide migration in the host rock. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Thermal Analysis of a 100 kW Polyphase Wireless Power Transfer System

Charging Electric Vehicles (EVs) fast and safely has a crucial role in the future of the EV technology. High-power Wireless Power Transfer (WPT) helps to significantly decrease the charging time. However, when the power transfer levels increase, thermal management becomes a significant challenge. The thermal design of the WPT systems needs more consideration in the design and implementation steps. This paper presents a thermal analysis of a 100 kW high-power WPT system. The thermal performance of the proposed design was evaluated at different power levels by considering the magnetic design and loss analysis. Finite Element Analysis (FEA) of the proposed design was performed and the thermal images of the implemented system were taken to prove the simulation results. The results show that, a liquid cooling design is needed for a high-power WPT systems for the long-time continuous operations of the charging pads.

Aydin, Emrullah↗
DOE OSTI2021

Ultrafast X-ray scattering offers a structural view of excited-state charge transfer

Intramolecular charge transfer and the associated changes in molecular structure in N,N'-dimethylpiperazine are tracked using femtosecond gas-phase X-ray scattering. The molecules are optically excited to the 3p state at 200 nm. Following rapid relaxation to the 3s state, distinct charge-localized and charge-delocalized species related by charge transfer are observed. The experiment determines the molecular structure of the two species, with the redistribution of electron density accounted for by a scattering correction factor. The initially dominant charge-localized state has a weakened carbon–carbon bond and reorients one methyl group compared with the ground state. Subsequent charge transfer to the charge-delocalized state elongates the carbon–carbon bond further, creating an extended 1.634 Å bond, and also reorients the second methyl group. At the same time, the bond lengths between the nitrogen and the ring-carbon atoms contract from an average of 1.505 to 1.465 Å. The experiment determines the overall charge transfer time constant for approaching the equilibrium between charge-localized and charge-delocalized species to 3.0 ps.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Adaptive Angular Quadrature Scheme for a backwards-in-time Method of Characteristics Solution to the Radiative Transfer Equation [Slides]

Radiative transfer/radiation transport are important problems to solve in astrophysics and high energy density physics. Various methods exist to solve radiation transport, such as Monte Carlo (MC), Discrete Ordinates (S N ), Method of Characteristics (MOC), and the spherical harmonics (P N ) method. Method of Characteristics requires “launching” of rays in discrete directions. Unresolved details of angular mesh create ray effects and can miss sources in the domain. Ray effects can lead to unphysical “stepping” in solution and incorrect energy deposition. Adaptive quadrature schemes can be used to detect and mitigate these effects. The Method of Characteristics (MOC) is a common method for solving hyperbolic PDEs in radiation transport and supersonic flow problems. Generally in MOC for radiation transport, virtual particles are tracked from birth to the end of a timestep. This requires interpolation to go from final location to cell averaged or corner values of angular intensity. Backwards-in-Time (BIT) particle tracking avoids this by prescribing the final position of the virtual particle at the cell nodes/corners. Angular intensities are computed at time k + 1 by launching ray back to previous timestep(s), or t = 0. Scheme allows solution to be computed as the characteristic ray is traced backwards in time.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Ultrafast spectroscopic studies of vibrational energy transfer in energetic materials

Shock-induced detonation is a key property of energetic materials (EM) that remains empirically understood. One proposed mechanism of shock-initiation in EM is “phonon up-pumping” to initiate chemical reactions, where excitation of lattice phonon modes rapidly transfers energy into intramolecular vibrations, ultimately resulting in the breaking of chemical bonds. We are developing novel ultrafast laser spectroscopy techniques to study vibrational energy transfer from phonon modes to intramolecular vibrations (phonon up-pumping), as well as competing energy transfer pathways from intramolecular vibrations to phonon modes (vibrational cooling). Through combinations of plasma- generated supercontinuum infrared, tunable near- and mid-infrared, and terahertz pulses in pump-probe spectroscopy, supplemented with ab inito simulations, we can explore the energy transfer processes on a sub-picosecond time scale to elucidate vibrational energy transfer pathways and lifetimes in EM. In this work, we highlight recent progress, including the spectral and temporal characteristics of the infrared and THz sources as well as preliminary results on select EM.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spin-optomechanical cavity interfaces by deep subwavelength phonon-photon confinement

A central goal of quantum information science is transferring qubits between space, time, and modality. Spin-based systems in solids are promising quantum memories, but high-fidelity transfer of their quantum states to telecom optical fields remains challenging. Here, we introduce a phonon-mediated interface between spins in a diamond nanobeam optomechanical crystal and telecom optical fields by a simultaneous deep-subwavelength confinement of optical and acoustic fields with mode volumes $V_{\textrm{mech}}$$/Λ^3_\textrm{p} ~ 10^{-5}$ and $V_{\textrm{opt}}$$/λ^3 ~ 10^{−3}$, respectively. This confinement boosts the spin-mechanical coupling rate of Group-IV silicon vacancy (SiV − ) centers by an order of magnitude to ~ 32 MHz while retaining high acousto-optical couplings. The optical cavity couples to the spin irrespective of the emitter’s native excited states, avoiding spectral diffusion. Using Quantum Monte Carlo simulations, we estimate heralded entanglement fidelities exceeding 0.96 between two such interfaces. We anticipate broad utility beyond diamond emitter-telecom systems to most solid-state quantum memories.

Raniwala, Hamza [Massachusetts Inst. of Technology↗
DOE OSTI2024

Real-Time Ethernet Interface for NSTX-U’s Thomson Scattering Diagnostic (2023)

Here, the multipoint Thomson scattering (MPTS) diagnostic system at the National Spherical Torus Experiment Upgrade (NSTX-U) facility is undergoing an upgrade to operate in real-time and interface with the plasma control system (PCS) for NSTX-U. Previous prototyping efforts have shown that spectral analysis and rapid calculations of electron temperature and density are possible on a real-time Linux machine when using up to a 100-Hz laser pulse repetition rate. A remaining challenge was transferring the real-time data to NSTX-U’s PCS, which utilizes the front panel data port (FPDP) protocol. The original proposed method was to convert the real-time data into analog values, but a new solution was developed to keep the output format digital by using an Ethernet controller with a field-programmable gate array (FPGA). This article focuses on a new input module that has been developed to accept incoming user datagram protocol (UDP) packets sent over Ethernet, convert into FPDP format, and integrate into the existing data stream under NSTX-U’s real-time framework.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Dynamics of carrier injection through V-defects in long wavelength GaN LEDs

The efficiency of high-power operation of multiple quantum well (QW) light emitting diodes (LEDs) to a large degree depends on the realization of uniform hole distribution between the QWs. In long wavelength InGaN/GaN QW LEDs, the thermionic interwell hole transport is hindered by high GaN barriers. However, in polar LED structures, these barriers may be circumvented by the lateral hole injection via semipolar {$10\bar{1}1$} QWs that form on the facets of V-defects. The efficiency of such carrier transfer depends on the transport time since transport in the semipolar QWs is competed by recombination. In this work, we study the carrier transfer from the semipolar to polar QWs by time-resolved photoluminescence in long wavelength (green to red) LEDs. We find that the carrier transfer through the semipolar QWs is fast, a few tens of picoseconds with the estimated room temperature ambipolar diffusion coefficient of ∼5.5 cm 2 /s. With diffusion much faster than recombination, the hole transport from the p-side of the structure to the polar QWs should proceed without a substantial loss, contributing to the high efficiency of long wavelength GaN LEDs.

Crystallographic defects↗
DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗
DOE OSTI2022

Mechanisms of shape transfer and preheating in indirect-drive double shell collisions

Implosions of Hohlraum-driven double shell targets as an alternative inertial confinement fusion concept are underway at the National Ignition Facility. The double shell system relies on a series of energy transfer processes starting from thermal x-ray absorption by the outer shell, followed by collisional transfer of kinetic energy to a heavy metal inner shell, and finally, conversion to the internal energy of the deuterium-tritium fuel. During each of these energy transfer stages, low-mode asymmetries can act to reduce the ideal transfer efficiency degrading double shell performance. Mechanisms, such as hard x-ray preheat from the Hohlraum, not only decrease the efficiency of kinetic energy transfer but may also be a source of low-mode asymmetry. In this article, we evaluate the shape transfer processes through the time of shell collision using two-dimensional integrated Hohlraum and capsule computations. We find that the dominant mode of the shape transfer is well described using a “radial impulse” model from the shape of the foam pressure reservoir. To evaluate the importance of preheat on inner shell shape, we also report on first measurements of Au L-shell preheat asymmetry in a double shell with a tungsten pusher. These measurements showed a 65% higher preheat velocity at the pole of the capsule relative to the equator. We also found that the experiments provided rigorous constraints by which to test the Hohlraum model settings that impact the amount and symmetry of Au L-shell preheat via the plasma conditions inside the outer cone Au bubble.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Coupling Subsurface and Above-Surface Models for Optimizing the Design of Borefields and District Heating and Cooling Systems

Accurate dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange for seasonal energy storage. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes without flowing groundwater. While detailed simulation tools for subsurface heat and mass transfer exist, these fall short in simulating above-surface energy systems. To support the design and operation of such systems, the study developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the subsurface. For the first, it uses the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, it uses the TOUGH simulator. The TOUGH family of codes can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. The study validated the coupled modeling approach by comparing the simulation results with one from the g-function based ground response model. It then looked into effects when the water table and the regional groundwater flow are considered in the ground, from the perspective of heat exchange between borehole and ground, and the electrical consumption of the district heating and cooling systems. To access the simulation models, please find the links in the submission: -- For coupled approach validation: see model Buildings.Fluid.Geothermal.Borefields.Examples.BorefieldsWithTough and Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from the "Modelica Building Library" resource, branch issue1495_tough_interface, commit a2667c0. -- For the study of the effect of water table: see model Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO, commit 760de49. -- For the study of the effect of regional groundwater flow: see Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO_3D, commit c2a2d2a. The coupling interface script "GrounResponse.py" can be found from the above links in the folder Buildings/Resources/Python-Sources. Also, the needed files for TOUGH simulation are in the folder Buildings/Resources/Python-Sources/ToughFiles that can be accessed through the above links. A brief description of these files is given below; detailed specifications for the first three files may be found in the TOUGH3 Users Guide (Jung et al., 2018) https://tough.lbl.gov/documentation/tough-manuals/. (1) INCON - initial conditions for each grid block (2) INFILE - main input file with material properties and control parameters (3) MESH - description of the computational grid (4) readsave - Modelica/TOUGH interface program: read the final output of TOUGH simulation after TOUGH time step and prepare for transfer to Modelica for next Modelica time step (5) readsave.inp - input parameters for program readsave (6) writeincon - Modelica/TOUGH interface program: write the output of Modelica after Modelica time step and prepare for transfer to TOUGH as initial conditions for the next TOUGH step (7) writeincon.inp - input parameters for program writeincon

15 GEOTHERMAL ENERGY↗