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
DOE OSTI2023

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Exascale Multiphysics Nuclear Reactor Simulations for Advanced Designs

ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.

Hamilton, Steven P.↗
DOE OSTI2022

A high accuracy/resolution spectral element/Fourier–Galerkin method for the simulation of shoaling non-linear internal waves and turbulence in long domains with variable bathymetry

A high-order hybrid continuous-Galerkin numerical method, designed for the simulation of non-linear, non -hydrostatic internal waves and turbulence in long computational domains with complex bathymetry, is presented. The spatial discretization in the non-periodic wave-propagating directions, utilizes the nodal spectral element method. Such a high-order element-based discretization allows the highly accurate representation of complex domain geometry along with the flexibility of concentrating resolution in areas of interest. Under the assumption of the normal-to-isobath propagation of non-linear internal waves, a third periodic direction is incorporated via a Fourier-Galerkin discretization. The distinct non-hydrostatic nature of non-linear internal waves and, any instabilities and turbulence therein, necessitates the numerically challenging solution of the pressure Poisson problem. A defining feature of this work is the application of a domain decomposition approach, combined with block-Jacobi/deflation-based preconditioning to the pressure Poisson problem. Such a combined approach is particularly suitable for the long high aspect-ratio complex domains of interest and enables the efficient high-accuracy reproduction of the non-hydrostatic dynamics of non-linear internal waves. Implementation details are also described in the context of the stability of the solver and its parallelization strategy. A series of benchmarks of increasing complexity demonstrate the robustness of the flow solver. The benchmarks culminate with the three-dimensional simulation of a convectively breaking mode-one non-linear internal wave over a realistic South-China-Sea bathymetric transect and background current/stratification profiles.

Deflation↗
DOE OSTI2021

Constraining the nonanalytic terms in the isospin-asymmetry expansion of the nuclear equation of state

In this report we examine the properties of the isospin-asymmetry expansion of the nuclear equation of state from chiral two- and three-body forces. We focus on extracting the high-order symmetry energy coefficients that consist of both normal terms (occurring with even powers of the isospin asymmetry) as well as terms involving the logarithm of the isospin asymmetry that are formally nonanalytic around the expansion point of isospin-symmetric nuclear matter. These coefficients are extracted from numerically precise perturbation theory calculations of the equation of state coupled with a new set of finite difference formulas that achieve stability by explicitly removing the effects of higher-order terms in the expansion. We consider contributions to the symmetry energy coefficients from both two- and three-body interactions. It is found that the coefficients of the logarithmic terms are generically larger in magnitude than those of the normal terms from second-order perturbation theory diagrams, but overall the normal terms give larger contributions to the ground state energy. The high-order isospin-asymmetry terms are especially relevant at large densities where they affect the proton fraction in β-equilibrium matter, and in particular we find that at twice saturation density they can reduce the proton fraction by up to 0.02.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Effects of low-frequency voltage on nonlinear standing wave excitation, plasma uniformity, and ion dynamics in dual-frequency asymmetric capacitive discharges

Abstract It is known that in very-high-frequency (VHF) capacitively coupled plasmas, the higher harmonics generated by nonlinear sheath motion can enhance the standing wave effect (SWE), which can lead to center-peaked plasma density profiles. In this work, an improved nonlinear electromagnetic model incorporating a transmission line model, an electron momentum balance model, a bulk plasma model, a collisionless nonlinear numerical sheath model, and an ion Monte-Carlo collision (MCC) model is developed to study the effects of low-frequency (LF) voltage V L on the nonlinear standing wave excitation, plasma uniformity, and ion energy and angular distribution functions (IEDFs and IADFs) in dual-frequency (DF) asymmetric capacitive argon discharges at relatively low pressure of 3 Pa. The plasma diffusion in the radial direction and ion dynamics within the LF oscillating sheath are self-consistently considered. The LF voltage V L at 2 MHz varies from 0 to 700 V while the HF voltage V H at 60 MHz is fixed at 100 V. Simulation results indicate that without the addition of an LF source (i.e. V L = 0 V), there are a considerable number of high-order harmonics with short wavelengths, leading to significant SWE and central peak in the radial plasma density profile. Nevertheless, the high-order harmonic excitations tend to be weakened and merely occur around the phase of the full LF sheath collapse due to a shorter characteristic damping time of the surface waves as V L increases. This, combined with increased surface wavelengths of both the driving frequency and the higher harmonics at a higher V L , leads to suppressed standing waves and improved plasma uniformity. Meanwhile, the simulations show that both the low and the high energy peaks of IEDF move towards higher energies, and the energy peak separation width ΔEbecomes wider with the increase of V L . The IEDF at the radial center of the powered electrode exhibits a broader ΔEthan that at the edge. For the IADF, an increased V L results in more ions incident on the electrode with a smaller deflection angle. Because of a thinner sheath and a higher sheath voltage at the electrode center, the peak value of IADF at the electrode center is greater than that at the edge.

Physics↗
DOE OSTI2020

Printing 2D Conjugated Polymer Monolayers and Their Distinct Electronic Properties

Abstract Recently, 2D monolayer films of conjugated polymers have gained increasing attention owing to the preeminence of 2D inorganic films that exhibit unique optoelectronic and mechanical properties compared to their bulk analogs. Despite numerous efforts, crystallization of semiconducting polymers into highly ordered 2D monolayer films still remains challenging. Herein, a dynamic‐template‐assisted meniscus‐guided coating is utilized to fabricate continuous, highly ordered 2D monolayer films of conjugated polymers over a centimeter scale with enhanced backbone π–π stacking. In contrast, monolayer films printed on solid substrates confer upon the 1D fiber networks strong alkyl side‐chain stacking at the expense of backbone packing. From single‐layers to multilayers, the polymer π‐stacks change from edge‐on to bimodal orientation as the film thickness reaches ≈20 nm. Spectroscopic and cyclic voltammetry analysis reveals an abrupt increase in J‐aggregation and absorption coefficient and a decrease in bandgap and highest occupied molecular orbital level until critical thickness, possibly arising from the straightened polymer backbone. This is corroborated by an abrupt increase in hole mobility with film thickness, reaching a maximum of 0.7 cm 2 V −1 s −1 near the critical thickness. Finally, fabrication of chemical sensors incorporating polymer films of various thicknesses is demonstrated, and an ultrahigh sensitivity of the ≈7 nm thick ultrathin film (bilayers) to 1 ppb ammonia is shown.

Kafle, Prapti↗
DOE OSTI2024

Large‐Area, High‐Numerical‐Aperture, Freeform Metasurfaces

Abstract Nanophotonic devices are optical platforms capable of unprecedented wavefront control. To push the limits of experimental device performance, scalable design methodologies that combine the simplicity and fabricability of conventional design paradigms with the extended capabilities of freeform optimization are required. A novel gradient‐based design framework for large‐area freeform metasurfaces is introduced in which nonlocal interactions between simply shaped nanostructures, placed on an irregular lattice, are tailored to produce high‐order hybridized modes that support customizable large‐angle scattering profiles. Utilizing this approach, multifunctional super‐dispersive metalenses are designed and experimentally demonstrated. The approach to high‐numerical‐aperture radial metalenses capable of diffraction limited focusing and the generation of donut‐shaped point spread functions is also extended. It is anticipated that these concepts will have utility in super‐resolution microscopy, particle trapping, additive manufacturing, and metrology applications that require ultra‐high numerical apertures.

Optics↗
DOE OSTI2020

Beating absorption in solid-state high harmonics

Since the new millennium coherent extreme ultra-violet and soft x-ray radiation has revolutionized the understanding of dynamical physical, chemical and biological systems at the electron’s natural timescale. Unfortunately, coherent laser-based upconversion of infrared photons to vacuum-ultraviolet and soft x-ray high-order harmonics in gaseous, liquid and solid targets is notoriously inefficient. In dense nonlinear media, the limiting factor is strong re-absorption of the generated high-energy photons. Here we overcome this limitation by generating high-order harmonics from a periodic array of thin one-dimensional crystalline silicon ridge waveguides. Adding vacuum gaps between the ridges avoids the high absorption loss of the bulk and results in a ~ 100-fold increase of the extraction depth. As the grating period is varied, each high harmonic shows a different and marked modulation, indicating their waveguiding in the vacuum slots with reduced absorption. Looking ahead, our results enable bright on-chip coherent short-wavelength sources and may extend the usable spectral range of traditional nonlinear crystals to their absorption windows. Potential applications include on-chip chemically-sensitive spectro-nanoscopy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The Paschen Jump as a Diagnostic of the Diffuse Nebular Continuum Emission in Active Galactic Nuclei*

Photoionization modeling of active galactic nuclei (AGN) predicts that diffuse continuum (DC) emission from the broad-line region makes a substantial contribution to the total continuum emission from ultraviolet through near-infrared wavelengths. Evidence for this DC component is present in the strong Balmer jump feature in AGN spectra, and possibly from reverberation measurements that find longer lags than expected from disk emission alone. However, the Balmer jump region contains numerous blended emission features, making it difficult to isolate the DC emission strength. In contrast, the Paschen jump region near 8200 Å is relatively uncontaminated by other strong emission features. Here, we examine whether the Paschen jump can aid in constraining the DC contribution, using Hubble Space Telescope Space Telescope Imaging Spectrograph spectra of six nearby Seyfert 1 nuclei. The spectra appear smooth across the Paschen edge, and we find no evidence of a Paschen spectral break or jump in total flux. We fit multicomponent spectral models over the range 6800–9700 Å and find that the spectra can still be compatible with a significant DC contribution if the DC Paschen jump is offset by an opposite spectral break resulting from blended high-order Paschen emission lines. The fits imply DC contributions ranging from ~10% to 50% at 8000 Å, but the fitting results are highly dependent on assumptions made about other model components. These degeneracies can potentially be alleviated by carrying out fits over a broader wavelength range, provided that models can accurately represent the disk continuum shape, Fe ii emission, high-order Balmer line emission, and other components.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Surface Atomic Rearrangement with High Cation Ordering for Ultra-Stable Single-Crystal Ni-Rich Co-Less Cathode Materials

It is crucial to minimize cobalt content in Ni-rich layered single-crystal cathodes due to their high price and limited availability, yet it will inevitably lead to cation disordering, capacity degradation, and thermal issues. Herein, to overcome the intrinsic trade-off between performance and composition of Ni-rich Co-less single-crystal cathodes, a precursor engineering strategy with an epitaxially grown cobalt enrichment on the surface is innovatively proposed. In contrast to traditional coating modifications with random orientation and rigid surface-bulk boundary, the epitaxially enriched surface cobalt layer on the precursor undergoes rapid interdiffusion with the internal Ni 3+ during the optimized sintering process. This interdiffusion eliminates the surface-bulk boundary, promoting the uniform distribution of cobalt and synergistically addressing the Li/Ni intermixing. Moreover, an enhanced surface Li + diffusion is obtained, thereby suppressing the Li + concentration gradient and intragranular cracks generation. Consequently, the modified LiNi 0.7 Co 0.07 Mn 0.23 O 2 exhibits impressive cycling stability with increased capacity retention in both coin-type half-cells and pouch-type full-cells (91% after 1000 cycles), even under the harsh condition of high-temperature, surpassing the majority of previously reported Ni-rich cathodes. Finally, this work opens new avenues toward the low cost, high energy density, thermal stability, and long cyclic life for Ni-rich Co-less cathodes and sheds light on large-scale commercial production.

25 ENERGY STORAGE↗
DOE OSTI2023

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Grey-Box System Identification of Grid-Forming Inverters

This paper demonstrates the use of grey-box system identification methods for simplifying and understanding the nonlinear power dynamics of grid-forming inverters (GFMs). The power and frequency outputs of complex high-order GFM models are fed into system identification software in order to fit them to a predetermined LTI system and learn system parameters such as (synthetic) inertia and droop constants. The same process is then run for a high-order synchronous generator model, and the outputs are fit to the same set of LTI equations. Simulation of a network of GFM inverters with diverse control architecture is also performed for the same process. The intent is threefold: first, to demonstrate the appropriateness of unified LTI models for describing the power and frequency dynamics of individual resources and connected networks, in order to facilitate analysis of larger heterogeneous networked systems; second, to discover the relationship between internal control parameters of GFMs and their externally observed values; and third, to validate that grey-box data-driven system identification techniques can be a valuable tool to discover the values of important parameters in the absence of explicit vendor models.

analytical models↗
DOE OSTI2023

Extending the predictive power of perturbative QCD using the principle of maximum conformality and the Bayesian analysis

In addition to the evaluation of high-order loop contributions, the precision and predictive power of perturbative QCD (pQCD) predictions depends on two important issues: (1) how to achieve a reliable, convergent fixed-order series, and (2) how to reliably estimate the contributions of unknown higher-order terms. The recursive use of renormalization group equation, together with the Principle of Maximum Conformality (PMC), eliminates the renormalization scheme-and-scale ambiguities of the conventional pQCD series. The result is a conformal, scale-invariant series of finite order which also satisfies all of the principles of the renormalization group. In this paper we propose a novel Bayesian-based approach to estimate the size of the unknown higher order contributions based on an optimized analysis of probability distributions. We show that by using the PMC conformal series, in combination with the Bayesian analysis, one can consistently achieve high degree of reliability estimates for the unknown high order terms. Thus the predictive power of pQCD can be greatly improved. We illustrate this procedure for two pQCD observables: R e+e– and R τ , which are each known up to four loops in pQCD. Numerical analyses confirm that by using the scale-independent and more convergent PMC conformal series, one can achieve reliable Bayesian probability estimates for the unknown higher-order contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Heralded Multiplexed High-Efficiency Cascaded Source of Dual-Rail Entangled Photon Pairs Using Spontaneous Parametric Down-Conversion

Deterministic sources of high-fidelity entangled qubit pairs encoded in the dual-rail photonic basis, i.e., presence of a single photon in one of two orthogonal modes, are a key enabling technology of many applications of quantum information processing, including high-rate, high-fidelity quantum communications over long distances. The most popular and mature sources of such photonic entanglement, e.g., those that leverage spontaneous parametric down-conversion (SPDC) or spontaneous four-wave mixing, generate an entangled (so-called continuous-variable) quantum state that contains contributions from high-order photon terms that lie outside the span of the dual-rail basis, which is detrimental to most applications. One often uses low pump power to mitigate the effects of those high-order terms. However, that reduces the pair generation rate, and the source becomes inherently probabilistic. We investigate a cascaded source that performs a linear-optical entanglement swap between two SPDC sources, to generate a heralded photonic entangled state that has a higher fidelity (to the ideal Bell state) compared to a free-running SPDC source. Furthermore, with the Bell swap providing a heralding trigger, we show how to build a multiplexed source, which despite reasonable switching losses and detector loss and noise, yields a fidelity versus success probability trade-off of a high-efficiency source of high-fidelity dual-rail photonic entanglement. We find, however, that there is a threshold of 1.5 dB of loss per switch, beyond which multiplexing hurts the fidelity versus success probability trade-off.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Multirate Exponential Rosenbrock Methods

In this paper we propose a novel class of methods for high-order accurate integration of multirate systems of ordinary differential equation initial-value problems. The proposed methods construct multirate schemes by approximating the action of matrix φ functions within explicit exponential Rosenbrock (ExpRB) methods, thereby called multirate ExpRB (MERB) methods. They consist of the solution to a sequence of modified “fast” initial-value problems, which may themselves be approximated through subcycling any desired initial-value problem solver. In addition to proving how to construct MERB methods from certain classes of ExpRB methods, we provide rigorous convergence analysis of these methods and derive efficient MERB schemes of orders 2 through 6 (the highest-order infinitesimal multirate methods to date). Lastly, we then present numerical simulations to confirm these theoretical convergence rates and to compare the efficiency of MERB methods against other recently introduced high-order multirate methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Synthesis and characterization of properties of (Tb1/3Mo2/3)2AlC polycrystalline

When molecules, atoms, or ions are grouped in a highly ordered microscopic structure to form a crystal lattice that stretches in all directions, the result is a solid material known as a crystal or crystalline solid. They are arranged in a highly ordered microscopic structure to form a crystal lattice. The smallest group of particles in material that constitutes this repeating pattern is unit cell of the structure. The unit cell is a repeating unit formed by the vectors spanning the points of a lattice.

Bretana, Alex↗
DOE OSTI2020

FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials

This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.

42 ENGINEERING↗
DOE OSTI2022

High harmonic generation in mixed XUV and NIR fields at a free-electron laser

Abstract We present the results of an experiment investigating the generation of high-order harmonics by a femtosecond near-infrared (NIR) laser pulse in the presence of an extreme ultraviolet (XUV) field provided by a free-electron laser (FEL), a process referred to as XUV-assisted high-order harmonic generation (HHG). Our experimental findings show that the XUV field can lead to a small enhancement in the harmonic yield when the XUV and NIR pulses overlap in time, while a strong decrease of the HHG yield and a red shift of the HHG spectrum is observed when the XUV precedes the NIR pulse. The latter observations are in qualitative agreement with model calculations that consider the effect of a decreased number of neutral emitters but are at odds with the predicted effect of the correspondingly increased ionization fraction on the phase matching. Our study demonstrates the technical feasibility of XUV-assisted HHG experiments at FELs, which may provide new avenues to investigate correlation-driven electron dynamics as well as novel ways to study and control propagation effects and phase matching in HHG.

47 OTHER INSTRUMENTATION↗