Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compressible flows”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Compressible pairwise interaction extended point-particle model for force prediction of shock-particle bed interaction

We propose a pairwise influence framework for the complex unsteady compressible particle-laden flow problem by accounting for the scattered hydrodynamic waves emitting from neighboring particles in a Euler-Lagrange simulation. It has been observed from particle-resolved (PR) simulations of randomly dispersed particle beds under a loading shock that the compressible pseudoturbulence dominates the flow system even after the primary shock has passed, which causes fluctuations observed in the forces experienced by the particle. Moreover, the fact that each particle exists in the vicinity of a random arrangement of other particles modifies the time history of the drag force experienced by each particle during and after the passage of the shock. First, the scattering flow field due to an incoming shock interacting with a single sphere is constructed using an analysis of the flow in the acoustic limit. Then we examine the validity of the compressible Maxey-Riley-Gatignol force model by comparing the force prediction against a PR simulation of two interacting particles for various particle arrangements and incoming shock strength. Subsequently, the neighboring influences are stored as a library of maps that can be used readily in the calculation of the perturbation force. Lastly, the pairwise interaction assumption is evaluated by comparing the force predicted with the model with PR simulations of a randomly packed particle bed of 10% volume fraction for both water and air as the fluid medium for an incoming shock Mach number 1.22. With a considerably lower cost for the implementation of the model compared to PR simulations, it is verified that the model is reasonably accurate in pinpointing particles whose peak force is significantly larger or smaller than the mean drag but also to capture the prolonged fluctuations after the initial shock.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Electrostatic Relativistic Fluid Models of Electron Emission in a Warm Diode

A semi-analytic fluid model has been developed for characterizing relativistic electron emission across a warm diode gap. Here we demonstrate the use of this model in (i) verifying multi-fluid codes in modeling compressible relativistic electron flows (the EMPIRE-Fluid code is used as an example; see also Ref. 1), (ii) elucidating key physics mechanisms characterizing the influence of compressibility and relativistic injection speed of the electron flow, and (iii) characterizing the regimes over which a fluid model recovers physically reasonable solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

PythonFOAM: In-situ data analyses with OpenFOAM and Python

Here, we outline the development of a general-purpose Python-based data analysis tool for OpenFOAM. Our implementation relies on the construction of OpenFOAM applications that have bindings to data analysis libraries in Python. Double precision data in OpenFOAM is cast to a NumPy array using the NumPy C-API and Python modules may then be used for arbitrary data analysis and manipulation on flow-field information. We highlight how the proposed wrapper may be used for an in-situ online singular value decomposition (SVD) implemented in Python and accessed from the OpenFOAM solver PimpleFOAM. Here, 'in-situ' refers to a programming paradigm that allows for a concurrent computation of the data analysis on the same computational resources utilized for the partial differential equation solver. In addition, to demonstrate parallel deployments, we deploy a distributed SVD, which collects snapshot data across the ranks of a distributed simulation to compute the global left singular vectors. Crucially, both OpenFOAM and Python share the same message passing interface (MPI) communicator for this deployment which allows Python objects and functions to exchange NumPy arrays across ranks. Subsequently, we provide scaling assessments of this distributed SVD on multiple nodes of Intel Broadwell and KNL architectures for canonical test cases such as the large eddy simulations of a backward facing step and a channel flow at friction Reynolds number of 395. Finally, we demonstrate the deployment of a deep neural network for compressing the flow-field information using an autoencoder to demonstrate an ability to use state-of-the-art machine learning tools in the Python ecosystem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Observation of Hydrodynamic Flows in Imploding Fusion Plasmas on the National Ignition Facility [plus supplemental information]

Inertial confinement fusion implosions designed to have minimal fluid motion at peak compression often show significant linear flows in the laboratory, attributable per simulations to percent-level imbalances in the laser drive illumination symmetry. In this work, we present experimental results which intentionally varied the mode 1 drive imbalance by up to 4% to test hydrodynamic predictions of flows and the resultant imploded core asymmetries and performance, as measured by a combination of DT neutron spectroscopy and high-resolution x-ray core imaging. Neutron yields decrease by up to 50%, and anisotropic neutron Doppler broadening increases by 20%, in agreement with simulations. Furthermore, a tracer jet from the capsule fill-tube perturbation that is entrained by the hot-spot flow confirms the average flow speeds deduced from neutron spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗
DOE OSTI2023

Switchable system for high-salinity brine desalination and fractional precipitation

A saline feed stream flows into a liquid-liquid extraction system; and a volatile organic solvent flows through a main compressor. The compressed volatile organic solvent then flows through a solvent regenerator, which can be a heat exchanger or a combination of a vaporization device and a condenser, to cool the volatile organic solvent. The cooled volatile organic solvent in liquid phase then flows into the liquid-liquid extraction system, where the saline feed stream contacts the volatile organic solvent to selectively extract water from the saline feed stream into the volatile organic solvent, producing a concentrated brine and an organic-rich mixture of water and the volatile organic solvent. The organic-rich mixture flows from the liquid-liquid extraction system into the solvent regenerator, where the organic-rich mixture is heated to produce an organic-rich vapor and desalinated water; and the organic-rich vapor is recycled as volatile organic solvent back into the liquid-liquid extraction system.

Wilson, Aaron D.↗
DOE OSTI2020

Method for Comparing Porous Carbon Electrode Performance in Redox Flow Batteries

In this paper, we present a method for comparing porous electrodes for use in redox flow batteries. By comparing physical and electrochemical properties of six carbon cloth electrodes (CCEs) including hydraulic permeability, Ohmic resistance, and mass transport overpotential, we show how optimal engineering conditions such as compression and electrolyte flow rates can be derived. These results are combined and applied to a quinone-based redox flow battery, leading to a peak power density around 0.4 W cm –2 for several different CCEs. This method can be applied broadly to other electrode architectures and flow battery chemistries, and provides a comprehensive framework for their comparison.

25 ENERGY STORAGE↗
DOE OSTI2022

Variational encoder geostatistical analysis (VEGAS) with an application to large scale riverine bathymetry

Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sinuosity‐Driven Hyporheic Exchange: Hydrodynamics and Biogeochemical Potentials

Abstract Hydrologic exchange processes are critical for ecosystem services along river corridors. Meandering contributes to this exchange by driving channel water, solutes, and energy through the surrounding alluvium, a process called sinuosity‐driven hyporheic exchange. This exchange is embedded within and modulated by the regional groundwater flow (RGF), which compresses the hyporheic zone and potentially diminishes its overall impact. Quantifying the role of sinuosity‐driven hyporheic exchange at the reach‐to‐watershed scale requires a mechanistic understanding of the interplay between drivers (meander planform) and modulators (RGF) and its implications for biogeochemical transformations. Here, we use a 2D, vertically integrated numerical model for flow, transport, and reaction to analyze sinuosity‐driven hyporheic exchange systematically. Using this model, we propose a dimensionless framework to explore the role of meander planform and RGF in hydrodynamics and how they constrain nitrogen cycling. Our results highlight the importance of meander topology for water flow and age. We demonstrate how the meander neck induces a shielding effect that protects the hyporheic zone against RGF, imposing a physical constraint on biogeochemical transformations. Furthermore, we explore the conditions when a meander acts as a net nitrogen source or sink. This transition in the net biogeochemical potential is described by a handful of dimensionless physical and biogeochemical parameters that can be measured or constrained from literature and remote sensing. This work provides a new physically based model that quantifies sinuosity‐driven hyporheic exchange and biogeochemical reactions, a critical step toward their representation in water quality models and the design and assessment of river restoration strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Crankcase ventilation system with dead space alignment sleeves

A cylinder aligning sleeve and adapter plate for use within an internal combustion engine. The sleeve includes apertures placed contiguously with holes within the adapter plate to form a flow passage through the dead space between compression and crosshead pistons. The flow passage exits to a vent line where gas detection may occur along with safe venting of any gas detected. The sleeve serves to ensure concentric alignment of the compression cylinder with the engine block bore.

Alvarado, Caleb A.↗
DOE OSTI2020

Stress effects on flow and transport in three-dimensional fracture networks

We investigate the effects of various external stress regimes on fracture apertures, fluid flow, and solute transport in three-dimensional fracture networks. We use well-established geomechanics equations coupled with discrete fracture network modeling to characterize changes in primary flow paths within a complex network as a function of stress magnitude and orientation. These changes manifest in the alterations of the fluid flow field and are measured in terms of Eulerian and Lagrangian flow observables including solute transport, which is a key problem in many hydrologic applications. Changes in primary flow paths affect the solute transport in the network by promoting anomalously early arrival or long tailing behavior. However, early time arrival is not ubiquitous in anisotropically stressed networks and in most cases there is a delayed arrival of solute, which is attributed to (i) the presence of low-velocity zones normal to the flow direction, (ii) the angle between flow direction and major compressive principal stress directions, and (iii) changes in the primary flow paths (i.e., increases in tortuosity and active network structure). Overall, flows become more channelized in anisotropically stressed fracture networks than in unstressed and isotropically stressed networks.

58 GEOSCIENCES↗
DOE OSTI2021

Time-resolved particle image velocimetry

Time-resolved particle image velocimetry (TR-PIV) has become widespread in fluid dynamics. Essentially a velocity field movie, the dynamic content provides temporal as well as spatial information, in contrast to conventional PIV offering only statistical ensembles of flow quantities. From these time series arise further analyses such as accelerometry, space-time correlations, frequency spectra of turbulence including spatial variability, and derivation of pressure fields and forces. The historical development of TR-PIV is chronicled, culminating in an assessment of the current state of technology in high-repetition-rate lasers and high-speed cameras. Commercialization of pulse-burst lasers has expanded TR-PIV into more flows, including the compressible regime, and has achieved MHz rates. Particle response times and peak locking during image interrogation require attention but generally are not impediments to success. Accuracy considerations are discussed, including the risks of noise and aliasing in spectral content. Oversampled TR-PIV measurements allow use of multi-frame image interrogation methods, which improve the precision of the correlation and raise the velocity dynamic range of PIV. In combination with volumetric methods and data assimilation, a full four-dimensional description of a flow is not only achievable but becoming standardized. Lastly, a survey of exemplary applications is followed by a few predictions concerning the future of TR-PIV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

A new re-redistribution scheme for weighted state redistribution with adaptive mesh refinement

State redistribution (SRD) is a recently developed technique for stabilizing cut cells that result from finite-volume embedded boundary methods. SRD has been successfully applied to a variety of compressible and incompressible flow problems. When used in conjunction with adaptive mesh refinement (AMR), additional steps are needed to preserve the accuracy and conservation properties of the solution if the embedded boundary is not restricted to a single level of the mesh hierarchy. In this work, we extend the weighted state redistribution algorithm to cases where cut cells live at or near a coarse-fine interface within the domain. Here, we present numerical results that demonstrate that the algorithm is conservative when the coarse-fine interface intersects the embedded boundary. Additionally we compare the numerical solution of the Sod shock tube problem in an inclined cylinder with the analytic solution, and we compare the simulation of a shock hitting a cylindrical obstacle with experimental data. Finally we demonstrate the methodology for simulation of the multicomponent compressible Navier-Stokes equations in a piston-bowl geometry, and discuss the computational efficiency gained by not requiring the entire embedded boundary to be defined at the finest level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

System and method for rapid, high throughput, high pressure synthesis of materials from a liquid precursor

The present disclosure relates to a system and method for synthesis of condensed nano-materials to at least one of create nanoparticles or modify existing nanoparticles. In one embodiment the system may have a source of liquid precursor, with the liquid precursor including a compound therein. A flow control element and a compression wave generating subsystem are also included. The flow control element is in communication with the source of the liquid precursor and creates a jet of liquid precursor. The compression wave generating subsystem drives a compression wave through at least a substantial portion of a thickness of the jet of liquid precursor to sufficiently compress the jet of liquid precursor, and to increase pressure and temperature of the jet of liquid precursor, to at least one of create nanoparticles or modify existing nanoparticles.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Scalability of OpenFOAM Density-Based Solver with Runge–Kutta Temporal Discretization Scheme

Compressible density-based solvers are widely used in OpenFOAM, and the parallel scalability of these solvers is crucial for large-scale simulations. In this paper, we report our experiences with the scalability of OpenFOAM’s native rhoCentralFoam solver, and by making a small number of modifications to it, we show the degree to which the scalability of the solver can be improved. The main modification made is to replace the first-order accurate Euler scheme in rhoCentralFoam with a third-order accurate, four-stage Runge-Kutta or RK4 scheme for the time integration. The scaling test we used is the transonic flow over the ONERA M6 wing. This is a common validation test for compressible flows solvers in aerospace and other engineering applications. Numerical experiments show that our modified solver, referred to as rhoCentralRK4Foam, for the same spatial discretization, achieves as much as a 123.2% improvement in scalability over the rhoCentralFoam solver. As expected, the better time resolution of the Runge–Kutta scheme makes it more suitable for unsteady problems such as the Taylor–Green vortex decay where the new solver showed a 50% decrease in the overall time-to-solution compared to rhoCentralFoam to get to the final solution with the same numerical accuracy. Finally, the improved scalability can be traced to the improvement of the computation to communication ratio obtained by substituting the RK4 scheme in place of the Euler scheme. All numerical tests were conducted on a Cray XC40 parallel system, Theta, at Argonne National Laboratory.

Li, Sibo↗
DOE OSTI2024

Internal measurements of electromagnetic geodesic acoustic mode (GAM) in EAST plasmas

Velocity, density, and magnetic fluctuations of the geodesic acoustic mode (GAM) have been measured using the Doppler backscattering system, Faraday-effect polarimeter-interferometer, and external pick-up coils in the Experimental Advanced Superconducting Tokamak. Simultaneous measurements of density and velocity fluctuations at the midplane and top of plasmas demonstrate that m = 1 density fluctuations are quantitatively balanced by the compression of perpendicular flow fluctuations. Furthermore, internal magnetic fluctuations associated with GAM have now been directly measured by laser-based Faraday-effect polarimetry for the first time. Line-averaged magnetic fluctuations (up to 16 Gauss, B̃¯R,GAMBT∼0.066%) are significantly larger than those extrapolated from edge coils (a few Gauss) and that magnetic fluctuations increase with β. The observed discrepancy between finite β theory and experimental data indicates the need for further theoretical investigations.

Physics↗
DOE OSTI2021

A multispecies turbulence model for the mixing and de-mixing of miscible fluids

A de-mix capable model for turbulence in compressible, variable density flows is proposed. The model is based on the Besnard-Harlow-Rauenzahn (BHR) family of models (Besnard D, Harlow F, Rauenzahn R, et al. Turbulence transport equations for variable-density turbulence and their relationship to two-field models. NM (United States): Los Alamos National Laboratory; 1992 (Technical Report LA-12303-MS), but is extended to track the evolution of the turbulent fluxes and fluctuations in the material mass fractions for each species present. The new evolution equations are introduced without requiring additional closures or new empirically tuned coefficients relative to previous BHR models, and are shown to improve the model’s ability to reproduce the behaviour of simulations containing mixing layers that are constrained by a stabilising force. As a result, the model is tested in a range of canonical flows including Rayleigh-Taylor driven, shock driven, and shear driven turbulence, and is shown to produce reasonable agreement with simulations and experiments in these scenarios.

42 ENGINEERING↗