Engineering Papers⌕ Search

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.

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

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Non-equilibrium low-temperature plasma-assisted combustion of iso-octane: Perturbing pyrolysis and oxidation kinetics

Here, in this study, a plasma-coupled flow reactor facility is used to examine the effects of non-equilibrium low-temperature plasmas on perturbing the pyrolysis and oxidation kinetics of iso-octane. Experiments were performed in highly dilute reactive mixtures of nitrogen, at near isothermal conditions for temperatures ranging from 523 K to 1203 K. Experiments cumulatively demonstrated enhanced chemical reactivity with the plasma for temperatures below 900 K, and a lowering of the hot-ignition temperature. Detailed kinetic insight was derived from a 0D plasma-coupled kinetic model, utilizing a constructed mechanism that combined both plasma-specific chemistry and the neutral combustion chemistry. For pyrolysis conditions, the model displayed relatively good agreement with fuel consumption and the formation of most intermediates compared to the experimental data, demonstrating the model is able to accurately predict primary radical formation from the plasma directly interacting with the fuel. Enhanced reactivity was attributed to collisional quenching of excited-states of N 2 with fuel, which led to efficient fuel fragmentation and enhancement of the H-radical flux. For oxidation conditions, the model displayed satisfactory agreement with the experiments. Model predictions were able to accurately predict fuel consumption and most intermediate speciation data for T > 800 K, but most discrepancies were towards T < 800 K in particular with oxygenated intermediates. In the presence of oxygen, plasma effects were predominantly spent on efficient enhancement of O- and H-radical fluxes, leading to further fuel fragmentation and initiation of both the OH- and HO 2 -radical pools. Subsequent reactivity of iso-octane was then dictated by the response of the temperature-dependent neutral chemistry. At low-temperatures (T = 643 K), enhanced fuel radicals and O 2 -additon chemistry lead to the formation of oxygenated species, while at intermediate temperatures (T = 843 K) net decrease in OH-radical reactivity led to an increase in hydrocarbon speciation. Near the self-ignition threshold (T = 1163 K), radicals generated by high-temperature branching reactions dominate the oxidation process and effectively ignition. This study ultimately demonstrated that the enhancement of radicals afforded by the plasma causes a deviation in known understanding of iso-octane kinetics in some regards and warrants future studies to reconcile these discrepancies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1994

Axisymmetric modeling of cometary mass loading on an adaptively refined grid: MHD results

The first results of an axisymmetric magnetohydrodynamic (MHD) model of the interaction of an expanding cometary atmosphere with the solar wind are presented. The model assumes that far upstream the plasma flow lines are parallel to the magnetic field vector. The effects of mass loading and ion-neutral friction are taken into account by the governing equations, whcih are solved on an adaptively refined unstructured grid using a Monotone Upstream Centered Schemes for Conservative Laws (MUSCL)-type numerical technique. The combination of the adaptive refinement with the MUSCL-scheme allows the entire cometary atmosphere to be modeled, while still resolving both the shock and the near nucleus of the comet. The main findingsare the following: (1) A shock is formed approximately = 0.45 Mkm upstream of the comet (its location is controlled by the sonic and Alfvenic Mach numbers of the ambient solar wind flow and by the cometary mass addition rate). (2) A contact surface is formed approximately = 5,600 km upstream of the nucleus separating an outward expanding cometary ionosphere from the nearly stagnating solar wind flow. The location of the contact surface is controlled by the upstream flow conditions, the mass loading rate and the ion-neutral drag. The contact surface is also the boundary of the diamagnetic cavity. (3) A closed inner shock terminates the supersonic expansion of the cometary ionosphere. This inner shock is closer to the nucleus on dayside than on the nightside.

Gombosi, Tamas I.↗
NASA NTRS1994

Magnetocentrifugally driven flows from young stars and disks. 2: Formulation of the dynamical problem

We formulate the dynamical problem of a cool wind centrifugally driven from the magnetic interface of a young star and an adjoining Keplerian disk. We examine the situation for mildly accreting T Tauri stars that rotate slowly as well as rapidly accreting protostars that rotate near break-up. In both cases a wind can be driven from a small X-region just outside the stellar magnetopause, where the field lines assume an open geometry and are rooted to material that rotates at an angular speed equal both to the local Keplerian value and to the stellar angular speed. Assuming axial symmetry for the ideal magnetohydrodynamic flow, which requires us to postpone asking how the (lightly ionized) gas is loaded onto field lines, we can formally integrate all the governing equations analytically except for a partial equation that describes how streamlines spread in the meridional plane. Apart from the difficulty of dealing with PDEs of mixed type, finding the functional forms of the conserved quantities along streamlines - the ratio beta of magnetic field to mass flux, the specific energy H of the fluid in the rotating frame, and the total specific angular momentum J carried in the matter and the field - constitutes a standard difficulty in this kind of (Grad-Shafranov) formalism. Fortunately, because the ratio of the thermal speed of the mass-loss regions to the Keplerian speed of rotation of the interface constitutes a small parameter epsilon, we can attack the overall problem by the method of matched asymptotic expansions. This procedure leads to a natural and systematic technique for obtaining the relevant functional dependences of beta, H, and J. Moreover, we are able to solve analytically for the properties of the flow emergent from the small transsonic region driven by gas pressure without having to specify the detailed form of any of the conserved functions, beta, H, and J. This analytical solution provides inner boundary conditions for the numerical computation in a companion paper by Najita & Shu of the larger region where the main acceleration to terminal speeds occurs.

Shu, Frank H.↗
DOE OSTI2021

2019 NETL Related Energy Diagrams

This document illustrates the energy production, consumption, use, flow, and CO 2 emission in the U.S. in 2019.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

CFD Simulation of Biomass Pyrolysis Vapor Upgrading over a Pt/TiO 2 Catalyst in Fixed and Moving Beds

This report presents a comprehensive computational fluid dynamics (CFD) study of the biomass pyrolysis vapor phase upgrading and the catalyst regeneration in different reactors under various operational conditions. This study used the open-source software Multiphase Flow with Interphase eXchanges (MFiX), developed at National Energy Technology Laboratory (NETL) for simulating hydrodynamics, heat transfer, and chemical reactions in multiphase systems. This simulation solved the conservation of mass, momentum, energy, species for gas and solid phases with chemical reactions.

09 BIOMASS FUELS↗
DOE OSTI2023

Bipolar Membranes with an Electrospun 3D Junction

Freestanding bipolar membranes (BPMs) with an entirely new and transformative morphology were fabricated and characterized. The key disruptive design element was a water-splitting/water-generation junction layer of finite thickness composed of interwoven and interlocking electrospun anion-exchange polymer (AEP) fibers and cation-exchange polymer (CEP) fibers. Methods were developed to sandwich the 3D junction layer between two pre-formed dense AEP and CEP films., with the closure of all interfiber voids, where the processing steps are easily incorporated into a commercial roll-to-roll membrane manufacturing scheme. A series of membranes were made using different polymers and junction layers to identify the optimum morphology and composition for water-splitting and water-generation applications. For water splitting, the 3D junction BPMs worked remarkably well with operating current densities at/above 1.0 A/cm2 (10X greater than commercial BPMs) at a transmembrane voltage drop of only 1.1 V. In water generation mode, the 3D junction BPMs operated at 0.5 A/cm2, a world-record current density. The membranes were found to exhibit outstanding durability and can be manufactured at scale for low cost. The high operating current densities and stable long-term operation are due to the 3D junction layer design, where there is a high interfacial area for water splitting or water generation reactions and where the interlocking fibers prevent delamination of the outer films. These new BPMs are ideal candidate materials for both existing water-splitting electrodialysis separations and new electrochemical processes, such as large-scale direct air capture of CO 2 , reactors for CO 2 reduction, self-hydrating fuel cells, and redox flow batteries.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS1973

Introduction to gasdynamics of explosions

Questions concerning the genesis and sustenance of an explosion are investigated, giving attention to the mechanics of explosions, the gasdynamics of explosions, aspects of technological significance, and future prospects. The dynamics of exothermic processes is discussed together with the most prominent effects of explosions. Blast waves are considered, taking into account conservation principles, blast wave transformation, conservative equations in nondimensional form, the equation of state, Eulerian space profiles, Eulerian time profiles, Lagrangian time profiles, boundary conditions and integral relations, and self-similar flow fields.

Oppenheim, A. K.↗
NASA NTRS1980

Numerical calculation of transonic axial turbomachinery flows

This paper presents a numerical method and the results of a computer program for solving an exact, three-dimensional, full-potential equation that models rotating and nonrotating inviscid, absolutely irrotational, homentropic flows. Besides calculating the flows through an arbitrarily shaped rotor or stator blade row mounted on an axisymmetric hub and confined in an axisymmetric duct, the computer program is also capable of analysing flow fields about arbitrarily shaped wing-body combinations, propellers, helicopter rotors in hover, and wind turbine rotors. The governing equation is solved numerically in a fully conservative form by using an artificial time concept, a finite volume technique, rotated type-dependent differencing, successive line overrelaxation, and sequential boundary-conforming grid refinement. An artificial viscosity is added in fully conservative form; and an initial guess for the potential field is applied, as determined by a two-dimensional cascade analysis.

Dulikravich, D. S.↗
NASA NTRS1978

Recent experiences with three-dimensional transonic potential flow calculations

Some recent experiences with computer programs capable of solving finitie-difference approximations to the full potential equation for the transonic flow past three dimensional swept wings and simple wing-fuselage combinations are discussed. The programs used are a nonconservative program for swept wings, a quasi-conservative finite-volume program capable of treating swept wings mounted on fuselages of slowly varying circular cross section, and a fully conservative finite volume scheme capable of treating swept wings and wing-cylinder combinations. The present capabilities of these codes are reviewed. The relative merits of the conservative and nonconservative formulations are discussed, and the results of calculations including corrections for the boundary-layer displacement effect are presented.

Caughey, D. A.↗
NASA NTRS1983

Numerical studies of the formation and destruction of vortices in a motored four-stroke piston-cylinder configuration

A finite-difference procedure which solves the conservation equations of mass, momentum, and energy is used to investigate the effects of the compression ratio, engine speed, bore-to-stroke ratio, and air intake flow angle on the turbulent flow field within an axisymmetric piston-cylinder configuration. It is shown that in a four-stroke piston-cylinder configuration, the intake stroke is characterized by the formation of a piston vortex. The piston vortex is stretched during the intake stroke, and the head vortex has an almost constant diameter. For a 0-deg air intake flow angle, both vortices disappear by the end of the compression stroke; for an air intake flow angle of 45 deg, the flow field within the cylinder shows three elongated vortices which persist into the compression stroke and then break up and merge. It is also shown that larger bore-to-stroke ratios give rise to lower turbulent levels than smaller bore-to-stroke ratios and that the turbulent intensity is almost independent of the rpm.

Schock, H. J.↗
NASA NTRS1984

Fresh attack on laminar flow

Laminar flow control is a technology with great potential for aircraft drag reduction. Stabilization of laminar boundary layers became known as natural laminar flow (NLF) and research led to the development of NLF airfoils. Research was also conducted on stabilization by suction, referred to as laminar flow control (LFC). Experiments demonstrated that extensive laminar flow could be achieved in flight. However, there remained doubts regarding the practicality of producing, with the technology then available, wing surfaces sufficiently smooth and wavefree to meet laminar-flow criteria and maintaining the wing surface quality in normal service. In 1976, the Aircraft Energy Efficiency (ACEE) program was begun by NASA to develop fuel-conservative technology for commercial transports. The progress of the ACEE program is discussed. Attention is given to LFC wing structures, and LFC leading-edge systems.

Wagner, R. D.↗
NASA NTRS1985

Drag characteristics of unsteady, perturbed boundary layer flows

A series of time-dependent numerical computations have been performed for the flow of an unsteady boundary layer over an embedded cavity, which is aligned normal to the flow direction. The unsteady flow is the result of the imposition of a Stuart vortex, at the inflow boundary of the computational domain, onto an otherwise unperturbed two-dimensional laminar boundary layer flow. This produces a strongly vortical, but spatially monochromatic, motion at inflow which is allowed to evolve downstream through the time-dependent numerical solution of appropriate conservation equations. As the imposed vortex evolves over the solid boundary, a weaker induced vortical motion is produced and together these perturbing motions interact with the fluid motion in the embedded cavity. An analysis of the pertinent dynamic variables as well as the time-dependent drag characteristics for this type of non-planar geometry is performed; with particular emphasis on the relative contribution of both pressure drag and frictional drag at various times in the evolution process.

Gatski, T. B.↗
NASA NTRS1988

Numerical modeling of multidimensional flow in seals and bearings used in rotating machinery

The rotordynamic behavior of turbomachinery is critically dependent on fluid dynamic rotor forces developed by various types of seals and bearings. The occurrence of self-excited vibrations often depends on the rotor speed and load. Misalignment and rotor wobbling motion associated with differential clearance were often attributed to stability problems. In general, the rotative character of the flowfield is a complex three dimensional system with secondary flow patterns that significantly alter the average fluid circumferential velocity. A multidimensional, nonorthogonal, body-fitted-grid fluid flow model is presented that describes the fluid dynamic forces and the secondary flow pattern development in seals and bearings. Several numerical experiments were carried out to demonstrate the characteristics of this complex flowfield. Analyses were performed by solving a conservation form of the three dimensional Navier-Stokes equations transformed to those for a rotating observer and using the general-purpose computer code PHOENICS with the assumptions that the rotor orbit is circular and that static eccentricity is zero. These assumptions have enabled a precise steady-state analysis to be used. Fluid injection from ports near the seal or bearing center increased fluid-film direct dynamic stiffness and, in some cases, significantly increased quadrature dynamic stiffness. Injection angle and velocity could be used for active rotordynamic control; for example, injection, when compared with no injection, increased direct dynamic stiffness, which is an important factor for hydrostatic bearings.

Hendricks, R. C.↗
NASA NTRS1987

Numerical and analytical study of fluid dynamic forces in seals and bearings

Study results for rotor/bearing/seal system models are compared with both a numerical model based on a transformed, conservative version of the three-dimensional Navier-Stokes equation, and an analytical model based on 'lumped' fluid parameters. Rotor-destabilizing factors are related to flow field rotation; it is shown that these factors can be reduced through a decrease in the fluid average circumferential velocity. A transformation of the Navier-Stokes equations to those for a rotating observer and a numerical code with a nonorthogonal body-fitted grid are used to carry out numerical experiments demonstrating the character of the complex flow field.

Tam, L. T.↗
NASA NTRS1991

The effect of surface radiation on flame spread in a quiescent, microgravity environment

In most theories of laminar flame spread over solid combustibles, radiation is neglected. However, when gas motion is completely absent, which may occur in the quiescent, microgravity environment of a spacecraft, the importance of radiation compared to convection is enhanced. Here, a theoretical model is developed of flame spread over a thin solid fuel into an opposing flow of oxidizer, the opposing flow being present in the quiescent environment in flame-fixed coordinates, including the effects of radiative heat transfer from the fuel surface. Numerical solutions to the conservation equations in the gas and the solid phase that describe the spreading flame are presented for a variety of ambient conditions and surface emittances. When surface radiation is significant, the solid surface acts as a heat sink, causing the flame to cool, shrink in size, and spread more slowly. Better agreement with experiment for net heat transfer to the surface is obtained when surface radiation is included. For low ambient oxygen levels, radiative effects lead to flame extinction. A radiation/conduction parameter is identified that adequately describes the importance of surface radiation.

Bhattacharjee, S.↗
NASA NTRS2000

Computational Aeroacoustics by the Space-Time CE/SE Method

The Space-Time Conservation Element and Solution Element Method, or CE/SE Method for short, is a newly developed numerical method for conservation laws. Despite its second order accuracy, it possesses low dispersion errors and low dissipation. The method is robust enough to cover a wide spectrum of compressible flows: from weak linear acoustic waves to strong discontinuous waves (shocks). An outstanding feature of the CE/SE scheme is its novel, simple but effective non-reflecting boundary condition (NRBC), which is particularly valuable for CAA (computational aeroacoustics). In this seminar, the 1-D and 2-D unstructured version of the CE/SE schemes are first briefly described. Secondly, some discussions on the NRBC are given. Then, various examples for linear, nonlinear aeroacoustics are presented.

Loh, Ching Y.↗
NASA NTRS2017

Numerical Modeling of Saturated Boiling in a Heated Tube

This paper describes a mathematical formulation and numerical solution of boiling in a heated tube. The mathematical formulation involves a discretization of the tube into a flow network consisting of fluid nodes and branches and a thermal network consisting of solid nodes and conductors. In the fluid network, the mass, momentum and energy conservation equations are solved and in the thermal network, the energy conservation equation of solids is solved. A pressure-based, finite-volume formulation has been used to solve the equations in the fluid network. The system of equations is solved by a hybrid numerical scheme which solves the mass and momentum conservation equations by a simultaneous Newton-Raphson method and the energy conservation equation by a successive substitution method. The fluid network and thermal network are coupled through heat transfer between the solid and fluid nodes which is computed by Chen's correlation of saturated boiling heat transfer. The computer model is developed using the Generalized Fluid System Simulation Program and the numerical predictions are compared with test data.

Majumdar, Alok↗
NASA NTRS2019

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗