Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ordering transformation”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Adaptive clipping‐and‐redistribution algorithms for bounded and conservative high‐order interpolations applied to discontinuous and reactive flows

Abstract A new adaptive clipping‐and‐redistribution method is presented which provides bounds‐preservation for multidimensional interpolation in the context of high‐order finite‐volume discretizations with adaptive mesh refinement (AMR). The underlying finite‐volume method (FVM) for the computational fluid dynamics applications is fourth‐order accurate for smooth solutions and utilizes AMR for computational efficiency in solving multiscale problems involving turbulence and combustion. High‐order interpolation between different AMR levels is required. However, this operation often leads to numerical issues because combustion species must have physical bounds preserved. The present study overcomes two major challenges in the development of the high‐order interpolation method. First, the method needs to be bound‐preserving near extrema or discontinuities to prevent the emergence of unphysical oscillations while maintaining fourth‐order accuracy in smooth flows. Second, the method needs to satisfy the conservation requirement in multiple dimensions, particularly in the context of curvilinear coordinate transformations. Additionally, the method is designed to be localized and computationally inexpensive. The new interpolation scheme is demonstrated by solving reacting flows, which are extremely sensitive to unphysical overshoots in conserved quantities. The test problems are shock‐induced ‐ combustion and a ‐air flame in a practical bluff‐body combustor. Results show the method prevents new extrema near discontinuities while maintaining high‐order accuracy in smooth regions. In particular, the method is extremely beneficial for combustion with stiff chemistry. With the proposed new method, even if flame fronts cross AMR interfaces or new grids are created in the vicinity of the flame, solution stability is retained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Local structural ordering determines the mechanical damage tolerance of amorphous grain boundary complexions

Amorphous grain boundary complexions act as toughening features within a microstructure because they can absorb dislocations more efficiently than traditional grain boundaries. This toughening effect should be a strong function of the local internal structure of the complexion, which has recently been shown to be determined by grain boundary crystallography. To test this hypothesis, molecular dynamics are used here to simulate dislocation absorption and damage nucleation for complexions with different distributions of structural short-range order. Here, the complexion with a more disordered structure away from the dislocation absorption site is actually found to better resist crack nucleation, as damage tolerance requires delocalized deformation and the operation of shear-transformation zones through the complexion thickness. The more damage tolerant complexion accommodates plastic strain efficiently within the entire complexion, providing the key mechanistic insight that local patterning and asymmetry of structural short-range order controls the toughening effect of amorphous complexions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Efficient Fourier transforms for transverse momentum dependent distributions

Hadron production at low transverse momenta in semi-inclusive deep inelastic scattering can be described by transverse momentum dependent (TMD) factorization. This formalism has also been widely used to study the Drell-Yan process and back-to-back hadron pair production in $e^+e^-$ collisions. These processes are the main ones for extractions of TMD parton distribution functions and TMD fragmentation functions, which encode important information about nucleon structure and hadronization. One of the most widely used TMD factorization formalism in phenomenology formulates TMD observables in coordinate $b_\perp$-space, the conjugate space of the transverse momentum. The Fourier transform from $b_\perp$-space back into transverse momentum space is sufficiently complicated due to oscillatory integrands that it requires a careful and computationally intensive numerical treatment in order to avoid potentially large numerical errors. Within the TMD formalism, the azimuthal angular dependence is analytically integrated and the two-dimensional $b_\perp$ integration reduces to a one-dimensional integration over the magnitude $b_\perp$. In this paper we develop a fast numerical Hankel transform algorithm for such a $b_\perp$-integration that improves the numerical accuracy of TMD calculations in all standard processes. Libraries for this algorithm are implemented in Python 2.7 and 3, C++, as well as FORTRAN77. All packages are made available open source.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Quantifying the Dynamics of Protein Self-Organization Using Deep Learning Analysis of Atomic Force Microscopy Data

The dynamics of protein self-assembly on the inorganic surface and the resultant geometric patterns are visualized using high-speed atomic force microscopy. The time dynamics of the classical macroscopic descriptors such as 2D fast Fourier transforms, correlation, and pair distribution functions are explored using the unsupervised linear unmixing, demonstrating the presence of static ordered and dynamic disordered phases and establishing their time dynamics. Here, the deep learning (DL)-based workflow is developed to analyze detailed particle dynamics and explore the evolution of local geometries. Finally, we use a combination of DL feature extraction and mixture modeling to define particle neighborhoods free of physics constraints, allowing for a separation of possible classes of particle behavior and identification of the associated transitions. Overall, this work establishes the workflow for the analysis of the self-organization processes in complex systems from observational data and provides insight into the fundamental mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Investigation of the background in coherent J/ψ production at the EIC

Understanding various fundamental properties of nucleons and nuclei is among the most important scientific goals at the upcoming Electron-Ion Collider (EIC). With the unprecedented opportunity provided by the next-generation machine, the EIC might provide definitive answers to many standing puzzles and open questions in modern nuclear physics. We investigate one of the golden measurements proposed at the EIC, which is to obtain the spatial gluon density distribution within a lead (Pb) nucleus. The proposed experimental process is the exclusive J/ψ vector-meson production off the Pb nucleus: e + Pb → e' + J/ψ + Pb'. The Fourier transformation of the momentum transfer |t| distribution of the coherent diffraction is the transverse gluon spatial distribution. In order to measure it, the experiment has to overcome an overwhelmingly large background arising from the incoherent diffractive production, where the nucleus Pb' mostly breaks up into fragments of particles in the far-forward direction close to the hadron-going beam rapidity. We systematically study the rejection of incoherent J/ψ production by vetoing products from these nuclear breakups—protons, neutrons, and photons—which is based on the BeAGLE event generator and the most up-to-date EIC Far-forward Interaction Region design. The achieved vetoing efficiency, the ratio between the numbers of vetoed events and total incoherent events, ranges from about 80% to 99% depending on |t|. Assuming a 5% smearing applied to the reconstructed |t| resolution in the Sar t re model, this vetoing efficiency can suppress the incoherent background to at least the first minimum of the coherent |t| distribution. Experimental and accelerator machine challenges as well as potential improvements are discussed herein.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A Semi-Analytical Approach for State-Space Electromagnetic Transient Simulation

Here, this paper proposes a semi-analytical approach for efficient and accurate electromagnetic transient (EMT) simulation of a power grid. The approach first derives a high-order semi-analytical solution (SAS) of the grid’s state-space EMT model using the differential transformation (DT), and then evaluates the solution over enlarged, variable time steps to significantly accelerate the simulations while maintaining its high accuracy on detailed fast EMT dynamics. The approach also addresses switches during large time steps by using a limit violation detection algorithm with a binary search-enhanced quadratic interpolation. Case studies are conducted on EMT models of the IEEE 39-bus system and large-scale systems to demonstrate the merits of the new simulation approach against traditional numerical methods.

electromagnetic transient↗
DOE OSTI2021

Coupling morphological and magnetic anisotropy for assembling tetragonal colloidal crystals

Morphological and magnetic anisotropy can be combined in colloidal assembly to create unconventional secondary structures. We show here that magnetite nanorods interact along a critical angle, depending on their aspect ratios and assemble into body-centered tetragonal colloidal crystals. Under a magnetic field, size-dependent attractive and repulsive domains develop on the ends and center of the nanorods, respectively. Our joint experiment-computational multiscale study demonstrates the presence of a critical angle in the attractive domain, which defines the equilibrium bonding states of interacting rods and leads to the formation of non–close-packed yet hard-contact tetragonal crystals. Small-angle x-ray scattering measurement attributes the perfect tetragonal phase to the slow assembly kinetics. The crystals exhibit brilliant structural colors, which can be actively tuned by changing the magnetic field direction. These highly ordered frameworks and well-defined three-dimensional nanochannels may offer new opportunities for manipulating nanoscale chemical transformation, mass transportation, and wave propagation.

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

Interplay of broken symmetry and delocalized excitations in the insulating state of 1 T –Ta⁢S 2

Coexistence of localized and extended excitations is central to the macroscopic properties of correlated materials. For 5⁢d transition-metal compounds, electron correlations alone generally do not lead to a metal-insulator (Mott) transition, with insulating behavior usually resulting from their coupling with magnetic ordering and/or structural distortions. 1⁢T–Ta⁢S 2 is a prototypical example of such correlated insulating behavior, with a high-symmetry metallic phase transforming into a distorted, charge-density wave (CDW) insulating state at low temperatures. The nature of the insulating phase as well as the existence and relevance of the localized electron physics remains debated. We resolve this standing controversy in 1⁢T–Ta⁢S 2 combining resonant inelastic x-ray spectroscopy and first-principles calculations. We observe five electronic excitations arising from the interband transitions of the Ta 5d orbitals and the S 3⁢p ligand state, with none of the excitations on the order of the Mott gap. These excitations cannot be explained within the framework of standard multiplet calculations that assume a localized wave function, but instead, are captured by a band-theory framework accounting for the low symmetry of the crystal field in the CDW state. Our findings suggest that the electronic properties of 1T–Ta⁢S 2 are dominated in the visible by both plasmonic quasiparticles and interband transitions, with no resonance associated with a putative Mott transition observed in the 0–3 eV energy range. Finally, our discovery provides insights into the electron localization and the Mott vs band-insulator debate in 1T–Ta⁢S 2 and other transition-metal materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Coupling of organic cation and inorganic lattice in methylammonium lead halide perovskites: Insights into a pressure-induced isostructural phase transition

Organic-inorganic hybrid metal-halide perovskite (OIHP) materials provide a tunable platform for engineering their optoelectronic properties. Although several high-pressure studies have been conducted from the OIHP family of single crystals and films, the exact nature of the dynamic coupling of the CH 3 NH 3 (MA) cation with the octahedral lattice framework and the mechanisms responsible for the structural phase transformation under pressure are not well captured. By combined photoluminescence (PL), synchrotron-based x-ray diffraction, and Raman-scattering studies as a function of pressure from methylammonium lead bromide (MAPbBr 3 ), we shed light on an isostructural phase transition due to the coupling of the MA cation and the PbBr 6 lattice through hydrogen bonding. The sharp discontinuities at ~ 1 GPa and ~ 3 GPa in the PL peak positions correlate with the structural changes observed in high-pressure XRD and Raman-scattering studies. The electronic band edge as a function of pressure is calculated within density-functional theory. The PL peak position, intensity and width of the excitonic peak show significant changes at 2 GPa, which corroborate the changes observed in high-pressure Raman-scattering studies. The frequencies of the lattice modes and the C-H/N-H bending and stretching modes of the MA cation show anomalous changes and other nuances at 2 GPa. The suppression of rotational and orientational disorder of the organic moiety is initiated at 2 GPa and the ordering is completed by 3.0 GPa, leading to an order-disorder type cubic II (Im$\overline{3}$) to orthorhombic ( Pnma ) phase transition. Along with the revelation of an isostructural transformation at 2 GPa, this paper highlights the impact of molecular vibrations on the electronic properties of MAPbBr 3 under pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effects of dimensionality on the electronic structure of Ruddlesden-Popper chromates Sr n+1 ⁢Cr n ⁢O 3⁢n+1

Transition-metal oxides host a wide variety of electronic phenomena that can be significantly influenced by the effective dimensionality of the system under consideration. These include charge, spin, and orbital orderings, as well as unconventional superconductivity. In this context, the Ruddlesden-Popper chromates Sr n+1 Cr n O 3n+1 emerge as a particularly intriguing series of materials. Formally, the chromium atom displays a rather special 4+ oxidation state throughout the entire series. However, the effective dimensionality changes from quasi-2D to 3D as n increases from 1 to ∞. As a result, the insulating antiferromagnetic behavior observed for n = 1, 2, 3 transforms into itinerant antiferromagnetism with reduced transition temperature for the n = ∞ end member of the series, i.e., the perovskite SrCrO 3 . Further, distinct orbital orderings with exotic singlet states have been predicted for these systems. However, the lack of single-crystal bulk or thin-film samples has made experimental progress difficult. Here we demonstrate the synthesis of thin films of the perovskite SrCrO 3 and the associated layered chromates via oxide molecular beam epitaxy for n = 1 to n = 5. Our electrical transport measurements reveal a gradual evolution from a strongly insulating state in Sr 2 CrO 4 to a metallic state in the end member SrCrO 3 . X-ray absorption spectroscopy measurements demonstrate a varying hybridization strength of the Cr 4+ valence electrons across the series, helping to explain the trend in conduction. Finally, density functional theory calculations further confirm the observed transport trend and identify additional distortions present in the system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Deep learning model for fast, science-based forecasting of fluid migration along faults in geologic carbon storage scenarios

Effective long-term geologic storage depends on robust site selection and credible, science-based forecasting of subsurface behavior to ensure storage integrity. For this work, we develop a deep learning–based reduced-order model (ROM) to quantify potential carbon dioxide (CO₂) and brine migration through geological faults. The ROM combines a Transformer model for binary classification and a Stacked Ensemble for regression, trained on a comprehensive dataset generated from 1400 physics-based reservoir simulations. Key geologic and operational parameters—including fault geometry, reservoir structure, and injection conditions—were systematically varied to capture a wide range of fluid migration scenarios. The ROM accurately predicts the onset of migration, cumulative migration volumes of both CO₂ and brine, and associated migration rates, as compared to an independent set of validation simulations, while significantly reducing computational cost compared to traditional simulation methods. Model performance was evaluated across diverse fault configurations, revealing that shallow reservoir geometry and fault angle are among the most influential factors governing migration behavior. Sensitivity analysis using SHapley Additive exPlanations (SHAP) provided interpretability, revealing distinct patterns in how geological and operational features drive transient versus cumulative migration outcomes. The ROM’s ability to rapidly simulate fault migration scenarios enables efficient sensitivity analyses, scenario evaluations, and decision support for site selection and monitoring design. This approach enhances the safety, scalability, and long-term operational performance of geologic carbon storage (GCS) systems by providing a robust, interpretable tool for predicting subsurface fluid migration and assessing fault-related migration potential.

42 ENGINEERING↗
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress–strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress-strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

Levitas, Valery↗
DOE OSTI2021

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg (1997). This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple- ζ basis set predicts the spin splittings of [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ to be –35.7 and –17.1 kcal mol –1 , respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

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

Hydroxyl radical‐driven transformations of bisphenol A and 2,4‐dinitroanisole: Experimental and computational analysis

Abstract This study used experimental and computational analysis to investigate the advanced oxidation of bisphenol A (BPA) and 2,4‐dinitroanisole (DNAN). The pseudo first‐order reaction rate constants depended on the molar peroxide ratio and were between 0.13 and 0.28 min −1 for BPA and between 0.018 and 0.032 min −1 for DNAN. The kinetic differences appear to be due in part to the energy requirements for oxidation, which depended on the reaction mechanism but were typically lower for BPA than they were for DNAN. Density functional theory (DFT) was used to develop transformation pathways that included experimentally‐detected byproducts. The most energetically favored pathway for BPA oxidation begins with the formation of hydroxylated derivatives, while for DNAN, the most energetically favorable degradation pathway begins with the substitution of the methoxy group. Overall, these findings demonstrate the power of combining experimental and computational tools to reveal transformation mechanisms during water treatment. Practitioner Points Advanced oxidation transformations for two emerging water pollutants, bisphenol A and dinitroanisole, was investigated. The observed reaction kinetics depended on molar peroxide ratio in a manner that is in keeping with previous findings. Density functional theory‐based analysis revealed reaction energy requirements and degradation pathways.

Engineering↗
DOE OSTI2021

GPU acceleration of rank-reduced coupled-cluster singles and doubles

Here, we have developed a graphical processing unit (GPU) accelerated implementation of our recently introduced rank-reduced coupled-cluster singles and doubles (RR-CCSD) method. RR-CCSD introduces a low-rank approximation of the doubles amplitudes. This is combined with a low-rank approximation of the electron repulsion integrals via Cholesky decomposition. The result of these two low-rank approximations is the replacement of the usual fourth-order CCSD tensors with products of second- and third-order tensors. In our implementation, only a single fourth-order tensor must be constructed as an intermediate during the solution of the amplitude equations. Owing in large part to the compression of the doubles amplitudes, the GPU-accelerated implementation shows excellent parallel efficiency (95% on eight GPUs). Our implementation can solve the RR-CCSD equations for up to 400 electrons and 1550 basis functions—roughly 50% larger than the largest canonical CCSD computations that have been performed on any hardware. In addition to increased scalability, the RR-CCSD computations are faster than the corresponding CCSD computations for all but the smallest molecules. We test the accuracy of RR-CCSD for a variety of chemical systems including up to 1000 basis functions and determine that accuracy to better than 0.1% error in the correlation energy can be achieved with roughly 95% compression of the ov space for the largest systems considered. We also demonstrate that conformational energies can be predicted to be within 0.1 kcal mol -1 with efficient compression applied to the wavefunction. Finally, we find that low-rank approximations of the CCSD doubles amplitudes used in the similarity transformation of the Hamiltonian prior to a conventional equation-of-motion CCSD computation will not introduce significant errors (on the order of a few hundredths of an electronvolt) into the resulting excitation energies.

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

Customizable wave tailoring nonlinear materials enabled by bilevel inverse design

Abstract Passive wave transformation via nonlinearity is ubiquitous in settings from acoustics to optics and electromagnetics. It is well known that different nonlinearities yield different effects on propagating signals, which raises the question of “what precise nonlinearity is the best for a given wave tailoring application?” In this work, considering a one-dimensional spring-mass chain connected by polynomial springs (a variant of the Fermi-Pasta-Ulam-Tsingou system), we introduce a bilevel inverse design method which couples the shape optimization of structures for tailored constitutive responses with reduced-order nonlinear dynamical inverse design. We apply it to two qualitatively distinct problems—minimization of peak transmitted kinetic energy from impact, and pulse shape transformation—demonstrating our method’s breadth of applicability. For the impact problem, we obtain two fundamental insights. First, small differences in nonlinearity can drastically change the dynamic response of the system, from severely under- to outperforming a comparative linear system. Second, the oft-used strategy of impact mitigation via “energy locking” bistability can be significantly outperformed by our optimal nonlinearity. We validate this case with impact experiments and find excellent agreement. This study establishes a framework for broader passive nonlinear mechanical wave tailoring material design, with applications to computing, signal processing, shock mitigation, and autonomous materials.

Science & Technology - Other Topics↗
DOE OSTI2025

Multi-modal characterization of the B2 phase in the Ta-Re binary system

The energy and transportation industries demand materials that retain their mechanical property at high temperatures. Refractory complex concentrated alloys (RCCAs) with a BCC + B2 microstructure offer a potential solution, where maintaining the high temperature mechanical properties can be achieved by precipitation strengthening. This depends on the B2 phase in RCCAs being thermodynamically stable with a high solvus temperature. Recently, we predicted the high temperature stability of the B2 structure in the Ta-Re binary system, using density functional theory. Here, in this study, we provide experimental evidence for the existence of this phase for the first time, using a Ta 65 Re 35 alloy. Despite Ta-Re binary phase diagrams predicting a single-phase BCC microstructure for Ta 65 Re 35 , we show that a high Z nanoscale secondary phase appears after heat treatment at 1550 °C and 1100 °C. Scanning transmission electron microscopy (STEM) revealed that this phase has a cubic structure and is equiatomic TaRe though B2 superlattice reflections were absent in fast Fourier transforms (FFT) and diffraction patterns (DPs). DP simulations indicate that the B2 TaRe superlattice reflections are up to two orders of magnitude weaker than their fundamental reflections, making their detection challenging via electron microscopy. Neutron diffraction confirmed the second phase had a B2 structure. This study identified a previously unobserved high temperature stable B2 phase in the Ta-Re system, enabling the development of new high temperature BCC + B2 RCCAs.

36 MATERIALS SCIENCE↗