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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A-type antiferromagnetic order in MnBi 4 Te 7 and MnBi 6 Te 10 single crystals

MnBi 4 Te 7 and MnBi 6 Te 10 are two members with n = 2 and 3 in the family of MnBi 2 n Te 3 n + 1 where the n = 1 member, MnBi 2 Te 4 , has been intensively investigated as the first intrinsic antiferromagnetic topological insulator. In this work, we report the A-type antiferromagnetic order in these two compounds by measuring magnetic properties, electrical and thermal transport, specific heat, and single-crystal neutron diffraction. Both compounds order into an A-type antiferromagnetic structure as does MnBi 2 Te 4 with ferromagnetic planes coupled antiferromagnetically along the crystallographic c axis. While no evidence for any in-plane ordered moment is found for MnBi 2 Te 4 or MnBi 6 Te 10 , weak reflections at half- L positions along the [0 0 L ] direction are observed for MnBi 4 Te 7 suggesting an in-plane ordered moment around 0.15 μ B / Mn . The ordering temperature, T N , is 13 K for MnBi 4 Te 7 and 11 K for MnBi 6 Te 10 . The magnetic order is also manifested in the anisotropic magnetic properties. Furthermore, for both compounds, the interlayer coupling is weak and a spin-flip transition occurs when a magnetic field of around 1.6 kOe is applied along the c axis at 2 K. As observed in MnBi 2 Te 4 , when cooling across T N , no anomaly was observed in the temperature dependence of thermopower. On the other hand, critical scattering effects are observed in thermal conductivity although the effect is less pronounced than that in MnBi 2 Te 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Unconventional magnetic order emerging from competing energy scales in the new R Rh 3 Si 7 intermetallics ( R = Gd-Yb)

The competition between Ruderman-Kittel-Kasuya-Yosida (RKKY), crystal electric field (CEF), and Kondo energy scales has recently emerged at the heart of complex magnetism in several Ce- or Yb-based intermetallics. Hard axis magnetic order has been observed in a handful of these compounds, independent of the crystal symmetry, size of the ordered moment, or the relative scale of the Kondo and magnetic ordering temperatures. This raises the question of the role of each energy scale in driving the ground state properties. In focusing on a single class of compounds, the rhombohedral RRh 3 Si 7 , we compare the anisotropy and magnetic ground states in members of this series with only RKKY interactions (R = Gd), or RKKY and CEF effects (R = Tb-Tm), with the behavior of the R = Yb compound, where all three energy scales (RKKY, CEF, Kondo) are at play. Moreover, we extend the comparison to two other isostructural Kondo systems YbIr 3 Si 7 and YbIr 3 Ge 7 , where hard axis magnetic order is also observed. The non-Kondo compounds RRh 3 Si 7 (R = Tb-Tm) lack the complexity of magnetic order along the hard CEF axis, pointing to the dominant role of the Kondo effect in driving this magnetic order. Furthermore, the CEF-RKKY competition is still responsible for complex magnetic ground states, and it appears that the electronic and magnetic degrees of freedom are entangled in all magnetic members of this series of compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2020

One-component order parameter in URu 2 Si 2 uncovered by resonant ultrasound spectroscopy and machine learning

The unusual correlated state that emerges in URu 2 Si 2 below T HO = 17.5 K is known as “hidden order” because even basic characteristics of the order parameter, such as its dimensionality (whether it has one component or two), are “hidden.” We use resonant ultrasound spectroscopy to measure the symmetry-resolved elastic anomalies across T HO . We observe no anomalies in the shear elastic moduli, providing strong thermodynamic evidence for a one-component order parameter. We develop a machine learning framework that reaches this conclusion directly from the raw data, even in a crystal that is too small for traditional resonant ultrasound. Our result rules out a broad class of theories of hidden order based on two-component order parameters, and constrains the nature of the fluctuations from which unconventional superconductivity emerges at lower temperature. Our machine learning framework is a powerful new tool for classifying the ubiquitous competing orders in correlated electron systems.

42 ENGINEERING↗
DOE OSTI2021

Elastic Wave Propagation in Curvilinear Coordinates with Mesh Refinement Interfaces by a Fourth Order Finite Difference Method

In this work, we develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinement interfaces with hanging nodes, physical interface conditions are imposed by using ghost points and interpolation. With a fourth order predictor-corrector time integrator, the fully discrete scheme is energy conserving. Numerical experiments are presented to verify the fourth order convergence rate and the energy conserving property.

58 GEOSCIENCES↗
DOE OSTI2020

Smoothers for Matrix-Free Algebraic Multigrid Preconditioning of High-Order Finite Elements

We investigate smoothers for use in matrix-free algebraic multigrid (AMG) preconditioning of high-order finite element problems. These AMG preconditioners are matrix-free in the sense that they are built from a related low-order refined finite element problem whose system matrix can be much more rapidly assembled than the high-order problem. Our proposed smoother, which we call distributive relaxation, is more robust to the anisotropy present in many low-order refined meshes which feature a clustering of nodes near the boundaries between high-order finite elements. For solving the low-order refined problem, we show that this new distributive relaxation smoother possesses significantly improved performance compared to more traditional smoothers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

High-Order Wall-Modeled Large-Eddy Simulation of High-Lift Configuration

This paper presents the assessment of several recent enhancements for a high-order wall-modeled large-eddy simulation (WMLES) approach and demonstrates order independence with a fixed data exchange location in the wall model. The two enhancements include the use of isotropic tetrahedral elements to improve accuracy and an explicit subgrid-scale model, the Vreman model, to improve accuracy and robustness. The [Formula: see text] study focused on the high-lift Common Research Model (HL-CRM) at the angle of attack of 19.57 deg, a benchmark problem from the 4th AIAA High-Lift Prediction Workshop. Solution polynomial orders of [Formula: see text], and 5 were used in the study. The study demonstrated [Formula: see text] independence in integrated forces, pitch moment, velocity profile in the wall-normal direction, and surface flow topology. It also showed that a [Formula: see text] order of at least 3 ([Formula: see text]) was needed to correctly predict the external inviscid flow and the surface flow topology. Thereafter, [Formula: see text] simulations over several other angles of attack demonstrated that the high-order WMLES approach can correctly predict the maximum lift and flow separation regions for HL-CRM with about 40 million degrees of freedom (DOF) compared to at least 250 million DOF required by second-order methods.

Engineering↗
DOE OSTI2025

Toward Higher-order Accuracy in Self-gravitating Hydrodynamics

High-order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution Δ x ) with only a moderate increase to computational expense. Significant effort has been placed in the development of finite-volume algorithms for (magneto)hydrodynamics; however, state-of-the-art astrophysical simulations tightly couple a plenitude of physics, additionally including gravity, photon transport, cosmic-ray transport, chemistry, and/or diffusion, to name a few. Algorithms frequently operator-split this additional physics (often a first-order error in time) and/or adopt a model wherein their evaluation is limited to second-order accuracy in space. In this work, we present a fourth-order-accurate finite-volume scheme for self-gravitating hydrodynamics on a uniform Cartesian grid. The method supplies source terms for the gravitational acceleration ( ρ g ) and gravitational energy release ( ρ v · g ) associated with fourth-order-accurate solutions to the Poisson equation. Our scheme (1) guarantees the conservation of total linear momentum while (2) decreasing (in proportion to Δ x 4 ) the effects of spurious heating and/or cooling associated with truncation error in the gravity. We demonstrate expected convergence rates for the algorithm by measuring errors in test problems evolving self-gravity modified linear waves and 3D polytropic equilibria. We test robustness of the algorithm by integrating an induced “inside-out” adiabatic collapse. We also discuss a method to smoothly downgrade the solution to second-order spatial accuracy to avoid spurious overshoots near steep density and/or pressure gradients.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Novel Charge-Ordered States in 1 T -IrTe 2

1T-IrTe 2 exhibits an intriguing series of charge-ordered states in both bulk and atomically thin layers. The charge orders in IrTe 2 emerge with consecutive first-order structural transitions involving in-plane Ir–Ir dimerization and interlayer Te–Te depolymerization, resulting in stripe patterns with various periods. Upon chemical doping, 1T-IrTe 2 can also be driven into a superconductor. Despite the intense research effort, a comprehensive understanding of the nature of the charge ordering and its relation to the superconductivity is yet to be reached. This review provides an overview of the novel charge orders and the emergence of superconductivity in IrTe 2 . Finally, we further introduce the recent studies on IrTe 2 nanoflake, in which novel charge orders attainable only in two-dimensional limits lay the grounds for an exotic phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A second order yield-temperature relation for accurate inference of burn-averaged quantities in multi-species plasmas

Measured yields and ion temperatures inferred from the fusion product energy spectra can be used as metrics for the performance of an ICF implosion. This can be to infer species separation, thermal decoupling, flows or other effects that can cause the inferred ion temperatures to deviate from the true underlying thermal temperature and the yield ratio to deviate from the expected value. Direct inference of the impact of these effects on observed temperatures and yields can be difficult to uncover due to underlying dependence on the shape and time evolution of the temperature and density proles of the fusing plasma. Due to differences in the temperature dependence of the reactivities, different fusion products are emitted from different regions and times within the implosion. In order to properly account for this, a second order analytic expression relating the apparent temperatures and yield ratios is developed. This expression can be coupled to models of yield and/or temperature altering effects to infer their burn-averaged impact on an implosion. The second order expression shows significant improvement over lower order expressions in synthetic data studies. Demonstrations of its applications to synthetic data coupled with models of ion thermal decoupling and radial flows are presented. In the case of thermal decoupling both first and second order expressions show reasonable levels of accuracy. To consistently infer the amplitude of radial flow with <10% error the second order equation is required.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗