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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2013

Radial Profiles of Plasma Electron Characteristics in a Low-Power Arcjet

Since 1960, the Arc Jet Complex at NASA Ames Research Center has been a source of long-duration, high-enthalpy flow for materials testing with application to the thermal protection of aerospace vehicle components. From their inception the facilities have played an integral role supporting many of NASA's space flight programs and numerous DoD projects. In recent years advancements in computational fluid dynamics (CFO) have made the resultant models a valuable tool for assessing and predicting performance, however, the inherent limitation of models to compensate for a dissociated, transitionally high temperature nonequilibrated flowfield have made further measurements necessary. The use of electrostatic probe diagnostics within similarly harsh plasma environments in previous studies have been met with much success. In this study, the use of a single Langmuir probe was implemented to characterize the plasma parameters of interest as they vary radially within a large volume of the plume. Classical Langmuir probe theory was applied to achieve first order estimates of the heavy particle temperature, the ratio T(sub i)/T(sub e), and the ionization fraction. As expected, both the electron temperature and electron density measurements show a dependence on radial distance from the plume centerline, with electron density profiles showing the largest dependence. This paper aims to validate and strengthen modeling work conducted at the NASA Ames Research Center by measuring the critical plasma parameters (electron temperature, electron density, and plasma potential) within the arc plume of a subscale arc jet. These parameters are intended to give physical insight into the flow characteristics while providing the necessary boundary conditions to validate full scale simulations.

Codron, Douglas↗
DOE OSTI2021

Statistical uncertainties of the v n { 2 k } harmonics from Q cumulants

Analytic formulas to calculate statistical uncertainties of v n { 2 k } harmonics extracted from the Q cumulants are presented. The Q cumulants are multivariate polynomial functions of the weighted means of 2m-particle azimuthal correlations, $\langle\langle 2 m \rangle\rangle$. Variances and covariances of the $\langle\langle 2 m \rangle\rangle$ are included in the analytic formulas of the uncertainties that can be calculated simultaneously with the calculations of the v n { 2 k } harmonics. The calculations are performed using a simple toy model, which roughly simulates elliptic flow azimuthal anisotropy with magnitudes around 0.05. The results are compared with the results obtained by the many data subsets, and by the bootstrapping method. The first one is a common way of estimation of the statistical uncertainties of the v n { 2 k } harmonics in a real experiment. In order to increase precision in the measurement of the v n { 2 k } harmonics, a large number of 15 000 subsets and the same number of the resampling in the bootstrap method is used. Unlike the other ways of the analytic calculation of the v n { 2 k } statistical uncertainties, our proposal that includes the use of squared weights in the calculation of both the variances and covariances, gives the best agreement with the results obtained from the subsets and bootstrap method. Additionally, a recurrence equation between Q cumulants of any order is also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Plasma Imaging, LOcal Measurement, and Tomographic Experiment (PILOT): A Mission Concept for Transformational Multi-Scale Observations of Mass and Energy Flow Dynamics in Earth’s Magnetosphere

We currently do not understand the fundamental physical processes that govern mass and energy flow through the Earth’s magnetosphere. Knowledge of these processes is critical to understanding the mass loss rate of Earth’s atmosphere, as well as for determining the role that a planetary magnetic field plays in atmospheric retention, and therefore habitability, for Earth-like planets beyond the solar system. Mass and energy flow processes are challenging to determine at Earth in part because Earth’s planetary magnetic field creates a complex “system of systems” composed of interdependent plasma populations and overlapping spatial regions that perpetually exchange mass and energy across a broad range of temporal and spatial scales. Further, the primary mass carrier in the magnetosphere is cold plasma (as cold as ~0.1 eV), which is invisible to many space-borne instruments that operate in the inner magnetosphere. The Plasma Imaging LOcal and Tomographic experiment (PILOT) mission concept, described here, provides the transformational multi-scale observations required to answer fundamental open questions about mass and energy flow dynamics in the Earth’s magnetosphere. PILOT uses a constellation of spacecraft to make radio tomographic, remote sensing, and in-situ measurements simultaneously, fully capturing cold plasma mass dynamics and its impact on magnetospheric systems over an unprecedented range of spatial and temporal scales. This article details the scientific motivation for the PILOT mission concept as well as a potential mission implementation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Treating network junctions in finite volume solution of transient gas flow models

A finite volume scheme for the numerical solution of a non-isothermal non-adiabatic compressible flow model for gas transportation networks on non-flat topography is introduced. Unlike standard Euler equations, the model takes into account wall friction, variable height and heat transfer between the pipe and the environment which are source terms. The case of one single pipe was considered in a previous reference by the authors, , where a finite volume method with upwind discretization of the flux and source terms has been proposed in order to get a well-balanced scheme. The main goal of the present paper is to go a step further by considering a network of pipes. The main issue is the treatment of junctions for which container-like 2D finite volumes are introduced. The couplings between pipes (1D) and containers (2D) are carefully described and the conservation properties are analyzed. Numerical tests including real gas networks are solved showing the performance of the proposed methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Imaging nuclear shape through anisotropic and radial flow in high-energy heavy-ion collisions

Most atomic nuclei exhibit ellipsoidal shapes characterized by quadrupole deformation β 2 and triaxiality γ, and sometimes even a pear-like octupole deformation β 3 . The STAR experiment introduced a new ‘imaging-by-smashing’ technique to image the nuclear global shape by colliding nuclei at ultra-relativistic speeds and analyzing outgoing debris. Features of nuclear shape manifest in collective observables like anisotropic flow vn and radial flow via mean transverse momentum [p T ]. We present new measurements of the variances of v n (n = 2, 3, and 4) and [p T ], and the covariance of $v$$^{2}_{n}$ with [p T ], in collisions of highly deformed 238 U and nearly spherical 197 Au. Ratios of these observables between the two systems effectively suppress common final-state effects, isolating the strong impact of uranium’s deformation. By comparing results with state-of-the-art hydrodynamic model calculations, we extract β 2U and γ U values consistent with those deduced from low-energy nuclear structure measurements. Measurements of v 3 and its correlation with [p T ] also provide the first experimental suggestion of a possible octupole deformation for 238 U. These findings provide significant support for using high-energy collisions to explore nuclear shapes on femtosecond timescales, with implications for both nuclear structure and quark-gluon plasma studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A physics-informed and hierarchically regularized data-driven model for predicting fluid flow through porous media

This paper presents a new deep learning data-driven model for predicting structure dependent pore-fluid velocity fields in rock. The model is based on a Convolutional Auto-Encoder (CAE) artificial neural network capable of learning from image data generated by direct numerical simulations of fluid flow through pore-structures, such as by Lattice Boltzmann or molecular dynamics methods. The main novelty of the model in comparison to previous CAE-based data-driven approaches consists of three parts. The first is a methodology for decomposing the full-domain of the porous media into sub-regions, or “sub-domains”, in order to reduce the overall size of the CAE, batch process the sub-domains in parallel, and enable the CAE to learn local and generalizable nonlinear mappings of pore-fluid velocities. The second consists of embedding the finite difference solutions of the incompressible Navier-Stokes and continuity equations into convolutional layers prior to the CAE in order to provide the CAE with knowledge of fluid dynamics physics (PhyFlow). The third main novelty is that the training of the CAE is regularized with a hierarchical loss function that encourages the learning of fluid flow patterns (in a way similar to ranked modes in principal component analysis), ranking from most to least important. This is shown to increase the stability in learning, reduce over-fitting, and promote interpretability of the CAE neural network layers (HierCAE). The comprehensive new data-driven model, which we call the PhyFlow-HierCAE model, is shown to exhibit improved accuracy and generalizability of flow field predictions over conventional CAE models, attributable to the embedded physical knowledge and the hierarchical regularization, as well as realize orders of magnitude speed-ups in computation times as a surrogate for the direct numerical simulations. Examples of training and forward predictions on unseen pore-structures are provided and evaluated for data from Lattice Boltzmann and molecular dynamics simulations of pore-fluid flow. The model is shown to be a fast and accurate emulator (or “surrogate”) for predicting effective permeability of unseen pore-structures based on learning from relatively small direct numerical simulation datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

The dynamics of the Venus ionosphere

Data from the Pioneer-Venus orbiter has demonstrated the importance of understanding ion dynamics in the Venus ionosphere. The analysis of the data has shown that during solar maximum the topside Venus ionosphere in the dark hemisphere is generated almost entirely on the dayside of the planet during solar maximum, and flows with supersonic velocities across the terminator into the nightside. The flow field in the ionosphere is mainly axially-symmetric about the sun-Venus axis, as are most measured ionospheric quantities. The primary data base used consisted of the ion velocity measurements made by the RPA during three years that periapsis of the orbiter was maintained in the Venus ionosphere. Examples of ion velocities were published and modeled. This research examined the planetary flow patterns measured in the Venus ionosphere, and the physical implications of departures from the mean flow.

Miller, K. L.↗
NASA NTRS1992

Using electron scattering to probe the environment of cluster cooling flows

The electron-scattered surface brightness profiles are calculated for a strong central source embedded in a cluster cooling flow, and the profiles are used to discuss the physical state of the gas in the flow and the behavior of the central source. The effects of temporal variations in the luminosity of the central source are calculated for a burst model and a source whose luminosity increases linearly with time. If the emission from the central AGN is polarized, the surface brightness profiles have cusps along the direction of the projection of the normal to the polarization vector on the sky. The 2D surface brightness profiles produced by beamed radiation are calculated.

Wise, Michael W.↗
NASA NTRS1998

Models of Plumes: Their Flow, Their Geometric Spreading, and Their Mixing with Interplume Flow

There are two types of plume flow models: (1) 1D models using ad hoc spreading functions, f(r); (2) MagnetoHydroDynamics (MHD) models. 1D models can be multifluid, time dependent, and incorporate very general descriptions of the energetics. They confirm empirical results that plume flow is slow relative to requirements for high speed wind. But, no published 1 D model incorporates the rapid local spreading at the base (fl(r)) which has an important effect on mass flux. The one published MHD model is isothermal, but confirms that if b=8*pi*p/absolute value(B)2<<l then the field is nearly potential below -70,000 km. Building on the MHD result, we apply a two scale approximation to calculate fl(r). We also compute the global spreading (fg(r)) out to 5.0 RSUN imposed by coronal hole geometry. Global MHD models provide a potent method of calculating fg(r). Unambiguous plume signatures have not yet been found in the solar wind. This is probably due to strong mixing of plume and interplume flows near the Sun. We describe a physical source for strong mixing due to the observed flows being unstable to shear instabilities that lead to rapid disruption.

Suess, Steven T.↗
NASA NTRS2012

Advances in Rotor Performance and Turbulent Wake Simulation Using DES and Adaptive Mesh Refinement

Time-dependent Navier-Stokes simulations have been carried out for a rigid V22 rotor in hover, and a flexible UH-60A rotor in forward flight. Emphasis is placed on understanding and characterizing the effects of high-order spatial differencing, grid resolution, and Spalart-Allmaras (SA) detached eddy simulation (DES) in predicting the rotor figure of merit (FM) and resolving the turbulent rotor wake. The FM was accurately predicted within experimental error using SA-DES. Moreover, a new adaptive mesh refinement (AMR) procedure revealed a complex and more realistic turbulent rotor wake, including the formation of turbulent structures resembling vortical worms. Time-dependent flow visualization played a crucial role in understanding the physical mechanisms involved in these complex viscous flows. The predicted vortex core growth with wake age was in good agreement with experiment. High-resolution wakes for the UH-60A in forward flight exhibited complex turbulent interactions and turbulent worms, similar to the V22. The normal force and pitching moment coefficients were in good agreement with flight-test data.

Chaderjian, Neal M.↗
DOE OSTI2016

Mixing model with multi-particle interactions for Lagrangian simulations of turbulent mixing

We report on the numerical study of the mixing volume model (MVM) for molecular diffusion in Lagrangian simulations of turbulent mixing problems. The MVM is based on the multi-particle interaction in a finite volume (mixing volume). A priori test of the MVM, based on the direct numerical simulations of planar jets, is conducted in the turbulent region and the interfacial layer between the turbulent and non-turbulent fluids. The results show that the MVM predicts well the mean effects of the molecular diffusion under various numerical and flow parameters. The number of the mixing particles should be large for predicting a value of the molecular diffusion term positively correlated to the exact value. The size of the mixing volume relative to the Kolmogorov scale η is important in the performance of the MVM. The scalar transfer across the turbulent/non-turbulent interface is well captured by the MVM especially with the small mixing volume. Furthermore, the MVM with multiple mixing particles is tested in the hybrid implicit large-eddy-simulation/Lagrangian-particle-simulation (LES–LPS) of the planar jet with the characteristic length of the mixing volume of O(100η). Despite the large mixing volume, the MVM works well and decays the scalar variance in a rate close to the reference LES. The statistics in the LPS are very robust to the number of the particles used in the simulations and the computational grid size of the LES. Both in the turbulent core region and the intermittent region, the LPS predicts a scalar field well correlated to the LES.

42 ENGINEERING↗
DOE OSTI2021

Drift waves enstrophy, zonal flow, and nonlinear evolution of the modulational instability

The interaction of the drift wave (DW) turbulence and zonal flow (ZF) is investigated with the modified Hasegawa–Mima equation taking into account the backreaction of ZF velocity on DW turbulence. It is shown that the y-averaged enstrophy of DW turbulence and the velocity of ZF are intrinsically related. By utilizing this feature, a nonlinear stage of DW modulational instability is considered within the framework of the wave kinetic equation. It is shown that in this approximation, the nonlinear stage of the modulational instability results in the collapsing solutions, accompanied by the “wave breaking” phenomenon. Numerical simulations based on the Hasegawa–Mima equation show that for a weak DW turbulence, Φ ̃ = ( e φ ̃ / T e ) ( L n / ρ s ) ⪝ 1 , the collapsing-like features on both ZF and y-averaged enstrophy of DW turbulence decay in time and then re-emerge again at different locations. For the case of a strong DW turbulence, Φ ̃ > 1 , where nonlinear interactions of DW harmonics dominate, stable spatial structures of ZF and y-averaged enstrophy of DW turbulence emerge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Anisotropic Flow in Fixed-Target 208 Pb + 20 Ne Collisions as a Probe of Quark-Gluon Plasma

The System for Measuring Overlap with Gas (SMOG2) at the LHCb detector enables the study of fixed-target ion-ion collisions at relativistic energies ($\sqrt{𝑠_{NN}}$ ∼ 100 GeV in the center of mass). Here, with input from ab initio calculations of the structure of 16 O and 20 Ne , we compute 3+1⁢D hydrodynamic predictions for the anisotropic flow of Pb+Ne and Pb+O collisions to be tested with upcoming LHCb data. This will allow the detailed study of quark-gluon plasma formation as well as experimental tests of the predicted nuclear shapes. Elliptic flow (𝑣 2 ) in Pb + Ne collisions is greatly enhanced compared to the Pb + O baseline due to the shape of 20 Ne , which is deformed in a bowling-pin geometry. Owing to the large 208 Pb radius, this effect is seen in a broad centrality range, a unique feature of this collision configuration. Larger elliptic flow further enhances the quadrangular flow (𝑣 4 ) of Pb + Ne collisions via nonlinear coupling, and impacts the sign of the kurtosis of the elliptic flow vector distribution (𝑐 2 ⁡{4}). Exploiting the shape of 20 Ne proves thus an ideal method to investigate the formation of quark-gluon plasma in fixed-target experiments at LHCb, and demonstrates the power of System for Measuring Overlap with Gas as a tool to image nuclear ground states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Durability test of a flowing-water target for isotope harvesting

A high intensity proton irradiation was performed with the flowing-water isotope harvesting target at the University of Wisconsin-Madison Cyclotron Laboratory to measure the rate of degradation of the target shell during irradiation conditions. The beam reached an intensity of 34 µA by the end of the irradiation and covered an area of 0.7 cm 2 on the target. Radiolysis products, such as H 2 O 2 , H 2 , and O 2 , were measured in the bulk water of the system and found to be present at much lower levels than predicted by literature escape yields. Radionuclides formed in the target shell were measured in the system water as a radiotracer for target degradation. Using a simple, beam intensity dependent model, a corrosion rate of 1.5E-6 μm/(μA*s) was found to match the measured radiotracer activities at various points in the irradiation. Furthermore, this rate was used to extrapolate the lifetime of future isotope harvesting targets at the NSCL and FRIB, using the areal power density of different ion beams to scale the corrosion rate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

A cut-cell finite volume – finite element coupling approach for fluid–structure interaction in compressible flow

We present a loosely coupled approach for the solution of fluid–structure interaction problems between a compressible flow and a deformable structure. The method is based on staggered Dirichlet–Neumann partitioning. The interface motion in the Eulerian frame is accounted for by a conservative cut-cell Immersed Boundary method. The present approach enables sub-cell resolution by considering individual cut-elements within a single fluid cell, which guarantees an accurate representation of the time-varying solid interface. The cut-cell procedure inevitably leads to non-matching interfaces, demanding for a special treatment. A Mortar method is chosen in order to obtain a conservative and consistent load transfer. We validate our method by investigating two-dimensional test cases comprising a shock-loaded rigid cylinder and a deformable panel. Moreover, the aeroelastic instability of a thin plate structure is studied with a focus on the prediction of flutter onset. Finally, we propose a three-dimensional fluid–structure interaction test case of a flexible inflated thin shell interacting with a shock wave involving large and complex structural deformations.

42 ENGINEERING↗