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

STOMPX

STOMPX is an OpenMPI implementation of selected operational modes of the Subsurface Transport Over Multiple Phases numerical simulator. STOMPX is designed to execute with both shared- and distributed-memory computer architectures. This implementation of the simulator is designed to solve the same problems as the STOMP simulator, but taking advantage of multiple processor execution through the OpenMPI language and libraries. The STOMP simulator is written in Fortran 90 and operates on single processors or shared-memory computer architectures. The STOMPX simulator extends execution to distributed-memory computer architectures, including super computers

White, Mark↗
DOE OSTI2024

Probing New Physics with Tau Leptons using the CMS Detector

The Standard Model (SM) of particle physics is a thoroughly well-tested theory of all electroweak scale phenomena to date—the most recent test being the discovery of the long-predicted Higgs boson. However, there is much evidence to suggest that it is incomplete. On the formal side, the mass of the Higgs particle itself is perplexingly small, while the largest theoretical energy scale of the SM is 17 orders of magnitude larger, a puzzle known as the hierarchy problem. On the experimental side, among other things, the SM lacks a candidate for dark matter, which is now known to play an indispensable role in astrophysical large scale structure formation. Many theories that address the fundamental puzzles of the SM predict the existence of multiple Higgs bosons, including a very light boson that couples preferentially to third-generation fermions. Observation of this new light boson, either directly produced or in the decay of the recently discovered SM Higgs, would provide an unmistakable sign of new dynamics in nature. The research supported by this award used data collected by the Compact Muon Solenoid (CMS) detector at the Large Hadron Collider (LHC) to explore boosted di-tau signatures of a new light boson as a tool in the quest to understand the shortcomings of the SM. With a proton-proton collision energy of 13 trillion electron-volts (TeV) and peak instantaneous proton luminosity of over 10^34 collisions per square centimeter per second (cm –2 s –1 ), the LHC is the only facility in the world that could enable this line of research. The research sought to address the hierarchy problem and other shortcomings of the SM through a mixture of LHC data analysis, maintenance of the CMS detector to ensure high quality data, and contributions to the Phase 2 upgrade detector for the upcoming High Luminosity phase of the LHC (HL-LHC).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

International collaboration in disposal research: comparative modeling of coupled processes in the DECOVALEX project

Abstract. This presentation gives an overview of an international research collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems. The creation of the international DECOVALEX project, now running for more than 25 years, was motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. DECOVALEX emphasizes joint analysis and comparative modeling of state-of-the-art field and laboratory experiments, across a range of host rock options and repository designs. Participating research teams are from radioactive waste management organizations, national research institutes, regulatory agencies, universities, and consulting groups, providing a wide range of perspectives and solutions to these complex problems. The presentation provides examples of the research contributions made collectively in past DECOVALEX tasks and also touches on the unique modeling challenges tackled in the ongoing project phase, referred to as DECOVALEX-2023. The current phase comprises 17 partner organizations, about 50 modeling teams, and 7 modeling tasks, which cover a broad portfolio from fundamental studies on gas migration to full-scale in situ heater experiments in different host rocks to performance assessment studies. Together, these examples illustrate that the insight and scientific knowledge gained within the DECOVALEX project would not have been possible if one group had studied these complex THMC modeling challenges alone rather than within a truly collaborative setting.

Birkholzer, Jens T.↗
DOE OSTI2023

Elastic functional changepoint detection of climate impacts from localized sources

Detecting changepoints in functional data has become an important problem as interest in monitoring of climate phenomenon has increased, where the data is functional in nature. Here, the observed data often contains both amplitude (y-axis) and phase (x-axis) variability. If not accounted for properly, true changepoints may be undetected, and the estimated underlying mean change functions will be incorrect. In this article, an elastic functional changepoint method is developed which properly accounts for these types of variability. The method can detect amplitude and phase changepoints which current methods in the literature do not, as they focus solely on the amplitude changepoint. This method can easily be implemented using the functions directly or can be computed via functional principal component analysis to ease the computational burden. We apply the method and its nonelastic competitors to both simulated data and observed data to show its efficiency in handling data with phase variation with both amplitude and phase changepoints. We use the method to evaluate potential changes in stratospheric temperature due to the eruption of Mt. Pinatubo in the Philippines in June 1991. Using an epidemic changepoint model, we find evidence of a increase in stratospheric temperature during a period that contains the immediate aftermath of Mt. Pinatubo, with most detected changepoints occurring in the tropics as expected.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Efficient numerical integration of thermal interaction rates

In many problems in particle cosmology, interaction rates are dominated by 2 ↔ 2 scatterings, or get a substantial contribution from them, given that 1 ↔ 2 and 1 ↔ 3 reactions are phase-space suppressed. We describe an algorithm to represent, regularize, and evaluate a class of thermal 2 ↔ 2 and 1 ↔ 3 interaction rates for general momenta, masses, chemical potentials, and helicity projections. A key ingredient is an automated inclusion of virtual corrections to 1 ↔ 2 scatterings, which eliminate logarithmic and double-logarithmic IR divergences from the real 2 ↔ 2 and 1 ↔ 3 processes. We also review thermal and chemical potential induced contributions that require resummation if plasma particles are ultrarelativistic.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Finding wombling boundaries in LHC data with Voronoi and Delaunay tessellations

We address the problem of finding a wombling boundary in point data generated by a general Poisson point process, a specific example of which is an LHC event sample distributed in the phase space of a final state signature, with the wombling boundary created by some new physics. We discuss the use of Voronoi and Delaunay tessellations of the point data for estimating the local gradients and investigate methods for sharpening the boundaries by reducing the statistical noise. The outcome from traditional wombling algorithms is a set of boundary cell candidates with relatively large gradients, whose spatial properties must then be scrutinized in order to construct the boundary and evaluate its significance. Here we propose an alternative approach where we simultaneously form and evaluate the significance of all possible boundaries in terms of the total gradient flux. We illustrate our method with several toy examples of both straight and curved boundaries with varying amounts of signal present in the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

X-ray multiphoton-induced Coulomb explosion images complex single molecules

Following structural dynamics in real time is a fundamental goal towards a better understanding of chemical reactions. Recording snapshots of individual molecules with ultrashort exposure times is a key ingredient towards this goal, as atoms move on femtosecond (10 -15 s) timescales. For condensed-phase samples, ultrafast, atomically resolved structure determination has been demonstrated using X-ray and electron diffraction. Pioneering experiments have also started addressing gaseous samples. However, they face the problem of low target densities, low scattering cross sections and random spatial orientation of the molecules. Therefore, obtaining images of entire, isolated molecules capturing all constituents, including hydrogen atoms, remains challenging. Here we demonstrate that intense femtosecond pulses from an X-ray free-electron laser trigger rapid and complete Coulomb explosions of 2-iodopyridine and 2-iodopyrazine molecules. We obtain intriguingly clear momentum images depicting ten or eleven atoms, including all the hydrogens, and thus overcome a so-far impregnable barrier for complete Coulomb explosion imaging—its limitation on molecules consisting of three to five atoms. In combination with state-of-the-art multi-coincidence techniques and elaborate theoretical modelling, this allows tracing ultrafast hydrogen emission and obtaining information on the result of intramolecular electron rearrangement. Our work represents an important step towards imaging femtosecond chemistry via Coulomb explosion.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Advection algorithms for quantum neutrino moment transport

Neutrino transport in compact objects is an inherently challenging multidimensional problem. Here, this difficulty is compounded if one includes flavor transformation—an intrinsically quantum phenomenon requiring one to follow the coherence between flavors and thus necessitating the introduction of complex numbers. To reduce the computational burden, simulations of compact objects that include neutrino transport often make use of momentum-angle-integrated moments (the lowest order ones being commonly referred to as the energy density and flux) and these quantities can be generalized to include neutrino flavor, i.e., they become quantum moments. Numerous finite-volume approaches to solving the moment evolution equations for classical neutrino transport have been developed based on solving a Riemann problem at cell interfaces. In this paper we describe our generalization of a Riemann solver for quantum moments, specifically decomposing complex numbers in terms of a (signed) magnitude and phase instead of real and imaginary parts. We then test our new algorithm in numerous cases showing a neutrino fast flavor instability, varying from toy models with analytic solutions to snapshots from neutron star merger simulations. Compared to previous algorithms for neutrino transport with flavor mixing, we find uniformly smaller growth rates of the flavor transformation along with concomitantly larger length-scales, and that the results are a better match with the growth rates seen from multiangle codes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Fierro Version 2.x

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three-dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2023

Fierro

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2023

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗
DOE OSTI2021

On the Use of Field RR Lyrae as Galactic Probes. V. Optical and Radial Velocity Curve Templates

We collected the largest spectroscopic catalog of RR Lyrae (RRLs) including ≈20,000 high-, medium-, and low-resolution spectra for ≈10,000 RRLs. We provide the analytical forms of radial velocity curve (RVC) templates. These were built using 36 RRLs (31 fundamental—split into three period bins—and five first-overtone pulsators) with well-sampled RVCs based on three groups of metallic lines (Fe, Mg, Na) and four Balmer lines (H$α$, H$β$, H$γ$, H$δ$). We tackled the long-standing problem of the reference epoch to anchor light-curve and RVC templates. For the V-band, we found that the residuals of the templates anchored to the phase of the mean magnitude along the rising branch are ~35% to ~45% smaller than those anchored to the phase of maximum light. For the RVC, we used two independent reference epochs for metallic and Balmer lines and we verified that the residuals of the RVC templates anchored to the phase of mean RV are from 30% (metallic lines) up to 45% (Balmer lines) smaller than those anchored to the phase of minimum RV. We validated our RVC templates by using both the single-point and the three phase point approaches. We found that barycentric velocities based on our RVC templates are two to three times more accurate than those available in the literature. We applied the current RVC templates to Balmer lines RVs of RRLs in the globular NGC 3201 collected with MUSE at VLT. We found the cluster barycentric RV of $V$$γ$ = 496.89 ± 8.37(error) ± 3.43 (standard deviation) km s -1 , which agrees well with literature estimates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Phase Transformation Driven by Oxygen Vacancy Redistribution as the Mechanism of Ferroelectric Hf 0.5 Zr 0.5 O 2 Fatigue

As a promising candidate for nonvolatile memory devices, the hafnia-based ferroelectric system has recently been a hot research topic. Although significant progress has been made over the past decade, the endurance problem is still an obstacle to its final application. In perovskite-based ferroelectrics, such as the well-studied Pb[Zr x Ti 1-x ]O 3 (PZT) family, polarization fatigue has been discussed within the framework of the interaction of charged defects (such as oxygen vacancies) with the moving domains during the switching process, particularly at the electrode-ferroelectric interface. Armed with this background, a hypothesis is set out to test that a similar mechanism can be in play with the hafnia-based ferroelectrics. The conducting perovskite La-Sr-Mn-O is used as the contact electrode to create La 0.67 Sr 0.33 MnO 3 / Hf 0.5 Zr 0.5 O 2 (HZO)/ La 0.67 Sr 0.33 MnO 3 capacitor structures deposited on SrTiO 3 -Si substrates. Nanoscale X-ray diffraction is performed on single capacitors, and a structural phase transition from polar o-phase toward non-polar m-phase is demonstrated during the bipolar switching process. The energy landscape of multiphase HZO has been calculated at varying oxygen vacancy concentrations. Based on both theoretical and experimental results, it is found that a polar to non-polar phase transformation caused by oxygen vacancy redistribution during electric cycling is a likely explanation for fatigue in HZO.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Applying Bayesian inference and deterministic anisotropy to retrieve the molecular structure ∣Ψ(R)∣2 distribution from gas-phase diffraction experiments

Abstract Currently, our general approach to retrieving molecular structures from ultrafast gas-phase diffraction heavily relies on complex ab initio electronic or vibrational excited state simulations to make conclusive interpretations. Without such simulations, inverting this measurement for the structural probability distribution is typically intractable. This creates a so-called inverse problem. Here we address this inverse problem by developing a broadly applicable method that approximates the molecular frame structure ∣Ψ( R , t )∣ 2 distribution independent of these complex simulations. We retrieve the vibronic ground state ∣Ψ( R )∣ 2 for both simulated stretched NO 2 and measured N 2 O. From measured N 2 O, we observe 40 mÅ coordinate-space resolution from 3.75 Å −1 reciprocal space range and poor signal-to-noise, a 50X improvement over traditional Fourier transform methods. In simulated NO 2 diffraction experiments, typical to high signal-to-noise levels predict 100–1000X resolution improvements, down to 0.1 mÅ. By directly measuring the width of ∣Ψ( R )∣ 2 , we open ultrafast gas-phase diffraction capabilities to measurements beyond current analysis approaches. This method has the potential to effectively turn gas-phase ultrafast diffraction into a discovery-oriented technique to probe systems that are prohibitively difficult to simulate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Reducing the cost of energy estimation in the variational quantum eigensolver algorithm with robust amplitude estimation

Quantum chemistry and materials is one of the most promising applications of quantum computing. Yet much work is still to be done in matching industry-relevant problems in these areas with quantum algorithms that can solve them. Most previous efforts have carried out resource estimations for quantum algorithms run on large-scale fault-tolerant architectures, which include the quantum phase estimation algorithm. In contrast, few have assessed the performance of near-term quantum algorithms, which include the variational quantum eigensolver (VQE) algorithm. Recently, a large-scale benchmark study [Gonthier et al. 2020] found evidence that the performance of the variational quantum eigensolver for a set of industry-relevant molecules may be too inefficient to be of practical use. This motivates the need for developing and assessing methods that improve the efficiency of VQE. In this work, we predict the runtime of the energy estimation subroutine of VQE when using robust amplitude estimation (RAE) to estimate Pauli expectation values. Under conservative assumptions, our resource estimation predicts that RAE can reduce the runtime over the standard estimation method in VQE by one to two orders of magnitude. Despite this improvement, we find that the runtimes are still too large to be practical. These findings motivate two complementary efforts towards quantum advantage: 1) the investigation of more efficient near-term methods for ground state energy estimation and 2) the development of problem instances that are of industrial value and classically challenging, but better suited to quantum computation.

Johnson, Peter D.↗
DOE OSTI2024

Analytical Modeling of Metal Foam Composite Phase Change Materials (PCM) in Thermal Energy Storage Using Asymptotic Analysis

The use of phase change materials (PCMs) for thermal energy storage can release or absorb a significant amount of latent heat during the freezing or melting process, offering a higher energy storage density. One of the main drawbacks of PCMs is their low thermal conductivity, resulting in poor thermal performance. Recent research has attempted to enhance heat transfer and increase the thermal conductivity of PCMs, including the use of metal foams. However, modeling the metal foam composite PCM using conventional methods is computationally expensive. This paper proposes an asymptotic solution for a Stefan-like problem subject to a convective boundary for outward solidification in a hollow cylinder, capable of predicting the freeze-melt cycle of the metal foam composite PCM. Specifically, three temporal regimes and four spatial layers are considered in the asymptotic analysis for each phase change process. The thermal conductivity is calculated by a theoretical three-dimensional tetrakaidecahedron model, while other thermophysical properties are obtained using the method of volume averaging. The results are verified with numerical data and validated against experimental data in the literature. The presented analytical modeling framework could have the potential to be applied to other types of composite PCMs with considerably lower computational costs compared with conventional methods.

analytical model↗
DOE OSTI2021

Interactive Mesh Generation and Comparison Tools for Nuclear Energy Advanced Modeling & Simulation Phase I. Final Report

The promise of computer simulations for nuclear reactor designs is the opportunity to speed up the design process while also improving safety, providing more detail and higher fidelity, reducing cost, minimizing schedule risk, and avoiding many other potential problems; because simulations have the potential to be much less expensive and time-consuming than traditional development – where scale models and physical tests require construction and long lead times – it is possible to consider more alternatives and complete more comprehensive assessments before the design is finalized and constructed. As more design scenarios are studied by simulation, the speed that simulations promise is only available as long as bottlenecks are addressed for these scenarios. A frequent bottleneck is the amount of human labor required to describe the shape of the objects being designed and decompose the shape into simple elements, such as hexahedra, that are small enough to capture the physical phenomena of interest without being so small that even fast computers cannot perform the simulation in a reasonable amount of time. Recently, techniques for creating all-hexahedral decompositions have been developed, but they require some human input. Many simulation codes prefer all-hexahedral decompositions, so we proposed to evaluate the feasibility of these techniques on nuclear reactor geometries to see whether they required too much human input to be commercially viable. Our study concluded that they can be made viable with some additional software tools to reduce the amount of user input required.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Developing aqueous solubilizing agents as an alternative to solvent extraction

Here, advancing separations science is important for the entire field of chemistry. One partitioning technique that would benefit from improvement is solvent extraction. Despite its effective and widespread use, the method suffers from some problems: generation of flammable organic waste, lengthy process times, and safety concerns associated with contacting organic solvents with acidic aqueous solutions. We developed an alternative separation method inspired by solvent extraction that side-steps those issues. Toward this end, we identified that the functionality of an extractant (an agent used in solvent extraction to pull analytes from the aqueous phase into the organic phase) would change if it was modified for water solubility. In this alternative scenario, the extractant transforms into an “aqueous solubilizing agent.” We discovered that adding this aqueous solubilizing agent alongside a precipitating agent caused the contaminants to precipitate, but not the analyte. This separation concept was demonstrated within the bounds of one of the most difficult partitioning problems, separating minor actinides (Am 3+ ) from lanthanides (Ln 3+ ). We discovered that the (HSO 3 Ph) 4 BTP (aq) aqueous solubilizing agent prevented Am 3+ (aq) from precipitating with Ln 3+ (aq) when f-element precipitating agents (NaF (aq) or HF (aq) ) were added. This separation boasts impressive Am 3+ (aq) recovery yield (90 ± 2 %), near quantitative Ln 3+ (aq) removal, and high separation factors [>3000, Am 3+ (aq) vs. Nd 3+ (aq) ]. It seems likely – given the large number of candidate extractants that could be modified for aqueous solubility and the numerous combinations of existing solubilizing and precipitating agents – that this alternative approach could be used broadly in place of solvent extraction and solve other important separation problems.

(HSO3Ph)4BTP(aq)↗