Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite difference 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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Design sensitivity analysis for an aeroelastic optimization of a helicopter blade

The sensitivity of vibratory hub loads of a four-bladed hingeless rotor with respect to blade design parameters is investigated using a finite element formulation in space and time. Design parameters include nonstructural mass distribution (spanwise and chordwise), chordwise offset of center of gravity from aerodynamic center, blade bending stiffnesses (flap, lag and torsion). Hub loads selected are 4/rev vertical hub shear and 3/rev hub moment in the rotating reference frame. The sensitivity derivatives of vertical hub loads with respect to blade design parameters are compared using two approaches, finite difference scheme and analytical approach using chain rule differentiation. The analytical derivative approach developed as an integral part of response solution (finite element in time) is a powerful method for an aeroelastic optimization of a helicopter rotor.

Lim, Joon↗
NASA NTRS1990

Numerical methods for hypersonic boundary layer stability

Four different schemes for solving compressible boundary layer stability equations are developed and compared, considering both the temporal and spatial stability for a global eigenvalue spectrum and a local eigenvalue search. The discretizations considered encompass: (1) a second-order-staggered finite-difference scheme; (2) a fourth-order accurate, two-point compact scheme; (3) a single-domain Chebychev spectral collocation scheme; and (4) a multidomain spectral collocation scheme. As Mach number increases, the performance of the single-domain collocation scheme deteriorates due to the outward movement of the critical layer; a multidomain spectral method is accordingly designed to furnish superior resolution of the critical layer.

Malik, M. R.↗
NASA NTRS1985

Shock capturing

Recent developments which have improved the understanding of how finite difference methods resolve discontinuous solutions to hyperbolic partial differential equations are discussed. As a result of this understanding improved shock capturing methods are currently being developed and tested. Some of these methods are described and numerical results are presented showing their performance on problems containing shocks in one and two dimensions. A conservative difference scheme is defined. Conservation implies that, except in very special circumstances, shocks must be spread over at least two grid intervals. These two interval shocks are actually attained in one dimension if the shock is steady and an upwind scheme is used. By analyzing this case, the reason for this excellent shock resolution can be determined. This result is used to provide a mechanism for improving the resolution of two dimensional steady shocks. Unfortunately, this same analysis shows that these results cannot be extended to shocks which move relative to the computing grid. Total variation diminishing (TVD) finite difference schemes and flux limiters are introduced to deal with money shocks and contact discontinuities.

Davis, S. F.↗
NASA NTRS1990

Energy conserving and potential-enstrophy dissipating schemes for the shallow water equations

To incorporate potential enstrophy dissipation into discrete shallow water equations with no or arbitrarily small energy dissipation, a family of finite-difference schemes have been derived with which potential enstrophy is guaranteed to decrease while energy is conserved (when the mass flux is nondivergent and time is continuous). Among this family of schemes, there is a member that minimizes the spurious impact of infinite potential vorticities associated with infinitesimal fluid depth. The scheme is, therefore, useful for problems in which the free surface may intersect with the lower boundary.

Arakawa, Akio↗
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↗
NASA NTRS2004

Nonlinear filtering and limiting in high order methods for ideal and non-ideal MHD

The various filtering mechanisms and base scheme options of the newly developed adaptive numerical dissipation control in spatially high order filter schemes for the ideal and non-ideal magnetohydrodynamics (MHD) equations are investigated. These filter schemes are applicable to complex unsteady MHD high-speed shock/shear/turbulence problems. They also provide a natural and efficient way for the minimization of Div(B) numerical error. The type of spatial base scheme to be used in conjunction with our filter idea is very general. For example, spectral, compact and non-compact spatially central finite difference schemes are possible candidates. The adaptive numerical dissipation mechanism 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 to leave the rest of the region free from numerical dissipation contamination. The numerical dissipation considered consists of high order linear dissipation for the suppression of high frequency oscillation and the nonlinear dissipative portion of high-resolution shock-capturing methods for discontinuity capturing. The applicable nonlinear dissipative portion of high-resolution shock-capturing methods is also very general. The objective of this paper is to investigate the performance of using compact and non-compact central base schemes in conjunction with three commonly used types of nonlinear numerical dissipation for both the ideal and non-ideal MHD. This extended abstract shows the performance of three nonlinear filters in conjunction with a sixth-order non-compact spatial central base scheme. In the final paper, the high order compact spatial central base scheme will be illustrated and compared with the non-compact base scheme. The reason for the investigation of the high order compact spatial central base scheme over the non-compact base scheme is to evaluate if additional accuracy can be gained in regions of fine scale turbulence that are away from shocks/shears.

Yee,H. C.↗
NASA NTRS1988

Numerical study of the vortex burst phenomenon for delta wings

A flux-difference splitting scheme is employed to compute low-speed flows over a delta wing for angles of attack from 0 to 40 deg as steady-state solutions to the three-dimensional, Reynolds-averaged Navier-Stokes equations in their thin-layer approximation. The finite-difference scheme is made spatially second-order accurate by applying a total variation diminishing-like discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. Using first-order accurate Euler backward-time differencing, an efficient implicit algorithm is constructed, which combines approximate factorization in cross planes with a symmetric planar Gauss-Seidel relaxation in the remaining third spatial direction. The geometry of the thin (maximum thickness is 0.021), slender (aspect ratio is unity), sharp-edged delta wing is taken from Hummel's (1967, 1978) wind tunnel model. Over the entire angle-of-attack range, the computed values of lift and pitching moment are in good agreement with the experimental data. Also details of the flow-fieldlike spanwise surface pressure distributions compare well with the experiment. Computed flow-field results with a bubble-type vortex burst are analyzed in detail.

Hartwich, PETER-M.↗
NASA NTRS1996

Inhomogeneous Radiation Boundary Conditions Simulating Incoming Acoustic Waves for Computational Aeroacoustics

A set of nonhomogeneous radiation and outflow conditions which automatically generate prescribed incoming acoustic or vorticity waves and, at the same time, are transparent to outgoing sound waves produced internally in a finite computation domain is proposed. This type of boundary condition is needed for the numerical solution of many exterior aeroacoustics problems. In computational aeroacoustics, the computation scheme must be as nondispersive ans nondissipative as possible. It must also support waves with wave speeds which are nearly the same as those of the original linearized Euler equations. To meet these requirements, a high-order/large-stencil scheme is necessary The proposed nonhomogeneous radiation and outflow boundary conditions are designed primarily for use in conjunction with such high-order/large-stencil finite difference schemes.

Tam, Christopher K. W.↗
NASA NTRS1981

A split explicit reformulation of the regional numerical weather prediction model of the Japan Meteorological Agency

The split explicit integration scheme for numerical weather prediction models is employed in a version of the regional numerical weather prediction model of the Japan Meteorological Agency. The finite-difference scheme of the model is designed in the manner proposed by Okamura (1975). The horizontal advection terms in the governing equations are integrated with a time step limited by the wind speed while the terms which describe inertial-gravity oscillations are integrated in a succession of shorter time steps. The physical processes included within the model are precipitation, small-scale convection, surface exchanges of sensible and latent heat, and radiative heating and cooling. An example of a surface pressure forecast over Europe is shown for initial data observed at 0000 GMT 29 December 1979. Quantitative precipitation forecasts over Europe and North America for the 24 h period beginning at 0000 GMT 30 December 1979 are also shown. It is concluded that the model is capable of realistically depicting the evolution of synoptic-scale systems.

Duffy, D. G.↗
NASA NTRS1993

Studies on nonequilibrium phenomena in supersonic chemically reacting flows

This study deals with a systematic investigation of nonequilibrium processes in supersonic combustion. The two-dimensional, elliptic Navier-Stokes equations are used to investigate supersonic flows with nonequilibrium chemistry and thermodynamics, coupled with radiation, for hydrogen-air systems. The explicit, unsplit MacCormack finite-difference scheme is used to advance the governing equations in time, until convergence is achieved. For a basic understanding of the flow physics, premixed flows undergoing finite rate chemical reactions are investigated. Results obtained for specific conditions indicate that the radiative interactions vary substantially, depending on reactions involving HO2 and NO species, and that this can have a noticeable influence on the flowfield. The second part of this study deals with premixed reacting flows under thermal nonequilibrium conditions. Here, the critical problem is coupling of the vibrational relaxation process with the radiative heat transfer. The specific problem considered is a premixed expanding flow in a supersonic nozzle. Results indicate the presence of nonequilibrium conditions in the expansion region of the nozzle. This results in reduction of the radiative interactions in the flowfield. Next, the present study focuses on investigation of non-premixed flows under chemical nonequilibrium conditions. In this case, the main problem is the coupled turbulence-chemistry interaction. The resulting formulation is validated by comparison with experimental data on reacting supersonic coflowing jets. Results indicate that the effect of heat release is to lower the turbulent shear stress and the mean density. The last part of this study proposes a new theoretical formulation for the coupled turbulence-radiation interactions. Results obtained for the coflowing jets experiment indicate that the effect of turbulence is to enhance the radiative interactions.

Tiwari, S. N.↗
NASA NTRS1993

A new sensitivity analysis for structural optimization of composite rotor blades

This paper presents a detailed mathematical derivation of the sensitivity derivatives for the structural dynamic, aeroelastic stability and response characteristics of a rotor blade in hover and forward flight. The formulation is denoted by the term semianalytical approach, because certain derivatives have to be evaluated by a finite difference scheme. Using the present formulation, sensitivity derivatives for the structural dynamic and aeroelastic stability characteristics, were evaluated for both isotropic and composite rotor blades. Based on the results, useful conclusions are obtained regarding the relative merits of the semi-analytical approach, for calculating sensitivity derivatives, when compared to a pure finite difference approach.

Venkatesan, C.↗
NASA NTRS1984

Error reduction program: A progress report

Five finite differences schemes were evaluated for minimum numerical diffusion in an effort to identify and incorporate the best error reduction scheme into a 3D combustor performance code. Based on this evaluated, two finite volume method schemes were selected for further study. Both the quadratic upstream differencing scheme (QUDS) and the bounded skew upstream differencing scheme two (BSUDS2) were coded into a two dimensional computer code and their accuracy and stability determined by running several test cases. It was found that BSUDS2 was more stable than QUDS. It was also found that the accuracy of both schemes is dependent on the angle that the streamline make with the mesh with QUDS being more accurate at smaller angles and BSUDS2 more accurate at larger angles. The BSUDS2 scheme was selected for extension into three dimensions.

Syed, S. A.↗
NASA NTRS1985

Higher order methods for convection-diffusion problems

This paper applies C1 cubic Hermite polynomials embedded in an orthogonal collocation scheme to the spatial discretization of the unsteady nonlinear Burgers equation as a model of the equations of fluid mechanics. The temporal discretization is carried out by means of either a noniterative finite difference or an iterative finite difference procedure. Results of this method are compared with those of a second-order finite difference scheme and a splined-cubic Taylor's series scheme. Stability limits are derived and the matrix structure of the several schemes are compared.

Murphy, J. D.↗
NASA NTRS1986

Unsteady hybrid vortex technique for transonic vortex flows and flutter applications

The well developed surface panel method is extended by adding a volume integral term allowing calculation of the full effect of compressibility. The full effect of compressibility is calculated by using mixed type finite difference scheme consistent with the mixed nature of transonic flow. The solution is obtained through successive iteration cycles for subcritical flows and for critical flows the solution is obtained through satisfying the Rankine-Hugoniot relations across the captured shock in addition to the captured shock. The method is applied to airfoils in subcritical and critical flows and the results are in good agreement with the experimental data and finite-difference solutions of the full-potential and Euler equations.

Kandil, Osama A.↗
DOE OSTI2024

2-D seismic wave propagation using the distributional finite-difference method: further developments and potential for global seismology

SUMMARY We present a time-domain distributional finite-difference scheme based on the Lebedev staggered grid for the numerical simulation of wave propagation in acoustic and elastic media. The central aspect of the proposed method is the representation of the stresses and displacements with different sets of B-splines functions organized according to the staggered grid. The distributional finite-difference approach allows domain-decomposition, heterogeneity of the medium, curvilinear mesh, anisotropy, non-conformal interfaces, discontinuous grid and fluid–solid interfaces. Numerical examples show that the proposed scheme is suitable to model wave propagation through the Earth, where sharp interfaces separate large, relatively homogeneous layers. A few domains or elements are sufficient to represent the Earth’s internal structure without relying on advanced meshing techniques. We compare seismograms obtained with the proposed scheme and the spectral element method, and we show that our approach offers superior accuracy, reduced memory usage, and comparable efficiency.

Geochemistry & Geophysics↗
NASA NTRS1991

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models

Turbulent combustion can not be simulated adequately by conventional moment closure turbulence models. The difficulty lies in the fact that the reaction rate is in general an exponential function of the temperature, and the higher order correlations in the conventional moment closure models of the chemical source term can not be neglected, making the applications of such models impractical. The probability density function (pdf) method offers an attractive alternative: in a pdf model, the chemical source terms are closed and do not require additional models. A grid dependent Monte Carlo scheme was studied, since it is a logical alternative, wherein the number of computer operations increases only linearly with the increase of number of independent variables, as compared to the exponential increase in a conventional finite difference scheme. A new algorithm was devised that satisfies a restriction in the case of pure diffusion or uniform flow problems. Although for nonuniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

Hsu, Andrew T.↗
NASA NTRS1996

A Mixed Finite Volume Element Method for Flow Calculations in Porous Media

A key ingredient in the simulation of flow in porous media is the accurate determination of the velocities that drive the flow. The large scale irregularities of the geology, such as faults, fractures, and layers suggest the use of irregular grids in the simulation. Work has been done in applying the finite volume element (FVE) methodology as developed by McCormick in conjunction with mixed methods which were developed by Raviart and Thomas. The resulting mixed finite volume element discretization scheme has the potential to generate more accurate solutions than standard approaches. The focus of this paper is on a multilevel algorithm for solving the discrete mixed FVE equations. The algorithm uses a standard cell centered finite difference scheme as the 'coarse' level and the more accurate mixed FVE scheme as the 'fine' level. The algorithm appears to have potential as a fast solver for large size simulations of flow in porous media.

Jim E Jones↗
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.↗