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

Observation of a C -type short-range antiferromagnetic order in layer spacing expanded FeS

Here, we report neutron diffraction studies of FeS single crystals obtained from Rb x Fe 2-y S 2 single crystals via a hydrothermal method. While no $\sqrt{5}$×$\sqrt{5}$ iron vacancy order or block antiferromagnetic order typical of Rb x Fe 2-y S 2 is found in our samples, we observe C-type short-range antiferromagnetic order with moments pointed along the c axis hosted by a different phase of FeS with an expanded interlayer spacing. The Néel temperature for this magnetic order is determined to be 170 ± 4 K. Our finding of a variant FeS structure hosting this C-type antiferromagnetic order demonstrates that the known FeS phase synthesized in this method is in the vicinity of a magnetically ordered ground state, providing insights into understanding a variety of phenomena observed in FeS and the related FeSe 1-x S x iron chalcogenide system.

36 MATERIALS SCIENCE↗
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↗
NASA NTRS1975

Higher order accurate difference solutions of fluid mechanics problems by a compact differencing technique

A general fourth order differencing scheme is developed and applied to three viscous test problems to verify the accuracy and applicability of the technique. The procedure is atypical since only three nodes are necessary to obtain the desired fourth order accuracy. This is accomplished by a differencing technique which considers the function and all necessary derivatives as unknowns. The relations for these derivatives yield simple tridiagonal equations which can be easily solved. Comparisons of the fourth order results with those computed using second order methods are presented for the test problems and clearly indicate that the accuracy achieved by these fourth order computations is always significantly better than current second order procedures.

Hirsh, R. S.↗
NASA NTRS1978

The stability of numerical methods for second order ordinary differential equations

An important characterization of a numerical method for first order ODE's is the region of absolute stability. If all eigenvalues of the linear problem dy/dt = Ay are inside this region, the numerical method is stable. If the second order system d/dt(dy/dt) = 2Ady/dt - By is solved as a first order system, the same result applies to the eigenvalues of the generalized eigenvalue problem (lambda-squared)I 2(lambda)A + B. No such region exists for general methods for second order equations, but in some cases a region of absolute stability can be defined for methods for the single second order equation d/dt(dy/dt) = 2ady/dt - by. The absence of a region of absolute stability can occur when different members of a system of first order equations are solved by different methods.

Gear, C. W.↗
NASA NTRS1984

Second-order-accurate spatial differencing for the transonic small-disturbance equation

Current methods for calculating transonic flows with the small-disturbance potential equation are typically only first-order accurate in the supersonic regions of the flow. However, calculations using the full-potential equation show significant improvements in accuracy when second-order methods are used instead of first-order methods. In this paper, algorithms for the small-disturbance equations are presented, for both steady and unsteady flows, with spatial differencing that is second-order accurate in both the subsonic and supersonic regions of the flow. These algorithms are stable, simple extensions of implicit monotone algorithms; the former are only first-order accurate spatially in the supersonic regions. Several calculations in one and two dimensions have been made, and the first- and second-order calculations are compared. In one dimension, the calculations include the four types of shock-wave motion. In two dimensions, the comparisons include steady flow over a Korn airfoil and unsteady flow over a pitching airfoil. All the comparisons in two dimensions show that the new methods are just as stable numerically as the old methods, but improve the accuracy of the solutions. The algorithm improvements can be implemented in present computer codes by making minor coding modifications.

Goorjian, P. M.↗
NASA NTRS1984

Substrate-induced orientational order in the isotropic phase of liquid crystals

Nematic order induced near a solid boundary in an otherwise isotropic liquid crystal is studied theoretically, at temperatures just above the bulk nematic-isotropic phase transition. Three distinct regimes are found, depending on the strength of orientational torques at the boundary: (1) strong orientational order, (2) strong orientational order followed by a first-order transition to a state of weak orientational order as temperature is raised, and (3) a state of weak orientational order.

Mauger, A.↗
NASA NTRS1987

Fixed order compensator design based on frequency shaped cost functionals

This paper treats the design of fixed order compensators using frequency shaped cost functionals. The formulation is done in an output feedback setting which exploits an observer canonical form to represent the compensator dynamics. The major advantages of the design approach are that the order of the compensator is fixed by the design process, and it is not necessary to realize the frequency shaping dynamics as an integral part of the compensator design. This avoids the two-step approach of full state feedback design followed by the design of an observer, which in the case of a full order observer requires a compensator of order n + nf where n is the order of the plant and nf is the order of the frequency shaping realization. The formulation also precludes the use of direct feedback of the plant output, thus reducing the effect of sensor noise and improving the robustness to high frequency unmodeled dynamics. A structural vibration model describing the fast dynamics of a lightweight flexible arm is used as an example to illustrate the design procedure.

Calise, Anthony J.↗
NASA NTRS1988

Finite-dimensional approximation for optimal fixed-order compensation of distributed parameter systems

In controlling distributed parameter systems it is often desirable to obtain low-order, finite-dimensional controllers in order to minimize real-time computational requirements. Standard approaches to this problem employ model/controller reduction techniques in conjunction with LQG theory. In this paper we consider the finite-dimensional approximation of the infinite-dimensional Bernstein/Hyland optimal projection theory. This approach yields fixed-finite-order controllers which are optimal with respect to high-order, approximating, finite-dimensional plant models. The technique is illustrated by computing a sequence of first-order controllers for one-dimensional, single-input/single-output, parabolic (heat/diffusion) and hereditary systems using spline-based, Ritz-Galerkin, finite element approximation. Numerical studies indicate convergence of the feedback gains with less than 2 percent performance degradation over full-order LQG controllers for the parabolic system and 10 percent degradation for the hereditary system.

Bernstein, Dennis S.↗
NASA NTRS1988

Third-order multi-dimensional Euler/Navier-Stokes solver

A nonoscillatory multidimensional advection code is presented which is applicable to the study of very high Reynolds number steady flow, and results are obtained for the oblique advection of a step profile. The scheme is based on an exponential upwinding or linear extrapolation refinement of a previous QUICK scheme, allowing sharp monotonic resolution of the step while retaining third-order accuracy. The resulting simple flux-limiting criteria have been applied to a higher order method and a formally third-order supercompressive algorithm. Results are presented for first-order, second-order, and third-order upwinding.

Leonard, B. P.↗
NASA NTRS1990

Robust eigensystem assignment for state estimators using second-order models

A novel design of a state estimator is presented using second-order dynamic equations of mechanical systems. The eigenvalues and eigenvectors of the state estimator are assigned by solving the second-order eigenvalue problem of the structural system. Three design method for the state estimator are given in this paper. The first design method uses collocated sensors to measure the desired signals and their derivatives. The second design method uses prefilters to shift signal phases to obtain estimates of the signal derivatives. These two methods are used to build a second-order state estimator model. The third design method is the conventional one which converts a typical second-order dynamic model to a first-order model, and then builds a state estimator based on the first-order model. It is shown that all the three designs for state estimation are similar. A numerical example representing a large space structure is given for illustration of the design methods presented in this paper.

Juang, Jer-Nan↗