SEARCH · Engineering Papers

Results for “consistent scheme”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Consistent boundary conditions for Reduced Navier-Stokes (RNS) scheme applied to 3-dimensional internal viscous flows

A consistent and efficient set of boundary conditions is developed for the multi-sweep space-marching pressure-elliptic Reduced Navier-Stokes (RNS) scheme as applied for 3-D internal viscous flow problems. No-slip boundary conditions are directly imposed on the solid walls. There is no iteration procedure required in the cross plane to ensure mass conservation across each marching plane. The finite difference equations forming the coefficient matrix are ordered such that the surface normal velocity is specified on all the solid walls; unlike external flows, a pressure boundary condition in the cross plane is not required. Since continuity is directly satisfied at all points in the flow domain, the first order momentum equations can be solved directly for the pressure without the need for a Poisson pressure correction equation. The procedure developed herein can also be applied with periodic boundary conditions. The analysis is given for general compressible flows. Incompressible flow solutions are obtained, for straight and curved ducts of square cross section, to validate the procedure. These solutions are used to demonstrate the applicability of the RNS scheme, with the improved boundary conditions for internal flows with strong interaction, as would be encountered in ducts and turbomachinery geometries.

Reddy, D. R.↗

NASA NTRS1978

An old integration scheme for compressible flows revisited, refurbished and put to work

A scheme for integrating the Euler equations of compressible flow in any hyperbolic case is presented. The scheme relies on the concept of characteristics but is strictly a finite difference scheme. Improvements in accuracy and physical consistence due to the scheme are discussed and results of its application to complex flows are shown.

Moretti, G.↗

NASA NTRS2016

Comparison of NACA 0012 Laminar Flow Solutions: Structured and Unstructured Grid Methods

In this paper we consider the solution of the compressible Navier-Stokes equations for a class of laminar airfoil flows. The principal objective of this paper is to demonstrate that members of this class of laminar flows have steady-state solutions. These laminar airfoil flow cases are often used to evaluate accuracy, stability and convergence of numerical solution algorithms for the Navier-Stokes equations. In recent years, such flows have also been used as test cases for high-order numerical schemes. While generally consistent steady-state solutions have been obtained for these flows using higher order schemes, a number of results have been published with various solutions, including unsteady ones. We demonstrate with two different numerical methods and a range of meshes with a maximum density that exceeds 8 × 106 grid points that steady-state solutions are obtained. Furthermore, numerical evidence is presented that even when solving the equations with an unsteady algorithm, one obtains steady-state solutions.

Swanson, R. C.↗

NASA NTRS1983

An implicit, transonic, full-potential code for cascade flow on H-grid topology

A transonic, full-potential code is developed for computing the flow through two-dimensional cascades using an H-type grid topology that employs an implicit approximate-factorization scheme. The body-conforming H-grid is generated numerically by solving Poisson's equation. The flow-solution algorithm at the coordinate mapping singularity associated with this grid is investigated using two different types of finite-difference schemes. The grid-geometry effect on these schemes is also studied by noting free-stream capturing properties. It is found that by implementing a consistent spatial differencing scheme, the mapping singularities can be resolved numerically, and the grid-geometry-induced error minimized. The code is verified by computing model cascade flow problems.

Kwak, D.↗

NASA NTRS2017

Sensitivity of CONUS Summer Rainfall to the Selection of Cumulus Parameterization Schemes in NU-WRF Seasonal Simulations

This study investigates the sensitivity of daily rainfall rates in regional seasonal simulations over the contiguous United States (CONUS) to different cumulus parameterization schemes. Daily rainfall fields were simulated at 24-km resolution using the NASA-Unified Weather Research and Forecasting (NU-WRF) Model for June-August 2000. Four cumulus parameterization schemes and two options for shallow cumulus components in a specific scheme were tested. The spread in the domain-mean rainfall rates across the parameterization schemes was generally consistent between the entire CONUS and most subregions. The selection of the shallow cumulus component in a specific scheme had more impact than that of the four cumulus parameterization schemes. Regional variability in the performance of each scheme was assessed by calculating optimally weighted ensembles that minimize full root-mean-square errors against reference datasets. The spatial pattern of the seasonally averaged rainfall was insensitive to the selection of cumulus parameterization over mountainous regions because of the topographical pattern constraint, so that the simulation errors were mostly attributed to the overall bias there. In contrast, the spatial patterns over the Great Plains regions as well as the temporal variation over most parts of the CONUS were relatively sensitive to cumulus parameterization selection. Overall, adopting a single simulation result was preferable to generating a better ensemble for the seasonally averaged daily rainfall simulation, as long as their overall biases had the same positive or negative sign. However, an ensemble of multiple simulation results was more effective in reducing errors in the case of also considering temporal variation.

Iguchi, Takamichi↗

NASA NTRS1997

The Implicit and Explicit alpha-mu Schemes

Artificial numerical dissipation is an important issue in large Reynolds number computations. In such computations, the artificial dissipation inherent in traditional numerical schemes can overwhelm the physical dissipation and yield inaccurate results on meshes of practical size. In the present work, the space-time conservation element and solution element method is used to construct new and accurate numerical schemes such that artificial numerical dissipation will not overwhelm physical dissipation. Specifically, these schemes have the property that numerical dissipation vanishes when the physical viscosity goes to zero. These new schemes therefore accurately model the physical dissipation even when it is extremely small. The method of space-time conservation element and solution element, currently under development, is a nontraditional numerical method for solving conservation laws. The method is developed on the basis of local and global flux conservation in a space-time domain, in which space and time are treated in a unified manner. Explicit solvers for model and fluid dynamic conservation laws have previously been investigated. In this paper, we introduce a new concept in the design of implicit schemes, and use it to construct two highly accurate solvers for a convection-diffusion equation. The two schemes become identical in the pure convection case, and in the pure diffusion case. The implicit schemes are applicable over the whole Reynolds number range, from purely diffusive equations to purely inviscid (convective) equations. The stability and consistency of the schemes are analyzed, and some numerical results are presented. It is shown that, in the inviscid case, the new schemes become explicit and their amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, their principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme. We also construct an explicit solver with the treatment of diffusion being based on that in the implicit solvers. The explicit solver has only a CFL stability limitation on the Courant number, yet it retains the second-order spatial accuracy of the implicit schemes.

Chang, Sin-Chung↗

NASA NTRS1999

Implicit Space-Time Conservation Element and Solution Element Schemes

Artificial numerical dissipation is in important issue in large Reynolds number computations. In such computations, the artificial dissipation inherent in traditional numerical schemes can overwhelm the physical dissipation and yield inaccurate results on meshes of practical size. In the present work, the space-time conservation element and solution element method is used to construct new and accurate implicit numerical schemes such that artificial numerical dissipation will not overwhelm physical dissipation. Specifically, these schemes have the property that numerical dissipation vanishes when the physical viscosity goes to zero. These new schemes therefore accurately model the physical dissipation even when it is extremely small. The new schemes presented are two highly accurate implicit solvers for a convection-diffusion equation. The two schemes become identical in the pure convection case, and in the pure diffusion case. The implicit schemes are applicable over the whole Reynolds number range, from purely diffusive equations to convection-dominated equations with very small viscosity. The stability and consistency of the schemes are analysed, and some numerical results are presented. It is shown that, in the inviscid case, the new schemes become explicit and their amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, their principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme.

Chang, Sin-Chung↗

NASA NTRS1981

Basic advances in the finite-volume method for transonic potential flow calculations

The finite-volume method of Jameson and Caughey provides a framework within which it is possible to calculate transonic potential flows about essentially arbitrary geometrical configurations. Improvements designed to increase the accuracy of the basic scheme and its consistency in the far field will be described. These include the incorporation of an artificial viscosity which maintains the formal second-order accuracy of the scheme in supersonic zones, and a modification of the flux balances to allow the free-stream conditions to satisfy the difference equations identically. Results of calculations illustrating the importance of these effects will be presented.

Caughey, D. A.↗

NASA NTRS1979

Formulation of coupled rotor/fuselage equations of motion

The governing equations of motion of a helicopter rotor coupled to a rigid body fuselage are derived. A consistent formulation is used to derive nonlinear periodic coefficient equations of motion which can be used to study coupled rotor/fuselage dynamics in forward flight. The methodology of rotor/fuselage coupling is clearly described and the importance of an ordering scheme in deriving consistent nonlinear equations of motion is reviewed. The final equations which are presented in partial differential form can be used to model coupled rotor/fuselage aeroelastic response or stability problems.

Warmbrodt, W.↗

NASA NTRS1981

Influence of boundary approximations and conditions on finite difference solutions

Numerical representations of boundary approximations and conditions for three problems are investigated to determine the resulting global accuracy of the steady state solution. Numerical accuracy with various boundary approximations is determined for quasi-one-dimensional inviscid flow in a duct with the interior grid points evaluated using the MacCormack scheme. When an extrapolation approximation with first order local truncation error is used, the global second order accuracy of the difference scheme can be destroyed. For one dimensional flow in a porous medium, an implicit midpoint difference scheme which is consistent with the boundary conditions is developed without the need of boundary approximations. A dissipative model problem is solved with the boundary conditions discretized with first and second order accuracy. The overall second order accuracy of the difference scheme is destroyed if first order numerical representation of one of the boundary conditions is used. With a boundary approximation, the second order global accuracy of the model problem is retained if either second order extrapolation or first order representation of the governing equation is used.

Blottner, F. G.↗

NASA NTRS1988

Consistent boundary conditions for reduced Navier-Stokes (RNS) scheme applied to three-dimensional internal viscous flows

Three-dimensional internal viscous flow problems are presently addressed by a consistent and efficient set of boundary conditions for the multisweep space-marching, pressure-elliptic Reduced Navier-Stokes scheme. Since continuity is by this means directly satisfied at all points in the flow domain, the first-order momentum equations are directly solvable for pressure without the requirement for a Poisson pressure-correction equation. Incompressible flow solutions are obtained for straight and curved ducts of square cross section, in order to validate the procedure. Usefulness is demonstrated for internal flows with strong interactions, as would be found in turbomachine geometries.

Reddy, D. R.↗

NASA NTRS1984

Annihilation of angular momentum drift during spinning-up and thrusting maneuvers of rigid bodies

A very simple, yet accurate, heuristic solution for the spiral path of the angular momentum vector during spin-up and spin-down maneuvers of rigid body spacecraft is presented. A two-burn scheme is proposed consisting of a burn, a coast, and a second burn of the spin thruster. The appropriate burn times are found by a transcendental equation similar to Kepler's equation. Numerical results verify the accuracy of the burn, coast, and burn times given by closed form expressions. The scheme can also be applied to the problem of axial thrusting during constant spin.

Longuski, J. M.↗

NASA NTRS1995

Self consistent MHD modeling of the solar wind from coronal holes with distinct geometries

Utilizing an iterative scheme, a self-consistent axisymmetric MHD model for the solar wind has been developed. We use this model to evaluate the properties of the solar wind issuing from the open polar coronal hole regions of the Sun, during solar minimum. We explore the variation of solar wind parameters across the extent of the hole and we investigate how these variations are affected by the geometry of the hole and the strength of the field at the coronal base.

Stewart, G. A.↗

NASA NTRS2004

Adaptive Numerical Dissipative Control in High Order Schemes for Multi-D Non-Ideal MHD

The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free of numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multi-resolution wavelets (WAV) (for the above types of flow feature). These filter approaches also provide a natural and efficient way for the minimization of Div(B) numerical error. The filter scheme consists of spatially sixth order or higher non-dissipative spatial difference operators as the base scheme for the inviscid flux derivatives. If necessary, a small amount of high order linear dissipation is used to remove spurious high frequency oscillations. For example, an eighth-order centered linear dissipation (AD8) might be included in conjunction with a spatially sixth-order base scheme. The inviscid difference operator is applied twice for the viscous flux derivatives. After the completion of a full time step of the base scheme step, the solution is adaptively filtered by the product of a 'flow detector' and the 'nonlinear dissipative portion' of a high-resolution shock-capturing scheme. In addition, the scheme independent wavelet flow detector can be used in conjunction with spatially compact, spectral or spectral element type of base schemes. The ACM and wavelet filter schemes using the dissipative portion of a second-order shock-capturing scheme with sixth-order spatial central base scheme for both the inviscid and viscous MHD flux derivatives and a fourth-order Runge-Kutta method are denoted.

Yee, H. C.↗

NASA NTRS1979

A viscous flow analysis for the tip vortex generation process

A three dimensional, forward-marching, viscous flow analysis is applied to the tip vortex generation problem. The equations include a streamwise momentum equation, a streamwise vorticity equation, a continuity equation, and a secondary flow stream function equation. The numerical method used combines a consistently split linearized scheme for parabolic equations with a scalar iterative ADI scheme for elliptic equations. The analysis is used to identify the source of the tip vortex generation process, as well as to obtain detailed flow results for a rectangular planform wing immersed in a high Reynolds number free stream at 6 degree incidence.

Shamroth, S. J.↗

NASA NTRS1980

Nonlinear Aeroelastic Equations of Motion of Twisted, Nonuniform, Flexible Horizontal-Axis Wind Turbine Blades

The second-degree nonlinear equations of motion for a flexible, twisted, nonuniform, horizontal axis wind turbine blade were developed using Hamilton's principle. A mathematical ordering scheme which was consistent with the assumption of a slender beam was used to discard some higher-order elastic and inertial terms in the second-degree nonlinear equations. The blade aerodynamic loading which was employed accounted for both wind shear and tower shadow and was obtained from strip theory based on a quasi-steady approximation of two-dimensional, incompressible, unsteady, airfoil theory. The resulting equations had periodic coefficients and were suitable for determining the aeroelastic stability and response of large horizontal-axis wind turbine blades.

Kaza, K. R. V.↗

NASA NTRS1980

Grid generation for two-dimensional finite element flowfield computation

The finite element method for fluid dynamics was used to develop a two dimensional mesh generation scheme. The method consists of shearing and conformal maps with upper and lower surfaces handled independently to allow sharp leading edges. The method also generates meshes of triangular or quadrilateral elements.

Tatum, K. E.↗

NASA NTRS1991

Homogeneity of spatial correlation statistics of tropical oceanic rainfall

The possibility of uniform horizontal correlation scales for tropical oceanic rainfall has been examined by a study of satellite-observed microwave data as a proxy measure of rain rates. From the brightness temperatures from the electrically scanning microwave radiometer on Nimbus 5 near nadir during the year 1974, the mean spatial autocorrelation function as a function of simultaneous pixel separation was calculated in each 5 by 5 deg grid box over the tropical Pacific and Atlantic for each season. The equal-time spatially lagged correlations were compared for geographical dependence to investigate the hypothesis of homogeneity. A simple model of the spatial statistics of the microwave brightness temperatures was used, consisting of a mixture of uncorrelated spatial white noise incoherently superimposed on a spatially correlated field (spatial red noise). The red noise signals are presumed to be generated by convective activity in the tropical atmosphere. The parameters of the red noise are consistent with this scheme over the tropical oceans, yielding a uniform spatial scale of about 50 km throughout.

Shin, Kyung-Sup↗