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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A generalized two-fluid picture of non-driven collisionless reconnection and its relation to whistler waves

A generalized, intuitive two-fluid picture of 2D non-driven collisionless magnetic reconnection is described using results from a full-3D numerical simulation. The relevant two-fluid equations simplify to the condition that the flux associated with canonical circulation Q=m e ∇×u e +q e B is perfectly frozen into the electron fluid. In the reconnection geometry, flux tubes defined by Q are convected with the central electron current, effectively stretching the tubes and increasing the magnitude of Q exponentially. This, coupled with the fact that Q is a sum of two quantities, explains how the magnetic fields in the reconnection region reconnect and give rise to strong electron acceleration. The Q motion provides an interpretation for other phenomena as well, such as spiked central electron current filaments. The simulated reconnection rate was found to agree with a previous analytical calculation having the same geometry. Energy analysis shows that the magnetic energy is converted and propagated mainly in the form of the Poynting flux, and helicity analysis shows that the canonical helicity ∫P·Q dV as a whole must be considered when analyzing reconnection. A mechanism for whistler wave generation and propagation is also described, with comparisons to recent spacecraft observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

A simulation method for the phase diagram of complex fluid mixtures

We present that the phase behavior of complex fluid mixtures is of continuing interest, but obtaining the phase diagram from computer simulations can be challenging. In the Gibbs ensemble method, for example, each of the coexisting phases is simulated in a different cell, and ensuring the equality of chemical potentials of all components requires the transfer of molecules from one cell to the other. For complex fluids such as polymers, successful insertions are rare. An alternative method is to simulate both coexisting phases in a single simulation cell, with an interface between them. The challenge here is that the interface position moves during the simulation, making it difficult to determine the concentration profile and coexisting concentrations. In this work, we propose a new method for single cell simulations that uses a spatial concentration autocorrelation function to (spatially) align instantaneous concentration profiles from different snapshots. This allows one to obtain average concentration profiles and hence the coexisting concentrations. We test the method by calculating the phase diagrams of two systems: the Widom-Rowlinson model and the symmetric blends of freely jointed polymer molecules for which phase diagrams from conventional methods are available. Excellent agreement is found, except in the neighborhood of the critical point where the interface is broad and finite size effects are important. Finally, the method is easy to implement and readily applied to any mixture of complex fluids.

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

A Landau-fluid closure for arbitrary frequency response

The kinetic Landau-fluid (LF) closure which can be regarded as the exact closure is derived. For Maxwellian plasma, the kinetic closure is the same as Hammett-Perkins closure in static limit and totally the same with Chang-Callen closure. A new LF closure for arbitrary frequency response constructed with harmonic average technique is presented in this paper. This new LF closure bridges the existing LF closures in the low and high frequency limits: it recovers Hammett-Perkins closure when weight coefficient κ = 0 and converges to Chang-Callen closure at high frequency when weight coefficient κ = 1. By picking an appropriate κ, the harmonic average closure contains both nonlocal transport and local transport and the resulting fluid response function of a three moment fluid model well matches the exact kinetic response function within the entire frequency range. On the computational side, a sum of diffusion-convection solves (SDCS) method is developed to facilitate numerical implementation of the harmonic average LF closure. By using SDCS method, good agreement is achieved for the response functions between driven initial-value simulations and matrix eigenvalue calculations within the BOUT++ framework. The harmonic average closure of shifted-Maxwellian is also outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

High-order two-fluid plasma solver for direct numerical simulations of plasma flows with full transport phenomena

The two-fluid plasma equations for a single ion species, with full transport terms, including temperature and magnetic field dependent ion and electron viscous stresses and heat fluxes, frictional drag force, and ohmic heating terms, have been implemented in the CFDNS code and solved by using sixth-order non-dissipative compact finite differences for plasma flows in several different regimes. In order to be able to fully resolve all the dynamically relevant time and length scales, while maintaining computational feasibility, the assumptions of infinite speed of light and negligible electron inertia have been made. Non-dimensional analysis of the two-fluid plasma equations shows that, by varying the characteristic/background number density, length scale, temperature, and magnetic strength, the corresponding Hall, resistive, and ideal magnetohydrodynamic equations can be recovered as limiting cases. The accuracy and robustness of this two-fluid plasma solver in handling plasma flows in different regimes have been validated against four canonical problems: Alfven and whistler dispersion relations, electromagnetic plasma shock, and magnetic reconnection. For all test cases, by using physical dissipation and diffusion, with negligible numerical dissipation/diffusion, fully converged Direct Numerical Simulation (DNS)-like solutions are obtained when the ion Reynolds number based on the grid size is smaller than a threshold value which is about 2.3. For the magnetic reconnection problem, the results show that the magnetic flux saturation time and value converge when the ion and magnetic Reynolds numbers are large enough. Thus, the DNS-like results become relevant to practical problems with much larger Reynolds numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Deep learning surrogate model for kinetic Landau-fluid closure with collision

In this work, the kinetic Landau-fluid (LF) closure with collision and periodic boundary condition is used in the development of the deep learning (DL) surrogate model. A classical neural network, namely, feedforward neural network or sometimes termed multilayer perceptron, is constructed and trained to learn the kinetic LF closure in the static limit and arbitrary mean free path in configuration space. The preliminary relation between best hyperparameters and critical parameters for data generation is found. Compared with the numerical approach (non-Fourier method) of the LF closure, the deep learning surrogate model shows an order of magnitude of improvement in terms of accuracy. Perhaps most importantly, the surrogate model closure has been integrated for the first time with fluid simulations. Our DL-enabled fluid simulations, for the first time, give the correct Landau damping rate for a wide range of wave vectors, while the Hammett–Perkins closure cannot produce the correct damping rate. We correctly connect the collisionless Hammett–Perkins closure and collisional Braginskii closure to reproduce the intrinsic nonlocal feature of the heat flux with DL techniques. We address the most concerning error accumulation problem and find that simulations with the deep learning surrogate model are as good as, if not better than, simulations with the analytic closure in terms of long-term numerical stability in the linear Landau damping test.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fluid-like shock compression of dilute polymer nanocomposites

The shock Hugoniot of heterogenous mixtures of discrete particles has been experimentally investigated for porous agglomerates and fluids, while the study of full-density solids has been primarily limited to compressed powders and high-dimensional composites. By dispersing ceria nanoparticles in a polyethylene matrix, we are able to examine the hydrodynamic behavior of a nonporous, heterogenous solid in thermal equilibrium during weak shock compression. Phase-driven discontinuities in the Hugoniot particle velocity–shock velocity (u−D) relationship of pure polyethylene are replicated in the nanocomposites but are shifted to lower velocity and to higher pressure with higher particle concentration. The results are explained using an isothermal, two-velocity fluid model under the hydrodynamic approximation. The model, which assumes a theoretical equation-of-state for ceria and either a low-order or high-order fit to the measured polyethylene Hugoniot, reasonably predicts the Hugoniot for two different polyethylene/ceria nanocomposites. Using the model, the mixture Hugoniot is shown to be insensitive to the Hugoniot of the stiffer constituent when the moduli are sufficiently disparate, while dependence on particle density and volume fraction is preserved through fluid-like motion.

Moore, Nathan W. [Sandia National Laboratories (SN↗
DOE OSTI2025

Coupled 2D MHD and runaway electron fluid simulations of SPARC disruptions

Runaway electrons (REs) generated during disruption events in tokamaks can carry mega-Ampère level currents, potentially causing damage to plasma-facing components. Understanding RE evolution during disruption events is important for evaluating strategies to mitigate RE damage. Using two-dimensional toroidally symmetric magnetohydrodynamic (MHD) simulations in M3D-C1, which incorporates a fluid RE model evolved self-consistently with the bulk MHD fluid, we examine the seeding and avalanching of REs during disruptions in the SPARC tokamak – a compact, high-field, high-current device designed to achieve a fusion gain Q > 2 in deuterium–tritium plasmas. The M3D-C1 simulations of unmitigated disruptions demonstrate RE plateau formation and peaking of the final current density, which agree well with the results of lower-fidelity reduced RE fluid models. This work provides the first systematic comparison and benchmarking of different primary sources, including activated tritium beta decay and Compton scattering, in SPARC disruption simulations with self-consistent MHD and RE coupling.

Compton scattering↗
DOE OSTI2025

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Implications of a two-fluid heat diffusion model for RF condensation effects in magnetic islands

Theoretical calculations have predicted an “RF condensation” effect through which the temperature perturbation and current drive may be nonlinearly enhanced in electron cyclotron (EC) heated magnetic islands. Previously, two-fluid effects had been neglected in the OCCAMI code used to model RF condensation. For realistic experimental conditions, we show that two-fluid effects in an EC heated plasma may have a major impact on RF condensation. Less collisional plasmas, where electron and ion temperatures are decoupled, exhibit an enhanced linear temperature perturbation in the island, and an even stronger nonlinear enhancement of the perturbed temperature. This can occur even when nonlinear effects are small in the single-fluid limit. We also calculate corresponding hysteresis effects. Finally, we demonstrate that other physical effects, including the rotation of the magnetic island and the Gaussian spread of the beam from a gyrotron, can significantly decrease the predicted observable RF condensation effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Control of optical transport parameters of 'porous medium – supercritical fluid' systems

The possibility of controlling optical transport parameters (in particular, transport scattering coefficient) of porous systems based on polymer fibres, saturated with carbon dioxide in different phase states (gaseous, liquid and supercritical) has been experimentally studied. An increase in the pressure of the saturating medium leads to a rise of its refractive index and, correspondingly, the diffuse-transmission coefficient of the system due to the decrease in the transport scattering coefficient. It is shown that, in the case of subcritical saturating carbon dioxide, the small-angle diffuse transmission of probed porous layers at pressures close to the saturated vapour pressure is determined by the effect of capillary condensation in pores. The immersion effect in 'porous medium – supercritical fluid' systems, where the fluid pressure is used as a control parameter, is considered. The results of reconstructing the values of transport scattering coefficient of probed layers for different refractive indices of a saturating fluid are presented. (radiation scattering)

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

Tolerable ranges of fluid shear for early life-stage fishes: implications for safe fish passage at hydropower and irrigation infrastructure

Egg and larval fish drifting downstream are likely to encounter river infrastructure, increasing mortality. Elevated fluid shear is one likely cause of this. To confirm this and determine tolerable levels, egg and larvae of three Australian species were exposed to elevated shear in a laboratory flume. Mortality was modelled over a broad range of strain rates, allowing critical thresholds to be estimated. Eggs were very susceptible to mortality at low strain rates and 100% of golden and silver perch died once strain rate exceeded 629 and 148 cm s-1 cm-1, respectively. Larvae were less vulnerable than eggs, but mortality increased at higher strain rates and at younger ages. Murray cod larvae were the most susceptible, particularly around 9 days post hatch. Younger golden perch and silver perch were also susceptible but little mortality occurred once maturity reached ~25 days post hatch (nearing juvenile metamorphosis). The results indicate that fluid shear will contribute to the large mortalities of drifting eggs and larvae entrained at river infrastructure, with the impact reducing with age. Minimising turbulent flows and fluid shear should be the objective when designing and operating structures like spillways, turbines and irrigation weirs in areas of significant egg and larval drift.

Navarro, Anna↗
DOE OSTI2020

Complex subsurface hydrothermal fluid mixing at a submarine arc volcano supports distinct and highly diverse microbial communities

Hydrothermally active submarine volcanoes are mineral-rich biological oases contributing significantly to chemical fluxes in the deep sea, yet little is known about the microbial communities inhabiting these systems. Here we investigate the diversity of microbial life in hydrothermal deposits and their metagenomics-inferred physiology in light of the geological history and resulting hydrothermal fluid paths in the subsurface of Brothers submarine volcano north of New Zealand on the southern Kermadec arc. From metagenome-assembled genomes we identified over 90 putative bacterial and archaeal genomic families and nearly 300 previously unknown genera, many potentially endemic to this submarine volcanic environment. While magmatically influenced hydrothermal systems on the volcanic resurgent cones of Brothers volcano harbor communities of thermoacidophiles and diverse members of the superphylum “DPANN,” two distinct communities are associated with the caldera wall, likely shaped by two different types of hydrothermal circulation. The communities whose phylogenetic diversity primarily aligns with that of the cone sites and magmatically influenced hydrothermal systems elsewhere are characterized predominately by anaerobic metabolisms. These populations are probably maintained by fluids with greater magmatic inputs that have interacted with different (deeper) previously altered mineral assemblages. However, proximal (a few meters distant) communities with gene-inferred aerobic, microaerophilic, and anaerobic metabolisms are likely supported by shallower seawater-dominated circulation. Furthermore, mixing of fluids from these two distinct hydrothermal circulation systems may have an underlying imprint on the high microbial phylogenomic diversity. Collectively our results highlight the importance of considering geologic evolution and history of subsurface processes in studying microbial colonization and community dynamics in volcanic environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Fluid inertia controls mineral precipitation and clogging in pore to network-scale flows

Mineral precipitation caused by fluid mixing presents complex control and predictability challenges in a variety of natural and engineering processes, including carbon mineralization, geothermal energy, and microfluidics. Precipitation dynamics, particularly under the influence of fluid flow, remain poorly understood. Combining microfluidic experiments and three-dimensional reactive transport simulations, we demonstrate that fluid inertia controls mineral precipitation and clogging at flow intersections, even in laminar flows. We observe distinct precipitation regimes as a function of Reynolds number (Re). At low Reynolds numbers (Re < 10), precipitates form a thin, dense layer along the mixing interface, which shuts precipitation off, while at high Reynolds numbers (Re > 50), strong three-dimensional flows significantly enhance precipitation over the entire intersection, resulting in rapid clogging. When injection rates from two inlets are uneven, flow symmetry-breaking leads to unexpected flow bifurcation phenomena, which result in enhanced concurrent precipitation in both downstream channels. Finally, we extend our findings to rough channel networks and demonstrate that the identified inertial effects on precipitation at the intersection scale are also present and even more dramatic at the network scale. This study sheds light on the fundamental mechanisms underlying mixing-induced mineral precipitation and provides a framework for designing and optimizing processes involving mineral precipitation.

Science & Technology - Other Topics↗
DOE OSTI2018

Modified Point-Kinetics Model for Neutron Precursors and Fission Product Behavior for Fluid-Fueled Molten Salt Reactors

Fluid-fueled nuclear reactors, such as molten salt reactors (MSRs), have recently gained significant interest. These advanced reactors represent a potential revolutionary shift in the implementation of nuclear power, and as a broad class of reactors, they have the potential to directly address many U.S. energy policy objectives. Fuel that is dissolved in the coolant requires methods to account for the birth, decay, and transport of fission products not only in the core but also throughout the loop and any auxiliary systems, such as off-gas, to which liquid fuel flows, gaseous products are carried, or solid particulates plate out. System models are particularly well suited to explore the wide range of phenomena that are associated with fluid-fueled systems, especially for safeguards analysis. However, before system dynamics can be explored, the compositions of fission products of the salt throughout the loop must be determined as they drive the dynamic behavior of a reactor.This paper describes the derivation of a modified point-kinetics model for obtaining a first-order approximation of the behavior of a salt-fueled system in which neutron precursors and fission products are born in the fuel-salt and transported outside the core. This paper also provides verification of the model using a steady-state analytic solution and provides additional cases exploring the response under transient cases. This model establishes a baseline model that can be used to explore the dynamic response of fluid-fueled reactors and to investigate important safeguards issues such as mass accountability of source terms. Here, the model is implemented in the Oak Ridge National Laboratory–developed, Modelica-based TRANSFORM library that was developed to investigate various aspects of advanced energy systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

A modified many-body dissipative particle dynamics model for mesoscopic fluid simulation: methodology, calibration, and application for hydrocarbon and water

The many-body dissipative particle dynamics (mDPD) is a prominent mesoscopic multiphase model for fluid transport in mesoconfinement. However, it has been a long-standing challenge for mDPD (and other multiphase-enabled DPD models) to accurately predict real-fluid static and dynamic properties simultaneously. We have developed a modified mDPD model to overcome the issue and a rigorous calibration approach that uses reference data, including experimental and/or molecular dynamics (MD) simulations to parameterise the modified mDPD for real fluids. We choose heptane as a representative example of hydrocarbon in source rocks to demonstrate the model's capability to accurately predict the equation of state (EOS), free surface tension, diffusivity, and viscosity. Our timing test shows that the modified mDPD is 400–500 times faster than its MD counterpart for simulating bulk heptane in equivalent volumes. To further demonstrate the robustness of the model, we revisited the benchmark problem of mesoscopic modelling of water, in which all the previous works on DPD reported only a limited portion of the water properties. Here, we show that the modified mDPD can provide accurate modelling of water static and dynamic properties and an EOS that matches the experimental data to a large range of confinement pressure.

42 ENGINEERING↗
DOE OSTI2014

Cosmological fluctuations of a random field and radiation fluid

A generalization of the random fluid hydrodynamic fluctuation theory due to Landau and Lifshitz is applied to describe cosmological fluctuations in systems with radiation and scalar fields. The viscous pressures, parametrized in terms of the bulk and shear viscosity coefficients, and the respective random fluctuations in the radiation fluid are combined with the stochastic and dissipative scalar evolution equation. This results in a complete set of equations describing the perturbations in both scalar and radiation fluids. These derived equations are then studied, as an example, in the context of warm inflation. Similar treatments can be done for other cosmological early universe scenarios involving thermal or statistical fluctuations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗
DOE OSTI2019

Two-fluid dynamics of one-dimensional quantum liquids in the absence of Galilean invariance

Luttinger liquid theory of one-dimensional quantum systems ignores exponentially weak backscattering of particles. This endows Luttinger liquids with superfluid properties. The corresponding two-fluid hydrodynamic description available at present applies only to Galilean-invariant systems, whereas most experimental realizations of one-dimensional quantum liquids lack Galilean invariance. In this work, we develop the two-fluid theory of such quantum liquids. In the low-frequency limit the theory reduces to single-fluid hydrodynamics. However, the absence of Galilean invariance brings about three new transport coefficients. Lastly, we obtain expressions for these coefficients in terms of the backscattering rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗