Search NASA⌕ Search

SEARCH · Search NASA

Results for “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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A cut-cell finite volume – finite element coupling approach for fluid–structure interaction in compressible flow

We present a loosely coupled approach for the solution of fluid–structure interaction problems between a compressible flow and a deformable structure. The method is based on staggered Dirichlet–Neumann partitioning. The interface motion in the Eulerian frame is accounted for by a conservative cut-cell Immersed Boundary method. The present approach enables sub-cell resolution by considering individual cut-elements within a single fluid cell, which guarantees an accurate representation of the time-varying solid interface. The cut-cell procedure inevitably leads to non-matching interfaces, demanding for a special treatment. A Mortar method is chosen in order to obtain a conservative and consistent load transfer. We validate our method by investigating two-dimensional test cases comprising a shock-loaded rigid cylinder and a deformable panel. Moreover, the aeroelastic instability of a thin plate structure is studied with a focus on the prediction of flutter onset. Finally, we propose a three-dimensional fluid–structure interaction test case of a flexible inflated thin shell interacting with a shock wave involving large and complex structural deformations.

42 ENGINEERING↗
DOE OSTI2017

A scalable nonlinear fluid–structure interaction solver based on a Schwarz preconditioner with isogeometric unstructured coarse spaces in 3D

Nonlinear fluid-structure interaction (FSI) problems on unstructured meshes in 3D appear many applications in science and engineering, such as vibration analysis of aircrafts and patient-specific diagnosis of cardiovascular diseases. In this work, we develop a highly scalable, parallel algorithmic and software framework for FSI problems consisting of a nonlinear fluid system and a nonlinear solid system, that are coupled monolithically. The FSI system is discretized by a stabilized finite element method in space and a fully implicit backward difference scheme in time. To solve the large, sparse system of nonlinear algebraic equations at each time step, we propose an inexact Newton-Krylov method together with a multilevel, smoothed Schwarz preconditioner with isogeometric coarse meshes generated by a geometry preserving coarsening algorithm. Here ''geometry'' includes the boundary of the computational domain and the wet interface between the fluid and the solid. We show numerically that the proposed algorithm and implementation are highly scalable in terms of the number of linear and nonlinear iterations and the total compute time on a supercomputer with more than 10,000 processor cores for several problems with hundreds of millions of unknowns.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Fabrication of nanotube arrays on commercially pure titanium and their apatite-forming ability in a simulated body fluid

In this study, we investigated self-organized TiO{sub 2} nanotubes that were grown using anodization of commercially pure titanium at 5 V or 10 V in NH{sub 4}F/NaCl electrolyte. The nanotube arrays were annealed at 450 °C for 3 h to convert the amorphous nanotubes to anatase and then they were immersed in simulated body fluid at 37 °C for 0.5, 1, and 14 days. The purpose of this experiment was to evaluate the apatite-formation abilities of anodized Ti nanotubes with different tube diameters and lengths. The nanotubes that formed on the surfaces of Ti were examined using a field emission scanning electron microscope, X-ray diffraction, and X-ray photoelectron spectroscope. When the anodizing potential was increased from 5 V to 10 V, the pore diameter of the nanotube increased from approximately 24–30 nm to 35–53 nm, and the tube length increased from approximately 590 nm to 730 nm. In vitro testing of the heat-treated nanotube arrays indicated that Ca-P formation occurred after only 1 day of immersion in simulated body fluid. This result was particularly apparent in the samples that were anodized at 10 V. It was also found that the thickness of the Ca-P layer increases as the applied potential for anodized c.p. Ti increases. The average thickness of the Ca-P layer on Ti that was anodized at 5 V and 10 V was approximately 170 nm and 190 nm, respectively, after immersion in simulated body fluid for 14 days. - Highlights: • TiO{sub 2} nanotube on Ti surface was formed by anodic oxidation in a NaCl/NH{sub 4}F solution. • TiO{sub 2} layers show a tube length of 590 nm and 730 nm at 5 V and 10 V, respectively. • After soaking in SBF, Ca-P layer completely covered the entire nanotubular surfaces. • The Ca-P layer was thicker on the Ti surface anodized at 10 V.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Approaches to multi-scale analyses of mechanically and thermally-driven migration of fluid inclusions in salt rocks

Fluid inclusions are found within mineral crystals or along grain boundaries in many sedimentary rocks, notably in evaporite formations, and can migrate along a thermal or hydro-mechanical gradient. Shale and salt rocks have been considered potential host rocks for radioactive waste disposal, due to their low permeability. Previously stagnant inclusions may become mobilised by a perturbation of the in situ state by a geotechnical installation or the emplacement of heat-generating waste. Here, the migration of fluid inclusions can thus have important impacts on the long-term performance of a geologic repository for high-level radioactive waste disposal. As a part of the international research project DECOVALEX-2019, two aspects of fluid inclusion migration in rock salt are currently investigated under different boundary conditions: a) altered hydro-mechanical conditions as a consequence of tunnel excavation or borehole drilling and b) coupled thermo-hydro-mechanical-chemical conditions during the heating period of the post-closure phase of a repository.

58 GEOSCIENCES↗
DOE OSTI2020

Fluid model for the piezothermal effect

When a gas in an externally imposed potential field is compressed, temperature gradients appear. This has been called the piezothermal effect. It is possible to analytically calculate the time-dependent behavior of the piezothermal effect using a linearized fluid model. Quantitative differences between the fluid-model results and previous numerical calculations can be explained by the effects of viscosity and heat conductivity. Here, the fluid model casts the piezothermal effect as a spectrum of buoyancy oscillations, which yields new physical insights into the effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Vanadium Geochemistry of Oil Sands Fluid Petroleum Coke

Vanadium has previously been linked to elevated toxicity of leachates derived from oil sands petroleum coke. However, geochemical controls on V mobility within coke deposits remain poorly constrained. Detailed examinations of porewater and solid-phase V geochemistry were therefore performed on oil sands fluid petroleum coke deposits in Alberta, Canada. Sample collection focused on both active and reclaimed deposits, which contained more than 3 × 10 7 m 3 of fluid petroleum coke. Dissolved V concentrations were highest (up to 3.0 mg L –1 ) immediately below the water table but decreased rapidly with increasing depth. This trend corresponded to a transition from mildly acidic (pH 6–7) and oxic conditions to mildly alkaline (pH 7–8.5) and anoxic conditions. Scanning electron microscopy (SEM), electron microprobe analysis (EMPA), and micro-X-ray fluorescence (μXRF) mapping revealed coke particles exhibited an internal structure characterized by successive concentric layers. Here, the outer margins of these layers were characterized by elevated V, Fe, Si, and Al concentrations, indicating the presence of inorganic phases. Micro-X-ray absorption near-edge structure (μXANES) spectroscopy revealed that V speciation was dominated by V(IV) porphyrins except at outer margins of layers, where octahedrally coordinated V(III) was a major component. Minor to trace V(V) was also detected within fluid petroleum coke particles.

58 GEOSCIENCES↗
DOE OSTI2021

Experiments on the Capillary Condensation/Evaporation Hysteresis of Pure Fluids and Binary Mixtures in Cylindrical Nanopores

Here, a new experimental method, i.e., isochoric cooling/heating cycle measurement using differential scanning calorimeter (DSC), is applied to investigate the hysteretic phenomenon of pure components and mixtures confined in two different types of cylindrical nanopores (SBA-15 and KIT-6). The capillary condensation/evaporation hysteresis end point (T he ) for a specific system is determined by performing a series of measurements with different densities by increasing initial temperature and pressure. We confirm that for a pure component, T he is lower when confined in smaller nanopores, whereas T he for CH 4 /CO 2 and CH 4 /C 2 H 6 mixtures at the specific composition investigated are higher than that of pure CO 2 and C 2 H 6 , respectively, when confined in the same nanopores. To the best of our knowledge, this is the first time that T he of fluid mixture is experimentally studied, which further demonstrates that binary mixtures behave similarly to pure fluids when confined in nanoporous media. At low temperatures where hysteresis is observed, the shape as well as the duration of the exothermic peak (due to capillary condensation) and endothermic peak (due to capillary evaporation) indicate that capillary condensation, rather than capillary evaporation, takes place at/near a thermodynamic equilibrium. Within the accuracy of the measurements, we find that the total heat of capillary condensation is identical to that of capillary evaporation for all systems examined, regardless of the presence of hysteresis. By comparing the heat per unit volume among different systems, we also find that the heat per unit volume of the phase transition of confined pure fluid is lower than that of bulk, which further demonstrates that the molar density of confined pure liquid is less than that of bulk pure liquid at saturation.

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

Complex Fluid‐Driven Fractures Caused by Crack‐Parallel Stress

Abstract Managing fluid‐driven fracture networks is crucial for subsurface resource utilization, yet the current understanding of the key controlling factors remains insufficient. While geologic discontinuities have been shown to significantly influence fracture network complexity, this study identifies another major contributor. We conducted a new set of experiments using a transparent true triaxial cell, which enabled video recording of the temporal evolution of fluid‐driven fracture paths. Using pseudo‐2D samples without macroscale structural discontinuities, we observed multiple occurrences of hydraulic fracture curving and branching under anisotropic boundary stresses. We proposed a theoretical model demonstrating that the stress parallel to the crack line in the solid matrix near the crack tip (i.e., the T ‐stress) accounts for the observed fracture curving behavior. This finding suggests that T ‐stress is an additional mechanism contributing to the complexity of fluid‐driven fracture networks in the subsurface, besides the geologic discontinuities.

58 GEOSCIENCES↗
DOE OSTI2016

Modulated phases of graphene quantum Hall polariton fluids

There is a growing experimental interest in coupling cavity photons to the cyclotron resonance excitations of electron liquids in high-mobility semiconductor quantum wells or graphene sheets. These media offer unique platforms to carry out fundamental studies of exciton-polariton condensation and cavity quantum electrodynamics in a regime, in which electron–electron interactions are expected to play a pivotal role. Here, focusing on graphene, we present a theoretical study of the impact of electron–electron interactions on a quantum Hall polariton fluid, that is a fluid of magneto-excitons resonantly coupled to cavity photons. We show that electron–electron interactions are responsible for an instability of graphene integer quantum Hall polariton fluids towards a modulated phase. We demonstrate that this phase can be detected by measuring the collective excitation spectra, which is often at a characteristic wave vector of the order of the inverse magnetic length.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Methods to capture proteomic and metabolomic signatures from cerebrospinal fluid and serum of healthy individuals

Discovery of reliable signatures for the empirical diagnosis of neurological diseases—both infectious and non-infectious—remains unrealized. One of the primary challenges encountered in such studies is the lack of a comprehensive database representative of a signature background that exists in healthy individuals, and against which an aberrant event can be assessed. For neurological insults and injuries, it is important to understand the normal profile in the neuronal (cerebrospinal fluid) and systemic fluids (e.g., blood). Here, we present the first comparative multi-omic human database of signatures derived from a population of 30 individuals (15 males, 15 females, 23–74 years) of serum and cerebrospinal fluid. In addition to empirical signatures, we also assigned common pathways between serum and CSF. Together, our findings provide a cohort against which aberrant signature profiles in individuals with neurological injuries/disease can be assessed—providing a pathway for comprehensive diagnostics and therapeutics discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

A Hamiltonian fluid-kinetic model for a two-species non-neutral plasma

A model for describing the dynamics of a pure electron plasma in the presence of a population of massive charged particles is presented. The model couples the fluid dynamics of the pure electron plasma with the dynamics of the massive particle population, the latter being treated kinetically. The model is shown to possess a noncanonical Hamiltonian structure and to preserve invariants analogous to those of the two-dimensional (2D) Euler equation for an incompressible inviscid fluid, and of the Vlasov equation. The Hamiltonian structure of the model is used to derive a set of stability conditions for rotating coherent structures of the two-species system, in the case of negatively charged massive particles. According to these conditions, stability is attained if both the equilibrium distribution function of the kinetic species and the equilibrium density of the electron fluid are monotonically decreasing functions of the corresponding single-particle energies in the rotating frame. For radially confined equilibria near the axis, the stability condition corresponds to the existence of a finite interval of rotation frequencies for the reference frame, with the upper bound determined by the presence of the kinetic population.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Two-fluid description of wave-particle interactions in strong Buneman turbulence

To understand the nature of anomalous resistivity in magnetic reconnection, we investigate turbulence-induced momentum transport and energy dissipation while a plasma is unstable to the Buneman instability in force-free current sheets. Using 3D particle-in-cell simulations, we find that the macroscopic effects generated by wave-particle interactions in Buneman instability can be approximately described by a set of electron fluid equations. We show that both energy dissipation and momentum transport along electric current in the current layer are locally quasi-static, but globally dynamic and irreversible. Turbulent drag dissipates both the streaming energy of the current sheet and the associated magnetic energy. The net loss of streaming energy is converted into the electron component heat conduction parallel to the magnetic field and increases the electron Boltzmann entropy. The growth of self-sustained Buneman waves satisfies a Bernoulli-like equation that relates the turbulence-induced convective momentum transport and thermal momentum transport. Electron trapping and de-trapping drive local momentum transports, while phase mixing converts convective momentum into thermal momentum. The drag acts like a micro-macro link in the anomalous heating processes. The decrease of magnetic field maintains an inductive electric field that re-accelerates electrons, but most of the magnetic energy is dissipated and converted into the component heat of electrons perpendicular to the magnetic field. This heating process is decoupled from the heating of Buneman instability in the current sheets. Ion heating is weak but ions play an important role in assisting energy exchanges between waves and electrons. Cold ion fluid equations together with our electron fluid equations form a complete set of equations that describes the occurrence, growth, saturation and decay of the Buneman instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Dynamics of a confined dusty fluid in a sheared ion flow

Dynamics of an isothermally driven dust fluid is analyzed which is confined in an azimuthally symmetric cylindrical setup by an effective potential and is in equilibrium with an unconfined sheared flow of a streaming plasma. Cases are analyzed where the confining potential constitutes a barrier for the driven fluid, limiting its spatial extension and boundary velocity. The boundary effects entering the formulation are characterized by applying the appropriate boundary conditions and a range of solutions exhibiting single and multiple vortex are obtained. The equilibrium solutions considered in the cylindrical setup feature a transition from single to multiple vortex state of the driven flow. Effects of (i) the variation in dust viscosity, (ii) coupling between the driving and the driven fluid, and (iii) a friction determining the equilibrium dynamics of the driven system are characterized.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Radiation induces turbulence in particle-laden fluids

When a transparent fluid laden with solid particles is subject to radiative heating, non-uniformities in particle distribution result in local fluid temperature fluctuations. Under the influence of gravity, buoyancy induces vortical fluid motion which can lead to strong preferential concentration, enhancing the local heating and more non-uniformities in particle distribution. By employing direct numerical simulations this study shows that the described feedback loop can create and sustain turbulence. The velocity and length scale of the resulting turbulence is not known a priori, and is set by balance between viscous forces and buoyancy effects. When the particle response time is comparable to a viscous time scale, introduced in our analysis, the system exhibits intense fluctuations of turbulent kinetic energy and strong preferential concentration of particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Communication: Superstabilization of fluids in nanocontainers

One of the main challenges of thermodynamics is to predict and measure accurately the properties of metastable fluids. Investigation of these fluids is hindered by their spontaneous transformation by nucleation into a more stable phase. We show how small closed containers can be used to completely prevent nucleation, achieving infinitely long-lived metastable states. Using a general thermodynamic framework, we derive simple formulas to predict accurately the conditions (container sizes) at which this superstabilization takes place and it becomes impossible to form a new stable phase. This phenomenon opens the door to control nucleation of deeply metastable fluids at experimentally feasible conditions, having important implications in a wide variety of fields.

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

The mold integration method for the calculation of the crystal-fluid interfacial free energy from simulations

The interfacial free energy between a crystal and a fluid, γ{sub cf}, is a highly relevant parameter in phenomena such as wetting or crystal nucleation and growth. Due to the difficulty of measuring γ{sub cf} experimentally, computer simulations are often used to study the crystal-fluid interface. Here, we present a novel simulation methodology for the calculation of γ{sub cf}. The methodology consists in using a mold composed of potential energy wells to induce the formation of a crystal slab in the fluid at coexistence conditions. This induction is done along a reversible pathway along which the free energy difference between the initial and the final states is obtained by means of thermodynamic integration. The structure of the mold is given by that of the crystal lattice planes, which allows to easily obtain the free energy for different crystal orientations. The method is validated by calculating γ{sub cf} for previously studied systems, namely, the hard spheres and the Lennard-Jones systems. Our results for the latter show that the method is accurate enough to deal with the anisotropy of γ{sub cf} with respect to the crystal orientation. We also calculate γ{sub cf} for a recently proposed continuous version of the hard sphere potential and obtain the same γ{sub cf} as for the pure hard sphere system. The method can be implemented both in Monte Carlo and Molecular Dynamics. In fact, we show that it can be easily used in combination with the popular Molecular Dynamics package GROMACS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗