SEARCH · Engineering Papers

Results for “conservation form”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS2018

Uncertainty in Inverted Pendulum Thrust Measurements

An uncertainty analysis of a common configuration of electric propulsion thrust stand is presented. The analysis applies to inverted pendulum thrust stands operating in a null-coil configuration with in-situ calibration. Several sources of bias and precision uncertainty are discussed, propagated, and combined to form conservative estimates of the relative and absolute thrust uncertainties. A case study of the NASA Glenn Research Center Vacuum Facility 6 thrust stand is presented. For the thruster investigated, the uncertainty was estimated to be Â±6.9mN over the entire span of thrust. This uncertainty represents a maximum instrument bias introduced by the thrust stand. The paper does not address repeatability of actual thrust measurements, as this is generally beyond the influence of the thrust stand and can be dependent on a large number of factors.

Mackey, Jonathan A.↗

NASA NTRS2011

Discretely Conservative Finite-Difference Formulations for Nonlinear Conservation Laws in Split Form: Theory and Boundary Conditions

Simulations of nonlinear conservation laws that admit discontinuous solutions are typically restricted to discretizations of equations that are explicitly written in divergence form. This restriction is, however, unnecessary. Herein, linear combinations of divergence and product rule forms that have been discretized using diagonal-norm skew-symmetric summation-by-parts (SBP) operators, are shown to satisfy the sufficient conditions of the Lax-Wendroff theorem and thus are appropriate for simulations of discontinuous physical phenomena. Furthermore, special treatments are not required at the points that are near physical boundaries (i.e., discrete conservation is achieved throughout the entire computational domain, including the boundaries). Examples are presented of a fourth-order, SBP finite-difference operator with second-order boundary closures. Sixth- and eighth-order constructions are derived, and included in E. Narrow-stencil difference operators for linear viscous terms are also derived; these guarantee the conservative form of the combined operator.

Fisher, Travis C.↗

NASA NTRS1986

On the Maxwellian distribution, symmetric form, and entropy conservation for the Euler equations

The Euler equations of gas dynamics have some very interesting properties in that the flux vector is a homogeneous function of the unknowns and the equations can be cast in symmetric hyperbolic form and satisfy the entropy conservation. The Euler equations are the moments of the Boltzmann equation of the kinetic theory of gases when the velocity distribution function is a Maxwellian. The present paper shows the relationship between the symmetrizability and the Maxwellian velocity distribution. The entropy conservation is in terms of the H-function, which is a slight modification of the H-function first introduced by Boltzmann in his famous H-theorem. In view of the H-theorem, it is suggested that the development of total H-diminishing (THD) numerical methods may be more profitable than the usual total variation diminishing (TVD) methods for obtaining wiggle-free solutions.

Deshpande, S. M.↗

NASA NTRS1978

Coefficient matrices for implicit finite difference solution of the inviscid fluid conservation law equations

Although the Navier-Stokes equations describe most flows of interest in aerodynamics, the inviscid conservation law equations may be used for small regions with viscous forces. Thus, Euler equations and several time-accurate finite difference procedures, explicit and implicit, are discussed. Although implicit techniques require more computational work, they permit larger time steps to be taken without instability. It is noted that the Jacobian matrices for Euler equations in conservation-law form have certain eigenvalue-eigenvector properties which may be used to construct conservative-form coefficient matrices. This reduces the computation time of several implicit and semiimplicit schemes. Extensions of the basic approach to other areas are suggested.

Steger, J. L.↗

NASA NTRS1991

Choice of velocity variables for complex flow computation

The issue of adopting the velocity components as dependent velocity variables for the Navier-Stokes flow computations is investigated. The viewpoint advocated is that a numerical algorithm should preferably honor both the physical conservation law in differential form and the geometric conservation law in discrete form. With the use of Cartesian velocity vector, the momentum equations in curvilinear coordinates can retain the full conservation-law form and satisfy the physical conservation laws. With the curvilinear velocity components, source terms appear in differential equations and hence the full conservation law form can not be retained. In discrete expressions, algorithms based on the Cartesian components can satisfy the geometric conservation-law form for convection terms but not for viscous terms; those based on the curvilinear components, on the other hand, cannot satisfy the geometric conservation-law form for either convection or viscous terms. Several flow solutions for domain with 90 and 360 degree turnings are presented to illustrate the issues of using the Cartesian velocity components and the staggered grid arrangement.

Shyy, W.↗

NASA NTRS1977

An implicit factored scheme for the compressible Navier-Stokes equations

An implicit finite difference scheme is developed for the numerical solution of the compressible Navier-Stokes equations in conservation-law form. The algorithm is second-order-time accurate, noniterative, and spatially factored. In order to obtain an efficient factored algorithm, the spatial cross-derivatives are evaluated explicitly. However, the algorithm is unconditionally stable and, although a three-time-level scheme, requires only two-time-levels of data storage. The algorithm is constructed in a 'delta' form (i.e., increments of the conserved variables and fluxes) that provides a direct derivation of the scheme and leads to an efficient computational algorithm. In addition, the delta form has the advantageous property of a steady-state (if one exists) independent of the size of the time step. Numerical results are presented for a two-dimensional shock boundary-layer interaction problem.

Beam, R. M.↗

NASA NTRS1981

Computation of wing-vortex interaction in transonic flow using implicit finite difference algorithm

An implicit delta form finite difference algorithm for Euler equations in conservation law form was used in preliminary calculations of three dimensional wing vortex interaction. Both steady and unsteady transonic flow wing vortex interactions are computed. The computations themselves are meant to guide upcoming wind tunnel experiments of the same flow field. Various modifications to the numerical method that are intended to improve computational efficiency are also described and tested in both two and three dimensions. Combination of these methods can reduce the overall computational time by a factor of 4.

Srinivasan, G.↗

NASA NTRS1981

Computation of simple three-dimensional wing-vortex interaction in transonic flow

An implicit delta form finite-difference algorithm for Euler equations in conservation law form has been used in preliminary calculations of three-dimensional wing-vortex interactions. Both steady and unsteady transonic flow wing-vortex interactions are computed. The computations themselves are meant to guide upcoming wind tunnel experiments of the same flow field. Various modifications to the numerical method that are intended to improve computational efficiency are also described and tested in both two- and three-dimensions.

Srinivasan, G. R.↗

NASA NTRS2001

High Order Difference Method for Low Mach Number Aeroacoustics

A high order finite difference method with improved accuracy and stability properties for computational aeroacoustics (CAA) at low Mach numbers is proposed. The Euler equations are split into a conservative and a symmetric non- conservative portion to allow the derivation of a generalized energy estimate. Since the symmetrization is based on entropy variables, that splitting of the flux derivatives is referred to as entropy splitting. Its discretization by high order central differences was found to need less numerical dissipation than conventional conservative schemes. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the split Euler equations are formulated in perturbation form. The unknowns are the small changes of the conservative variables with respect to their large stagnation values. All nonlinearities and the conservation form of the conservative portion of the split flux derivatives can be retained, while cancellation errors are avoided with its discretization opposed to the conventional conservative form. The finite difference method is third-order accurate at the boundary and the conventional central sixth-order accurate stencil in the interior. The difference operator satisfies the summation by parts property analogous to the integration by parts in the continuous energy estimate. Thus, strict stability of the difference method follows automatically. Spurious high frequency oscillations are suppressed by a characteristic-based filter similar to but without limiter. The time derivative is approximated by a 4-stage low-storage second-order explicit Runge-Kutta method. The method has been applied to simulate vortex sound at low Mach numbers. We consider the Kirchhoff vortex, which is an elliptical patch of constant vorticity rotating with constant angular frequency in irrotational flow. The acoustic pressure generated by the Kirchhoff vortex is governed by the 2D Helmholtz equation, which can be solved analytically using separation of variables.

Mueller, B.↗

NASA NTRS1994

Theoretical studies of a molecular beam generator

We consider the problem of modeling a high temperature diatomic gas N2 flowing through a converging-diverging high expansion nozzle. The problem is modeled in two ways. The first model uses a single temperature with variable specific heats as functions of this temperature. For the second model we assume that the various degrees of freedom all have a Boltzmann population distribution which means that each degree of freedom has its own temperature and consequently each system state can be characterized by these temperatures. This suggests the formulation of a second model with a vibrational degree of freedom as having its own temperature along with a rotational-translational degree of freedom with its own temperature. Initially the vibrational degree of freedom is excited by heating the gas to a high temperature. As the high temperature gas expands through the nozzle throat there is a sudden drop in the rotational-translational temperature along with a finite relaxation time for the vibrational degree of freedom to achieve equilibrium with the rotational-translational degree of freedom. That is, the temperature change that occurs when the N2 gas passes through the nozzle throat is so great that the changes in the vibrational degree of freedom lags behind the rotational-translational energy changes. The resulting relaxation time is finite. It is in this context that the term nonequilibrium is used. That is, the term nonequilibrium denotes the fact that the energy content of the various degrees of freedom are characterized by two temperatures. We neglect any chemical reactions resulting from the high temperatures which could also add nonequilibrium effects. We develop the basic equations for the two models in various forms in order to check the derivations with other sources. The final form which is solved numerically consists of the scaled equations in a conservative dimensionless form.

Heinbockel, John H.↗

NASA NTRS1984

Approximate-factorization algorithms - Theory and applications in viscous-flow computations

A systematic development of implicit approximate-factorization algorithms in delta form for both unsteady and steady viscous flow is presented. The algorithms are cast in conservation-law form and simplified by using a thin-layer approximation to the governing equations. The implementation of implicit surface viscous boundary conditions is discussed in detail, and an example is presented illustrating the advantage of using the implicit boundary conditions. Three-dimensional results from the steady form of the algorithm are presented and compared with experimental data.

Chaussee, D. S.↗

NASA NTRSDate not supplied

Impact of Anthropogenic Activity, Climate Change, and Urbanization on Wetland Habitat in the Platte River Basin

The Platte River Basin (PRB) is a dynamic ecosystem where wetlands play a pivotal role as essential habitats for various flora and fauna, including local and migratory birds. It provides many crucial ecosystem services that benefit humans directly and indirectly. However, anthropogenic activity, climate change, and urbanization have resulted in decline in wildlife habitat, elevated flood risk, and wetland loss. To address this issue, NASA DEVELOP partnered with Audubon Great Plains (AGP) to address the vital habitats within urban areas to protect bird species, reduce flood hazards, and analyze the potential impact of future development on wetlands. We utilized remotely sensed data from Landsat 8 Operational Land Imager (OLI), Sentinel-2 Multispectral Instrument (MSI), and Sentinel-1 Synthetic Aperture Radar (SAR) to assess land use and land cover (LULC) change. Nighttime lights data from Suomi National Polar-orbiting Partnership (NPP) Visible Infrared Imaging Radiometer Suite (VIIRS) as well as NASA Socioeconomic Data and Applications Center (SEDAC) population data were also used as inputs to simulate urban growth potential up to 2050 using the open-source FUTure Urban-Regional Environment Simulation (FUTURES) model. A broad scale analysis across 13 focal cities showed varied changes in land use patterns across the PRB, with the most notable being a decrease in agricultural land coverage and an increase in vegetation and grassland coverage. We overlaid a flood extent map with the LULC classifications in Grand Island to identify possible restoration sites under AGP’s Urban Woods and Prairies Initiative. The results for two proposed scenarios showed that at least 51 counties out of 81 in the PRB would experience growth by year 2050. The first scenario (all wetlands are protected) showed that there will be no loss of wetlands by 2050. However, the second scenario (no wetlands are protected) showed a decrease in wetland area and loss of habitat for bird conservation. The results will help AGP to lead awareness workshops for communities about wetland protection and to form impactful conservation strategies in the future.

Nancee Uniyal↗

NASA NTRS2024

Platte River Basin Water Resources II: Predicting Land Cover Change in the Platte River Basin to Select Wetland Protection Sites Vulnerable to Urban Encroachment

The Platte River Basin (PRB) represents a critical dynamic ecosystem where wetlands play a pivotal role as essential habitats for local and migratory birds, as well as various flora and fauna. It also provides many crucial ecosystem services that directly and indirectly benefit human welfare. However, it is threatened by anthropogenic activity, climate change, and urbanization. Our team partnered with Audubon Great Plains, the regional office of a national non-profit organization (NGO), to analyze the future potential development of the region and its potential impact on the wetlands. We utilized remotely sensed data, such as Landsat 8 Operational Land Imager (OLI) and Suomi National Polar Orbiting Partnership (NPP) Visible Infrared Imaging Radiometer Suite (VIIRS), as well as NASA Socioeconomic Data and Application Center (SEDAC) data, to simulate urban growth potential up to 2050 using the open-source model FUTure Urban-Regional Environment Simulation (FUTURES) model. Our results for two proposed scenarios (all wetlands are protected, and no wetlands are protected) revealed that at least 56 counties out of 81 in the PRB would experience growth by the year 2050. The results for the first scenario show that there will be no loss of wetlands in the future. However, the results for the second scenario indicate a basin-wide reduction of 43.22 square kilometers (0.58%) in wetland areas, leading to a significant loss of bird habitat, critical for conservation. The results will help Audubon Great Plains' Urban Woods and Prairies Initiative to lead awareness workshops for communities about wetland protection and to form impactful conservation strategies.

FUTURES↗

NASA NTRS1976

Implicit finite-difference procedures for the computation of vortex wakes

Implicit finite-difference procedures for the primitive form of the incompressible Navier-Stokes and the compressible Euler equations are used to compute vortex wake flows. The partial differential equations in strong conservation-law form are transformed to cluster grid points in regions with large changes in vorticity. In addition to clustering, fourth-order accurate, spatial difference operators are used to help resolve the flow-field gradients. The use of implicit time-differencing permits large time steps to be taken since temporal variations are typically small. Computational efficiency is achieved by approximate factorization. Both two-dimensional and preliminary three-dimensional calculations are described and qualitatively compared with existing experimental data.

Steger, J. L.↗

NASA NTRS1983

Convair 990 transonic flow-field simulation about the forward fuselage

A three-dimensional, viscous flow code was used to calculate the transonic flow about the forebody of the Convair CV-990 (Galileo II) research aircraft stationed at NASA Ames Research Center. The computations were used to determine the location for a differential pressure system. In addition, attitude sensor placements were verified. These instruments comprise a meteorological measurement system, which will be used for global determination of three-dimensional wind data. The code solves the thin layer form of the Reynolds-averaged Navier-Stokes equations using an implicit numerical procedure. The governing equations are written in a generalized, nonorthogonal coordinate system, and are cast in a strong conservation law form. Laminar boundary layer results are presented for free stream Mach number of 0.8 and angles of attack of zero and 2 deg. Use of this computational tool reduced the development time for the location of the sensors and aided in the optimal placement on the aircraft of these devices.

Chaussee, D. S.↗

NASA NTRS1989

Algorithm For Hypersonic Flow In Chemical Equilibrium

Implicit, finite-difference, shock-capturing algorithm calculates inviscid, hypersonic flows in chemical equilibrium. Implicit formulation chosen because overcomes limitation on mathematical stability encountered in explicit formulations. For dynamical portion of problem, Euler equations written in conservation-law form in Cartesian coordinate system for two-dimensional or axisymmetric flow. For chemical portion of problem, equilibrium state of gas at each point in computational grid determined by minimizing local Gibbs free energy, subject to local conservation of molecules, atoms, ions, and total enthalpy. Major advantage: resulting algorithm naturally stable and captures strong shocks without help of artificial-dissipation terms to damp out spurious numerical oscillations.

Palmer, Grant↗

NASA NTRS1992

Numerical simulation of conservation laws

A new numerical framework for solving conservation laws is being developed. This new approach differs substantially from the well established methods, i.e., finite difference, finite volume, finite element and spectral methods, in both concept and methodology. The key features of the current scheme include: (1) direct discretization of the integral forms of conservation laws, (2) treating space and time on the same footing, (3) flux conservation in space and time, and (4) unified treatment of the convection and diffusion fluxes. The model equation considered in the initial study is the standard one dimensional unsteady constant-coefficient convection-diffusion equation. In a stability study, it is shown that the principal and spurious amplification factors of the current scheme, respectively, are structurally similar to those of the leapfrog/DuFort-Frankel scheme. As a result, the current scheme has no numerical diffusion in the special case of pure convection and is unconditionally stable in the special case of pure diffusion. Assuming smooth initial data, it will be shown theoretically and numerically that, by using an easily determined optimal time step, the accuracy of the current scheme may reach a level which is several orders of magnitude higher than that of the MacCormack scheme, with virtually identical operation count.

Chang, Sin-Chung↗

NASA NTRS1981

Numerical approximation of boundary conditions with applications to inviscid equations of gas dynamics

A comprehensive overview of the state of the art of well-posedness and stability analysis of difference approximations for initial boundary value problems of the hyperbolic type is presented. The applicability of recent theoretical development to practical calculations for nonlinear gas dynamics is examined. The one dimensional inviscid gas dynamics equations in conservation law form are selected for numerical experiments. The class of implicit schemes developed from linear multistep methods in ordinary differential equations is chosen and the use of linear extrapolation as an explicit or implicit boundary scheme is emphasized. Specification of boundary data in the primitive variables and computation in terms of the conservative variables in the interior is discussed. Some numerical examples for the quasi-one-dimensional nozzle are given.

Yee, H. C.↗