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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗
DOE OSTI2024

Assessment of Accelerated Stress Testing Data for Silicon Photovoltaics Using Tensor Decomposition Methods

In this work, we examine the use of high-order tensor decompositions to analyze degradation pathways emerging from accelerated stress testing of silicon photovoltaic (PV) modules. Matrix-based decompositions are powerful tools for studying two-dimensional data arrays and form the foundation of a host of classical data analysis techniques. Tensors are high-order extrapolations of matrices that are able to account for more parameter dimensions, and a variety of tensor decomposition methods have been developed that similarly seek to extend insights from matrix decompositions to higher dimensions. Applying and interpreting tensor decomposition methods to sequences of PV module image data, we seek to uncover and isolate different degradation modes occurring from accelerated stress testing procedures. Further, we consider the contributions of different modes to PV module performance degradations.

data analysis↗
DOE OSTI2023

Activation of human STING by a molecular glue-like compound

Stimulator of interferon genes (STING) is a dimeric transmembrane adapter protein that plays a key role in the human innate immune response to infection and has been therapeutically exploited for its antitumor activity. The activation of STING requires its high-order oligomerization, which could be induced by binding of the endogenous ligand, cGAMP, to the cytosolic ligand-binding domain. Here we report the discovery through functional screens of a class of compounds, named NVS-STGs, that activate human STING. Our cryo-EM structures show that NVS-STG2 induces the high-order oligomerization of human STING by binding to a pocket between the transmembrane domains of the neighboring STING dimers, effectively acting as a molecular glue. Our functional assays showed that NVS-STG2 could elicit potent STING-mediated immune responses in cells and antitumor activities in animal models.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

First-passage time statistics on surfaces of general shape: Surface PDE solvers using Generalized Moving Least Squares (GMLS)

Here, we develop numerical methods for computing statistics of stochastic processes on surfaces of general shape with drift-diffusion dynamics d X t = a (X t ) dt + b(X t ) d W t . We formulate descriptions of Brownian motion and general drift-diffusion processes on surfaces. We consider statistics of the form u (x) = E x [$∫^{τ}_{0}$ g (X t ) dt ] + E x [ f (X τ )] for a domain Ω and the exit stopping time τ = inf t { t >0 | X i Ω}, where f , g are general smooth functions. For computing these statistics, we develop high-order Generalized Moving Least Squares (GMLS) solvers for associated surface PDE boundary-value problems based on Backward- Kolmogorov equations. We focus particularly on the mean First Passage Times (FPTs) given by the case f = 0, g = 1 where u (x) = E x [τ]. We perform studies for a variety of shapes showing our methods converge with high-order accuracy both in capturing the geometry and the surface PDE solutions. We then perform studies showing how statistics are influenced by the surface geometry, drift dynamics, and spatially dependent diffusivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Quantum coherence tomography of light-controlled superconductivity

The coupling between superconductors and oscillation cycles of light pulses, i.e., lightwave engineering, is an emerging control concept for superconducting quantum electronics. Although progress has been made towards terahertz-driven superconductivity and supercurrents, the interactions able to drive non-equilibrium pairing are still poorly understood, partially due to the lack of measurements of high-order correlation functions. In particular, the sensing of exotic collective modes that would uniquely characterize light-driven superconducting coherence, in a way analogous to the Meissner effect, is very challenging but much needed. Here we report the discovery of parametrically driven superconductivity by light-induced order-parameter collective oscillations in iron-based superconductors. The time-periodic relative phase dynamics between the coupled electron and hole bands drives the transition to a distinct parametric superconducting state out-of-equalibrium. This light-induced emergent coherence is characterized by a unique phase–amplitude collective mode with Floquet-like sidebands at twice the Higgs frequency. We measure non-perturbative, high-order correlations of this parametrically driven superconductivity by separating the terahertz-frequency multidimensional coherent spectra into pump–probe, Higgs mode and bi-Higgs frequency sideband peaks. We find that the higher-order bi-Higgs sidebands dominate above the critical field, which indicates the breakdown of susceptibility perturbative expansion in this parametric quantum matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part I: the Linear Setting

Fully implicit Runge--Kutta (IRK) methods have many desirable properties as time integration schemes in terms of accuracy and stability, but high-order IRK methods are not commonly used in practice with numerical PDEs due to the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic preconditioning framework for solving the systems of equations that arise from IRK methods applied to linear numerical PDEs (without algebraic constraints). Additionally, this framework also naturally applies to discontinuous Galerkin discretizations in time. Under quite general assumptions on the spatial discretization that yield stable time integration, the preconditioned operator is proven to have condition number bounded by a small, order-one constant, independent of the spatial mesh and time-step size, and with only weak dependence on number of stages/polynomial order; for example, the preconditioned operator for 10th-order Gauss IRK has condition number less than two, independent of the spatial discretization and time step. The new method can be used with arbitrary existing preconditioners for backward Euler-type time-stepping schemes and is amenable to the use of three-term recursion Krylov methods when the underlying spatial discretization is symmetric. The new method is demonstrated to be effective on various high-order finite-difference and finite element discretizations of linear parabolic and hyperbolic problems, demonstrating fast, scalable solution of up to 10th-order accuracy. The new method consistently outperforms existing block preconditioning approaches, and in several cases, the new method can achieve 4th-order accuracy using Gauss integration with roughly half the number of preconditioner applications and wallclock time as required using standard diagonally IRK methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Cleavable quaternary oxychlorides with high magnetic ordering temperatures

Quaternary oxychlorides derived from Ruddlesden–Popper 3d transition metal oxides offer a route to cleavable crystals with bulk antiferromagnetic ordering temperatures reaching at least 550 K. Here, we study the magnetic, optical, and mechanical behavior of Sr 2 FeO 3 Cl, Ca 2 FeO 3 Cl, Ca 3 Fe 2 O 5 Cl 2 , and Sr 3 Fe 2 O 5 Cl 2 . Through optical absorption measurements, we show that these antiferromagnetic semiconductors have optical band gaps of ≈2.1(1) eV. The magnetic ordering symmetries and temperatures were probed by neutron powder diffraction and Mössbauer spectroscopy on polycrystalline samples, demonstrating Néel temperatures (T N ) near room temperature in the single layer Sr 2 FeO 3 Cl (T N ≈ 311 K) and Ca 2 FeO 3 Cl (T N ≈ 360 K), and the double-layer compound Sr 3 Fe 2 O 5 Cl 2 has T N ≈ 545 K. The high-spin moments of Fe 3+ lie within the basal plane and the magnetic structures are compensated within each magnetic layer and characterized by magnetic propagation vectors k = ($\frac{1}{2}$ $\frac{1}{2}$ 0). Magnetization results demonstrate the quasi-2D nature of the magnetism, with a broad maximum in the susceptibility near 2T N for Sr 2 FeO 3 Cl. Scotch tape tests and mechanical exfoliation onto SiO 2 confirm the micaceous nature of these crystals with cleavage down to a single unit cell (two magnetic layers) achieved for Sr 3 Fe 2 O 5 Cl 2 . In conclusion, this paper highlights strong antiferromagnetic interactions, semiconducting band gaps, and cleavability of quaternary Fe-based oxychlorides and motivates future work on crystals and exfoliated flakes of these and related oxyhalide systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sylvester-preconditioned adaptive-rank implicit time integrators for advection-diffusion equations with variable coefficients

Here, we consider the adaptive-rank integration of multi-dimensional time-dependent advection-diffusion partial differential equations (PDEs) with variable coefficients. We employ a standard finite-difference method for spatial discretization coupled with high-order diagonally implicit Runge-Kutta temporal schemes. The discrete equation is a generalized Sylvester equation (GSE), which we solve with a projection-based adaptive-rank algorithm structured around two key strategies: (i) constructing dimension-wise subspaces using a novel atypical extended Krylov strategy, and (ii) efficiently solving the basis coefficient matrix with a preconditioned GMRES solver. The low-rank decomposition is performed in 2D using SVD and with high-order SVD (HOSVD) in 3D to represent the tensor in a compressed Tucker format. For d-dimensional problems (here, d = 2 or 3), the computational complexity and memory storage of the approach are found numerically to scale as and $\mathscr{O}(Nr^2) + \mathscr{O} (r^{d+1})$ and $\mathscr{O}(Nr) + \mathscr{O} (r^{d})$, respectively, with the one-dimensional resolution and the maximal rank during the Krylov iteration (which we find to be largely independent of on our numerical examples). We present numerical examples that illustrate the advertised properties of the algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Spectral Mesh-Free Quadrature for Planar Regions Bounded by Rational Parametric Curves

This article presents spectral, mesh-free, Green’s theorem-based numerical quadrature schemes for integrating functions over planar regions bounded by rational parametric curves. Our algorithm proceeds in two steps: (1) We first find intermediate quadrature rules for line integrals along the region’s boundary curves corresponding to Green’s theorem. (2) We then use a high-order quadrature rule to compute the numerical antiderivative of the integrand along a coordinate axis, which is used to evaluate the Green’s theorem line integral. We present two methods to compute the intermediate quadrature rule. The first is spectrally accurate (it converges faster than any algebraic order with respect to number of quadrature points) and is relatively easy to implement, but has no guarantee of polynomial exactness. The second guarantees exactness for polynomial integrands up to a pre-specified degree with an a priori-known number of quadrature points and retains the convergence properties of the first, but is slightly more complicated. The quadrature schemes have applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order meshes, and initialization of multi-material simulations with rational geometry. We compare the quadrature schemes produced using our method to other methods in the literature and show that they are much more efficient both in terms of number of quadrature points and computational time. We provide an open-source implementation of the algorithm in MATLAB.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Hierarchical Truncations for Many-Body Expansion Potentials

In this work, a new strategy to truncate high-order terms in the many-body expansion (MBE) is proposed. This new approach, which we call a hierarchical many-body expansion (HMBE), is based on a hierarchical partition of the system into multitier fragments and can in principle be applied to any large molecular system. Numerical tests on a series of (H 2 O) 64 structures are presented, demonstrating satisfactory relative energies between the clusters and binding energies of individual clusters compared with full-cluster calculations, with significantly fewer high-order terms computed than conventional MBE. The hierarchical truncation can be augmented by certain many-body terms for fragments at the interface between the partitions (called “Schengen terms”) to further improve accuracy. This work establishes the HMBE scheme as a promising framework to model very large systems (e.g., proteins), which are naturally built on a hierarchical structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A holographic perspective on the axion quality problem

The axion provides a compelling solution to the strong CP problem as well as a candidate for the dark matter of the universe. However, the axion solution relies on the spontaneous breaking of a global U(1)P Q symmetry, which is also explicitly violated by quantum gravity. To preserve the axion solution, gravitational violations of the U(1)P Q symmetry must be suppressed to sufficiently high order. We present a simple, geometric solution of the axion quality problem by modelling the axion with a bulk complex scalar field in a slice of AdS5, where the U(1)P Q symmetry is spontaneously broken in the bulk but explicitly broken on the UV brane. By localising the axion field towards the IR brane, gravitational violations of the PQ symmetry on the UV brane are sufficiently sequestered. This geometric solution is holographically dual to 4D strong dynamics where the global U(1)P Q is an accidental symmetry to sufficiently high order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Energy conserving and well-balanced discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry

ABSTRACT This paper presents high-order Runge–Kutta (RK) discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry. The scheme can preserve a general polytropic equilibrium state and achieve total energy conservation up to machine precision with carefully designed spatial and temporal discretizations. To achieve the well-balanced property, the numerical solutions are decomposed into equilibrium and fluctuation components that are treated differently in the source term approximation. One non-trivial challenge encountered in the procedure is the complexity of the equilibrium state, which is governed by the Lane–Emden equation. For total energy conservation, we present second- and third-order RK time discretization, where different source term approximations are introduced in each stage of the RK method to ensure the conservation of total energy. A carefully designed slope limiter for spherical symmetry is also introduced to eliminate oscillations near discontinuities while maintaining the well-balanced and total-energy-conserving properties. Extensive numerical examples – including a toy model of stellar core collapse with a phenomenological equation of state that results in core bounce and shock formation – are provided to demonstrate the desired properties of the proposed methods, including the well-balanced property, high-order accuracy, shock-capturing capability, and total energy conservation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

MFA-DVR: direct volume rendering of MFA models

3D volume rendering is widely used to reveal insightful intrinsic patterns of volumetric datasets across many domains. However, the complex structures and varying scales of volumetric data can make efficiently generating high-quality volume rendering results a challenging task. Multivariate functional approximation (MFA) is a new data model that addresses some of the critical challenges: high-order evaluation of both value and derivative anywhere in the spatial domain, compact representation for large-scale volumetric data, and uniform representation of both structured and unstructured data. Here, in this paper, we present MFA-DVR, the first direct volume rendering pipeline utilizing the MFA model, for both structured and unstructured volumetric datasets. We demonstrate improved rendering quality using MFA-DVR on both synthetic and real datasets through a comparative study. We show that MFA-DVR not only generates more faithful volume rendering than using local filters but also performs faster on high-order interpolations on structured and unstructured datasets. MFA-DVR is implemented in the existing volume rendering pipeline of the Visualization Toolkit (VTK) to be accessible by the scientific visualization community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Hybridized Discontinuous Galerkin Methods for Computational Fluid Dynamics

Hybridizable Discontinuous Galerkin (HDG) methods hold promise for any applications with significant advection character, including thermal hydraulics in light water reactors and advanced reactor concepts and fluid models of plasmas in magnetic confinement fusion. Its features include natural upwinding, local element conservation, and extensibility to arbitrarily high order accuracy. In the last fiscal year we have implemented HDG in the Multiphysics Object-Oriented Simulation Environment (MOOSE). We developed a first-of-its-kind automatic static condensation system in MOOSE’s underlying finite element library libMesh which can condense out arbitrarily many internal variables. Finally, we developed the first preconditioner for HDG discretizations of the Navier-Stokes equations which shows robust performance across a wide range of problem sizes and Reynolds numbers. This preconditioner yields solution times that are equivalent to the fastest developed for industry standard finite volume methods. Moreover, the arbitrarily high-order nature of HDG makes it a prime candidate for acceleration via graphical processing units (GPUs). We believe these developments will hold significant importance for future DOE Nuclear Energy (NE) and Fusion Energy Science (FES) programs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Directional Variations in Tidal Flow Multifractality and Intermittency

We explored distinct directional variations of multifractal and intermittent characteristics of ebb and flood flow velocities at nodule point, WA, tidal energy site, and complementary inspection on the East River, NY using scaling exponents of the structure function, distribution flatness, detrending moving average (DMA) analysis, multifractal detrended fluctuation analysis (MF-DFA), and high-order spectral moments. Our findings reveal that tidal flow presents higher long-range dependence (LRD) and intermittent levels for the ebb flow sections, whereas the flood flow sections exhibit a higher degree of multifractality and greater sensitivity to larger magnitude of turbulent fluctuations. We demonstrate that long-range dependence predominantly contributes to multifractal behavior in both ebb and flood flows, as evidenced by the significantly reduced multifractal spectrum width for temporally randomly permuted time series. Moreover, spectral kurtosis analysis uncovers a higher intermittent level across all frequency scales for ebb flow sections and reveals a distinct pattern of tidal flow intermittency differing from the monotonically increasing intermittent level observed in wall-bounded and grid turbulence. Finally, we demonstrate that DMA, MF-DFA, and high-order spectral moments provide more comprehensive insights than structure function scaling exponents and PDF flatness methods.

Cheng, Shyuan [Univ. of Illinois at Urbana-Champai↗
DOE OSTI2021

Improved optical standing-wave beam splitters for dilute Bose–Einstein condensates

Bose–Einstein condensate (BEC)-based atom interferometry exploits low temperatures and long coherence lengths to facilitate high-precision measurements. Progress in atom interferometry promises improvements in navigational devices like gyroscopes and accelerometers, as well as applications in fundamental physics such as accurate determination of physical constants. Previous work demonstrates that beam splitters and mirrors for coherent manipulation of dilute BEC momentum in atom interferometers can be implemented with sequences of non-resonant standing-wave light pulses. While previous work focuses on the optimization of the optical pulses’ amplitude and duration to produce high-order momentum states with high fidelity, we explore how varying the shape of the optical pulses affects optimal beam-splitter performance, as well as the effect of pulse shape on the sensitivity of optimized parameters in achieving high fidelity in high-momentum states. Here, in simulations of two-pulse beam splitters utilizing optimized square, triangle, and sinc-squared pulse shapes applied to dilute BECs, we, in some cases, reduce parameter sensitivity by an order of magnitude while maintaining fidelity.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

A Refinement-by-Superposition -Method for (curl)- and (div)-Conforming Discretizations

Here, we present refinement-by-superposition (RBS) hp-refinement infrastructure for computational electromagnetics (CEMs), which permits exponential rates of convergence. In contrast to dominant approaches to hp-refinement for continuous Galerkin methods, which rely on explicit constraint equations, the multilevel strategy presented drastically reduces the implementation complexity. Through the RBS methodology, enforcement of continuity occurs by construction, enabling arbitrary levels of refinement with ease, and without the practical (but not theoretical) limitations of constrained-node refinement. We outline the construction of the RBS hp-method for refinement with H (curl)- and H (div)-conforming finite cells. Numerical simulations for the 2-D finite element method (FEM) solution of the Maxwell eigenvalue problem demonstrate the effectiveness of RBS hp-refinement. As an additional goal of this work, we aim to promote the use of mixed-order (low- and high-order) elements in practical CEM applications.

42 ENGINEERING↗
DOE OSTI2022

In-Situ Nanoscale Focusing of Extreme Ultraviolet Solid-State High Harmonics

Extreme ultraviolet light delivering radiation with a wavelength shorter than approximately 100 nm is now available from solid-state sources. However, despite exceptional progress, efficient focusing of extreme ultraviolet photons to their ultimate diffraction limit remains a formidable challenge because of the precision of the focusing by curved, optical surfaces. Here we integrate coherent short-wavelength high-order harmonics from a MgO crystal, with a high-numerical-aperture, nanostructured, focusing element etched onto the surface of the crystal itself. We focus extreme ultraviolet light, seventh harmonic of an 800-nm laser, with a zone plate of numerical aperture 0.35, down to a waist radius of 150 nm, with 18% focusing efficiency. The estimated intensity approaches 10 7 W/cm 2 . Future developments may demonstrate nanoscale laser ablation and miniaturization of extreme ultraviolet coherent sources on a chip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗