Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Sampling and resolution characteristics in reduced order models of shallow water equations: Intrusive vs nonintrusive

We investigate the sensitivity of reduced order models (ROMs) to training data spatial resolution as well as sampling rate. In particular, we consider proper orthogonal decomposition (POD), coupled with Galerkin projection (POD-GP), as an intrusive reduced order modeling technique. For non-intrusive ROMs, we consider two frameworks. The first is using dynamic mode decomposition (DMD), and the second is based on artificial neural networks (ANNs). For ANN, we utilized a residual deep neural network, and for DMD we have studied two versions of DMD approaches; one with hard thresholding and the other with sorted bases selection. Also, we highlight the differences between mean-subtracting the data (centering) and using the data without mean-subtraction. We tested these ROMs using a system of 2D shallow water equations for four different numerical experiments, adopting combinations of sampling rates and spatial resolutions. For these cases, we found that the DMD basis obtained with hard threshodling is sensitive to sampling rate. The sorted DMD algorithm helps to mitigate this problem and yields more stabilized and converging solution. Furthermore, we demonstrate that both DMD approaches without mean subtraction provide significantly more accurate results than DMD with mean-subtracting the data. On the other hand, POD is relatively insensitive to sampling rate and yields better representation of the flow field. Meanwhile, spatial resolution has little effect on both POD and DMD performances. Here, the numerical results reveal that an ANN on POD subspace (POD-ANN) performs remarkably better than POD-GP and DMD in capturing system dynamics, even with a small number of modes.

42 ENGINEERING↗
DOE OSTI2021

Nitride-Oxide-Metal Heterostructure with Self-Assembled Core–Shell Nanopillar Arrays: Effect of Ordering on Magneto-Optical Properties

Magneto-optical (MO) coupling incorporates photon-induced change of magnetic polarization that can be adopted in ultrafast switching, optical isolators, mode convertors, and optical data storage components for advanced optical integrated circuits. However, integrating plasmonic, magnetic, and dielectric properties in one single material system poses challenges since one natural material can hardly possess all these functionalities. Here, co-deposition of a three-phase heterostructure composed of a durable conductive nitride matrix with embedded core–shell vertically aligned nanopillars, is demonstrated. Additionally, the unique coupling between ferromagnetic NiO core and atomically sharp plasmonic Au shell enables strong MO activity out-of-plane at room temperature. Further, a template growth process is applied, which significantly enhances the ordering of the nanopillar array. The ordered nanostructure offers two schemes of spin polarization which result in stronger antisymmetry of Kerr rotation. The presented complex hybrid metamaterial platform with strong magnetic and optical anisotropies is promising for tunable and modulated all-optical-based nanodevices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High-Order Mesh hr-adaptivity for Surface Fitting to Implicit Geometries

We present an ℎ𝑟-adaptivity framework for morphing a given mesh to fit a target surface prescribed as the zero isocontour of a discrete function. In this framework, high-order meshing is posed as a variational minimization problem that depends on the mesh quality prescribed via the target matrix optimization paradigm (TMOP) and position of a subset of mesh nodes with respect to the target surface. The proposed formulation ensures that the variational problem is converged and mesh quality degradation near the surface is limited, even when the mesh topology is incompatible with the target surface. Additionally, a mesh subset-based approach and ℎ-refinement is introduced to efficiently increase fitting accuracy while reducing the computational cost of the mesh morphing problem. The ℎ𝑟-adaptivity technique extends to different element types in two- and three-dimensions, and can be used in existing finite element and spectral element frameworks to obtain high-order body-fitted meshes. Various numerical experiments demonstrate the robustness and accuracy of the fitting approach for problems of practical interest such as Lagrangian hydrodynamics and topology optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Third order QCD predictions for fiducial W -boson production

Measurements of W-boson production at the LHC have reached percent-level precision and impose challenging demands on theoretical predictions. Such predictions directly limit the precision of measurements of fundamental quantities such as the W-boson mass and the weak mixing angle. A dominant source of uncertainty in predictions is from higher-order QCD effects. We present a calculation of W-boson production at the level of $α^3_s$ at fixed order and including transverse-momentum resummation. We further show predictions for a direct comparison with low-pileup ATLAS transverse-momentum and fiducial cross-section measurements at $\sqrt{s}$ = 5.02 TeV. We discuss in detail the impact of modern PDFs. Our calculation including the matching to W+jet production at NNLO will be publicly available the upcoming CuTe-MCFM release and allows for theory-data comparison at the state-of-the-art level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Factorial growth at low orders in perturbative QCD: control over truncation uncertainties

A method, known as “minimal renormalon subtraction”, relates the factorial growth of a perturbative series (in QCD) to the power p of a power correction Λ p /Q p . (Λ is the QCD scale, Q some hard scale.) Here, the derivation is simplified and generalized to any p , more than one such correction, and cases with anomalous dimensions. Strikingly, the well-known factorial growth is seen to emerge already at low or medium orders, as a consequence of constraints on the Q dependence from the renormalization group. The effectiveness of the method is studied with the gluonic energy between a static quark and static antiquark (the “static energy”). Truncation uncertainties are found to be under control after next-to-leading order, despite the small exponent of the power correction ( p = 1) and associated rapid growth seen in the first four coefficients of the perturbative series.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

hr-Adaptivity for nonconforming high-order meshes with the target matrix optimization paradigm

In this study, we present an \(hr\) -adaptivity framework for optimization of high-order meshes. This work extends the r -adaptivity method by Dobrev et al. (Comput Fluids, 2020), where we utilized the Target-Matrix Optimization Paradigm (TMOP) to minimize a functional that depends on each element’s current and target geometric parameters: element aspect-ratio , size , skew , and rotation . Since fixed mesh topology limits the ability to achieve the target size and aspect-ratio at each position, in this paper, we augment the r -adaptivity framework with nonconforming adaptive mesh refinement to further reduce the error with respect to the target geometric parameters. The proposed formulation, referred to as \(hr\) -adaptivity, introduces TMOP-based quality estimators to satisfy the aspect-ratio target via anisotropic refinements and size target via isotropic refinements in each element of the mesh. The methodology presented is purely algebraic, extends to both simplices and hexahedra/quadrilaterals of any order, and supports nonconforming isotropic and anisotropic refinements in 2D and 3D. Using a problem with a known exact solution, we demonstrate the effectiveness of \(hr\) -adaptivity over both r - and \(h\) -adaptivity in obtaining similar accuracy in the solution with significantly fewer mesh nodes. We also present several examples that show that \(hr\) -adaptivity can help satisfy geometric targets even when \(r\) -adaptivity fails to do so, due to the topology of the initial mesh.

42 ENGINEERING↗
DOE OSTI2021

On the Order Parameter of the Continuous Phase Transition in the Classical and Quantum Mechanical Limits

The mean field theory is revisited in the classical and quantum mechanical limits. Taking into account the boundary conditions at the phase transition and the third law of the thermodynamics the physical properties of the ordered and disordered phases were reported. The equation for the order parameter predicts the occurrence of a saturation of Ψ 2 = 1 near Θ S , the temperature below the quantum mechanical ground state is reached. The theoretical predictions are also compared with high resolution thermal expansion data of SrTiO 3 monocrystalline samples and other some previous results. An excellent agreement has been found suggesting a universal behavior of the theoretical model to describe continuous structural phase transitions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A sweeping positivity-preserving high-order finite difference WENO scheme for Euler equations

We develop a simple, high-order, conservative and robust positivity-preserving sweeping procedure for the density and the nonlinear pressure function in the compressible Euler equations. Using the scaling limiter in Zhang and Shu (J Comput Phys 229:3091–3120, 2010), we obtain a non-trivial extension of the scalar sweeping technique in Liu et al. (J Sci Comput 73:1028–1071, 2017) for the positivity of pressure. The sweeping procedure developed in this paper is a post-processing technique, which can be applied to any concave functions of the conserved variables in hyperbolic conservation law systems. Thus, it has applications beyond the Euler equations. This procedure preserves positivity and conservation of physical quantities without destroying the accuracy of the underlying scheme. The algorithm works for general schemes including finite difference, finite volume and discontinuous Galerkin methods; however, in this paper we focus on finite difference weighted essentially non-oscillatory (WENO) methods. As a result, we provide numerical tests of the fifth-order finite difference WENO scheme to demonstrate the accuracy and robustness of the technique.

Compressible Euler equations↗
DOE OSTI2022

Well-Balanced Second-Order Convex Limiting Technique for Solving the Serre–Green–Naghdi Equations

In this article, we introduce a numerical method for approximating the dispersive Serre–Green–Naghdi equations with topography using continuous finite elements. The method is an extension of the hyperbolic relaxation technique introduced in Guermond et al. (J Comput Phys 450:110809, 2022). It is explicit, second-order accurate in space, third-order accurate in time, and is invariant-domain preserving. It is also well balanced and parameter free. Special attention is given to the convex limiting technique when physical source terms are added in the equations. The method is verified with academic benchmarks and validated by comparison with laboratory experimental data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

First-principles and cluster expansion study of the effect of magnetism on short-range order in Fe–Ni–Cr austenitic stainless steels

Short-range order (SRO), the regular and predictable arrangement of atoms over short distances, alters the mechanical properties of technologically relevant structural materials such as medium/high entropy alloys and austenitic stainless steels. In this study, we present a generalized spin cluster expansion (CE) model and show that magnetism is a primary factor influencing the level of SRO present in austenitic Fe-Ni-Cr alloys. The spin CE consists of a chemical cluster expansion combined with an Ising model for Fe-Ni-Cr austenitic alloys. It explicitly accounts for local magnetic exchange interactions, thereby capturing the effects of finite temperature magnetism on SRO. Model parameters are obtained by fitting to a first-principles data set comprising both chemically and magnetically diverse FCC configurations. The magnitude of the magnetic exchange interactions are found to be comparable to the chemical interactions. Compared to a conventional implicit magnetism CE built from only magnetic ground state configurations, the spin CE shows improved performance on several experimental benchmarks over a broad spectrum of compositions, particularly at higher temperatures due to the explicit treatment of magnetic disorder. We find that SRO is strongly influenced by alloy Cr content, since Cr atoms prefer to align antiferromagnetically with nearest neighbors but become magnetically frustrated with increasing Cr concentration. Using the spin CE, we predict that increasing the Cr concentration in typical austenitic stainless steels promotes the formation of SRO and increases order-disorder transition temperatures. Furthermore, this study underscores the significance of considering magnetic interactions explicitly when exploring the thermodynamic properties of complex transition metal alloys. It also highlights guidelines for customizing SRO through adjustments of alloy composition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Local chemical ordering of a neutron-irradiated CrFeMnNi compositionally complex alloy

While ion-irradiation studies are a critical first step in studying compositionally complex alloys (CCAs) for nuclear applications, they do not capture all the microstructural changes occurring under the low irradiation dose rates and different particles’ scattering patterns experienced in a nuclear reactor setting. To explore these phenomena in reactor-relevant conditions for the first time in CCA, the single-phase solid-solution Cr 10 Fe 30 Mn 30 Ni 30 was neutron irradiated up to 6.61 displacements per atom at 395 and 579 °C. Irradiation-enhanced local chemical ordering (LCO) well beyond the range of short range ordering was observed, and is predicted to be the precursor to the precipitation of a coherent Ni-Mn L1 0 phase and a Cr-rich α’ phase, though TEM analysis did not indicate the presence of either in any irradiation condition. The line density of faulted dislocation loops decreased from 6.47 to 1.69 ∙ 10 15 m -2 from 3.43 to 6.61 dpa at 579 °C despite no appreciable faulted loop content in the unirradiated material. LCO is expected to increase the complexity of the energy landscape within this alloy, restricting interstitial point defect mobility and creating local regions of greater stacking fault energy. These contribute to the negative correlation between irradiation dose and faulted dislocation loop density in this study, as well as the lack of void swelling observed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Driving rapid atomic order in MnAl via low-magnitude magnetic field annealing

Application of a mild (60 mT), uniform magnetic field during short-term thermal treatment of kinetically retained, atomically disordered (paramagnetic) ε-MnAl was found to deliver a significant ~50 % increase in the formation of L1 0 atomically ordered (ferromagnetic) τ-MnAl product phase, compared to that produced by conventional (i.e., zero-field) annealing under identical thermal conditions. The magnetic field, applied in a passive closed-circuit configuration during annealing, induced significant changes in the structural, magnetic, and phase evolution of the material. Computational results based on electronic structure calculations demonstrate that the effective magnetic susceptibility of τ-MnAl is sensitive to the orientation, rather than the magnitude, of an applied magnetic field in the vicinity of the Curie temperature. The uniaxial magnetocrystalline anisotropy of the L1 0 structure is proposed to act as a filter for selective propagation of the population of τ-MnAl variants that are favorably aligned with the applied field. In this manner, crystallographic “gridlock” is alleviated that would otherwise arise from the coexistence of multiple, energetically equivalent τ-phase variants within the parent ε-phase matrix. These results confirm that static, low-magnitude magnetic field annealing is able to accelerate L1 0 atomic ordering in the MnAl system and likely can exert similar influences in relevant magnetic systems, facilitating efficient tailoring of structure-sensitive magnetic properties for the manufacture of magnetic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluation of dual-weighted residual and machine learning error estimation for projection-based reduced-order models of steady partial differential equations

Projection-based reduced-order models (pROMs) show great promise as a means to accelerate many-query applications such as forward error propagation, solving inverse problems, and design optimization. In order to deploy pROMs in the context of high-consequence decision making, accurate error estimates are required to determine the region(s) of applicability in the parameter space. The following paper considers the dual-weighted residual (DWR) error estimate for pROMs and compares it to another promising pROM error estimate, machine learned error models (MLEM). Here, we show how DWR can be applied to ROMs and then evaluate DWR on two partial differential equations (PDEs): a two-dimensional linear convection–reaction–diffusion equation, and a three-dimensional static hyper-elastic beam. It is shown that DWR is able to estimate errors for pROMs extrapolating outside of their training set while MLEM is best suited for pROMs used to interpolate within the pROM training set.

42 ENGINEERING↗
DOE OSTI2021

Kink-antikink collisions and multi-bounce resonance windows in higher-order field theories

In this work, we study collisions of coherent structures in higher-order field-theoretic models, such as the Φ 8 , Φ 10 and Φ 12 ones. The main distinguishing feature, of the example models considered herein, is that the collision arises due to the long-range interacting algebraic tails of these solitary waves. We extend the approach to suitably initialize the relevant kinks, in the additional presence of finite initial velocity, in order to minimize the dispersive wave radiation potentially created by their slow spatial decay. We find that, when suitably initialized, these models still feature the multi-bounce resonance windows earlier found in models in which the kinks bear exponential tails, such as the Φ 4 and Φ 6 field theories among others. Also present is the self-similar structure of the associated windows with three- and more-bounce windows at the edges of two- and lower-bounce ones. Moreover, phenomenological but highly accurate (and predictive), scaling relations are derived for the dependence of the time between consecutive collisions and, e.g., the difference in kinetic energy between the incoming one and the critical one for one-bounces. Such scalings are traced extensively over two-bounce collision windows throughout the three models, hinting at the possibility of an analytical theory in this direction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A fast matrix-free approach to the high-order control volume finite element method with application to low-Mach flow

Here, a fast matrix-free formulation of the control volume finite element method is presented, requiring much less memory and computational work than previous efforts. The method is implemented and evaluated as a solver for low-Mach flow, including the evaluation of a preconditioning strategy for the pressure Poisson equation. The efficiency and scaling with polynomial order is evaluated on simple turbulent flows of interest, with appropriate solution quality metrics, and compared with a reference node-centered finite volume discretization. For a turbulent channel flow test, we show improvement in computational work for a given accuracy with the high-order scheme. The performance on a GPU accelerated platform is also investigated, with benefit shown for the matrix-free discretization.

42 ENGINEERING↗
DOE OSTI2025

General field evaluation in high-order meshes on GPUs

Robust and scalable function evaluation at any arbitrary point in the finite/spectral element mesh is required for querying the partial differential equation solution at points of interest, comparison of solution between different meshes, and Lagrangian particle tracking. This is a challenging problem, particularly for high-order unstructured meshes partitioned in parallel with MPI, as it requires identifying the element that overlaps a given point and computing the corresponding reference space coordinates. Here, we present a robust and efficient technique for general field evaluation in large-scale high-order meshes with quadrilaterals and hexahedra. In the proposed method, a combination of globally partitioned and processor-local maps are used to first determine a list of candidate MPI ranks, and then locally candidate elements that could contain a given point. Next, element-wise bounding boxes further reduce the list of candidate elements. Finally, Newton’s method with trust region is used to determine the overlapping element and corresponding reference space coordinates. Since GPU-based architectures have become popular for accelerating computational analyses using meshes with tensor-product elements, specialized kernels have been developed to utilize the proposed methodology on GPUs. The method is also extended to enable general field evaluation on surface meshes. The paper concludes by demonstrating the use of the proposed method in various applications ranging from mesh-to-mesh transfer during r-adaptivity to Lagrangian particle tracking.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering

Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Detecting anomalous sequences in electronic health records using higher-order tensor networks

Detecting anomalous sequences is an integral part of building and protecting modern large-scale health information technology (HIT) systems. These HIT systems generate a large volume of records of patients’ state and significant events, which provide a valuable resource to help improve clinical decisions, patient care processes, and other issues. However, detecting anomalous sequences in electronic health records (EHR) remains a challenge in healthcare applications for several reasons, including imbalances in the data, complexity of relationships between events in the sequence, and the curse of dimensionality. Conventional anomaly detection methods use the finite sequence of events to discriminate sequences. They fail to incorporate salient event details under variable higher-order dependencies (e.g., duration between events) that can provide better discrimination of sequences in their models. To address this problem, we propose event sequence and subsequence anomaly detection algorithms that (1) use network-based representations of interactions in the data, (2) account for variable higher-order dependencies in the data, and (3) incorporate events duration for adequate discrimination of the data. The proposed approach identifies anomalies by monitoring the change in the graph after the test sequence is removed from the network. The change is quantified using graph distance metrics so that dramatic changes in the network can be attributed to the removed sequence. Furthermore, the proposed subsequence algorithm recommends plausible paths and salient information for the detected anomalous subsequences. Our results show that the proposed event sequence anomaly detection algorithm outperforms the baseline methods for both synthetic data and real-world EHR data.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗