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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Timing and Multi-Channel: Novel Method for Determining the Neutrino Mass Ordering from Supernovae

One of the few remaining unknowns in the standard three-flavor neutrino oscillation paradigm is the ordering of neutrino masses. In this work we propose a novel method for determining neutrino mass ordering using the time information on early supernova neutrino events. In a core-collapse supernova, neutrinos are produced earlier than antineutrinos and, depending on the mass ordering which affects the adiabatic flavor evolution, may cause earlier observable signals in $\nu_e$ detection channels than in others. Hence, the time differences are sensitive to the mass ordering. We find that using the time information on the detection of the first galactic supernova events at future detectors like DUNE, JUNO and Hyper-Kamiokande, the mass ordering can already be determined at 2$\sigma$ CL, while $\mathcal{O}(10)$ events suffice for the discovery. Our method does not require high-statistics and could be used within the supernova early warning system (SNEWS) which will have access to the time information on early supernova neutrino events recorded in a number of detectors. The method proposed in this paper also implies a crucial interplay between the mass ordering and the triangulation method for locating supernovae.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

High-order accurate solutions of the point kinetics equations with the spectral deferred correction method

Solving initial value problems with higher-order methods can improve the accuracy of the simulation results or the efficiency of the calculation. In this paper, we apply the spectral deferred correction (SDC) method to solve the initial value problem of the point kinetics equations (PKE). SDC is a stable, robust, and efficient high-order time-integration method capable of an arbitrary order of accuracy. For our implementation we show that it is A-stable for orders up to 8 and the order of accuracy is verified for PKE problems with a range of different reactivities. A 5.-order SDC method was then implemented to solve the exact PKE (EPKE) in the Transient Multilevel (TML) method of MPACT. The error from solutions of the EPKE is shown to be negligible. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Composition-dependent ordering transformations in Pt–Fe nanoalloys

Significance Dynamically understanding the microscopic processes governing ordering transformations has rarely been attained. The situation becomes even more challenging for nanoscale alloys, where the significantly increased surface-area-to-volume ratio not only opens up a variety of additional freedoms to initiate an ordering transformation but also allows for kinetic interplay between the surface and bulk due to their close proximity. We provide direct evidence of the microscopic processes controlling the ordering transformation through the surface–bulk interplay in Pt–Fe nanoalloys and new features rendered by variations in alloy composition and chemical stimuli. These results provide a mechanistic detail of ordering transformation phenomena which are widely relevant to nanoalloys as chemical ordering occurs in most multicomponent materials under suitable environmental bias.

25 ENERGY STORAGE↗
DOE OSTI2022

Testing for pair density wave order in La 1.875 Ba 0.125 CuO 4

We report charge order is commonly believed to compete with superconducting order. An intertwined form of a superconducting wave function, known as pair density wave (PDW) order, has been proposed; however, direct evidence, theoretical or experimental, that it forms the ground state of any cuprate superconductor is lacking. As a test case, we consider La 2-x Ba x CuO 4 with x = 1/8, where charge and spin stripe orders within the CuO 2 planes compete with three-dimensional superconducting order. We report measurements of the superconducting critical current perpendicular to the planes in the presence of an in-plane magnetic field. The variation of the critical current with orientation of the field is inconsistent with a theoretical prediction specific to the PDW model. It appears, instead, that the orientation dependence of the critical-current density might be determined by a minority phase of d-wave superconductivity that is present as a consequence of doped-charge inhomogeneity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Simulation-driven optimization of high-order meshes in ALE hydrodynamics

Here we propose tools for high-order mesh optimization and demonstrate their benefits in the context of multi-material Arbitrary Lagrangian-Eulerian (ALE) compressible shock hydrodynamic applications. The mesh optimization process is driven by information provided by the simulation which uses the optimized mesh, such as shock positions, material regions, known error estimates, etc. These simulation features are usually represented discretely, for instance, as finite element functions on the Lagrangian mesh. The discrete nature of the input is critical for the practical applicability of the algorithms we propose and distinguishes this work from approaches that strictly require analytical information. Our methods are based on node movement through a high-order extension of the Target-Matrix Optimization Paradigm (TMOP). The proposed formulation is fully algebraic and relies only on local Jacobian matrices, so it is applicable to all types of mesh elements, in 2D and 3D, and any order of the mesh. We discuss the notions of constructing adaptive target matrices and obtaining their derivatives, reconstructing discrete data in intermediate meshes, node limiting that enables improvement of global mesh quality while preserving space-dependent local mesh features, and appropriate normalization of the objective function. The adaptivity methods are combined with automatic ALE triggers that can provide robustness of the mesh evolution and avoid excessive remap procedures. The benefits of the new high-order TMOP technology are illustrated on several simulations performed in the high-order ALE application BLAST.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An arbitrarily high-order three-dimensional Cartesian-grid method for reconstructing interfaces from volume fraction fields

Here Tthis work describes a newly developed, arbitrarily high-order Cartesian-grid method for reconstructing material interfaces from a volume fraction field. The method begins by identifying all of the grid cells in the volume fraction field that are intersected by the interface and need to be approximated by the reconstruction scheme. Finite-differences are used to calculate the gradient of the volume fraction field and provide an estimate of the surface normal in all of the interfacial grid cells. Groups of connected grid cells are then identified which all have the same dominant component of the normal vector. This grouping by orientation determines the proper dependent variable to use in the surface reconstruction (e.g. for a 2D curve, this step determines if the surface will be approximated by a function of x or y). A cumulative integral over the surface is constructed and fit using b-splines for two-dimensional problems or tensor-product b-splines for three-dimensional problems. This construction allows for the interface to be recovered through application of the second fundamental theorem of calculus. Fitting the cumulative integral with $\mathscr{N}$ th-order b-splines (or tensor-product b-splines) yields an ($\mathscr{N}$-1) th-order convergence rate of the interface shape. Differentiation of the b-spline interface function(s) allows for the high-order approximation of the normal vector and curvature to be obtained directly anywhere along b-spline. Together, the proposed reconstruction technique can achieve arbitrarily high mesh convergence rates. Validation tests are presented with mesh convergence rates ranging from fourth- to tenth-order.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

High-order dimensionally-split Cartesian embedded boundary method for non-dissipative schemes

Centered finite-difference schemes are commonly used for high-fidelity turbulent flow simulations in canonical configurations because of their non-dissipative property and computational efficiency. However, their use in flow simulations over complex geometries is limited by the requirements of a structured grid and a stable boundary treatment in the absence of artificial (numerical) dissipation. Cartesian embedded boundary (EB) approaches provide an efficient structured-grid framework to apply difference schemes over complex domains. However, they are often restricted to low orders of accuracy because of numerical instabilities at the embedded boundaries and the issues of small-cell problem that are difficult to address with high-order accuracy. The present work discusses a systematic approach to obtain high-order EB methods with non-dissipative centered schemes in the interior. This approach, based on satisfying the primary and secondary conservation conditions, is employed to derive EB schemes that are up to sixth-order accurate in the interior and fourth-order accurate globally for hyperbolic, parabolic as well as incompletely parabolic problems. The proposed finite-difference discretization is, by construction, dimensionally split and addresses the small-cell problem without any cell/geometry transformations, thus, highly simplifying implementation in a flow solver. Various linear and non-linear numerical tests are performed to evaluate the stability and the accuracy of the proposed EB schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A high-order finite difference method for moving immersed domain boundaries and material interfaces

Here, we present a high-order sharp treatment of immersed moving domain boundaries and material interfaces, and apply it to the advection-diffusion equation in two and three dimensions. The spatial discretization combines dimension-split finite difference schemes with an immersed boundary treatment based on a weighted least-squares reconstruction of the solution, providing stable discretizations with up to sixth order accuracy for diffusion terms and third order accuracy for advection terms. The temporal discretization relies on a novel strategy for maintaining high-order temporal accuracy in problems with moving boundaries that minimizes implementation complexity and allows arbitrary explicit or diagonally-implicit Runge-Kutta schemes. The approach is broadly compatible with popular PDE-specialized Runge-Kutta time integrators, including low-storage, strong stability preserving, and diagonally implicit schemes. Through numerical experiments we demonstrate that the full discretization maintains high-order spatial and temporal accuracy in the presence of complex 3D geometries and for a range of boundary conditions, including Dirichlet, Neumann, and flux conditions with large jumps in coefficients.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Features of relaxation modes in soft magnetic metallic glasses and their correlation with magnetic ordering

Although the anneal-induced relaxation process of Fe-based metallic glasses has been widely reported, there is a lack of a systematic description of the relaxation modes of atoms. Here, in this work, the relaxation modes of FeSiBPC metallic glasses and their relationship with magnetic ordering have been investigated. A two-stage relaxation behavior has been observed, and the pinning effect of Curie temperature is found in the first stage. The short-range ordering involved in the β relaxation in the first stage is confined and facilitated by magnetic ordering. It causes a smaller E β compared to that of non-magnetic metallic glasses. In the second stage, the pinning effect disappears. It is mainly the medium-range rearrangements involved in α relaxation that begins to yield the influence on magnetic ordering. This work uncovers the inner magnet-influenced β relaxation in Fe-based metallic glasses and provides insights into the interaction between structural arrangement during relaxation and magnetic ordering.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Impact of chemical ordering on thermodynamic properties of point defects and Xe substitutional in U-10Mo

The accurate knowledge of defect energetics is critical to understanding the aging and irradiation behavior of U-10Mo nuclear fuel, which is selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). Furthermore, using hybrid molecular dynamics and Monte Carlo (MDMC) simulation, we studied the impact of chemical ordering on the formation energies of vacancies, interstitials, and the solution energy of the Xe substitutionals. Instead of forming a random solid solution (RSS), substantial short-range-order (SRO) develops in U-10Mo, particularly at low temperatures. Mo atoms are found to repel each other and prefer U-rich local atomic environments within the 1st nearest neighbor (1NN) cutoff. Compared to the case of a RSS, the state with equilibrated Mo ordering shifts the distributions of vacancy and interstitial formation energies due to the dependence of defect energies on the local atomic environment, without a clear effect on Xe solution energy. In the operation temperature range (100–250 °C) of U-Mo fuels, neglecting SRO can lead to an inaccurate estimate of thermal equilibrium point defect concentrations by over an order of magnitude and incorrectly predict the preference among different types of dumbbells, highlighting the critical importance of accounting for the impact of chemical ordering for accurate atomistic calculations of defect properties.

Molecular dynamics and Monte Carlo (MDMC)↗
DOE OSTI2022

Proper orthogonal decomposition based reduced-order modeling of flux-Limited gray thermal radiation

Here, in this work, a proper orthogonal decomposition (POD) based reduced-order model (ROM) is developed to solve gray, flux-limited thermal radiation diffusion. We focus on the variable opacity radiation penetration benchmark posed by Olson, Auer, and Hall. The T -3 relationship for opacity in conjunction with high-temperature radiation penetrating an initially cold material produces a strong thermal radiation shock. This class of problems is particularly challenging for standard POD-based reduced-order modeling due to the nonlinearities presented by 1) the T 4 source term, and 2) flux-limited diffusion operator. To address these challenges and develop a cost competitive ROM, we employ a “hyper-reduction” technique through discrete empirical interpolation (DEIM) and allow for adaptive reduced-order projections through principal interval decomposition (PID). Performance of the proposed methodology is quantified by comparing the cost savings and accuracy relative to a full-order computation. Reference solutions and snapshot data are obtained through a full-order calculation performed by the University of Chicago maintained astrophysics code, FLASH. For consistency and potential extensibility, the developed ROM is also implemented in FLASH. We find that in the initialization regime, where the thermal radiation wave is initially created by the warming of the material, this class of problems is highly reducible and suitable for POD-based ROMs. However, the strong convective nature of the wave propagation regime is less reducible and more challenging to create an efficient ROM.

42 ENGINEERING↗
DOE OSTI2025

Slender-body approach for computing second-order wave loads in the frequency domain

This work presents a slender-body approach to evaluate the second-order wave loads acting on a floating structure in the frequency domain. The approach is in the same spirit as the common use of Morison’s equation to approximate the wave loads without solving the radiation/diffraction problem. To do so, we employ Rainey’s equation, which can be seen as an extension of the inertial part of Morison’s equation to include nonlinear effects. We introduce modifications to Rainey’s formulation in order to evaluate wave kinematics at the mean body position instead of the original approach of considering instantaneous displacements. We also propose a simple approximation to partially account for wave scattering effects on the second-order loads based on the analytical solution of a surface-piercing bottom-mounted vertical circular cylinder. Though limited to structures composed of cylinders, this slender-body approach is orders of magnitude faster than computing second-order wave coefficients with a radiation/diffraction code. We implemented this approach for difference-frequency (slow drift) loads in an open-source frequency-domain floating wind turbine model. We present comparisons against results obtained with radiation/diffraction theory for three reference floating wind turbine designs: the OC3-Hywind spar, the OC4-DeepCwind semisubmersible, and the VolturnUS-S semisubmersible. In general, the results show that the proposed slender-body approach with the correction to approximate wave scattering effects provides useful estimations of the difference-frequency wave loads and the resulting motions of the floater.

17 WIND ENERGY↗
DOE OSTI2026

Second-order wave excitation forces in WEC-Sim/MOST: Implementation, experimental validation, and code-to-code comparison

Accurate prediction of second-order hydrodynamic loads is essential for floating bodies, including floating offshore wind turbines, wave energy converters, and hybrid wind–wave platforms. These nonlinear effects, arising from both sum- and difference-frequency forcing, are critical for capturing key response characteristics but remain challenging to model efficiently. In this work, we extend the open-source Wave Energy Converter Simulator / MATLAB for Offshore Simulation Tool by implementing second-order wave excitation forces, supporting both the full Quadratic Transfer Function formulation and the Newman approximation. The full Quadratic Transfer Function method is used for all code-to-code comparisons and experimental validation, while the Newman approximation is provided as a computationally lighter alternative. To benchmark the new capability, we perform a code-to-code comparison with OpenFAST and OrcaFlex. We then validate the enhanced model using wave-tank measurements of a 1:96 scale DeepCwind semi-submersible, showing that second-order effects are required to reproduce platform motions. The implementation employs a computationally efficient pre-computation strategy for second-order wave excitation forces, reducing simulation cost while maintaining engineering accuracy. Overall, this work advances the tool as an open-source and versatile tool for modelling floating offshore renewable-energy systems requiring second-order hydrodynamic fidelity.

17 WIND ENERGY↗
DOE OSTI2021

Rapid Atomic Ordering Transformation toward Intermetallic Nanoparticles

Chemically ordered intermetallic nanoparticles are promising candidates for energy-related applications such as electrocatalysis. However, the synthesis of intermetallics generally requires long annealing (several hours) to achieve the ordered structure, which causes nanoparticles agglomeration and diminished performance, particularly for catalysis. Herein, we demonstrate a new rapid Joule heating approach that can synthesize highly ordered and well-dispersed intermetallic nanoparticles. As a proof-of-concept, we synthesized fully ordered Pd 3 Pb intermetallic nanoparticles that feature small size distribution (~6 nm). Computational analysis of the L1 2 Pd 3 Pb material suggests that this rapid atomic ordering transformation can be attributed to a vacancy-mediated diffusion mechanism. Moreover, the nanoparticles demonstrate excellent electrocatalytic activity and exceptional stability for the oxygen reduction reaction (ORR), retaining >95% of the current density over 10 h of chronoamperometry test with negligible structural and compositional changes. Furthermore, this study demonstrates a new strategy of providing a new direction for intermetallic synthesis and catalyst discovery.

25 ENERGY STORAGE↗
DOE OSTI2021

First-order antiferromagnetic transitions of SrMn 2 P 2 and CaMn 2 P 2 single crystals containing corrugated-honeycomb Mn sublattices

Significance With rare exceptions, an antiferromagnetic (AFM) transition in zero magnetic field is thermodynamically of second order where the thermal-average magnetic moments of the magnetic atoms (ordered moments) vary continuously on cooling through the AFM ordering temperature T N with no latent heat at the transition. Such materials include the AFM pnictides CaMn 2 As 2 , SrMn 2 As 2 , CaMn 2 Sb 2 , SrMn 2 Sb 2 , and CaMn 2 Bi 2 . Here we demonstrate that the closely related SrMn 2 P 2 and CaMn 2 P 2 insulators instead exhibit first-order AFM transitions at T N = 53 and 70 K, respectively, where the heat capacity exhibits a latent heat at T N . The mechanism causing these first-order transitions remains to be explained, but its understanding may lead to the development of novel magnetic materials of technological interest.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Competing magnetic orders in quantum critical Sr 3 Ru 2 O 7

We investigated Sr 3 RuO 7 , a quantum critical metal that shows a metamagnetic quantum phase transition and electronic nematicity, through density functional calculations. These predict a ferromagnetic ground state in contrast to the experimentally observed paramagnetism, raising the question of competing magnetic states and associated fluctuations that may suppress magnetic order. We did a search to identify such low-energy antiferromagnetically ordered metastable states. We find that the lowest-energy antiferromagnetic state has a striped order. This corresponds to the E-type order that has been shown to be induced by Mn alloying. We also note significant transport anisotropy in this E-type ordered state. Furthermore, these results are discussed in relation to experimental observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Spin-density-wave order controlled by uniaxial stress in CeAuSb 2

We report the tetragonal heavy-fermion compound CeAuSb 2 (space group P4/nmm ) exhibits incommensurate spin-density wave (SDW) order below T N ≈ 6.5 K with the propagation vector q A = (δ A , δ A , 1/2 ). The application of uniaxial stress along the [010] direction induces a sudden change in the resistivity ratio ρ a /ρ b at a compressive strain of ε ≈ - 0.5 %. Here we use neutron scattering to show that the uniaxial stress induces a first-order transition to a SDW state with a different propagation vector (0, δ B ,1/2) with δ B = 0.25.The magnetic structure of the new (B) phase consists of Ce layers with ordered moments alternating with layers with zero moment stacked along the c axis. The ordered layers have an up-up-down-down configuration along the b axis. This is an unusual situation in which the loss of spatial inversion in a metallic system is driven by the magnetic order. We argue that the change in SDW wave vector leads to Fermi-surface reconstruction and a concomitant change in the transport properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗