Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conservation equations”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

A finite difference scheme for three-dimensional steady laminar incompressible flow

A finite difference scheme for three-dimensional steady laminar incompressible flows is presented. The Navier-Stokes equations are expressed conservatively in terms of velocity and pressure increments (delta form). First order upwind differences are used for first order partial derivatives of velocity increments resulting in a diagonally dominant matrix system. Central differences are applied to all other terms for second order accuracy. The SIMPLE pressure correction algorithm is used to satisfy the continuity equation. Numerical results are presented for cubic cavity flow problems for Reynolds numbers up to 2000 and are in good agreement with other numerical results.

Hwang, Danny P.↗
DOE OSTI2024

A data-driven sensor placement approach for detecting voltage violations in distribution systems

Stochastic fluctuations in power injections from distributed energy resources (DERs) combined with load variability can cause constraint violations (e.g., exceeded voltage limits) in electric distribution systems. To monitor grid operations, sensors are placed to measure important quantities such as the voltage magnitudes. Here, in this paper, we consider a sensor placement problem which seeks to identify locations for installing sensors that can capture all possible violations of voltage magnitude limits. We formulate a bilevel optimization problem that minimizes the number of sensors and avoids false sensor alarms in the upper level while ensuring detection of any voltage violations in the lower level. This problem is challenging due to the nonlinearity of the power flow equations and the presence of binary variables. Accordingly, we employ recently developed conservative linear approximations of the power flow equations that overestimate or underestimate the voltage magnitudes. By replacing the nonlinear power flow equations with conservative linear approximations, we can ensure that the resulting sensor locations and thresholds are sufficient to identify any constraint violations. Additionally, we apply various problem reformulations to significantly improve computational tractability while simultaneously ensuring an appropriate placement of sensors. Lastly, we improve the quality of the results via an approximate gradient descent method that adjusts the sensor thresholds. We demonstrate the effectiveness of our proposed method for several test cases, including a system with multiple switching configurations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1979

Description and preliminary results of the 9-level UCLA general circulation model

The 9-level general circulation model, with the upper boundary of 50 mb, is a modification of the earlier 6- and 12-level versions used at UCLA. Two major design changes have been made in an effort to improve the model: (1) a potential enstrophy conserving scheme in the equation of motion has replaced the previous scheme which only conserved enstrophy for nondivergent flow, and (2) the treatment of the planetary boundary layer has been modified by making the predicted planetary boundary layer top a coordinate surface. Experiments with idealized steep mountains have been made to test each of the major modifications, and the model is now ready for general circulation simulation experiments.

Suarez, M. J.↗
NASA NTRS1991

Sensing and actuation in a liquid propellant rocket engine

Equations describing the conservation of mass, momentum, and energy, and multicomponent transport and chemical reaction in a liquid propellant rocket (LPR) engine are given, along with appropriate boundary conditions. The formulation includes a detailed description of the two-phase (gas/liquid) flow in terms of a moving boundary problem, with droplets and vapor occupying time-dependent subdomains of the LPR chamber. Several sensor and actuator types are considered in detail, and their mathematical descriptions are incorporated into the differential equations and boundary conditions.

Bentsman, Joseph↗
NASA NTRS1994

Particle Distribution in Suspension Shear Flow

An earlier proposed constitutive relation for normal stresses originated by random particle fluctuations is used to describe a joint effect of thermal and shear-induced fluctuations on concentrational distributions in suspension flow. Averaged products of components of the fluctuation velocity are evaluated on a basis of the rational mechanics approach combined with a simple kinematic consideration. The equation of momentum conservation of the dispersed phase of a suspension closed with this constitutive relation is applied to unidirectional shear flow in the gravity field and to rotational Couette flow. Coupling of the thermal and shear-induced fluctuations results in that the ability of shear flow to suspend particles has a minimum at a certain particle size, all other things being equal. The developed model provides also for a reasonable explanation of particle distributions observed in Couette flow. The approach based on the consideration of momentum balance for the dispersed phase is proved to lead to an effective equation of convective diffusion of the suspended particles. Coefficients of mutual diffusion due to both thermal and shear-induced fluctuations are drastically different from corresponding self-diffusivities as regards both their scaling and their concentrational dependence.

Buyevich, Yu A.↗
DOE OSTI2022

Energy conserving and well-balanced discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry

ABSTRACT This paper presents high-order Runge–Kutta (RK) discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry. The scheme can preserve a general polytropic equilibrium state and achieve total energy conservation up to machine precision with carefully designed spatial and temporal discretizations. To achieve the well-balanced property, the numerical solutions are decomposed into equilibrium and fluctuation components that are treated differently in the source term approximation. One non-trivial challenge encountered in the procedure is the complexity of the equilibrium state, which is governed by the Lane–Emden equation. For total energy conservation, we present second- and third-order RK time discretization, where different source term approximations are introduced in each stage of the RK method to ensure the conservation of total energy. A carefully designed slope limiter for spherical symmetry is also introduced to eliminate oscillations near discontinuities while maintaining the well-balanced and total-energy-conserving properties. Extensive numerical examples – including a toy model of stellar core collapse with a phenomenological equation of state that results in core bounce and shock formation – are provided to demonstrate the desired properties of the proposed methods, including the well-balanced property, high-order accuracy, shock-capturing capability, and total energy conservation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1979

Cropping management using color and color infrared aerial photographs

The Universal Soil Loss Equation (USLE) is a widely accepted tool for erosion prediction and conservation planning. Solving this equation yields the long-term average annual soil loss that can be expected from rill and inter-rill erosion. In this study, manual interpretation of color and color infrared 70 mm photography at the scale of 1:60,000 is used to determine the cropping management factor in the USLE. Accurate information was collected about plowing practices and crop residue cover (unharvested vegetation) for the winter season on agricultural land in Pheasant Branch Creek watershed in Dane County, Wisconsin.

Morgan, K. M.↗
NASA NTRS2009

Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions

A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.

plasma↗
NASA NTRS1993

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗
DOE OSTI2021

A two-phase three-field modeling framework for heat pipe application in nuclear reactors

Heat pipes and two-phase thermosyphons are highly efficient heat transfer devices utilizing continuous evaporation and condensation of working fluid for two-phase heat transport in closed systems. Because of the nearly isothermal and fully passive phase-change heat transfer mechanism, heat pipes and thermosyphons have found many applications in nuclear engineering, space technologies, and other energy systems. High-temperature heat pipes are used in nuclear microreactors to remove fission power from the primary system and are coupled with power conversion systems or process heat applications. Modeling of the two-phase flow phenomena inside a heat pipe is essential to its design and safety analysis. In this study, a comprehensive one-dimensional two-phase three-field flow model has been developed for the analysis of heat pipes in normal operation conditions and transients. The conservation or field equations of mass, momentum, and energy were developed for the liquid film, vapor, and droplet. In addition, constitutive models or correlations were reviewed thoroughly and provided for the closure of the three-field equations. Specific constitutive equations regarding interfacial mass and heat transfer at two interfaces, namely film-gas interface and gas-droplet interface, were reviewed for droplet entrainment and deposition rates as well as film and droplet evaporation rates. Furthermore, mechanistic correlations of annular flow film thickness were recommended for the modeling of the thermosyphons without a wick as a critical constitutive correlation. Furthermore, experimental data needs from new experiments using a prototype working fluid or surrogate fluids for the model validation of high-temperature heat pipes in microreactors were recommended for future research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A multi-phase-field model of topological pattern formation during electrochemical dealloying of binary alloys

Porous metal structures produced by electrochemical dealloying have been extensively studied for use in various applications that require high surface areas. A new comprehensive multi-phase-field (MPF) model is proposed to study topological porous patterns formed by spontaneously etching a bulk binary alloy that involves electrochemical reactions, bulk and surface diffusion, ion transport, applied electrode potential, and charge conservation. The governing equations for the alloy-porous cluster-electrolyte system account for a generalized Butler-Volmer electrochemical reaction and are in accordance with the classical nucleation theory. Based on a quantitative examination of the effects of electrode potential and precursor composition, the simulation results reproduce typical phenomena including passive surface dealloying, active porosity evolution, critical potential, and characteristic length scale in two-dimension (2-D) and three-dimension (3-D). To consider more complicated dealloying systems, the evolutions of a bimodal porous metal structure, nanocomposite, and nested porous network with a structured hierarchy are investigated. Finally, the proposed model can be a useful tool for understanding and predicting the morphology evolution of diverse porous structures during electrochemical dealloying.

Hierarchical porous structure↗
DOE OSTI2023

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1963

CONVECTION IN PLANETARY INTERIORS

Formulation of linearized equations safeguarding the conservation of mass, energy and momentum of viscous flow inside a fluid sphere, and application to the problem of convection in planetary interiors

VISCOUS FLOW↗
NASA NTRS1972

Numerical computation of multishocked, three-dimensional supersonic flow fields with real gas effects.

A computational procedure is presented which is capable of determining the supersonic flow field surrounding three-dimensional wing-body configurations such as a delta-wing space shuttle. The governing equations in conservation-law form are solved by a finite difference method using a second-order noncentered algorithm between the body and the outermost shock wave, which is treated as a sharp discontinuity. Secondary shocks which form between these boundaries are captured automatically, and the intersection of these shocks with the bow shock posed no difficulty. Resulting flow fields about typical blunt nose shuttle-like configurations at angle of attack are presented. The differences between perfect and real gas effects for high Mach number flows are shown.

Kutler, P.↗
NASA NTRS1974

A comparison of numerical solutions of the advective equation

Second- and third-order finite-difference methods recently applied to problems in high-speed fluid flow are applied to the model advection equation cast in conservative form. The differencing methods considered use forward time differencing with both centered and preferential space differences employed in the predictor-corrector sequences. The free parameter required for stability in the third-order method is adjusted to cause the solution to be either minimum dispersive or minimum dissipative in nature. Results indicate that the third-order method using minimum dispersion is the most accurate method tested. Computer time requirements are approximately twice those needed for second-order techniques.

Anderson, D.↗
NASA NTRS1975

Three-dimensional, shock-on-shock interaction problem

The unsteady, three-dimensional flowfield resulting from the interaction of a plane shock with a cone-shaped vehicle traveling supersonically is determined, using a second-order, shock-capturing, finite-difference approach. The time-dependent, inviscid gasdynamic equations are transformed to include the self-similar property of the flow, to align various coordinate surfaces with known shock waves, and to cluster points in the vicinity of the intersection of the transmitted incident shock and the surface of the vehicle. The governing partial differential equations in conservation-law form are then solved iteratively using MacCormack's (1969) algorithm.

Kutler, P.↗
NASA NTRS1976

A BBGKY framework for fluid turbulence

A framework is presented for a systematic kinetic theory of turbulence originating from the Liouville equation for the Fourier coefficients of fluid variables. The real and imaginary parts of these Fourier coefficients play the role that particle coordinates (positions and momenta) play in the BBGKY theory. The basic relations of the problem are the incompressible Navier-Stokes equations in two dimensions with zero viscosity, with the probability distributions of Fourier coefficients rather than moments being the basic variables of the theory. A kinetic equation is derived and shown to possess a number of requirements that any reasonable kinetic equation must have: conservation laws, positive-definite spectral densities, and an H-theorem. The major lack in the theory is any reliable information on the relaxation predicted by the complicated linear operator H. Closure of the hierarchy is achieved by the hypothesis that the five-coefficient correlation function is negligible. Problems associated with inclusion of viscosity and external driving forces are discussed.

Montgomery, D.↗
NASA NTRS1976

Computation of the inviscid supersonic flow over an external axial corner

A second-order finite-difference procedure is used to evaluate the inviscid supersonic flowfield surrounding an external axial corner composed of swept planar compression surfaces and representing the inlets on existing high-speed aircraft. The governing partial differential equations in conservation-law form are hyperbolic with respect to the axial coordinate and are solved iteratively by means of MacCormack's algorithm. The procedure treats both the peripheral shock wave and vortical singularities as discontinuities. Numerical results are presented for two parametric studies regarding the effects on the flowfield of varying the free-stream Mach number and the leading edge sweep of the horizontal wedge. Results of parametric Mach number study agree with the Mach number independence principle in that as the Mach number increases, such characteristics as shock shape, cross-flow sonic line location, and vortical singularity position approach an asymptote.

Kutler, P.↗