Search NASA⌕ Search

SEARCH · Search NASA

Results for “high order 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Low-dissipation and -dispersion Runge-Kutta schemes for computational acoustics

In this paper, we investigate accurate and efficient time advancing methods for computational acoustics, where non-dissipative and non-dispersive properties are of critical importance. Our analysis pertains to the application of Runge-Kutta methods to high-order finite difference discretization. In many CFD applications multi-stage Runge-Kutta schemes have often been favored for their low storage requirements and relatively large stability limits. For computing acoustic waves, however, the stability consideration alone is not sufficient, since the Runge-Kutta schemes entail both dissipation and dispersion errors. The time step is now limited by the tolerable dissipation and dispersion errors in the computation. In the present paper, it is shown that if the traditional Runge-Kutta schemes are used for time advancing in acoustic problems, time steps greatly smaller than that allowed by the stability limit are necessary. Low-Dissipation and -Dispersion Runge-Kutta (LDDRE) schemes are proposed, based on an optimization that minimizes the dissipation and dispersion errors for wave propagation. Order optimizations of both single-step and two-step alternating schemes are considered. The proposed LDDRK schemes are remarkably more efficient than the classical Runge-Kutta schemes for acoustic computations. Moreover, low storage implementations of the optimized schemes are discussed. Special issues of implementing numerical boundary conditions in the LDDRK schemes are also addressed.

Hu, F. Q.↗
NASA NTRS1995

Application of low dissipation and dispersion Runge-Kutta schemes to benchmark problems in computational aeroacoustics

We investigate accurate and efficient time advancing methods for computational aeroacoustics, where non-dissipative and non-dispersive properties are of critical importance. Our analysis pertains to the application of Runge-Kutta methods to high-order finite difference discretization. In many CFD applications, multi-stage Runge-Kutta schemes have often been favored for their low storage requirements and relatively large stability limits. For computing acoustic waves, however, the stability consideration alone is not sufficient, since the Runge-Kutta schemes entail both dissipation and dispersion errors. The time step is now limited by the tolerable dissipation and dispersion errors in the computation. In the present paper, it is shown that if the traditional Runge-Kutta schemes are used for time advancing in acoustic problems, time steps greatly smaller than that allowed by the stability limit are necessary. Low Dissipation and Dispersion Runge-Kutta (LDDRK) schemes are proposed, based on an optimization that minimizes the dissipation and dispersion errors for wave propagation. Optimizations of both single-step and two-step alternating schemes are considered. The proposed LDDRK schemes are remarkably more efficient than the classical Runge-Kutta schemes for acoustic computations. Numerical results of each Category of the Benchmark Problems are presented. Moreover, low storage implementations of the optimized schemes are discussed. Special issues of implementing numerical boundary conditions in the LDDRK schemes are also addressed.

Hu, F. Q.↗
NASA NTRS2012

Comparative Study of Three High Order Schemes for LES of Temporally Evolving Mixing Layers

Three high order shock-capturing schemes are compared for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (Mc) ranging from the quasi-incompressible regime to highly compressible supersonic regime. The considered high order schemes are fifth-order WENO (WENO5), seventh-order WENO (WENO7) and the associated eighth-order central spatial base scheme with the dissipative portion of WENO7 as a nonlinear post-processing filter step (WENO7fi). This high order nonlinear filter method (H.C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22-26, 2009, Trondheim, Norway) is designed for accurate and efficient simulations of shock-free compressible turbulence, turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters. The LES results by WENO7fi using the same scheme parameter agree well with experimental results of Barone et al. (2006), and published direct numerical simulations (DNS) work of Rogers & Moser (1994) and Pantano & Sarkar (2002), whereas results by WENO5 and WENO7 compare poorly with experimental data and DNS computations.

Yee, Helen M. C.↗
NASA NTRS2016

Improving Fidelity of Launch Vehicle Liftoff Acoustic Simulations

Launch vehicles experience high acoustic loads during ignition and liftoff affected by the interaction of rocket plume generated acoustic waves with launch pad structures. Application of highly parallelized Computational Fluid Dynamics (CFD) analysis tools optimized for application on the NAS computer systems such as the Loci/CHEM program now enable simulation of time-accurate, turbulent, multi-species plume formation and interaction with launch pad geometry and capture the generation of acoustic noise at the source regions in the plume shear layers and impingement regions. These CFD solvers are robust in capturing the acoustic fluctuations, but they are too dissipative to accurately resolve the propagation of the acoustic waves throughout the launch environment domain along the vehicle. A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed to improve such liftoff acoustic environment predictions. The framework combines the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate discontinuous Galerkin (DG) solver, Loci/THRUST, developed in the same computational framework. Loci/THRUST employs a low dissipation, high-order, unstructured DG method to accurately propagate acoustic waves away from the source regions across large distances. The DG solver is currently capable of solving up to 4th order solutions for non-linear, conservative acoustic field propagation. Higher order boundary conditions are implemented to accurately model the reflection and refraction of acoustic waves on launch pad components. The DG solver accepts generalized unstructured meshes, enabling efficient application of common mesh generation tools for CHEM and THRUST simulations. The DG solution is coupled with the CFD solution at interface boundaries placed near the CFD acoustic source regions. Both simulations are executed simultaneously with coordinated boundary condition data exchange.

Liever, Peter↗
NASA NTRSDate not supplied

Efficient Basis Derivatives Evaluations for High-order Discontinuous Galerkin Schemes

Computational methods of evaluating a basis of Lagrange polynomials are developed for the purpose of implementing efficient discontinuous Galerkin conservation laws solvers. Special attention is payed to the computation of higher-order partial derivatives which may be required for certain applications. Two different approaches are considered; one involves hardcoding explicit one-line formulae into the source code in order to make it as simple as possible, whereas the other uses algorithms designed to minimize the asymptotic order of growth with respect to the order of the scheme. Timing experiments show that either approach can perform well if implemented effectively, and identify advantages and disadvantages of each.

Micaiah Smith-Pierce↗
NASA NTRSDate not supplied

Enabling Supersonic Over-Land Flight Using Computational Modeling

NASA's QueSST mission is interested in exploring the potential of commercial, over-land supersonic flight, which has been banned since 1973 due to the noise intensity of the sonic booms that supersonic vehicles create. NASA has partnered with Lockheed Martin to design and construct an experimental aircraft, the X-59, that will shape the sonic boom to achieve a perceived loudness on the ground that is lower than 75dB. In this talk, we present new developments in the Launch Ascent and Vehicle Aerodynamics (LAVA) CFD solver framework which are tailored to the analysis of supersonic vehicles, like the X-59, and sonic boom propagation. Utilizing a high-order space marching method, adjoint-based grid adaptation, and robust database capabilities have enabled significant cost savings over traditional CFD methods for the same level of accuracy in perceived boom loudness. Example cases demonstrating the toolset are presented, along with comparisons to experimental wind tunnel data. These developments will enable LAVA to quickly and accurately provide simulation data to the project during X-59's first and subsequent flights in 2024 and beyond.

CST↗
NASA NTRSDate not supplied

Enhanced Simulation Techniques in Predicting Sonic Boom Loudness Using CFD

This paper outlines advancements in predicting sonic boom loudness within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework. Traditionally, a two step process consisting of a steady state computational fluid dynamics problem for near-field analysis and a far-field propagation solver for calculation of loudness metrics has been used. Improvements to this process made in this work include utilizing a high-order space marching method for mid-field computations, developing a novel output-based mesh adaptation method targeting error in near-field pressure sig-natures, and developing a robust scripting system using curvilinear grids to increase robustness and simplify the process of running large databases of simulation cases. These advancements are detailed and applied to the simulation of the X-59, presenting comparative cost and timing analyses between the prior two step workflow and the current three step procedure. We achieve increased accuracy and robustness for loudness predictions with at least a50%computational cost reduction.

CST↗
NASA NTRS1989

Universal limiter for high order explicit conservative advection schemes

A simple method for construction of nonoscillatory, explicit, conservative advection schemes of arbitrarily high accuracy is described. Four test profiles are considered: an isolated sine-squared wave, a unit step function, a semiellipse, and a narrow Gaussian. A universal limiter guarantees monotonic resolution of the step. Sharpness increases uniformly with the order of the base scheme.

Leonard, B. P.↗
NASA NTRS1993

Protostellar hydrodynamics: Constructing and testing a spacially and temporally second-order accurate method. 2: Cartesian coordinates

In Boss & Myhill (1992) we described the derivation and testing of a spherical coordinate-based scheme for solving the hydrodynamic equations governing the gravitational collapse of nonisothermal, nonmagnetic, inviscid, radiative, three-dimensional protostellar clouds. Here we discuss a Cartesian coordinate-based scheme based on the same set of hydrodynamic equations. As with the spherical coorrdinate-based code, the Cartesian coordinate-based scheme employs explicit Eulerian methods which are both spatially and temporally second-order accurate. We begin by describing the hydrodynamic equations in Cartesian coordinates and the numerical methods used in this particular code. Following Finn & Hawley (1989), we pay special attention to the proper implementations of high-order accuracy, finite difference methods. We evaluate the ability of the Cartesian scheme to handle shock propagation problems, and through convergence testing, we show that the code is indeed second-order accurate. To compare the Cartesian scheme discussed here with the spherical coordinate-based scheme discussed in Boss & Myhill (1992), the two codes are used to calculate the standard isothermal collapse test case described by Bodenheimer & Boss (1981). We find that with the improved codes, the intermediate bar-configuration found previously disappears, and the cloud fragments directly into a binary protostellar system. Finally, we present the results from both codes of a new test for nonisothermal protostellar collapse.

Myhill, Elizabeth A.↗
NASA NTRSDate not supplied

Predicting SLS Launch Environment using a Novel Multiphase Formulation

Powerful acoustic waves generated during ignition of launch vehicles may be dangerous to the vehicle, its payload, or the surrounding structures. The water-based Ignition Overpressure and Sound Suppression (IOP/SS) system at Kennedy Space Center’s (KSC) Launch Complex 39B (LC-39B) will be used to protect the Space Launch System (SLS) from the acoustic vibrations generated during launch. The IOP/SS system uses enormous amounts of water to dampen and attenuate these sound waves. To better understand the launch environment risks and to study the effectiveness of the IOP/SS system it is desirable to have time-accurate unsteady simulations of the vehicle ignition with water-based sound suppression. This paper presents results obtained with a novel, high-order accurate, and robust numerical method designed for simulating compressible multiphase flows. A positivity-preserving finite difference scheme is utilized which is formally high-order accurate and also provably robust. Robustness is critical due to the extreme nature of the flow which exhibits highly nonlinear shock and rarefaction waves interacting with liquid-gas interfaces with density ratios of the order of 1000:1. Furthermore, the high-order accuracy (and the high resolution property) is desirable for predicting wave phenomena like IOP waves since the signal can be resolved accurately and propagated long distances with fewer grid points. This finite-difference method was developed using NASA’s Launch, Ascent, and Vehicle Aerodynamics (LAVA) Cartesian immersed boundary framework. We present a validation case by applying our solver to the SLS Scale Model Acoustic Test (SMAT). The SLS SMAT is a well-instrumented 5% scale model test meant to represent the SLS at NASA KSC’s LC-39B pad. Scale IOP tests were performed with and without the sound suppression water and included many sensors which recorded the pressure waves produced during ignition. For this validation case we conduct two simulations, likewise with and without sound suppression water, and compare the SLS SMAT pressure sensor signals with our numerical signals at identical locations. Following this validation case we present a study of the SLS launch environment to examine engineering safety concerns about the mobile launch pad. Engineers at KSC redesigned the main flame deflector at LC-39B anticipating the increased loads from the SLS and to repair damage from prior Shuttle missions. This deflector redesign made use of surface pressure and temperature data from LAVA full-scale SLS simulations without the sound suppression system. The engineers were questioning the possibility of increased pressure loads on the underside of the mobile launcher due to the water in the flame trench. Based on the results established in our simulations of the SLS SMAT, we performed updated calculations for SLS at LC-39B with and without water systems active to assess the readiness of the launch pad for Artemis I launch. Our results show that the IOP/SS system is effective at reducing the overpressure signal and overall sound pressure levels felt by the vehicle and additionally that the pressure loads experienced by the mobile launcher (ML) during engine startup is not increased by the presence of water.

EGS↗
NASA NTRS1993

Nonrecursive formulations of multibody dynamics and concurrent multiprocessing

Since the late 1980's, research in recursive formulations of multibody dynamics has flourished. Historically, much of this research can be traced to applications of low dimensionality in mechanism and vehicle dynamics. Indeed, there is little doubt that recursive order N methods are the method of choice for this class of systems. This approach has the advantage that a minimal number of coordinates are utilized, parallelism can be induced for certain system topologies, and the method is of order N computational cost for systems of N rigid bodies. Despite the fact that many authors have dismissed redundant coordinate formulations as being of order N(exp 3), and hence less attractive than recursive formulations, we present recent research that demonstrates that at least three distinct classes of redundant, nonrecursive multibody formulations consistently achieve order N computational cost for systems of rigid and/or flexible bodies. These formulations are as follows: (1) the preconditioned range space formulation; (2) penalty methods; and (3) augmented Lagrangian methods for nonlinear multibody dynamics. The first method can be traced to its foundation in equality constrained quadratic optimization, while the last two methods have been studied extensively in the context of coercive variational boundary value problems in computational mechanics. Until recently, however, they have not been investigated in the context of multibody simulation, and present theoretical questions unique to nonlinear dynamics. All of these nonrecursive methods have additional advantages with respect to recursive order N methods: (1) the formalisms retain the highly desirable order N computational cost; (2) the techniques are amenable to concurrent simulation strategies; (3) the approaches do not depend upon system topology to induce concurrency; and (4) the methods can be derived to balance the computational load automatically on concurrent multiprocessors. In addition to the presentation of the fundamental formulations, this paper presents new theoretical results regarding the rate of convergence of order N constraint stabilization schemes associated with the newly introduced class of methods.

Kurdila, Andrew J.↗
NASA NTRS2013

Estimate Low and High Order Wavefront Using P1640 Calibrator Measurements

P1640 high contrast imaging system on the Palomar 200 inch Telescope consists of an apodized-pupil Lyot coronagraph, the PALM-3000 adaptive optics (P3K-AO), and P1640 Calibrator (CAL). Science images are recorded by an integral field spectrograph covering J-H bands for detecting and characterizing stellar companions. With aberrations from atmosphere corrected by the P3K-AO, instrument performance is limited mainly by the quasi-static speckles due to noncommon path wavefront aberrations for the light to propagate to the P3K-AO wavefront sensor and to the coronagraph mask. The non-common path wavefront aberrations are sensed by CAL, which measures the post-coronagraph E-field using interferometry, and can be effectively corrected by offsetting the P3K-AO deformable mirror target position accordingly. Previously, we have demonstrated using CAL measurements to correct high order wavefront aberrations, which is directly connected to the static speckles in the image plane. Low order wavefront, on the other hand, usually of larger amplitudes, causes light to leak through the coronagraph making the whole image plane brighter. Knowledge error in low order wavefront aberrations can also affect the estimation of the high order wavefront. Even though, CAL is designed to sense efficiently high order wavefront aberrations, the low order wavefront front can be inferred with less sensitivity. Here, we describe our method for estimating both low and high order wavefront aberrations using CAL measurements by propagating the post-coronagraph E-field to a pupil before the coronagraph. We present the results from applying this method to both simulated and experiment data.

aberrations↗
NASA NTRS2002

Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids II: Extension to Two Dimensional Scalar Equation

The framework for constructing a high-order, conservative Spectral (Finite) Volume (SV) method is presented for two-dimensional scalar hyperbolic conservation laws on unstructured triangular grids. Each triangular grid cell forms a spectral volume (SV), and the SV is further subdivided into polygonal control volumes (CVs) to supported high-order data reconstructions. Cell-averaged solutions from these CVs are used to reconstruct a high order polynomial approximation in the SV. Each CV is then updated independently with a Godunov-type finite volume method and a high-order Runge-Kutta time integration scheme. A universal reconstruction is obtained by partitioning all SVs in a geometrically similar manner. The convergence of the SV method is shown to depend on how a SV is partitioned. A criterion based on the Lebesgue constant has been developed and used successfully to determine the quality of various partitions. Symmetric, stable, and convergent linear, quadratic, and cubic SVs have been obtained, and many different types of partitions have been evaluated. The SV method is tested for both linear and non-linear model problems with and without discontinuities.

Wang, Z. J.↗
DOE OSTI2025

Explicit block encodings of boundary value problems for many-body elliptic operators

Simulation of physical systems is one of the most promising use cases of future digital quantum computers. In this work we systematically analyze the quantum circuit complexities of block encoding the discretized elliptic operators that arise extensively in numerical simulations for partial differential equations, including high-dimensional instances for many-body simulations. When restricted to rectangular domains with separable boundary conditions, we provide explicit circuits to block encode the many-body Laplacian with separable periodic, Dirichlet, Neumann, and Robin boundary conditions, using standard discretization techniques from low-order finite difference methods. To obtain high-precision, we introduce a scheme based on periodic extensions to solve Dirichlet and Neumann boundary value problems using a high-order finite difference method, with only a constant increase in total circuit depth and subnormalization factor. We then present a scheme to implement block encodings of differential operators acting on more arbitrary domains, inspired by Cartesian immersed boundary methods. We then block encode the many-body convective operator, which describes interacting particles experiencing a force generated by a pair-wise potential given as an inverse power law of the interparticle distance. This work provides concrete recipes that are readily translated into quantum circuits, with depth logarithmic in the total Hilbert space dimension, that block encode operators arising broadly in applications involving the quantum simulation of quantum and classical many-body mechanics.

Kharazi, Tyler [University of California, Berkeley↗
NASA NTRS1997

Computational Considerations for the Simulation of Discontinuous Flows

The numerical study of aeroacoustic problems places stringent demands on the choice of a computational algorithm, because it requires the ability to propagate disturbances of small amplitude and short wavelength. The demands are particularly high when shock waves are involved, because the chosen algorithm must also resolve discontinuities in the solution. In a previous work the capabilities and deficiencies of shock-capturing methods for aeroacoustic problems were demonstrated using a high-order essentially nonoscillatory (ENO) numerical method. It was shown that first-order results are obtained when simulating time-dependent flows with discontinuities. The present study reaffirms this conclusion by comparing the ENO results with those obtained using a conventional linear scheme. A sixth-order-accurate compact implicit finite difference scheme is used to investigate various discontinuous flows. The design order of accuracy is achieved in the smooth regions of a steady-state, quasi-one-dimensional Euler test case, as well as in the time-dependent Burgers' equation. However, in the unsteady Euler sound-shock interaction, first-order results are obtained downstream of the shock. A comparison is made between the linear and nonlinear results, noting the advantages of each method. A discontinuous linear model problem is then used to identify the cause of the first-order results. Here, the nature of the solution error is quantified as being predominantly a numerical phase shift, and a post-processing procedure is demonstrated which increases the solution accuracy downstream of the discontinuity to second-order.

Carpenter Mark H.↗
NASA NTRS2018

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗