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,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

The quantum mechanics of high-order kinematic values

In this paper, within the framework of a unified mathematical model, the new formulation of quantum mechanics – quantum mechanics of higher order kinematic values – is proposed. In contrast to the well-known formulations of quantum mechanics in the new formulation, the wave function has not only a coordinate or momentum representation, but also representations through acceleration and accelerations of higher orders. The representations of the wave function in terms of higher order kinematic values make it possible to consider the dissipative systems problem and the variable entropy systems one within a single mathematical apparatus. The new formulation is not built phenomenologically, but from first principles and in a particular case contains the formulation of the de Broglie–Bohm ⟨wave-pilot⟩ quantum mechanics. The new formulation is based on the Vlasov equations chain. This approach allows us to obtain a natural connection between classical and quantum systems. Examples of such systems are considered in this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2018

Ultrahigh-order Maxwell solver with extreme scalability for electromagnetic PIC simulations of plasmas

The advent of massively parallel supercomputers, with their distributed-memory technology using many processing units, has favored the development of highly-scalable local low-order solvers at the expense of harder-to-scale global very high-order spectral methods. Indeed, FFT-based methods, which were very popular on shared memory computers, have been largely replaced by finite-difference (FD) methods for the solution of many problems, including plasmas simulations with electromagnetic Particle-In-Cell methods. For some problems, such as the modeling of so-called “plasma mirrors” for the generation of high-energy particles and ultra-short radiations, we have shown that the inaccuracies of standard FD-based PIC methods prevent the modeling on present supercomputers at sufficient accuracy. Here, we demonstrate here that a new method, based on the use of local FFTs, enables ultrahigh-order accuracy with unprecedented scalability, and thus for the first time the accurate modeling of plasma mirrors in 3D.

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

The nature of the first order isostructural transition in GdRhSn

We present structural, magnetic, thermal, magnetocaloric, and electrical transport properties of polycrystalline GdRhSn. Magnetization data show that it orders antiferromagnetically at TN = 16.2 K. The compound has the ZrNiAl type hexagonal crystal structure at room temperature and undergoes a first order iso-structural transition in the paramagnetic state at 245 K. The unit cell volume change at the transition is small (-0.07 %) but discontinuous, in agreement with the first-order nature of the transition observed by magnetic, transport, and heat capacity measurements. The anisotropic changes of the lattice parameters are Δa/a = 0.28 % and Δc/c = -0.64 % on cooling. A substantial change in the 4f and conduction electron hybridization, giving rise to an increased integrated DOS, occurs when the high temperature phase transforms to the low temperature phase. A moderate magnetocaloric effect at TN (ΔSM = -6.5 J/kg K and ΔTad = 4.5 K for ΔH= 50 kOe) has been measured using both magnetization and heat capacity data

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 OSTI2013

A second-order semi-implicit δf method for hybrid simulation

A second-order accurate semi-implicit Lorentz force ions, fluid electrons δf hybrid model has been developed using a “current closure” scheme. The model assumes quasi-neutrality and is fully electromagnetic. The implicit field solver improves numerical accuracy by separating the equilibrium terms in the presence of small perturbations. The equilibrium part of the generalized Ohm’s law is solved by direct matrix inversion along the direction of gradients for every Fourier mode in the other two directions, while the nonlinear part is solved iteratively. The simulation has been benchmarked on Alfvén waves, ion sound waves and whistler waves against analytical dispersion relation in a slab. In particular, the first-order and second-order schemes are compared by studying the numerical damping of whistler waves. The full evolution of the resistive tearing mode using the Harris sheet equilibrium is also investigated. The linear growth rate and mode structure are compared with the resistive MHD theory. Important tearing mode nonlinear phenomena such as the Rutherford regime and saturation are demonstrated. We also presented systematic study of Rutherford growth rates and saturation island width, which is consistent with previous MHD studies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

A high-order discontinuous Galerkin method for fluid–structure interaction with efficient implicit–explicit time stepping

We present a high-order accurate scheme for coupled fluid–structure interaction problems. The fluid is discretized using a discontinuous Galerkin method on unstructured tetrahedral meshes, and the structure uses a high-order volumetric continuous Galerkin finite element method. Standard radial basis functions are used for the mesh deformation. The time integration is performed using a partitioned approach based on implicit–explicit Runge–Kutta methods. The resulting scheme fully decouples the implicit solution procedures for the fluid and the solid parts, which we perform using two separate efficient parallel solvers. We demonstrate up to fifth order accuracy in time on a non-trivial test problem, on which we also show that additional subiterations are not required. We solve a benchmark problem of a cantilever beam in a shedding flow, and show good agreement with other results in the literature. Finally, we solve for the flow around a thin membrane at a high angle of attack in both 2D and 3D, and compare with the results obtained with a rigid plate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗

Refine your search to explore more results.