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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Relaxation solution of the full Euler equations

A numerical procedure for the relaxation solution of the full steady Euler equations is described. By embedding the Euler system in a second order surrogate system, central differencing may be used in subsonic regions while retaining matrix forms well suited to iterative solution procedures and convergence acceleration techniques. Hence, this method allows the development of stable, fully conservative differencing schemes for the solution of quite general inviscid flow problems. Results are presented for both subcritical and shocked supercritical internal flows. Comparisons are made with a standard time dependent solution algorithm. Previously announced in STAR as N82-24859

Johnson, G. M.↗
NASA NTRS1985

Breakdown of the conservative potential equation

The conservative full-potential equation is used to study transonic flow over five airfoil sections. The results of the study indicate that once shock waves are present in the flow, the qualitative approximation is different from that observed with the Euler equations. The difference in behavior of the potential eventually leads to multiple solutions.

Salas, M. D.↗
NASA NTRS1985

A fast efficient implicit scheme for the gasdynamic equations using a matrix reduction technique

An efficient implicit finite-difference algorithm for the gasdynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved without loss of the stability characteristics generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows and used to show that the method offers computational advantages over the conventional Beam and Warming scheme. Existing Beam and Warming codes can be retrofit with minimal effort. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations.

Barth, T. J.↗
NASA NTRS1984

Dynamical relations for left ventricular ejection - Flow rate, momentum, force and impulse

An investigation was carried out to quantitatively evaluate left ventricular volume flow rate, momentum, force and impulse derived from application of conservation principles for mass and momentum of blood within the ventricle during the ejection phase. An automated digital image processing system was developed and applied to left ventricular angiograms which are computer processed and analyzed frame by frame to determine the dynamical relations by numerical methods. The initial experience with force and impulse has indicated that neither quantity seemed to be a sensitive indicator of coronary artery disease as evaluated by qualitative angiography for the particular patient group studied. Utilization of the dynamical relations in evaluating human left ventricular performance requires improved means of measurement and interpretation of clinical studies.

Back, L. H.↗
NASA NTRS1986

Breakdown of the conservative potential equation

The conservative full-potential equation is used to study transonic flow over five airfoil sections. The results of the study indicate that once shock are present in the flow, the qualitative approximation is different from that observed with the Euler equations. The difference in behavior of the potential eventually leads to multiple solutions.

Salas, M. D.↗
NASA NTRS1988

Preshock region acceleration of implanted cometary H(+) and O(+)

A self-consistent, three-fluid model of plasma transport and implanted ion acceleration in the unshocked solar wind is presented. The solar wind plasma is depleted by charge exchange with the expanding cometary exosphere, while implanted protons and heavy ions are produced by photoionization and charge transfer and lost by charge exchange. A generalized transport equation describing convection, adiabatic and diffusive velocity change, and the appropriate production terms is used to describe the evolution of the two cometary ion components, while the moments of the Boltzmann equation are used to calculate the solar wind density and pressure. The flow velocity is obtained self-consistently by combining the conservation equations of the three ion species. The results imply that second-order Fermi acceleration can explain the implanted spectra observed in the unshocked solar wind. Comparison of measured and calculated distribution indicates that spatial diffusion of implanted ions probably plays an important role in forming the energetic particle environment in the shock vicinity.

Gombosi, Tamas I.↗
NASA NTRS1988

Second order accurate finite difference approximations for the transonic small disturbance equation and the full potential equation

New shock-capturing finite difference approximations for solving two scalar conservation law nonlinear partial differential equations describing inviscid, isentropic, compressible flows of aerodynamics at transonic speeds are presented. A global linear stability theorem is applied to these schemes in order to derive a necessary and sufficient condition for the finite element method. A technique is proposed to render the described approximations total variation-stable by applying the flux limiters to the nonlinear terms of the difference equation dimension by dimension. An entropy theorem applying to the approximations is proved, and an implicit, forward Euler-type time discretization of the approximation is presented. Results of some numerical experiments using the approximations are reported.

Mostrel, M. M.↗
NASA NTRS1996

High Resolution Euler Solvers Based on the Space-Time Conservation Element and Solution Element Method

The I-D, quasi I-D and 2-D Euler solvers based on the method of space-time conservation element and solution element are used to simulate various flow phenomena including shock waves, Mach stem, contact surface, expansion waves, and their intersections and reflections. Seven test problems are solved to demonstrate the capability of this method for handling unsteady compressible flows in various configurations. Numerical results so obtained are compared with exact solutions and/or numerical solutions obtained by schemes based on other established computational techniques. Comparisons show that the present Euler solvers can generate highly accurate numerical solutions to complex flow problems in a straightforward manner without using any ad hoc techniques in the scheme.

Wang, Xiao-Yen↗
NASA NTRS2012

Ice Accretion Modeling using an Eulerian Approach for Droplet Impingement

A three-dimensional Eulerian analysis has been developed for modeling droplet impingement on lifting bodes. The Eulerian model solves the conservation equations of mass and momentum to obtain the droplet flow field properties on the same mesh used in CFD simulations. For complex configurations such as a full rotorcraft, the Eulerian approach is more efficient because the Lagrangian approach would require a significant amount of seeding for accurate estimates of collection efficiency. Simulations are done for various benchmark cases such as NACA0012 airfoil, MS317 airfoil and oscillating SC2110 airfoil to illustrate its use. The present results are compared with results from the Lagrangian approach used in an industry standard analysis called LEWICE.

Kim, Joe Woong↗
DOE OSTI2023

Operation Optimization Using Reinforcement Learning with Integrated Artificial Reasoning Framework

In large and complex systems, operational decision-making requires a systematic analysis with a vast amount of data from both process parameters and component status monitoring. In this paper, we present an integrated artificial reasoning approach for system state transition models that can help operational decision-making with explainable and traceable reasoning. The integrated artificial reasoning framework is a physics-based approach of defining the system structure in a Bayesian network, so we leveraged it in a Markov decision process (MDP) for finding optimal operational solutions. In our proposed framework, the MDP is implemented on a dynamic Bayesian network (DBN), which represents causalities in a system. The multilevel flow modeling was utilized in order to extract these causalities in a more efficient and objective manner. Since multilevel flow modeling is based on the fundamental energy and mass conservation laws, the target system is decomposed into several mass, energy, and information structures, which serve as the basis for a DBN. The MDP consists of the processes of finding a solution for the Bellman equation, which can be derived from the conditional probability equations of the constructed DBN. System operators can capture stochastic system dynamics as multiple subsystem state transitions based on their physical relations and uncertainties coming from the component degradation process or random failures. We analyzed a simplified example system to illustrate finding an optimal operational policy with this approach.

Kim, Junyung↗
DOE OSTI2022

Implicit shock tracking for unsteady flows by the method of lines

Here, a recently developed high-order implicit shock tracking (HOIST) framework for resolving discontinuous solutions of inviscid, steady conservation laws is extended to the unsteady case. Central to the framework is an optimization problem which simultaneously computes a discontinuity-aligned mesh and the corresponding high-order approximation to the flow, which provides nonlinear stabilization and a high-order approximation to the solution. This work extends the implicit shock tracking framework to the case of unsteady conservation laws using a method of lines discretization via a diagonally implicit Runge-Kutta method by “solving a steady problem at each timestep”. We formulate and solve an optimization problem that produces a feature-aligned mesh and solution at each Runge-Kutta stage of each timestep, and advance this solution in time by standard Runge-Kutta update formulas. A Rankine-Hugoniot based prediction of the shock location together with a high-order, untangling mesh smoothing procedure provides a high-quality initial guess for the optimization problem at each time, which results in rapid convergence of the sequential quadratic programing (SQP) optimization solver. This method is shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover the design accuracy of the Runge-Kutta scheme. We demonstrate this framework on a series of inviscid, unsteady conservation laws in both one- and two- dimensions. We also verify that our method is able to recover the design order of accuracy of our time integrator in the presence of a strong discontinuity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1996

Multigrid, Fractional-Step Computation Of Flow

Speed of computer code solving three-dimensional Navier-Stokes equations of flow of incompressible fluid by fractional-step method increased significantly by use of multigrid procedures. In method, equations solved on general nonorthogonal curvilinear coordinate grid, using volume fluxes. At each time step, computations performed in two fractional steps. In first step, equations of conservation of momentum solved by use of gradient of pressure from previous time step via explicit approximate-factorization method, yielding approximate flow field that does not satisfy equation of conservation of mass. In second step, discrete Poisson-like equation with Neumann-type boundary conditions, formed by combining equations of conservation of momentum and mass, solved iteratively.

Kwak, Dochan↗
DOE OSTI2023

Faster-than-real-time Simulation with Demonstration for Resilient DER Integration

The US electric grid is facing operational, stability, and security challenges. Transmission system operators need some measure of visibility into distribution system renewable generation. Distribution system generation needs to support transmission system voltage. The grid is experiencing an expansion in measurement systems. How to take full advantage of this expansion and defend against attacks, both cyber and physical, poses additional challenges. The Faster-than-real-time Simulation with demonstration for Resilient DER Integration project set out to do the following: a. Flatten the voltage profile through the feeders and system for cost saving and voltage stabilization needs. b. Increase the amount of intermittent distributed energy resources (IDERs) that could be deployed on a utility feeder and provide 100% or more energy needed for the demands on that feeder, and c. based on an accurate model (Digital Twin) of the utilities system, be able to detect any abnormalities on the utilities distribution system. To manage the voltage and increase IDER penetration (a,b), Graph Trace Analysis is employed in a time-series, optimal power flow to coordinate the time-varying feedback control setpoints of a distribution feeder’s utility control devices. Under the coordinated control are a Load Tap Changing Transformer, a voltage regulator, and five switched capacitor banks. The feeder serves over 2000 customers, the feeder secondaries are modeled, and the feeder has 2.3 MW of PV generation, corresponding to a 17.4% penetration of PV generation. The feeder model has over 12,000 components, where every customer load bus and PV generator are modeled. The accuracy of the power flow solution is compared against historical meter voltage measurements, the improvement in conservation voltage reduction energy savings as a function of the coordinated control desired voltage profile is investigated, and the increase in PV penetration of the coordinated control over the existing control is presented. To achieve improved control performance while observing system operation constraints, bellwether Advanced Metering Infrastructure (AMI) voltage measurements are used to adjust the desired voltage profile used by the optimal power flow analysis. To detect and alleviate or negate attacks or failures on the distribution and transmission utility grids (c) the grid needs to be resilient and self-healing. In this project software was designed to do just that. At the center of the software is an Integrated System Model (ISM) that spans from transmission to secondary distribution. The ISM is employed in real-time abnormality detection, voltage stability forecasting, and multi-mode control. Testing results are presented for: 1—attacks on utility infrastructure; 2—energy savings from optimal control; 3—distribution system control response during a low voltage transmission system event; 4—cyber-attacks on PV inverters, where physical inverters are used in hard-ware-in-the-simulation-loop studies. Contributions of this work include real-time analysis that spans from three-phase transmission through secondary distribution; an approach for detecting abnormalities that employs measurements from three independent measurement systems; and a multi-mode distribution system control that responds to cyber-attacks, physical attacks, equipment failures, and transmission system needs.

Integrated System Model, Graph Trace Analysis, Adv↗
NASA NTRS2017

Numerical Analysis of Mixed-Phase Icing Cloud Simulations in the NASA Propulsion Systems Laboratory

This presentation describes the development of a numerical model that couples the thermal interaction between ice particles, water droplets, and the flowing gas of an icing wind tunnel for simulation of NASA Glenn Research Centers Propulsion Systems Laboratory (PSL). The ultimate goal of the model is to better understand the complex interactions between the test parameters and have greater confidence in the conditions at the test section of the PSL tunnel. The model attempts to explain the observed changes in test conditions by coupling the conservation of mass and energy equations for both the cloud particles and flowing gas mass. Model predictions were compared to measurements taken during May 2015 testing at PSL, where test conditions varied gas temperature, pressure, velocity and humidity levels, as well as the cloud total water content, particle initial temperature, and particle size distribution.

NASA PSL↗
NASA NTRS2017

Numerical Analysis of Mixed-Phase Icing Cloud Simulations in the NASA Propulsion Systems Laboratory

This paper describes the development of a numerical model that couples the thermal interaction between ice particles, water droplets, and the flowing gas of an icing wind tunnel for simulation of NASA Glenn Research Centers Propulsion Systems Laboratory (PSL). The ultimate goal of the model is to better understand the complex interactions between the test parameters and have greater confidence in the conditions at the test section of the PSL tunnel. The model attempts to explain the observed changes in test conditions by coupling the conservation of mass and energy equations for both the cloud particles and flowing gas mass. Model predictions were compared to measurements taken during May 2015 testing at PSL, where test conditions varied gas temperature, pressure, velocity and humidity levels, as well as the cloud total water content, particle initial temperature, and particle size distribution.

NASA PSL↗
NASA NTRS1975

ECASTAR: Energy Conservation; an Assessment of Systems, Technologies and Requirements

A methodology for a systems approach display and assessment of the potential for energy conservation actions and the impacts of those actions was presented. The U.S. economy is divided into four sectors: energy industry, industry, residential/commercial and transportation. Each sector is assessed with respect to energy conservation actions and impacts. The four sectors are combined and three strategies for energy conservation actions for the combined sectors are assessed. The three strategies (national energy conservation, electrification and diversification) represent energy conservation actions for the near term (now to 1985), the mid term (1985 to 2000) and the far term (2000 and beyond). The assessment procedure includes input/output analysis to bridge the flows between the sectors, and net economics and net energetics as performance criteria for the conservation actions. Targets of opportunity for large net energy net energy savings and the application of technology to achieve these savings are discussed.

Source record↗
NASA NTRS2013

Effects of Nonequilibrium Chemistry and Darcy-Forchheimer Pyrolysis Flow for Charring Ablator

The fully implicit ablation and thermal response code simulates pyrolysis and ablation of thermal protection materials and systems. The governing equations, which include energy conservation, a three-component decomposition model, and a surface energy balance, are solved with a moving grid.This work describes new modeling capabilities that are added to a special version of code. These capabilities include a time-dependent pyrolysis gas flow momentum equation with Darcy-Forchheimer terms and pyrolysis gas species conservation equations with finite rate homogeneous chemical reactions. The total energy conservation equation is also enhanced for consistency with these new additions. Two groups of parametric studies of the phenolic impregnated carbon ablator are performed. In the first group, an Orion flight environment for a proposed lunar-return trajectory is considered. In the second group, various test conditions for arcjet models are examined. The central focus of these parametric studies is to understand the effect of pyrolysis gas momentum transfer on material in-depth thermal responses with finite-rate, equilibrium, or frozen homogeneous gas chemistry. Results indicate that the presence of chemical nonequilibrium pyrolysis gas flow does not significantly alter the in-depth thermal response performance predicted using the chemical equilibrium gas model.

Chen, Yih-Kanq↗
NASA NTRS2006

Cassini Observations of Saturn's Dawn-Magnetotail Region and their Relation to Models of Saturn's Aurora: Preliminary Results

Using Cassini plasma and magnetic field observations from the dawn meridian of Saturn s outer magnetosphere to Saturn s magnetotail region, we investigate the applicability of the centrifugal instability model by Sittler et al. [2006] for Saturn s auroral response to the solar wind, versus the reconnection model of Saturn s aurora by Cowley et al. [2005]. We use Cassini Plasma Spectrometer (CAPS) Ion Mass Spectrometer (IMS) and Electron Plasma Spectrometer (ELS) observations to characterize the plasma environment. ELS and magnetometer observations are used to map out the morphology of the outer magnetosphere from dawn to midnight local time. IMS observations are used to measure plasma flow velocities from which one can infer rotation versus convective flows. IMS composition measurements are used to trace the source of plasma from the inner magnetosphere (protons, H2 and water group ions) versus an external solar wind source (protons and Heff ions). A critical parameter for both models is the strength of the convection electric field with respect to the rotational electric field for the large scale magnetosphere. Is there a significant return flow from the magnetotail? Pitch angle distributions also play an important role as a discriminator. If the magnetosphere tends to conserve angular momentum as suggested by Sittler et al. [2006], then we expect to see an anti-correlation between rotational flow component and radial flow velocities. All will be investigated.

Sittler, E. C.↗