Search NASA⌕ Search

SEARCH · Search NASA

Results for “general fluid mechanics”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Identifying and tracking bubbles and drops in simulations: A toolbox for obtaining sizes, lineages, and breakup and coalescence statistics

Knowledge of bubble and drop size distributions in two-phase flows is important for characterizing a wide range of phenomena, including combustor ignition, sonar communication, and cloud formation. The physical mechanisms driving the background flow also drive the time evolution of these distributions. Accurate and robust identification and tracking algorithms for the dispersed phase are necessary to reliably measure this evolution and thereby quantify the underlying mechanisms in interface-resolving flow simulations. The identification of individual bubbles and drops traditionally relies on an algorithm used to identify connected regions. This traditional algorithm can be sensitive to the presence of spurious structures. A cost-effective refinement is proposed to maximize volume accuracy while minimizing the identification of spurious bubbles and drops. An accurate identification scheme is crucial for distinguishing bubble and drop pairs with large size ratios. The identified bubbles and drops need to be tracked in time to obtain breakup and coalescence statistics that characterize the evolution of the size distribution, including breakup and coalescence frequencies, and the probability distributions of parent and child bubble and drop sizes. An algorithm based on mass conservation is proposed to construct bubble and drop lineages using simulation snapshots that are not necessarily from consecutive time steps. These lineages are then used to detect breakup and coalescence events, and obtain the desired statistics. Accurate identification of large-size-ratio bubble and drop pairs enables accurate detection of breakup and coalescence events over a large size range. Accurate detection of successive breakup and coalescence events requires that the snapshot interval be an order of magnitude smaller than the characteristic breakup and coalescence times to capture these successive events while minimizing the identification of repeated confounding events. Together, these algorithms serve as a toolbox for detailed analysis of two-phase simulations, and enable insights into the mechanisms behind bubble and drop formation and evolution in flows of practical importance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Transferability of data-driven, many-body models for CO 2 simulations in the vapor and liquid phases

Here, extending on the previous work by Riera et al. [J. Chem. Theory Comput. 16, 2246–2257 (2020)], we introduce a second generation family of data-driven many-body MB-nrg models for CO 2 and systematically assess how the strength and anisotropy of the CO 2 –CO 2 interactions affect the models’ ability to predict vapor, liquid, and vapor–liquid equilibrium properties. Building upon the many-body expansion formalism, we construct a series of MB-nrg models by fitting one-body and two-body reference energies calculated at the coupled cluster level of theory for large monomer and dimer training sets. Advancing from the first generation models, we employ the charge model 5 scheme to determine the atomic charges and systematically scale the two-body energies to obtain more accurate descriptions of vapor, liquid, and vapor–liquid equilibrium properties. Challenges in model construction arise due to the anisotropic nature and small magnitude of the interaction energies in CO 2 , calling for the necessity of highly accurate descriptions of the multidimensional energy landscape of liquid CO 2 . These findings emphasize the key role played by the training set quality in the development of transferable, data-driven models, which, accurately representing high-dimensional many-body effects, can enable predictive computer simulations of molecular fluids across the entire phase diagram.

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

The oscillatory anisotropy in the spatially flat cosmological models

Highlights: • Relations among asymptotic anisotropy and matter content. • A new scheme for obtaining physically motivated exact solutions. • The process of isotropization/anisotropization in cosmology. • The conversion among the anisotropy and the energy density into one another. In this manuscript, we investigate the oscillatory behaviour of the anisotropy in the diagonal Bianchi-I spacetimes. Our starting point is a simplification of Einstein’s equations using only observable or physical variables. As a consequence, we are able to: (a) Prove general results concerning the existence of oscillations of the anisotropy in the primordial and the late-time universe. For instance, in the expanding scenario, we show that a past weakly mixmaster behaviour (oscillations as we approach the Kasner solutions) might appear even with no violation of the usual energy conditions, while in the future, the pulsation (oscillations around isotropic solutions) seems to be most favoured; (b) Determine a large scheme for deriving classes of physically motivated exact solutions, and we give some (including the general barotropic perfect fluid and the magnetic one); (c) Understand the physical conditions for the occurrence of the isotropization or anisotropization during the cosmological evolution; (d) Understand how anisotropy and energy density are converted one into another. In particular, we call attention to the presence of a residue in the energy density in a late-time isotropic universe coming from its past anisotropic behaviour.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Clogging, dynamics, and reentrant fluid for active matter on periodic substrates

Here, we examine the collective states of run-and-tumble active matter disks driven over a periodic obstacle array. When the drive is applied along a symmetry direction of the array, we find a clog-free uniform liquid state for low activity, while at higher activity, the density becomes increasingly heterogeneous and an active clogged state emerges in which the mobility is strongly reduced. For driving along nonsymmetry or incommensurate directions, there are two different clogging behaviors consisting of a drive-dependent clogged state in the low activity thermal limit and a drive-independent clogged state at high activity. These regimes are separated by a uniform flowing liquid at intermediate activity. There is a critical activity level above which the thermal clogged state does not occur, as well as an optimal activity level that maximizes the disk mobility. Thermal clogged states are dependent on the driving direction while active clogged states are not. In the low activity regime, diluting the obstacles produces a monotonic increase in the mobility; however, for large activities, the mobility is more robust against obstacle dilution. We also examine the velocity-force curves for driving along nonsymmetry directions and find that they are linear when the activity is low or intermediate but become nonlinear at high activity and show behavior similar to that found for the plastic depinning of solids. At higher drives, the active clustering is lost. For low activity, we also find a reentrant fluid phase, where the system transitions from a high mobility fluid at low drives to a clogged state at higher drives and then back into another fluid phase at very high drives. We map the regions in which the thermally clogged, partially clogged, active uniform fluid, clustered fluid, active clogged, and directionally locked states occur as a function of disk density, drift force, and activity.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Implementation of High Temperature and Pressure fluid property interpolation tables

This document describes the use of module property_interpolate by John Doherty (Doherty,2006). The document “Fast Lookup of CO 2 Properties” describes the numerical details associated with lookup table. Rajeh Pawar modified property_interpolate for easier implementation in the FEHM code. George Zyvoloski modified the package to produce tables for water and air at very high temperatures and pressures. He also modified the auxiliary codes(drivers) to use newer datasets from the National Institute of Standards and Technology. We note here that within the files associated with FEHM there three property_interpolate modules. They are interpolate_2a.f90 (water), interpolate_2b.f90 (air) , and interpolate_2c.f90(CO 2 ). The interpolate_2c module is different than the earlier module produced by Rajeh Pawar in that it includes high temperature and pressure data. While FEHM has a number of fluid physics modules, this document will describe only the software and algorithms associated with water, water vapor, and heat (WH).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The role of plasma instabilities in relativistic radiation-mediated shocks: stability analysis and particle-in-cell simulations

ABSTRACT Relativistic radiation-mediated shocks are likely formed in prodigious cosmic explosions. The structure and emission of such shocks are regulated by copious production of electron–positron pairs inside the shock-transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the non-linear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced-transport equation for the particles via pitch-angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to non-adiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of non-thermal distributions, arising from the microturbulence, can affect the observed shock-breakout signal in ways unaccounted for by current single-fluid models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Spatiotemporal control of structure and dynamics in a polar active fluid

We apply optimal control theory to drive a polar active fluid into new behaviors: relocating asters, reorienting waves, and on-demand switching between states. This study reveals general principles to program active matter for useful functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Deterministic Calibration of MFiX-PIC, Part 1: Settling Bed

The Particle-in-cell (PIC) numerical approach for modeling granular solids in fluid flow has gained significant interest in recent years. Valued for its often shorter time-to-solution, the PIC formulation relies on modeling statistical groupings of particles called parcels in cooperation with a solids stress model to affect local solids velocity. This is in contrast to the discrete element model (DEM) where every particle in a system is modelled individually and directly coupled to local solids velocity through Newtonian mechanics. The U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) develops and maintains Multiphase Flow with Interphase eXchanges (MFiX), a collection of open-source computational fluid dynamics (CFD) solvers. Included in the MFiX suite are traditional two-fluid model (TFM) and DEM solvers, and a recently added PIC solver (NETL, 2021). In general, PIC methodologies offer an accuracy trade-off in lieu of computational speed; and therefore, it is important to assess the credibility of MFiX-PIC simulations. For this purpose, a systematic verification, validation and uncertainty quantification (VVUQ) effort was initiated at NETL to assess the new PIC solver

42 ENGINEERING↗
DOE OSTI2022

Hydrodynamization times of a holographic fluid far from equilibrium

Here we investigate several hydrodynamization times for an ensemble of different far-from-equilibrium solutions of the strongly coupled $\mathscr{N}$ = 4 supersymmetric Yang-Mills plasma undergoing Bjorken flow. For the ensemble of initial data analyzed in the present work, we find that, with typical tolerances between 3% to 5%, the average hydrodynamization time associated with the late time convergence of the pressure anisotropy to the corresponding Borel resummed hydrodynamic attractor is approximately equal to the average hydrodynamization time associated with the Navier-Stokes result, while both are shorter than the average hydrodynamization time associated with second-order hydrodynamics. On the other hand, we find that the entropy density of the different solutions coalesces to second-order hydrodynamics long before entering in the Navier-Stokes regime. A clear hierarchy between the different average hydrodynamization times of the Bjorken expanding fluid is established for the set of analyzed initial data, comprising also some solutions which, whilst satisfying the dominant and the weak energy conditions at the initial time, evolve such as to transiently violate one or both conditions when the fluid is still far from equilibrium. In particular, solutions violating the weak energy condition are generally found to take a longer time to enter in the hydrodynamic regime than the other solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Characterization of flowing liquid films as a regenerating plasma mirror for high repetition-rate laser contrast enhancement

In this paper, we characterize a high repetition-rate regenerating plasma mirror produced by the thin film of liquid formed when two laminar streams collide. The use of a flowing liquid film is inexpensive and the interaction surface refreshes automatically, avoiding buildup of on-target debris. The composition of the liquid material and the relative angle of the film-generating nozzles was optimized for this application. Spectra measured in reflection from a water-based plasma mirror showed a blue shift but an optical reflectivity of up to 30%. The thickness of the film was found to be of the order of 2 μm, and the stability of the reflected spot was ≈1 mrad. The reflected beam profile was highly distorted but stable. Further optimization of the nozzles to affect the fluid flow should enable significant improvements in control of the fluid films and increase in the reflectivity of these mirrors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

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

Chemical evolution in nitrogen shocked beyond the molecular stability limit

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. Finally, this suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects.

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

Facile polymerization in a bicellar template to produce polymer nano-rings

We report a well-defined discoidal bicelle composed of three lipids, specifically zwitterionic long-chain 1,2-dipalmitoyl phosphocholine (DPPC) and short-chain 1,2-dihexanoyl phosphocholine (DHPC) doped with anionic 1,2-dipalmitoyl phosphoglycerol (DPPG) provides a generalized template for the synthesis of hydrophobic polymer nano-rings. The lipid molar ratio of DPPC/DHPC/DPPG is 0.71/0.25/0.04. The detailed investigation and discussion were based on styrene but tested on three other vinyl monomers. The structure of nano-rings is identified through the detailed analysis of small angle X-ray/neutron scattering (SAXS and SANS) data and transmission electron micrographs (TEM), supported by the differential scanning calorimetric (DSC) data before and after polymerization. The investigation covers samples with a styrene-to-lipid ratio ranged varied from 1:50 to 1:10.The styrene monomers are initially located at both the discoidal planar (long-chain lipid rich) and rim (short-chain lipid rich) regions. During polymerization, they migrate to the more fluid rim region-section. The formation mechanism involves the interplay of hydrophobic interaction, mismatched miscibility of polystyrene between the ordered and disordered phases, and crystallinity of the long lipid acyl chains. This facile synthesis is proven applicable for several hydrophobic monomers. The well-defined nano-rings greatly enhance the interfacial area and have the potential to be the building blocks for functional materials, if monomers are incorporated with desirable functions, for future applications. Crosslinked polymer nano-rings are produced in well-defined discoidal bicelles. The enhanced immiscibility in the planar region drives the hydrophobic polymers to the fluidic rim. The platform can be generalized for synthesis of other hydrophobic polymer nano-rings.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Planar turbulent wakes under pressure gradient: Integral and self-similarity analyses

By using a combination of integral and self-similarity analyses, the generalized analytical solutions for the mean transverse velocity and Reynolds shear stress are rigorously derived for the first time for the far field of planar turbulent wakes under arbitrary pressure gradients. Specifically, by assuming self-similarity for the mean axial velocity, the analytical formulation for the mean transverse velocity is obtained from the integral of the mean continuity equation, and the analytical formulation for the Reynolds shear stress is obtained from the integral of the momentum equation. The generalized analytical formulations for the mean transverse velocity and Reynolds shear stress consist of multiple components, each with its unique scale and physical mechanism. In the zero pressure gradient limit, the generalized formulations recover the single-scale equations reported by Wei, Liu, and Livescu. Furthermore, simpler approximate formulations for the mean transverse velocity and Reynolds shear stress are also obtained, and show excellent agreement with the experimental measurements. Here, the findings provide new insights into the properties of planar turbulent wakes under pressure gradients, filling some long-standing gaps in the existing literature.

42 ENGINEERING↗
DOE OSTI2023

Ponderomotive recoil for electromagnetic waves

When waves damp or amplify on resonant particles in a plasma, nonresonant particles experience a recoil force that conserves the total momentum between particles and electromagnetic fields. This force is important to understand, as it can completely negate current drive and rotation drive mechanisms that are predicted on the basis of only resonant particles. Here, the existing electrostatic theory of this recoil force is extended to electromagnetic waves. While the result bears close similarity to historical fluid theories of laser–plasma interactions, it now incorporates both resonant and nonresonant particles, allowing momentum conservation to be self-consistently proven. Furthermore, the result is shown to be generally valid for kinetic plasmas, which is verified through single-particle hot-plasma simulations. Here, the new form of the force provides physical insight into the nature of the generalized Minkowski (plasmon) momentum of geometrical optics, which is shown to correspond to the momentum gained by the field and nonresonant particles as the wave is self-consistently ramped up from vanishing amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

3D polycatenated architected materials

Architected materials derive their properties from the geometric arrangement of their internal structural elements. Their designs rely on continuous networks of members to control the global mechanical behavior of the bulk. Here, in this study, we introduce a class of materials that consist of discrete concatenated rings or cage particles interlocked in three-dimensional networks, forming polycatenated architected materials (PAMs). We propose a general design framework that translates arbitrary crystalline networks into particle concatenations and geometries. In response to small external loads, PAMs behave like non-Newtonian fluids, showing both shear-thinning and shear-thickening responses, which can be controlled by their catenation topologies. At larger strains, PAMs behave like lattices and foams, with a nonlinear stress-strain relation. At microscale, we demonstrate that PAMs can change their shapes in response to applied electrostatic charges. The distinctive properties of PAMs pave the path for developing stimuli-responsive materials, energy-absorbing systems, and morphing architectures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

New systems in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) serves as a common library of classes between applications developed for advanced reactor analysis, fusion device engineering, spent fuel cask analysis, geochemistry studies, among other fields. These applications drive the development of the framework to meet their needs. Systems in MOOSE group capabilities that share a common purpose and generally common code. They can be leveraged by all downstream applications, providing extensive code re-use and shared maintenance. They facilitate the discovery by new users of the classes meeting at least partially their needs, and offer the same opportunities for customization as other systems. The addition of a new system to MOOSE opens new ways of solving or discretizing nonlinear problems, of performing distributed postprocessing, and a plethora of other needs. While new systems can be introduced in downstream applications rather than at the framework level, the framework team monitors common needs across the community and often triggers their addition. Documentation, training material, development needs can be centralized, limiting duplicated work across the community. The last three years have seen a large expansion in the capabilities of MOOSE. The supporting role of the framework in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has created numerous feature requests to support neutronics, thermal hydraulics, computational fluid dynamics and thermo-mechanics simulations in the Griffin, SAM, Pronghorn and Bison applications respectively. Similarly, laboratory-directed research and development (LDRD) projects in additive manufacturing, high-Reynolds flow simulations, structure optimization also necessitate an expansion of the framework capabilities. This summary reports on the new systems created in MOOSE, their design, their capabilities and some of the relevant interfaces.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

System and method for measuring fluid properties using magnetic field techniques via magnetic tracer

A method and apparatus for determining a property of a fluid in a vessel. The method uses magnetic tracer particles and an externally applied magnetic field which orients the particles. When the fluid moves, it changes the orientation of the tracer particles, thus changing the magnetic fields. These changes are detected by external magnetic field sensors. By using mathematical models, the property of the fluid in the vessel is determined from the detected magnetic field. In this manner, fluid vorticity, velocity, strain and stress may be estimated.

42 ENGINEERING↗