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

Influence of Ordered Mesoporous Oxides in Plasma-Assisted Ammonia Synthesis

Widespread implementation of dielectric barrier discharge (DBD)-assisted NH 3 synthesis, a nascent technology operating under sustainable, ambient conditions, is hindered by low energy yields due to, in part, poor fundamental understanding. Porous oxides used to support metal nanoparticle catalysts have shown significant energy yield contributions for DBD-assisted NH 3 synthesis even without metal. Using an AC-powered, coaxial, single-stage reactor at 16 kV with equimolar (N 2 /H 2 ) feed, we measured NH 3 synthesis rates in the presence of different nonordered oxides, ordered SiO 2 structures (SBA-15 and MCM-41), and ordered Al-incorporated analogues (γ-Al 2 O 3 -coated with varying Al-loadings and Al-substitution, respectively: Al 2 O 3 -SBA-15 and Al-MCM-41). We systematically quantified NH 3 energy yield dependence on pore structures and material identities (i.e., ordered pores and Al incorporation) known to facilitate higher DBD-assisted NH 3 synthesis rates. SBA-15 displayed a higher steady-state energy yield than MCM-41, indicating that framework type is a crucial factor, with both ordered porous systems outperforming fumed SiO 2 . 10 wt % Al maximized in situ NH 3 uptake among the various Al loadings, exhibiting a higher steady-state energy yield and similar power to SBA-15. However, Al-MCM-41 had a similar steady-state energy yield and lower power than MCM-41, likely due to the extended γ-Al 2 O 3 surface that has a dielectric constant higher than that of SiO 2 . Both Al-incorporated analogues benefit from surface acid sites that can adsorb NH 3 in situ, resulting in higher overall NH 3 energy yields than that of their parent ordered SiO 2 . Al 2 O 3 -SBA-15 shielded more NH 3 than Al-MCM-41, likely due to a higher acid site density than the acid site identity. Furthermore, Al incorporation via γ-Al 2 O 3 coating more successfully improves the NH 3 energy yield; together with the high-performing ordered framework, these analogues are potential metal catalyst supports with promising energy yields for DBD-assisted synthesis of NH 3 and other chemicals.

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

The virtual element method for linear elastodynamics models: Design, analysis, and implementation

We design the conforming virtual element method for the numerical simulation of two dimensional time-dependent elastodynamics problems. We investigate the performance of the method both theoretically and numerically. We prove the stability and the convergence of the semi-discrete approximation in the energy norm and derive optimal error estimates. We also show the convergence in the L 2 norm. The performance of the virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally seen in the p-refinement setting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Fourth-Order Embedded Boundary Finite Volume Method for the Unsteady Stokes Equations with Complex Geometries

A fourth-order finite volume embedded boundary (EB) method is presented for the unsteady Stokes equations. The algorithm represents complex geometries on a Cartesian grid using EB, employing a technique to mitigate the ``small cut-cell"" problem without mesh modifications, cell merging, or state redistribution. Spatial discretizations are based on a weighted least-squares technique that has been extended to fourth-order operators and boundary conditions, including an approximate projection to enforce the divergence-free constraint. Solutions are advanced in time using a fourth-order additive implicit-explicit Runge-Kutta method, with the viscous and source terms treated implicitly and explicitly, respectively. Formal accuracy of the method is demonstrated with several grid convergence studies, and results are shown for an application with a complex bio-inspired material. In conclusion, the developed method achieves fourth-order accuracy and is stable despite the pervasive small cells arising from complex geometries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

High Resolution Polar Kerr Effect Studies of CsV 3 Sb 5 : Tests for Time-Reversal Symmetry Breaking below the Charge-Order Transition

We report high resolution polar Kerr effect measurements on CsV 3 Sb 5 single crystals in search for signatures of spontaneous time reversal symmetry breaking below the charge order transition at T * = 94 K. Utilizing two different versions of zero-area loop Sagnac interferometers operating at 1550 nm wavelength, each with the fundamental attribute that without a time reversal symmetry breaking sample at its path, the interferometer is perfectly reciprocal, we find no observable Kerr effect to within the noise floor limit of the apparatus at 30 nanoradians. Simultaneous coherent reflection ratio measurements confirm the sharpness of the charge order transition in the same optical volume as the Kerr measurements. At finite magnetic field we observe a sharp onset of a diamagnetic shift in the Kerr signal at T *, which persists down to the lowest temperature without change in trend. Since 1550 nm is an energy that was shown to capture all features of the optical properties of the material that interact with the charge order transition, we are led to conclude that it is highly unlikely that time reversal symmetry is broken in the charge ordered state in CsV 3 Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The influence of chemical short-range order on the nanoindentation properties of high-entropy alloys prepared via laser powder bed fusion

Chemical short-range order (CSRO) plays an instrumental role in determining the mechanical properties of high-entropy alloys (HEAs). Current methods for controlling CSRO mainly focus on cast HEAs, which can not solve the problem of grain coarsening caused by homogenization treatments. Laser Powder Bed Fusion (LPBF), which has extremely high cooling rates, can suppress compositional segregation, refine grain structures, and achieve low CSRO. This offers a promising avenue for quantitatively controlling CSRO without the need for homogenization treatments. In this study, we investigate the CSRO of HEAs fabricated via LPBF and elucidate their influence on nanoindentation behavior. The results show that, at the low cooling rate, the average size of CSRO is 0.66 nm and they occupy 2.8% of the cross-sectional area, while at high cooling rate, the corresponding values are 0.92 nm and 9.9%, respectively, both of which are significantly smaller than those of cast HEAs. Furthermore, by analyzing samples with different CSROs, we observe that the hardness increases and then decreases with the increase of CSRO, primarily due to the preferential formation of dislocations in regions enriched with Cr-Mn-Ni. With the improvement of CSRO, both the size and quantity of Cr-Mn-Ni rich region gradually increased; subsequently, as multiple Cr-Mn-Ni rich regions tend to merge into one during the CSRO improvement, the size of the Cr-Mn-Ni rich region continued to increase while the quantity gradually decreased; when the size of the Cr-Mn-Ni rich region is moderate and widely dispersed, the dislocation density is higher with much more uniform distributions, which favorably improves nanoindentation hardness.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An adaptive discontinuous Petrov-Galerkin method for the Grad-Shafranov equation

In this work, we propose and develop an arbitrary-order adaptive discontinuous Petrov--Galerkin (DPG) method for the nonlinear Grad--Shafranov equation. An ultraweak formulation of the DPG scheme for the equation is given based on a minimal residual method. The DPG scheme has the advantage of providing more accurate gradients compared to conventional finite element methods, which is desired for numerical solutions to the Grad--Shafranov equation. The numerical scheme is augmented with an adaptive mesh refinement approach, and a criterion based on the residual norm in the minimal residual method is developed to achieve dynamic refinement. Nonlinear solvers for the resulting system are explored and a Picard iteration with Anderson acceleration is found to be efficient to solve the system. Finally, the proposed algorithm is implemented in parallel on MFEM using a domain-decomposition approach, and our implementation is general, supporting arbitrary order of accuracy and general meshes. Furthermore, numerical results are presented to demonstrate the efficiency and accuracy of the proposed algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

High-Fidelity CFD Assessments of Flow Resistance in a 61-Pin Wire-Wrapped Assembly with Partially Blocked Channels

The examination of thermal-hydraulic behaviors in wire-wrapped rod bundles continues to be an active area of research. The sodium fast reactor, a prominent candidate in next-generation nuclear designs, utilizes a hexagonal configuration of wire-wrapped fuel pins. Here, the potential for channel blockage within this compact arrangement poses a significant safety challenge, spurring a number of recent experimental and computational investigations to evaluate its impact on flow and heat transfer. The present work aims to benchmark the high-fidelity NekRS computational fluid dynamics (CFD) solver in predicting the pressure drops associated with substantial blockages, using available experimental data as a reference. A 61-pin wire-wrapped fuel assembly with two flow blockage configurations has been simulated and investigated at a range of low to moderate Reynolds numbers (487 ≤ Re ≤ 14 600). The NekRS solver demonstrates an exponential reduction of spatial discretization error with increasing polynomial order. The high level of agreement between the numerical results and measured data confirms the accuracy and consistency of the present numerical approach. This benchmark study establishes the capability of NekRS to perform reliable hydrodynamic simulations for sodium fast reactor applications and supports its use in design, licensing, and safety analyses.

CFD Benchmarking↗
DOE OSTI2021

Novel niobium-doped titanium oxide towards electrochemical destruction of forever chemicals

Electrochemical advanced oxidative processes (EAOP) are a promising route to destroy recalcitrant organic contaminants such as per- and polyfluoroalkyl substances (PFAS) in drinking water. Central to EAOP are catalysis-induced reactive free radicals for breaking the carbon fluorine bonds in PFAS. Generating these reactive species electrochemically at electrodes provides an advantage over other oxidation processes that rely on chemicals or other harsh conditions. Herein, we report on the performance of niobium (Nb) doped rutile titanium oxide (TiO 2 ) as a novel EAOP catalytic material, combining theoretical modeling with experimental synthesis and characterization. Calculations based on density functional theory are used to predict the overpotential for oxygen evolution at these candidate electrodes, which must be high in order to oxidize PFAS. The results indicate a non-monotonic trend in which Nb doping below 6.25 at.% is expected to reduce performance relative to TiO 2 , while higher concentrations up to 12.5 at.% lead to increased performance, approaching that of state-of-the-art Magnéli Ti 4 O 7 . TiO 2 samples were synthesized with Nb doping concentration at 10 at.%, heat treated at temperatures from 800 to 1100 °C, and found to exhibit high oxidative stability and high generation of reactive oxygen free radical species. The capability of Nb-doped TiO 2 to destroy two common species of PFAS in challenge water was tested, and moderate reduction by ~ 30% was observed, comparable to that of Ti 4 O 7 using a simple three-electrode configuration. We conclude that Nb-doped TiO 2 is a promising alternative EAOP catalytic material with increased activity towards generating reactive oxygen species and warrants further development for electrochemically destroying PFAS contaminants.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Interaction-induced topological phase transition and Majorana edge states in low-dimensional orbital-selective Mott insulators

Topological phases of matter are among the most intriguing research directions in Condensed Matter Physics. It is known that superconductivity induced on a topological insulator’s surface can lead to exotic Majorana modes, the main ingredient of many proposed quantum computation schemes. In this context, the iron-based high critical temperature superconductors are a promising platform to host such an exotic phenomenon in real condensed-matter compounds. The Coulomb interaction is commonly believed to be vital for the magnetic and superconducting properties of these systems. This work bridges these two perspectives and shows that the Coulomb interaction can also drive a canonical superconductor with orbital degrees of freedom into the topological state. Namely, we show that above a critical value of the Hubbard interaction the system simultaneously develops spiral spin order, a highly unusual triplet amplitude in superconductivity, and, remarkably, Majorana fermions at the edges of the system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Degradation of Ta 2 O 5 / SiO 2 dielectric cavity mirrors in ultra-high vacuum

In order for optical cavities to enable strong light-matter interactions for quantum metrology, networking, and scalability in quantum computing systems, their mirrors must have minimal losses. However, high-finesse dielectric cavity mirrors can degrade in ultra-high vacuum (UHV), increasing the challenges of upgrading to cavity-coupled quantum systems. We observe the optical degradation of high-finesse dielectric optical cavity mirrors after high-temperature UHV bake in the form of a substantial increase in surface roughness. We provide an explanation of the degradation through atomic force microscopy (AFM), X-ray fluorescence (XRF), selective wet etching, and optical measurements. We find the degradation is explained by oxygen reduction in Ta 2 O 5 followed by growth of tantalum sub-oxide defects with height to width aspect ratios near ten. We discuss the dependence of mirror loss on surface roughness and finally give recommendations to avoid degradation to allow for quick adoption of cavity-coupled systems.

Optics↗
DOE OSTI2024

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Second-order invariant domain preserving approximation of the compressible Navier–Stokes equations

Here, we present a fully discrete approximation technique for the compressible Navier–Stokes equations that is second-order accurate in time and space, semi-implicit, and guaranteed to be invariant domain preserving. The restriction on the time step is the standard hyperbolic CFL condition, i.e. τ ≲ O ( h ) / V where V is some reference velocity scale and h the typical meshsize.

42 ENGINEERING↗
DOE OSTI2020

The arbitrary-order virtual element method for linear elastodynamics models. Convergence, stability and dispersion-dissipation analysis.

We design the conforming virtual element method for the numerical approximation of the two dimensional elastodynamics problem. We prove stability and convergence of the semi-discrete approximation and derive optimal error estimates under $\textit{h}$-refinement in both the energy and the $L^2$ norms, and optimal error estimates under $\textit{p}$-refinement in the energy norm. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally observed under p-refinement. Finally, we present a dispersion-dissipation analysis for both the semi-discrete and fully-discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion-dissipation properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

X-Ray Performance of Critical-angle Transmission Grating Prototypes for the Arcus Mission

Arcus is a proposed Explorer Class soft X-ray grating spectrometer. It aims to explore cosmic feedback by mapping hot gases within and between galaxies and galaxy clusters and characterizing jets and winds from supermassive black holes and to investigate the dynamics of protoplanetary disks and stellar accretion. Arcus features 12 m focal-length grazing-incidence silicon pore optics (SPO) developed for the Athena mission. Critical-angle transmission (CAT) gratings efficiently disperse high diffraction orders onto CCDs. We report new and improved X-ray performance results for Arcus-like CAT gratings, including a record resolving power for two coaligned CAT gratings. Multiple Arcus prototype grating facets were illuminated by an SPO at the PANTER facility. The facets consist of 32 × 32.5 mm 2 patterned silicon membranes, bonded to metal frames. The bonding angle is adjusted according to the measured average tilt angle of the grating bars in the membrane. Two simultaneously illuminated facets show a minor broadening of the Al-K α doublet in the 18th and 21st orders with the best-fit record effective resolving power of ${R}_{G}\approx {1.3}_{-0.5}^{+\infty }\times {10}^{4}$ (3σ), about three to four times the Arcus requirement. We measured the diffraction efficiency of quasi-fully illuminated gratings at O-K wavelengths in orders 4–7 in an Arcus-like configuration and compare results with synchrotron spot measurements. After corrections for geometrical effects and bremsstrahlung continuum, we find agreement between full and spot illumination at the two different facilities, as well as with the models used for Arcus effective area predictions. We find that these flight-like gratings meet the diffraction efficiency and greatly exceed the resolving power Arcus requires.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Coupling Approaches with Non-matching Grids for Classical Linear Elasticity and Bond-based Peridynamic Models in 1D

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. To facilitate the coupling of the two models, non-matching grids are often desirable as nonlocal grids usually require a finer resolution than local grids. In that case, it is often convenient to resort to interpolation operators so that models can exchange information in the overlap regions when nodes from the two grids do not coincide. This paper studies three existing coupling approaches, namely 1) a method that enforces matching displacements in an overlap region, 2) a variant that enforces a constraint on the stresses instead, and 3) a method that considers a variable horizon in the vicinity of the interfaces. Further, the effect of the interpolation order and of the grid ratio on the performance of the three coupling methods with non-matching grids is carefully studied on one-dimensional examples using polynomial manufactured solutions. The numerical results show that the degree of the interpolants should be chosen with care to avoid introducing additional modeling errors, or simply minimize these errors, in the coupling approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A space-time parallel algorithm with adaptive mesh refinement for computational fluid dynamics

This work describes a space-time parallel algorithm with space-time adaptive mesh refinement (AMR). AMR with subcycling is added to multigrid reduction-in-time (MGRIT) in order to provide solution efficient adaptive grids with a reduction in work performed on coarser grids. This algorithm is achieved by integrating two software libraries: XBraid (Parallel time integration with multigrid. https://computation.llnl.gov/projects/parallel-timeintegration-multigrid) and Chombo (Chombo software package for AMR applications—design document, 2014). The former is a parallel time integration library using multigrid and the latter is a massively parallel structured AMR library. Employing this adaptive space-time parallel algorithm is Chord (Comput Fluids 123:202–217, 2015), a computational fluid dynamics (CFD) application code for solving compressible fluid dynamics problems. For the same solution accuracy, speedups are demonstrated from the use of space-time parallelization over the time-sequential integration on Couette flow and Stokes’ second problem. On a transient Couette flow case, at least a 1.5× speedup is achieved, and with a time periodic problem, a speedup of up to 13.7× over the time-sequential case is obtained. In both cases, the speedup is achieved by adding processors and exploring additional parallelization in time. The numerical experiments show the algorithm is promising for CFD applications that can take advantage of the time parallelism. Future work will focus on improving the parallel performance and providing more tests with complex fluid dynamics to demonstrate the full potential of the algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Long-range magnetic order with disordered spin orientations in a high-entropy antiferromagnet

Disorder in magnetic systems typically suppresses long-range order, promoting short-range states such as spin glasses and magnetic clusters. This is particularly prominent in high-entropy materials, characterized by the random distributions of local magnetic entities and exchange interactions. However, in rare exceptions, long-range magnetic order can persist in high-entropy systems, while the microscopic characters and underlying mechanisms remain elusive, especially the magnetic behaviors of individual elements. Here, combining neutron diffraction and resonant soft x-ray scattering, we have conducted an element-specific investigation into the magnetic order of a high-entropy honeycomb-lattice van der Waals material (Mn 1/4 Fe 1/4 Co 1/4 Ni 1/4 )PS 3 . Despite significant atomic disorder, long-range zigzag antiferromagnetic order is observed below 72 K, with all four transition-metal elements participating in a unified phase transition. However, the spin orientations of various elements are distinct, attributed to the competition between single-ion anisotropies and exchange interactions. Our findings showcase a novel form of long-range magnetic order with disordered spin orientations, which is synergically stabilized by distinct magnetic elements in a high entropy magnet, offering a new paradigm for understanding complex magnetic systems.

Shen, Yao [Chinese Academy of Sciences (CAS), Beij↗
DOE OSTI2024

ASb 3 Mn 9 O 19 (A = K or Rb): New Mn-Based 2D Magnetoplumbites with Geometric and Magnetic Frustration

Magnetoplumbites are one of the most broadly studied families of hexagonal ferrites, typically with high magnetic ordering temperatures, making them excellent candidates for permanent magnets. However, magnetic frustration is rarely observed in magnetoplumbites. Herein, the discovery, synthesis, and characterization of the first Mn-based magnetoplumbite, as well as the first magnetoplumbite involving pnictogens (Sb), ASb 3 Mn 9 O 19 (A = K or Rb) are reported. The Mn 3+ ( S = 2) cations, further confirmed by DC magnetic susceptibility and X-ray photoelectron spectroscopy, construct three geometrically frustrated sublattices, including Kagome, triangular, and puckered honeycomb lattices. Magnetic properties measurements revealed strong antiferromagnetic spin–spin coupling as well as multiple low-temperature magnetic features. Heat capacity data does not show any prominent λ-anomaly, suggesting minimal associated magnetic entropy. Moreover, neutron powder diffraction (NPD) implied the absence of long-range magnetic ordering in KSb 3 Mn 9 O 19 down to 3 K. However, several magnetic peaks are observed in RbSb 3 Mn 9 O 19 at 3 K, corresponding to an incommensurate magnetic structure. Interestingly, strong diffuse scattering is seen in the NPD patterns of both compounds at low angles and is analyzed by reverse Monte Carlo refinements, indicating short-range spin ordering related to frustrated magnetism as well as 2D magnetic correlations in ASb 3 Mn 9 O 19 (A = K or Rb).

2-D magnetic correlation↗