Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conservative”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

An adaptive mesh-moving and refinement procedure for one-dimensional conservation laws

We examine the performance of an adaptive mesh-moving and /or local mesh refinement procedure for the finite difference solution of one-dimensional hyperbolic systems of conservation laws. Adaptive motion of a base mesh is designed to isolate spatially distinct phenomena, and recursive local refinement of the time step and cells of the stationary or moving base mesh is performed in regions where a refinement indicator exceeds a prescribed tolerance. These adaptive procedures are incorporated into a computer code that includes a MacCormack finite difference scheme wih Davis' artificial viscosity model and a discretization error estimate based on Richardson's extrapolation. Experiments are conducted on three problems in order to qualify the advantages of adaptive techniques relative to uniform mesh computations and the relative benefits of mesh moving and refinement. Key results indicate that local mesh refinement, with and without mesh moving, can provide reliable solutions at much lower computational cost than possible on uniform meshes; that mesh motion can be used to improve the results of uniform mesh solutions for a modest computational effort; that the cost of managing the tree data structure associated with refinement is small; and that a combination of mesh motion and refinement reliably produces solutions for the least cost per unit accuracy.

Biswas, Rupak↗
NASA NTRS1994

UARS MILS O3 soundings compared with lidar measurements using the conservative coordinates reconstruction technique

A technique based on conservative properties of certain meteorological fields is used to compare ozone measurements from the Microwave Limb Sounder (MLS) on the Upper Atmosphere Research Satellite (UARS) with soundings from a lidar system operated at midlatitudes by the University of L'Aquila, Italy. A few typical cases are analyzed in connection with the position of the vortex relative to the observing station, and it is shown that in general lidar observations taken within the vortex compare well with the UARS data, regardless of whether they are coincident with a satellite overpass. It is shown that such analysis may be useful for comparing measurements of the same quantity taken at different sites using different measurement techniques.

Redaelli, G.↗
NASA NTRS1995

Optimization methods, flux conserving methods for steady state Navier-Stokes equation

Navier-Stokes equation as discretized by new flux conserving method proposed by Chang and Scott results in the system: vector F(vector x) = 0, where F is a vector valued function. The Optimization method we use is based on Quasi-Newton methods: given a nonlinear function vector F(vector x) = 0, we solve, Delta(vector x) = -BF(vector x), where Delta(vector x) is the correction term and B is the inverse Jacobian of F(x). Then, iteratively, vector(x(sub (i+1))) = vector(x (sub i)) + alpha.Delta(vector x(sub i)), where alpha is a line search correction term determined by a line search routine. We use the BFCG's update the Jacobian matrix B(sub k) at each iteration. It is well known that B(sub k) approaches B(*) at the solution X(*). This algorithm has several advantages over the Newton-Raphson method. For example, we do not need to calculate the Jacobian matrix at each iteration which is computationally very expensive.

Adeyeye, John↗
NASA NTRS1995

Nonlinear (time domain) and linearized (time and frequency domain) solutions to the compressible Euler equations in conservation law form

Two linearized solvers (time and frequency domain) based on a high resolution numerical scheme are presented. The basic approach is to linearize the flux vector by expressing it as a sum of a mean and a perturbation. This allows the governing equations to be maintained in conservation law form. A key difference between the time and frequency domain computations is that the frequency domain computations require only one grid block irrespective of the interblade phase angle for which the flow is being computed. As a result of this and due to the fact that the governing equations for this case are steady, frequency domain computations are substantially faster than the corresponding time domain computations. The linearized equations are used to compute flows in turbomachinery blade rows (cascades) arising due to blade vibrations. Numerical solutions are compared to linear theory (where available) and to numerical solutions of the nonlinear Euler equations.

Sreenivas, Kidambi↗
NASA NTRS1995

An hp-adaptivity and error estimation for hyperbolic conservation laws

This paper presents an hp-adaptive discontinuous Galerkin method for linear hyperbolic conservation laws. A priori and a posteriori error estimates are derived in mesh-dependent norms which reflect the dependence of the approximate solution on the element size (h) and the degree (p) of the local polynomial approximation. The a posteriori error estimate, based on the element residual method, provides bounds on the actual global error in the approximate solution. The adaptive strategy is designed to deliver an approximate solution with the specified level of error in three steps. The a posteriori estimate is used to assess the accuracy of a given approximate solution and the a priori estimate is used to predict the mesh refinements and polynomial enrichment needed to deliver the desired solution. Numerical examples demonstrate the reliability of the a posteriori error estimates and the effectiveness of the hp-adaptive strategy.

Bey, Kim S.↗
NASA NTRS1995

The method of space-time and conservation element and solution element: A new approach for solving the Navier-Stokes and Euler equations

A new numerical framework for solving conservation laws is being developed. This new framework differs substantially in both concept and methodology from the well-established methods, i.e., finite difference, finite volume, finite element, and spectral methods. It is conceptually simple and designed to overcome several key limitations of the above traditional methods. A two-level scheme for solving the convection-diffusion equation is constructed and used to illuminate the major differences between the present method and those previously mentioned. This explicit scheme, referred to as the a-mu scheme, has two independent marching variables.

Chang, Sin-Chung↗
NASA NTRS1995

Conservative properties of finite difference schemes for incompressible flow

The purpose of this research is to construct accurate finite difference schemes for incompressible unsteady flow simulations such as LES (large-eddy simulation) or DNS (direct numerical simulation). In this report, conservation properties of the continuity, momentum, and kinetic energy equations for incompressible flow are specified as analytical requirements for a proper set of discretized equations. Existing finite difference schemes in staggered grid systems are checked for satisfaction of the requirements. Proper higher order accurate finite difference schemes in a staggered grid system are then proposed. Plane channel flow is simulated using the proposed fourth order accurate finite difference scheme and the results compared with those of the second order accurate Harlow and Welch algorithm.

Morinishi, Youhei↗
NASA NTRS1996

A Conservation Strategy for the Florida Scrub-Jay on John F. Kennedy Space Center/Merritt Island National Wildlife Refuge: An Initial Scientific Basis for Recovery

The Florida Scrub-Jay (Aphelocoma coerulescens) is an indicator of ecosystem integrity of Florida scrub, an endangered ecosystem that requires frequent fire. One of the largest populations of this federally threatened species occurs on John F. Kennedy Space Center/Merritt Island National Wildlife Refuge. Population trends were predicted using population modeling and field data on reproduction and survival of Florida Scrub-Jays collected from 1988 - 1995. Analyses of historical photography indicated that habitat suitability has been declining for 30 years. Field data and computer simulations suggested that the population declined by at least 40% and will decline by another 40% in 1 0 years, if habitat management is not greatly intensified. Data and computer simulations suggest that habitat suitability cannot deviate greatly from optimal for the jay population to persist. Landscape trajectories of vegetation structure, responsible for declining habitat suitability, are associated with the disruption of natural fire regimes. Prescribed fire alone can not reverse the trajectories. A recovery strategy was developed, based on studies of Florida Scrub-Jays and scrub vegetation. A reserve design was formulated based on conservation science principles for scrub ecosystems. The strategy emphasizes frequent fire to restore habitat, but includes mechanical tree cutting for severely degraded areas. Pine thinning across large areas can produce rapid increases in habitat quality. Site-specific strategies will need to be developed, monitored, and modified to achieve conditions suitable for population persistence.

Breininger, D. R.↗
NASA NTRS1993

Theory of Metastable State Relaxation in a Gravitational Field for Non-Critical Binary Systems with Non-Conserved Order Parameter

A new mathematical ansatz is developed for solution of the time-dependent Ginzburg-Landau nonlinear partial differential equation describing metastable state relaxation in binary (solute+solvent) non-critical solutions with non-conserved scalar order parameter in presence of a gravitational field. It has been demonstrated analytically that in such systems metastability initiates heterogeneous solute redistribution which results in the formation of a non-equilibrium singly-periodic spatial solute structure in the new solute-rich phase. The critical radius of nucleation and the induction time in these systems are gravity-dependent. It has also been proved that metastable state relaxation in vertical columns of supersaturated non-critical binary solutions leads to formation of the solute concentration gradient. Analytical expression for this concentration gradient is found and analysed. It is concluded that gravity can initiate phase separation (nucleation or spinodal decomposition).

Izmailov, Alexander F.↗
NASA NTRS1993

Theory of Metastable State Relaxation for Non-Critical Binary Systems with Non-Conserved Order Parameter

A new mathematical ansatz for a solution of the time-dependent Ginzburg-Landau non-linear partial differential equation is developed for non-critical systems such as non-critical binary solutions (solute + solvent) described by the non-conserved scalar order parameter. It is demonstrated that in such systems metastability initiates heterogeneous solute redistribution which results in formation of the non-equilibrium singly-periodic spatial solute structure. It is found how the time-dependent period of this structure evolves in time. In addition, the critical radius r(sub c) for solute embryo of the new solute rich phase together with the metastable state lifetime t(sub c) are determined analytically and analyzed.

Izmailov, Alexander↗
NASA NTRS1997

Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws

In these lecture notes we describe the construction, analysis, and application of ENO (Essentially Non-Oscillatory) and WENO (Weighted Essentially Non-Oscillatory) schemes for hyperbolic conservation laws and related Hamilton- Jacobi equations. ENO and WENO schemes are high order accurate finite difference schemes designed for problems with piecewise smooth solutions containing discontinuities. The key idea lies at the approximation level, where a nonlinear adaptive procedure is used to automatically choose the locally smoothest stencil, hence avoiding crossing discontinuities in the interpolation procedure as much as possible. ENO and WENO schemes have been quite successful in applications, especially for problems containing both shocks and complicated smooth solution structures, such as compressible turbulence simulations and aeroacoustics. These lecture notes are basically self-contained. It is our hope that with these notes and with the help of the quoted references, the reader can understand the algorithms and code them up for applications.

Shu, Chi-Wang↗
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 NTRS1998

A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy

Stable and accurate interface conditions are derived for the linear advection-diffusion equation. The conditions are functionally independent of the spatial order of accuracy and rely only on the form of the discrete operator. We focus on high-order finite-difference operators that satisfy the summation-by-parts (SBP) property. We prove that stability is a natural consequence of the SBP operators used in conjunction with the new boundary conditions. In addition, we show that the interface treatments are conservative. New finite-difference operators of spatial accuracy up to sixth order are constructed: these operators satisfy the SBP property. Finite-difference operators are shown to admit design accuracy (p(sup th)-order global accuracy) when (p - 1)(sup th)-order stencil closures are used near the boundaries if the physical boundary conditions are implemented to at least p(sup th)-order accuracy. Stability and accuracy are demonstrated on the nonlinear Burgers' equation for an twelve-subdomain problem with randomly distributed interfaces.

Carpenter, Mark H.↗
NASA NTRS1996

Evidence of Momentum Conservation at a Nonepitaxial Metal/Semiconductor Interface Using Ballistic Electron Emission Microscopy

Ballistic-Electron-Emission Microscopy (BEEM) spectroscopy has been performed on Au/Si(111) structures as a function of Au thickness and temperature. At 77 K a direct signature of parallel momentum conservation at the Au/Si interface is observed in the BEEM spectra. The variation in spectral shape with both Au thickness and temperature places restrictions on allowable values of inelastic and elastic mean-free paths in the metal, and also requires the presence of multiple electron passes within the Au layer. An independent indication of multiple reflections is directly observed in the attenuation of BEEM current with Au thickness.

Bell, L. D.↗
NASA NTRS1981

Fuel-Conservation Guidance System for Powered-Lift Aircraft

A technique is described for the design of fuel-conservative guidance systems and is applied to a system that was flight tested on board NASA's sugmentor wing jet STOL research aircraft. An important operational feature of the system is its ability to rapidly synthesize fuel-efficient trajectories for a large set of initial aircraft positions, altitudes, and headings. This feature allows the aircraft to be flown efficiently under conditions of changing winds and air traffic control vectors. Rapid synthesis of fuel-efficient trajectories is accomplished in the airborne computer by fast-time trajectory integration using a simplified dynamic performance model of the aircraft. This technique also ensures optimum flap deployment and, for powered-lift STOL aircraft, optimum transition to low-speed flight. Also included in the design is accurate prediction of touchdown time for use in four-dimensional guidance applications. Flight test results have demonstrated that the automatically synthesized trajectories produce significant fuel savings relative to manually flown conventional approaches.

Erzberger, Heinz↗
NASA NTRS1980

Fast, Conservative Algorithm for Solving the Transonic Full-Potential Equation

A fast, fully implicit approximate factorization algorithm designed to solve the conservative, transonic, full-potential equation in either two or three dimensions is described. The algorithm uses an upwind bias of the density coefficient for stability in supersonic regions. This provides an effective upwind difference of the streamwise terms for any orientation of the velocity vector (i.e., rotated differencing), thereby greatly enhancing the reliability of the present algorithm. A numerical transformation is used to establish an arbitrary body-fitted, finite-difference mesh. Computed results for both airfoils and simplified wings demonstrate substantial improvement in convergence speed for the new algorithm relative to standard successive-line over-relaxation algorithms.

Holst, Terry L.↗
NASA NTRS1999

Application of the Space-Time Conservation Element and Solution Element Method to One-Dimensional Advection-Diffusion Problems

Test problems are used to examine the performance of several one-dimensional numerical schemes based on the space-time conservation and solution element (CE/SE) method. Investigated in this paper are the CE/SE schemes constructed previously for solving the linear unsteady advection-diffusion equation and the schemes derived here for solving the nonlinear viscous and inviscid Burgers equations. In comparison with the numerical solutions obtained using several traditional finite-difference schemes with similar accuracy, the CE/SE solutions display much lower numerical dissipation and dispersion errors.

Wang, Xiao-Yen↗