Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Runge Kutta Methods”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Three-dimensional unstructured method for flows past bodies in 6-DOF relative motion

A three dimensional, unstructured-mesh methodology was developed to simulate unsteady flows past bodies in relative motion, where the trajectory was determined from the instantaneous aerodynamics. The method coupled the equations of fluid flow and those of rigid-body dynamics, and captured the time-dependent interference between stationary and moving boundaries. The unsteady, compressible Euler equations were solved on dynamic, unstructured meshes by an explicit, finite-volume, upwind method. The grid adaptation was performed within a window placed around the moving body. The Euler equations of dynamics were solved by a Runge-Kutta integration scheme. The flow solver and the adaptation scheme were validated by simulating the transonic, unsteady flow around a wing undergoing a forced, periodic pitching motion, then comparing the results with the experimental data. To validate the trajectory code, the six-degrees-of-freedom (DOF) motion of a store separating from a wing was computed using the experimentally determined force and moment fields, then comparing with an independently generated trajectory. Finally, the overall methodology was demonstrated by simulating the unsteady flowfield and the trajectory of a store dropped from a wing. The methodology, its computational cost notwithstanding, has proven to be accurate, automated, easy for dynamic gridding, and relatively efficient for the required man-hours.

Singh, K. P.↗
NASA NTRS1999

Low-Storage, Explicit Runge-Kutta Schemes for the Compressible Navier-Stokes Equations

The derivation of storage explicit Runge-Kutta (ERK) schemes has been performed in the context of integrating the compressible Navier-Stokes equations via direct numerical simulation. Optimization of ERK methods is done across the broad range of properties, such as stability and accuracy efficiency, linear and nonlinear stability, error control reliability, step change stability, and dissipation/dispersion accuracy, subject to varying degrees of memory economization. Following van der Houwen and Wray, 16 ERK pairs are presented using from two to five registers of memory per equation, per grid point and having accuracies from third- to fifth-order. Methods have been assessed using the differential equation testing code DETEST, and with the 1D wave equation. Two of the methods have been applied to the DNS of a compressible jet as well as methane-air and hydrogen-air flames. Derived 3(2) and 4(3) pairs are competitive with existing full-storage methods. Although a substantial efficiency penalty accompanies use of two- and three-register, fifth-order methods, the best contemporary full-storage methods can be pearl), matched while still saving two to three registers of memory.

Kennedy, Chistopher A.↗
NASA NTRS2001

Textbook Multigrid Efficiency for Computational Fluid Dynamics Simulations

Considerable progress over the past thirty years has been made in the development of large-scale computational fluid dynamics (CFD) solvers for the Euler and Navier-Stokes equations. Computations are used routinely to design the cruise shapes of transport aircraft through complex-geometry simulations involving the solution of 25-100 million equations; in this arena the number of wind-tunnel tests for a new design has been substantially reduced. However, simulations of the entire flight envelope of the vehicle, including maximum lift, buffet onset, flutter, and control effectiveness have not been as successful in eliminating the reliance on wind-tunnel testing. These simulations involve unsteady flows with more separation and stronger shock waves than at cruise. The main reasons limiting further inroads of CFD into the design process are: (1) the reliability of turbulence models; and (2) the time and expense of the numerical simulation. Because of the prohibitive resolution requirements of direct simulations at high Reynolds numbers, transition and turbulence modeling is expected to remain an issue for the near term. The focus of this paper addresses the latter problem by attempting to attain optimal efficiencies in solving the governing equations. Typically current CFD codes based on the use of multigrid acceleration techniques and multistage Runge-Kutta time-stepping schemes are able to converge lift and drag values for cruise configurations within approximately 1000 residual evaluations. An optimally convergent method is defined as having textbook multigrid efficiency (TME), meaning the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in the discretized system of equations (residual equations). In this paper, a distributed relaxation approach to achieving TME for Reynolds-averaged Navier-Stokes (RNAS) equations are discussed along with the foundations that form the basis of this approach. Because the governing equations are a set of coupled nonlinear conservation equations with discontinuities (shocks, slip lines, etc.) and singularities (flow- or grid-induced), the difficulties are many. This paper summarizes recent progress towards the attainment of TME in basic CFD simulations.

Brandt, Achi↗
NASA NTRS2002

Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think of a mechanism that activates the split form of the equations only at some parts of the domain. Another issue is how to define good sensors for determining in which parts of the computational domain a certain feature should be filtered by the appropriate numerical dissipation. For the present study we employ a wavelet technique introduced in as sensors. Here, the method is briefly described with selected numerical experiments.

Sjoegreen, Bjoern↗
NASA NTRS1989

Design of optimally smoothing multi-stage schemes for the Euler equations

In this paper, a method is developed for designing multi-stage schemes that give optimal damping of high-frequencies for a given spatial-differencing operator. The objective of the method is to design schemes that combine well with multi-grid acceleration. The schemes are tested on a nonlinear scalar equation, and compared to Runge-Kutta schemes with the maximum stable time-step. The optimally smoothing schemes perform better than the Runge-Kutta schemes, even on a single grid. The analysis is extended to the Euler equations in one space-dimension by use of 'characteristic time-stepping', which preconditions the equations, removing stiffness due to variations among characteristic speeds. Convergence rates independent of the number of cells in the finest grid are achieved for transonic flow with and without a shock. Characteristic time-stepping is shown to be preferable to local time-stepping, although use of the optimally damping schemes appears to enhance the performance of local time-stepping. The extension of the analysis to the two-dimensional Euler equations is hampered by the lack of a model for characteristic time-stepping in two dimensions. Some results for local time-stepping are presented.

Van Leer, Bram↗
NASA NTRS1992

Numerical simulations using a dynamic solution-adaptive grid algorithm, with applications to unsteady internal flows

An investigation into the numerical simulation of unsteady flows is undertaken using a two-stage Runge-Kutta scheme coupled with the dynamic solution-adaptive grid algorithm developed by the authors. The inviscid fluxes are described by a modified Advective Upwind Split Method to eliminate the need for artificial dissipation. A well-documented numerical example containing moving discontinuities is presented that demonstrates the ability of the coupled grid/solver scheme to accurately capture unsteady flowfield phenomena. Applications are to a typical inlet diffuser configuration at Mach 3.0 with excessive back pressure inducing inlet unstart.

Benson, Rusty A.↗
NASA NTRS1993

Time-accurate simulation of a self-excited oscillatory supersonic external flow with a multi-block solution-adaptive mesh algorithm

Results are presented of an investigation of the time-accurate simulation of supersonic unsteady flow oscillations over spike-tipped bodies using the multistage Runge-Kutta scheme coupled with a dynamic solution-adaptive grid algorithm modified for multiblock capabilities. The inviscid fluxes are described by a modified advective upwind split method to obviate the need for artificial dissipation. If a time-varying, solution-adaptive mesh algorithm is incorporated, resolution of the details of the unsteady spike-tipped body flow is improved. The adaptive algorithm is also shown to resolve multiple and diverse features of the flow simultaneously, with the adapted regions in the mesh convecting with these features as they translate.

Ingram, Clint L.↗
NASA NTRS1986

Computation of three-dimensional, rotational flow through turbomachinery blade rows for improved aerodynamic design studies

The performance of a three dimensional computer code developed for predicting the flowfield in stationary and rotating turbomachinery blade rows is described in this study. The four stage Runge-Kutta numerical integration scheme is used for solving the governing flow equations and yields solution to the full, three dimensional, unsteady Euler equations in cylindrical coordinates. This method is fully explicit and uses the finite volume, time marching procedure. In order to demonstrate the accuracy and efficiency of the code, steady solutions were obtained for several cascade geometries under widely varying flow conditions. Computed flowfield results are presented for a fully subsonic turbine stator and a low aspect ratio, transonic compressor rotor blade under maximum flow and peak efficiency design conditions. Comparisons with Laser Anemometer measurements and other numerical predictions are also provided to illustrate that the present method predicts important flow features with good accuracy and can be used for cost effective aerodynamic design studies.

Subramanian, S. V.↗
NASA NTRS1991

An effective multigrid method for high-speed flows

The use is considered of a multigrid method with central differencing to solve the Navier-Stokes equations for high speed flows. The time dependent form of the equations is integrated with a Runge-Kutta scheme accelerated by local time stepping and variable coefficient implicit residual smoothing. Of particular importance are the details of the numerical dissipation formulation, especially the switch between the second and fourth difference terms. Solutions are given for 2-D laminar flow over a circular cylinder and a 15 deg compression ramp.

Swanson, R. C.↗
NASA NTRS1992

An effective multigrid method for high-speed flows

The use is considered of a multigrid method with central differencing to solve the Navier-Stokes equations for high speed flows. The time dependent form of the equations is integrated with a Runge-Kutta scheme accelerated by local time stepping and variable coefficient implicit residual smoothing. Of particular importance are the details of the numerical dissipation formulation, especially the switch between the second and fourth difference terms. Solutions are given for 2-D laminar flow over a circular cylinder and a 15 deg compression ramp.

Swanson, R. C.↗
NASA NTRS1997

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

An investigation is conducted on several numerical schemes for use in the computation of two-dimensional, spatially evolving, laminar variable-density compressible shear layers. Schemes with various temporal accuracies and arbitrary spatial accuracy for both inviscid and viscous terms are presented and analyzed. All integration schemes use explicit or compact finite-difference derivative operators. Three classes of schemes are considered: an extension of MacCormack's original second-order temporally accurate method, a new third-order variant of the schemes proposed by Rusanov and by Kutier, Lomax, and Warming (RKLW), and third- and fourth-order Runge-Kutta schemes. In each scheme, stability and formal accuracy are considered for the interior operators on the convection-diffusion equation U(sub t) + aU(sub x) = alpha U(sub xx). Accuracy is also verified on the nonlinear problem, U(sub t) + F(sub x) = 0. Numerical treatments of various orders of accuracy are chosen and evaluated for asymptotic stability. Formally accurate boundary conditions are derived for several sixth- and eighth-order central-difference schemes. Damping of high wave-number data is accomplished with explicit filters of arbitrary order. Several schemes are used to compute variable-density compressible shear layers, where regions of large gradients exist.

Kennedy, Christopher A.↗
NASA NTRSDate not supplied

Explicit Discontinuous Galerkin Methods for Conservation Laws

The two explicit DG methods in this study are based on a ‘predictor-corrector’ formulation, the first introduced by Lörcher, Gassner, and Munz (2007, 2008) called space–time expansion discontinuous Galerkin or STE-DG scheme, and the second, introduced independently by the author (Huynh 2006, 2013) called the upwind moment scheme. The predictor step of the two methods is essentially identical using a Cauchy-Kovalevsky (CK) procedure, which involves no interaction of the data among neighboring cells. The corrector step also shares the same space-time integration formulation and is where interaction of the data among neighboring cells takes place; the difference, however, is in how the resulting space-time volume integral is estimated. As a consequence of the different estimates, for the case of advection in one spatial dimension (1D), the moment scheme has a CFL (Courant-Friedrichs-Lewy) condition of 1 for all p and is accurate to order 2p+1, i.e., it possesses the super accuracy property, whereas the STE-DG method has a more restrictive CFL condition and is accurate to the expected order of p+1. For 1D advection, compared with the CFL conditions of 1/(2p+1) of standard RK-DG (Runge-Kutta) scheme where space and time discretization are of the same order, the moment scheme allows a significantly larger time step size. It also turns out that the scheme yields a result identical to Van Leer’s scheme III (1977), which amounts to shifting the data a distance of advection corresponding to the time step and projecting the result onto the space of polynomial solutions. Contrary to Van Leer’s approach, however, the space-time ‘predictor-corrector’ formulation facilitates extensions to the case of systems of equations. Concerning 2D extensions, in the case of advection, when the flow is along the diagonal direction, the CFL conditions for the moment schemes become restrictive as will be shown by Fourier (Von Neumann) stability and accuracy analyses. Since the moment scheme employs the right Radau points as collocation points in time, the method is closely related to the implicit Radau IIA scheme, which is stable for any time step size. The role of Radau IIA in relieving stability restriction for these explicit DG schemes remains to be explored

Discontinuous Galerkin↗
NASA NTRS1994

Some considerations on velocity vector accuracy in dust trajectory analysis

The relative contributions of comets and asteroids to the reservoir of dust in the interplanetary medium is not known. There are direct observations of dust released from comets and there is evidence to associate the IRAS dust bands with possible collisions of asteroids in the main belt. A means towards sorting out the parent sources has been proposed in the establishment of a dust collector in orbit about the Earth. The purpose of such a facility would be to collect not only cosmic dust particles intact but also the state vectors, as they arrive at the detector, the idea being that one may combine analytical laboratory analysis of the physics and chemistry of the captured particles with orbital data in order to help distinguish between bodies and identify parent bodies. The theoretical study of dust particle orbits in the solar system takes on greatly more importance if we use collected trajectory data. The orbital motion of dust when radiation and forces alone are acting is well understood. When gravitational forces due to the planets are included, the motion can become quite complex. In order to characterize the orbits of particles as they crossed the Earth's orbits, a study of the long-time dust orbital evolution was undertaken. We have considered various parameters associated with these dust orbits to see if one may in a general way discriminate between particles evolved from comets and asteroids. We proceed in this study as we have done previously. That is, we considered the dust particles as ideal black bodies, of density 1 gm/cc, spherical, with radii 10-100 microns. Particles of this size are affected by radiation forces, photon pressure, and Poynting-Robertson drag. Account was also taken of solar wind drag, which amounts to about 30 percent of the Poynting-Robertson drag negligible. The gravitational forces due to the planets are included, unlike in our previous study; the planetary orbits are those of true n-body interaction so that the possibility of secular resonance is included. Our method was to calculate explicitly by a numerical procedure the orbits of dust particles after they left their parent bodies. The motion is determined numerically with the implicit Runge-Kutta integrator using Gauss-Radau spacings.

Jackson, A. A.↗
DOE OSTI2023

Exponential Runge-Kutta Parareal for non-diffusive equations

Parareal is a well-known parallel-in-time algorithm that combines a coarse and fine propagator within a parallel iteration. It allows for large-scale parallelism that leads to significantly reduced computational time compared to serial time-stepping methods. However, like many parallel-in-time methods it can fail to converge when applied to non-diffusive equations such as hyperbolic systems or dispersive nonlinear wave equations. Here, this paper explores the use of exponential integrators within the Parareal iteration. Exponential integrators are particularly interesting candidates for Parareal because of their ability to resolve fast-moving waves, even at the large stepsizes used by coarse propagators. This work begins with an introduction to exponential Parareal integrators followed by several motivating numerical experiments involving the nonlinear Schrödinger equation. These experiments are then analyzed using linear analysis that approximates the stability and convergence properties of the exponential Parareal iteration on nonlinear problems. The paper concludes with two additional numerical experiments involving the dispersive Kadomtsev-Petviashvili equation and the hyperbolic Vlasov-Poisson equation. These experiments demonstrate that exponential Parareal methods offer improved time-to-solution compared to serial exponential integrators when solving certain non-diffusive equations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Computational techniques for high-speed flows with viscous and chemical effects

Algorithms for solving the Euler and the Navier-Stokes equations in conjunction with chemical kinetic equations are presented. The convective flux is estimated from a quasi one dimensional interpolation procedure. Shock, contact, and expansion waves and thermochemical nonequilibrium phenomena are captured by the Lax-Friedrichs technique. Relaxation techniques were developed to enhance their effectiveness in dealing with spatial and temporal stiffness associated with the physical problems. Both explicit and implicit smoothers were implemented into the standard multigrid time stepping method. Unsteady and steady scalar problems are discussed. A perfect gas and equilibrium air shock tube problem is investigated. Numerical schemes and techniques are compared for the problems of shock and boundary layer interaction and three dimensional viscous, nonequilibrium flow encompassing an aerobrake. The results are comparable in accuracy against other high order non-oscillatory techniques. The multigrid methods are assessed using a Mach 8 flow over a complete planar body, a sphere, and a blunt delta wing at 20 deg incidence. Applying an implicit multigrid method on a nested grid of 128 by 64 nodes the reduction factor is 0.25. The central processing unit reduction factor is 2.2 after both the single and multigrid Runge-Kutta solutions converged to machine zero on a grid of 37 by 41 by 73 nodes.

Li, C. P.↗
NASA NTRS1993

Computation of propulsion-related flowfields using unstructured adaptive meshes

In this paper, we describe a computational fluid dynamics (CFD) technique based on unstructured triangular/tetrahedral meshes. A finite-volume scheme is used in conjunction with a multi-stage Runge-Kutta algorithm. Convergence enhancements in the form of dual time-stepping and time-derivative preconditioning are used to overcome the limitations of conventional multi-stage schemes. The method is applied to propulsion-related flows and shown to perform satisfactorily.

Weiss, Jonathan M.↗
NASA NTRSDate not supplied

Discretization Error Estimation and Control for Farfield Acoustic Signatures

We investigate the utility of adjoint-based error estimates for sonic boom farfield simulations governed by solutions of the augmented Burgers’ equation. Solution of this nonlinear system uses operator splitting with a second-order finite volume discretization in space and second-order Runge-Kutta time marching, while the absorption and molecular relaxation are solved using second-order central differencing. The discretization error in selected ground sonic boom cost functionals is estimated using the method of adjoint-weighted residuals. Key elements of the implementation process are emphasized with details provided on the practical aspects as appliedto the sonic boom farfield propagation. We establish the accuracy of the adjoint solutions usingcomplex step and finite difference approaches, and examine the accuracy of the error estimates using analytical N-wave solutions. We then apply it to a pressure waveform corresponding to the X-59 research aircraft. The investigations demonstrate that the method of adjoint-weighted residuals accurately predicts the level of discretization error present in sonic boom farfield simulations while offering insight into which features of the near field signal are the primary drivers of ground noise metrics. The numerical results indicate that at sampling frequencies as low as50kHz, discretization error in the propagation is under 0.01 dB[A] for realistically complex examples.

CST↗