Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex fluid”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
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 OSTI2019

Investigative Studies toward the Extraction and Recovery of Praseodymium From Ionic Liquids Using Supercritical Carbon Dioxide

Recycle of rare earth elements requires the extraction and recovery of these metals from consumer products such as rechargeable batteries, electric motor magnets, and aircraft engines once the rare earths have been oxidatively dissolved. One approach to this is to combine the metal extraction and recovery steps into a continuous process. Such a scheme would involve direct-dissolution and ligand-assisted supercritical fluid extraction of the metal, followed by deposition of the extracted metal complex into an ionic liquid where it can be electrochemically reduced to its metallic state. Whereas supercritical fluid extraction of rare earth metals has been recently studied, deposition of rare earth metal-ligand complexes into an ionic liquid is not well known for many rare earth metal-ligand systems. In these studies we measured the solubility of several Praseodymium complexes in scCO2. We also measured the solubility of the Praseodymium complexes in ionic liquids with different physical properties of interest including polarizability, viscosity, and conductivity. The ionic liquids studied are from the pyrrolidinium family and include n-butyl-n-methylpyrrolidinium bis(trifluoromethylsulfonyl) imide (BMPy TFSI). Upon completion of these experiments, we added an additional modifier to the ionic liquid solution in the form of a chelating agent and carbon dioxide. Overall, the additions of the various chelating agents and BMPy TFSI result in large changes to viscosity, conductivity, and metal complexation properties of the solutions. Here we report the results of our findings and discuss the possibility of tuning the extractive and reductive properties of the solutions by picking ionic liquids and chelating agents with physical properties conducive to the desired metal recovery scheme.

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 OSTI2014

Code Coupling via Jacobian-Free Newton-Krylov Algorithms with Application to Magnetized Fluid Plasma and Kinetic Neutral Models

Jacobian-free Newton-Krylov (JFNK) algorithms are a potentially powerful class of methods for solving the problem of coupling codes that address dfferent physics models. As communication capability between individual submodules varies, different choices of coupling algorithms are required. The more communication that is available, the more possible it becomes to exploit the simple sparsity pattern of the Jacobian, albeit of a large system. The less communication that is available, the more dense the Jacobian matrices become and new types of preconditioners must be sought to efficiently take large time steps. In general, methods that use constrained or reduced subsystems can offer a compromise in complexity. The specific problem of coupling a fluid plasma code to a kinetic neutrals code is discussed as an example.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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.↗
DOE OSTI2016

The Nature and Origin of the REE Mineralization in the Wicheeda Carbonatite, British Columbia, Canada

The Wicheeda carbonatite is a deformed plug or sill that hosts relatively high grade light rare earth elements (LREE) mineralization in the British Columbia alkaline province. It was emplaced within metasedimentary rocks belonging to the Kechika Group, which have been altered to potassic fenite near the intrusion and sodic fenite at greater distances from it. The intrusion comprises a ferroan dolomite carbonatite core, which passes gradationally outward into calcite carbonatite. The potentially economic REE mineralization is hosted by the dolomite carbonatite. We recognized three types of dolomite. Dolomite constitutes the bulk of the dolomite carbonatite, dolomite replaced dolomite near veins and vugs, and dolomite occurs in veins and vugs together with the REE mineralization. Carbon and oxygen isotope ratios indicate that the calcite carbonatite crystallized from a magma of mantle origin, that dolomite is of primary igneous origin, that dolomite has a largely igneous signature with a small hydrothermal component, and that dolomite is of hydrothermal origin. Furthermore, the REE minerals comprise REE fluorocarbonates, ancylite-(Ce), and monazite-(Ce). In addition to dolomite, they occur with barite, molybdenite, pyrite, and thorite. Minor concentrations of niobium are present as magmatic pyrochlore in the calcite carbonatite. model is proposed in which crystallization of calcite carbonatite preceded that of dolomite carbonatite. During crystallization of the latter, an aqueous-carbonic fluid was exsolved, which mobilized the REE as chloride complexes into vugs and fractures in the dolomite carbonatite, where they precipitated mainly in response to the increase in pH that accompanied fluid-rock interaction and, in the case of the REE fluorocarbonates, decreasing temperature. These fluids altered the host metasedimentary rock to potassic fenite adjacent to the carbonatite and, distal to it, they mixed with formational waters to produce sodic fenite.

58 GEOSCIENCES↗
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 OSTI2016

Directing the phase behavior of polyelectrolyte complexes using chiral patterned peptides

Polyelectrolyte complexes (PECs) have a broad range of promising applications as soft materials due to their self-assembly and diversity of structure and chemical composition. Peptide polymer PECs are highly biocompatible and biodegradable, making them particularly useful for encapsulation of food additives and flavors, micellar drug delivery, medical and underwater adhesives, fetal membrane patches, and scaffolds for cell growth in tissue engineering. While parameters affecting PEC formation and stability in regards to charge effects are well researched, little is known about the effects of van der Waals interactions, hydrogen bonding, and secondary structure in these materials. Peptide chirality provides a unique opportunity to manipulate PEC phase to modulate the amount of solid-like (precipitate) or liquid-like (coacervate) character by influencing hydrogen bonding interactions among peptide chains. In previous work, we showed that chiral peptides form solid complexes, while complexes with even one racemic peptide were fluid. This raised the interesting question of how long a homochiral sequence must be to result in solid phase formation. In this work, we designed chiral patterned peptides of polyglutamic acid and polylysine ranging from 50 to 90% L-chiral residues with increasing numbers of sequential L-chiral residues before a chirality change. These polymers were mixed together to form PECs. We observed that 8 or more sequential L-chiral residues are necessary to achieve both the appearance of a precipitate phase and sustained beta-sheets in the complex, as determined by optical imaging and FTIR Spectroscopy. Less homochiral content results in formation of a coacervate phase. Thus, we show that chiral sequence can be used to control the phase transition of PECs. Understanding how to manipulate PEC phase using chiral sequence as presented here may enable tuning of the material properties to achieve the desired mechanical strength for coatings and polymer brushes, or the most effective molecular release kinetics for drug delivery applications, for example.

Pacalin, Naomi M.↗
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↗
NASA NTRS1986

The method of complex characteristics for design of transonic blade sections

A variety of computational methods were developed to obtain shockless or near shockless flow past two-dimensional airfoils. The approach used was the method of complex characteristics, which determines smooth solutions to the transonic flow equations based on an input speed distribution. General results from fluid mechanics are presented. An account of the method of complex characteristics is given including a description of the particular spaces and coordinates, conformal transformations, and numerical procedures that are used. The operation of the computer program COMPRES is presented along with examples of blade sections designed with the code. A user manual is included with a glossary to provide additional information which may be helpful. The computer program in Fortran, including numerous comment cards is listed.

Bledsoe, M. R.↗