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

Polarization Dependent Excitation and High Harmonic Generation from Intense Mid-IR Laser Pulses in ZnO

The generation of high order harmonics from femtosecond mid-IR laser pulses in ZnO has shown great potential to reveal new insight into the ultrafast electron dynamics on a few femtosecond timescale. In this work we report on the experimental investigation of photoluminescence and high-order harmonic generation (HHG) in a ZnO single crystal and polycrystalline thin film irradiated with intense femtosecond mid-IR laser pulses. The ellipticity dependence of the HHG process is experimentally studied up to the 17th harmonic order for various driving laser wavelengths in the spectral range 3–4 µm. Interband Zener tunneling is found to exhibit a significant excitation efficiency drop for circularly polarized strong-field pump pulses. For higher harmonics with energies larger than the bandgap, the measured ellipticity dependence can be quantitatively described by numerical simulations based on the density matrix equations. The ellipticity dependence of the below and above ZnO band gap harmonics as a function of the laser wavelength provides an efficient method for distinguishing the dominant HHG mechanism for different harmonic orders.

tunneling excitation↗
DOE OSTI2025

Extreme-ultraviolet spatiotemporal vortices via high harmonic generation

Spatiotemporal optical vortices (STOVs) are space–time structured light pulses with a unique topology that couples spatial and temporal domains and carry transverse orbital angular momentum (OAM). Up to now, their generation has been limited to the visible and infrared regions of the spectrum. During the last decade, it was shown that through the process of high-order harmonic generation, it is possible to upconvert spatial optical vortices that carry longitudinal OAM from the near-infrared into the extreme-ultraviolet (EUV), thereby producing vortices with distinct femtosecond and attosecond structure. In this work, we demonstrate theoretically and experimentally the generation of EUV spatiotemporal and spatiospectral vortices using near-infrared STOV driving laser pulses. Here, we use analytical expressions for focused STOVs to perform macroscopic calculations of high-order harmonic generation that are directly compared to the experimental results. As STOV beams are not eigenmodes of propagation, we characterize the highly charged EUV STOVs in both the near and far fields to show that they represent conjugated spatiotemporal and spatiospectral vortex pairs. Our work provides high-frequency light beams topologically coupled at the nanometre/attosecond scales domains with transverse OAM that could be suitable to explore electronic dynamics in magnetic materials, chiral media and nanostructures.

High-harmonic generation↗
DOE OSTI2021

A new method to compute periodic orbits in general symplectic maps

The search of high-order periodic orbits has been typically restricted to problems with symmetries that help to reduce the dimension of the search space. Well-known examples include reversible maps with symmetry lines. The present work proposes a new method to compute high-order periodic orbits in twist maps without the use of symmetries. The method is a combination of the parameterization method in Fourier space and a Newton–Gauss multiple shooting scheme. The parameterization method has been successfully used in the past to compute quasi-periodic invariant circles. However, this is the first time that this method is used in the context of periodic orbits. Numerical examples are presented showing the accuracy and efficiency of the proposed method. Furthermore, the method is also applied to verify the renormalization prediction of the residues’ convergence at criticality (extensively studied in reversible maps) in the relatively unexplored case of maps without symmetries.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Science & Technology Review: The Road to Exascale Computing

At Lawrence Livermore National Laboratory, we focus on science and technology research to ensure our nation’s security. We also apply that expertise to solve other important national problems in energy, bioscience, and the environment. Science & Technology Review is published eight times a year to communicate, to a broad audience, the Laboratory’s scientific and technological accomplishments in fulfilling its primary missions. The publication’s goal is to help readers understand these accomplishments and appreciate their value to the individual citizen, the nation, and the world. The Department of Energy’s Exascale Computing Project (ECP) and Lawrence Livermore’s RADIUSS (Rapid Application Development via an Institutional Universal Software Stack) initiative benefit from strategically developed software tools. The front cover shows a simulation of advection under twisting rotation that uses high-order finite elements from Livermore’s Modular Finite Element Methods (MFEM) software library and GLVis visualization tool. On the back cover, the logo (also created with GLVis) for the MFEM project illustrates the curved mesh and sub-element resolution used in high-order simulations. MFEM and GLVis are key components of the ECP’s co-design Center for Efficient Exascale Discretizations (CEED) and RADIUSS. MFEM is also part of ECP’s Extreme-Scale Scientific Software Development Kit (xSDK).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Attosecond pulse shaping using a seeded free-electron laser

Attosecond pulses are central to the investigation of valence- and core-electron dynamics on their natural timescales. The reproducible generation and characterization of attosecond waveforms has been demonstrated so far only through the process of high-order harmonic generation. Several methods for shaping attosecond waveforms have been proposed, including the use of metallic filters, multilayer mirrors and manipulation of the driving field. However, none of these approaches allows the flexible manipulation of the temporal characteristics of the attosecond waveforms, and they suffer from the low conversion efficiency of the high-order harmonic generation process. Free-electron lasers, by contrast, deliver femtosecond, extreme-ultraviolet and X-ray pulses with energies ranging from tens of microjoules to a few millijoules. Recent experiments have shown that they can generate subfemtosecond spikes, but with temporal characteristics that change shot-to-shot. In this paper we report reproducible generation of high-energy (microjoule level) attosecond waveforms using a seeded free-electron laser. We demonstrate amplitude and phase manipulation of the harmonic components of an attosecond pulse train in combination with an approach for its temporal reconstruction. The results presented here open the way to performing attosecond time-resolved experiments with free-electron lasers.

47 OTHER INSTRUMENTATION↗
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↗