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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
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 OSTI2015

A POD reduced order model for resolving angular direction in neutron/photon transport problems

This article presents the first Reduced Order Model (ROM) that efficiently resolves the angular dimension of the time independent, mono-energetic Boltzmann Transport Equation (BTE). It is based on Proper Orthogonal Decomposition (POD) and uses the method of snapshots to form optimal basis functions for resolving the direction of particle travel in neutron/photon transport problems. A unique element of this work is that the snapshots are formed from the vector of angular coefficients relating to a high resolution expansion of the BTE's angular dimension. In addition, the individual snapshots are not recorded through time, as in standard POD, but instead they are recorded through space. In essence this work swaps the roles of the dimensions space and time in standard POD methods, with angle and space respectively. It is shown here how the POD model can be formed from the POD basis functions in a highly efficient manner. The model is then applied to two radiation problems; one involving the transport of radiation through a shield and the other through an infinite array of pins. Both problems are selected for their complex angular flux solutions in order to provide an appropriate demonstration of the model's capabilities. It is shown that the POD model can resolve these fluxes efficiently and accurately. In comparison to high resolution models this POD model can reduce the size of a problem by up to two orders of magnitude without compromising accuracy. Solving times are also reduced by similar factors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Linearly first- and second-order, unconditionally energy stable schemes for the phase field crystal model

In this paper, we develop a series of linear, unconditionally energy stable numerical schemes for solving the classical phase field crystal model. The temporal discretizations are based on the first order Euler method, the second order backward differentiation formulas (BDF2) and the second order Crank–Nicolson method, respectively. The schemes lead to linear elliptic equations to be solved at each time step, and the induced linear systems are symmetric positive definite. We prove that all three schemes are unconditionally energy stable rigorously. Various classical numerical experiments in 2D and 3D are performed to validate the accuracy and efficiency of the proposed schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes–Brinkman–Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes–Brinkman–Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies the effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2016

Dependence of short and intermediate-range order on preparation in experimental and modeled pure a-Si

Variability in the short-to-intermediate range order of pure amorphous silicon prepared by different experimental and computational techniques is probed by measuring mass density, atomic coordination, bond-angle deviation, and dihedral angle deviation. It is found that there is significant variability in order parameters at these length scales in this archetypal covalently bonded, monoatomic system. This diversity strongly reflects preparation technique and thermal history in both experimental and simulated systems. Experiment and simulation do not fully quantitatively agree, partly due to differences in the way parameters are accessed. However, qualitative agreement in the trends is identified. Relaxed forms of amorphous silicon closely resemble continuous random networks generated by a hybrid method of bond-switching Monte Carlo and molecular dynamics simulation. As-prepared ion implanted amorphous silicon can be adequately modeled using a structure generated from amorphization via ion bombardement using energetic recoils. Preparation methods which narrowly avoid crystallization such as experimental pressure-induced amorphization or simulated melt-quenching result in inhomogeneous structures that contain regions with significant variations in atomic ordering. Ad hoc simulated structures containing small (1 nm) diamond cubic crystal inclusions were found to possess relatively high bond-angle deviations and low dihedral angle deviations, a trend that could not be reconciled with any experimental material.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Highly conducting divalent Mg{sup 2+} cation solid electrolytes with well-ordered three-dimensional network structure

A three-dimensionally well-ordered NASICON-type Mg{sup 2+} cation conductor, (Mg{sub x}Hf{sub 1−x}){sub 4/(4−2x)}Nb(PO{sub 4}){sub 3}, was firstly developed by partial substitution of lower valent Mg{sup 2+} cation onto the Hf{sup 4+} sites in a HfNb(PO{sub 4}){sub 3} solid to realize high Mg{sup 2+} cation conductivity even at moderate temperatures. Due to the formation of well-ordered NASICON-type structure, both the high Mg{sup 2+} cation conductivity below 450 °C and the low activation energy for Mg{sup 2+} cation migration was successfully realized for the (Mg{sub 0.1}Hf{sub 0.9}){sub 4/3.8}Nb(PO{sub 4}){sub 3} solid. Pure Mg{sup 2+} cation conduction in the NASICON-type (Mg{sub 0.1}Hf{sub 0.9}){sub 4/3.8}Nb(PO{sub 4}){sub 3} solid was directly and quantitatively demonstrated by means of two kinds of dc electrolysis. - Graphical abstract: Image of the Mg{sup 2+} cation conduction in NASICON-type (Mg{sub 0.1}Hf{sub 0.9}){sub 4/3.8}Nb(PO{sub 4}){sub 3} and its Mg{sup 2+} conductivity. - Highlights: • We develop a three-dimensionally well-ordered NASICON-type Mg{sup 2+} cation conductor. • A high magnesium cation conductivity is realized even at moderate temperatures. • Divalent magnesium cation conduction is demonstrated directly and quantitatively.

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

Oriented and ordered mesoporous ZrO{sub 2}/TiO{sub 2} fibers with well-organized linear and spring structure

Graphical abstract: The ultra-stable order mesoporous ZrO{sub 2}/TiO{sub 2} fibers with well-organized linear and spring structure and large surface area under higher temperatures were prepared by a simple EISA process. - Highlights: • The ZrO{sub 2}/TiO{sub 2} fibers were prepared by EISA process combined with steam heat-treatment. • The mesoporous ZrO{sub 2}/TiO{sub 2} fibers have well-organized linear and spring structure. • The fibers were composed of oval rod nanocrystals of ZrTiO{sub 4}. - Abstract: The ultra-stable order mesoporous ZrO{sub 2}/TiO{sub 2} fibers with well-organized linear and spring structure and large surface areas under higher temperatures were prepared by a (simple evaporation-induced assembly) EISA process. The preparation, microstructures and formation processes were characterized by Fourier transformation infrared (FTIR), X-ray diffraction (XRD), scanning electron microscopy (SEM), transmission electron microscopy (TEM) and N{sub 2} adsorption–absorption measurements. The fibers take on pinstripe configuration which is very orderly along or perpendicular to the axial direction of the fibers. The diameters of the pinstripe are in the region of 200–400 nm and arranges regularly, which are composed of oval rod nanocrystals of ZrTiO{sub 4}.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Structural, optical and femtosecond third-order nonlinear optical properties of LiNbO{sub 3} thin films

Highlights: • Fs third order NLO properties of LiNbO{sub 3} thin films were investigated. • Al:ZnO was used as a buffer layer to control the growth of LiNbO{sub 3}. • The magnitude of the χ{sup (3)} of all films was ∼10{sup −11} e.s.u. • χ{sup (3)} obtained here is much higher that of single crystal Z-cut LiNbO{sub 3}. - Abstract: LiNbO{sub 3} thin films were deposited by means of Al:ZnO as intermediate layer using pulsed laser deposition and their structural, linear and nonlinear optical (NLO) properties were investigated. The films were found to be polycrystalline, (110) orientation for Al:ZnO/fused silica and (300) orientation for Al:ZnO/quartz substrates. The films demonstrated a transparency of >80% in the visible and near-infrared region. The NLO studies were performed with the Z-scan technique using ∼150 fs, 80 MHz pulses at 800 nm. The studies revealed that the polycrystalline LiNbO{sub 3}films and (300) LiNbO{sub 3}films demonstrated reverse saturable absorption while a transition to saturable absorption was observed in the case of (110) LiNbO{sub 3}. The presence of two-photon absorption is attributed to the indirect transition via intermediate defect levels which is related to the structural asymmetry. The large nonlinearity in refractive index could be a result of the strain induced ferroelectric polarization. The magnitude of the third order optical nonlinearity of all films was ∼10{sup −11} e.s.u. which is orders of magnitude higher than the value obtained for single crystal Z-cut LiNbO{sub 3}.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Room temperature magnetic ordering and analysis by bound magnetic polaron model of Yb3+ doped nanocrystalline zinc oxide (Zn0.98Yb0.02O)

Highlights: • ZnO and Yb doped ZnO are synthesized by co-precipitation method. • Rietvelt anlysis of XRD patterns confirms the complete substitution of Yb{sup 3+} in ZnO. • Vacancy and defect mediated RTFM is introduced in doped ZnO. • BMP model successfully analyze the magnetic property of Zn{sub 0.98}Yb{sub 0.02}O. - Abstract: Nanocrystalline ZnO and Yb-doped ZnO (Zn{sub 0.98}Yb{sub 0.02}O) are prepared by simple co-precipitation method. The pure crystallographic phases of ZnO and Zn{sub 0.98}Yb{sub 0.02}O are confirmed by Rietveld analyses of XRD patterns. Raman and photoluminescence spectra recorded at room temperature (RT) indicates that lattice defects and oxygen vacancy in Zn{sub 0.98}Yb{sub 0.02}O is more than that of ZnO. Magnetic measurements of ZnO and Zn{sub 0.98}Yb{sub 0.02}O are carried out from 300 to 5 K. Ferromagnetic ordering is observed in Zn{sub 0.98}Yb{sub 0.02}O at RT and the strength of the ordering increases with the lowering of temperature. Ferromagnetism of Zn{sub 0.98}Yb{sub 0.02}O is attributed to the oxygen vacancy introduced by annealing in vacuum. Magnetization vs. temperature curve is well fitted by the combination of 3D spin wave and Curie-Weiss law, which suggest simultaneous presence of para- and ferro-magnetic phases in Zn{sub 0.98}Yb{sub 0.02}O. Ferromagnetic ordering is successfully explained by vacancy mediated bound-magnetic polaron model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2016

Collapse for the higher-order nonlinear Schrödinger equation

We examine conditions for finite-time collapse of the solutions of the higher-order nonlinear Schr odinger (NLS) equation incorporating third-order dispersion, self-steepening, linear and nonlinear gain and loss, and Raman scattering; this is a system that appears in many physical contexts as a more realistic generalization of the integrable NLS. By using energy arguments, it is found that the collapse dynamics is chiefly controlled by the linear/nonlinear gain/loss strengths. We identify a critical value of the linear gain, separating the possible decay of solutions to the trivial zero-state, from collapse. The numerical simulations, performed for a wide class of initial data, are found to be in very good agreement with the analytical results, and reveal long-time stability properties of localized solutions. The role of the higher-order effects to the transient dynamics is also revealed in these simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗