Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,261 records · Page 70

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗

Localization of deformation in fluid-saturated sandstone

Herein we discuss a laboratory study of inelastic deformation and localization was conducted to evaluate microseismicity and fracture in fluid-saturated rock. Plane strain compression experiments were performed on oil- and water-saturated Berea sandstone under dry, drained, undrained, and quasi-unjacketed conditions. Parameters associated with the elastic and inelastic material response were accurately measured and used in a constitutive model to predict localization of deformation under the drained and undrained conditions. Acoustic emission (AE) activity was recorded and the onset of inelastic response in fluid-saturated rock, similar to dry rock, coincided with an increase in AE rate. Fracture initiation and propagation were successfully monitored by AE locations. Clustering of AE events and the start of inhomogeneous deformation inferred from displacement measurements are shown to be well correlated with the model prediction of the onset of localized deformation under various rock-fluid interaction conditions.

42 ENGINEERING↗

Classification of machine learning frameworks for data-driven thermal fluid models

Thermal fluid processes are inherently multi-physics and multi-scale, involving mass-momentum-energy transport phenomena at multiple scales. Thermal fluid simulation (TFS) is based on solving conservative equations, for which – except for “first-principles” direct numerical simulation – closure relations (CRs) are required to provide microscopic interactions or so-called sub-grid-scale physics. In practice, TFS is realized through reduced-order modeling, and its CRs as low-fidelity models can be informed by observations and data from relevant and adequately evaluated experiments and high-fidelity simulations. This paper is focused on data-driven TFS models, specifically on their development using machine learning (ML). Five ML frameworks are introduced including physics-separated ML (PSML or Type I ML), physics-evaluated ML (PEML or Type II ML), physics-integrated ML (PIML or Type III ML), physics-recovered (PRML or Type IV ML), and physics-discovered ML (PDML or Type V ML). The frameworks vary in their performance for different applications depending on the level of knowledge of governing physics, source, type, amount and quality of available data for training. Notably, outlined for the first time in this paper, Type III models present stringent requirements on modeling, substantial computing resources for training, and high potential in extracting value from “big data” in thermal fluid research. Here the current paper demonstrates and investigates ML frameworks in three examples. First, we utilize the heat diffusion equation with a nonlinear conductivity model formulated by convolutional neural networks (CNNs) and feedforward neural networks (FNNs) to illustrate the applications of Type I, Type II, Type III, and Type V ML. The results indicate a preference for Type II ML under deficient data support. Type III ML can effectively utilize field data, potentially generating more robust predictions than Type I and Type II ML. CNN-based closures exhibit more predictability than FNN-based closures, but CNN-based closures require more training data to obtain accurate predictions. Second, we illustrate how to employ Type I ML and Type II ML frameworks for data-driven turbulence modeling using reference works. Third, we demonstrate Type I ML by building a deep FNN-based slip closure for two-phase flow modeling. The results show that deep FNN-based closures exhibit a bounded error in the prediction domain.

42 ENGINEERING↗

Multi-solid and multi-fluid diffuse interface model: Applications to dynamic fracture and fragmentation

We extend the model of diffuse solid–fluid interfaces developed earlier by authors of this paper to the case of arbitrary number of interacting hyperelastic solids. Plastic transformations of solids are taken into account through a Maxwell type model. The specific energy of each solid is given in separable form: it is the sum of a hydrodynamic part of the energy depending only on the density and the entropy, and an elastic part of the energy which is unaffected by the volume change. It allows us to naturally pass to the fluid description in the limit of vanishing shear modulus. In spite of a large number of governing equations, the model has a quite simple mathematical structure: it is a duplication of a single visco-elastic model. The model is well posed both mathematically and thermodynamically: it is hyperbolic and compatible with the second law of thermodynamics. The resulting model can be applied in the situations involving an arbitrary number of fluids and solids. In particular, we show the ability of the model to describe spallation and penetration phenomena occurring during high velocity impacts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Wave resonances and the partition of energy in ideal compressible magnetohydrodynamic fluids

Phase mixing and resonant absorption are two processes that have been under scrutiny for some time because of their role in wave damping and in providing a mechanism for heating space and laboratory plasmas. The accumulation or absorption of energy that develops within a resonant layer is usually attributed to a logarithmic singularity, but it will be shown that this build up of energy is inextricably tied to a discontinuity in the fluid displacement at the resonant point. This change in the dynamics of the problem will be examined by establishing a partition of energy that identifies and isolates the individual resonances within the fluid. The partition is based on a variational description of the Fourier transformed equations and is guided by an electrical model of the MHD system that not only illustrates the resonant structure threading the fluid but also exposes the mechanism driving the resonant absorption process. A simplified version of this model is then constructed to help determine the approximate rate of energy absorption.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computer simulation of liquid-vapor coexistence of confined quantum fluids

The liquid-vapor coexistence (LV) of bulk and confined quantum fluids has been studied by Monte Carlo computer simulation for particles interacting via a semiclassical effective pair potential V{sub eff}(r) = V{sub LJ} + V{sub Q}, where V{sub LJ} is the Lennard-Jones 12-6 potential (LJ) and V{sub Q} is the first-order Wigner-Kirkwood (WK-1) quantum potential, that depends on β = 1/kT and de Boer's quantumness parameter Λ=h/σ√(mε), where k and h are the Boltzmann's and Planck's constants, respectively, m is the particle's mass, T is the temperature of the system, and σ and ε are the LJ potential parameters. The non-conformal properties of the system of particles interacting via the effective pair potential V{sub eff}(r) are due to Λ, since the LV phase diagram is modified by varying Λ. We found that the WK-1 system gives an accurate description of the LV coexistence for bulk phases of several quantum fluids, obtained by the Gibbs Ensemble Monte Carlo method (GEMC). Confinement effects were introduced using the Canonical Ensemble (NVT) to simulate quantum fluids contained within parallel hard walls separated by a distance L{sub p}, within the range 2σ ⩽ L{sub p} ⩽ 6σ. The critical temperature of the system is reduced by decreasing L{sub p} and increasing Λ, and the liquid-vapor transition is not longer observed for L{sub p}/σ < 2, in contrast to what has been observed for the classical system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗