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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2025

Ab initio computations of the fourth-order charge density moments of 48 Ca and 208 Pb

Neutron skins of neutron-rich nuclei connect nuclei with the matter in neutron stars. High-precision measurements of nuclear charge densities to extract higher-order moments are proposed to be sensitive to neutron radii and skin thicknesses. We investigate the charge density of 48 Ca and 208 Pb, leading candidates for such studies, with ab initio nuclear structure calculations. We find strong correlations between the fourth-order charge density moment R$^{4}_{ch}$ and the charge and neutron radii, allowing us to predict R$^{4}_{ch}$ for 48 Ca and 208 Pb. We find a substantially weaker correlation between the fourth-order charge density moment and the neutron skin, limiting the ability of high-precision electron scattering to determine the neutron skin in a model-independent manner.

Ab initio nuclear theory↗
DOE OSTI2026

Long-range magnetic order and relaxor ferroelectricity in a hexagonal high-entropy ferrite

Multiferroics that combine ferroelectricity and magnetic order are attractive for electronic and spintronic technologies, yet chemical disorder that promotes relaxor ferroelectricity usually suppresses long-range magnetic order. Here, we report entropy-stabilized relaxor multiferroicity in epitaxial hexagonal (Tb0.2Dy0.2Ho0.2Lu0.2Yb0.2)FeO3 thin films. Structural, magnetic, dielectric, and synchrotron spectroscopic measurements show the coexistence of relaxor ferroelectricity and long-range ferromagnetic order. We find that improper ferroelectricity remains robust against A-site configurational disorder, while the Fe sublattice preserves magnetic exchange. This separation of the microscopic origins of the polar and magnetic responses enables chemically disordered multiferroicity. Our results establish entropy engineering in hexagonal ferrites as a route toward multifunctional oxide thin films and provide a general design strategy for high-entropy multiferroics.

Miertschin, Duncan [Baylor University]↗
DOE OSTI2023

Genome-scale phylogeny and comparative genomics of the fungal order Sordariales

The order Sordariales is taxonomically diverse, and harbours many species with different lifestyles and large economic importance. Despite its importance, a robust genome-scale phylogeny, and associated comparative genomic analysis of the order is lacking. In this study, we examined whole-genome data from 99 Sordariales, including 52 newly sequenced genomes, and seven outgroup taxa. We inferred a comprehensive phylogeny that resolved several contentious relationships amongst families in the order, and cleared-up intrafamily relationships within the Podosporaceae. Extensive comparative genomics showed that genomes from the three largest families in the dataset (Chaetomiaceae, Podosporaceae and Sordariaceae) differ greatly in GC content, genome size, gene number, repeat percentage, evolutionary rate, and genome content affected by repeat-induced point mutations (RIP). All genomic traits showed phylogenetic signal, and ancestral state reconstruction revealed that the variation of the properties stems primarily from within-family evolution. Together, the results provide a thorough framework for understanding genome evolution in this important group of fungi.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Role of Short-Range Order on Diffusion Coefficients in the Li–Mg Alloy

Li–Mg alloys are important because of their beneficial role in fostering uniform plating and stripping of lithium in all-solid-state batteries. The alloy Li x Mg 1–x forms a solid solution on the BCC crystal structure when the lithium content is greater than x ≈ 0.3. The activation barriers of lithium and magnesium exchanges with a vacancy, crucial for substitutional diffusion, are predicted to be exceptionally low and almost identical, with negligible dependence on the alloy composition. The equilibrium vacancy concentration at room temperature is predicted to be very low, and it also remains almost constant with no dependence on Mg content in the alloy (for x Li ≥ 0.5). Nevertheless, both experiments and kinetic Monte Carlo simulations indicate that the tracer diffusion coefficients decrease by almost an order of magnitude with the addition of Mg to the alloy. In this contribution, the crucial role that chemical short-range order plays in affecting the diffusion coefficients is studied. Chemical short-range order is found to increase the effective activation barrier for lithium and magnesium diffusion by making successive atomic hops with vacancies correlated.

36 MATERIALS SCIENCE↗