Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow physics”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Numerical Simulations of High-Speed Chemically Reacting Flow

The Essentially NonOscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Ton, V. T.↗
NASA NTRS1994

Numerical Simulations of High-Speed Chemically Reacting Flow

The essentially nonoscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Ton, V. T.↗
DOE OSTI2025

Measurement of directed flow in Au + Au collisions at $\sqrt{s{NN}}$ = 19.6 and 27 GeV with the STAR event plane detector

In heavy-ion collision experiments, the global collectivity of final-state particles can be quantified by anisotropic flow coefficients (𝑣 𝑛 ) . The first-order flow coefficient, also referred to as the directed flow (𝑣 1 ), describes the collective sideward motion of produced particles and nuclear fragments in heavy-ion collisions. It carries information on the very early stage of the collision, especially at large pseudorapidity (𝜂), where it is believed to be generated during the nuclear passage time. Directed flow therefore probes the onset of bulk collective dynamics during thermalization, providing valuable experimental guidance to models of the pre-equilibrium stage. In 2018, the Event Plane Detector (EPD) was installed in STAR and used for the Beam Energy Scan phase-II (BES-II) data taking. The combination of EPD (2.1 < |𝜂| < 5.1) and high-statistics BES-II data enables us to extend the 𝑣 1 measurement to the forward and backward 𝜂 regions. In this paper, we present the measurement of 𝑣 1 over a wide 𝜂 range in Au + Au collisions at $\sqrt{s_{NN}}$ = 19.6 and 27 GeV using the STAR EPD. The results of the analysis at $\sqrt{s_{NN}}$ = 19.6 GeV exhibit excellent consistency with the previous PHOBOS measurement, while elevating the precision of the overall measurement. The increased precision of the measurement also revealed finer structures in heavy-ion collisions, including a potential observation of the first-order event-plane decorrelation. Multiple physics models were compared to the experimental results. Only a transport model and a three-fluid hybrid model can reproduce a sizable 𝑣 1 at large 𝜂 as was observed experimentally. The model comparison also indicates 𝑣 1 at large 𝜂 might be sensitive to the QGP phase transition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Assessing the ultracentral flow puzzle in hydrodynamic modeling of heavy-ion collisions

An outstanding problem in heavy-ion collisions is the inability for models to accurately describe ultra-central experimental flow data, despite that being precisely the regime where a hydrodynamic description should be most applicable. Here, we reassess the status of this puzzle by computing the flow in ultra-central collisions obtained from multiple recent Bayesian models that were tuned to various observables in different collision systems at typical centralities. While central data can now be described with better accuracy than in previous calculations, tension with experimental observation remains and worsens as one goes to ultra-central collisions. Tuning the model parameters cannot remove this tension without destroying the fit at other centralities. As such, new elements are likely needed in the standard modeling of heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Observations of phase changes in monoolein during high viscous injection

Serial crystallography of membrane proteins often employs high-viscosity injectors (HVIs) to deliver micrometre-sized crystals to the X-ray beam. Typically, the carrier medium is a lipidic cubic phase (LCP) media, which can also be used to nucleate and grow the crystals. However, despite the fact that the LCP is widely used with HVIs, the potential impact of the injection process on the LCP structure has not been reported and hence is not yet well understood. The self-assembled structure of the LCP can be affected by pressure, dehydration and temperature changes, all of which occur during continuous flow injection. These changes to the LCP structure may in turn impact the results of X-ray diffraction measurements from membrane protein crystals. To investigate the influence of HVIs on the structure of the LCP we conducted a study of the phase changes in monoolein/water and monoolein/buffer mixtures during continuous flow injection, at both atmospheric pressure and under vacuum. The reservoir pressure in the HVI was tracked to determine if there is any correlation with the phase behaviour of the LCP. The results indicated that, even though the reservoir pressure underwent (at times) significant variation, this did not appear to correlate with observed phase changes in the sample stream or correspond to shifts in the LCP lattice parameter. During vacuum injection, there was a three-way coexistence of the gyroid cubic phase, diamond cubic phase and lamellar phase. During injection at atmospheric pressure, the coexistence of a cubic phase and lamellar phase in the monoolein/water mixtures was also observed. The degree to which the lamellar phase is formed was found to be strongly dependent on the co-flowing gas conditions used to stabilize the LCP stream. A combination of laboratory-based optical polarization microscopy and simulation studies was used to investigate these observations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Direct Analysis of Organic Compounds in Liquid Using a Miniature Photoionization Ion Trap Mass Spectrometer with Pulsed Carrier-Gas Capillary Inlet

A miniature ion trap mass spectrometer with capillary direct sampling and vacuum ultraviolet photoionization source was developed to conduct trace analysis of organic compounds in liquids. Self-aspiration sampling is available where the samples are drawn into the vacuum chamber through a capillary with an extremely low flow rate (less than 1 μL/min), which minimizes sample consumption in each analysis to tens of micrograms. A pulsed gas-assisted inlet was designed and optimized to promote sample transmission in the tube and facilitate the cooling of ions, thereby improving instrument sensitivity. A limit of detection of 2 ppb could be achieved for 2,4-dimethylaniline in a methanol solution. The sampling system described in the present study is specifically suitable for a miniature photoionization ion trap mass spectrometer that can perform rapid and online analysis for liquid samples. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Investigation of boundary-layer wind predictions during nocturnal low-level jet events using the Weather Research and Forecasting model

Simulations of two periods featuring three consecutive low level jet (LLJ) events in the US Upper Great Plains during the autumn of 2011 were conducted to explore the impacts of various setup configurations and physical process models on simulated flow parameters within the lowest 200 m above the surface, using the Weather Research and Forecasting (WRF) model. Sensitivities of simulated flow parameters to the horizontal and vertical grid spacing, planetary boundary layer (PBL) and land surface model (LSM) physics options, were assessed. Data from a Light Detection and Ranging (lidar) system, deployed to the Weather Forecast Improvement Project (WFIP; Finley et al. 2013) were used to evaluate the accuracy of simulated wind speed and direction at 80 m above the surface, as well as their vertical distributions between 120 and 40 m, covering the typical span of contemporary tall wind turbines. All of the simulations qualitatively captured the overall diurnal cycle of wind speed and stratification, producing LLJs during each overnight period, however large discrepancies occurred at certain times for each simulation in relation to the observations. 54-member ensembles encompassing changes of the above discussed configuration parameters displayed a wide range of simulated vertical distributions of wind speed and direction, and potential temperature, reflecting highly variable representations of stratification during the weakly stable overnight conditions. Root mean square error (RMSE) statistics show that different ensemble members performed better and worse in various simulated parameters at different times, with no clearly superior configuration . Simulations using a PBL parameterization designed specifically for the stable conditions investigated herein provided superior overall simulations of wind speed at 80 m, demonstrating the efficacy of targeting improvements of physical process models in areas of known deficiencies. However, the considerable magnitudes of the RMSE values of even the best performing simulations indicate ample opportunities for further improvements.

Mirocha, Jeff D.↗
DOE OSTI2022

Reactive Transport Modeling of Mineral Precipitation and Carbon Trapping in Discrete Fracture Networks

In this study we use numerical experiments to analyze reactive flow and transport behavior in discrete intersecting fracture networks, focusing on (a) how reaction-induced changes in physical and chemical properties affect flow connectivity and (b) how fracture networks developed in the Earth's critical zone contribute to carbon sequestration via mineral weathering reactions. In the first part of the study, we used two-dimensional reactive flow and transport simulations to analyze the impacts of mixing in a natural discrete fracture network. We concluded that reaction-induced changes can substantially alter the flow connectivity, especially at fracture intersections. The second set of simulations considered the problem of natural weathering of fractured mafic and ultramafic rocks in the partially saturated Earth's critical zone as a function of infiltration rates, fracture permeability, and partially saturated flow parameters. As a model system, we considered an incongruent reaction network with dissolution of forsterite and precipitation of magnesite. The behavior is complex in terms of the rate-controlling processes because of the multicomponent nature of the system as shown by the grid Peclet number: the CO 2 behavior is gas diffusion-controlled in the partially saturated zone, while the rate of water flow via the Damkӧhler number controls Mg 2+ transport through the fracture network. The amounts of carbon that can be trapped are modest, but the naturally fractured domain considered here provides a useful “base case” against which various engineered solutions can be compared.

58 GEOSCIENCES↗
DOE OSTI2018

Hg isotopes reveal in-stream processing and legacy inputs in East Fork Poplar Creek, Oak Ridge, Tennessee, USA

In this paper, natural abundance stable Hg isotope measurements were used to place new constraints on sources, transport, and transformations of Hg along the flow path of East Fork Poplar Creek (EFPC), a point-source contaminated headwater stream in Oak Ridge, Tennessee. Particulate-bound Hg in the water column of EFPC within the Y-12 National Security Complex, was isotopically similar to average metallic Hg(0) used in industry, having a mean δ 202 Hg value of -0.42 ± 0.09‰ (1SD) and near-zero Δ 199 Hg. On average, particulate fraction δ 202 Hg values increased downstream by 0.53‰, while Δ 199 Hg decreased by -0.10‰, converging with the Hg isotopic composition of the fine fraction of streambed sediment along the 26 km flow path. The dissolved fraction behaved differently. Although initial Δ 199 Hg values of the dissolved fraction were also near-zero, these values increased transiently along the flow path. Initial δ 202 Hg values of the dissolved fraction were more variable than in the particulate fraction, ranging from -0.44 to 0.18‰ among three seasonal sampling campaigns, but converged to an average δ 202 Hg value of 0.01 ± 0.10‰ (1SD) downstream. Dissolved Hg in the hyporheic and riparian pore water had higher and lower δ 202 Hg values, respectively, compared to dissolved Hg in stream water. Finally, variations in Hg isotopic composition of the dissolved and suspended fractions along the flow path suggest that: (1) physical processes such as dilution and sedimentation do not fully explain decreases in total mercury concentrations along the flow path; (2) in-stream processes include photochemical reduction, but microbial reduction is likely more dominant; and (3) additional sources of dissolved mercury inputs to EFPC at baseflow during this study predominantly arise from the hyporheic zone.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

Influence of zonal flow and density on resistive drift wave turbulent transport

The generations of zonal flow (ZF) and density (ZD) and their feedback on the resistive drift wave turbulent transport are investigated within the modified Hasegawa-Wakatani model. With proper normalization, the system is only controlled by an effective adiabatic parameter, a, where the ZF dominates the collisional drift wave (DW) turbulence in the adiabatic limit a > 1. By conducting direct numerical simulations, we found that the ZF can significantly reduce the transport by trapping the DWs in the vicinities of its extrema for a > 1, whereas the ZD itself has little impact on the turbulence but can only assist ZF to further decrease the transport by flattening the local plasma density gradient.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Molecular gas dynamics applied to low-thrust propulsion

The Direct Simulation Monte Carlo method is currently being applied to study flowfields of small thrusters, including both the internal nozzle and the external plume flow. The DSMC method is employed because of its inherent ability to capture nonequilibrium effects and proper boundary physics in low-density flow that are not readily obtained by continuum methods. Accurate prediction of both the internal and external nozzle flow is important in determining plume expansion which, in turn, bears directly on impingement and contamination effects.

Zelesnik, Donna↗
NASA NTRS1996

Simulations of Active Vortex Generators

We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. we consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise vorticity is generated and ejected due to the oscillatory motion of the lid. The present simulations complement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin (1994). The viscous effects are taken into account by modifying the strength Of the particles, whereas fast multipole schemes employing hundreds of thousands of particles allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies.

Mansour, N. N.↗
NASA NTRS1996

Simulations of Active Vortex Generators

We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. We consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise voracity is generated and ejected due to the oscillatory motion of the lid. The present simulations c Implement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University. Jacobson and Reynolds used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds. More recently, Lachowiez and Wlezien are investigating the flow generated by an electro-mechanically driven lid to be used for assertion control in aerodynamic applications. We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin. The viscous effects are taken into account by modifying the strength of the particles, whereas fast multipole schemes employing hundreds of thousands ol'particle's allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies. Estimates of the forces needed to drive the devices will also be presented.

Mansour, N. N.↗
DOE OSTI2023

Scaling approach to nuclear structure in high-energy heavy-ion collisions

In high-energy heavy-ion collisions, the initial condition of the produced quark-gluon plasma (QGP) and its evolution are sensitive to collective nuclear structure parameters describing the shape and radial profiles of the nuclei. Here, we find a general scaling relation between these parameters and many experimental observables such as elliptic flow, triangular flow, and particle multiplicity distribution. In particular, the ratios of observables between two isobar systems depend only on the differences of these parameters, but not on the details of the final state interactions, hence offering a new way to constrain the QGP initial condition. Using this scaling relation, we show how the structure parameters of $^{96}_{44}$Ru and $^{96}_{40}$Zr conspire to produce the rich centrality dependences of these ratios, as measured by the STAR Collaboration. Our scaling approach demonstrates that isobar collisions are a precision tool to probe the initial condition of heavy-ion collisions, as well as the collective nuclear structures, including the neutron skin, of the atomic nuclei across energy scales.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Probing triaxial deformation of atomic nuclei in high-energy heavy ion collisions

Most atomic nuclei are deformed with a quadrupole shape described by its overall strength β 2 and triaxiality γ. The deformation can be accessed in high-energy heavy ion collisions by measuring the collective flow response of the produced quark-gluon plasma to the eccentricity ε 2 and the density gradient d ⊥ in the initial state. Using an analytical estimate and a Glauber model, I show that the variances $\langleε^{2}_{2}\rangle$ or $\langle$(δd ⊥ /d ⊥ ) 2 $\rangle$ and skewnesses $\langleε^{2}_{2}δd_{⊥}/d_{⊥}\rangle$ or $\langle(δd_{⊥}/d_{⊥})^{3}\rangle$ have a simple analytical form of a'+b'β$^{2}_{2}$ and a'+[b'+c' cos(3γ)]β$^{3}_{2}$, respectively. From these, I constructed several normalized skewnesses to isolate the γ dependence from that of β 2 and show that the correlations between a normalized skewness and a variance can constrain simultaneously β 2 and γ. Assuming a linear relation with elliptic flow v 2 and mean-transverse momentum [p T ] of final-state particles, v 2 ∝ ε 2 and δ [p T ] / [p T ] ∝ δd ⊥ /d ⊥ , similar conclusions are also expected for the variances and skewnesses of v 2 and [p T ], i.e., a + bβ$^{2}_{2}$ for $\langle v^{2}_{2}\rangle$ and $\langle(δ[p_{T}] / [p_{T}])^{2}\rangle$ and a + [b + c cos (3γ)] β$^{3}_{2}$ for $\langle v^{2}_{2}δ [p_{T}] / [p_{T}]\rangle$ or $\langle(δ[p_{T}] / [p_{T}])^{3}\rangle$. My findings motivate a dedicated system scan of high-energy heavy ion collisions at RHIC and LHC to measure triaxiality of atomic nuclei: one first determines the coefficients b and c by collisions of isobaric near prolate nuclei, cos (3γ) ≈ 1, and near oblate nuclei, cos (3γ) ≈ –1 , with known β 2 values, followed by collisions of other species of interest with similar mass number. The (β 2 , γ) values for this species can be inferred directly from the measured variance and skewness observables from these collisions. Finally, the results demonstrate the unique opportunities offered by high-energy collisions as a tool to perform interdisciplinary nuclear physics studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Electric charge and strangeness-dependent directed flow splitting of produced quarks in Au+Au collisions

We report directed flow (v 1 ) of multistrange baryons (Ξ and Ω) and improved v1 data for $K¯$, $\bar{p}$, $\bar{Λ}$ and φ in Au+Au collisions at $\sqrt{s_{NN}} =27$ and 200 GeV from the STAR at the Relativistic Heavy Ion Collider (RHIC). We focus on particles whose constituent quarks are not transported from beam rapidity rather produced in the collisions. In midcentral collisions, we observe a coalescence sum rule for hadron combinations with identical quark content and a difference (“splitting”) in the slope of v 1 vs. rapidity for combinations having nonidentical quark content. The splitting strength appears to increase with the electric charge difference and strangeness content difference of the constituent quarks in the combinations, consistent with an electromagnetic effect. The peripheral collision statistics are insufficient to draw firm conclusions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Fluid dynamic mathematical aspects of supernova remnants

Supernovae—explosions of stars—are a central problem in astrophysics since they contain information on the entire process of stellar evolution and nucleosynthesis. Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities, developing during the supernova blast, lead to intense interfacial RT/RM mixing of the star's materials and couple astrophysical to atomic scales. This work analyzes some fluid dynamic mathematical aspects of the titanic task of supernova's blast. We handle mathematical challenges of RT/RM dynamics in supernova relevant conditions by directly linking the conservation laws governing RT/RM dynamics to symmetry-based momentum model, by exactly deriving the model parameters in the scale-dependent and scale-invariant regimes, and by exploring the special self-similar class for RT/RM interfacial mixing with variable accelerations. Here we reveal that RT/RM dynamics is strongly influenced by deterministic (the initial and the flow) conditions in the scale-dependent linear and nonlinear regimes and in the self-similar mixing regime. The theory outcomes are consistent with the observations of supernova remnants, explain the results of the scaled laboratory experiments in high energy density plasmas, and yield the design of future experiments for the accurate quantification of RT/RM dynamics in supernova relevant conditions. We find that from fluid dynamic mathematical perspectives, supernovae can be regarded as an astrophysical initial value problem. Along with the guidance of what explodes at microscopic scales, supernova remnants encapsulate information on the explosion hydrodynamics and the associated deterministic conditions at macroscopic scales. We urge such effects be considered in interpretations of the observational data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Efficient numerical schemes for viscoplastic avalanches. Part 1: The 1D case

This paper deals with the numerical resolution of a shallow water viscoplastic flow model. Viscoplastic materials are characterized by the existence of a yield stress: below a certain critical threshold in the imposed stress, there is no deformation and the material behaves like a rigid solid, but when that yield value is exceeded, the material flows like a fluid. In the context of avalanches, it means that after going down a slope, the material can stop and its free surface has a non-trivial shape, as opposed to the case of water (Newtonian fluid). The model involves variational inequalities associated with the yield threshold: finite-volume schemes are used together with duality methods (namely Augmented Lagrangian and Bermúdez–Moreno) to discretize the problem. To be able to accurately simulate the stopping behavior of the avalanche, new schemes need to be designed, involving the classical notion of well-balancing. In the present context, it needs to be extended to take into account the viscoplastic nature of the material as well as general bottoms with wet/dry fronts which are encountered in geophysical geometries. We derived such schemes and numerical experiments are presented to show their performances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗