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

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

High-order multirate explicit time-stepping schemes for the baroclinic-barotropic split dynamics in primitive equations

In order to treat the multiple time scales of ocean dynamics in an efficient manner, the baroclinic-barotropic splitting technique has been widely used for solving the primitive equations for ocean modeling. Based on the framework of strong stability-preserving Runge-Kutta approach, we propose two high-order multirate explicit time-stepping schemes (SSPRK2-SE and SSPRK3-SE) for the resulting split system in this paper. The proposed schemes allow for a large time step to be used for the three-dimensional baroclinic (slow) mode and a small time step for the two-dimensional barotropic (fast) mode, in which each of the two mode solves just need to satisfy their respective CFL conditions for numerical stability. Specifically, at each time step, the baroclinic velocity is first computed by advancing the baroclinic mode and fluid thickness of the system with the large time-step and the assistance of some intermediate approximations of the barotropicmode obtained by substepping with the small time step; then the barotropic velocity is corrected by using the small time step to re-advance the barotropic mode un-der an improved barotropic forcing produced by interpolation of the forcing terms from the preceding baroclinic mode solves; lastly, the fluid thickness is updated by coupling the baroclinic and barotropic velocities. Additionally, numerical inconsistencies on the discretized sea surface height caused by the mode splitting are relieved via a reconciliation process with carefully calculated flux deficits. Here, two benchmark tests from the “MPAS-Ocean” platform are carried out to numerically demonstrate the performance and parallel scalability of the proposed SSPRK-SE schemes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Accelerating high-order mesh optimization using finite element partial assembly on GPUs

In this paper we present a new GPU-oriented mesh optimization method based on high order finite elements. Our approach relies on node movement with fixed topology, through the Target-Matrix Optimization Paradigm (TMOP) and uses a global nonlinear solve over the whole computational mesh, i.e., all mesh nodes are moved together. A key property of the method is that the mesh optimization process is recast in terms of finite element operations, which allows us to utilize recent advances in the field of GPU-accelerated high order finite element algorithms. For example, we reduce data motion by using tensor factorization and matrix-free methods, which have superior performance characteristics compared to traditional full finite element matrix assembly and offer advantages for GPU based HPC hardware. Furthermore, we describe the major mathematical components of the method along with their efficient GPU-oriented implementation. In addition, we propose an easily reproducible mesh optimization test that can serve as a performance benchmark for the mesh optimization community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Metamagnetism and kinetic arrest in a long-range ferromagnetically ordered multicaloric double perovskite Y 2 CoMnO 6

A systematic magnetic study of the double perovskite oxide Y 2 CoMnO 6 (YCMO) has been performed. A monoclinic P21/n phase of YCMO was synthesized using a sol-gel method. Neutron diffraction (ND) measurements evidence the onset of long-range ferromagnetic (FM) ordering at TC ~ 76 K, which persists down to 5 K. The presence of 25% antisite disorder, estimated from the ND data, leads to the appearance of short-range antiferromagnetic (AFM) interactions. The existence of thermal hysteresis due to competing interactions between FM and AFM phases is observed in the M vs. T measurements. The pinning of magnetic domain walls at the Co/Mn antiphase boundaries results in a metamagnetic-like behavior. The field dependence of thermomagnetic irreversibility and the nature of virgin curves indicate the occurrence of a kinetic arrest phenomenon, which is further verified via cooling and heating of the system in an unequal fields (CHUF) protocol. In full agreement with the ND data, a detailed analysis of critical exponents near the paramagnetic (PM)-FM phase transition also establishes YCMO as a long-range interacting mean-field system. A close examination of the temperature- and field-dependent magnetic entropy change yields an in-depth understanding of coexisting magnetically ordered and disordered phases in YCMO, leading to a comprehensive magnetic phase diagram of this multifunctional double perovskite system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Neutron diffraction study of magnetic ordering in high pressure phases of rare earth metal holmium

The magnetic ordering in rare earth metals is well established for ambient pressure crystal structures, however, little is known about the magnetic ordering in their corresponding high-pressure crystalline modifications. Holmium (Ho) was studied in a large-volume diamond anvil cell at the Spallation Neutron Source to high-pressure up to 20 GPa and to low-temperature to 10 K. We have conducted two independent high-pressure low-temperature experiments under non-hydrostatic and quasi-hydrostatic pressure conditions respectively. The ambient pressure hexagonal close packed (hcp) phase of holmium shows two magnetic transitions below 10 GPa one to an incommensurate Antiferromagnetic (AFM) phase and another to a conical-Ferromagnetic (c-FM) phase. In contrast, alpha-Samarium-type (α-Sm) phase above 10 GPa and the double hexagonal close packed (dhcp) phase above 19 GPa show only one FM transition marked by the appearance of a magnetic peak at 3 Å and the concurrent enhancement of nuclear peaks below 30 K. Finally, these new transitions observed by neutron diffraction can be accounted by a commensurate superlattice formation along c-axis in both the α-Sm-type phase and the dhcp phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Magnetic ordering in GdAuAl 4 Ge 2 and TbAuAl 4 Ge 2 : Layered compounds with triangular lanthanide nets

In this work, we report the synthesis and magnetic properties of GdAuAl 4 Ge 2 and TbAuAl 4 Ge 2 , where temperature and magnetic field dependent magnetization, heat capacity, and electrical resistivity measurements reveal that both compounds exhibit several magnetically ordered states at low temperatures, with evidence for magnetic fluctuations extending into the paramagnetic temperature region. For magnetic fields applied in the ab-plane there is particularly rich behavior, with several ordered state regions that are separated by metamagnetic phase transitions. Despite Gd being an isotropic S-state ion and Tb having an anisotropic J-state, there are similarities in the phase diagrams for the two compounds, suggesting that factors such as the symmetry of the crystalline lattice, which features well separated triangular planes of lanthanide ions, or the Ruderman–Kittel–Kasuya– Yosida interaction control the magnetism. We also point out similarities to other centrosymmetric compounds that host skyrmion lattices such as Gd 2 PdSi 3 , and propose that the LnAuAl 4 Ge 2 family of compounds are of interest as reservoirs for complex magnetism and electronic behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗
DOE OSTI2023

Search for octupolar order in NpO 2 by neutron powder diffraction

The magnetic structure of neptunium dioxide (NpO 2 ) remains a mystery despite decades of research. We revisited the search for magnetic ordering in NpO 2 by performing a powder neutron diffraction experiment at low temperatures and relatively large values of Q. Diffraction data were collected for 300 MK, 15 K, and 35 K to understand the phase transition near 25 K. However, no significant changes in the neutron diffraction pattern were identified. Using the octupolar magnetic form factor j 6 , the maximum theoretical expected value of the magnetic moment (3.27 μ B ), the calculated magnetic Bragg peak intensities would not be observable in neutron powder diffraction data. This value is much larger than the previously estimated magnetic moment of ~0.1 μ B . In conclusion, supported by experimental results, these calculations suggest that unpolarized neutron powder diffraction is unsuitable for measuring octupolar ordering in NpO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Architected mesoporous crystalline magnesium silicates with ordered pore structures

Novel approaches to harness earth abundant silicates as building blocks for carbon dioxide removal, capture, utilization, and storage are gaining increasing attention in the context of sustainable and low carbon energy and resource recovery. Advancing a calibrated understanding of these fluid-silicate interactions is essential for developing scalable processes. One of the challenges in developing predictive controls over these interactions is the compositional and morphological heterogeneity of naturally occurring, heterogeneous magnesium silicate minerals. To address this challenge, the synthesis of architected mesoporous crystalline magnesium silicate (Mg 2 SiO 4 ) is proposed. While synthesis routes for producing amorphous mesoporous magnesium silicates have been developed via sol-gel methods, approaches to synthesize crystalline magnesium silicates with well-controlled pore size distributions have not been explored. The conventional approaches of converting matter that is amorphous to crystalline states at elevated temperatures results in a heterogeneous pore size distribution. To develop controls on pore size distribution, amorphous mesoporous magnesium silicates are coated with carbon. This approach preserves the pore size distributions during the amorphous to crystalline transformations of Mg-silicates at elevated temperatures. The carbon coating is removed on heating. Magnesium silicate particles produced using this approach have highly ordered pores around 2.58 nm and a specific surface area of 124.25 m 2 /g. In this study, we report the chemical compositions, morphologies, phase transitions, and pore structures of the intermediate materials produced during the synthesis of crystalline mesoporous magnesium silicates. The transitions, and pore structures of the intermediate materials produced during the synthesis of crystalline mesoporous magnesium silicates. Furthermore, the synthesis routes discussed in this study can be applied translationally to produce metal silicates with ordered mesoporous structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multi-modal imaging of a single mouse brain over five orders of magnitude of resolution

Mammalian neurons operate at length scales spanning six orders of magnitude; they project millimeters to centimeters across brain regions, are composed of micrometer-scale-diameter myelinated axons, and ultimately form nanometer scale synapses. Capturing these anatomical features across that breadth of scale has required imaging samples with multiple independent imaging modalities. Translating between the different modalities, however, requires imaging the same brain with each. Here, we imaged the same postmortem mouse brain over five orders of spatial resolution using MRI, whole brain micrometer-scale synchrotron x-ray tomography ($\mu$CT), and large volume automated serial electron microscopy. Using this pipeline, we can track individual myelinated axons previously relegated to axon bundles in diffusion tensor MRI or arbitrarily trace neurons and their processes brain-wide and identify individual synapses on them. This pipeline provides both an unprecedented look across a single brain's multi-scaled organization as well as a vehicle for studying the brain's multi-scale pathologies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A long short-term memory embedding for hybrid uplifted reduced order models

In this paper, we introduce an uplifted reduced order modeling (UROM) approach through the integration of standard projection based methods with long short-term memory (LSTM) embedding. Our approach has three modeling layers or components. In the first layer, we utilize an intrusive projection approach to model dynamics represented by the largest modes. The second layer consists of an LSTM model to account for residuals beyond this truncation. This closure layer refers to the process of including the residual effect of the discarded modes into the dynamics of the largest scales. However, the feasibility of generating a low rank approximation tails off for higher Kolmogorov n -width systems due to the underlying nonlinear processes. The third uplifting layer, called super-resolution, addresses this limited representation issue by expanding the span into a larger number of modes utilizing the versatility of LSTM. Therefore, our model integrates a physics-based projection model with a memory embedded LSTM closure and an LSTM based super-resolution model. In several applications, we exploit the use of Grassmann manifold to construct UROM for unseen conditions. We performed numerical experiments by using the Burgers and Navier-Stokes equations with quadratic nonlinearity. Finally, our results show robustness of the proposed approach in building reduced order models for parameterized systems and confirm the improved trade-off between accuracy and efficiency.

42 ENGINEERING↗
DOE OSTI2020

Dilute Fermi gas at fourth order in effective field theory

Using effective field theory methods, we calculate for the first time the complete fourth-order term in the Fermi-momentum or k F a s expansion for the ground-state energy of a dilute Fermi gas. The convergence behavior of the expansion is examined for the case of spin one-half fermions and compared against quantum Monte-Carlo results, showing that the Fermi-momentum expansion is well-converged at this order for |k F a s | ≲ 0.5.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A simple model for short-range ordering kinetics in multi-principal element alloys

Short-range ordering (SRO) in multi-principal element alloys influences material properties such as strength and corrosion. While some degree of SRO is expected at equilibrium, predicting the kinetics of its formation is challenging. We present a simplified isothermal concentration-wave (CW) model to estimate an effective relaxation time of SRO formation. Estimates from the CW model agree to within a factor of five with relaxation times obtained from kinetic Monte Carlo (kMC) simulations when above the highest ordering instability temperature. Further, the advantage of the CW model is that it only requires mobility and thermodynamic parameters, which are readily obtained from alloy mobility databases and Metropolis Monte Carlo simulations, respectively. The simple parameterization of the CW model and its analytical nature makes it an attractive tool for the design of processing conditions to promote or suppress SRO in multicomponent alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ultrafast and selective gas transport through highly ordered black phosphorene nanochannels

Two-dimensional (2D) materials bring a great opportunity to fabricate molecular sieving membranes that can potentially break the permeability-selectivity trade-off. Although 2D laminar membranes with interlayer nanochannels as molecular sieving channel were widely studied, for most of reported 2D laminar membranes, it is of a great challenge to fabricate highly ordered interlayer nanochannels for mass transport. Herein, we report a novel kind of black phosphorene membrane which is made from the stacking of highly ordered 2D black phosphorene nanoflakes. The as prepared black phosphorene membrane shows H 2 permeance > 1000 GPU and H 2 /CO 2 selectivity > 100 for H 2 /CO 2 mixed gas, demonstrating an extremely high gas separation performance. The DFT calculation results demonstrate that the interlayer galleries in the black phosphorene membrane allow the H 2 passing through easily while block the other gases with bigger kinetic diameters, matching well with the experimental findings. In conclusion, the present results indicate that the interlayer galleries in the black phosphorene membrane can be applied as molecular sieving channels for gas separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Geometric control of emergent antiferromagnetic order in coupled artificial spin ices

Artificial spin ices (ASIs) composed of coupled nanomagnets offer the possibility to create designer geometrical frustration and manipulate inter-nanomagnet interactions. In particular, by using a dimer motif consisting of two strongly coupled single-domain nano magnets as a building block, we can control and realize intriguing antiferromagnetic physical states in which the magnetic charge is not conserved. Here we create a dimer kagome ASI system within which the antiferromagnetic order is controlled by tuning lattice geometry. Different disordered antiferromagnetic phases and the inter-nanomagnet interactions are visualized directly in demagnetized lattices with different lattice parameters. Monte Carlo simulations establish that the collective ground state consists of disordered antiferromagnetic dimers across the lattice when the intra-dimer interaction is dominant. However, for lattices governed by intra-triad-unit interactions, the ground state exhibits a longrange spin-ordered state in which the vertex magnetic charge is uniform across all triads.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a five-dimensional continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular two-dimensional velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. Here we describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma using first- and second-order elements with the pseudo-inverse method and comparing with the bilinear mapping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Enumeration as a Tool for Structure Solution: A Materials Genomic Approach to Solving the Cation-Ordered Structure of Na 3 V 2 (PO 4 ) 2 F 3

While powder diffraction methods are routinely utilized to optimize structural models for compounds whose crystal structures are known, the determination of unknown structures is far more challenging. When the unknown structure is large, structure solution can become a virtually intractable problem using standard structure solution methodologies, especially when the space group cannot be unambiguously resolved. One such system is the promising Na-ion battery cathode material Na 3 V 2 (PO 4 ) 2 F 3 whose high temperature and room temperature structures were previously solved, but whose more complex low-temperature structure could not be determined. Here, a novel materials genomic approach is demonstrated for the solution of the unknown 100 K structure of Na 3 V 2 (PO 4 ) 2 F 3 in which enumeration methods are first used to generate a large number (~3,000) of trial structures based on plausible orderings of Na ions and then automated Rietveld refinements are carried out to optimize each of these trial structures. Based on both the analysis of the ensemble of optimized trial structures and the density functional theory energy minimization of selected trial structures, the 100 K structure of Na 3 V 2 (PO 4 ) 2 F 3 is best described as belonging to the space group A2 1 am with unit cell dimensions of a = 9.01928(4), b = 27.1379(1), c = 10.73307(5). The 100 K unit cell has a large volume of 2627.07(2) Å 3 with Z = 12 and 33 independent crystallographic sites (9 Na, 3 V, 3 P, 12 O, and 6 F) that is 3x and 6x larger than the room- and high-temperature polymorphs of this phase, respectively. Finally, the novel methods described here will be generally applicable for the solution of the complex cation-ordered structures that commonly occur for battery materials.

36 MATERIALS SCIENCE↗