Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space phase model”

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 595 records · Page 33

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↗

Toward more accurate adiabatic connection approach for multireference wavefunctions

A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the “fixed-RDM” approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DC-Link Capacitor Design for a Neutral-Point-Less Three-Level Dual-Phase Inverter for Traction Application

Conventional multilevel inverters, such as neutral￾point-clamped and T-type inverters, have gained popularity in electric vehicle applications due to their advantages, including high voltage range, high power capability, low switching losses, low total harmonic distortion, and low electromagnetic interference. However, these traditional multilevel inverters require a neutral point connection to generate a zero-voltage vector. The neutral point current oscillates at three times the fundamental frequency, leading to voltage imbalance and overvoltage stress on the power modules. Additionally, the use of two stacked DC-link capacitors increases the volume required for the same overall capacitance and complicates packaging due to ripple current and heat dissipation from separate components. This is a significant concern for traction drive units, where space is limited. In this paper, a neutral-point-less multilevel dual three￾phase inverter topology is investigated for traction inverter applications. Simulation results demonstrate that the proposed topology effectively retains the benefits of multilevel operation while utilizing a single DC-link capacitor. The inverter model was simulated in conjunction with an industry-standard battery model to evaluate the potential for capacitor size reduction compared to a conventional three-level inverter.

33 ADVANCED PROPULSION SYSTEMS↗

Control of Large-Amplitude Combustion Oscillations Using Nanosecond Repetitively Pulsed Plasmas

This paper details the use of nanosecond repetitively pulsed discharges to attenuate combustion instabilities in a 14 kW swirl-stabilized methane/air combustor. The combustor exhibits large-amplitude pressure oscillations ranging from 1 to 4% of the mean pressure during which the flame exhibits bulk motion in each instability cycle, upstream and downstream, as revealed by high-speed chemiluminescence. Control is accomplished with an electrode comprising a pin anode at the centerline of the combustor, allowing a nanosecond spark to be generated in a region spanning close to the flame base, through the shear layers of the swirling flow and ending at the metallic combustor wall. The discharges are generated using 20 kV, 9 kHz pulses; and they correspond to about 120 W of mean power. This results in a suppression of the peak amplitude of the pressure oscillations by a factor of two to four, and 5 dB in the rms value. Using phase-averaged visualizations of the flame with and without plasma, we detail the sequence of flame motion in the course of the instability. With the plasma active, this reveals significant interactions between the flame and the plasma during the suppression. Finally, we present a state-space model of the thermoacoustic system, and we demonstrate open-loop control of the instabilities.

Engineering↗

High-Pressure Investigation of 2,4,6-Trinitro-3-bromoanisole (TNBA): Structural Determination and Piezochromism

Understanding phase transitions in energetic materials is crucial for developing predictive models of detonation. 2,4,6-Trinitro-3-bromoanisole (TNBA), an energetic material, was studied in its single-crystal form up to pressures of 45 GPa in a diamond anvil cell. The material was characterized by using X-ray, Raman, and optical transmission measurements. From single-crystal X-ray diffraction, the ambient structure of TNBA was determined which crystallizes in the P2 1 /c space group having four molecular units per unit cell. The X-ray data up to 9.2 GPa were fitted to a third-order Birch–Murnaghan equation of state by using the parameters K o = 13.2(2.4) GPa and K p = 5.1(1.4). Between 6.8 and 7.3 GPa, a phase transition was inferred in TNBA from concurrent fading of X-ray diffraction, disappearance of Raman peaks, increase in sample fluorescence, and discontinuous color change. The new phase was consistent with an amorphous state of at least partially intact molecules judging from the presence of higher-order Raman modes and irreversibility of the Raman spectra upon release. Piezochromism was observed with the translucent yellow TNBA gradually darkening and becoming opaque black at ~25 GPa. This correlated to the absorption edge gradually shifting to the red in the visible spectrum. Signs of two possible additional structural transitions were detected in the 32.4–41.0 GPa range as suggested by a jump in the absorption edge, the irreversible changes in the absorption spectrum upon release to ambient pressure, and by the lack of Raman modes in recovered samples. Here, the crystal and electronic structures of TNBA were also investigated up to 10 GPa by using DFT calculations and crystal structure prediction (CSP) simulations. In agreement with the experimentally observed transition at 7 GPa, the simulations at 10 GPa found a bevy of polymorphs lower in enthalpy and higher in density than P2 1 /c. The lowest calculated enthalpy structure was determined to be P2 1 2 1 2 1 , being in a different space group than the ambient experimental result.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Validation of MCNP Critical Benchmark Models of PU-MET-FAST-016

A new centralized repository of high-quality MCNP models of critical benchmark experiments is currently under development at Los Alamos National Laboratory (LANL). The repository is the combined effort of the Nuclear Criticality Safety, Nuclear Data, and Monte Carlo code development/application organizations at LANL. The initial set of benchmark models in the repository are derived from the Whisper Suite provided with MCNP6.2 and are evaluated against the benchmark experiments described in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook. The goal is to build a single LANL benchmark collection that is up to date with the latest ICSBEP revision, has a formal review and revision process, is contained in an open-source repository, and utilizes new Python tools for improved input and output file review. This paper describes the validation of the models associated with PU-MET-FAST-016, “Flooded 3X3X3 Arrays of 3-kg Plutonium (Pu) Metal Cylinders – Phase I”. In this critical experiment, twenty-seven 3-kg Pu metal cylinders with a height of 46.33 ± 0.15 mm and diameter of 65.25 ± 0.05 mm were arranged in a cubic array. The Pu cylinders were sealed in aluminum cans, and three cans were placed in each of nine perforated aluminum sleeves, as shown in Fig. 1. The nine sleeves were supported by an aluminum frame that allowed for control of the lateral spacing of the cylinders, both in the X and Y directions. The benchmark experiment measured thirteen configurations. Of these configurations, six were modeled using MCNP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

Null Space Monte Carlo Evaluation of the Plateau to River Model

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model provides a computational framework to simulate the F&T of contaminants in groundwater associated with the 200-PO-1, 200-UP-1, 200-BP-5, and 200-ZP-1 Groundwater Operable Units (OUs) in the Hanford Site Central Plateau. In addition, the model includes adjacent areas and facilities (e.g., the State Approved Land Disposal Site). Intended and anticipated uses of the model include calculating water levels, hydraulic gradients, and groundwater flows throughout the model domain (encompassing the 200 West and 200 East Areas) for use in subsequent F&T calculations for contaminants of concern and developing scale-appropriate, telescopic-mesh refinement models for detailed evaluation of areas within the model domain where required. The overall objective of the modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within the OU boundaries. The objective for the model development phase is to create a common modeling platform that can be used for investigations of the Central Plateau groundwater OUs and areas downgradient toward the Columbia River. The P2R Model calibration to historical data observed at the Hanford Site is documented in CP-57037, Model Package Report for the Plateau to River Model Version 8.3. The purpose of this environmental calculation is to describe a null space Monte Carlo (NSMC) evaluation was conducted with the historic calibration of the P2R model. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. Use of the NSMC approach results in 100 groundwater flow models that are variants of the calibrated P2R Model. These variant models can be used to evaluate uncertainty in model predictions made by the calibrated P2R Model for other analyses. A secondary purpose of the environmental calculation is to establish these variant models for use with other applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Digital polycrystalline microstructure generation using diffusion probabilistic models

Accurate micromechanical simulation of polycrystalline materials requires a realistic digital representation of the grain scale microstructure. Here, this work demonstrates the use of a generative diffusion probabilistic model for synthesizing single phase polycrystalline realizations. The model performs well and is capable of producing realistic microstructures consisting of not just simple equiaxed structures but also structures exhibiting more complex spatial arrangements. Masked microstructure generation reveals that the model is context aware of morphological descriptors which may be encoded in the latent space. Training on more diverse data sets, with scaled up architectures, may enable development of future models capable of synthesizing even more complex microstructural features.

36 MATERIALS SCIENCE↗

Distilling physical origins of hardness in multi-principal element alloys directly from ensemble neural network models

Abstract Despite a plethora of data being generated on the mechanical behavior of multi-principal element alloys, a systematic assessment remains inaccessible via Edisonian approaches. We approach this challenge by considering the specific case of alloy hardness, and present a machine-learning framework that captures the essential physical features contributing to hardness and allows high-throughput exploration of multi-dimensional compositional space. The model, tested on diverse datasets, was used to explore and successfully predict hardness in Al x Ti y (CrFeNi) 1- x - y , Hf x Co y (CrFeNi) 1- x - y and Al x (TiZrHf) 1- x systems supported by data from density-functional theory predicted phase stability and ordering behavior. The experimental validation of hardness was done on TiZrHfAl x . The selected systems pose diverse challenges due to the presence of ordering and clustering pairs, as well as vacancy-stabilized novel structures. We also present a detailed model analysis that integrates local partial-dependencies with a compositional-stimulus and model-response study to derive material-specific insights from the decision-making process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DiffLense: a conditional diffusion model for super-resolution of gravitational lensing data

Abstract Gravitational lensing data is frequently collected at low resolution due to instrumental limitations and observing conditions. Machine learning-based super-resolution techniques offer a method to enhance the resolution of these images, enabling more precise measurements of lensing effects and a better understanding of the matter distribution in the lensing system. This enhancement can significantly improve our knowledge of the distribution of mass within the lensing galaxy and its environment, as well as the properties of the background source being lensed. Traditional super-resolution techniques typically learn a mapping function from lower-resolution to higher-resolution samples. However, these methods are often constrained by their dependence on optimizing a fixed distance function, which can result in the loss of intricate details crucial for astrophysical analysis. In this work, we introduce DiffLense , a novel super-resolution pipeline based on a conditional diffusion model specifically designed to enhance the resolution of gravitational lensing images obtained from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Our approach adopts a generative model, leveraging the detailed structural information present in Hubble space telescope (HST) counterparts. The diffusion model, trained to generate HST data, is conditioned on HSC data pre-processed with denoising techniques and thresholding to significantly reduce noise and background interference. This process leads to a more distinct and less overlapping conditional distribution during the model’s training phase. We demonstrate that DiffLense outperforms existing state-of-the-art single-image super-resolution techniques, particularly in retaining the fine details necessary for astrophysical analyses.

Computer Science↗

Exact Fock-State Preparation with $n^{1/4}$ Circuit Depth

Efficient, deterministic, and high-fidelity preparation of large Fock states is essential for scaling bosonic quantum technologies and exploring quantum phenomena at large excitation energies. We introduce a deterministic one-parameter (D1p) protocol that maps Fock-state preparation in an infinite-dimensional Hilbert space onto two-dimensional amplitude amplification. Starting from a coherent state with $|α|\simeq\sqrt{n}$, the initial target-state population scales as $n^{-1/2}$, yielding an iteration count and circuit depth of $\mathcal{O}(n^{1/4})$. Phase matching guarantees unit fidelity in the ideal model; remarkably, preparing $|{10^6}\rangle$ requires only 39 iterations. The protocol uses only displacements and number-selective phase operations, requires no numerical optimization, and further extends to state transfer, general superpositions, finite-dimensional systems, and multipartite entangled states. In the large-amplitude regime, its multi-target form prepares $L$-legged cat states with an iteration count determined only by $L$; cats with up to ten legs require only two iterations, independent of the coherent-state amplitude. This framework provides a broadly applicable route to highly excited bosonic states on platforms supporting these elementary controls.

Roy, Tanay [Fermilab] (ORCID:000000019442862X)↗

Filamentation and focusing of electron beams due to interactions with plasma waves

Results of numerical modeling of the interaction of an electron beam propagating across relativistic plasma waves indicate that electron beam filamentation and focusing may occur under certain conditions. The model is based on solving the relativistic equation of motion in three dimensions for the individual electrons in a tenuous Gaussian beam, as they pass through a relativistic plasma wave. Several electron beam and plasma wave parameters were varied, and the results are summarized. One of the results is that the spacing of the electron beam filaments correlates with the wavelength of the plasma wave. The electron beam filaments appear as vertical slabs after the beam exits the plasma waves. The electron beam also compresses to a focus after it exits the plasma, and the focal distance depends on several parameters including the electron beam energy and phase velocity of the relativistic plasma wave. It is suggested that these focusing and filamentation phenomena may be the basis for diagnostics schemes for laser plasma interactions. The parameters used in the model were electron beam energies in the 5–50 keV range and plasma wave properties typical for the beat-wave produced by CO2 lasers, which correspond to the facilities available in our laboratory. The limitations of these results to lower energy density beam and plasma regimes and to higher energy density regimes will be discussed.

47 OTHER INSTRUMENTATION↗

Autonomous phase mapping of gold nanoparticles synthesis with differentiable models of spectral shape

Autonomous experimentation–or self-driving labs–offers a systematic approach to accelerate materials discovery by integrating automated synthesis, characterization, and data-driven decision-making. We present a closed-loop workflow for the on-demand synthesis and structural characterization of colloidal gold nanoparticles, enabling direct mapping from composition to nanoscale structure. Our framework leverages differentiable models of spectral shape to address two central tasks in self-driving labs: (a) phase mapping, or identifying compositional regions with distinct structural behavior; and (b) material retrosynthesis, or optimizing compositions for target structure. Using functional data analysis, we develop a data-driven model with generative pre-training, active learning, and high-throughput experiments to predict spectral responses across composition space. We demonstrate the approach on seed-mediated growth of gold nanoparticles, showcasing its ability to extract design rules, reveal secondary interactions, and efficiently navigate morphology space. Gradient-based optimization of the models enables inverse design, making this a unified platform.

36 MATERIALS SCIENCE↗

Ferromagnets, a new anomaly, instantons, and (noninvertible) continuous translations

We discuss a large class of classical field theories with continuous translation symmetry. In the quantum theory, a new anomaly explicitly breaks this translation symmetry to a discrete symmetry. Furthermore, this discrete translation symmetry is extended by a d – 2-form global symmetry. All these theories can be described as U(1) gauge theories where Gauss law states that the system has nonzero charge density. Special cases of such systems can be phrased as theories with a compact phase space. Examples are ferromagnets and lattices in the lowest Landau level. In some cases, the broken continuous translation symmetry can be resurrected as a noninvertible symmetry. We clarify the relation between the discrete translation symmetry of the continuum theory and the discrete translation symmetry of an underlying lattice model. Our treatment unifies, clarifies, and extends earlier works on the same subject.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Choosing Transport Events for Initiating Splitting and Rouletting

A study was performed to determine which transport events should be used to initiate a weight window lookup to achieve the best variance reduction performance. A weight window lookup potentially triggers particle splitting (in important regions of phase space) or rouletting (in unimportant regions), thereby optimizing computational effort. Potential initiating transport events include collisions (both pre- and post-collision), geometry surface crossings, traversing a mean-free path, and streaming across a weight window boundary. Permutations of these initiating events were tested on an urban model with background radiation sources and a spent fuel cask with a neutron dose mesh tally. Generally, all methods perform better with finer weight window meshes. Tracking on weight windows performs well for coarse weight window meshes, while a combination of splitting each mean-free path, geometric surface crossing, and before collisions performs well for fine weight window meshes.

42 ENGINEERING↗

Micro-Pulse DIAL for Atmospheric Water Vapor Profiling (Phase I Final Report)

Water vapor plays a key role in Earth’s climate and weather patterns. Accurate, high resolution measurements of the water vapor vertical profiles in a region have been found to significantly improve the numerical weather prediction models’ short-term forecasts, precipitation forecasts, as well as severe weather forecasts. However, current measurement techniques provide very coarse resolution in space, time, and/or altitude. A network of continuously sampling water vapor profilers would greatly expand the scientific knowledge of the atmosphere and serve as a valuable input into numerical weather models. This would be of interest not only to the Department of Energy’s Environmental Sciences Division but to other institutions such as the National Weather Service, the National Center for Atmospheric Research, the Federal Emergency Management Agency, and the military.

54 ENVIRONMENTAL SCIENCES↗

Phase-field modeling of the interactions between an edge dislocation and an array of obstacles

Obstacles, such as voids and precipitates, are prevalent in crystalline materials. They strengthen crystals by serving as barriers to dislocation glide. Here in this work, we develop a phase-field dislocation dynamics (PFDD) technique for investigating the interactions between dislocations and second-phase obstacles, which can be either voids or precipitates. The PFDD technique is constructed to account for elastic heterogeneity, elastic anisotropy, dissociation of the dislocation, and dislocation transmission across bicrystalline interfaces. Within the framework, we present a model for “pseudo-voids”, which are voids shearable by dislocations, in contrast to unphysical, unshearable voids in conventional phase-field dislocation formulations. We employ the PFDD technique to investigate the in-plane interactions between an edge dislocation and an array of nano-scale obstacles with different spacings. In this application, the interactions take place in glide planes of either a face-centered cubic (FCC) Cu or a body-centered cubic (BCC) Nb matrix, while the precipitates have a Cu 1-x Nb x composition, with x varying from 0.1 to 0.9. Our atomistic simulations find that the alloy precipitates can have an FCC, an amorphous, or a BCC phase, depending on the compositional ratio between Cu and Nb, i.e., value of x. Among all types of obstacles, the critical stresses for dislocation bypass are the highest for unshearable amorphous precipitates, followed by shearable crystalline precipitates, and then the pseudo-voids.

42 ENGINEERING↗