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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High-Order Hybrid RANS-LES Study of NACA0012 Wing Sections

We develop hybrid RANS-LES strategies within Nek5000 for application to airfoil sections at small flight configurations. We present a validation and verification study of $k \space – \space \tau$ SST applied to a NACA 0012 wing section in a pure RANS and in a hybrid RANS-LES setup. The study shows good corroboration with existing experimental and numerical datasets. We also analyze some of the observed discrepancies with the experiments by evaluating the side wall “blocking” effect. We demonstrate that for the hybrid turbulence modeling approach a high-order spectral- element discretization converges faster (i.e., with less resolution) than a representative low-order finite-volume-based approach.

42 ENGINEERING↗
DOE OSTI2024

Highly Ordered Eutectic Mesostructures via Template‐Directed Solidification within Thermally Engineered Templates

Template-directed self-assembly of solidifying eutectics results in emergence of unique microstructures due to diffusion constraints and thermal gradients imposed by the template. Here, the importance of selecting the template material based on its conductivity to control heat transfer between the template and the solidifying eutectic, and thus the thermal gradients near the solidification front, is demonstrated. Simulations elucidate the relationship between the thermal properties of the eutectic and template and the resultant microstructure. The overarching finding is that templates with low thermal conductivities are generally advantageous for forming highly organized microstructures. When electrochemically porosified silicon pillars (thermal conductivity < 0.3 Wm −1 K −1 ) are used as the template into which an AgCl-KCl eutectic is solidified, 99% of the unit cells in the solidified structure exhibit the same pattern. In contrast, when higher thermal conductivity crystalline silicon pillars (≈100 Wm −1 K −1 ) are utilized, the expected pattern is only present in 50% of the unit cells. The thermally engineered template results in mesostructures with tunable optical properties and reflectances nearly identical to the simulated reflectances of perfect structures, indicating highly ordered patterns are formed over large areas. This work highlights the importance of controlling heat flows in template-directed self-assembly of eutectics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Scalable generalized screening for high-order terms in the many-body expansion: Algorithm, open-source implementation, and demonstration

The many-body expansion lies at the heart of numerous fragment-based methods that are intended to sidestep the nonlinear scaling of ab initio quantum chemistry, making electronic structure calculations feasible in large systems. In principle, inclusion of higher-order n-body terms ought to improve the accuracy in a controllable way, but unfavorable combinatorics often defeats this in practice and applications with n ≥ 4 are rare. Here, we outline an algorithm to overcome this combinatorial bottleneck, based on a bottom-up approach to energy-based screening. This is implemented within a new open-source software application (“Fragme∩t”), which is integrated with a lightweight semi-empirical method that is used to cull subsystems, attenuating the combinatorial growth of higher-order terms in the graph that is used to manage the calculations. This facilitates applications of unprecedented size, and we report four-body calculations in (H2O)64 clusters that afford relative energies within 0.1 kcal/mol/monomer of the supersystem result using less than 10% of the unique subsystems. We also report n-body calculations in (H2O)20 clusters up to n = 8, at which point the expansion terminates naturally due to screening. These are the largest n-body calculations reported to date using ab initio electronic structure theory, and they confirm that high-order n-body terms are mostly artifacts of basis-set superposition error.

Chemistry↗
DOE OSTI2020

A novel sector mass spectrograph design for high-order coded aperture Mass Spectrometry with stigmatic aberration correction

Previous works experimentally demonstrated that the application of spatially-coded apertures to simple sector mass spectrographs can result in increased signal intensity while maintaining spectral resolution comparable to that of a conventional slit aperture.[Chen, E.X., et al.: J Am Soc Mass Spectr. 26, 1633 (2015)][Russell, Z.E., et al.: J Am Soc Mass Spectr. 26, 248 (2015)] This spatial coding technique was later applied to a double-focusing Mattauch-Herzog mass spectrograph, and shown to be generally compatible with the architecture.[Russell, Z.E., et al.: J Am Soc Mass Spectr. 27, 1 (2016)] However, the achieved signal gain was limited by the spatial size of the electric sector gap. In this work, we introduce a modified version of the Mattauch-Herzog mass spectrograph with a wide-gap electric sector that is enabled by a novel, segmented-electrode approach, and show via simulation that this modified architecture alone (without aperture coding) achieves a 1.8x improvement in mass resolving power compared to a conventional Mattauch-Herzog spectrograph of equivalent size. Further, simulation of an aperture-coded version of the instrument supports a predicted signal strength increase of 50x with minimal loss of aperture clarity on the sensor due to the increased open area of the high-order coded aperture which can be passed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Residual estimation for grid modification in wall-modeled large eddy simulation using unstructured high-order methods

Here, the accuracy and computational cost of a large eddy simulation are highly dependent on the computational grid. Building optimal grids manually from a priori knowledge is not feasible in most practical use cases; instead, solution-adaptive strategies can provide a robust and cost-efficient method to generate a grid with the desired accuracy. We adapt the residual estimation algorithm developed by Toosi and Larsson for Discontinuous Galerkin Spectral Elements Methods (DGSEM) to guide the grid-adaptation process. The core of the method is the computation of the estimated modeling residual using the polynomial basis functions used in DGSEM, and the averaging of the estimated residual over each element. The final method is assessed in multiple channel flow test cases and for the transonic flow over an airfoil, in both cases making use of mortar interfaces between elements with hanging nodes. The method is found to be robust and reliable, and to provide solutions on grids with significantly fewer elements at comparable accuracy compared to when using human-generated grids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A high order Cartesian grid, finite volume method for elliptic interface problems

We present a higher-order finite volume method for solving elliptic PDEs with jump conditions on interfaces embedded in a 2D Cartesian grid. Second, fourth, and sixth order accuracy is demonstrated on a variety of tests including problems with high-contrast and spatially varying coefficients, large discontinuities in the source term, and complex interface geometries. We include a generalized truncation error analysis based on cell-centered Taylor series expansions, which then define stencils in terms of local discrete solution data and geometric information. In the process, we develop a simple method based on Green's theorem for computing exact geometric moments directly from an implicit function definition of the embedded interface. This approach produces stencils with a simple bilinear representation, where spatially-varying coefficients and jump conditions can be easily included and finite volume conservation can be enforced.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Novel Approach to Fast Calculation of High-Order Q-Cumulants

The method of Q-cumulants has been shown as a powerful tool to study the fine details of the azimuthal anisotropies in high-energy nucleus-nucleus collisions. A new method for the fast calculation of arbitrary order Q-cumulant $v_{n}\{2k\}$ values, based on the partition of a non-negative integer l $\le$ m for calculation of the 2m-particle azimuthal correlations is presented in this paper. Unlike the standard Q-cumulants method in which the calculation of high-order multi-particle calculations is impractical, the newly proposed method enables easy calculation. The validity of the method is proven via a toy model that uses the elliptic power distribution to simulate anisotropic emission of particles. The method enables the study of fine details of the $v_{2}$ distribution, such as higher-order central moments of the $v_{2}$ distribution, as well as the hydrodynamic behavior of the Quark-Gluon Plasma.

FOS: Physical sciences↗
DOE OSTI2023

An Efficient High-Order Solver for Diffusion Equations with Strong Anisotropy on Non-Anisotropy-Aligned Meshes

This paper concerns numerical solution of the diffusion equation with strong anisotropy on meshes not aligned with the anisotropic vector field. In order to resolve the numerical pollution for simulations on a non-anisotropy-aligned mesh and reduce the associated high computational cost we propose an effective preconditioner, extending our previous work. Similar to the anisotropy-aligned mesh case, we apply the auxiliary space preconditioning framework to design a preconditioner where a continuous finite element space is used as the auxiliary space for the discontinuous finite element space. The key component is an effective line smoother that can mitigate the high-frequency errors perpendicular to the magnetic field. We design a graph-based approach to find such a line smoother that is approximately perpendicular to the vector fields when the mesh does not align with the anisotropy. Finally, numerical experiments for several benchmark problems are presented, demonstrating the effectiveness and robustness of the proposed preconditioner when applied to Krylov iterative methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Extended quantitative rescattering model for simulating high-order harmonic streaking spectra by synchronization of an intense IR laser and a time-delayed attosecond XUV pulse

In this work, we theoretically investigate the modulated high-harmonic generation (HHG) driven by an intense few-cycle infrared (IR) laser field and a weak extreme-ultraviolet (XUV) pulse at a delayed time. We establish an extended quantitative rescattering (EQRS) model to simulate the HHG streaking spectra, with the ideas of correcting the IR ionization and the transition from the ground to continuum states in the strong-field approximation. The EQRS model has an accuracy comparable to that from “exactly” solving the time-dependent Schrödinger equation (TDSE). We reveal that the fringes in the streaking spectra are caused by the interference between the attosecond XUV pulse and harmonics resulting from different recombination pathways under the intense IR laser. We then demonstrate that the XUV pulse can be accurately retrieved by treating the single-atom TDSE results or macroscopic propagation results as the “input” data. This work provides with a tool for efficiently simulating and thoroughly analyzing the XUV-assisted HHG, which could also enhance its capability for tracing the electron dynamics involved in the strong-field phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Efficient Parallel Sparse Symmetric Tucker Decomposition for High-Order Tensors

Tensor based methods are receiving renewed attention in recent years due to their prevalence in diverse real-world applications. There is considerable literature on tensor representations and algorithms for tensor decompositions, both for dense and sparse tensors. Many applications in hypergraph analytics, machine learning, psychometry, and signal processing result in tensors that are both sparse and symmetric, making it an important class for further study. Similar to the critical Tensor Times Matrix chain operation (TTMc) in general sparse tensors, the Sparse Symmetric Tensor Times Same Matrix chain (S3TTMc) operation is compute and memory intensive due to high tensor order and the associated factorial explosion in the number of non-zeros. In this work, we present a novel compressed storage format CSS for sparse symmetric tensors, along with an efficient parallel algorithm for the S3TTMc operation. We theoretically establish that S3TTMc on CSS achieves a better memory versus run-time trade-off compared to state-of-the-art implementations. We demonstrate experimental findings that confirm these results and achieve up to 2.9× speedup on synthetic and real datasets.

Shivakumar, Shruti↗
DOE OSTI2026

Polynomial range estimation as a troubled-cell indicator for high-order methods

Two troubled-cell indicators based on polynomial range estimation methods are used to flag cells that may violate positivity constraints. One method uses interval extension, and the second uses the range enclosure property of the Bernstein polynomial basis. Furthermore, both methods reduce compute time for the positivity preserver by limiting its application to a subset of cells. The Bernstein polynomial method remains effective as the problem dimensionality increases. Interval extension applied to the internal energy equation permits the use of the troubled-cell indicators for rational functions, though performance suffers compared to directly applying the indicators to polynomial functions.

42 ENGINEERING↗
DOE OSTI2020

Two-color-driven enhanced high-order harmonic generation in solids

We theoretically investigate the emission of high-harmonic (HH) radiation in model crystals by bichromatic few-cycle driving pulses that are composed as the phase-coherent superposition of a mid-infrared fundamental pulse and its second harmonic. Adjusting the model-crystal parameters to reproduce the lowest band gap of MgO, we examine the extent to which distinct domains of the HH spectrum can be controlled and enhanced by tuning the temporal profile of the bichromatic driving laser electric field. We change the driving-pulse shape by varying its fundamental-versus-second-harmonic pulse amplitude ratio and delay, while keeping the energy of the driving laser pulse fixed. For suitable amplitude ratios and delays, we find an up to fivefold enhancement of the spectral HH yield and significant shifts of the HH cutoff frequency.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Foundations for high-order, conservative cut-cell methods: Stable discretizations on degenerate meshes

Cut-cell methods for unsteady flow problems can greatly simplify the grid generation process and allow for high-fidelity simulations on complex geometries. However, cut-cell methods have been limited to low orders of accuracy. This is driven, largely, by the variety of procedures typically introduced to evaluate derivatives in a stable manner near the highly irregular embedded geometry. Here, a completely new approach, termed TEMO (truncation error matching and optimization), is taken to solve this problem. The approach is based on two simple and intuitive design principles. These principles directly allow for the construction of stable 8th To the best of the authors' knowledge, these are the highest orders ever achieved for a cut-cell discretization by a significant margin. This is done for both explicit and compact finite differences and is accomplished without any geometric transformations or artificial stabilization procedures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗