Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high order”

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 NTRS2001

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible; (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred; (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used; and (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for nonlinear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows. The minimization of employing very fine grids to overcome the production of spurious numerical solutions and/or instability due to under-resolved grids is also sought. The incremental studies to illustrate the performance of the approach are summarized. Extensive testing and full implementation of the approach is forthcoming. The results shown so far are very encouraging.

Yee, Helen C.↗
NASA NTRS2004

High-Order Shock-Capturing Methods for Modeling Dynamics of the Solar Atmosphere

We use one-dimensional high-order central shock capturing numerical methods to study the response of various model solar atmospheres to forcing at the solar surface. The dynamics of the atmosphere is modeled with the Euler equations in a variable-sized flux tube in the presence of gravity. We study dynamics of the atmosphere suggestive of spicule formation and coronal oscillations. These studies are performed on observationally-derived model atmospheres above the quiet sun and above sunspots. To perform these simulations, we provide a new extension of existing second- and third- order shock-capturing methods to irregular grids. We also solve the problem of numerically maintaining initial hydrostatic balance via the introduction of new variables in the model equations and a careful initialization mechanism. We find several striking results: all model atmospheres respond to a single impulsive perturbation with several strong shock waves consistent with the rebound-shock model. These shock waves lift material and the transition region well into the initial corona, and the sensitivity of this lift to the initial impulse depends non-linearly on the details of the atmosphere model. We also reproduce an observed 3-minute coronal oscillation above sunspots compared to 5-minute oscillations above the quiet sun.

Bryson, Steve↗
NASA NTRS2005

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible. (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred. (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used. (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for non-linear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows.

Yee, H. C.↗
NASA NTRS2019

High-Order Finite-Difference Nonlinear Filter Methods for Subsonic Turbulence Simulation with Stochastic Forcing

Numerical stability of high-order filter schemes developed by Yee & Sjogreen is tested on three-dimensional turbulence simulations with stochastic forcing and their performance is compared with that of TVD and WENO schemes. The best­ performing filter method employs an eighth-order central base scheme with the Kennedy & Gruber skew-symmetric splitting of the inviscid flux derivative, a wavelet-based local flow sensor, a nonlinear filter utilizing the dissipative portion of seventh-order \VENO scheme, and an explicit third - or fourth-order Runge­-Kutta time integration. We show that the filter scheme is more computational]y efficient and provides a wider spectral bandwidth compared to the seventh-order WENO scheme. The method also demonstrates robust long-time integration for moderately compressible turbulence. In contrast, the fifth- and seventh­ order WENO schemes show non-trivial evolution of the velocity and density power spectra. over a. few dozen dynamical times, where both TVD and filter schemes recover a so lid statistically stationary turbulent state

Kritsuk, Alexei G.↗
NASA NTRS1990

High-order state space simulation models of helicopter flight mechanics

This paper describes the formulation and validation of a high-order linearized mathematical model of helicopter flight mechanics, which includes rotor flap and lag degrees of freedom as well as inflow dynamics. The model is extracted numerically from an existing nonlinear, blade element, time simulation model. Extensive modifications in the formulation and solution process of the nonlinear model, required for a theoretically rigorous linearization, are described in detail. The validation results show that the linearized model successfully captures the coupled rotor-fuselage dynamics in the frequency band most critical for the design of advanced flight control systems. Additional results quantify the extent to which the order of the model can be reduced without loss of fidelity.

Kim, Frederick D.↗
NASA NTRS1993

High-order state space simulation models of helicopter flight mechanics

This paper describes the formulation and validation of a high-order linearized mathematical model of helicopter flight mechanics, which includes rotor flap and lag degrees of freedom as well as inflow dynamics. The model is extracted numerically from an existing nonlinear, blade element, real-time simulation model. Extensive modifications in the formulation and solution process of the nonlinear model, required for a theoetically rigorous linearization, are described in detail. The validation results show that the linearized model successfully captures the coupled rotor-fuselage dynamics in the frequency band most critical for the design of advanced flight control systems. Additional results quantify the extent to which the order of the model can be reduced without loss of fidelity.

Kim, Frederick D.↗
NASA NTRS1996

A Review of High-Order and Optimized Finite-Difference Methods for Simulating Linear Wave Phenomena

This paper presents a review of high-order and optimized finite-difference methods for numerically simulating the propagation and scattering of linear waves, such as electromagnetic, acoustic, or elastic waves. The spatial operators reviewed include compact schemes, non-compact schemes, schemes on staggered grids, and schemes which are optimized to produce specific characteristics. The time-marching methods discussed include Runge-Kutta methods, Adams-Bashforth methods, and the leapfrog method. In addition, the following fourth-order fully-discrete finite-difference methods are considered: a one-step implicit scheme with a three-point spatial stencil, a one-step explicit scheme with a five-point spatial stencil, and a two-step explicit scheme with a five-point spatial stencil. For each method studied, the number of grid points per wavelength required for accurate simulation of wave propagation over large distances is presented. Recommendations are made with respect to the suitability of the methods for specific problems and practical aspects of their use, such as appropriate Courant numbers and grid densities. Avenues for future research are suggested.

Zingg, David W.↗
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↗
DOE OSTI2025

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

High Order Filter Methods for the Non-ideal Compressible MHD Equations

The generalization of a class of low-dissipative high order filter finite difference methods for long time wave propagation of shock/turbulence/combustion compressible viscous gas dynamic flows to compressible MHD equations for structured curvilinear grids has been achieved. The new scheme is shown to provide a natural and efficient way for the minimization of the divergence of the magnetic field numerical error. Standard divergence cleaning is not required by the present filter approach. For certain non-ideal MHD test cases, divergence free preservation of the magnetic fields has been achieved.

Yee, H. C.↗
NASA NTRS1986

A study of the very high order natural user language (with AI capabilities) for the NASA space station common module

The requirements are identified for a very high order natural language to be used by crew members on board the Space Station. The hardware facilities, databases, realtime processes, and software support are discussed. The operations and capabilities that will be required in both normal (routine) and abnormal (nonroutine) situations are evaluated. A structure and syntax for an interface (front-end) language to satisfy the above requirements are recommended.

Gill, E. N.↗
NASA NTRS2002

High-Order Central WENO Schemes for 1D Hamilton-Jacobi Equations

In this paper we derive fully-discrete Central WENO (CWENO) schemes for approximating solutions of one dimensional Hamilton-Jacobi (HJ) equations, which combine our previous works. We introduce third and fifth-order accurate schemes, which are the first central schemes for the HJ equations of order higher than two. The core ingredient is the derivation of our schemes is a high-order CWENO reconstructions in space.

Bryson, Steve↗
DOE OSTI2025

Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods

The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.

Vargas, Arturo [Lawrence Livermore National Labora↗
DOE OSTI2022

High-order ( N = 4–6) multiphoton absorption and mid-infrared Kerr nonlinearity in GaP, ZnSe, GaSe, and ZGP crystals

We report a study of high-order multiphoton absorption, nonlinear refraction, and their anisotropy in four notable mid-infrared χ (2) crystals: GaP, ZnSe, GaSe, and ZGP using femtosecond pulses at 2.35 μm with peak intensity in excess of 200 GW/cm 2 . We found that the nonlinear absorption obeys a perturbation model with multiphoton absorption (MPA) orders from N = 4–6, in agreement with the bandgaps of the crystals. A study of the role of free carrier absorption, performed by changing the pulse duration from 30 to 70 fs while maintaining a constant peak intensity, showed that, at our intensity levels, free carriers generated during the MPA process absorb much more strongly than would be expected from their known linear absorption cross section. Possible mechanisms for this anomalous behavior are high-field effects, such as intravalley scattering in the conduction band and absorption to higher lying bands. Nonlinear refractive indices were measured using (i) closed aperture Z-scan and (ii) spectral broadening due to self-phase modulation, both methods agreeing with each other.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Analysis of single-mode Richtmyer–Meshkov instability using high-order incompressible vorticity-streamfunction and shock-capturing simulations

In this study, two- and three-dimensional simulation results obtained using a new high-order incompressible, variable-density vorticity–streamfunction (VS) method and data from previous ninth-order weighted essentially nonoscillatory (WENO) shock-capturing simulations are used to investigate the nonlinear dynamics of single-mode Richtmyer–Meshkov instability using a model of a Mach 1.3 air(acetone)/SF6 shock tube experiment. A comparison of the density fields from both simulations with the experimental images demonstrates very good agreement in the large-scale structure with both methods but differences in the small-scale structure. The WENO method captures the small-scale disordered structure observed in the experiment, while the VS method partially captures such structure and yields a strong rotating core. The perturbation amplitude growth from the simulations generally agrees well with the experiment. The simulation bubble and spike amplitudes agree well at early times. At later times, the WENO bubble amplitude is smaller than the VS amplitude and vice versa for the spike amplitude. The predictions of nonlinear single-mode instability growth models are shown to agree with the simulation amplitudes at early-to-intermediate times but underpredict the amplitudes at later times in the nonlinear regime. Visualizations of the mass fraction and enstrophy isosurfaces, velocity and vorticity fields, and baroclinic vorticity production and vortex stretching terms from the three-dimensional simulations indicate that, with the exception of the small-scale structure within the rollups, the VS and WENO results are in good agreement.

42 ENGINEERING↗
NASA NTRS1993

Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes

We present a systematic method for constructing boundary conditions (numerical and physical) of the required accuracy, for compact (Pade-like) high-order finite-difference schemes for hyperbolic systems. First, a roper summation-by-parts formula is found for the approximate derivative. A 'simultaneous approximation term' (SAT) is then introduced to treat the boundary conditions. This procedure leads to time-stable schemes even in the system case. An explicit construction of the fourth-order compact case is given. Numerical studies are presented to verify the efficacy of the approach.

Carpenter, Mark H.↗