Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase problem”

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

Material Model Parameters Optimization in Liquid Mercury Target Dynamics Simulation With Machine Learning Surrogates

A pulsed spallation target is subjected to very short (∼0.7μs) but intense loads (23.3 kJ) from repeated proton pulses, which knock away neutrons from the mercury atoms’ nuclei for a wide range application in physics, engineering, medicine, petroleum exploration, biology, chemistry, etc. The effect of this pulsed loading on the stainless-steel target module which contains the flowing mercury target material is difficult to predict not only due to its short but intense explosive-like physical reaction, but also the nonlinear material behavior of the liquid mercury in the structure. Injecting small helium bubbles in the mercury has been an efficient method of mitigating the pressure wave at high power level stage. However, prediction of the resultant loading on the target is more difficult when helium gas is intentionally injected into the mercury. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. A parameter sensitivity study was firstly employed to understand their impact on the simulation strains. The investigated parameters included E, μ, γ, σ, n, VFgas, and gas cumulative volume curve control parameters a and b. Verification and validation results from sparse polynomial expansions (SPE) method and directional Gaussian smoothing (DGS) optimization show that the surrogate model had training error of ∼7% and validation error of ∼15%, indicating that machine learning methods and surrogate models can help optimize the uncertain parameters in the complex 2-phase material model. This approach is expected to fill the knowledge gap between unknown liquid-gas mixture material model and measured vessel strain responses.

Lin, Lianshan↗
DOE OSTI2024

Learnability Transitions in Monitored Quantum Dynamics via Eavesdropper’s Classical Shadows

Monitored quantum dynamics—unitary evolution interspersed with measurements—has recently emerged as a rich domain for phase structure in quantum many-body systems away from equilibrium. Here we study monitored dynamics from the point of view of an eavesdropper who has access to the classical measurement outcomes, but not to the quantum many-body system. We show that a measure of information flow from the quantum system to the classical measurement record—the —undergoes a phase transition in correspondence with the measurement-induced phase transition (MIPT). This transition determines the eavesdropper’s (in)ability to learn properties of an unknown initial quantum state of the system, given a complete classical description of the monitored dynamics and arbitrary classical computational resources. We make this learnability transition concrete by defining classical shadow protocols that the eavesdropper may apply to this problem, and show that the MIPT manifests as a transition in the sample complexity of various shadow-estimation tasks, which become harder in the low-measurement phase. We focus on three applications of interest: Pauli expectation values (where we find the MIPT appears as a point of optimal learnability for typical Pauli operators), many-body fidelity, and global charge in U ( 1 ) -symmetric dynamics. Our work unifies different manifestations of the MIPT under the umbrella of and gives this notion a general operational meaning via classical shadows. Published by the American Physical Society 2024

Physics↗
DOE OSTI2025

Quench process modeling & simulation in the heat-treatment of critical aerospace components

To manufacture light-weight, advanced metal alloy components for gas turbine engines, quench heat-treatment processes are typically used. By quenching the component from elevated temperatures, the alloy sometimes undergoes a solid-state phase transformation which produces special microstructures with the required, enhanced mechanical properties. However, the quenching can also lead to cracks forming in the component. Addressing the quench cracking problems adds a significant burden to the cost, schedule, and energy demand of manufacture. Currently, optimizing the quench process to mitigate or avoid the cracking is performed largely by trial-and-error, relying heavily on costly experimental (thermocouple)trials to understand the local thermal gradients which cause the cracks to form. In this first part (Phase 1) of the work, high-performance computing is employed to establish the ability of modern CFD (computational fluid dynamics) to alleviate or wholly replace the experimental quenching trials by virtual testing. A Baseline CFD model is defined and its accuracy established to be comparable to(and which usually exceeds) the accuracy of existing HTC (heat-transfer coefficient) based simulation methods of quenching. As a first-principles based approach, “calibration” of the Baseline CFD model is independent of the quench process itself, but instead relies on the accuracy of the underlying (modeled),generic two-phase fluid processes which cannot be currently resolved by CFD for large, industrial-scale cases. A novel, high-fidelity DNS capability has been developed and verified to examine and further improve upon the mean-field closure submodels on which the Baseline CFD approach is based.1

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Quench process modeling & simulation in the heat-treatment of critical aerospace components

To manufacture light-weight, advanced metal alloy components for gas turbine engines, quench heat-treatment processes are typically used. By quenching the component from elevated tempera-tures, the alloy sometimes undergoes a solid-state phase transformation which produces special microstructures with the required, enhanced mechanical properties. However, the quenching can also lead to cracks forming in the component. Addressing the quench cracking problems adds a significant burden to the cost, schedule, and energy demand of manufacture. Currently, optimizing the quench process to mitigate or avoid the cracking is performed largely by trial-and-error, relying heavily on costly experimental (thermocouple) trials to understand the local thermal gradients which cause the cracks to form. In this first part (Phase 1) of the work, high-performance computing is employed to establish the ability of modern CFD (computational fluid dynamics) to alleviate or wholly replace the experimental quenching trials by virtual testing. A baseline CFD model is defined and its accuracy established to be comparable to (and which usually exceeds) the accuracy of existing HTC (heat-transfer correlation) based simulation methods of quenching. As a first-principles based approach, “calibration” of the Baseline CFD model is independent of the quench process itself, but instead relies on the accuracy of the underlying (modeled), generic two-phase fluid processes which cannot be currently resolved by CFD for large, industrial-scale cases. A novel, high-fidelity DNS capability has been developed and verified to examine and further improve upon the mean-field closure submodels on which the Baseline CFD approach is based.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Alamo multiphysics solver for phase field simulations with strong-form mechanics and block structured adaptive mesh refinement

Alamo is a high-performance scientific code that uses block-structured adaptive mesh refinement to solve such problems as: the ignition and burn of solid rocket propellant, plasticity, damage and fracture in materials undergoing loading, and the interaction of compressible flow with eroding solid materials. Alamo is powered by AMReX, and provides a set of unique methods, models, and algorithms that enable it to solve solid-mechanics problems (coupled to other physical behavior such as fluid flow or thermal diffusion) using the power of block-structured adaptive mesh refinement.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Towards Quantum Computing Phase Diagrams of Gauge Theories with Thermal Pure Quantum States

The phase diagram of strong interactions in nature at finite temperature and chemical potential remains largely theoretically unexplored due to inadequacy of Monte-Carlo–based computational techniques in overcoming a sign problem. Quantum computing offers a sign-problem-free approach, but evaluating thermal expectation values is generally resource intensive on quantum computers. To facilitate thermodynamic studies of gauge theories, we propose a generalization of the thermal-pure-quantum-state formulation of statistical mechanics applied to constrained gauge-theory dynamics, and numerically demonstrate that the phase diagram of a simple low-dimensional gauge theory is robustly determined using this approach, including mapping a chiral phase transition in the model at finite temperature and chemical potential. Quantum algorithms, resource requirements, and algorithmic and hardware error analysis are further discussed to motivate future implementations. Thermal pure quantum states, therefore, may present a suitable candidate for efficient thermal simulations of gauge theories in the era of quantum computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nonlinear kinetic simulation study of the ion–ion streaming instability in single- and multi-ion species plasmas

The nonlinear evolution of the ion–ion streaming instability (IISI) is studied using numerical techniques novel to this problem that afford direct insight into the evolution of the particle distributions of each species. Here, during the linear phase of the instability, we demonstrate quantitative agreement with linear kinetic theory. Subsequently, the electrostatic field generated by the IISI causes ring-like velocity distributions of ions to form that are both heated and slowed to varying degrees relative to their initial flows. Due to variation in the trapping conditions for ion species of differing charge-to-mass ratio, when flows of multiple species interact, the nonlinear evolution of each species can be starkly different: we show a case where a lighter ion species is completely stopped by a heavier ion species via the IISI alone (i.e., without collisions) and, for the first time, demonstrate how the IISI can introduce a relative flow between ion species that initially have the same flow velocities, thereby separating them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Diffusion bonding of tungsten-vanadium-zirconium using vacuum hot pressing for the development of a low decay heat cladding solution for tungsten spallation targets

Tantalum has been used as cladding material for water-cooled solid tungsten targets at many leading spallation neutron production facilities thanks to its high neutron yield, manageable radiation damage behavior, and excellent corrosion/erosion resistance in radiation environments. However, from a safety hazard perspective, thermal neutron capture of tantalum in spallation environments causes a high specific decay heat in the target volume, which often becomes a limiting factor in increasing the beam power on the target. In this paper, we studied vacuum hot pressing (VHP) parameters to diffusion bond zirconium to tungsten to explore the feasibility of using zirconium alloys as an alternative cladding material to tantalum. Zirconium alloys have long been used as cladding material for early generation solid spallation targets, and nuclear fuel rods. In spallation environments zirconium has significantly lower decay heat with shorter decay time compared to tantalum. The hot isostatic pressing (HIP) of zirconium and tungsten is known to produce limited bonding quality due to the formation of the brittle ZrW2 intermetallic layer. To overcome this problem, placing a vanadium interlayer between tungsten and zirconium has been proposed by exploring parameter space in binary alloy phase diagrams. Under the VHP conditions, 860 ◦C at 70 MPa for 4 hours, Zr-V and V-W showed good diffusion bonding, which demonstrates the feasibility of a single step HIP process to make the zirconium alloy clad tungsten spallation volumes.

Lee, Yong Joong↗
DOE OSTI2023

High temperature lightweight Al—Fe—Si based alloys

Described herein are approaches to stabilizing AlFeSi ternary intermetallic compounds while destabilizing competing phases. The inclusion of metals such as Mn, Ni, Co, Cu, or Zn to produce quaternary systems accomplishes this problem associated with AlFeSi ternary intermetallic compounds.

Soto Medina, Sujeily↗
DOE OSTI2020

A conservative phase-space moving-grid strategy for a 1D-2V Vlasov–Fokker–Planck Solver

In this work, we develop a conservative configuration- and velocity-space (i.e., phase-space) moving-grid strategy for the Vlasov–Fokker–Planck (VFP) equation in a planar geometry. The velocity-space grid is normalized and shifted in terms of the thermal speed and the bulk-fluid velocity, respectively. The configuration-space grid is moved according to a mesh-motion-partial-differential equation (MMPDE), which equidistributes a monitor function that is inversely proportional to the gradient-length scales of the macroscopic plasma quantities. The resulting inertial terms in the transformed VFP equations are discretized to ensure the discrete conservation of mass, momentum, and energy. To satisfy the discrete conservation theorems in the presence of phase-space mesh motion, we employ the method of discrete nonlinear constraints – explored in previous studies – but the underlying symmetries are determined in a much more efficient manner than before. The conservative grid-adaptivity strategy provides an efficient scheme that resolves important physical structures in the phase-space while controlling the computational complexity at all times. We demonstrate the favorable features of the proposed algorithm through a set of test cases of increasing complexity. The problems test independent components of the algorithms, as well as the integrated capability on settings relevant to inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Gradient-based constrained optimization using a database of linear reduced-order models

A methodology grounded in model reduction is presented for accelerating the gradient-based solution of a family of linear or nonlinear constrained optimization problems where the constraints include at least one linear Partial Differential Equation (PDE). A key component of this methodology is the construction, during an offline phase, of a database of pointwise, linear, Projection-based Reduced-Order Models (PROM)s associated with a design parameter space and the linear PDE(s). A parameter sampling procedure based on an appropriate saturation assumption is proposed to maximize the efficiency of such a database of PROMs. A real-time method is also presented for interpolating at any queried but unsampled parameter vector in the design parameter space the relevant sensitivities of a PROM. The practical feasibility, computational advantages, and performance of the proposed methodology are demonstrated for several realistic, nonlinear, aerodynamic shape optimization problems governed by linear aeroelastic constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A numerical-homogenization based phase-field fracture modeling of linear elastic heterogeneous porous media

Most porous media, such as geomaterials and biomaterials are highly heterogeneous in nature, and they contain large variations of microscopic pore structures, such as pore sizes, pore distribution, and pore shapes. The oscillation of microscopic structures is a substantial challenge in theoretical characterization and is usually ignored in continuous modeling. However, mechanical behavior of porous media such as deformation and failure, are essentially impacted by the microscopic heterogeneity which needs to be considered in modeling a porous media. Here, this research proposes a numerical modeling framework with a capability to investigate the effect of microscopic heterogeneity on the macroscopic fracture behavior in porous media by using a numerical homogenization technique, combined with the phase-field fracture modeling method. This numerical modeling strategy computes a homogenized elasticity tensor based on microscopic heterogeneous pore structures heterogenous porous domain by solving boundary value problems at microscopic domain. The strain energy and subsequent propagation of macroscopic fractures will be updated using homogenized stiffness information. Using this numerical scheme, the microscopic pore structure’s impact on the fracture behavior through the homogenized elastic tensor will be taken into account. This multiscale technique is benchmarked against classical problems. Finally, the results highlight the importance of the underlying pore structure and reveal that both fracture strength and propagation path can be influenced by the microscopic heterogeneity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spatial Mode Correction of Single Photons Using Machine Learning

Abstract Spatial modes of light constitute valuable resources for a variety of quantum technologies ranging from quantum communication and quantum imaging to remote sensing. Nevertheless, their vulnerabilities to phase distortions, induced by random media, impose significant limitations on the realistic implementation of numerous quantum‐photonic technologies. Unfortunately, this problem is exacerbated at the single‐photon level. Over the last two decades, this challenging problem has been tackled through conventional schemes that utilize optical nonlinearities, quantum correlations, and adaptive optics. In this article, the self‐learning and self‐evolving features of artificial neural networks are exploited to correct the complex spatial profile of distorted Laguerre–Gaussian modes at the single‐photon level. Furthermore, the potential of this technique is used to improve the channel capacity of an optical communication protocol that relies on structured single photons. The results have important implications for real‐time turbulence correction of structured photons and single‐photon images.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Disentangling observable dependence in SCET I and SCET II anomalous dimensions: angularities at two loops

The resummation of radiative corrections to collider jet observables using soft collinear effective theory is encoded in differential renormalization group equations (RGEs), with anomalous dimensions depending on the observable under consideration. This observable dependence arises from the ultraviolet (UV) singular structure of real phase space integrals in the effective field theory. We show that the observable dependence of anomalous dimensions in SCET I problems can be disentangled by introducing a suitable UV regulator in real radiation integrals. Resummation in the presence of the new regulator can be performed by solving a two-dimensional system of RGEs in the collinear and soft sectors, and resembles many features of resummation in SCET II theories by means of the rapidity renormalization group. We study the properties of SCET I with the additional regulator and explore the connection with the system of RGEs in SCET II theories, highlighting some universal patterns that can be exploited in perturbative calculations. As an application, we compute the two-loop soft and jet anomalous dimensions for a family of recoil-free angularities and give new analytic results. This allows us to study the relations between the SCET I and SCET II limits for these observables. We also discuss how the extra UV regulator can be exploited to calculate anomalous dimensions numerically, and the prospects for numerical resummation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Atomic-Scale Structural Mapping of Active Sites in Monolayer PGM-Free Catalysts by Low-Voltage 4D-STEM

Two-dimensional (2D) materials have attracted a large amount of attention in both basic and applied fields, and scanning transmission electron microscopy (STEM) is often uniquely well-suited for characterizing the atomic-scale structure of these materials [1-4]. As a result, STEM is poised to significantly impact progress on platinum group metal (PGM)-free catalysts, which are currently under intense development to enable low-cost, commercially viable hydrogen fuel cells [5]. While recent advancements have resulted in fuel cell performance comparable to Pt catalysts by some measures [6], cell durability remains a significant challenge, limiting practical applications [7]. Catalytically active sites in PGM-free materials are proposed to be FeN4 complexes embedded in a graphene lattice (Fig. 1b) within layered or other larger materials, but this is still under debate largely due to the range of potential actives sites predicted by computational methods and lack of methods for directly validating these models [5]. Fundamental insights into the atomic structure and resulting degradation pathways of proposed active sites are therefore needed to fully understand and control cell performance and durability [6].2D materials typically make ideal samples for STEM, but those within PGM-free catalysts present additional challenges since these materials are often defect-rich, with a high density of edges, dopant atoms, etc., which significantly increase susceptibility to beam damage at standard operating voltages. This makes analysis of potential FeN4 active sites particularly challenging, since a large proportion of Fe exists at edge sites where beam-induced atomic displacements can prohibit high-resolution structural characterization [6]. Conventional dark-field imaging compounds this problem by producing less signal for a given dose and being less sensitive to light elements than dose-efficient phase contrast imaging techniques such as those enabled by four-dimensional (4D)-STEM [8-10] (Fig. 1a). Consequently, active site structural analysis is often left to methods such as low-resolution imaging combined with quantum chemical calculations [6], which hinders accurate determination of reaction and degradation mechanisms.Here, we demonstrate direct atomic-scale structural mapping of FeN4 sites by performing low-voltage 4D-STEM on a model PGM-free catalyst system with many exposed monolayer regions. To accomplish this, we pair a 30 keV aberration-corrected probe with a fast pixelated detector that has optimal performance at low beam voltages [11]. This enables us to simultaneously image light and heavy elements with high signal-to-noise by center-of-mass analysis (Fig. 1c) while minimizing beam-induced atomic displacements at sensitive sites. The monolayer nature of these materials additionally allows for experimental validation by direct comparison with multislice simulations [12] of model structures (Fig. 1d-e). This work demonstrates how low-voltage 4D-STEM will provide new insights into the atomic-scale structure and degradation mechanisms of active sites in PGM-free catalysts, facilitating the development of low-cost hydrogen fuel cells and other energy conversion technologies in the future [13].

Zachman, Michael↗
DOE OSTI2021

Reversing the Irreversible: Thermodynamic Stabilization of LiAlH 4 Nanoconfined Within a Nitrogen-Doped Carbon Host

A general problem when designing functional nanomaterials for energy storage is the lack of control over the stability and reactivity of metastable phases. Using the high-capacity hydrogen storage candidate LiAlH4 as an exemplar, we demonstrate an alternative approach to the thermodynamic stabilization of metastable metal hydrides by coordination to nitrogen binding sites within the nanopores of N-doped CMK-3 carbon (NCMK-3). The resulting LiAlH 4 @NCMK-3 material releases H 2 at temperatures as low as 126 °C with full decomposition below 240 °C, bypassing the usual Li 3 AlH 6 intermediate observed in bulk. Moreover, >80% of LiAlH 4 can be regenerated under 100 MPa H 2 , a feat previously thought to be impossible. Nitrogen sites are critical to these improvements, as no reversibility is observed with undoped CMK-3. Density functional theory predicts a drastically reduced Al–H bond dissociation energy and supports the observed change in the reaction pathway. Finally, the calculations also provide a rationale for the solid-state reversibility, which derives from the combined effects of nanoconfinement, Li adatom formation, and charge redistribution between the metal hydride and the host.

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

Unified origin of inflation, baryon asymmetry, and neutrino mass

In this work, we present a unified theoretical framework that simultaneously addresses some of the most intriguing puzzles in particle physics and cosmology, namely the origins of neutrino mass, baryon asymmetry, and cosmic inflation. In our model, inflation is driven by a combination of the Standard Model Higgs, the type II seesaw Higgs responsible for neutrino mass generation, and the unified symmetry-breaking Higgs field. During inflation, nonzero values of the latter field ensure the absence of the monopole problem. The baryon asymmetry is generated through the Affleck-Dine mechanism, facilitated by the nonzero angular motion in the phase of a complex scalar field, which is part of the inflaton. We find that the successful parameter region for generating baryon asymmetry through a renormalizable term in the scalar potential requires a rather heavy type II seesaw triplet, with a mass well beyond the TeV scale. Inflationary observables, in particular, the spectral index is in excellent agree with experimental observation, whereas tensor-to scalar ratio is expected to be probed by the future LiteBIRD and CMB-S4 missions. Published by the American Physical Society 2024

Kaladharan, Ajay (ORCID:0000000165096427)↗