Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field 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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Higgs boson production cross-section measurements and their EFT interpretation in the $4\ell $ decay channel at $\sqrt{s}=$13 TeV with the ATLAS detector

Higgs boson properties are studied in the four-lepton decay channel (where lepton = e, μ) using 139 fb -1 of proton–proton collision data recorded at $\sqrt{s}=$13 TeV by the ATLAS experiment at the Large Hadron Collider. The inclusive cross-section times branching ratio for H→ZZ* decay is measured to be 1.34±0.12 pb for a Higgs boson with absolute rapidity below 2.5, in good agreement with the Standard Model prediction of 1.33±0.08 pb. Cross-sections times branching ratio are measured for the main Higgs boson production modes in several exclusive phase-space regions. The measurements are interpreted in terms of coupling modifiers and of the tensor structure of Higgs boson interactions using an effective field theory approach. Exclusion limits are set on the CP-even and CP-odd ‘beyond the Standard Model’ couplings of the Higgs boson to vector bosons, gluons and top quarks.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Scale-bridging with the extended/generalized finite element method for linear elastodynamics

This paper presents an extended/generalized finite element method for bridging scales in linear elastodynamics in the absence of scale separation. More precisely, the GFEMgl framework is expanded to enable the numerical solution of multiscale problems through the automated construction of specially-tailored shape functions, thereby enabling high-fidelity finite element modeling on simple, fixed finite element meshes. Furthermore, this introduces time-dependencies in the shape functions in that they are subject to continuous adaptation with time. The temporal aspects of the formulation are investigated by considering the Newmark-β time integration scheme, and the efficacy of mass lumping strategies is explored in an explicit time-stepping scheme. This method is demonstrated on representative wave propagation examples as well as a dynamic fracture problem to assess its accuracy and flexibility.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ultrafast jet classification at the HL-LHC

Abstract Three machine learning models are used to perform jet origin classification. These models are optimized for deployment on a field-programmable gate array device. In this context, we demonstrate how latency and resource consumption scale with the input size and choice of algorithm. Moreover, the models proposed here are designed to work on the type of data and under the foreseen conditions at the CERN large hadron collider during its high-luminosity phase. Through quantization-aware training and efficient synthetization for a specific field programmable gate array, we show that O ( 100 ) ns inference of complex architectures such as Deep Sets and Interaction Networks is feasible at a relatively low computational resource cost.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Error field measurements with rotating RMP fields for DIII-D H-mode

3D magnetic sensors are employed to identify the amplitude and toroidal phase of error fields (EF) by analyzing the torque balance for magnetic islands entrained by rotating resonant magnetic perturbations (RMPs) in DIII-D H-mode plasmas. This technique of torque balance allows for efficient error field identification, offering a valuable tool for scenario-specific and optimized error field compensation (EFC) and requires only magnetic diagnostics. The torque balance used in this work includes the contribution from electromagnetic torque due to error fields, wall response, RMP fields, and viscous. Results show that viscous torque plays a crucial role, particularly during locked modes and H-mode plasmas, ensuring accurate data fits with lower residuals. The torque balance technique reveals that the L- and H-mode plasmas have distinct EF configurations, and consistent with a model-based EF analysis including MHD response in IPEC and the SURFMN EF simulation. This technique shows great robustness in measuring the intrinsic EF amplitude regardless of its amplitude or toroidal phase. Repeated discharges with EFC disparities exhibit consistent results of intrinsic error field within a reasonable range near the “standard” error field compensation. Additionally, the use of a rotating n = 1 resonant magnetic perturbation offers the advantage of reducing disruption risks by entraining saturated magnetic islands. These findings are instrumental for optimizing EF correction in fusion devices, thereby enhancing tearing mode suppression and overall plasma stability.

3D magnetic sensors↗
DOE OSTI2020

ELM control optimization for various ITER scenarios based on linear and quasi-linear figures of merit

For the purpose of better understanding type-I edge localized mode (ELM) control in ITER with resonant magnetic perturbation (RMP) fields, the plasma response to RMP is computed by a resistive full magneto-hydrodynamic model in toroidal geometry. Five scenarios designed for ITER are considered, ranging from the pre-nuclear to nuclear phases. The plasma response to RMP is quantified by the plasma surface displacement near the X-point of the divertor plasma and at the outboard mid-plane. The optimal coil configurations between two high- Q deuterium-tritium (DT) scenarios (at the same plasma current of 15 MA and the same magnetic field of 5.3 T but different fusion gains, Q = 5 and 10) are predicted to be similar. For the other ITER scenarios with similar edge safety factor q 95 ~ 3 to that of the baseline scenario, the optimal coil phasing is also similar. The optimization results are different for a half-current full-field (7.5 MA/5.3 T) scenario, largely due to the difference in q 95 . The RMP coil currents are also optimized to tailor the core vs edge toroidal torques exerted by the 3D RMP fields on the plasma column. Torque optimization, with various objective functions proposed in the study, is useful for minimizing the side effects of RMP on the plasma core flow in ITER, while still maintaining the ELM control capability. Full utilization of three rows of ELM control coils in ITER is found to be essential to ensure both flexibility and robustness of ELM control, in terms of both linear and quasilinear plasma responses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Non-oscillatory no-scale inflation

We propose a non-oscillatory no-scale supergravity model of inflation (NO-NO inflation) in which the inflaton does not oscillate at the end of the inflationary era. Instead, the Universe is then dominated by the inflaton kinetic energy density (kination). During the transition from inflation to kination, the Universe preheats instantly through a coupling to Higgs-like fields. These rapidly annihilate and scatter into ultra-relativistic matter particles, which subsequently dominate the energy density, and reheating occurs at a temperature far above that of Big Bang Nucleosynthesis. After the electroweak transition, the inflaton enters a tracking phase as in some models of quintessential inflation. The model predictions for cosmic microwave background observables are consistent with Planck 2018 data, and the density of gravitational waves is below the upper bound from Big Bang Nucleosynthesis. Furthermore, we also find that the density of supersymmetric cold dark matter produced by gravitino decay is consistent with Planck 2018 data over the expected range of supersymmetric particle masses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Electric and magnetic black holes in a new nonlinear electrodynamics model

Highlights: • We introduce a new NED model which is comparable with the BI model in the weak-field limit. • We find an electric black hole solution in the context of Einstein’s gravity minimally coupled with the new NED model. • We study the thermal stability of the electric black hole solution. • Modified Smarr’s formula consistent with the first law of black hole thermodynamics is obtained. A new nonlinear electrodynamics (NED) model is introduced in the form of a nonpolynomial Lagrangian which admits static spherical and also plane wave solutions in a flat space. Upon coupling with gravity the electric field is finite and comparable with the Born–Infeld counterpart. The electric and magnetic black hole solutions in the Einstein’s gravity coupled with this NED model are presented. The solutions give both asymptotically and in the weak field limit Reissner–Nordström (RN) black hole and unlike the other known models our electric solution is expressed in terms of elementary functions in a closed form. We study the first law and derive the modified Smarr’s formula for the electric type extension of our model. Considerable rich structure, especially thermodynamic ones, ranging from first to the second order phase transitions are added to the RN black hole of linear electrodynamics with this NED model. Having the exact solution for the metric function at our disposal we investigate the stability of the electric black hole from both the thermodynamical and causal points of view.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Transient Pressure Interference During CO2 Injection in Saline Aquifers

Abstract CO2 injection in subsurface geological formations (e.g. deep saline aquifers) causes pressure perturbations over a large area surrounding the injection well. Observation wells are widely considered in geologic CO2 storage (GCS) projects where the pressure perturbation induced by CO2 injection is measured. In this work, we use analytical and numerical modeling tools along with field data to examine the pressure behavior in GCS projects before and after CO2 arrival at an observation well. Prior to CO2 arrival, a baseline pressure trend is established which corresponds to single-phase brine flow across the observation well (approximated by Theis solution). Therefore, analysis of early-time pressure data is straightforward, provides the single-phase flow characteristics (mobility and storativity), and helps establishing a baseline pressure change that can be extended beyond the single-phase flow period at the observation well. Upon CO2 arrival, a departure from this baseline trend is expected. For the pressure to detect the CO2 arrival at an observation well, the departure from baseline pressure behavior must be significant and well above the background noise levels. We use existing analytical models to determine the strength of the expected pressure departure signal from the baseline trend upon CO2 arrival. The strength of the expected pressure departure is found to be directly proportional to the mobility ratio. Accordingly, we establish a criterion to determine whether the pressure at an observation well can detect the CO2 arrival. We present an analysis approach through application to synthetic and field data and show the characteristic pressure behavior before and after CO2 arrival. We show that while generally the pressure can be either above or below the expected baseline pressure trend, it would be likely above the baseline upon CO2 arrival. This is because the mobility ratio becomes less than unity after CO2 arrival. We show that depending on the reservoir characteristics, changes in the pressure trend may or may not be sufficient to detect the CO2 arrival.

Engineering↗
DOE OSTI2022

Dynamical topological phase realized in a trapped-ion quantum simulator

Nascent platforms for programmable quantum simulation offer unprecedented access to new regimes of far-from-equilibrium quantum many-body dynamics in almost isolated systems. Here achieving precise control over quantum many-body entanglement is an essential task for quantum sensing and computation. Extensive theoretical work indicates that these capabilities can enable dynamical phases and critical phenomena that show topologically robust methods to create, protect and manipulate quantum entanglement that self-correct against large classes of errors. However, so far, experimental realizations have been confined to classical (non-entangled) symmetry-breaking orders. In this work, we demonstrate an emergent dynamical symmetry-protected topological phase, in a quasiperiodically driven array of ten 171 Yb + hyperfine qubits in Quantinuum’s System Model H1 trapped-ion quantum processor. This phase shows edge qubits that are dynamically protected from control errors, cross-talk and stray fields. Crucially, this edge protection relies purely on emergent dynamical symmetries that are absolutely stable to generic coherent perturbations. Here, this property is special to quasiperiodically driven systems: as we demonstrate, the analogous edge states of a periodically driven qubit array are vulnerable to symmetry-breaking errors and quickly decohere. Our work paves the way for implementation of more complex dynamical topological orders that would enable error-resilient manipulation of quantum information.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Magnetized laser–plasma interactions in high-energy-density systems: Parallel propagation

In this study, we investigate parametric processes in magnetized plasmas, driven by a large-amplitude pump light wave. Our focus is on laser–plasma interactions relevant to high-energy-density (HED) systems, such as the National Ignition Facility and the Sandia MagLIF concept. We present a self-contained derivation of a “parametric” dispersion relation for magnetized three-wave interactions, meaning the pump wave is included in the equilibrium, similar to the unmagnetized work of Drake et al., Phys. Fluids 17, 778 (1974). For this, we use a multi-species plasma fluid model and Maxwell's equations. The application of an external B field causes right- and left-polarized light waves to propagate with differing phase velocities. This leads to Faraday rotation of the polarization, which can be significant in HED conditions. Phase-matching and linear wave dispersion relations show that Raman and Brillouin scattering have modified spectra due to the background B field, though this effect is usually small in systems of current practical interest. We study a scattering process we call stimulated whistler scattering, where a light wave decays to an electromagnetic whistler wave (ω≲ω ce ) and a Langmuir wave. This only occurs in the presence of an external B field, which is required for the whistler wave to exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Intrinsic Thermal Hall Effect in Mott Insulators

Here, in light of recent experimental data indicating a substantial thermal Hall effect in square lattice antiferromagnetic Mott insulators, we investigate whether a simple Mott insulator can sustain a finite thermal Hall effect. We verify that the answer is “no” if one performs calculations within a spin-only low-energy effective spin model with noninteracting magnons. However, by performing determinant quantum Monte Carlo simulations, we show the single-band 𝑡−𝑡′−𝑈 Hubbard model coupled to an orbital magnetic field does support a finite thermal Hall effect when 𝑡′ ≠0 and 𝐵 ≠0 in the Mott insulating phase. We argue that the (carrier agnostic) necessary conditions for observing a finite thermal Hall effect are time-reversal and particle-hole symmetry breaking. By considering magnon-magnon scattering using a semiclassical Boltzmann analysis, we illustrate a physical mechanism by which finite transverse thermal conductivity may arise, consistent with our symmetry argument and numerical results. Our results indicate that square and triangular lattices with SU(2) symmetry can support a finite thermal Hall effect and call for a critical reexamination of thermal Hall effect data in insulating magnets, as the magnon contribution should not be excluded a priori .

Mott insulators↗
DOE OSTI2022

Deploying a Model Predictive Traffic Signal Control Algorithm - A Field Deployment Experiment Case Study: Preprint

This paper presents a field deployment experiment of a real-time traffic signal control algorithm. We implemented the Model Predictive Control (MPC) algorithm based on Virtual Phase-Link (VPL) model. We selected the deployment locations and times based on an energy saving potential concept. We developed a set of experiment systems to enable the deployment including sensing, processing, and actuating components. We tested the systems rigorously before the experiment days. We reported the key procedures on the experiment days including the steps took, real-time control procedure, and the monitoring of the experiment. We evaluated the impact of the deployment with the changes in delay and energy consumption.

deployment↗
DOE OSTI2025

Quantum stabilization of unexpected ordered phases on the honeycomb lattice

In this work, I discuss and showcase the utility of the method called minimally-augmented spin-wave theory (MAGSWT) as a relatively simple semi-analytical approach to the phase diagrams of quantum magnets. It complements numerical methods by providing physical insight into which states are competitive and by yielding approximate phase boundaries that agree well with much more numerically intensive calculations. This approach enabled me to construct the faithful phase diagrams of two paradigmatic honeycomb-lattice models in the quantum S=1/2 limit, the J 1 –J 3 FM-AF and J 1 –J 2 AF models, for the collinear quantum phases that replace or extend the classical ones. The results are in good qualitative and semi-quantitative agreement with state-of-the-art numerical studies for these models, correctly capturing the emergence of unexpected quantum phases and the suppression of classically favored spiral orders by quantum fluctuations. This study provides a much-needed important guidance to the ongoing theoretical and experimental searches of the unconventional quantum states. This work will be of significant and timely interest to both theorists and experimentalists in the field of quantum magnetism, broadly defined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗
DOE OSTI2025

Parametric raytracing modeling for NSTX-U scenario development with high harmonic fast waves and neutral beam injection

High harmonic fast waves (HHFW) are a versatile heating and current drive tool for scenario development. Extensive modeling scans were performed to find optimal parameters for different uses of HHFW in National Spherical Tokamak Experiment (NSTX-U). Scans of plasma density, temperature, magnetic field, and antenna phasing were performed both with and without neutral beam injection. For speed of calculation, the ray-tracing code GENRAY coupled to the quasilinear Fokker–Planck code CQL3D was used. CQL3D allows for a more accurate description of the fast ion population, as well as for quasilinear effects such as HHFW-induced modifications of the distribution function. Best current drive results are obtained at elevated electron temperatures and with the lowest k φ phasing. Adding neutral beams however typically strongly reduces the HHFW current drive efficiency at the low density cases due to HHFW absorption on beam ions. Results of this parametric study will feed into scenario development and predict-first whole-shot modeling of NSTX-U discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Trends in spatial correlations, two-phase coexistence, and criticality in a class of type-2 Schloegl models for autocatalysis

A class of type-2 Schloegl models is considered for particles on a square lattice with variable-range cooperativity. These models involve: (i) spontaneous particle annihilation at rate p; (ii) autocatalytic particle creation at unoccupied sites (i, j) with n ⩾ 2 particles within a specified neighborhood, Ω 𝑁 (i, j), of sites at rate $k_n$ = $\frac{^{(^n_2)}}{_{(^N_2)}}$ = $\frac{n{(n-1)}}{_{N(N-1)}}$; and (iii) possible spontaneous particle creation at unoccupied sites at “small” rate ɛ ⩾ 0. In some cases, Ω 𝑁 just includes all symmetry-equivalent sites at a single specific distance 𝑑 (in units of lattice constants) from the unoccupied site, e.g., 𝑑 = 1 (nearest-neighbor sites) where 𝑁 = 4, or 𝑑 = √5 (or √13 or…) where 𝑁 = 8. In other cases, Ω 𝑁 includes sites multiple distances from the unoccupied site, e.g., 𝑑 = {1,√2}, where 𝑁 = 8. Kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition between high- and low-density states below a critical point, ɛ < ɛ c , with generic two-phase coexistence (2PC) at least for smaller 𝑁. With some exceptions, there is an approach toward mean-field behavior with increasing 𝑁 (so the regime of generic 2PC shrinks, and ɛ c approaches the mean-field value of 1/27). Additional insight into trends is provided by analysis of the exact master equations for the models via hierarchical truncation. These truncations utilize suitably tailored pair approximations which reflect the dominant nonequilibrium spatial correlations. These correlations in turn are shown to reflect the details of the autocatalytic particle creation process. For spatially heterogeneous states, the truncations produce coupled sets of lattice differential equations (LDE) which can describe orientation-dependent propagation of an interface between high- and low-density steady states for ɛ < ɛ c . Pair approximation values of 𝑝 = 𝑝 eq where the interface is stationary, and its orientation-dependence, are in semiquantitative agreement with KMC results. In conclusion, this comparison accounts for propagation failure in the LDE which complicates interpretation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Measurements and Model Improvement: Insight into NWP Model Error Using Doppler Lidar and Other WFIP2 Measurement Systems

Abstract Doppler-lidar wind-profile measurements at three sites were used to evaluate NWP model errors from two versions of NOAA’s 3-km-grid HRRR model, to see whether updates in the latest version 4 reduced errors when compared against the original version 1. Nested (750-m grid) versions of each were also tested to see how grid spacing affected forecast skill. The measurements were part of the field phase of the Second Wind Forecasting Improvement Project (WFIP2), an 18-month deployment into central Oregon–Washington, a major wind-energy-producing region. This study focuses on errors in simulating marine intrusions, a summertime, 600–800-m-deep, regional sea-breeze flow found to generate large errors. HRRR errors proved to be complex and site dependent. The most prominent error resulted from a premature drop in modeled marine-intrusion wind speeds after local midnight, when lidar-measured winds of greater than 8 m s −1 persisted through the next morning. These large negative errors were offset at low levels by positive errors due to excessive mixing, complicating the interpretation of model “improvement,” such that the updates to the full-scale versions produced mixed results, sometimes enhancing but sometimes degrading model skill. Nesting consistently improved model performance, with version 1’s nest producing the smallest errors overall. HRRR’s ability to represent the stages of sea-breeze forcing was evaluated using radiation budget, surface-energy balance, and near-surface temperature measurements available during WFIP2. The significant site-to-site differences in model error and the complex nature of these errors mean that field-measurement campaigns having dense arrays of profiling sensors are necessary to properly diagnose and characterize model errors, as part of a systematic approach to NWP model improvement. Significance Statement Dramatic increases in NWP model skill will be required over the coming decades. This paper describes the role of major deployments of accurate profiling sensors in achieving that goal and presents an example from the Second Wind Forecast Improvement Program (WFIP2). Wind-profile data from scanning Doppler lidars were used to evaluate two versions of HRRR, the original and an updated version, and nested versions of each. This study focuses on the ability of updated HRRR versions to improve upon predicting a regional sea-breeze flow, which was found to generate large errors by the original HRRR. Updates to the full-scale HRRR versions produced mixed results, but the finer-mesh versions consistently reduced model errors.

Meteorology & Atmospheric Sciences↗