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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Effects of spatial order of accuracy on the computation of vortical flowfields

The effect of the order-of-accuracy, used for the spatial discretization, on the resolution of the leading edge vortices over sharp-edged delta wings is investigated. The flowfield is computed using a viscous/inviscid zonal approach. The viscous flow in the vicinity of the wing is computed using the conservative formulation of the compressible, thin-layer Navier-Stokes equations. The leeward-side vortical flowfield and the other flow regions away from the surface are computed as inviscid. The time integration is performed with both an explicit fourth-order Runge-Kutta scheme and an implicit, factorized, iterative scheme. High-order-accurate inviscid fluxes are computed using both a conservative and a non-conservative (primitive variable) formulation. The nonlinear, inviscid terms of the primitive variable form of the governing equations are evaluated with a finite-difference numerical scheme based on the sign of the eigenvalues. High-order, upwind-biased, finite difference formulas are used to evaluate the derivatives of the nonlinear convective terms. Computed results are compared with available experimental data, and comparisons of the flowfield in the vicinity of the vortex cores are presented.

Ekaterinaris, J. A.↗

NASA NTRS2013

Numerical Modeling of Ablation Heat Transfer

A unique numerical method has been developed for solving one-dimensional ablation heat transfer problems. This paper provides a comprehensive description of the method, along with detailed derivations of the governing equations. This methodology supports solutions for traditional ablation modeling including such effects as heat transfer, material decomposition, pyrolysis gas permeation and heat exchange, and thermochemical surface erosion. The numerical scheme utilizes a control-volume approach with a variable grid to account for surface movement. This method directly supports implementation of nontraditional models such as material swelling and mechanical erosion, extending capabilities for modeling complex ablation phenomena. Verifications of the numerical implementation are provided using analytical solutions, code comparisons, and the method of manufactured solutions. These verifications are used to demonstrate solution accuracy and proper error convergence rates. A simple demonstration of a mechanical erosion (spallation) model is also provided to illustrate the unique capabilities of the method.

Ewing, Mark E.↗

NASA NTRS1985

Assessment of discretization schemes to reduce numerical diffusion in the calculation of complex flows

Comprehensive criteria for evaluating finite-difference schemes to reduce numerical diffusion in the solution of the Reynolds averaged Navier-Stokes equations have been established in this study. Several candidate schemes have been examined, based on these criteria. Taylor series analyses, model problem studies and examination of the difference equations have been used to assess the schemes. Two schemes, which showed the most promise, and satisfied a majority of the criteria, are recommended for incorporation into two dimensional and three dimensional computer codes for further evaluation.

Syed, S. A.↗

NASA NTRS1982

Application of a finite element algorithm for the solution of steady transonic Euler equations

A finite element algorithm for the solution of two-dimensional, steady Euler equations is presented which, through a Clebsch-type transformation for the velocity vector, solves the conservation of mass equation with one primary variable and two additional equations for the convection of two new variables. The accuracy and efficiency of this scheme is discussed, and the second-order accuracy attained in the analysis of the convection of the vorticity is demonstrated together with the efficient treatment of rotational and irrotational flow subregions. A sample problem is used to show the accuracy of the numerical scheme and its convergence characteristics.

Akay, H. U.↗

NASA NTRS1988

Simulation of 2-dimensional viscous flow through cascades using a semi-elliptic analysis and hybrid C-H grids

A semi-elliptic formulation, termed the interacting parabolized Navier-Stokes (IPNS) formulation, is developed for the analysis of a class of subsonic viscous flows for which streamwise diffusion is neglible but which are significantly influenced by upstream interactions. The IPNS equations are obtained from the Navier-Stokes equations by dropping the streamwise viscous-diffusion terms but retaining upstream influence via the streamwise pressure-gradient. A two-step alternating-direction-explicit numerical scheme is developed to solve these equations. The quasi-linearization and discretization of the equations are carefully examined so that no artificial viscosity is added externally to the scheme. Also, solutions to compressible as well as nearly compressible flows are obtained without any modification either in the analysis or in the solution process. The procedure is applied to constricted channels and cascade passages formed by airfoils of various shapes. These geometries are represented using numerically generated curilinear boundary-oriented coordinates forming an H-grid. A hybrid C-H grid, more appropriate for cascade of airfoils with rounded leading edges, was also developed. Satisfactory results are obtained for flows through cascades of Joukowski airfoils.

Ramamurti, R.↗

NASA NTRS2015

Structured Overlapping Grid Simulations of Contra-rotating Open Rotor Noise

Computational simulations using structured overlapping grids with the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are presented for predicting tonal noise generated by a contra-rotating open rotor (CROR) propulsion system. A coupled Computational Fluid Dynamics (CFD) and Computational AeroAcoustics (CAA) numerical approach is applied. Three-dimensional time-accurate hybrid Reynolds Averaged Navier-Stokes/Large Eddy Simulation (RANS/LES) CFD simulations are performed in the inertial frame, including dynamic moving grids, using a higher-order accurate finite difference discretization on structured overlapping grids. A higher-order accurate free-stream preserving metric discretization with discrete enforcement of the Geometric Conservation Law (GCL) on moving curvilinear grids is used to create an accurate, efficient, and stable numerical scheme. The aeroacoustic analysis is based on a permeable surface Ffowcs Williams-Hawkings (FW-H) approach, evaluated in the frequency domain. A time-step sensitivity study was performed using only the forward row of blades to determine an adequate time-step. The numerical approach is validated against existing wind tunnel measurements.

Grid↗

NASA NTRS1992

Protostellar hydrodynamics: Constructing and testing a spatially and temporally second-order-accurate method. I - Spherical coordinates

Two related numerical schemes for calculating the 3D collapse of protostellar clouds are defined, developed, and checked on a wide variety of test problems in spherical symmetry and multiple dimensions. One scheme is first-order accurate in time (code S), and the other second-order accurate in time (code ST). Through convergence testing, the codes are shown to be second-order accurate in spatial differences. Compared with the previous 3D code, the combination of reduced numerical dissipation through second-order accuracy and of removing the systematic bias toward central concentrations implies that the tendency for fragmentation into binary or multiple protostars should increase. A reinvestigation of fragmentation as a mechanism for forming binary stars is expected to yield an even more favorable evaluation.

Boss, Alan P.↗

NASA NTRS1981

Measurements and predictions of turbulence in sound excited jets

The modifications to the structure of turbulence in a round jet excited by acoustic tones upstream of the jet exit are measured and compared with predictions. The prediction scheme assumes that the acoustic waves excite instability waves or large turbulent structures in the jet. The presence of these waves at finite amplitude causes an increase in the levels of the random turbulence. Models are presented for the interaction between the instability waves and the turbulence. The numerical scheme permits the enhanced turbulence and large turbulent structures to modify the shape of the mean velocity profile. The changes in the turbulence spectra with excitation are measured. The jet has a Reynolds number of 3.75 x 10 to the 5th and the initial boundary layer is turbulent.

Morris, P. J.↗

NASA NTRS1977

Turbulent plane Couette flow using probability distribution functions

A numerical scheme employing a combination of the discrete ordinate method and finite differences is developed for solving the one-dimensional form of Lundgren's (1967) model equation for turbulent plane Couette flow. The approach used requires no a priori assumption about the form of the turbulent distribution function, and the numerical solution is obtained directly from the governing differential equations. Two different types of boundary conditions (zero-gradient and Chapman-Enskog) for the distribution function are evaluated by comparing the numerical results with experimental data. It is found that: (1) the present approach gives convergent and stable results over a wide range of Reynolds numbers; (2) Lundgren's equation yields results that compare well with experimental data for mean velocity and skin friction in the case of simple Couette flow; (3) the zero-gradient boundary condition leads to a logarithmic flow profile; and (4) the Chapman-Enskog boundary condition provides very good agreement with experimental data when applied within the near-wall region.

Srinivasan, R.↗

NASA NTRS1973

The use of Levy-Lees variables in three-dimensional boundary-layer flows

A method for solving a general class of three-dimensional boundary layer flows is developed. In the development, Levy-Lees variables are extended to three dimensions and equations are placed in these similarity variables. An implicit finite difference scheme which is stable for negative transverse velocities is used to solve these equations. The method developed is applied to obtain solutions for sharp and spherically blunted circular cones at angle of attack. Longitudinal and transverse distributions are presented for these cases. Good agreement is found with the results obtained by other numerical schemes and the experimental data of Tracy, for sharp circular cones at angle of attack. For spherically blunted cones at angle of attack, the results are in good agreement with axisymmetric sphere results up to the region where spherical symmetry holds.

Vatsa, V. N.↗

NASA NTRS1990

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. Part 1: The ODE connection and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Yee, H. C.↗

NASA NTRS1991

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. I - The dynamics of time discretization and its implications for algorithm development in computational fluid dynamics

Yee, H. C.↗

NASA NTRS2001

Asymmetry Effects in Numerical Simulation of Supersonic Flows with Upstream Separated Regions

The present paper studies the numerical simulation of flows with shock/boundary-layer upstream interaction, under conditions of symmetry in geometry, boundary conditions, and grid. For this purpose, a series of two- and three-dimensional numerical test-cases were carried out. The tests showed that standard numerical schemes, which appear to be symmetry preserving under most flow configurations, produce nonsymmetric perturbations when large separated regions are present. These perturbations are amplified when the core flow is under compression. If the flow-blockage due to separation is sufficiently large, the symmetry of the flow may collapse altogether. Experimental evidence of this numerical behavior is also considered.

Riodriguez, C. G.↗

NASA NTRS2001

Asymmetry Effects in Numerical Simulation of Supersonic Flows with Upstream Separated Regions

The present paper studies the numerical simulation of flows with shock/boundary-layer upstream interaction, under conditions of symmetry in geometry, boundary conditions, and grid. For this purpose, a series of two- and three-dimensional numerical test-cases were carried out. The tests showed that standard numerical schemes, which appear to be symmetry-preserving under most flow configurations, produce nonsymmetric perturbations when large separated regions are present. These perturbations are amplified when the core flow is under compression. If the flow-blockage due to separation is sufficiently large, the symmetry of the flow may collapse altogether. Experimental evidence of this numerical behavior is also considered.

Rodriguez, C. G.↗

NASA NTRS2003

Computation of Pressurized Gas Bearings Using CE/SE Method

The space-time conservation element and solution element (CE/SE) method is extended to compute compressible viscous flows in pressurized thin fluid films. This numerical scheme has previously been used successfully to solve a wide variety of compressible flow problems, including flows with large and small discontinuities. In this paper, the method is applied to calculate the pressure distribution in a hybrid gas journal bearing. The formulation of the problem is presented, including the modeling of the feeding system. the numerical results obtained are compared with experimental data. Good agreement between the computed results and the test data were obtained, and thus validate the CE/SE method to solve such problems.

Cioc, Sorin↗

NASA NTRS1996

Progress Toward an Efficient and General CFD Tool for Propulsion Design/Analysis

The simulation of propulsive flows inherently involves chemical activity. Recent years have seen substantial strides made in the development of numerical schemes for reacting flowfields, in particular those involving finite-rate chemistry. However, finite-rate calculations are computationally intensive and require knowledge of the actual kinetics, which are not always known with sufficient accuracy. Alternatively, flow simulations based on the assumption of local chemical equilibrium are capable of obtaining physically reasonable results at far less computational cost. The present study summarizes the development of efficient numerical techniques for the simulation of flows in local chemical equilibrium, whereby a 'Black Box' chemical equilibrium solver is coupled to the usual gasdynamic equations. The generalization of the methods enables the modelling of any arbitrary mixture of thermally perfect gases, including air, combustion mixtures and plasmas. As demonstration of the potential of the methodologies, several solutions, involving reacting and perfect gas flows, will be presented. Included is a preliminary simulation of the SSME startup transient. Future enhancements to the proposed techniques will be discussed, including more efficient finite-rate and hybrid (partial equilibrium) schemes. The algorithms that have been developed and are being optimized provide for an efficient and general tool for the design and analysis of propulsion systems.

Cox, C. F.↗

NASA NTRS1991

Reliable numerical computation in an optimal output-feedback design

This paper presents a reliable algorithm for the evaluation of a quadratic performance index and its gradients with respect to the controller design parameters. The algorithm is part of a design algorithm for optimal linear dynamic output-feedback controller that minimizes a finite-time quadratic performance index. The numerical scheme is particularly robust when it is applied to the control-law synthesis for systems with densely packed modes and where there is a high likelihood of encountering degeneracies in the closed-loop eigensystem. The algorithm has been included in a control design package for optimal robust low-order controllers. Usefulness of the proposed numerical algorithm has been demonstrated using numerous practical design cases where degeneracies occur frequently in the closed-loop system under an arbitrary controller design initialization and during the numerical search.

Vansteenwyk, Brett↗

NASA NTRS1992

Bypass transition in compressible boundary layers

Transition to turbulence in aerospace applications usually occurs in a strongly disturbed environment. For instance, the effects of free-stream turbulence, roughness and obstacles in the boundary layer strongly influence transition. Proper understanding of the mechanisms leading to transition is crucial in the design of aircraft wings and gas turbine blades, because lift, drag and heat transfer strongly depend on the state of the boundary layer, laminar or turbulent. Unfortunately, most of the transition research, both theoretical and experimental, has focused on natural transition. Many practical flows, however, defy any theoretical analysis and are extremely difficult to measure. Morkovin introduced in his review paper the concept of bypass transition as those forms of transition which bypass the known mechanisms of linear and non-linear transition theories and are currently not understood by experiments. In an effort to better understand the mechanisms leading to transition in a disturbed environment, experiments are conducted studying simpler cases, viz. the effects of free stream turbulence on transition on a flat plate. It turns out that these experiments are very difficult to conduct, because generation of free stream turbulence with sufficiently high fluctuation levels and reasonable homogeneity is non trivial. For a discussion see Morkovin. Serious problems also appear due to the fact that at high Reynolds numbers the boundary layers are very thin, especially in the nose region of the plate where the transition occurs, which makes the use of very small probes necessary. The effects of free-stream turbulence on transition are the subject of this research and are especially important in a gas turbine environment, where turbulence intensities are measured between 5 and 20 percent, Wang et al. Due to the fact that the Reynolds number for turbine blades is considerably lower than for aircraft wings, generally a larger portion of the blade will be in a laminar transitional state. The effects of large free stream turbulence in compressible boundary layers at Mach numbers are examined both in the subsonic and transonic regime using direct numerical simulations. The flow is computed over a flat plate and curved surface. while many applications operate in the transonic regime. Due the nature of their numerical scheme, a non-conservation formulation of the Navier-Stokes equations, it is a non-trivial extension to compute flow fields in the transonic regime. This project aims at better understanding the effects of large free-stream turbulence in compressible boundary layers at mach number both in the subsonic and transonic regime using direct numerical simulations. The present project aims at computing the flow over a flat plate and curved surface. This research will provide data which can be used to clarify mechanisms leading to transition in an environment with high free stream turbulence. This information is useful for the development of turbulence models, which are of great importance for CFD applications, and are currently unreliable for more complex flows, such as transitional flows.

Vandervegt, J. J.↗