SEARCH · Engineering Papers

Results for “flow conservation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Electromagnetically driven relativistic jets - A class of self-consistent numerical solutions

We explore the generalized Grad-Shafranov (GS) equation (Lovelace et al. 1986), which describes the cross-field balance of the magnetic flux surfaces in a stationary, axisymmetric, cold relativistic MHD wind. We construct a family of self-consistent, jet-type, non-self-similar solutions in which flux surfaces thread the equatorial plane vertically and are eventually collimated into nested cylinders. These solutions conserve the total energy and angular momentum along each flux surface and satisfy the relativistic GS equation in the cross-field direction. We find that the total specific energy in the flow depends primarily on a dimensionless combination of the poloidal magnetic flux, the rotation, and the mass loading rate, which reduces to Michel's (1969) magnetization parameter for radial winds. The final width of the flow and its ratio of kinetic energy flux to Poynting flux depend sensitively on the pressure distribution of the ambient medium, which is required to contain the jet at its outer edge when no jet-confirming external toroidal magnetic field is present. Furthermore, we show that the magnetic pressure in the flow can be much higher near the axis than close to the edge, reflecting the 'pinching' effect of the toroidal magnetic field.

Li, Zhi-Yun↗

NASA NTRS2000

Aerothermodynamic Analysis of Stardust Sample Return Capsule with Coupled Radiation and Ablation

An aerothermodynamic analysis of the forebody aeroshell of the Stardust Sample Return Capsule is carried out by using the axisymmetric viscous shock-layer equations with and without fully coupled radiation and ablation. Formulation of the viscous shock-layer equations with shoulder radius as the length scale and implementation of the Vigneron pressure condition allow resolution of the flowfield over the shoulder. With a predominantly supersonic outflow over the shoulder, a globally iterated solution or viscous shock-layer equations can be obtained. The stagnation-point results are obtained along a specified trajectory, whereas detailed calculations along the body are provided at the peak-heating point. The equilibrium calculations with ablation injection are the focus of the present study because of the lack of a general chemical nonequilibrium analysis that accounts for both surface and flowfield effect. The equilibrium calculations also provide a simple way to conserve surface (and flowfield) elemental composition for the current small ablation injection rates, where the surface elemental composition is a mixture of freestream and ablator elements. Therefore, the coupled laminar and turbulent flow solutions with radiation and ablation are obtained by using the equilibrium flow chemistry, whereas a nonequilibrium chemistry model is used for solutions without ablation and turbulence. Various computed results are compared with those obtained by the other researchers.

Gupta, Roop N.↗

NASA NTRS2006

A Computational and Experimental Investigation of Shear Coaxial Jet Atomization

The instability and subsequent atomization of a viscous liquid jet emanated into a high-pressure gaseous surrounding is studied both computationally and experimentally. Liquid water issued into nitrogen gas at elevated pressures is used to simulate the flow conditions in a coaxial shear injector element relevant to liquid propellant rocket engines. The theoretical analysis is based on a simplified mathematical formulation of the continuity and momentum equations in their conservative form. Numerical solutions of the governing equations subject to appropriate initial and boundary conditions are obtained via a robust finite difference scheme. The computations yield real-time evolution and subsequent breakup characteristics of the liquid jet. The experimental investigation utilizes a digital imaging technique to measure resultant drop sizes. Data were collected for liquid Reynolds number between 2,500 and 25,000, aerodynamic Weber number range of 50-500 and ambient gas pressures from 150 to 1200 psia. Comparison of the model predictions and experimental data for drop sizes at gas pressures of 150 and 300 psia reveal satisfactory agreement particularly for lower values of investigated Weber number. The present model is intended as a component of a practical tool to facilitate design and optimization of coaxial shear atomizers.

Ibrahim, Essam A.↗

NASA NTRS2014

Icing Analysis of a Swept NACA 0012 Wing Using LEWICE3D Version 3.48

Icing calculations were performed for a NACA 0012 swept wing tip using LEWICE3D Version 3.48 coupled with the ANSYS CFX flow solver. The calculated ice shapes were compared to experimental data generated in the NASA Glenn Icing Research Tunnel (IRT). The IRT tests were designed to test the performance of the LEWICE3D ice void density model which was developed to improve the prediction of swept wing ice shapes. Icing tests were performed for a range of temperatures at two different droplet inertia parameters and two different sweep angles. The predicted mass agreed well with the experiment with an average difference of 12%. The LEWICE3D ice void density model under-predicted void density by an average of 30% for the large inertia parameter cases and by 63% for the small inertia parameter cases. This under-prediction in void density resulted in an over-prediction of ice area by an average of 115%. The LEWICE3D ice void density model produced a larger average area difference with experiment than the standard LEWICE density model, which doesn't account for the voids in the swept wing ice shape, (115% and 75% respectively) but it produced ice shapes which were deemed more appropriate because they were conservative (larger than experiment). Major contributors to the overly conservative ice shape predictions were deficiencies in the leading edge heat transfer and the sensitivity of the void ice density model to the particle inertia parameter. The scallop features present on the ice shapes were thought to generate interstitial flow and horse shoe vortices which enhance the leading edge heat transfer. A set of changes to improve the leading edge heat transfer and the void density model were tested. The changes improved the ice shape predictions considerably. More work needs to be done to evaluate the performance of these modifications for a wider range of geometries and icing conditions

Ice Accretion Software↗

NASA NTRS2014

Icing Analysis of a Swept NACA 0012 Wing Using LEWICE3D Version 3.48

Icing CFD↗

DOE OSTI2021

Divergent species-specific impacts of whole ecosystem warming and elevated CO 2 on vegetation water relations in an ombrotrophic peatland

Boreal peatland forests have relatively low species diversity and thus impacts of climate change on one or more dominant species could shift ecosystem function. Despite abundant soil water availability, shallowly rooted vascular plants within peatlands may not be able to meet foliar demand for water under drought or heat events that increase vapor pressure deficits while reducing near surface water availability, although concurrent increases in atmospheric CO 2 could buffer resultant hydraulic stress. We assessed plant water relations of co-occurring shrub (primarily Rhododendron groenlandicum and Chamaedaphne calyculata) and tree (Picea mariana and Larix laricina) species prior to, and in response to whole ecosystem warming (0 to +9°C) and elevated CO 2 using 12.8-m diameter open-top enclosures installed within an ombrotrophic bog. Water relations (water potential [Ψ], turgor loss point, foliar and root hydraulic conductivity) were assessed prior to treatment initiation, then Ψ and peak sap flow (trees only) assessed after 1 or 2 years of treatments. Under the higher temperature treatments, L. laricina Ψ exceeded its turgor loss point, increased its peak sap flow, and was not able to recover Ψ overnight. In contrast, P. mariana operated below its turgor loss point and maintained constant Ψ and sap flow across warming treatments. Similarly, C. calyculata Ψ stress increased with temperature while R. groenlandicum Ψ remained at pretreatment levels. The more anisohydric behavior of L. laricina and C. calyculata may provide greater net C uptake with warming, while the more conservative P. mariana and R. groenlandicum maintained greater hydraulic safety. These latter species also responded to elevated CO 2 by reduced Ψ stress, which may also help limit hydraulic failure during periods of extreme drought or heat in the future. Along with Sphagnum moss, the species-specific responses of peatland vascular communities to drier or hotter conditions will shape boreal peatland composition and function in the future.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Analytical closure to the spatially-filtered Euler equations for shock-dominated flows

To ensure numerical stability in the vicinity of shocks, a variety of methods have been used, including shock-capturing schemes such as weighted essentially non-oscillatory schemes, as well as the addition of artificial diffusivities to the governing equations. Centered finite difference schemes are often avoided near discontinuities due to the tendency for significant oscillations. However, such schemes have desirable conservation properties compared to many shock-capturing schemes. The objective of this work is to derive all necessary viscous/diffusion terms from first principles and then demonstrate the performance of these analytical terms within a centered differencing framework. The physical Euler equations are spatially-filtered with a Gaussian-like filter. Sub-filter scale (SFS) terms arise in the momentum and energy equations. Analytical closure is provided for each of them by leveraging the jump conditions for a shock. No SFS terms are present in the continuity or species equations. Here, this approach is tested for several problems involving shocks in one and two dimensions. Implemented within a centered difference code, the SFS terms perform well for a range of flow conditions without introducing excessive diffusion.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2013

MHD Simulations of the Plasma Flow in the Magnetic Nozzle

The magnetohydrodynamic (MHD) flow of plasma through a magnetic nozzle is simulated by solving the governing equations for the plasma flow in the presence of an static magnetic field representing the applied nozzle. This work will numerically investigate the flow and behavior of the plasma as the inlet plasma conditions and magnetic nozzle field strength are varied. The MHD simulations are useful for addressing issues such as plasma detachment and to can be used to gain insight into the physical processes present in plasma flows found in thrusters that use magnetic nozzles. In the model, the MHD equations for a plasma, with separate temperatures calculated for the electrons and ions, are integrated over a finite cell volume with flux through each face computed for each of the conserved variables (mass, momentum, magnetic flux, energy) [1]. Stokes theorem is used to convert the area integrals over the faces of each cell into line integrals around the boundaries of each face. The state of the plasma is described using models of the ionization level, ratio of specific heats, thermal conductivity, and plasma resistivity. Anisotropies in current conduction due to Hall effect are included, and the system is closed using a real-gas equation of state to describe the relationship between the plasma density, temperature, and pressure.A separate magnetostatic solver is used to calculate the applied magnetic field, which is assumed constant for these calculations. The total magnetic field is obtained through superposition of the solution for the applied magnetic field and the self-consistently computed induced magnetic fields that arise as the flowing plasma reacts to the presence of the applied field. A solution for the applied magnetic field is represented in Fig. 1 (from Ref. [2]), exhibiting the classic converging-diverging field pattern. Previous research was able to demonstrate effects such as back-emf at a super-Alfvenic flow, which significantly alters the shape of the magnetic field in both the near- and far-field regions. However, in that work the downstream domain was constrained to a channel of constant cross-sectional area. In the present work we seek to address this issue by modeling the downstream region with a domain that permits free expansion of the plasma, permitting a better evaluation of the downstream effects the applied field has on the plasma. The inlet boundary conditions and applied magnetic field values will also be varied to determine the effect the initial plasma energy content and applied magnetic field energy density have on the near- and far-field plasma properties on the MHD code. This will determine the effect of inlet boundary conditions on the results downstream and address issues related to the restrictive numerical domain previously used.

Smith, T. E. R.↗

NASA NTRS1994

Progress in Grid Generation: From Chimera to DRAGON Grids

Hybrid grids, composed of structured and unstructured grids, combines the best features of both. The chimera method is a major stepstone toward a hybrid grid from which the present approach is evolved. The chimera grid composes a set of overlapped structured grids which are independently generated and body-fitted, yielding a high quality grid readily accessible for efficient solution schemes. The chimera method has been shown to be efficient to generate a grid about complex geometries and has been demonstrated to deliver accurate aerodynamic prediction of complex flows. While its geometrical flexibility is attractive, interpolation of data in the overlapped regions - which in today's practice in 3D is done in a nonconservative fashion, is not. In the present paper we propose a hybrid grid scheme that maximizes the advantages of the chimera scheme and adapts the strengths of the unstructured grid while at the same time keeps its weaknesses minimal. Like the chimera method, we first divide up the physical domain by a set of structured body-fitted grids which are separately generated and overlaid throughout a complex configuration. To eliminate any pure data manipulation which does not necessarily follow governing equations, we use non-structured grids only to directly replace the region of the arbitrarily overlapped grids. This new adaptation to the chimera thinking is coined the DRAGON grid. The nonstructured grid region sandwiched between the structured grids is limited in size, resulting in only a small increase in memory and computational effort. The DRAGON method has three important advantages: (1) preserving strengths of the chimera grid; (2) eliminating difficulties sometimes encountered in the chimera scheme, such as the orphan points and bad quality of interpolation stencils; and (3) making grid communication in a fully conservative and consistent manner insofar as the governing equations are concerned. To demonstrate its use, the governing equations are discretized using the newly proposed flux scheme, AUSM+, which will be briefly described herein. Numerical tests on representative 2D inviscid flows are given for demonstration. Finally, extension to 3D is underway, only paced by the availability of the 3D unstructured grid generator.

Liou, Meng-Sing↗

DOE OSTI2024

"Hall" transport of liquid crystal solitons in Couette flow

Topology establishes a unifying framework for a diverse range of scientific areas including particle physics, cosmology, and condensed matter physics. One of the most fascinating manifestations of topology in the context of condensed matter is the topological Hall effect, and its relative: the Skyrmion Hall effect. Skyrmions are stable vortex-like spin configurations in certain chiral magnets, and when subject to external electric currents can drift in the transverse direction to the current. These quasi-particles are characterised by a conserved topological charge which in the Skyrmion Hall effect plays the role of electric charges in the ordinary Hall effect. Recently, it has been shown that liquid crystals endowed with chiral properties serve as an ideal testbed for the fundamental investigation of topological solitons, including their two- and three-dimensional realisations. Here, we show experimentally and numerically that three-dimensional solitons aka "torons" exhibit a Hall-like effect when driven by shear flows: the torons are deflected in the direction perpendicular to the shear plane. The experimental results are rationalised in terms of the dynamic Ericksen-Leslie equations, which predict the emergence of the transverse component of the net mass flow, the magnitude of which scales as the 3rd power of the shear rate. The perturbation analysis highlights an interplay of the viscous and chiral elastic torques as the mechanism for the emergence of net transverse currents. Numerical simulations demonstrate, however, that torons are not merely dragged by the flow but move with their own transverse speed, much larger than the average flow velocity in the transverse direction. Our findings may enable responsive microfluidic applications relying on soft topological solitons.

FOS: Physical sciences↗

DOE OSTI2020

Development of a CPU/GPU portable software library for Lagrangian–Eulerian simulations of liquid sprays

The Lagrangian–Eulerian method is widely used for simulations of fuel sprays in turbulent combustion because of the advantage of treating the spray droplets as discrete points. One challenge of the Lagrangian–Eulerian method is the intense computational requirement when tracking the large number of Lagrangian particles needed for high fidelity. We have developed a performance-portable library, Grit, to track the Lagrangian particles in parallel on central processing unit (CPU) and graphics processing unit (GPU) accelerated high performance computing (HPC) architectures. Grit is a C++ library which employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs. The parallel algorithms, key parallel kernels, and their performances on the pre-exascale supercomputer, Summit, are presented. Grit is coupled with a direct numerical simulation (DNS) solver, S3D, for multiphase simulations. A conservative formulation has been developed and implemented in Grit for phase coupling with thermodynamic consistency. The formulation separates the conservation of mass, momentum and energy from the physical models to prevent accidental violation of conservation laws due to inconsistent models. The formulation also enforces consistent definitions of enthalpies of the fuel for both phases and the latent heat of evaporation. Finally, simulations of turbulent particle-laden flow and the evaporation of dilute turbulent spray jet are performed to verify the software implementation, and to demonstrate the scalability of Grit for large-scale multiphase simulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Spatiotemporal Variability in Transport and Reactive Processes Across a First- to Fifth-Order Fluvial Network

Fluvial networks integrate and transport constituents from terrestrial and aquatic ecosystems. To date, most research on water quality dynamics has focused on process understanding at individual streams and, as a result, there is a lack of studies analyzing how small-scale, physical and biogeochemical drivers scale across fluvial networks. We performed tracer tests in five stream orders of the Jemez River continuum in New Mexico, USA, to quantify reach-scale hyporheic exchange during two different seasonal periods to address: how do hyporheic zone contributions to riverine processing change with space and time? And does the spatiotemporal variability of hyporheic exchange scale across fluvial networks? Combining conservative (i.e., bromide) and reactive (i.e., resazurin) tracer analyses with solute transport modeling, we found a dominance of reaction-limited transport conditions in space (for all stream orders) and time (two contrasting flow regimes). Moreover, we found a trend toward more reaction-limited conditions with increasing discharge across all stream orders. While our studied fluvial network did not consistently follow the expectations of decreasing hyporheic exchange contributions with increasing stream order, as it has been previously suggested from modeling studies, we found that processing rate coefficients and Damköhler numbers consistently decreased along the continuum with increasing discharge. Our findings suggest that knowledge transferability of hyporheic zone processing along fluvial networks may be possible if the scaling patterns that we found hold true in other catchments. Thus, more research on scaling patterns along fluvial networks is highly needed.

Gootman, Kaylyn S.↗

DOE OSTI2021

The effect of rigid electron rotation on the Grad–Shafranov equilibria of a class of FRC devices

Rigid electron rotation of a Rotamak-FRC produces a pressure flux function that is more peaked than the Solov'ev flux function. This paper explores the implications of this peaked pressure flux function, including the isothermal case, which appear when the temperature profile is broader than the density profile, creating both benefits and challenges to a Rotamak-FRC based fusion reactor. In this regime, the density distribution becomes very peaked, enhancing the fusion power. The separatrix has a tendency to become oblate, which can be mitigated by flux conserving current loops. Plasma extends outside the separatrix, notably in the open field line region. Furthermore, this model does not apply to very kinetic FRCs or FRCs in which there are significant ion flows, but it may have some applicability to their outer layers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Learning unknown physics of non-Newtonian fluids

Here, we use physics-informed neural networks (PINNs) to learn viscosity models of two non-Newtonian systems (polymer melts and suspensions of particles) using only velocity measurements. For synthetic velocity data generated with the power-law viscosity model, the PINN-inferred viscosity model agrees with the analytical model for shear rates with large absolute values but deviates for shear rates near zero where the analytical model has an unphysical singularity. Once the viscosity model is learned the PINN method can solve the momentum conservation equation using only the boundary conditions.

42 ENGINEERING↗

DOE OSTI2024

Semi–Analytical Modeling of Transient Stream Drawdown and Depletion in Response to Aquifer Pumping

Analytical and semi–analytical models for stream depletion with transient stream stage drawdown induced by groundwater pumping are developed to address a deficiency in existing models, namely, the use of a fixed stream stage condition at the stream–aquifer interface. Here field data are presented to demonstrate that stream stage drawdown does indeed occur in response to groundwater pumping near aquifer–connected streams. A model that predicts stream depletion with transient stream drawdown is developed based on stream channel mass conservation and finite stream channel storage. The resulting models are shown to reduce to existing fixed–stage models in the limit as stream channel storage becomes infinitely large, and to the confined aquifer flow with a no–flow boundary at the streambed in the limit as stream storage becomes vanishingly small. The model is applied to field measurements of aquifer and stream drawdown, giving estimates of aquifer hydraulic parameters, streambed conductance, and a measure of stream channel storage. The results of the modeling and data analysis presented herein have implications for sustainable groundwater management.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

M-Star ® Software Test and Verification

Savannah River Mission Completion (SRMC) currently manages the risk for retained hydrogen in the Defense Waste Processing Facility (DWPF) vessels by implementing a Retained Hydrogen Program. The current program relies on Sludge Batch (SB) 8 Gas Chromatograph data and on conservative assumptions concerning gas release and retention to determine allowable vessel Quiescent time (Q-time). The authors identified M-Star ® CFD as a software that could simulate processes such as fluid flow, heat transfer, species transport, chemical reactions, particle transport, and retained hydrogen gas release. Preliminary simulation results suggest that more realistic assumptions on gas retention and release may be feasible for the DWPF retained hydrogen program. Because of the desire to use the M-Star ® software to perform analyses that support nuclear safety, SRMC has requested Savannah River National Laboratory (SRNL) to upgrade the software classification level of M-Star ® CFD from level D to level A to perform analyses that support nuclear safety.

08 HYDROGEN↗

NASA NTRS1978

Numerical study of transonic flow over oscillating airfoils using the full potential equation

The behavior of unsteady aerodynamic loadings on airfoils oscillating in transonic flow has been investigated numerically with particular attention given to supercritical airfoil sections. A previously developed finite difference method, which is based on the full potential equation and which uses a quasi-conservative scheme for proper capture of a shock wave motion, was employed for the present study. The unsteady aerodynamic pressure and load distributions on several different airfoil sections are presented with particular emphasis on the effects of free-stream Mach number, reduced frequency, and mean angle of attack. These parameters are demonstrated to have a significant effect on the behavior of the unsteady aerodynamic loadings. Comparisons of the present calculations with the exact inviscid solution and with the experimental results are also presented.

Isogai, K.↗

NASA NTRS1980

Comparison of quick-look plume heating calculations and Monte Carlo direct simulation

A quick-look method has been developed for estimating convective heat flux due to plume impingement in the transition flow regime on a concave surface surrounding a rocket nozzle. Comparison with the flowfield and heat fluxes computed by the Monte Carlo Direct Simulation Method show the quick-look method to be conservative. Assumptions regarding the nature of the flowfield based on engineering judgment were shown to be essentially valid. Further, the assumption of free molecular impingement was shown to be non-conservative for portions of the surface. Potential refinements of the quick-look method are discussed.

Guernsey, C. S.↗