Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Porous flow”

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

Modeling Monthly Nitrate Concentration in a Karst Spring with and without Discrete Conduit Flow

Understanding the coupled continuum pipe-flow framework for modeling contaminant transport in karst systems is critical for protecting water resources therein. This study simulated point and non-point source nitrate concentration in a karst spring and investigated the results generated from the flow and transport model with and without discrete conduit flow. CFPv2 and CMT3D models were integrated to address the changes in nitrate concentration at a monthly scale in a karst spring, and the results were compared with that from an equivalent porous media (EPM) model with high-hydraulic conductivity (K) zones set in the main karstified area to represent conduits. The results show that the CFPv2+CMT3D model is able to describe well the recession of nitrate concentration in spring discharge, and the relatively larger deviation (slower nitrate recession) from the observed trend for the EPM model is probably a result of the limitation of utilizing high hydraulic conductivity cells to represent conduit. Moreover, simulated hydraulic heads in poorly karstified areas from the two models both show slight differences from the observations (the head RMSE values of calibration/validation for CFPv2 and MODFLOW models are 0.16 m/0.25 m and 0.26 m/0.17 m, respectively), indicating the inclusion of conduits may not affect the simulation considerably, and the lower the proportion of karstic area, the slight effects brought from the inclusion of conduits in the model. For highly karstified areas, the CFPv2+CMT3D model could provide more accurate results (head RMSE of calibration/validation for CFPv2 and MODFLOW are 0.22 m/0.06 and 0.52 m/0.47 m, respectively), showing the coupled continuum pipe-flow framework may be more appropriate for applying to highly and maturely karstified areas where the variations in the behavior of flow and contaminant transport are more affected by turbulent flow regime.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

The design of test-section inserts for higher speed aeroacoustic testing in the Ames 80- by 120-foot wind tunnel

An engineering feasibility study was made of aeroacoustic inserts designed for large-scale acoustic research on aircraft models in the 80- by 120-Foot Wind Tunnel at NASA Ames Research Center. The goal was to find test-section modifications that would allow improved aeroacoustic testing at airspeeds equal to and above the current 100 knots limit. Results indicate that the required maximum airspeed drives the design of a particular insert. Using goals of 200, 150, and 100 knots airspeed, the analysis led to a 30 x 60 ft open-jet test section, a 40 x 80 ft open-jet test section, and a 70 x 110 ft closed test section with enhanced wall lining respectively. The open-jet inserts would be composed of a nozzle, collector, diffuser, and acoustic wedges incorporated in the existing 80 x 120 ft test section. The closed test section would be composed of approximately 5-ft acoustic wedges covered by a porous plate attached to the test-section walls of the existing 80 x 120. All designs would require a double row of acoustic vanes between the test section and fan drive to attenuate fan noise and, in the case of the open-jet designs, to control flow separation at the diffuser downstream end. The inserts would allow virtually anechoic acoustics studies of large helicopter models, jets, and V/STOL aircraft models in simulated flight. Model scale studies would be necessary to optimize the aerodynamic and acoustic performance of any of the designs. Successful development of acoustically transparent walls, though not strictly necessary to the project, would lead to a porous-wall test section that could be substituted for any of the open-jet designs, and thereby eliminate many aerodynamic and acoustic problems characteristic of open-jet shear layers.

Soderman, Paul T.↗
DOE OSTI2022

Design optimization of gap distance for the capillary limitation of a heat pipe with annular-type wick structure

In this study, an experimental investigation was conducted on the rising height and contact angle of fluid in an annular wick-type heat pipe. The annular wick-type heat pipe was characterized by a small gap between the wick structure and tube wall, which compensated for the pressure drop along the porous media and created additional capillary force. To describe and model the advantage of this gap, the rising of a wetting liquid in the gap between a vertical solid plate and a mesh (with a small angle between them) was experimentally measured and analyzed. An additional experiment was performed to investigate the effect of curvature on the capillary rise using tubes and meshes of varying radii. Resultantly, we confirmed that the linear combination of the contact angles of the solid plate and mesh could be applied to calculate the rising height from the Laplace–Young equation. Furthermore, the effect of curvature on the rising height of the liquid was negligible. These results were extended to the investigation of finding the optimal gap distance for the annular wick-type heat pipe by referring to previous studies. We observed that a gap distance of 1.27 mm provided the largest permeability (K) over the effective pore radius (r eff ) value for a heat pipe with ethanol, which in turn resulted in the highest capillary limitation. For a sodium heat pipe, a gap distance of 0.84 mm resulted in the highest capillary limitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2014

NASA Bioculture System: From Experiment Definition to Flight Payload

Starting in 2015, the NASA Bioculture System will be available to the science community to conduct cell biology and microbiology experiments on ISS. The Bioculture System carries ten environmentally independent Cassettes, which house the experiments. The closed loop fluids flow path subsystem in each Cassette provides a perfusion-based method for maintain specimen cultures in a shear-free environment by using a biochamber based on porous hollow fiber bioreactor technology. Each Cassette contains an incubator and separate insulated refrigerator compartment for storage of media, samples, nutrients and additives. The hardware is capable of fully automated or manual specimen culturing and processing, including in-flight experiment initiation, sampling and fixation, up to BSL-2 specimen culturing, and the ability to up to 10 independent cultures in parallel for statistical analysis. The incubation and culturing of specimens in the Bioculture System is a departure from standard laboratory culturing methods. Therefore, it is critical that the PI has an understanding the pre-flight test required for successfully using the Bioculture System to conduct an on-orbit experiment. Overall, the PI will conduct a series of ground tests to define flight experiment and on-orbit implementation requirements, verify biocompatibility, and determine base bioreactor conditions. The ground test processes for the utilization of the Bioculture System, from experiment selection to flight, will be reviewed. Also, pre-flight test schedules and use of COTS ground test equipment (CellMax and FiberCell systems) and the Bioculture System will be discussed.

Bioculture System↗
DOE OSTI2021

Nanopore networks in colloidal silica assemblies characterized by XCT for confined fluid flow modeling

Fluid flow through nanopore structures has exhibited different behavior than that described by Darcy's law, derived from pore networks of micrometer scale. Here, pressure assemblies of colloidal silica spheres of specified diameters from 94 to 620 nanometers were prepared as model nano-porous media. The smaller silica spheres tended to form random packing, but some substantial ordering was found in the packing of the larger silica spheres. X-ray computed tomography (XCT), with a voxel resolution of 16 nm, was used to characterize the nanopore networks. The porosity of the Nano-XCT pore network (34%) was less than the total porosity of the assemblies as determined by X-ray attenuation (52%), possibly due to the limits of voxel resolution during segmentation. Fluid flow in the nanopore networks was simulated using the single-phase Lattice Boltzmann Method, and the simulated permeability was compared with the empirical and experimental values. Based on our characterization, although a well-ordered packing of silica spheres was not achieved, it was found that the nanopore networks in the colloidal silica assemblies had pore size distributions corresponding to the particle sizes. The simulated permeability was less than the experimental measurement for water flow, but the complex packing of silica spheres and surface chemistry issues need to be considered in future research.

02 PETROLEUM↗
NASA NTRS1987

Centrifuge impact cratering experiments: Scaling laws for non-porous targets

A geotechnical centrifuge was used to investigate large body impacts onto planetary surfaces. At elevated gravity, it is possible to match various dimensionless similarity parameters which were shown to govern large scale impacts. Observations of crater growth and target flow fields have provided detailed and critical tests of a complete and unified scaling theory for impact cratering. Scaling estimates were determined for nonporous targets. Scaling estimates for large scale cratering in rock proposed previously by others have assumed that the crater radius is proportional to powers of the impactor energy and gravity, with no additional dependence on impact velocity. The size scaling laws determined from ongoing centrifuge experiments differ from earlier ones in three respects. First, a distinct dependence of impact velocity is recognized, even for constant impactor energy. Second, the present energy exponent for low porosity targets, like competent rock, is lower than earlier estimates. Third, the gravity exponent is recognized here as being related to both the energy and the velocity exponents.

Schmidt, Robert M.↗
NASA NTRS2001

Numerical Simulation of One- and Two-Phase Flows in Propulsion Systems

In this report, we present some results of problems investigated during joint research between the Hampton University Fluid Mechanics and Acoustics Laboratory (HU/FM&AL), NASA GRC, and the LaRC Hyper-X Program. This work is supported by joint research between the NASA GRC and the Institute of Mechanics at Moscow State University (IM/MSU) in Russia under a CRDF grant. The main areas of current scientific interest of the HU/FM&AL include an investigation of the proposed and patented advanced methods for aircraft engine thrust and noise benefits. These methods are based on nontraditional 3D corrugated and composite nozzle, inlet, propeller and screw designs such as a Bluebell and Telescope nozzle, Mobius-shaped screw, etc. This is the main subject of our other projects, of which one is presented at the current conference. Here we analyze additional methods for exhaust jet noise reduction without essential thrust loss and even with thrust augmentation. Such additional approaches are: (1) to add some solid, fluid, or gas mass at discrete locations to the main supersonic gas stream to minimize the negative influence of strong shock waves formed in propulsion systems. This mass addition may be accompanied by heat addition to the main stream as a result of the fuel combustion or by cooling of this stream as a result of the liquid mass evaporation and boiling; (2) Use of porous or permeable nozzles and additional shells at the nozzle exit for preliminary cooling of the hot jet exhaust and pressure compensation for off-design conditions (so-called continuous ejector with small mass flow rate); and (3) to propose and analyze new effective methods of fuel injection into the flow stream in air-breathing engines. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. The FM&AL Team uses analytical methods, numerical simulations, and experimental tests at the Hampton University campus, NASA, and IM/MSU. The main outcomes during this reporting period are: (1) Publications: The AIAA Paper #01-1893 has been accepted for the AIAA/NAL-NASDA-ISAS 10th International Space Planes and Hypersonic Systems and Technologies Conference, 24-27 April 2001, Kyoto, Japan. The AIAA Paper #01-3204 has been accepted for presentation at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference being held on 08-11 July, in Salt Lake City, UT; (2) Grants and proposals: The HU/FM&AL was awarded the NASA grant NAG3-2495 in October 2000 and the CRDF award was granted to the NASA GRC-HU FM&AL and IM/MSU (Russia) in July 2000. A solicited proposal was submitted for the NASA NRA-01GRC-02 competition and two unsolicited proposals to NASA are in preparation; (3) Theory and numerical simulations: Analytical theory, numerical simulation, comparison of theoretical with experimental results, and modification of theoretical approaches, models, grids, etc. Such investigations have been conducted for three main problems: (a) Combustion efficiency optimization in the half-duct combustor system; (b) Drag reduction effects for blunt bodies with solid needles; and (c) Solid particle injection from the butt-end against a supersonic flow. The NASA CFL3D, HU/FM&AL, and IM/MSU GODUNOV codes were used and modified for solution of these problems. The codes are based on full Euler and Navier-Stokes solvers with and without nonequilibrium oxygen-nitrogen and air-hydrogen chemical reactions in laminar and turbulent gas flow regimes; (4) Experimental tests: In the IM/MSU supersonic wind tunnel, experimental tests were conducted with different number, location, and geometric parameters of solid needles mounted at the front of the butt-end in supersonic flow. Optimal parameters were determined that provide minimal butt-end drag in the stationary flow regime. Experimental and numerical simulation results are in good agreement; (5) Student Research Activity: Involvement of one graduate and two undergraduate students as research assistants in the current project.

Gilinsky, Mikhail↗
DOE OSTI2022

Prediction of local concentration fields in porous media with chemical reaction using a multi scale convolutional neural network

The study of solute transport in porous media is of interest in many chemical engineering systems. Some example applications include packed bed catalytic reactors, filtration devices, and batteries. The pore scale modeling of these systems is time consuming and may require large computing resources, for this reason computational fluid dynamics (CFD) simulations are not practical if a large number of simulations is required, like in multiscale modeling, where a model at a large scale calls for pore scale simulations. It has been shown that neural networks can be trained with a dataset of flow simulations and then predict fields orders of magnitude faster, and with less computational resources, in new domains. However, it is crucial to provide the neural network with an effective description of the domain and the undergoing operating conditions to be able to train models that generalize accurately in unseen samples. Therefore, research is needed to employ neural networks in new complex systems. The appropriate training of a network for predicting coupled flow and solute transport processes is an outstanding problem due to the complex interplay between geometry and operating conditions. In this work, we train a multi scale convolutional neural network (MSNet) with a diverse dataset of simulations of transport and chemical reaction in porous media to predict the local concentration fields in images of porous media. Our dataset contains a wide diversity of sphere pack arrangements under different operating conditions (Péclet and Reynolds numbers). Further, we train a robust model by employing different input descriptors that represent the medium and the different operating conditions of each system. Our trained model is able to provide nearly instantaneous predictions, compared to around twenty hours of the CFD workflow, with less than 3.5% error on new geometries and transport conditions. Thus the model could be easily integrated in a multiscale workflow where fast response is needed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Large-scale aeroacoustic research feasibility and conceptual design of test-section inserts for the Ames 80- by 120-foot wind tunnel

An engineering feasibility study was made of aeroacoustic inserts designed for large-scale acoustic research on aircraft models in the 80 by 120 foot Wind Tunnel at NASA Ames Research Center. The advantages and disadvantages of likely designs were analyzed. Results indicate that the required maximum airspeed leads to the design of a particular insert. Using goals of 200, 150, and 100 knots airspeed, the analysis indicated a 30 x 60 ft open-jet test section, a 40 x 80 ft open jet test section, and a 70 x 100 ft closed test section with enhanced wall lining, respectively. The open-jet inserts would be composed of a nozzle, collector, diffuser, and acoutic wedges incorporated in the existing 80 x 120 test section. The closed test section would be composed of approximately 5 ft acoustic wedges covered by a porous plate attached to the test section walls of the existing 80 x 120. All designs would require a double row of acoustic vanes between the test section and fan drive to attenuate fan noise and, in the case of the open-jet designs, to control flow separation at the diffuser downstream end. The inserts would allow virtually anechoic acoustic studies of large helicopter models, jets, and V/STOL aircraft models in simulated flight. Model scale studies would be necessary to optimize the aerodynamic and acoustic performance of any of the designs. In all designs studied, the existing structure would have to be reinforced. Successful development of acoustically transparent walls, though not strictly necessary to the project, would lead to a porous-wall test section that could be substituted for any of the open-jet designs, and thereby eliminate many aerodynamic and acoustic problems characteristic of open-jet shear layers. The larger size of the facility would make installation and removal of the insert components difficult. Consequently, scheduling of the existing 80 x 120 aerodynamic test section and scheduling of the open-jet test section would likely be made on an annual or longer basis. The enhanced wall-lining insert would likely be permanent. Although the modifications are technically feasible, the economic practicality of the project was not evaluated.

Soderman, Paul T.↗
NASA NTRS1995

On the effect of pressure, oxygen concentration, air flow and gravity on simulated pool fires

The initial development of a fire is characterized by the establishment of a diffusion flame over the surface of a the condensed fuel and is particularly influenced by gravity, with most of the gaseous flow induced by natural convection. Low initial momentum of the fuel vapor, strong buoyant flows induced by the hot post-combustion gases and consequently low values of the Froude number (inertia-gravity forces ratio) are typical of this kind of scenario. An experimental study is conducted by using a porous burner to simulate the burning of a horizontal combustible surface. Ethane is used as fuel and different mixtures of oxygen and nitrogen as oxidizer. The magnitude of the fuel injection velocities is restricted to values that will keep the Froude number on the order of 10-5, when calculated at normal gravity and pressure, which are characteristic of condensed fuel burning. Two different burners are used, a circular burner (62 mm diameter) placed inside a cylindrical chamber (0.3 m diameter and 1.0 m height) and a rectangular burner (50 mm wide by 200 mm long) placed in a wind tunnel (350 mm long) of rectangular cross section (120 mm wide and 90 mm height). The first burner is used to study the effect of pressure and gravity in the absence of a forced flow parallel to the surface. The second burner is used to study the effect of a forced flow parallel to the burner surface as well as the effect of oxygen concentration in the oxidizer flow. In this case experiments are also conducted at different gravity levels (micro-gravity, 0.2 g(sub 0), g(sub 0) and 1.8 g(sub 0)) to quantify the relative importance of buoyancy.

Torero, J. L.↗
NASA NTRS1983

Analysis of coolant entrance boundary shape of porous region to control cooling along exit boundary

A cooled porous region has a plane surface exposed to a specified spatially varying heat flux. The coolant leaves the region through this surface, and it is desired to control the flow distribution to maintain a specified uniform surface temperature. This is accomplished by having the coolant entrance surface shaped to provide in the region the necessary variation of path length and, hence, flow resistance. The surface shape at the coolant entrance is found by solving a Cauchy boundary value problem. An exact solution is obtained that will deal with a wide variety of heating distributions for both two- and three-dimensional shapes.

Siegel, R.↗
NASA NTRS1992

Crystallization, sublimation, and gas release in the interior of a porous comet nucleus

A numerical code is developed for evolutionary calculations of the thermal structure and composition of a porous comet nucleus made of water ice, in amorphous or crystalline form, other volatiles, dust, and gases trapped in amorphous ice. Bulk evaporation, crystallization, gas release, and free (Knudsen) flow of gases through the pores are taken into account. The numerical scheme yields exact conservation laws for mass and energy. The code is used to study the effect of bulk evaporation of ice in the interior of a comet nucleus during crystallization. It is found that evaporation controls the temperature distribution; the vapor prevents cooling of the crystallized layer of ice, by recondensation and release of latent heat. Thus high temperatures are maintained below the surface of the nucleus and down to depths of tens or hundreds of meters, even at large heliocentric distances, as long as crystallization goes on. Gas trapped in the ice and released during the phase transition flows both toward the interior and toward the surface and out of the nucleus. The progress of crystallization is largely determined by the contribution of gas fluxes to heat transfer.

Prialnik, Dina↗
NASA NTRS2009

Testing of Commercial Hollow Fiber Membranes for Space Suit Water Membrane Evaporator

Three commercial-off-the-shelf (COTS) hollow fiber (HoFi) membrane evaporators, modified for low pressure, were tested in a vacuum chamber at pressures below 33 pascals as potential space suit water membrane evaporator (SWME) heat rejection technologies. Water quality was controlled in a series of 25 tests, first simulating potable water reclaimed from waste water and then changing periodically to simulate the ever concentrating make-up of the circulating coolant over that is predicted over the course of 100 EVAs. Two of the systems, comprised of non-porous tubes with hydrophilic molecular channels as the water vapor transport mechanism, were severely impacted by the increasing concentrations of cations in the water. One of the systems, based on hydrophobic porous polypropylene tubes was not affected by the degrading water quality, or the presence of microbes. The polypropylene system, called SWME 1, was selected for further testing. An inverse flow configuration was also tested with SWME 1, with vacuum exposure on the inside of the tubes, provided only 20% of the performance of the standard configuration. SWME 1 was also modified to block 50% and 90% of the central tube layers, and tested to investigate performance efficiency. Performance curves were also developed in back-pressure regulation tests, and revealed important design considerations arising from the fully closed valve. SWME 1 was shown to be insensitive to air bubbles injected into the coolant loop. Development and testing of a full-scale prototype based on this technology and these test results is in progress.

Bue, Grant C.↗
NASA NTRS2019

Effect of Out-Gassing on the Onset of Hypersonic Boundary Layer Transition

Prediction and control of the onset of transition and the associated variation in aerothermodynamic parameters in high-speed flows is key to optimize the performance and design of Thermal Protection Systems (TPS) of next-generation aerospace vehicles [1]. Boundary Layer Transition (BLT) characteristics can influence the surface heating budget determining the TPS thickness and consequently its weight penalty. Ablative heatshields are designed to alleviate the high heat flux at the surface through pyrolysis of their polymeric matrix and subsequent fiber ablation [2]. Pyrolysis leads to out-gassing and non-uniform ablation lead to surface roughness, both of which are known to influence the transition process. An ablator impacts BLT through three main routes: gas injecting into the boundary layer from the wall, changing the surface heat transfer due to wall-flow chemical reactions, and modifying surface roughness [3]. In preparation to Mars 2020 mission post-flight analysis, the predictive transition capability has been initiated toward hard-coupling porous material response analysis and aerothermal environment calculation.

Karimi, Mona↗
DOE OSTI2023

Energy storage material and method of producing the same

The instant invention includes a spherical porous secondary silicon-based particle and methods for producing the same. The spherical porous secondary silicon-based particle is comprised of agglomerated primary silicon-based nanoparticles. The secondary particle comprises a carbon coating that reduces the effective exposed surface area of the primary particles to the electrolyte, thus improving first cycle efficiency. The secondary particle further comprises porous regions that enable the silicon nanoparticles to expand during lithiation. Advantages include ease of castability with micron-sized spherical particles, ease of mixing spherical particles, ease of flow for spherical particles in various processing steps, and ease with obtaining higher loading, which translates to higher areal capacity and overall energy density of the cell. A readily scalable process for producing the particles using low-cost materials and low-cost processing methods is disclosed.

Holt, Christopher T.↗
DOE OSTI2020

In Situ Chemical Monitoring and Imaging of Contents within Microfluidic Devices Having a Porous Membrane Wall Using Liquid Microjunction Surface Sampling Probe Mass Spectrometry

The ability to observe dynamic chemical processes (e.g., signaling, transport, etc.) in vivo or in situ using nondestructive chemical imaging opens a new door to understanding the complex dynamics of developing biological systems. With the advent of “biology-on-a-chip” devices has come the ability to monitor dynamic chemical processes in a controlled environment, using these engineered habitats to capture key features of natural systems while allowing visual observation of system development. Having the capability to spatially and temporally map the chemical signals within these devices may yield new insights into the forces that drive biosystem development. Here, a porous membrane sealed microfluidic device was designed to allow normal microfluidic operation while enabling continuous, location specific sampling and chemical characterization by liquid microjunction surface sampling probe mass spectrometry (LMJ-SSP MS). LMJ-SSP was used to extract fluids with nL-to-μL/min flow rates directly from selected areas of the microfluidic device without negatively impacting the device function. These extracts were subsequently characterized using MS. This technique was used to acquire MS images of the entirety of several multi-input microfluidic devices having different degrees of fluid mixing. LMJ-SSP MS imaging visualized the spatial distribution of chemical components within the microfluidic channels and could visualize chemical reactions occurring in the device. These microfluidic devices with a porous membrane wall are wholly compatible with the construction of biology-on-a-chip devices. This ultimately would enable correlation of biosystem physical structure with an evolving chemical environment

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Insights into the rapid two-phase transport dynamics in different structured porous transport layers of water electrolyzers through high-speed visualization

In proton exchange membrane electrolyzer cells (PEMECs), maintaining efficient two-phase transport is one of the most important functions of porous transport layers (PTLs). To enhance the two-phase transport in PTLs, thin/titanium liquid/gas diffusion layers (TT-LGDLs) are introduced in PEMECs, and their difference from the conventional Ti felt PTLs are analyzed in-situ through high-speed and microscale visualization and electrochemical characterizations. The visualization results show that unfavorable large slugs can be greatly reduced in the PEMEC with a TT-LGDL compared to the PEMEC with a Ti felt PTL. More importantly, the recovery capability of water starvation with different PTLs is studied. After water starvation, the PEMEC with the TT-LGDL can recover the water starvation much more rapidly than the Ti felt cell, benefiting from its short and straight-through flow paths. Furthermore, the TT-LGDL tends to generate oxygen bubbles that are almost six times smaller and 236 times more frequently than the Ti felt PTL, indicating significantly boosted removal efficiency of produced oxygen and PEMEC performance. Finally, this study offers new insights into the dynamic processes of two-phase transport and the recovery capability of water starvation for different PTLs, which will provide valuable guidance for further optimization of PTLs and performance enhancement of PEMECs.

08 HYDROGEN↗
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↗