Search NASA⌕ Search

SEARCH · Search NASA

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Energy spectra of finite temperature superfluid helium-4 turbulence

A mesoscopic model of finite temperature superfluid helium-4 based on coupled Langevin-Navier-Stokes dynamics is proposed. Drawing upon scaling arguments and available numerical results, a numerical method for designing well resolved, mesoscopic calculations of finite temperature superfluid turbulence is developed. The application of model and numerical method to the problem of fully developed turbulence decay in helium II, indicates that the spectral structure of normal-fluid and superfluid turbulence is significantly more complex than that of turbulence in simple-fluids. Analysis based on a forced flow of helium-4 at 1.3 K, where viscous dissipation in the normal-fluid is compensated by the Lundgren force, indicate three scaling regimes in the normal-fluid, that include the inertial, low wavenumber, Kolmogorov k{sup −5/3} regime, a sub-turbulence, low Reynolds number, fluctuating k{sup −2.2} regime, and an intermediate, viscous k{sup −6} range that connects the two. The k{sup −2.2} regime is due to normal-fluid forcing by superfluid vortices at high wavenumbers. There are also three scaling regimes in the superfluid, that include a k{sup −3} range that corresponds to the growth of superfluid vortex instabilities due to mutual-friction action, and an adjacent, low wavenumber, k{sup −5/3} regime that emerges during the termination of this growth, as superfluid vortices agglomerate between intense normal-fluid vorticity regions, and weakly polarized bundles are formed. There is also evidence of a high wavenumber k{sup −1} range that corresponds to the probing of individual-vortex velocity fields. The Kelvin waves cascade (the main dynamical effect in zero temperature superfluids) appears to be damped at the intervortex space scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

4D Proxy Imaging of Fracture Dilation and Stress Shadowing Using Electrical Resistivity Tomography During High Pressure Injections into a Dense Rock Formation

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

electrical resisitivity tomography, stress, 4D Tim↗
DOE OSTI2022

Data-driven, multi-moment fluid modeling of Landau damping

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. Here, the learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effect such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve the fluid closure and reduce the computational cost of multi-scale modeling of global systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Liquid-based gating mechanism with tunable multiphase selectivity and antifouling behaviour

Living organisms make extensive use of micro- and nanometre-sized pores as gatekeepers for controlling the movement of fluids, vapours and solids between complex environments. The ability of such pores to coordinate multiphase transport, in a highly selective and subtly triggered fashion and without clogging, has inspired interest in synthetic gated pores for applications ranging from fluid processing to 3D printing and lab-on-chip systems(1-10). But although specific gating and transport behaviours have been realized by precisely tailoring pore surface chemistries and pore geometries(6,11-17), a single system capable of controlling complex, selective multiphase transport has remained a distant prospect, and fouling is nearly inevitable(11,12). Here we introduce a gating mechanism that uses a capillary-stabilized liquid as a reversible, reconfigurable gate that fills and seals pores in the closed state, and creates a non-fouling, liquid-lined pore in the open state. Theoretical modelling and experiments demonstrate that for each transport substance, the gating threshold-the pressure needed to open the pores-can be rationally tuned over a wide pressure range. This enables us to realize in one system differential response profiles for a variety of liquids and gases, even letting liquids flow through the pore while preventing gas from escaping. These capabilities allow us to dynamically modulate gas-liquid sorting in a microfluidic flow and to separate a three-phase air-water-oil mixture, with the liquid lining ensuring sustained antifouling behaviour. Because the liquid gating strategy enables efficient long-term operation and can be applied to a variety of pore structures and membrane materials, and to micro- as well as macroscale fluid systems, we expect it to prove useful in a wide range of applications.

Hou, X↗
NASA NTRS1988

The role of boron and fluids in high temperature, shallow level metamorphism of the Chugach Metamorphic Complex, Alaska

The possible role of boron (B) involvement in granite equilibria and generation of melts during crustal metamorphism has been a focus of speculation in recent literature. Most of the evidence for such involvement derives from experimental data which implies that the addition of B will lower the temperature of the granite solidus. Also the presence of tourmaline has a minor effect on the temperature of the solidus. Further indirect evidence that B may be involved in partial melting processes is the observation that granulites are commonly depleted in B, whereas the B content of low grade metapelites can be high (up to 2000 ppm). Researchers' measurements of the whole-rock B contents of granulites from the Madras region, India are low, ranging from 0.4 to 2.6 ppm. Ahmad and Wilson suggest that B was mobilized in the fluid phase during granulite facies metamorphism of the Broken Hill Complex, Australia. Thus, it appears that during the amphibolite to granulite transition, B is systematically lost from metasediments. The B that is released will probably partition into the vapor phase and/or melt phase. Preliminary measurements imply that the boron content of rocks in the Chugach Metamorphic Complex is not sufficient to influence the processes of partial melting at low pressures.

Sisson, V. B.↗
NASA NTRS2005

Investigation of Compressibility Effect for Aeropropulsive Shear Flows

Rocket Based Combined Cycle (RBCC) engines operate within a wide range of Mach numbers and altitudes. Fundamental fluid dynamic mechanisms involve complex choking, mass entrainment, stream mixing and wall interactions. The Propulsion Research Center at the University of Alabama in Huntsville is involved in an on- going experimental and numerical modeling study of non-axisymmetric ejector-based combined cycle propulsion systems. This paper attempts to address the modeling issues related to mixing, shear layer/wall interaction in a supersonic Strutjet/ejector flow field. Reynolds Averaged Navier-Stokes (RANS) solutions incorporating turbulence models are sought and compared to experimental measurements to characterize detailed flow dynamics. The effect of compressibility on fluids mixing and wall interactions were investigated using an existing CFD methodology. The compressibility correction to conventional incompressible two- equation models is found to be necessary for the supersonic mixing aspect of the ejector flows based on 2-D simulation results. 3-D strut-base flows involving flow separations were also investigated.

Balasubramanyam, M. S.↗
DOE OSTI2024

Solving the Hele–Shaw flow using the Harrow–Hassidim–Lloyd algorithm on superconducting devices: A study of efficiency and challenges

The development of quantum processors for practical fluid flow problems is a promising yet distant goal. Recent advances in quantum linear solvers have highlighted their potential for classical fluid dynamics. In this study, we evaluate the Harrow–Hassidim–Lloyd (HHL) quantum linear systems algorithm (QLSA) for solving the idealized Hele–Shaw flow. Our focus is on the accuracy and computational cost of the HHL solver, which we find to be sensitive to the condition number, scaling exponentially with problem size. This emphasizes the need for preconditioning to enhance the practical use of QLSAs in fluid flow applications. Moreover, we perform shots-based simulations on quantum simulators and test the HHL solver on superconducting quantum devices, where noise, large circuit depths, and gate errors limit performance. Error suppression and mitigation techniques improve accuracy, suggesting that such fluid flow problems can benchmark noise mitigation efforts. Finally, our findings provide a foundation for future, more complex application of QLSAs in fluid flow simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A radioisotope - enabled reactive transport model for deep vadose zone carbon

In mountainous regions, which constitute the principle source of recharge to major rivers and regional aquifers, infiltration occurs through fractured, partially saturated, weathered bedrock that acts as a boundary layer between saturated aquifers and surface soil. Commonly this deep vadose zone (DVZ) is many meters thick, and yet its role in regulating the generation, retention and mobility of reactive solutes, including nutrients, contaminants and weathering products, is largely unknown. In particular, many of the key reactions that drive the formation of the weathered DVZ and the quality of water moving through it are redox processes, regulated by the availability of organic carbon and oxygen below the soil layer. The role of the DVZ is thus also poorly constrained in the context of carbon stocks and mobility, particularly in lithologies that are naturally high in organic carbon, such as shales. The overarching hypothesis of this study is that upland regions developed in geologic settings with abundant petrogenic carbon store and actively cycle carbon in the weathered DVZ below the soil and above the water table at rates that are significant and currently unconstrained. In order to quantify this cycling, the current study combines novel instrumentation techniques allowing new direct sampling of DVZ systems with advanced numerical reactive transport simulations of carbon transport and transformation. Critically, these simulations will explicitly treat the three isotopes of carbon (the abundant 12C, the stable rare 13C and the radioactive 14C) in a unified framework, thus clearly parsing between the contributions of modern surface derived carbon and lithologic carbon sources in integrated measurements of fluid and gas phase fluxes. This novel model capability will be applied to test the role of DVZ carbon cycling as a regulator of water quality and geological weathering in two complementary field sites both located in organic carbon rich shale lithologies. The first is the Eel River Critical Zone Observatory (ERCZO) in Mendocino County, California, and the second is the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area (SFA) in the East River watershed, near Crested Butte, Colorado. At the ERCZO, a novel vadose zone monitoring system has been installed in a 20 m thick, partially saturated, weathered shale hillslope, and preliminary data already indicate substantial CO2 flux generated many meters below the soil surface. At the SFA field site, an instrumented hillslope transect indicates a more complex multi-dimensional fluid and solute transport regime, which will serve as a key test of the calibrated models. Collectively, this project will advance understanding of the cycling of carbon belowground and in relation to transport pathways across the poorly constrained DVZ characteristic of primary water recharge areas. The key product of this work will be enhanced isotope simulation capabilities that are robust and publicly available for application across a broad diversity of systems.

58 GEOSCIENCES↗
NASA NTRS1988

Modeling of supersonic reacting flow fields

A detailed understanding of the scramjet combustor flow field is critical to the achievement of a successful design. Even though the combustor flow field is quite complex, it can be realistically viewed as a collection of spatially developing and reacting supersonic mixing layers that are initially discrete, but that ultimately merge into larger more complex zones. These mixing layers begin downstream of a set of fuel injectors that introduce gaseous hydrogen in both a parallel and transverse direction into a supersonic air stream entering from the engine inlet. The behavior of the initial portion of the combustor flow, in the mixing layers near the fuel injectors, appears to be most critical, since this is where the mechanism for efficient high speed mixing must be established to achieve the required degree of combustion downstream. Because of the structure of the flow field in this initial portion of the combustor, a single supersonic, spatially developing and reacting mixing layer serves as an excellent physical model for the overall flow field. Even though this reacting mixing layer flow is geometrically simple, it can still be made to retain all of the fluid mechanical and chemical complexities present in the actual combustor flow field.

Drummond, J. P.↗
NASA NTRS1995

A CFD study of complex missile and store configurations in relative motion

An investigation was conducted from May 16, 1990 to August 31, 1994 on the development of computational fluid dynamics (CFD) methodologies for complex missiles and the store separation problem. These flowfields involved multiple-component configurations, where at least one of the objects was engaged in relative motion. The two most important issues that had to be addressed were: (1) the unsteadiness of the flowfields (time-accurate and efficient CFD algorithms for the unsteady equations), and (2) the generation of grid systems which would permit multiple and moving bodies in the computational domain (dynamic domain decomposition). The study produced two competing and promising methodologies, and their proof-of-concept cases, which have been reported in the open literature: (1) Unsteady solutions on dynamic, overlapped grids, which may also be perceived as moving, locally-structured grids, and (2) Unsteady solutions on dynamic, unstructured grids.

Baysal, Oktay↗
NASA NTRS1988

Three dimensional grid generation for complex configurations: Recent progress

Some of the capabilities of the Computational Fluid Dynamics (CFD) community are surveyed for griding complex three dimensional configurations. The intent is to provide some insight as to the present state of grid generation for aircraft configurations in order to help assess whether this task presents a long term stumbling block to routine use of CFD in aerodynamic applications. A brief review is given of some of the techniques that are available for generating body conforming curvilinear grids. In order to assess capabilities in grid generation, colleagues at selected institutions were solicited to describe their experiences and difficulties in grid generation of complex configurations.

J.F. Thompson↗
DOE OSTI2024

A Review of Coupled Geochemical–Geomechanical Impacts in Subsurface CO 2 , H 2 , and Air Storage Systems

Increased demand for decarbonization and renewable energy has led to increasing interest in engineered subsurface storage systems for large-scale carbon reduction and energy storage. In these applications, a working fluid (CO 2 , H 2 , air, etc.) is injected into a deep formation for permanent sequestration or seasonal energy storage. The heterogeneous nature of the porous formation and the fluid–rock interactions introduce complexity and uncertainty in the fate of the injected component and host formations in these applications. Interactions between the working gas, native brine, and formation mineralogy must be adequately assessed to evaluate the efficiency, risk, and viability of a particular storage site and operational regime. This study reviews the current state of knowledge about coupled geochemical–geomechanical impacts in geologic carbon sequestration (GCS), underground hydrogen storage (UHS), and compressed air energy storage (CAES) systems involving the injection of CO 2 , H 2 , and air. Specific review topics include (1) existing injection induced geochemical reactions in these systems; (2) the impact of these reactions on the porosity and permeability of host formation; (3) the impact of these reactions on the mechanical properties of host formation; and (4) the investigation of geochemical-geomechanical process in pilot scale GCS. This study helps to facilitate an understanding of the potential geochemical–geomechanical risks involved in different subsurface energy storage systems and highlights future research needs.

08 HYDROGEN↗
DOE OSTI2020

Analysis of combustion acoustic phenomena in compression–ignition engines using large eddy simulation

As computational capabilities continue to grow, exploring the limits of computational fluid dynamics to capture complex and elusive phenomena, which are otherwise difficult to study by experimental techniques, is one of the main targets for the research community. This paper presents a detailed analysis of the physical processes that lead to combustion noise emissions in internal combustion engines. In particular, diesel combustion in a compression-ignition (CI) engine is studied in order to understand the singular behavior of the in-cylinder flow field responsible for the acoustic emissions. The main objective is, therefore, to improve the understanding of the phenomena involved in CI engine noise using large eddy simulations. Several visualization methods are employed to investigate the connection between combustion behavior and its effects on the pressure field. In addition, proper orthogonal decomposition is used to analyze the modal energy distribution among all the acoustic modes. The results show that the acoustic signature is fundamentally conditioned by the intensity of the premixed combustion rather than by the pressure oscillations generated by turbulent fluctuations in the flame surface established during the diffusion stage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗
DOE OSTI2026

Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media

Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.

CT scanning↗
DOE OSTI2019

Highly resolved Eulerian simulations of fuel spray transients in single and multi-hole injectors: Nozzle flow and near-exit dynamics

In high pressure fuel injectors, needle opening and closing transients cause complex off-design fluid dynamics behaviors that profoundly impact the spray and mixture formation processes. These dynamics are completely different from what is known to occur in steady state conditions. In this study, diesel spray transients have been investigated in single-hole and 3-hole nozzles, encompassing internal and external nozzle flow and including needle motion, performing highly resolved (2.5 μm) computational fluid dynamics (CFD) simulations. We focused on end-of-injection (EOI) and start-of-injection (SOI) processes, in order to provide insights in to the physics. The liquid fuel, vapor and gas species are modeled with a single-fluid multiphase mixture approach, with diffuse interface, and with large eddy simulations (LES) of the turbulence. Occurrence of phase change due to cavitation is accounted for, and the spray dispersion is described with a turbulent dispersion model. Detailed needle motion data and orifice internal surface are available from xray synchrotron source measurements carried out at Argonne National Laboratory, and shared through the Engine Combustion Network (ECN) community. Simulations are compared against x-ray phase contrast imaging and radiography of the internal and near-exit flow, in addition to optical microscopy data of the near-exit sprays. Simulation results are found to agree well with available experimental data, and are able to realistically capture local and global features. The simulations allow to gain insight into the physics of gas ingestion and dribbles at EOI, for different hole diameters, operating conditions and number of holes. At SOI, timing of liquid appearance out of the injector and spray tip penetration are adequately predicted, by using the EOI flow field as in-nozzle initialization, and by prescribing the measured tip needle displacement with an informed effective valve opening point inferred from the x-ray observations. Lastly, the variation of spreading angle over time is also discussed in detail for the multi-hole case, including hole-to-hole variations. Due to real geometry features and asymmetric needle motion with eccentric components, it is found that the three holes exhibit swirling flows of increasing intensity as the lift decreases, causing the near cone angle to open and spread, in a quasi-hollow cone structure. These features are not observed in axial single-hole injectors because of their relative simplicity and intrinsic symmetry.

42 ENGINEERING↗
DOE OSTI2018

Microseismicity cloud can be substantially larger than the associated stimulated fracture volume: the case of the Paralana Enhanced Geothermal System

Here, the goal of hydraulic stimulation is to increase formation permeability in the near vicinity of a well. However, there remain technical challenges around measuring the outcome of this operation. During two Enhanced Geothermal System (EGS) stimulations in South Australia, Paralana in 2011 and Habanero in 2003, extensive catalogs of microseismicity were recovered. It is often assumed that shear failure of existing fractures is the main mechanism behind both the induced earthquakes and any permeability enhancement. This underpins a common notion, that the seismically active volume is also the stimulated reservoir. In this paper, we compute the density of earthquake hypocenters and provide evidence that, under certain conditions, this spatiotemporal quantity is a reasonable proxy for pore pressure increase. We then apply an inverse modeling approach that uses the earthquake observations and a modified reservoir simulator to estimate the parameters of a permeability evolution relation. The regime implied by the data indicates that most permeability enhancement occurred very near to the wellbore and was not coincident with the bulk of the seismicity, whose volume was about two orders of magnitude larger. Thus, we conclude that, in some cases, it is possible for permeability enhancement and induced seismicity to be decoupled, in which case the seismically active volume is a poor indicator of the stimulated reservoir. Our results raise serious questions about the effectiveness of hydroshearing as a stimulation mechanism in EGS. Finally, this study extends our understanding of the complex processes linking earthquakes, fluid pressure, and permeability in both natural and engineered settings.

58 GEOSCIENCES↗
DOE OSTI2016

Rapid expulsion of microswimmers by a vortical flow

Interactions of microswimmers with their fluid environment are exceptionally complex. Macroscopic shear flow alters swimming trajectories in a highly nontrivial way and results in dramatic reduction of viscosity and heterogeneous bacterial distributions. Here we report on experimental and theoretical studies of rapid expulsion of microswimmers, such as motile bacteria, by a vortical flow created by a rotating microparticle. We observe a formation of a macroscopic depletion area in a high-shear region, in the vicinity of a microparticle. The rapid migration of bacteria from the shear-rich area is caused by a vortical structure of the flow rather than intrinsic random fluctuations of bacteria orientations, in stark contrast to planar shear flow. Our mathematical model reveals that expulsion is a combined effect of motility and alignment by a vortical flow. Our findings offer a novel approach for manipulation of motile microorganisms and shed light on bacteria-flow interactions.

36 MATERIALS SCIENCE↗