Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase space methods”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simple and Robust Carrier-Based PWM Technique for Single-Stage Three-phase Rectifier Indirect Matrix Converter

Indirect matrix converters (IMC) can offer high efficiency, high power density, and bidirectional power flow capability with a single-stage three-to single-phase AC-AC power conversion integrating a conventional three-phase active rectifier with a primary side full-bridge inverter of a dual active bridge converter. However, the IMCs generally require space vector modulation techniques, leading to complex sectors and dwell times calculations. A practical implementation of the modulation methods with a commonly used digital signal processor becomes even more challenging when considering power factor correction, harmonics reduction, and deadtime between sector transitions. This paper proposes a simple and robust carrier-based pulse width modulation (PWM) technique determined through a voltage-source inverter (VSI) to current-source inverter (CSI) to IMC conversion. This translation enables the adoption of all conventional VSI control methods with active and reactive power control in the IMCs through simple carrier and reference waveform comparison. The VSI vectors are first one-on-one matched with CSI vectors and translated to IMC modulation signals without calculating sectors and dwell times. Nine different carrier waveforms are proposed, and the converter performance is analyzed. Grid current THD and efficiency are estimated for each modulation method. The symmetric A and B methods had the lowest THD of the grid current, 2.80%, and the efficiency remained at approximately 98.3%, making it the most suitable carrier-based method proposed in this paper. The space vector PWM (SVPWM) modulation method was implemented to the DSP and generated the expected twelve gate signals.

Benson, Mikayla↗
DOE OSTI2020

Hierarchical Phenomena in Multicomponent Liquids: Simulation Methods, Analysis, Chemistry

Complex, multicomponent, solutions have often been studied solely through the lens of specific applications of interest. Yet advances to both simulation methodologies (enhanced sampling, etc.) and analysis techniques (network analysis algorithms and others), are creating a trove of data that reveal transcending characteristics across vast compositional phase space. This perspective discusses technical considerations of the reliable and accurate simulations of complex solutions, followed by the advances to analysis algorithms that elucidate coupling of different length and timescale behavior (hierarchical phenomena). The different manifestations of hierarchical phenomena are presented across an array of solution environments, emphasizing fundamental and ongoing science questions. With a more advanced molecular understanding in hand, a quintessential application (solvent extraction) is discussed, where significant opportunities exist to re-imagine the technical scope of an established technology.

Servis, Michael J.↗
DOE OSTI2021

Parameter Reduction of Composite Load Model Using Active Subspace Method

Over the past decades, the increasing penetration of distributed energy resources (DERs) has dramatically changed the power load composition in the distribution networks. The traditional static and dynamic load models can hardly capture the dynamic behavior of modern loads especially for fault-induced delayed voltage recovery (FIDVR) events. Thus, a more comprehensive composite load model with combination of static load, different types of induction motors, single-phase A/C motor, electronic load and DERs has been proposed by Western Electricity Coordinating Council (WECC). However, due to the large number of parameters and model complexity, the WECC composite load model (WECC CMLD) raises new challenges to power system studies. To overcome these challenges, in this paper, a cutting-edge parameter reduction (PR) approach for WECC CMLD based on active subspace method (ASM) is proposed. Firstly, the WECC CMLD is parameterized in a discrete-time manner for the application of the proposed method. Then, parameter sensitivities are calculated by discovering the active subspace, which is a lower-dimensional linear subspace of the parameter space of WECC CMLD in which the dynamic response is most sensitive. The interdependency among parameters can be taken into consideration by our approach. Finally, the numerical experiments validate the effectiveness and advantages of the proposed approach for WECC CMLD model.

active subspace↗
DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗
DOE OSTI2024

Dose-efficient automatic differentiation for ptychographic reconstruction

Ptychography, as a powerful lensless imaging method, has become a popular member of the coherent diffractive imaging family over decades of development. The ability to utilize low-dose X-rays and/or fast scans offers a big advantage in a ptychographic measurement (for example, when measuring radiation-sensitive samples), but results in low-photon statistics, making the subsequent phase retrieval challenging. Here, we demonstrate a dose-efficient automatic differentiation framework for ptychographic reconstruction (DAP) at low-photon statistics and low overlap ratio. As no reciprocal space constraint is required in this DAP framework, the framework, based on various forward models, shows superior performance under these conditions. It effectively suppresses potential artifacts in the reconstructed images, especially for the inherent periodic artifact in a raster scan. We validate the effectiveness and robustness of this method using both simulated and measured datasets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

An extended/generalized phase-field finite element method for crack growth with global-local enrichment

In this paper, an extended/generalized finite element method (XFEM/GFEM) for simulating quasistatic crack growth based on a phase-field method is presented. The method relies on approximations to solutions associated with two different scales: a global scale, that is, structural and discretized with a coarse mesh, and a local scale encapsulating the fractured region, that is, discretized with a fine mesh. A stable XFEM/GFEM is employed to embed the displacement and damage fields at the global scale. The proposed method accommodates approximation spaces that evolve between load steps, while preserving a fixed background mesh for the structural problem. In addition, a prediction-correction algorithm is employed to facilitate the dynamic evolution of the confined crack regions within a load step. Several numerical examples of benchmark problems in two- and three-dimensional quasistatic fracture are provided to demonstrate the approach.

multiscale↗
DOE OSTI2024

Data-driven discovery of dynamics from time-resolved coherent scattering

Coherent X-ray scattering (CXS) techniques are capable of interrogating dynamics of nano- to mesoscale materials systems at time scales spanning several orders of magnitude. However, obtaining accurate theoretical descriptions of complex dynamics is often limited by one or more factors—the ability to visualize dynamics in real space, computational cost of high-fidelity simulations, and effectiveness of approximate or phenomenological models. In this work, we develop a data-driven framework to uncover mechanistic models of dynamics directly from time-resolved CXS measurements without solving the phase reconstruction problem for the entire time series of diffraction patterns. Our approach uses neural differential equations to parameterize unknown real-space dynamics and implements a computational scattering forward model to relate real-space predictions to reciprocal-space observations. This method is shown to recover the dynamics of several computational model systems under various simulated conditions of measurement resolution and noise. Moreover, the trained model enables estimation of long-term dynamics well beyond the maximum observation time, which can be used to inform and refine experimental parameters in practice. Finally, we demonstrate an experimental proof-of-concept by applying our framework to recover the probe trajectory from a ptychographic scan. Our proposed framework bridges the wide existing gap between approximate models and complex data.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A new approach for simulating inhomogeneous chemical kinetics

Abstract In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling $$\alpha$$ α - and $$\gamma$$ γ -radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method’s efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach’s capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method’s suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Introduction to the Demonstration of Pyroshock Response Prediction Methods Using NASA ShockSat Data Virtual Collection

Mechanical shocks, particularly pyroshocks, are environments that all spacecraft and launch vehicles experience. Pyroshocks can be challenging design environments for space hardware because these events are often high frequency, high intensity transients. Furthermore, since design decisions must be made relatively early in the design phase, understanding, and accurately defining pyroshock excitation and the loads that it imparts on components is important to ensure the high reliability required of spacecraft and rockets. Here, general guidelines about shock transmission through spacecraft structures are provided in NASA and MIL handbooks and standards, but these are relatively broad guidelines, and recent studies have shown that more accurate and less-conservative predictions are needed. However, predicting shock response accurately is still a challenge because of the high frequencies in a pyroshock excitation.

Babuska, Vit [Sandia National Laboratories (SNL-NM↗
DOE OSTI2020

Bounded motion design in the Earth zonal problem using differential algebra based normal form methods

Establishing long-term relative bounded motion between orbits in perturbed dynamics is a key challenge in astrodynamics to enable cluster flight with minimum propellant expenditure. Herein, we present an approach that allows for the design of long-term relative bounded motion considering a zonal gravitational model. Entire sets of orbits are obtained via high-order Taylor expansions of Poincarè return maps about reference fixed points. The high-order normal form algorithm is used to determine a change in expansion variables of the map into normal form space, in which the phase space behavior is circular and can be easily parameterized by action–angle coordinates. The action–angle representation of the normal form coordinates is then used to parameterize the original Poincarè return map and average it over a full phase space revolution by a path integral along the angle parameterization. As a result, the averaged nodal period and drift in the ascending node are obtained, for which the bounded motion conditions are straightforwardly imposed. Sets of highly accurate bounded orbits are obtained, extending over several thousand kilometers, and valid for decades.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Synthesis, Structure, and Properties of the Complex Zintl Phase Eu 9 Zn 4.5 As 9 : A Candidate Topological Insulator and Thermoelectric Material

Reported are the synthesis and detailed analysis of the crystal and electronic structure of the novel Zintl phase Eu 9 Zn 4.5 As 9 . This material was identified in the densely populated Eu–Zn–As compositional space. For structure determination and for property measurements, suitable single crystals of this compound were grown from either Sn- or Pb-flux. Single-crystal X-ray diffraction methods indicate that Eu 9 Zn 4.5 As 9 crystallizes in the orthorhombic crystal system with the space group Pnma (a = 12.1953(7) Å, b = 4.3730(2) Å, c = 42.674(2) Å) and is formally isostructural to Ca 9 Mn 4+x Sb 9 , the less common “9–4–9” type. The structure is heavily disordered, with multiple partially occupied sites, yet, according to the Zintl-Klemm formalism, a charge-balanced composition (Eu 2+ ) 9 (Zn 2+ ) 4.5 (As 3− ) 9 is attained. Electronic structure calculations for a model, disorder-free structure indicate no energy gap between the valence and the conduction bands and suggest (semi)metallic behavior. Preliminary susceptibility measurements confirm the expected divalent nature of Eu 2+ ([Xe] 4 f 7 ground state).

Zintl phases↗
DOE OSTI2025

Data Reproducibility of Spin-Echo Small-Angle Neutron Scattering Instruments

Spin-echo small-angle neutron scattering (SESANS) is a unique method to measure structures of materials in real space with length scales from ∼ 30 μm to ∼20 µm. As shown in Figure 1, the accessible length scale of SESANS is given by its ability to encode the momentum transfer into the Larmor phase, namely Φ = $\vec{𝛿}$ ⋅ $\vec{𝑄}$, where $\vec{𝑄}$ is the momentum transfer and $\vec{𝛿}$ is the encoding vector of the setup and its projection along Q (δQ) is called spin-echo length (SEL). The spin echo length, which is synonymous with the spatial correlation distance probed, is defined as the following 𝛿𝑄 ∝ 𝜆 2 𝐵𝐿cot𝜃 where 𝜆 is neutron wavelength, B is magnetic field, L is length of the parallelogram magnetic field region, and 𝜃 is the angle between the inclined magnetic field boundary and the beam direction, as shown in Figure 1. The result of the SESANS experiment is a Hankel transformation of the SANS scattering function I(Q), which yields the correlation function of the sample in real space.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Femtosecond two-photon-resonant parametric four-wave mixing detection of atomic hydrogen

Here, we demonstrate an ultrafast parametric four-wave mixing method for the collisional-quenching-free detection of atomic hydrogen. Two photons of 243 nm light resonantly excite the n = 1 to n = 2 transition of atomic hydrogen, and a third probe photon scatters from the electronic coherence to generate the parametric emission of a degenerate signal photon. Laser-induced photofragmentation fluorescence from water, also excited with the same 243 nm pulses, is used as a convenient method to visualize beam overlap in space and time. The femtosecond time resolution enables the suppression of unwanted nonresonant contributions to the signal and separation from resonant mixing in water molecules because of the significantly different coherence lifetimes. A folded BOXCARS phase-matching configuration yields a spatially isolated, laser-like, signal beam, amenable to measurements in highly luminous backgrounds. Measurements in an H 2 /CH 4 /air flame demonstrate quantitative agreement with simulated H-atom concentration profiles. A coherence decay time constant of 42 ps is recorded in the atmospheric pressure flame.

Kliewer, Christopher Jesse [Sandia National Labora↗
DOE OSTI2025

Sparks in the dark

This study presents a novel method for the definition of signal regions in searches for new physics at collider experiments. By leveraging multi-dimensional histograms with precise arithmetic and utilizing the SparkDensityTree library, it is possible to identify high-density regions within the available phase space, potentially improving sensitivity to very small signals. Inspired by a search for dark mesons at the ATLAS experiment, CMS open data is used for this proof-of-concept intentionally targeting an already excluded signal. Signal regions are defined based on density estimates of signal and background. These preliminary regions align well with the physical properties of the signal while effectively rejecting background events.

Gudnadottir, Olga Sunneborn [Uppsala Univ. (Sweden↗
DOE OSTI2024

VOPO 4 cathode for sodium ion batteries

An electrode comprising a space group Pna2 1 VOPO 4 lattice, capable of electrochemical insertion and release of alkali metal ions, e.g., sodium ions. The VOPO 4 lattice may be formed by solid phase synthesis of KVOPO 4 , milled with carbon particles to increase conductivity. A method of forming an electrode is provided, comprising milling a mixture of ammonium metavanadate, ammonium phosphate monobasic, and potassium carbonate; heating the milled mixture to a reaction temperature, and holding the reaction temperature until a solid phase synthesis of KVOPO 4 occurs; milling the KVOPO 4 together with conductive particles to form a conductive mixture of fine particles; and adding binder material to form a conductive cathode. A sodium ion battery is provided having a conductive NaVOPO 4 cathode derived by replacement of potassium in KVOPO 4 , a sodium ion donor anode, and a sodium ion transport electrolyte. The VOPO 4 , preferably has a volume greater than 90 Å 3 per VOPO 4 .

Whittingham, M. Stanley↗
DOE OSTI2021

Influence of Wetting on Viscous Fingering Via 2D Lattice Boltzmann Simulations

We present simulations of two-phase flow using the Rothman and Keller colour gradient Lattice Boltzmann method to study viscous fingering when a “red fluid” invades a porous model initially filled with a “blue” fluid with different viscosity. We conducted eleven suites of 81 numerical experiments totalling 891 simulations, where each suite had a different random realization of the porous model and spanned viscosity ratios in the range $$M\in [0.01,100]$$ M ϵ [ 0.01 , 100 ] and wetting angles in the range $$\theta _w\in [180^\circ ,0^\circ ]$$ θ w ϵ [ 180 ° , 0 ° ] to allow us to study the effect of these parameters on the fluid-displacement morphology and saturation at breakthrough (sweep). Although sweep often increased with wettability, this was not always so and the sweep phase space landscape, defined as the difference in saturation at a given wetting angle relative to saturation for the non-wetting case, had hills, ridges and valleys. At low viscosity ratios, flow at breakthrough is localized through narrow fingers that span the model. After breakthrough, the flow field continues to evolve and the saturation continues to increase albeit at a reduced rate, and eventually exceeds 90% for both non-wetting and wetting cases. The existence of a complicated sweep phase space at breakthrough, and continued post-breakthrough evolution suggests the hydrodynamics and sweep is a complicated function of wetting angle, viscosity ratio and time, which has major potential implications to Enhanced Oil Recovery by water flooding, and hence, on estimates of global oil reserves. Validation of these results via experiments is required to ensure they translate to field studies.

Engineering↗
DOE OSTI2024

Advancing Heteroanionicity in Zintl Phases: Crystal Structures, Thermoelectric and Magnetic Properties of Two Quaternary Semiconducting Arsenide Oxides, Eu 8 Zn 2 As 6 O and Eu 14 Zn 5 As 12 O

Two novel quaternary oxyarsenides, Eu 8 Zn 2 As 6 O and Eu 14 Zn 5 As 12 O, were synthesized through metal flux reactions, and their crystal structures were established by single-crystal X-ray diffraction methods. Eu 8 Zn 2 As 6 O crystallizes in the orthorhombic space group Pbca, featuring polyanionic ribbons composed of corner-shared triangular [ZnAs 3 ] units, running along the [100] direction. The structure of Eu 14 Zn 5 As 12 O crystallizes in the monoclinic space group P2/m and its anionic substructure can be described as an infinite “ribbonlike” chain comprised of [ZnAs 3 ] trigonal-planar units, although the structural complexity here is greater and also amplified by disorder on multiple crystallographic positions. In both structures, the O 2– anion occupies an octahedral void with six neighboring Eu 2+ cations. Formal electron counting, electronic structure calculations, and transport properties reveal the charge-balanced semiconducting nature of these heteroanionic Zintl phases. High-temperature thermoelectric transport properties measurements on Eu 14 Zn 5 As 12 O reveal relatively high resistivity (ρ 500K = 8 Ω·cm) and Seebeck coefficient values (S 500K = 220 μV K –1 ), along with a low concentration and mobility of holes as the dominant charge-carriers (n 500K = 8.0 × 10 17 cm –3 , μ 500K = 6.4 cm 2 /V s). Magnetic studies indicate the presence of divalent Eu 2+ species in Eu 14 Zn 5 As 12 O and complex magnetic ordering, with two transitions observed at T 1 = 21.6 K and T 2 = 9 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗