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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

3D Higher Order Modeling in the BEM/FEM Hybrid Formulation

Higher order divergence- and curl-conforming bases have been shown to provide significant benefits, in both convergence rate and accuracy, in the 2D hybrid finite element/boundary element formulation (P. Fink and D. Wilton, National Radio Science Meeting, Boulder, CO, Jan. 2000). A critical issue in achieving the potential for accuracy of the approach is the accurate evaluation of all matrix elements. These involve products of high order polynomials and, in some instances, singular Green's functions. In the 2D formulation, the use of a generalized Gaussian quadrature method was found to greatly facilitate the computation and to improve the accuracy of the boundary integral equation self-terms. In this paper, a 3D, hybrid electric field formulation employing higher order bases and higher order elements is presented. The improvements in convergence rate and accuracy, compared to those resulting from lower order modeling, are established. Techniques developed to facilitate the computation of the boundary integral self-terms are also shown to improve the accuracy of these terms. Finally, simple preconditioning techniques are used in conjunction with iterative solution procedures to solve the resulting linear system efficiently. In order to handle the boundary integral singularities in the 3D formulation, the parent element- either a triangle or rectangle-is subdivided into a set of sub-triangles with a common vertex at the singularity. The contribution to the integral from each of the sub-triangles is computed using the Duffy transformation to remove the singularity. This method is shown to greatly facilitate t'pe self-term computation when the bases are of higher order. In addition, the sub-triangles can be further divided to achieve near arbitrary accuracy in the self-term computation. An efficient method for subdividing the parent element is presented. The accuracy obtained using higher order bases is compared to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained using higher order bases are compared to those obtained with lower order bases for selected sample

Fink, P. W.↗
DOE OSTI2020

Using resonant energy X-ray diffraction to extract chemical order parameters in ternary semiconductors

II–IV–V 2 materials, ternary analogs to III–V materials, are emerging for their potential applications in devices such as LEDs and solar cells. Controlling cation ordering in II–IV–V 2 materials offers the potential to tune properties at nearly fixed compositions and lattice parameters. While tuning properties at a fixed lattice constant through ordering has the potential to be a powerful tool used in device fabrication, cation ordering also creates challenges with characterization and quantification of ordering. In this work, we investigate two different methods to quantify cation ordering in ZnGeP 2 thin films: a stretching parameter calculated from lattice constants , and an order parameter determined from the cation site occupancies ( S ). We use high resolution X-ray diffraction (HRXRD) to determine and resonant energy X-ray diffraction (REXD) to extract S . REXD is critical to distinguish between elements with similar Z -number ( e.g. Zn and Ge). We found that samples with a corresponding to the ordered chalcopyrite structure had only partially ordered S values. The optical absorption onset for these films occurred at lower energy than expected for fully ordered ZnGeP 2 , indicating that S is a more accurate descriptor of cation order than the stretching parameter. Since disorder is complex and can occur on many length scales, metrics for quantifying disorder should be chosen that most accurately reflect the physical properties of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantifying multipoint ordering in alloys

A central problem in multicomponent lattice systems is to systematically quantify multipoint ordering. Ordering in such systems is often described in terms of pairs, even though this is not sufficient when three-point and higher-order interactions are included in the Hamiltonian. Current models and parameters for multipoint ordering are often only applicable for very specific cases or require approximating a subset of correlated occupational variables on a lattice as being uncorrelated. In this paper, cluster order parameters are introduced to systematically quantify arbitrary multipoint ordering motifs in substitutional systems through direct calculations of normalized cluster probabilities. These parameters can describe multipoint chemical ordering in crystal systems with multiple sublattices, multiple components, and systems with reduced symmetry. These are defined in this paper and applied to quantify four-point chemical ordering motifs in platinum/palladium alloy nanoparticles that are of practical interest to the synthesis of catalytic nanocages. Impacts of chemical ordering on nanocage stability are discussed. It is demonstrated that approximating four-point probabilities from superpositions of lower-order pair probabilities is not sufficient in cases where three- and four-body terms are included in the energy expression. Furthermore, conclusions about the formation mechanisms of nanocages may change significantly when using common pair approximations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Designing phase sensitive probes of monopole superconducting order

Distinct from familiar s-, p-, or d-wave pairings, the monopole superconducting order represents a novel class of pairing order arising from nontrivial monopole charge of the Cooper pair. In the weak-coupling regime, this order can emerge when pairing occurs between Fermi surfaces with different Chern numbers in, for example, doped Weyl semimetal systems. However, the phase of monopole pairing order is not well-defined over an entire Fermi surface, making it challenging to design experiments sensitive to both its symmetry and topology. To address this, we propose a scheme based on symmetry and topological principles to identify this elusive pairing order through a set of phase-sensitive Josephson experiments. By examining the discrepancy between global and local angular momentum of the pairing order, we can unveil the monopole charge of the pairing order, including for models with higher pair monopole charge |q p |=1,2, and 3. We demonstrate the proposed probe of monopole pairing order through analytic and numerical studies of Josephson coupling in models of monopole superconductor junctions. This work opens a promising avenue to uncover the unique topological properties of monopole pairing orders and to distinguish them from known pairing orders based on spherical harmonic symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Neural network-based order parameter for phase transitions and its applications in high-entropy alloys

Phase transition is one of the most important phenomena in nature and plays a central role in materials design. All phase transitions are characterized by suitable order parameters, including the order–disorder phase transition. However, finding a representative order parameter for complex systems is non-trivial, such as for high-entropy alloys. Given the strength of dimensionality reduction of a variational autoencoder (VAE), we introduce a VAE-based order parameter. Here, we propose that the Manhattan distance in the VAE latent space can serve as a generic order parameter for order–disorder phase transitions. The physical properties of our order parameter are quantitatively interpreted and demonstrated by multiple refractory high-entropy alloys. Using this order parameter, a generally applicable alloy design concept is proposed by mimicking the natural mixing process of elements. Our physically interpretable VAE-based order parameter provides a computational technique for understanding chemical ordering in alloys, which can facilitate the development of rational alloy design strategies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Theory of topological defects and textures in two-dimensional quantum orders with spontaneous symmetry breaking

In this article, we consider two-dimensional (2d) quantum many-body systems with long-range orders, where the only gapless excitations in the spectrum are Goldstone modes of spontaneously broken continuous symmetries. To understand the interplay between classical long-range order of local order parameters and quantum order of long-range entanglement in the ground states, we study the topological point defects and textures of order parameters in such systems. We show that the universal properties of point defects and textures are determined by the remnant symmetry enriched topological order in the symmetry-breaking ground states with a nonfluctuating order parameter, and provide a classification for their properties based on the inflation-restriction exact sequence. We highlight a few phenomena revealed by our theory framework. First, in the absence of intrinsic topological orders, we show a connection between the symmetry properties of point defects and textures to deconfined quantum criticality. Second, when the symmetry-breaking ground state has intrinsic topological orders, we show that the point defects can permute different anyons when braided around. They can also obey projective fusion rules in the sense that multiple vortices can fuse into an Abelian anyon, a phenomenon for which we coin “defect fractionalization.” Finally, we provide a formula to compute the fractional statistics and fractional quantum numbers carried by textures (skyrmions) in Abelian topological orders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1977

The 4th order GISS model of the global atmosphere

The new GISS 4th order model of the global atmosphere is described. It is based on 4th order quadratically conservative differences with the periodic application of a 16th order filter on the sea level pressure and potential temperature equations, a combination which is approximately enstrophy conserving. Several short range forecasts indicate a significant improvement over 2nd order forecasts with the same resolution (approximately 400 km). However the 4th order forecasts are somewhat inferior to 2nd order forecasts with double resolution. This is probably due to the presence of short waves in the range between 1000 km and 2000 km, which are computed more accurately by the 2nd order high resolution model. An operation count of the schemes indicates that with similar code optimization, the 4th order model will require approximately the same amount of computer time as the 2nd order model with the same resolution. It is estimated that the 4th order model with a grid size of 200 km provides enough accuracy to make horizontal truncation errors negligible over a period of a week for all synoptic scales (waves longer than 1000 km).

Kalnay-Rivas, E.↗
NASA NTRS2002

Control of Initialized Fractional-Order Systems

Due to the importance of historical effects in fractional-order systems, this paper presents a general fractional-order control theory that includes the time-varying initialization response. Previous studies have not properly accounted for these historical effects. The initialization response, along with the forced response, for fractional-order systems is determined. Stability properties of fractional-order systems are presented in the complex Airplane, which is a transformation of the s-plane. Time responses are discussed with respect to pole positions in the complex Airplane and frequency response behavior is included. A fractional-order vector space representation, which is a generalization of the state space concept, is presented including the initialization response. Control methods for vector representations of initialized fractional-order systems are shown. Nyquist, root-locus, and other input-output control methods are adapted to the control of fractional-order systems. Finally, the fractional-order differintegral is generalized to continuous order-distributions that have the possibility of including a continuum of fractional orders in a system element.

Hartly, Tom T.↗
NASA NTRS2002

Control of Initialized Fractional-Order Systems

Due to the importance of historical effects in fractional-order systems, this paper presents a general fractional-order control theory that includes the time-varying initialization response. Previous studies have not properly accounted for these historical effects. The initialization response, along with the forced response, for fractional-order systems is determined. Stability properties of fractional-order systems are presented in the complex w-plane, which is a transformation of the s-plane. Time responses are discussed with respect to pole positions in the complex w-plane and frequency response behavior is included. A fractional-order vector space representation, which is a generalization of the state space concept, is presented including the initialization response. Control methods for vector representations of initialized fractional-order systems are shown. Nyquist, root-locus, and other input-output control methods are adapted to the control of fractional-order systems. Finally, the fractional-order differintegral is generalized to continuous order-distributions that have the possibility of including a continuum of fractional orders in a system element.

Hartley, Tom T.↗
NASA NTRS2002

Higher Order, Hybrid BEM/FEM Methods Applied to Antenna Modeling

In this presentation, the authors address topics relevant to higher order modeling using hybrid BEM/FEM formulations. The first of these is the limitation on convergence rates imposed by geometric modeling errors in the analysis of scattering by a dielectric sphere. The second topic is the application of an Incomplete LU Threshold (ILUT) preconditioner to solve the linear system resulting from the BEM/FEM formulation. The final tOpic is the application of the higher order BEM/FEM formulation to antenna modeling problems. The authors have previously presented work on the benefits of higher order modeling. To achieve these benefits, special attention is required in the integration of singular and near-singular terms arising in the surface integral equation. Several methods for handling these terms have been presented. It is also well known that achieving ~he high rates of convergence afforded by higher order bases may als'o require the employment of higher order geometry models. A number of publications have described the use of quadratic elements to model curved surfaces. The authors have shown in an EFIE formulation, applied to scattering by a PEC .sphere, that quadratic order elements may be insufficient to prevent the domination of modeling errors. In fact, on a PEC sphere with radius r = 0.58 Lambda(sub 0), a quartic order geometry representation was required to obtain a convergence benefi.t from quadratic bases when compared to the convergence rate achieved with linear bases. Initial trials indicate that, for a dielectric sphere of the same radius, - requirements on the geometry model are not as severe as for the PEC sphere. The authors will present convergence results for higher order bases as a function of the geometry model order in the hybrid BEM/FEM formulation applied to dielectric spheres. It is well known that the system matrix resulting from the hybrid BEM/FEM formulation is ill -conditioned. For many real applications, a good preconditioner is required to obtain usable convergence from an iterative solver. The authors have examined the use of an Incomplete LU Threshold (ILUT) preconditioner . to solver linear systems stemming from higher order BEM/FEM formulations in 2D scattering problems. Although the resulting preconditioner provided aD excellent approximation to the system inverse, its size in terms of non-zero entries represented only a modest improvement when compared with the fill-in associated with a sparse direct solver. Furthermore, the fill-in of the preconditioner could not be substantially reduced without the occurrence of instabilities. In addition to the results for these 2D problems, the authors will present iterative solution data from the application of the ILUT preconditioner to 3D problems.

Fink, P. W.↗
DOE OSTI2025

Explicit Runge–Kutta Methods that Alleviate Order Reduction

Explicit Runge–Kutta (RK) methods are susceptible to a reduction in the observed order of convergence when applied to an initial boundary value problem with time-dependent boundary conditions. We study conditions on explicit RK methods that guarantee high order convergence for linear problems; we refer to these conditions as weak stage order conditions. We prove a general relationship between the method’s order, weak stage order, and number of stages. Furthermore, we derive explicit RK methods with high weak stage order and demonstrate, through numerical tests, that they avoid the order reduction phenomenon up to any order for linear problems and up to order three for nonlinear problems.

explicit Runge–Kutta↗
DOE OSTI2021

Spectral deferred correction methods for high-order accuracy in poroelastic problems

In this work, we investigate high-order accuracy in time integration by examining two operator splitting methods for poroelastic problems: the two-pass and the spectral deferred correction (SDC) methods. To enhance the order of accuracy, the two-pass method partitions a coupled operator symmetrically, whereas the SDC method corrects truncation errors by establishing an error equation. These high-order methods are applied to underlying solution strategies, i.e., monolithic, fixed-stress sequential, and undrained sequential methods. We observe that semi-discretized systems from spatial discretization have forms similar to those of index-1 differential algebraic equations (DAEs), causing order reduction against the two-pass method when it is used in conjunction with either the monolithic or sequential method. On the other hand, the SDC in conjunction with the monolithic method exhibits the desired second-order accuracy in poroelastic problems while increasing the order of accuracy for index-1 DAEs. However, the SDC in conjunction with either of the two sequential methods does not achieve the desired order of accuracy, and maintains first order because the flow equation for poroelasticity has an additional approximation associated with the volumetric strain rate term, which does not yield exactly the same forms as those of conventional DAEs. Thus, the monolithic SDC method can achieve higher-order accuracy, but may require higher computational costs because it involves solving matrix systems larger than those for the sequential methods.

02 PETROLEUM↗
DOE OSTI2020

Matrix-free preconditioning for high-order H (curl) discretizations

The greater arithmetic intensity of high-order finite element discretizations makes them attractive for implementation on next-generation hardware, but assembly of high-order finite element operators as matrices is prohibitively expensive. As a result, the development of general algebraic solvers for such operators has been an open research challenge. Fast matrix-free application of high-order operators has received significant attention in the literature in the context of Poisson-type problems, but preconditioners and solvers for inverting more general operators are not very well-developed. In this paper, we consider the problem of preconditioning a definite Maxwell operator at high polynomial order without assembling a matrix. We show that given efficient preconditioners for high-order H 1 finite element problems on the same mesh, efficient H(curl) preconditioners can be constructed in an auxiliary space framework. We demonstrate the resulting preconditioners in a practical setting with tensor-product basis functions on an unstructured mesh of quadrilaterals. Overall, our approach uses a sparsified H 1 solver constructed on a low-order mesh of the nodal points of the underlying high-order space, and we show that the resulting H(curl) preconditioner is effective at very high polynomial orders for two-dimensional model problems with complicated geometry, varying piecewise constant coefficients, and curved elements. The resulting preconditioner scales with nearly optimal O(p d+1 ) floating point operation count and optimal O(p d ) memory transfer requirements, outperforming existing Maxwell preconditioners in the high-order regime.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A second-order distributed memory parallel fast sweeping method for the Eikonal equation

The Eikonal equation is used to calculate wave propagation and distance fields, and due to its complexity requires numerical treatment for its solution. In this work, we present a second-order distributed memory parallel fast sweeping method. The second-order solution switches on a two-point stencil when two upwind points are available, and reverts to first-order otherwise. In all examples, the second-order method improves the solution over the first-order, allowing for significant savings in memory while achieving the same accuracy. Parallelization over distributed memory saw good weak scaling with optimal convergence. The computational time for second-order was approximately 2.5 times slower than first-order, where the largest amount of mesh points ran on 144 cores (512 GB) was ≈20 billion. The savings in memory from the second-order method combined with the distributed memory algorithm result in the ability to solve problems much larger than are possible with the serial first-order method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Absence of precursor incommensurate charge order in electronic nematic Ba 0.35 Sr 0.65 Ni 2 As 2

Recent discoveries of charge order and electronic nematic order in the iron-based superconductors and cuprates have pointed towards the possibility of nematic and charge fluctuations playing a role in the enhancement of superconductivity. The Ba 1–x Sr x Ni 2 As 2 system, closely related in structure to the BaFe 2 As 2 system, has recently been shown to exhibit both types of ordering without the presence of any magnetic order. We report single-crystal x-ray diffraction experiments on Ba 0.35 Sr 0.65 Ni 2 As 2 , providing evidence that the previously reported incommensurate charge order with wave-vector (0, 0.28, 0) tet in the tetragonal state of BaNi 2 As 2 vanishes by 65% Sr substitution together with nematic order. Furthermore, our measurements suggest that the nematic and incommensurate charge orders are closely tied in the tetragonal state and show that the (0, 0.33, 0) tri charge ordering in the triclinic phase of BaNi 2 As 2 evolves to become (0, 0.5, 0) tri charge ordering at x = 0.65 before vanishing at x = 0.71.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hierarchy among the crystal lattice, charge density wave, and superconducting orders in transition metal dichalcogenides

Using high-energy x-ray scattering and large-scale three-dimensional (3D) structure modeling, in this study we investigate the relationship between the crystal lattice, charge density wave (CDW), and superconducting (SC) orders in transition metal dichalcogenides (TMDs). In particular, we systematically substitute Te for Se in Ta-Se-Te solid solutions, determine changes in their crystal lattice, and relate them to changes in the CDW transition temperature, $T_{CDW}$, and SC critical temperature, $T_c$ . We find that strong lattice distortions such as buckling of Ta layers are detrimental to the CDW and SC orders. The presence of a perfect lattice order in two dimensions is a prerequisite to the emergence of CDWs but insufficient to achieve a SC ordered state. For the SC order to emerge, the Ta sublattice should also appear periodic in 3D. Local chemical disorder may promote the SC order, and the perfectness of Ta coordination polyhedra is a factor contributing to its strength. A hierarchical relationship among the crystal lattice, CDW, and SC orders thus appears to exist in TMDs in a sense that different degrees of crystal lattice order-disorder promote and maintain the CDW and SC orders to a different extent, offering an opportunity to control the latter through modifying the former by rational design. Our findings are a step towards a better understanding of the correlation between lattice and electronic degrees of freedom in TMDs. We also demonstrate an efficient experimental approach to study them in fine detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Landau-Ginzburg theory of charge density wave formation accompanying lattice and electronic long-range ordering

Here we propose an analytical Landau-Ginzburg (LG) theory of the charge density waves coupled with lattice and electronic long-range order parameters. Examples of long-range order include the electronic wave function of superconducting Cooper pairs, structural distortions, electric polarization, and magnetization. We formulate the LG free energy density as a power expansion with respect to the charge density and other long-range order parameters as well as their spatial gradients and biquadratic coupling terms. We introduced a biquadratic coupling between the charge density gradient and long-range order parameters as well as nonlinear higher gradients of the long-range order parameters. The biquadratic gradient coupling is critical to the appearance of different spatially modulated phases in charge-ordered ferroics and high-temperature superconductors. We derived the thermodynamic conditions for the stability of the spatially modulated phases, which are the intertwined spatial waves of charge density and lattice/electronic long-range order. The analytical expressions for the energies of different phases, corresponding order parameters, charge density waves amplitudes, and modulation periods obtained in this paper can be employed to guide the comprehensive physical explanation, deconvolution, and Bayesian analysis of experimental data on quantum materials ranging from charge-ordered ferroics to high-temperature superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Coarsening dynamics of Ising-nematic order in a frustrated Heisenberg antiferromagnet

We study the phase ordering dynamics of the classical antiferromagnetic 𝐽 1 −𝐽 2 (nearest-neighbor and next-nearest-neighbor couplings) Heisenberg model on the square lattice in the strong frustration regime (𝐽 2 /𝐽 1 > 1/2). While thermal fluctuations preclude any long-range magnetic order at finite temperatures, the system exhibits a long-range spin-driven nematic phase at low temperatures. The transition into the nematic phase is further shown to belong to the two-dimensional Ising universality class based on the critical exponents near the phase transition. Our large-scale stochastic Landau-Lifshitz-Gilbert simulations find a two-stage phase ordering when the system is quenched from a high-temperature paramagnetic state into the nematic phase. In the early stage, collinear alignments of spins lead to a locally saturated Ising-nematic order. Once domains of well-defined Ising order are developed, the late-stage relaxation is dominated by curvature-driven domain coarsening, as described by the Allen-Cahn equation. The characteristic size of Ising-nematic domains scales as the square root of time, similar to the kinetic Ising model described by the time-dependent Ginzburg-Landau theory. Our results confirm that the late-stage ordering kinetics of the spin-driven nematic, which is a vestigial order of the frustrated Heisenberg model, belongs to the dynamical universality class of a nonconserved Ising order. Interestingly, the system shows no violation of the superuniversality hypothesis under weak bond disorder. The dynamic scaling invariance is preserved in the presence of weak bond disorder. Here, we also discuss possible applications of our results to materials for which vestigial Ising-nematic order is realized.

Antiferromagnets↗