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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Out-of-time-order correlators and Lyapunov exponents in sparse SYK

We use a combination of analytical and numerical methods to study out-of-time order correlators (OTOCs) in the sparse Sachdev-Ye-Kitaev (SYK) model. We find that at a given order of N, the standard result for the q-local, all-to-all SYK, obtained through the sum over ladder diagrams, is corrected by a series in the sparsity parameter, k. We present an algorithm to sum the diagrams at any given order of 1/(kq) n . We also study OTOCs numerically as a function of the sparsity parameter and determine the Lyapunov exponent. We find that numerical stability when extracting the Lyapunov exponent requires averaging over a massive number of realizations. This trade-off between the efficiency of the sparse model and consistent behavior at finite N becomes more significant for larger values of N.

2D gravity↗
DOE OSTI2024

Calculation of machine precision second order derivatives using dual-complex numbers

It is well known that both complex and dual numbers can be employed to obtain machine precision first-order derivatives; however, neither, on their own, can compute machine precision 2nd order derivatives. To address this limitation, it is demonstrated in this paper that combined dual-complex numbers can be used to compute machine precision 1st and 2nd order derivatives. The dual-complex approach is simpler than utilizing multicomplex or hyper-dual numbers as existing dual libraries can be used as is or easily augmented to accept complex numbers, and the complexity of developing, integrating, and deploying multicomplex or hyper-dual libraries is avoided. The efficacy of this approach is demonstrated for both univariate and multivariate functions. Finally, source code examples using the Python, Julia, and Mathematica languages are provided as supplemental material.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Short- and medium-range orders in Al90Tb10 glass and their relation to the structures of competing crystalline phases

Molecular dynamics simulations using an interatomic potential developed by artificial neural network deep machine learning are performed to study the local structural order in Al90Tb10 metallic glass. We show that more than 80% of the Tb-centered clusters in Al90Tb10 glass have short-range order (SRO) with their 17 first coordination shell atoms stacked in a ‘3661’ or ‘15551’ sequence. Medium-range order (MRO) in Bergman-type packing extended out to the second and third coordination shells is also clearly observed. Analysis of the network formed by the ‘3661’ and ‘15551’ clusters show that ~82% of such SRO units share their faces or vertexes, while only ~6% of neighboring SRO pairs are interpenetrating. Such a network topology is consistent with the Bergman-type MRO around the Tb-centers. Moreover, crystal structure searches using genetic algorithm and the neural network interatomic potential reveal several low-energy metastable crystalline structures in the composition range close to Al 90 Tb 10 . Some of these crystalline structures have the ‘3661’ SRO while others have the ‘15551’ SRO. While the crystalline structures with the ‘3661’ SRO also exhibit the MRO very similar to that observed in the glass, the ones with the ‘15551’ SRO have very different atomic packing in the second and third shells around the Tb centers from that of the Bergman-type MRO observed in the glassy phase.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Critical cooling rates for amorphous-to-ordered complexion transitions in Cu-rich nanocrystalline alloys

Amorphous complexions in nanocrystalline metals have the potential to improve mechanical properties and radiation tolerance, as well as resistance to grain growth. In this study, the stability of amorphous complexions in binary and ternary Cu-based alloys is investigated by observing the effect of cooling rate from high temperature on the occurrence of amorphous-to-ordered complexion transitions. Bulk Cu-Zr and Cu-Zr-Hf alloy samples were annealed to induce boundary premelting and then quenched through a procedure that induces a gradient of local cooling rate through the sample height. Amorphous complexion thickness distributions were found to be invariant to local cooling rate in the Cu-Zr-Hf alloy, demonstrating enhanced stability of the amorphous complexion structure compared to the Cu-Zr alloy, which had thinner amorphous complexions in the regions that were slowly cooled. The experimental results are used to construct time-temperature-transformation diagrams for the amorphous-to-ordered complexion transition in both the binary and ternary alloys, enabling a deeper understanding of the influence of cooling rate and grain boundary chemistry on complexion transitions. In conclusion, the critical cooling rate necessary to avoid complexion transitions in the ternary alloy is found to be at least three orders of magnitude slower than that for the binary alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Topological order in Mott insulators

Here in this commemorative issue for Phil Anderson, I reproduce parts of a joint unpublished paper that first got me thinking about topological order in Mott insulators. With Phil, we explored aspects of the momentum distribution function and the non-classical rotational inertia in the Bose Hubbard model close to the quantum phase transition. My interest on this question continued, and later in the article I take up the question of topological order in the Kitaev Mott insulator on the honeycomb lattice with bond-dependent interactions in a finite magnetic field. Kitaev obtained the exact solution for zero field that demonstrated an exotic Z 2 quantum spin liquid phase with topological order. In a finite magnetic field, we discovered a new intermediate U(1) gapless quantum spin liquid phase, sandwiched between the known gapped Kitaev spin liquid and a polarized phase, that harbors fractionalized neutral fermionic excitations, dubbed spinons, that remarkably form a Fermi surface in a Mott insulator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Galerkin formulations of isogeometric shell analysis: Alleviating locking with Greville quadratures and higher-order elements

Here we propose new quadrature schemes that asymptotically require only four in-plane points for Reissner–Mindlin shell elements and nine in-plane points for Kirchhoff–Love shell elements in B-spline and NURBS-based isogeometric shell analysis, independent of the polynomial degree p of the elements. The quadrature points are Greville abscissae associated with pth-order B-spline basis functions whose continuities depend on the specific Galerkin formulations, and the quadrature weights are calculated by solving a linear moment fitting problem in each parametric direction. The proposed shell element formulations are shown through numerical studies to be rank sufficient and to be free of spurious modes. The studies reveal comparable accuracy, in terms of both displacement and stress, compared with fully integrated spline-based shell elements, while at the same time reducing storage and computational cost associated with forming element stiffness and mass matrices and force vectors. The high accuracy with low computational cost makes the proposed quadratures along with higher-order spline bases, in particular polynomial orders, p = 5 and 6, good choices for alleviating membrane and shear locking in shells.

42 ENGINEERING↗
DOE OSTI2022

Stress-constrained topology optimization of lattice-like structures using component-wise reduced order models

We report lattice-like structures can provide a combination of high stiffness with light weight that is useful in many applications, but a resolved finite element mesh of such structures results in a computationally expensive discretization. This computational expense may be particularly burdensome in many-query applications, such as optimization. We develop a stress-constrained topology optimization method for lattice-like structures that uses component-wise reduced order models as a cheap surrogate, providing accurate computation of stress fields while greatly reducing run time relative to a full order model. We demonstrate the ability of our method to produce large reductions in mass while respecting a constraint on the maximum stress in a pair of test problems. The ROM methodology provides a speedup of about 150x in forward solves compared to full order static condensation and provides a relative error of less than 5% in the relaxed stress.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

SIDIS-RC EvGen: A Monte-Carlo event generator of semi-inclusive deep inelastic scattering with the lowest-order QED radiative corrections

SIDIS-RC EvGen is a C++ standalone Monte-Carlo event generator for studies of semi-inclusive deep inelastic scattering (SIDIS) processes at medium to high lepton beam energies. In particular, the generator contains binary and library components for generating SIDIS events and calculating cross sections for unpolarized or longitudinally polarized beam and unpolarized, longitudinally or transversely polarized target. The structure of the generator incorporates transverse momentum-dependent parton distribution and fragmentation functions, whereby we obtain multi-dimensional binned simulation results, which will facilitate the extraction of important information about the three-dimensional nucleon structure from SIDIS measurements. In order to build this software, we have used recent elaborate QED calculations of the lowest-order radiative effects, applied to the leading order Born cross section in SIDIS. Here, in this paper, we provide details on the theoretical formalism as well as the construction and operation of SIDIS-RC EvGen, e.g., how we handle the event generation process and perform multi-dimensional integration. We also provide example programs, flowcharts, and numerical results on azimuthal transverse single-spin asymmetries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The magnetic order in multiferroic DyMnO 3

We report that with flexibility in tuning their electric and magnetic properties, multiferroics can be used in information exchange and storage in ways that are very different from the present electronic materials. Here we use resonant soft X-ray scattering spectroscopy to study the F-type (0,τ,0) and C-type (0, 1–2τ, 0) diffraction peaks from sinusoidal antiferromagnetic spin order in multiferroic DyMnO 3 . By comparing the temperature dependence of ordering wave vectors τ, peak intensities I, and correlation lengths λ measured at Mn L 2 - , O K-, and Dy M 5 -edges, we show that the nearly perfect locking between the ordering wave vectors from Dy 4f states and Mn 3d orbitals manifesting the second harmonic diffraction peak implies the notable orbital involvement in the coupling between Mn and Dy spins. Our DFT calculations further suggest that the lattice response to different antiferromagnetic ground states (A-type versus E-type) is much weaker in TbMnO 3 , in agreement with previous claim that the symmetric exchange interaction can be an important factor for understanding the ferroelectricity in DyMnO 3 than in TbMnO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bulk-like first-order magnetoelastic transition in FeRh particles

Near-equiatomic, chemically-ordered iron-rhodium (FeRh) alloy is a fundamentally interesting material that may become useful in niche applications making use of its unique magneto functional phenomena, for example, the giant inverse magnetocaloric effect near room temperature that is associated with a sharp first-order magnetic phase transition. The nearly discontinuous antiferromagnetic-ferromagnetic phase transformation in bulk FeRh is well-known; however, the transition broadens considerably in fine particles and films with thickness less than 50 nm, precluding their potential applications. Here, we report an abrupt, bulk-like first-order magnetoelastic transformation in powders consisting of sub-micron particles of nearly equiatomic FeRh compound synthesized via solid-state mechanochemical co-reduction of FeF 2 and RhCl3 and subsequent heat treatments. In this work, we demonstrate that annealing at temperatures ranging from 600 °C to 800 °C enables tailoring phase content, particle size, and magnetic properties of the powders. A maximum magnetic-field-induced entropy change of ~10 J/kg K at μ0ΔH = 1 T has been achieved in powders annealed at 800 °C. The retention of extraordinary responsiveness in sub-micron particles of FeRh is likely to open doors for system component fabrication using additive manufacturing methods, along with new opportunities to employ FeRh in theranostics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Implementing contact angle boundary conditions for second-order Phase-Field models of wall-bounded multiphase flows

In the present work, a general formulation is proposed to implement the contact angle boundary conditions for the second-order Phase-Field models, which is applicable to N-phase (N ≥ 2) moving contact line problems. To remedy the issue of mass change due to the contact angle boundary condition, a source term or Lagrange multiplier is added to the original second-order Phase-Field models, which is determined by the consistent and conservative volume distribution algorithm so that the summation of the order parameters and the consistency of reduction are not influenced. To physically couple the proposed formulation to the hydrodynamics, especially for large-density-ratio problems, the consistent formulation is employed. The reduction-consistent conservative Allen-Cahn models are chosen as examples to illustrate the application of the proposed formulation. The numerical scheme that preserves the consistency and conservation of the proposed formulation is employed to demonstrate its effectiveness. Results produced by the proposed formulation are in good agreement with the exact and/or asymptotic solutions. The proposed method captures complex dynamics of moving contact line problems having large density ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Improving the accuracy of discretisations of the vector transport equation on the lowest-order quadrilateral Raviart-Thomas finite elements

Within finite element models of fluids, vector-valued fields such as velocity or momentum variables are commonly discretised using the Raviart-Thomas elements. However, when using the lowest-order quadrilateral Raviart-Thomas elements, standard finite element discretisations of the vector transport equation typically have a low order of spatial accuracy. This paper describes two schemes that improve the accuracy of transporting such vector-valued fields on two-dimensional curved manifolds. The first scheme that is presented reconstructs the transported field in a higher-order function space, where the transport equation is then solved. The second scheme applies a mixed finite element formulation to the vector transport equation, simultaneously solving for the transported field and its vorticity. In this work, an approach to stabilising this mixed vector-vorticity formulation is presented that uses a Streamline Upwind Petrov-Galerkin (SUPG) method. These schemes are then demonstrated, along with their accuracy properties, through some numerical tests. Two new test cases are used to assess the transport of vector-valued fields on curved manifolds, solving the vector transport equation in isolation. The improvement of the schemes is also shown through two standard test cases for rotating shallow-water models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

In a computational fluid model of the atmosphere, the advective transport of trace species, or tracers, can be computationally expensive. For efficiency, models often use semi-Lagrangian advection methods. High-order interpolation semi-Lagrangian (ISL) methods, in particular, can be extremely efficient, if the problem of property preservation specific to them can be addressed. Atmosphere models often use geometrically and logically nonuniform grids for efficiency and, as a result, element-based discretizations. Such grids and discretizations make stability a particular problem for ISL methods. Generally, high-order, element-based ISL methods that use the natural polynomial interpolant associated with a nodal finite-element discretization are unstable. Here, we derive new bases having order of accuracy up to nine, with positive nodal weights, that stabilize the element-based ISL method. We use these bases to construct the linear advection operator in the property-preserving Interpolation Semi-Lagrangian Element-based Transport (Islet) method. Then we discuss key software implementation details. Finally, we show performance results for the Energy Exascale Earth System Model's atmosphere dynamical core, comparing the original and new transport methods. These simulations used up to 27,600 Graphical Processing Units (GPU) on the Oak Ridge Leadership Computing Facility's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Efficient data-driven regression for reduced-order modeling of spatial pattern formation

We present an efficient data-driven regression approach for constructing reduced-order models (ROMs) of reaction-diffusion systems exhibiting pattern formation. The ROMs are learned non-intrusively from available training data of physically accurate numerical simulations. The method can be applied to general nonlinear systems through the use of polynomial model form, while not requiring knowledge of the underlying physical model, governing equations, or numerical solvers. The process of learning ROMs is posed as a low-cost least-squares problem in a reduced-order subspace identified via Proper Orthogonal Decomposition (POD). Numerical experiments on classical pattern-forming systems–including the Schnakenberg and Mimura–Tsujikawa models–demonstrate that higher-order surrogate models significantly improve prediction accuracy while maintaining low computational cost. The proposed method provides a flexible, non-intrusive model reduction framework, well suited for the analysis of complex spatio-temporal pattern formation phenomena.

Data-driven modeling↗
DOE OSTI2023

Towards the design of nature-inspired materials: Impact of complex pore morphologies via higher-order homogenization

Even though the development of novel materials that mimic nature is widely used in a variety of engineering and scientific fields, the relationship between effective material properties and underlying, often complex pore morphology is still not fully understood. To address this knowledge gap and accelerate the development of novel nature-inspired materials, this paper adopts a higher-order asymptotic homogenization method to numerically investigate the effect of complex micropore morphology on the effective mechanical properties of a porous system. Specifically, we create unique pore morphologies with varying levels of complexity that serve as a more realistic representation of natural materials. Here, we then use the second-order homogenization method to capture the role of pore size, shape, orientation, and distribution on effective properties. By creating different pore morphologies, we systematically studied the relationship between morphology and effective mechanical properties. The results highlight the necessity of higher-order parameters to fully capture the role of realistic pore morphologies on effective mechanical properties and provide a path forward in the design of nature-inspired materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Identifying point defects and ordering in the high-entropy layered oxide Li 1.5 MO 3-δ (M=Mn, Al, Fe, Co, Ni) for energy storage applications

High-entropy layered oxides (HELOs) represent a very promising class of next-generation battery cathodes, combining the well-studied properties of layered cathode materials such as LiCoO 2 with the chemical tunability and stability of high-entropy materials. HELO materials often form particles with complex defects and domain structures, complicating accurate characterization of structure and cation mixing. Understanding disorder and order in HELO materials is necessary for understanding their performance and utility as cathodes. Here we demonstrate the characterization of the HELO Li 1.5 MO 3-δ (M = Mn, Al, Fe, Co, Ni), wherein X-ray powder diffraction, transmission electron microscopy imaging, and electron diffraction patterns are analyzed to reveal the presence of ordering in the HELO, with imaging and diffraction simulation employed to compare experimental results to atomic modeling. Finally, without rigorous characterization at the atomic scale, important features such as defect ordering can be easily overlooked and therefore remain unconsidered when interpreting experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Closing in on critical net-baryon fluctuations at LHC energies: Cumulants up to third order in Pb–Pb collisions

Fluctuation measurements are important sources of information on the mechanism of particle production at LHC energies. This article reports the first experimental results on third-order cumulants of the net-proton distributions in Pb–Pb collisions at a center-of-mass energy $\sqrt{s_{NN}}$ = 5.02 TeV recorded by the ALICE detector. The results on the second-order cumulants of net-proton distributions at $\sqrt{s_{NN}}$ = 2.76 and 5.02 TeV are also discussed in view of effects due to the global and local baryon number conservation. The results demonstrate the presence of long-range rapidity correlations between protons and antiprotons. Such correlations originate from the early phase of the collision. The experimental results are compared with HIJING and EPOS model calculations, and the dependence of the fluctuation measurements on the phase-space coverage is examined in the context of lattice quantum chromodynamics (LQCD) and hadron resonance gas (HRG) model estimations. The measured third-order cumulants are consistent with zero within experimental uncertainties of about 4% and are described well by LQCD and HRG predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗