SEARCH · Engineering Papers

Results for “complex fluids”

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 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Ground-Based Cryogenic Leak Test of Fittings for Cryogenic Fluid Management

Mechanically connected fluid joints are virtually unavoidable in complex cryogenic system designs. In the case of spacecraft, the performance of these joints is critical to mission success. This is especially true for long-duration space missions where even very low leak rates can eventually lead to significant propellant losses or failure of vital cryogenic cooling systems. NASA has recently undertaken an effort to quantify the leak rate of Vacuum Coupling Radiation (VCR) fittings over the temperature range from ambient to 20 K, both before and after exposure to a launch vibration profile. A test apparatus employing a cryocooler and a calibrated helium mass spectrometer was developed, validated, and used to obtain quantifiable leak rates at a fitting test pressure of 31 bar (450 psig). Three different fitting sizes were tested, 25.4 mm (¼ inch), 12.7 mm (½ inch), and 6.35 mm (1 inch) and two gasket materials, stainless steel, and silver-plated nickel. Each fitting configuration (size/seal material) was subjected to two consecutive cryogenic thermal cycles/measurement tests, followed by a launch vibration test profile at ambient temperature, and then two additional cryogenic thermal cycles/measurement tests. The design and development of the test apparatus, and test data are presented and discussed in detail.

G J Tamasy↗

NASA NTRS1985

Laser Doppler velocimetry primer

Advanced research in experimental fluid dynamics required a familiarity with sophisticated measurement techniques. In some cases, the development and application of new techniques is required for difficult measurements. Optical methods and in particular, the laser Doppler velocimeter (LDV) are now recognized as the most reliable means for performing measurements in complex turbulent flows. And such, the experimental fluid dynamicist should be familiar with the principles of operation of the method and the details associated with its application. Thus, the goals of this primer are to efficiently transmit the basic concepts of the LDV method to potential users and to provide references that describe the specific areas in greater detail.

Bachalo, William D.↗

DOE OSTI2020

Chemically controlled shape-morphing of elastic sheets

Two-dimensional responsive materials that change shape into complex three-dimensional structures are valuable for creating systems ranging from wearable electronics to soft robotics. Typically, the final 3D structure is unique and predetermined through the materials’ processing. Here, we use theory and simulation to devise a distinctive approach for driving shape changes of 2D elastic sheets in fluid-filled microchambers. The sheets are coated with catalyst to generate controllable fluid flows, which transform the sheets into complex 3D shapes. A given shape can be achieved by patterning the arrangement of the catalytic domains on the sheet and introducing the appropriate reactant to initiate a specific catalytic reaction. Moreover, a single sheet that encompasses multiple catalytic domains can be transformed into a variety of 3D shapes through the addition of one or more reactants. Materials systems that morph on-demand into a variety of distinct structures can simplify manufacturing processes and broaden the utility of soft materials.

36 MATERIALS SCIENCE↗

NASA NTRS2002

Stirring Up an Elastic Fluid: Critical Viscosity of Xenon-2 (CVX-2)

Whipped cream stays in place even when turned upside down. Yet it readily flows through the nozzle of a spray can to reach the dessert plate. This demonstrates the phenomenon of shear thinning that is important to many industrial and physical processes. Paints, film emulsions, and other complex solutions that are highly viscous under normal conditions but become thin and flow easily under shear forces. A simple fluid, such as water, does not exhibit shear thinning under normal conditions. Very close to the liquid-vapor critical point, where the distinction between liquid and vapor disappears, the fluid becomes more complex and is predicted to display shear thinning. At the critical point, xenon atoms interact over long distances in a classical model of cooperative phenomena. Physicists rely on this system to learn how long-range order arises. The Critical Viscosity of Xenon Experiment (CVX-2) will measure the viscous behavior of xenon, a heavy inert gas used in flash lamps and ion rocket engines, at its critical point. Although it does not easily combine with other chemicals, its viscosity at the critical point can be used as a model for a range of fluids. Viscosity originates from the interactions of individual molecules. It is so complicated that, except for the simplest gas, it cannot be calculated accurately from theory. Tests with critical fluids can provide key data, but are limited on Earth because critical fluids are highly compressed by gravity. CVX-2 employs a tiny metal screen vibrating between two electrodes in a bath of critical xenon. The vibrations and how they dampen are used to measure viscosity. CVX flew on STS-85 (1997), where it revealed that, close to the critical point, the xenon is partly elastic: it can 'stretch' as well as flow. For STS-107, the hardware has been enhanced to determine if critical xenon is a shear-thinning fluid.

Berg, Robert F.↗

DOE OSTI2023

Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales

An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Phenomenological Model of Nonlinear Dynamics and Deterministic Chaotic Gas Migration in Bentonite: Experimental Evidence and Diagnostic Parameters

Understanding gas migration in compacted clay materials, e.g., bentonite and claystone, is important for the design and performance assessment of an engineered barrier system of a radioactive waste repository system, as well as many practical applications. Existing field and laboratory data on gas migration processes in low-permeability clay materials demonstrate the complexity of flow and transport processes, including various types of instabilities, caused by nonlinear dynamics of coupled processes of liquid–gas exchange, dilation, fracturing, fracture healing, etc., which cannot be described by classical models of fluid dynamics in porous media. We here show that the complexity of gas migration processes can be explained using a phenomenological concept of nonlinear dynamics and deterministic chaos theory. To do so, we analyzed gas pressure and gas influx (i.e., input) and outflux (i.e., output), recorded during the gas injection experiment in the compact Mx80-D bentonite sample, and calculated a set of the diagnostic parameters of nonlinear dynamics and chaos, such a global embedding dimension, a correlation dimension, an information dimension, and a spectrum of Lyapunov exponents, as well as plotted 2D and 3D pseudo-phase-space strange attractors, based on the univariate influx and outflux time series data. These results indicate the presence of phenomena of low-dimensional deterministic chaotic behavior of gas migration in bentonite. In particular, during the onset of gas influx in the bentonite core, before the breakthrough, the development of gas flow pathways is characterized by the process of chaotic gas diffusion. After the breakthrough, with inlet-to-outlet movement of gas, the prevailing process is chaotic advection. During the final phase of the experiment, with no influx to the sample, the relaxation pattern of gas outflux is resumed back to a process of chaotic diffusion. The types of data analysis and a proposed phenomenological model can be used to establish the basic principles of experimental data-gathering, modeling predictions, and a research design.

36 MATERIALS SCIENCE↗

NASA NTRS1988

Gneiss-charnockite transformation at Kottavattam, Southern Kerala (India)

At Kottavattam, leucocratic granitic garnet-biotite gneisses (age less than 2 Ga) were partially transformed to coarse-grained charnockite along a system of conjugate fractures (N70E and N20W) and the foliation planes (N60 to 80W; dip 80 to 90 SW) about 550 m.y. ago. To examine and quantify changes in fabric, mineralogy, pore fluids and chemical composition associated with this process, large rock specimens showing gneiss-charnockite transition were studied in detail. The results of the present study corroborate the concept that charnockite formation at Kottavattam is an internally-generated phenomenon and was not triggered by the influx of carbonic fluids from a deep-seated source. It is suggested that charnockitization was caused by the following mechanism: (1) near-isothermal decompression during uplift of the gneiss complex led to an increase of the pore fluid pressure (P sub fluid greater than P sub lith) which - in a regime of anisotropic stress - triggered or at least promoted the development of conjugate fractures; (2) the simultaneous release of pore fluids from bursting fluid inclusions and their escape into the developing fracture system resulted in a drop of fluid pressure; and (3) the internal generation and buffering of the fluids and their, probably, limited migration in an entirely granitic rock system explains the absence of any significant metasomatic mass transfer.

Raith, M.↗

NASA NTRS1989

CFD for applications to aircraft aeroelasticity

Strong interactions of structures and fluids are common in many engineering environments. Such interactions can give rise to physically important phenomena such as those occurring for aircraft due to aeroelasticity. Aeroelasticity can significantly influence the safe performance of aircraft. At present exact methods are available for making aeroelastic computations when flows are in either the linear subsonic or supersonic range. However, for complex flows containing shock waves, vortices and flow separations, computational methods are still under development. Several phenomena that can be dangerous and limit the performance of an aircraft occur due to the interaction of these complex flows with flexible aircraft components such as wings. For example, aircraft with highly swept wings experience vortex induced aeroelastic oscillations. Correct understanding of these complex aeroelastic phenomena requires direct coupling of fluids and structural equations. Here, a summary is presented of the development of such coupled methods and applications to aeroelasticity since about 1978 to present. The successful use of the transonic small perturbation theory (TSP) coupled with structures is discussed. This served as a major stepping stone for the current stage of aeroelasticity using computational fluid dynamics. The need for the use of more exact Euler/Navier-Stokes (ENS) equations for aeroelastic problems is explained. The current development of unsteady aerodynamic and aeroelastic procedures based on the ENS equations are discussed. Aeroelastic results computed using both TSP and ENS equations are discussed.

Guruswamy, Guru P.↗

DOE OSTI2025

Spectrophotometric determination of the stability of La hydroxyl complexes at near neutral to alkaline pH from 25 to 75 °C

The hydrolysis of rare earth elements (REE) potentially controls their mobility during fluid-rock interaction in a broad range of pH and temperature conditions. However, there is still a lack of thermodynamic data for modeling accurately the stability of REE hydroxyl complexes in hydrothermal aqueous fluids. Here, in this study, UV–Vis spectrophotometric experiments were conducted from 25 to 75 °C in near-neutral to alkaline NaOH-bearing aqueous solutions with varying lanthanum (La) concentrations (0 to ∼0.23 mmol/kg). The color indicator m-cresol purple was used to determine in situ pH and derive the average OH− ligand number ($\overrightarrow{n}$) and formation constants for the La hydroxyl complexes (LaOH 2+ , La(OH) 2 + , and La(OH) 3 0 ). From 25 to 50 °C, $\overrightarrow{n}$ ranges between ∼1 and 2 at pH from 7.0 to 9.3. At 75 °C, $\overrightarrow{n}$ ranges between ∼1.5 and 3 at pH from 6.3 to 8.8. These results suggest the predominance of LaOH 2+ and La(OH) 2 + complexes from 25 to 50 °C, and an increased predominance of La(OH) 3 0 at 75 °C. The cumulative formation constants (β n °, n = 1 to 3) are derived for the reaction La 3+ + nOH − = La(OH) n 3-n , and fitted between 25 and 250 °C by combining the UV–Vis and literature solubility data. The resulting logβ n ° are expressed as function of temperature (T in Kelvin): logβ 1 ° = −1.786 + 0.0133 T + 1.049·10 3 /T; logβ 2 ° = −5.797 + 0.0267 T + 2.713·10 3 /T; logβ 3 ° = 6.435 + 0.0223 T + 512.7/T. A comparison between these new fits and existing extrapolations using the Helgeson-Kirkham-Flowers equation of state indicates significant differences in the predicted hydrolysis of La. The latter extrapolations should therefore be updated for the hydrolysis of REE.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Coupled heat pulse propagation in two-fluid plasmas

Because of the large mass differences between electrons and ions, the heat diffusion in electron-ion plasmas exhibits more complex behavior than simple heat diffusion found in typical gas mixtures. In particular, heat is diffused in two distinct, but coupled, channels. Conventional single fluid models neglect the resulting complexity, and can often inaccurately interpret the results of heat pulse experiments. However, by recognizing the sensitivity of the electron temperature evolution to the ion diffusivity, not only can previous experiments be interpreted correctly, but informative simultaneous measurements can be made of both ion and electron heat channels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Insight into the assembly of lipid-hyaluronan complexes in osteoarthritic conditions

Interactions between molecules in the synovial fluid and the cartilage surface may play a vital role in the formation of adsorbed films that contribute to the low friction of cartilage boundary lubrication. Osteoarthritis (OA) is the most common degenerative joint disease. Previous studies have shown that in OA-diseased joints, hyaluronan (HA) not only breaks down resulting in a much lower molecular weight (MW), but also its concentration is reduced ten times. Here, we have investigated the structural changes of lipid-HA complexes as a function of HA concentration and MW to simulate the physiologically relevant conditions that exist in healthy and diseased joints. Small angle neutron scattering and dynamic light scattering were used to determine the structure of HA-lipid vesicles in bulk solution, while a combination of atomic force microscopy and quartz crystal microbalance was applied to study their assembly on a gold surface. We infer a significant influence of both MW and HA concentrations on the structure of HA-lipid complexes in bulk and assembled on a gold surface. Our results suggest that low MW HA cannot form an amorphous layer on the gold surface, which is expected to negatively impact the mechanical integrity and longevity of the boundary layer and could contribute to the increased wear of the cartilage that has been reported in joints diseased with OA.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Origins of Clustering of Metalate–Extractant Complexes in Liquid–Liquid Extraction

Effective and energy-efficient separation of precious and rare metals is very important for a variety of advanced technologies. Liquid-liquid extraction (LLE) is a relatively less energy intensive separation technique, widely used in separation of lanthanides, actinides, and platinum group metals (PGMs). In LLE, the distribution of an ion between an aqueous phase and an organic phase is determined by enthalpic (coordination interactions) and entropic (fluid reorganization) contributions. The molecular scale details of these contributions are not well understood. Preferential extraction of an ion from the aqueous phase is usually correlated with the resulting fluid organization in the organic phase, as the longer-range organization increases with metal loading. However, it is difficult to determine the extent to which organic phase fluid organization causes, or is caused by, metal loading. In this study, we demonstrate that two systems with the same metal loading may impart very different organic phase organizations and investigate the underlying molecular scale mechanism. Small-angle X-ray scattering shows that the structure of a quaternary ammonium extractant solution in toluene is affected differently by the extraction of two metalates (octahedral PtCl 6 2- and square-planar PtCl 6 2- ), although both are completely transferred into the organic phase. The aggregates formed by the metalate-extractant complexes (approximated as reverse micelles) exhibit a more long-range order (clustering) with PtCl 6 2- compared to that with PtCl 6 2- . Vibrational sum frequency generation spectroscopy and complementary atomistic molecular dynamics simulations on model Langmuir monolayers indicate that the two metalates affect the interfacial hydration structures differently. Furthermore, the interfacial hydration is correlated with water extraction into the organic phase. Overall, these results support a strong relationship between the organic phase organizational structure and the different local hydration present within the aggregates of metalate-extractant complexes, which is independent of metalate concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

NASA NTRS2004

"Smart" Magnetic Fluids Experiment Operated on the International Space Station

InSPACE is a microgravity fluid physics experiment that was operated on the International Space Station (ISS) in the Microgravity Science Glovebox from late March 2003 through early July 2003. (InSPACE is an acronym for Investigating the Structure of Paramagnetic Aggregates From Colloidal Emulsions.) The purpose of the experiment is to obtain fundamental data of the complex properties of an exciting class of smart materials termed magnetorheological (MR) fluids. MR fluids are suspensions, or colloids, comprised of small (micrometer-sized) superparamagnetic particles in a nonmagnetic medium. Colloids are suspensions of very small particles suspended in a liquid. (Examples of other colloids are blood, milk, and paint.) These controllable fluids can quickly transition into a nearly solid state when exposed to a magnetic field and return to their original liquid state when the magnetic field is removed. Controlling the strength of the magnetic field can control the relative stiffness of these fluids. MR fluids can be used to improve or develop new seat suspensions, robotics, clutches, airplane landing gear, and vibration damping systems. The principal investigator for InSPACE is Professor Alice P. Gast of the Massachusetts Institute of Technology (MIT). The InSPACE hardware was developed at the NASA Glenn Research Center. The InSPACE samples were delivered to the ISS in November 2002, on the Space Shuttle Endeavour, on Space Station Utilization Flight UF-2/STS113. Operations began on March 31, 2003, with the processing of three different particle size samples at multiple test parameters. This investigation focused on determining the structural organization of MR colloidal aggregates when exposed to a pulsing magnetic field. On Earth, the aggregates take the shape of footballs with spiky tips. This characteristic shape may be influenced by the pull of gravity, which causes most particles initially suspended in the fluid to sediment, (i.e., settle and collect at the bottom of the cell). In the absence of sedimentation effects on the ISS, the behavior and shape of these MR aggregate structures are dominated exclusively by magnetic and surface tension forces. The microscopic detail of these structures was imaged under two orthogonal camera views. The video was downlinked to the InSPACE team at Glenn's Telescience Support Center and to MIT and also recorded onboard the ISS on videotapes that will be brought back to the ground by the space shuttles. The study examined the effect on the structure formation by varying the magnetic field strength and pulse frequency, and particle size. Fundamental data that characterized the structure formation were obtained. InSPACE completed its last planned test run on July 2, 2003. Operations occurred on 21 days over approximately a 3-month period. Forty-one test points were completed during 26 test runs. During the initial testing, the procedures followed by the crew were modified to maximize the observation of some unexpected and interesting aggregate behavior. As a result Dr. Gast has reported on the formation of aggregate shapes that are more extended and diverse than those observed on the ground. Sheets of magnetic material folded over in a labyrinth pattern and large columnar aggregates with complex interfaces with the surrounding fluid are examples of the interesting structures that have been observed on the ISS. In light of these early findings, the understanding of the fundamental properties of MR fluids on the basis of ground-based observations may need to be reconsidered.The experiments on the ISS have provided a vast amount of video data for analysis. While this analysis is ongoing, plans are being made for additional experimental runs. For this purpose, additional hardware and cells containing samples of different magnetic particles and sizes are being fabricated for a future launch to the ISS. The InSPACE hardware will remain on orbit until this testing is completed.

Agui, Juan H.↗

DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗

DOE OSTI2025

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗

DOE OSTI2020

Experimental and Numerical Investigation of a Miniature Additively Manufactured Vortex Tube

Ranque–Hilsch vortex tubes are simple devices that can produce a cooling effect using compressed air. A key advantage of vortex tubes is the lack of moving solid parts; however, their efficiencies are relatively low. The present study focuses on the development of a miniature variable-diameter tube using additive manufacturing. A metal-based 3D printing technique was utilized to fabricate this vortex tube monolithically. Computational fluid dynamics simulations employing software STAR-CCM + with a compressible Reynolds-Averaged Navier–Stokes (RANS) approach and the elliptic-blending lag k-epsilon turbulence model have been applied to model thermofluid processes inside the vortex tube, to good agreement with the experiment. In conclusion, a temperature decrease of 13.3 °C and a cooling power of approximately 4 W were experimentally achieved with a pressure ratio of 4 in the air at normal conditions. This result shows promise for the goal of utilizing additive manufacturing to design and build complex-geometry vortex tubes intended for use with cryogenic fluids.

08 HYDROGEN↗

DOE OSTI2024

Multiphysics Demonstration of Temperature-Driven Assembly Bowing in SFRs using MOOSE-Based Codes

Core bowing is an important passive safety mechanism in liquid metal cooled fast reactors. When the core restraint system is properly designed, temperature and flux gradients influence assemblies in the core to bow into less reactive configurations during accident scenarios, resulting in negative reactivity feedback. Prediction of core bowing involves complex interplay of radiation transport, impacts of fluid flow and heat transfer on duct temperature, and mechanical responses to the induced temperature and flux gradients. Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program [1], an integrated multiphysics approach is being developed to model the core bowing phenomena in liquid metal-cooled fast reactors with the Multiphysics Object Oriented Simulation Environment (MOOSE) [2]. In this methodology, the MOOSE-based reactor physics code Griffin [3] will solve the neutron transport equation and determine the power distribution. With the detailed power distribution from Griffin, the subchannel analysis codes MOOSE-Subchannel [4] and Pronghorn [5] are utilized to calculate the assembly temperature distribution. MOOSE’s Solid Mechanics [6] and Contact [7] Modules are leveraged to calculate the thermal expansion and duct bowing displacement with the duct wall temperature from thermal hydraulics calculation. In this work, an initial one-way coupling demonstration of the integrated multiphysics approach has been performed on a seven-assembly problem based on the sodium-cooled fast reactor ABR-1000 design [8]. The neutronics calculation with Griffin is not yet involved in the current simulation. MOOSE-Subchannel and Pronghorn evaluate fluid and solid temperature based on a fixed power distribution. In addition, one-way coupling is utilized in this coupled calculation, via Pronghorn passing the duct temperature data to the MOOSE Solid Mechanics calculation. An assessment of the Solid Mechanics module was performed in parallel to verify duct bowing behavior with duct-to-duct contact phenomenon [9]. The displacement from MOOSE Solid Mechanics is not yet transferred back and utilized in the Pronghorn and MOOSE-Subchannel calculation. This model will be available on the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB) repository [10]. Future stages of this work will involve solving problems of increasing complexity as well as adding more physics (e.g. reactor physics) to the integrated workflow to reach the end goal of modeling the core bowing phenomenon with an integrated multiphysics workflow.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗