Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical schemes”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Aerothermal modeling program, phase 2

The main objective of the NASA sponsored Aerothermal Modeling Program, Phase 2--Element A, is to develop an improved numerical scheme for predicting combustor flow fields. This effort consists of the following three technical tasks. Task 1 involves the selection and evaluation of various candidate numerical techniques. Task 2 involves an in-depth evaluation of the selected numerical schemes. Task 3 involves the convection-diffusion scheme and the direct solver that will be incorporated in the NASA 3-D elliptic code (COM3S).

Karki, K. C.↗
NASA NTRSDate not supplied

Numerical Investigation of Combustor-Turbine Interactions with a Two-Stage High-Pressure Turbine Including Cooling Airflows

Combustor–turbine interactions (CTIs) are investigated by performing three-dimensional unsteady simulations using a realistic combustor and high-pressure turbine (HPT) configurations from the Energy Efficient Engine (E3) program. To understand CTIs, we compare the numerically predicted flow fields from single-component simulations (Step 1: the combustor + the first-stage stator of turbine; Step 2: the two-stage HPT imposing the time-averaged flow solution from Step 1 as the inflow boundary condition) and a fully coupled combustor–turbine simulation (Step 3) at the sea-level takeoff (SLTO) condition. In addition to three previous simulations where the cooling airflows inside the HPT had been neglected for all Step 1, Step 2, and Step 3, two new simulations of Step 2 and Step 3 take into account the cooling airflows using the source team approach. In this approach, to mimic the cooling airflows, we impose the source term at a specific area of a cooling airflow hole at the solid surfaces without making a mesh of each hole. The objective of this study is twofold. One is to investigate the effect of the presence of the cooling airflows on the aerodynamics of the combustor and HPT as well as the HPT performance. The second is to perform a detailed comparison among the calculated flowfields by two different numerical schemes, the central-difference with the standard Jameson–Schmidt–Turkel (CD-JST) scheme and the AUSM scheme. There is a noticeable difference in the hot-streak distributions at the first-stage stator. In addition, depending on the choice of the numerical scheme and the presence of the cooling airflows, an occurrence of shock waves appearing at the first-stage stator is greatly influenced. Thus, this has a noticeable impact on the HPT performance. It is shown that the estimated turbine efficiencies from Step 3 are about 7 % less than the ones from Step 2.

Combustor-turbine interaction↗
NASA NTRS2008

Problems Associated with Grid Convergence of Functionals

The current use of functionals to evaluate order-of-convergence of a numerical scheme can lead to incorrect values. The problem comes about because of interplay between the errors from the evaluation of the functional, e.g., quadrature error, and from the numerical scheme discretization. Alternative procedures for deducing the order-property of a scheme are presented. The problem is studied within the context of the inviscid supersonic flow over a blunt body; however, the problem and solutions presented are not unique to this example.

Salas, Manuel D.↗
NASA NTRS2009

On Problems Associated with Grid Convergence of Functionals

The current use of functionals to evaluate order-of-convergence of a numerical scheme can lead to incorrect values. The problem comes about because of interplay between the errors from the evaluation of the functional, e.g., quadrature error, and from the numerical scheme discretization. Alternative procedures for deducing the order property of a scheme are presented. The problems are studied within the context of the inviscid supersonic flow over a blunt body; however, the problems and solutions presented are not unique to this example.

Salas, Manuael D.↗
NASA NTRS1998

Higher-Order Compact Schemes for Numerical Simulation of Incompressible Flows

A higher order accurate numerical procedure has been developed for solving incompressible Navier-Stokes equations for 2D or 3D fluid flow problems. It is based on low-storage Runge-Kutta schemes for temporal discretization and fourth and sixth order compact finite-difference schemes for spatial discretization. The particular difficulty of satisfying the divergence-free velocity field required in incompressible fluid flow is resolved by solving a Poisson equation for pressure. It is demonstrated that for consistent global accuracy, it is necessary to employ the same order of accuracy in the discretization of the Poisson equation. Special care is also required to achieve the formal temporal accuracy of the Runge-Kutta schemes. The accuracy of the present procedure is demonstrated by application to several pertinent benchmark problems.

Wilson, Robert V.↗
NASA NTRS2016

Computational Flow Field in Energy Efficient Engine (EEE)

In this paper, preliminary results for the recently-updated Open National Combustion Code (Open NCC) as applied to the EEE are presented. The comparison between two different numerical schemes, the standard Jameson-Schmidt-Turkel (JST) scheme and the advection upstream splitting method (AUSM), is performed for the cold flow and the reacting flow calculations using the RANS. In the cold flow calculation, the AUSM scheme predicts a much stronger reverse flow in the central recirculation zone. In the reacting flow calculation, we test two cases: gaseous fuel injection and liquid spray injection. In the gaseous fuel injection case, the overall flame structures of the two schemes are similar to one another, in the sense that the flame is attached to the main nozzle, but is detached from the pilot nozzle. However, in the exit temperature profile, the AUSM scheme shows a more uniform profile than that of the JST scheme, which is close to the experimental data. In the liquid spray injection case, we expect different flame structures in this scenario. We will give a brief discussion on how two numerical schemes predict the flame structures inside the EEE using different ways to introduce the fuel injection.

CFD↗
NASA NTRS1993

Numerical Simulation of a High Mach Number Jet Flow

The recent efforts to develop accurate numerical schemes for transition and turbulent flows are motivated, among other factors, by the need for accurate prediction of flow noise. The success of developing high speed civil transport plane (HSCT) is contingent upon our understanding and suppression of the jet exhaust noise. The radiated sound can be directly obtained by solving the full (time-dependent) compressible Navier-Stokes equations. However, this requires computational storage that is beyond currently available machines. This difficulty can be overcome by limiting the solution domain to the near field where the jet is nonlinear and then use acoustic analogy (e.g., Lighthill) to relate the far-field noise to the near-field sources. The later requires obtaining the time-dependent flow field. The other difficulty in aeroacoustics computations is that at high Reynolds numbers the turbulent flow has a large range of scales. Direct numerical simulations (DNS) cannot obtain all the scales of motion at high Reynolds number of technological interest. However, it is believed that the large scale structure is more efficient than the small-scale structure in radiating noise. Thus, one can model the small scales and calculate the acoustically active scales. The large scale structure in the noise-producing initial region of the jet can be viewed as a wavelike nature, the net radiated sound is the net cancellation after integration over space. As such, aeroacoustics computations are highly sensitive to errors in computing the sound sources. It is therefore essential to use a high-order numerical scheme to predict the flow field. The present paper presents the first step in a ongoing effort to predict jet noise. The emphasis here is in accurate prediction of the unsteady flow field. We solve the full time-dependent Navier-Stokes equations by a high order finite difference method. Time accurate spatial simulations of both plane and axisymmetric jet are presented. Jet Mach numbers of 1.5 and 2.1 are considered. Reynolds number in the simulations was about a million. Our numerical model is based on the 2-4 scheme by Gottlieb & Turkel. Bayliss et al. applied the 2-4 scheme in boundary layer computations. This scheme was also used by Ragab and Sheen to study the nonlinear development of supersonic instability waves in a mixing layer. In this study, we present two dimensional direct simulation results for both plane and axisymmetric jets. These results are compared with linear theory predictions. These computations were made for near nozzle exit region and velocity in spanwise/azimuthal direction was assumed to be zero.

Hayder, M. Ehtesham↗
NASA NTRS1990

Counterrotating prop-fan simulations which feature a relative-motion multiblock grid decomposition enabling arbitrary time-steps

Improvements are presented of a computer algorithm developed for the time-accurate flow analysis of rotating machines. The flow model is a finite volume method utilizing a high-resolution approximate Riemann solver for interface flux definitions. The numerical scheme is a block LU implicit iterative-refinement method which possesses apparent unconditional stability. Multiblock composite gridding is used to orderly partition the field into a specified arrangement of blocks exhibiting varying degrees of similarity. Block-block relative motion is achieved using local grid distortion to reduce grid skewness and accommodate arbitrary time step selection. A general high-order numerical scheme is applied to satisfy the geometric conservation law. An even-blade-count counterrotating unducted fan configuration is chosen for a computational study comparing solutions resulting from altering parameters such as time step size and iteration count. The solutions are compared with measured data.

Janus, J. Mark↗
NASA NTRS1991

An investigation of voids formation mechanisms and their effects on freeze and thaw processes of lithium and lithium fluoride

The mechanisms of void formation during the cooldown and freezing of lithium coolant within the primary loop of SP-100 type systems are investigated. These mechanisms are: (1) homogeneous nucleation; (2) heterogeneous nucleation; (3) normal segregation of helium gas dissolved in liquid lithium; and (4) shrinkage of lithium during freezing. To evaluate the void formation potential due to segregation, a numerical scheme that couples the freezing and mass diffusion processes in both the solid and liquid regions is developed. The results indicated that the formation of He bubbles is unlikely by either homogeneous or heterogeneous nucleation during the cooldown process. However, homogeneous nucleation of He bubbles following the segregation of dissolved He in liquid lithium ahead of the solid-liquid interface is likely to occur. Results also show that total volume of He void is insignificant when compared to that of shrinkage voids. In viewing this, the subsequent research focuses on the effects of shrinkage void forming during freezing of lithium on subsequent thaw processes are investigated using a numerical scheme that is based on a single (solid/liquid) cell approach. The cases of lithium-fluoride are also investigated to show the effect of larger volume shrinkage upon freezing on the freeze and thaw processes. Results show that a void forming at the wall appreciably reduces the solid-liquid interface velocity, during both freeze and thaw, and causes a substantial rise in the wall temperature during thaw. However, in the case of Li, the maximum wall temperature was much lower than the melting temperature of PWC-11, which is used as the structure material in the SP-100 system. Hence, it is included that a formation of hot spots is unlikely during the startup or restart of the SP-100 system.

El-Genk, Mohamed S.↗
NASA NTRS1993

Numerical vorticity capturing for vortex-solid body interaction problems

A new numerical method based on adding a term to the Euler/Navier-Stokes equations has been demonstrated it can effectively treat vortex-dominated flows using low-order numerical schemes and coarse grids. The numerical diffusion introduced by convection schemes can be eliminated by modifying the velocity field as a result of solving the modified flow governing equations. The modification of the velocity field, which conserves the total vorticity, essentially convects the vorticity toward its local extreme to offset the numerical diffusion. This method preserves the vortex structure even when vortices travel on coarse grid region. The method is implemented with two distinctive codes: a Navier-Stokes flow solver based on vorticity-velocity formulation and an Euler code based on velocity-pressure formulation. Problems presented in this paper include airfoil dynamic stall, vortex airfoil interaction, and vortex-fuselage interactions.

Wang, Clin M.↗
NASA NTRS2001

The CE/SE Method: a CFD Framework for the Challenges of the New Millennium

The space-time conservation element and solution element (CE/SE) method, which was originated and is continuously being developed at NASA Glenn Research Center, is a high-resolution, genuinely multidimensional and unstructured-mesh compatible numerical method for solving conservation laws. Since its inception in 1991, the CE/SE method has been used to obtain highly accurate numerical solutions for 1D, 2D and 3D flow problems involving shocks, contact discontinuities, acoustic waves, vortices, shock/acoustic waves/vortices interactions, shock/boundary layers interactions and chemical reactions. Without the aid of preconditioning or other special techniques, it has been applied to both steady and unsteady flows with speeds ranging from Mach number = 0.00288 to 10. In addition, the method has unique features that allow for (i) the use of very simple non-reflecting boundary conditions, and (ii) a unified wall boundary treatment for viscous and inviscid flows. The CE/SE method was developed with the conviction that, with a solid foundation in physics, a robust, coherent and accurate numerical framework can be built without involving overly complex mathematics. As a result, the method was constructed using a set of design principles that facilitate simplicity, robustness and accuracy. The most important among them are: (i) enforcing both local and global flux conservation in space and time, with flux evaluation at an interface being an integral part of the solution procedure and requiring no interpolation or extrapolation; (ii) unifying space and time and treating them as a single entity; and (iii) requiring that a numerical scheme be built from a nondissipative core scheme such that the numerical dissipation can be effectively controlled and, as a result, will not overwhelm the physical dissipation. Part I of the workshop will be devoted to a discussion of these principles along with a description of how the ID, 2D and 3D CE/SE schemes are constructed. In Part II, various applications of the CE/SE method, particularly those involving chemical reactions and acoustics, will be presented. The workshop will be concluded with a sketch of the future research directions.

Chang, Sin-Chung↗
NASA NTRS2003

A Model for the Oxidation of C/SiC Composite Structures

A mathematical theory and an accompanying numerical scheme have been developed for predicting the oxidation behavior of C/SiC composite structures. The theory is derived from the mechanics of the flow of ideal gases through a porous solid. Within the mathematical formulation, two diffusion mechanisms are possible: (1) the relative diffusion of one species with respect to the mixture, which is concentration gradient driven and (2) the diffusion associated with the average velocity of the gas mixture, which is total gas pressure gradient driven. The result of the theoretical formulation is a set of two coupled nonlinear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations must be solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of space and time. The local rate of carbon oxidation is determined as a function of space and time using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The nonlinear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual method, allowing for the solution of the differential equations numerically. The end result is a numerical scheme capable of determining the variation of the local carbon oxidation rates as a function of space and time for any arbitrary C/SiC composite structures.

Sullivan, Roy M.↗
NASA NTRS1985

Lubrication of nonconformal contacts

Minimum film thickness results for piezoviscous-rigid regime of lubrication are developed for a compressible Newtonian fluid with Roelands viscosity. The results provide a basis for the analysis and design of a wide range of machine elements operating in the piezoviscous-rigid regime of lubrication. A new numerical method of calculating elastic deformation in contact stresses is developed using a biquadratic polynomial to approximate the pressure distribution on the whole domain analyzed. The deformation of every node is expressed as a linear combination of the nodal pressures whose coefficients can be combined into an influence coefficient matrix. This approach has the advantages of improved numerical accuracy, less computing time and smaller storage size required for influence matrix. The ideal elastohydrodynamic lubrication is extended to real bearing systems in order to gain an understanding of failure mechanisms in machine elements. The improved elastic deformation calculation is successfully incorporated into the EHL numerical scheme. Using this revised numerical technique and the flow factor model developed by Patir and Cheng (1978) the surface roughness effects on the elastohydrodynamic lubrication of point contact is considered. Conditions typical of an EHL contact in the piezoviscous-elastic regime entrained in pure rolling are investigated. Results are compared with the smooth surface solutions. Experiments are conducted to study the transient EHL effects in instrument ball bearings.

Jeng, Y. R.↗
NASA NTRS1986

Theory of electrophoretic separations. II - Construction of a numerical simulation scheme and its applications

The implementation of a previously derived model describing electrophoretic transport processes is reformulated, and its usefulness is demonstrated by simulating moving boundary electrophoresis and isoelectric focusing with immobilized species. The resulting solution is stable, evolves smoothly when suitable small step sizes are used, and has a loosely defined convergence as the step sizes are made smaller. The major computational problem involves a need to use small spatial step sizes at high current levels.

Palusinski, O. A.↗
NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗
NASA NTRS1986

The discrete one-sided Lipschitz condition for convex scalar conservation laws

Physical solutions to convex scalar conservation laws satisfy a one-sided Lipschitz condition (OSLC) that enforces both the entropy condition and their variation boundedness. Consistency with this condition is therefore desirable for a numerical scheme and was proved for both the Godunov and the Lax-Friedrichs scheme--also, in a weakened version, for the Roe scheme, all of them being only first order accurate. A new, fully second order scheme is introduced here, which is consistent with the OSLC. The modified equation is considered and shows interesting features. Another second order scheme is then considered and numerical results are discussed.

Brenier, Yann↗
NASA NTRS1990

Parallelization of implicit finite difference schemes in computational fluid dynamics

Implicit finite difference schemes are often the preferred numerical schemes in computational fluid dynamics, requiring less stringent stability bounds than the explicit schemes. Each iteration in an implicit scheme involves global data dependencies in the form of second and higher order recurrences. Efficient parallel implementations of such iterative methods are considerably more difficult and non-intuitive. The parallelization of the implicit schemes that are used for solving the Euler and the thin layer Navier-Stokes equations and that require inversions of large linear systems in the form of block tri-diagonal and/or block penta-diagonal matrices is discussed. Three-dimensional cases are emphasized and schemes that minimize the total execution time are presented. Partitioning and scheduling schemes for alleviating the effects of the global data dependencies are described. An analysis of the communication and the computation aspects of these methods is presented. The effect of the boundary conditions on the parallel schemes is also discussed.

Decker, Naomi H.↗
NASA NTRS2002

High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations

We present high-order semi-discrete central-upwind numerical schemes for approximating solutions of multi-dimensional Hamilton-Jacobi (HJ) equations. This scheme is based on the use of fifth-order central interpolants like those developed in [1], in fluxes presented in [3]. These interpolants use the weighted essentially nonoscillatory (WENO) approach to avoid spurious oscillations near singularities, and become "central-upwind" in the semi-discrete limit. This scheme provides numerical approximations whose error is as much as an order of magnitude smaller than those in previous WENO-based fifth-order methods [2, 1]. Thee results are discussed via examples in one, two and three dimensions. We also pregnant explicit N-dimensional formulas for the fluxes, discuss their monotonicity and tl!e connection between this method and that in [2].

Bryson, Steve↗