Search NASA⌕ Search

SEARCH · Search NASA

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

High-Order Stabilized Finite Elements on Dynamic Meshes

The development of dynamic mesh capability for turbulent flow simulations using the Streamlined Upwind Petrov-Galerkin (SUPG) discretization is described. The current work extends previous research to include high-order spatial accuracy, including the satisfaction of the discrete geometric conservation law (GCL) on curved elements. Two closely-related schemes are described and the ability of these schemes to satisfy the GCL, while also maintaining temporal accuracy and conservation is assessed. Studies indicate that although one scheme discretizes the time derivative in conservative form, both schemes exhibit temporal conservation errors that decrease according to the expected design order of accuracy. The source of the temporal conservation errors is examined, and it is demonstrated that many finite-volume and finite-element schemes can also be expected to have difficulty strictly satisfying conservation in time. The effects on conservation are examined and, while present in the simulations, are seen to be negligible for the problems considered.

Anderson, W. Kyle↗
DOE OSTI2023

A high‐order discontinuous Galerkin approach for physics‐based thermospheric modeling

Abstract The accurate prediction of aerodynamic drag on satellites orbiting in the upper atmosphere is critical to the operational success of modern space technologies, such as satellite‐based communication or navigation systems, which have become increasingly popular in the last few years due to the deployment of constellations of satellites in low‐Earth orbit. As a result, physics‐based models of the ionosphere and thermosphere have emerged as a necessary tool for the prediction of atmospheric outputs under highly variable space weather conditions. This paper proposes a high‐fidelity approach for physics‐based space weather modeling based on the solution of the Navier–Stokes equations using a high‐order discontinuous Galerkin method, combined with a matrix‐free strategy suitable for high‐performance computing on GPU architectures. The approach consists of a thermospheric model that describes a chemically frozen neutral atmosphere in nonhydrostatic equilibrium driven by the external excitation of the Sun. A novel set of variables is considered to treat the low densities present in the upper atmosphere and to accommodate the wide range of scales present in the problem. At the same time, and unlike most existing approaches, radial and angular directions are treated in a nonsegregated approach. The study presents a set of numerical examples that demonstrate the accuracy of the approximation and validate the current approach against observational data along a satellite orbit, including estimates of established empirical and physics‐based models of the ionosphere‐thermosphere system. Finally, a one‐dimensional radial derivation of the physics‐based model is presented and utilized for conducting a parametric study of the main thermal quantities under various solar conditions.

Engineering↗
DOE OSTI2024

Switchable planar chirality and spin texture in highly ordered ferroelectric hybrid perovskite domains

Chiral organic-inorganic hybrid perovskites offer a promising platform for developing non-linear chiro-optical applications and chiral-induced spin selectivity. Here, we show that achiral hybrid perovskites that have highly ordered ferroelectric domains with orthogonal polarization exhibit planar chirality, as manifested by second harmonic generation with strong circular dichroism. Interestingly, the handedness of the second harmonic generation circular dichroism response can be alternatingly switched between orthogonally polarized domains and domain walls. To correlate the origin of planar chirality in these domains, atom-resolved atomic force microscopy is used to reveal distinct arrangements of atomic rows in these ferroelectric crystals that indicate planar symmetry breaking. Remarkably, a large SHG anisotropy factor of 0.41 is measured at the domain wall. Additionally, spatially resolved two-photon photoluminescence excitation measurement provide evidence of larger Rashba spin-splitting in these chiral ferroelectric domains compared to domain-free areas. Our discoveries of domain-dependent planar chirality and spin texture hold great promise for advancement of domain-specific chiro-optical applications.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Evaluation of high order schemes for nonlinear wave computations

We present results of the workshop's benchmark problems of the category 2. This category of problems is designed to test the nonlinear wave propagation properties of a computational scheme. We chose three high order spatially accurate algorithms for our computations. These are the Dispersion- Relation- Preserving (DRP) scheme proposed by Tam and his colleagues, a fourth order extension of the MacCormack scheme proposed by Gottlieb and Turkel and an Essentially Non Oscillatory (ENO) scheme proposed by Shu and Osher.

Hayder, M. Ehtesham↗
NASA NTRS1999

A High-Order Method Using Unstructured Grids for the Aeroacoustic Analysis of Realistic Aircraft Configurations

A method for the prediction of acoustic scatter from complex geometries is presented. The discontinuous Galerkin method provides a framework for the development of a high-order method using unstructured grids. The method's compact form contributes to its accuracy and efficiency, and makes the method well suited for distributed memory parallel computing platforms. Mesh refinement studies are presented to validate the expected convergence properties of the method, and to establish the absolute levels of a error one can expect at a given level of resolution. For a two-dimensional shear layer instability wave and for three-dimensional wave propagation, the method is demonstrated to be insensitive to mesh smoothness. Simulations of scatter from a two-dimensional slat configuration and a three-dimensional blended-wing-body demonstrate the capability of the method to efficiently treat realistic geometries.

Atkins, Harold L.↗
NASA NTRS1992

High-order wide-band frequency domain identification using composite curve fitting

A method is presented for curve fitting nonparametric frequency domain data so as to identify a parametric model composed of two models in parallel, where each model has dynamics in a specified portion of the frequency band. This decomposition overcomes the problem of numerical sensitivity since lower order polynomials can be used compared to existing methods which estimate the model as a single entity. Consequently, composite curve fitting is useful for frequency domain identification of high-order systems and/or systems whose dynamics are spread over a large bandwidth. The approach can be extended to identify an arbitrary number of parallel subsystems in specified frequency regimes.

Bayard, D. S.↗
NASA NTRS1998

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗
DOE OSTI2024

Results of High Order Modes Spectra Measurements in 1.3 GHz Cavities for LCLS-II

Fermilab recently completed production and testing of 1.3 GHz cryomodules for the LCLS-II project. Each cryomodule consists of eight TESLA-shaped superconducting elliptical cavities equipped with two High Order Mode (HOM) coupler ports. Measurement of the HOM spectrum is part of the incoming quality control of cavities at room temperature and the final qualification cold test of cryomodules at the Cryomodule Test Facility (CMTF). In this paper we describe the procedure for measuring the HOM spectrum along with further data processing. Finally, we present accumulated statistics of individual HOM frequencies and quality factors related to various cavity vendors and discuss the possible contribution of HOMs to heat loads and beam dynamics.

43 PARTICLE ACCELERATORS↗
NASA NTRS1998

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Abstract Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗
NASA NTRS1991

Forward flight trim calculation and frequency response validation of a high-order helicopter simulation model

This paper describes a new trim procedure, that includes the calculation of the steady-state response of the rotor blades, and that is applicable to straight flight and steady coordinated turns. This paper also describes the results of a validation study for a high order linearized model of helicopter flight dynamics, that includes rotor, inflow, and actuator dynamics. The model is obtained by numerical perturbations of a nonlinear, blade element type mathematical model. Predicted responses are compared with flight test data for two values of flight speed. The comparison is carried out in the frequency domain. Numerical simulations show that the trim algorithm is very accurate, and preserves the periodicity of the aircraft states. The results also indicate that the predictions of the linearized model are in good agreement with flight test data, especially at medium and high frequencies.

Kim, Frederick D.↗
DOE OSTI2022

Thermal feedback coupling in a transient Monte Carlo high-order/low-order scheme

Transient simulations of nuclear systems face the computational challenge of resolving both space and time during reactivity changes. A common strategy for tackling this issue is to split the neutron flux into shape and amplitude functions. This split can be solved with high- order/low-order methods. While this multi-fidelity approach has traditionally been reserved for deterministic methods, it is also possible to implement in Monte Carlo as an efficient alternative to Dynamic Monte Carlo. This work implements the frequency transform method with thermal feedback in high-order/low-order Monte Carlo by blending static coupling methods such as single-batch Monte Carlo, with a simple thermal-fluids calculation. While previous work focused solely on prescribed transients, the addition of time-dependent thermal-fluids allows transients to be self-propagating. Tests were run in a fluids-initiated transient to showcase the basic functionalities of this methodology. Preliminary results behave as expected, paving the way for studying more sophisticated thermal-fluids coupling methods. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

High-order replica bands in monolayer FeSe/SrTiO3 revealed by polarization-dependent photoemission spectroscopy

Abstract The mechanism of the enhanced superconductivity in monolayer FeSe/SrTiO 3 has been enthusiastically studied and debated over the past decade. One specific observation has been taken to be of central importance: the replica bands in the photoemission spectrum. Although suggestive of electron-phonon interaction in the material, the essence of these spectroscopic features remains highly controversial. In this work, we conduct angle-resolved photoemission spectroscopy measurements on monolayer FeSe/SrTiO 3 using linearly polarized photons. This configuration enables unambiguous characterization of the valence electronic structure with a suppression of the spectral background. We consistently observe high-order replica bands derived from various Fe 3 d bands, similar to those observed on bare SrTiO 3 . The intensity of the replica bands is unexpectedly high and different between d xy and d yz bands. Our results provide new insights on the electronic structure of this high-temperature superconductor and the physical origin of the photoemission replica bands.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Solution of the hydrodynamic device model using high-order non-oscillatory shock capturing algorithms

A micron n+ - n - n+ silicon diode is simulated via the hydrodynamic model for carrier transport. The numerical algorithms employed are for the non-steady case, and a limiting process is used to reach steady state. The simulation employs shock capturing algorithms, and indeed shocks, or very rapid transition regimes, are observed in the transient case for the coupled system, consisting of the potential equation and the conservation equations describing charge, momentum, and energy transfer for the electron carriers. These algorithms, termed essentially non-oscillatory, were successfully applied in other contexts to model the flow in gas dynamics, magnetohydrodynamics, and other physical situations involving the conservation laws in fluid mechanics. The method here is first order in time, but the use of small time steps allows for good accuracy. Runge-Kutta methods allow one to achieve higher accuracy in time if desired. The spatial accuracy is of high order in regions of smoothness.

Fatemi, Emad↗
NASA NTRS1991

Solution of the hydrodynamic device model using high-order non-oscillatory shock capturing algorithms

A micron n+ - n - n+ silicon diode is simulated via the hydrodynamic model for carrier transport. The numerical algorithms employed are for the non-steady case, and a limiting process is used to reach steady state. The simulation employs shock capturing algorithms, and indeed shocks, or very rapid transition regimes, are observed in the transient case for the coupled system, consisting of the potential equation and the conservation equations describing charge, momentum, and energy transfer for the electron carriers. These algorithms, termed essentially nonoscillatory, were successfully applied in other contexts to model the flow in gas dynamics, magnetohydrodynamics, and other physical situations involving the conservation laws in fluid mechanics. The method here is first order in time, but the use of small time steps allows for good accuracy. Runge-Kutta methods allow one to achieve higher accuracy in time if desired. The spatial accuracy is of high order in regions of smoothness.

Fatemi, Emad↗
NASA NTRS2001

Technique for Very High Order Nonlinear Simulation and Validation

Finding the sources of sound in large nonlinear fields via direct simulation currently requires excessive computational cost. This paper describes a simple technique for efficiently solving the multidimensional nonlinear Euler equations that significantly reduces this cost and demonstrates a useful approach for validating high order nonlinear methods. Up to 15th order accuracy in space and time methods were compared and it is shown that an algorithm with a fixed design accuracy approaches its maximal utility and then its usefulness exponentially decays unless higher accuracy is used. It is concluded that at least a 7th order method is required to efficiently propagate a harmonic wave using the nonlinear Euler equations to a distance of 5 wavelengths while maintaining an overall error tolerance that is low enough to capture both the mean flow and the acoustics.

Dyson, Rodger W.↗
DOE OSTI2022

High-order diamond differencing schemes for the Boltzmann Fokker-Planck equation in 3D Cartesian geometries

The Boltzmann Fokker-Planck, an approximate form of the linear Boltzmann equation is commonly used to treat efficiently the transport of charged particles in matter. This paper introduces the application of high-order diamond differencing schemes (HODD), specifically the DD1 and DD2 schemes which are 4- and 6-order accurate respectively, to handle the spatial discretization of that equation in 3D Cartesian geometries. The energy deposition solutions for the coupled transport of electrons and photons presented in this work shows that HODD, compared to classical DD scheme, provides correction to the oscillations and a reduced propensity to yield negative fluxes. They are useful tools to minimize local error, notably in regions with abrupt variations of the flux solution. They also can be used to reduced execution time by decreasing the needed number of voxels to obtain a fixed accuracy. On the tested benchmarks, the DD1 scheme is 87%- 92%-91% more accurate than the classical DD scheme for total, mean per-voxels and maximum deviation of energy deposition values respectively. For comparison, a calculation with 8 times more voxels, requiring roughly 2.5 times more time to execute, is 92%-90%-77% more accurate. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

High-order dimensionally-split Cartesian embedded boundary method for non-dissipative schemes

Centered finite-difference schemes are commonly used for high-fidelity turbulent flow simulations in canonical configurations because of their non-dissipative property and computational efficiency. However, their use in flow simulations over complex geometries is limited by the requirements of a structured grid and a stable boundary treatment in the absence of artificial (numerical) dissipation. Cartesian embedded boundary (EB) approaches provide an efficient structured-grid framework to apply difference schemes over complex domains. However, they are often restricted to low orders of accuracy because of numerical instabilities at the embedded boundaries and the issues of small-cell problem that are difficult to address with high-order accuracy. The present work discusses a systematic approach to obtain high-order EB methods with non-dissipative centered schemes in the interior. This approach, based on satisfying the primary and secondary conservation conditions, is employed to derive EB schemes that are up to sixth-order accurate in the interior and fourth-order accurate globally for hyperbolic, parabolic as well as incompletely parabolic problems. The proposed finite-difference discretization is, by construction, dimensionally split and addresses the small-cell problem without any cell/geometry transformations, thus, highly simplifying implementation in a flow solver. Various linear and non-linear numerical tests are performed to evaluate the stability and the accuracy of the proposed EB schemes.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

Nonlinear Analysis of Airfoil High-Intensity Gust Response Using a High-Order Prefactored Compact Code

The nonlinear response of symmetric and loaded airfoils to an impinging vortical gust is investigated in the parametric space of gust dimension, intensity, and frequency. The study, which was designed to investigate the validity limits for a linear analysis, is implemented by applying a nonlinear high-order prefactored compact code and comparing results with linear solutions from the GUST3D frequency-domain solver. Both the unsteady aerodynamic and acoustic gust responses are examined.

Crivellini, A.↗