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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2021

GR-Athena++: Puncture Evolutions on Vertex-centered Oct-tree Adaptive Mesh Refinement

Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and multimessenger astronomy call for highly performant and scalable codes on modern massively parallel architectures. We present GR-Athena++, a general-relativistic, high-order, vertex-centered solver that extends the oct-tree, adaptive mesh refinement capabilities of the astrophysical (radiation) magnetohydrodynamics code Athena++. To simulate dynamical spacetimes, GR-Athena++ uses the Z4c evolution scheme of numerical relativity coupled to the moving puncture gauge. We demonstrate stable and accurate binary black hole merger evolutions via extensive convergence testing, cross-code validation, and verification against state-of-the-art effective-one-body waveforms. GR-Athena++ leverages the task-based parallelism paradigm of Athena++ to achieve excellent scalability. We measure strong-scaling efficiencies above 95% for up to ~1.2 × 10 4 CPUs and excellent weak scaling is shown up to ~10 5 CPUs in a production binary black hole setup with adaptive mesh refinement. GR-Athena++ thus allows for the robust simulation of compact binary coalescences and offers a viable path toward numerical relativity at exascale.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗

DOE OSTI2026

Stochastic tensor contraction for quantum chemistry

Many computational methods in ab initio quantum chemistry are formulated in terms of high-order tensor contractions, whose cost determines the size of system that can be studied. We introduce stochastic tensor contraction to perform such operations with greatly reduced cost, and present its application to the gold-standard quantum chemistry method, coupled cluster theory with up to perturbative triples. For total energy errors more stringent than chemical accuracy, we reduce the computational scaling to that of mean-field theory, while starting to approach the mean-field absolute cost, thereby challenging the existing cost-to-accuracy landscape. Benchmarks against state-of-the-art local correlation approximations further show that we achieve an order-of-magnitude improvement in both total computation time and error, with significantly reduced sensitivity to system dimensionality and electron delocalization. We conclude that stochastic tensor contraction is a powerful computational primitive to accelerate a wide range of quantum chemistry.

Chemical Physics (physics.chem-ph)↗

DOE OSTI2026

Disentangling High Harmonic Generation from Surface and Bulk States of a Topological Insulator

The discovery of topological phases has introduced a new dimension to materials science. Three-dimensional (3D) topological insulators (TIs) are a remarkable class of matter that is insulating in the bulk while hosting conductive topological surface states (TSSs) with unique charge and spin properties. High-order harmonic generation (HHG) has emerged as a powerful tool to probe condensed matter systems by providing insights into their electronic structure and dynamic behavior. Here, we investigate HHG in the prototype 3D-TI Bi$_2$Se$_3$. We demonstrate that the contributions of bulk and surface states to the harmonic emission can be controlled by tuning the thickness of thin film samples. An ultrathin (6 nm) film substantially enhances HHG from the surface states, while the bulk states dominate HHG in a thicker (50 nm) film. By applying a quasi-static terahertz perturbing field, we disentangle the bulk and surface responses and reveal the significant impact of the surface states' shift vector and Berry curvature on HHG. Our study provides effective methods for isolating the optical responses of TSSs from those of the bulk, which opens the door to resolving an ongoing debate regarding whether it is possible to reliably extract topological signatures in HHG.

Atomic Physics (physics.atom-ph)↗

DOE OSTI2020

The QCD phase diagram and statistics friendly distributions

The preliminary STAR data for proton cumulants for central collisions at \sqrt{s}=7.7 GeV are consistent with a two-component proton multiplicity distribution. We show that this two-component distribution is statistics friendly in that factorial cumulants of surprisingly high orders may be extracted with a relatively small number of events. As a consequence the two-component model can be tested and verified right now with the presently available STAR data from the first phase of the RHIC beam energy scan.

Koch, Volker↗

DOE OSTI2020

Multigrid Algorithms with Projection and Prolongation over Elements of the Phase Space for K-Eigenvalue Transport Problems

This paper describes new multilevel acceleration methods for solving the multigroup neu- tron transport eigenvalue problems. These multilevel algorithms use different projection and prolongation operators in the phase space. The Nonlinear Diffusion Acceleration (NDA) method with multiple grids in energy is formulated with the prolongation oper- ator based on multiplication iterative correction and linear-in-energy mapping. Another multilevel NDA method uses the projection operator with coarsening in energy between the high-order transport and low-order NDA equations. The third algorithm is formu- lated with the partial-current based CMFD low-order equations and applies projection operators in space and energy. The numerical results are presented.

Cornejo, Luke↗

DOE OSTI2021

High-current conduction cooled superconducting radio-frequency cryomodule

A high-current, compact, conduction cooled superconducting radio-frequency cryomodule for particle accelerators. The cryomodule will accelerate an electron beam of average current up to 1 ampere in continuous wave (CW) mode or at high duty factor. The cryomodule consists of a single-cell superconducting radio-frequency cavity made of high-purity niobium, with an inner coating of Nb3Sn and an outer coating of pure copper. Conduction cooling is achieved by using multiple closed-cycle refrigerators. Power is fed into the cavity by two coaxial couplers. Damping of the high-order modes is achieved by a warm beam-pipe ferrite damper.

Ciovati, Gianluigi↗

DOE OSTI2021

Method to synthesize nanoparticle supercrystals

Highly ordered arrays of 3D faceted nanoparticle supercrystals are formed by self-assembly with controlled nanoparticle packing and unique facet dependent optical property by using a binary solvent diffusion method. The binary diffusion results in supercrystals whose size and quality are determined by initial nanoparticle concentration and diffusion speed. The supercrystal solids display unique facet-dependent surface plasmonic and surface-enhanced Raman characteristics. The supercrystals have potential applications in areas such as optics, electronics, and sensor platforms.

Fan, Hongyou↗

DOE OSTI2021

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗

DOE OSTI2021

FULL RANGE TUNE SCAN STUDIES USING GRAPHICS PROCESSING UNITS WITH CUDA IN EIC BEAM-BEAM SIMULATIONS

The hadron beam in the Electron-Ion Collider (EIC) suffers high order betatron and synchro-betatron resonances. In this paper, we present a weak-strong full range (0.0 ~ 0.5) fractional tune scan with a step size as small as 0.001. Multiple Graphics Processing Units (GPUs) are used to speed up the simulation. A code parallelized with MPI and CUDA is implemented. The good tune region from weak-strong scan is further checked by the self-consistent strong-strong simulation. This study provides beam dynamics guidance in choosing proper working points for the future EIC.

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

thornado-transport: Anderson- and GPU-accelerated nonlinear solvers for neutrino-matter coupling

Algorithms for neutrino-matter coupling in core-collapse supernovae (CCSNe) are investigated in the context of a spectral two-moment model, which is discretized in space with the discontinuous Galerkin method, integrated in time with implicit-explicit (IMEX) methods, and implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado). The model considers electron neutrinos and antineutrinos and tabulated opacities from Bruenn (1985), which includes neutrino-electron scattering and pair processes. The nonlinear system arising from implicit time discretization of the equations governing neutrino-matter coupling is iterated to convergence using Anderson-accelerated fixed-point methods, which avoid formation of Jacobians and inversion of dense linear systems. Numerical experiments show that, for a given tolerance, a nested iteration scheme which aims to reduce opacity evaluations can lower the computational cost. Our initial port to GPUs, using both OpenMP and OpenACC, shows an overall speedup of up to ~ 100× when compared to results using a single CPU core. These results indicate that the algorithms implemented in thornado are well-suited to GPU acceleration.

Laiu, Paul↗

DOE OSTI2022

Improvements to the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. This document summarizes the transport solver development efforts pursued during Fiscal Year 2022. We added the multiphysics transient capability for the coarse-mesh finite difference accelerated Richardson iteration for discontinuous finite element method discrete ordinates (DFEM-SN) scheme to support high-order heterogeneous transport simulations. HFEM (hybrid finite element method) - PN (spherical harmonics expansion) was completed and red-black iteration was added for solving the HFEM-PN system with both preconditioned Jacobian-free Newton Krylov and Richardson iteration solvers. The HFEM-PN scheme, as one of the low-order transport schemes, is expected for supporting routine design simulations. Pin power reconstruction capability was also designed and implemented with the Griffin ISOXML module to enhance all the low-order transport solvers for more accurate multiphysics simulations. Numerical results are presented for demonstrating the capabilities and verifying their performance, and future works are discussed.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Magnetic ordering at anomalously high temperatures in selected lanthanides: what about Pr?

Recent experiments have given evidence that for elemental Nd, Tb, and Dy the application of sufficient pressure pushes their magnetic ordering temperatures to anomalously high values. Pr metal has a dhcp structure and is trivalent with the 4f 2 configuration. Its singlet ground state suppresses magnetic order from the 15 K anticipated from de Gennes scaling to 50 mK. Four-point electrical resistivity measurements were carried out on Pr for temperatures 1.5–295 K under pressures to 48 GPa. Although no clear evidence for magnetic order (or superconductivity) is observed above 1.5 K, the temperature dependence of the resistivity gives evidence that Pr may enter a dense Kondo state above 10 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Magnetic Ordering in a High‐Spin Donor–Acceptor Conjugated Polymer

Abstract The development of open‐shell organic molecules that magnetically order at room temperature,which can be practically applied, remains a grand challenge in chemistry, physics, and materials science. Despite the exploration of vast chemical space, design paradigms for organic paramagnetic centers generally result in unpaired electron spins that are unstable or isotropic. Here, a high‐spin conjugated polymer is demonstrated, which is composed of alternating cyclopentadithiophene and benzo[1,2‐ c ;4,5‐ c ′]bis[1,2,5]thiadiazole heterocycles, in which macromolecular structure and topology coalesce to promote the spin center generation and intermolecular exchange coupling. Electron paramagnetic resonance (EPR) spectroscopy is consistent with spatially localized spins, while magnetic susceptibility measurements show clear anisotropic spin ordering and exchange interactions that persist at room temperature. The application of long‐range π‐correlations for spin center generation promotes remarkable stability. This work offers a fundamentally new approach to the implementation of this long‐sought‐after physical phenomenon within organic materials and the integration of manifold properties within emerging technologies.

Steelman, Michael E.↗

DOE OSTI2021

Designing order–disorder transformation in high-entropy ferritic steels

Order–disorder transformations hold an essential place in chemically complex high-entropy ferritic steels (HEFSs) due to their critical technological application. The chemical inhomogeneity arising from mixing of multi-principal elements of varying chemistry can drive property altering changes at the atomic scale, in particular short-range order. Using density-functional theory-based linear-response theory, we predict the effect of compositional tuning on the order–disorder transformation in ferritic steels—focusing on Cr–Ni–Al–Ti–Fe HEFSs. We show that Ti content in Cr–Ni–Al–Ti–Fe solid solutions can be tuned to modify short-range order that changes the order–disorder path from BCC-B2 (Ti atomic-fraction = 0) to BCC-B2-L2 1 (Ti atomic-fraction > 0) consistent with existing experiments. Our study suggests that tuning degree of SRO through compositional variation can be used as an effective means to optimize phase selection in technologically useful alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Grain Boundary Segregation Suppresses Local Short‐Range Ordering in Nanocrystalline High‐Entropy Alloys

Multi-principal-element alloys like high-entropy alloys (HEAs) have potential applications in many engineering fields due to their unique mechanical/functional properties. While HEAs are generally considered random solid solutions, recent studies revealed that they are prone to short-range-ordering (SRO) due to the complex multi-pair-wise interactions among the constituent elements. Meanwhile, SROs' evolution can sometimes be deleterious, and it is necessary to have control over their evolution. Examining the AlCoCrFe-Zr model alloy, long-range ordering occurs following the expectation of enthalpic predictions. Advanced characterization techniques—transmission electron microscopy, high-energy synchrotron X-ray diffraction/pair distribution function, and atom probe tomography, reveal that SRO is suppressed in as-milled and GB-decorated NC-(AlCoCrFe)100-xZrx (x = 0–1.5 atomic %). Warren-Cowley coefficient calculations are further used to validate the suppression of SRO. Besides the low segregation enthalpies of Cr, Fe, and Zr, and the high-mixing enthalpy of Cr and Fe, the short diffusion path to GBs due to high-GB density in the NC-HEAs and the higher energy state of the GBs than the matrix promotes GB-segregation that further alters the matrix chemistry and consequently disfavors SRO formation within the matrix. Despite the GB-segregation of Cr, Fe, and Zr, the matrices and GBs remain in a random solid solution.

36 MATERIALS SCIENCE↗

DOE OSTI2021

First saturation correction in high energy proton-nucleus collisions. Part I. Time evolution of classical Yang-Mills fields beyond leading order

In high energy proton-nucleus collisions, the single- and double-inclusive soft gluon productions at the leading order have been calculated and phenomenologically studied in various approaches for many years. These studies do not take into account the saturation and multiple rescatterings in the field of the proton. The first saturation correction to these leading order results (the terms that are enhanced by the combination $ {\alpha}_s^2{\mu}^2 $ , where μ 2 is the proton’s color charge squared per unit transverse area) has not been completely derived despite recent attempts using a diagrammatic approach. This paper is the first in a series of papers towards analytically completing the first saturation correction to physical observables in high energy proton-nucleus collisions. Our approach is to analytically solve the classical Yang-Mills equations in the dilute-dense regime using the Color Glass Condensate effective theory and compute physical observables constructed from classical gluon fields. In the current paper, the Yang-Mills equations are solved perturbatively in the field of the dilute object (the proton). Next-to-leading order and next-to-next-to-leading order analytic solutions are explicitly constructed. A systematic way to obtain all higher order analytic solutions is outlined.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

High-Fidelity Large-Signal Order Reduction Approach for Composite Load Model

With the increasing penetration of electronic loads and distributed energy resources, conventional load models cannot capture their dynamics. Therefore, a new comprehensive composite load model is developed by Western Electricity Coordinating Council (WECC). However, this model is a complex high-order non-linear system with multi-time-scale property, which poses challenges on power system studies with heavy computational burden. In order to reduce the model complexity, the authors firstly develop a large-signal order reduction (LSOR) method using singular perturbation theory. In this method, the fast dynamics are integrated into the slow ones to preserve transient characteristics of the former. Then, accuracy assessment conditions are proposed and embedded into the LSOR to improve and guarantee the accuracy of reduced-order model. Finally, the reduced-order WECC composite load model is derived by using the proposed algorithm. Overall, simulation results show that the reduced-order large-signal model significantly alleviates the computational burden while maintaining similar dynamic responses as the original composite load model.

24 POWER TRANSMISSION AND DISTRIBUTION↗