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 1,063 records · Page 59

Multimodal Approach Reveals the Symmetry-Breaking Pathway to the Broken Helix in EuIn 2 ⁢As 2

Understanding and manipulating emergent phases, which are themes at the forefront of quantum-materials research, rely on identifying their underlying symmetries. This general principle has been particularly prominent in materials with coupled electronic and magnetic degrees of freedom, in which magnetic order influences the electronic band structure and can lead to exotic topological effects. However, identifying symmetry of a magnetically ordered phase can pose a challenge, particularly in the presence of small domains. Here we introduce a multimodal approach for determining magnetic structures, which combines symmetry-sensitive optical probes, scattering, and group-theoretical analysis. We apply it to EuIn 2 ⁢As 2 , a material that has received attention as a candidate axion insulator. While first-principles calculations predict this state on the assumption of a simple collinear antiferromagnetic structure, subsequent neutron-scattering measurements reveal a much more intricate magnetic ground state characterized by two coexisting magnetic wave vectors reached by successive thermal phase transitions. The proposed high- and low-temperature phases are a spin helix and a state with interpenetrating helical and Néel antiferromagnetic order termed a “broken helix,” respectively. Employing a multimodal approach, we identify the magnetic structure associated with these two phases of EuIn 2 ⁢As 2 . We find that the higher-temperature phase is characterized by a variation of the magnetic moment amplitude from layer to layer, with the moment vanishing entirely in every third Eu layer. The lower-temperature structure is similar to the broken helix, with one important difference: Because of local strain, the relative orientation of the magnetic structure and the lattice is not fixed. Consequently, the symmetry required to protect the axion phase is not generically protected in EuIn 2 ⁢As 2 , but we show that it can be restored if the magnetic structure is tuned with uniaxial strain. Finally, we present a spin Hamiltonian that identifies the spin interactions that account for the complex magnetic order in EuIn 2 ⁢As 2 . Our work highlights the importance of a multimodal approach in determining the symmetry of complex order parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multiple critical points in effective quark models

We consider the two-flavor version of the linear sigma model as well as of the Nambu-Jona-Lasinio model, at finite temperature and quark chemical potential, beyond the mean field approximation. Using parameter values for the pion and quark current masses which weakly break chiral symmetry, we show that both models can present more than one critical end point. In particular, we explicitly show that the appearance of a new critical point associated with a first-order line at high temperature and low densities could help to conciliate some lattice results with model predictions. Using different techniques, we perform an extensive thermodynamical analysis to understand the physical nature of the different critical points. For both models, our results suggest that the new first-order line which starts at vanishing chemical potential has a more chiral character than the usual line which displays a character more reminiscent of a liquid-gas phase transition.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Gravitational waves from first-order phase transition in a simple axion-like particle model

We consider a gauge-singlet complex scalar field Φ with a global U(1) symmetry that is spontaneously broken at some high energy scale ƒ a . As a result, the angular part of the Φ-field becomes an axion-like particle (ALP). We show that if the Φ-field has a nonzero coupling κ to the Standard Model Higgs boson, there exists a certain region in the (ƒ a , κ) parameter space where the global U(1) symmetry-breaking induces a strongly first order phase transition, thereby producing stochastic gravitational waves that are potentially observable in current and future gravitational-wave detectors. In particular, we find that future gravitational-wave experiments such as TianQin, BBO and Cosmic Explorer could probe a broad range of the energy scale 10 3 GeV ≳ ƒ a ≳ 108 GeV, independent of the ALP mass. Since all the ALP couplings to the Standard Model particles are proportional to inverse powers of the energy scale ƒ a (up to model-dependent O (1) coefficients), the gravitationalwave detection prospects are largely complementary to the current laboratory, astrophysical and cosmological probes of the ALP scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗

Particle Production in Nucleus-Nucleus Collisions at the SPS and the QCD Phase Diagram

Results of the NA49 experiment on particle production in the energy regime of the CERN-SPS are reviewed. In order to collect information on the properties of the QCD phase diagram systematic studies of the energy dependence of particle production in heavy ion collisions have been performed. Net-baryon distributions and results on strangeness production are discussed. Recent attempts to search for a critical point, such as multiplicity fluctuations and the transverse mass dependence of anti-baryon/baryon ratios are reviewed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Superfluid, solid, and supersolid phases of dipolar bosons in a quasi-one-dimensional optical lattice

We discuss a model of dipolar bosons trapped in a weakly coupled planar array of one-dimensional tubes. We consider the situation where the dipolar moments are aligned by an external field, and we find a rich phase diagram as a function of the angle of this field exhibiting quantum phase transitions between solid, superfluid, and supersolid phases. In the low energy limit, the model turns out to be identical to one describing quasi-one-dimensional superconductivity in condensed matter systems. This opens the possibility of using bosons as a quantum analog simulator of electronic systems, a scenario arising from the intricate relation between statistics and interactions in quasi-one-dimensional systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

Combustion of 10–100 μm aluminum droplets in detonation products gases

Here we describe a two-phase model of combustion effects in aluminized high explosive (HE) charges. It is based on: (i) a Gas Dynamic Model of the expansion of the detonation product gases and their turbulent combustion with air; and (ii) a Heterogeneous Continuum Model of aluminum (Al) droplets and their combustion with the detonation product gases. Initial conditions are based on an analytical similarity solution for a cylindrical Chapman-Jouguet (CJ) detonation propagating at the CJ detonation velocity. The CJ jump conditions are computed at the thermodynamic equilibrium state by the Cheetah code, assuming the Al droplets are inert. We assume that the Al is 10 % of the charge mass and occurs as droplets at the CJ state. Different initial droplet diameters, ranging from 10 to 100 microns, are studied. A hydrodynamic combustion model based on large Damköhler numbers is employed in this study, -3σ(1 + 0.276 $\sqrt{Re}$)/($Kd^2_w$). It has a square-root dependence on the Reynolds number (Re) and inversed-squared dependence on the droplet diameter (d w ). The burnout time (t B ) of the Al droplets has a three-halves dependence on the droplet diameter, t B ~ ($d^0_w$) 3/2 . After burnout, the detonation products act as detonation products of the HE charge with active Al. They turbulently mix with air and form a combustion layer on the outer edge of the fireball. Details of the two-phase model, initial conditions and evolution of the flow field will be described.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inference on accretion flow dynamics using TCAF solution from the analysis of spectral evolution of H 1743-322 during the 2010 outburst

We study accretion flow dynamics of the Galactic transient black hole candidate (BHC) H 1743-322 during its 2010 outburst by analyzing spectral data using the two-component advective flow (TCAF; Keplerian and sub-Keplerian) solution after its inclusion in XSPEC as a local model. We compare our TCAF solution fitted results with combined disk blackbody (DBB) and power-law (PL) model fitted results and find a similar smooth variation of thermal (Keplerian or DBB) and non-thermal (PL or sub-Keplerian) fluxes/rates in two types of model fits. For a spectral analysis, 2.5-25 keV spectral data from the Rossi X-Ray Timing Explorer Proportional Counter Array instrument are used. From the TCAF solution fit, accretion flow parameters, such as Keplerian rate, sub-Keplerian rate, location of centrifugal pressure-supported shock, and strength of the shock, are extracted, providing a deeper understanding of the accretion process and properties of accretion disks around BHC H 1743-322 during its X-ray outburst. Based on the halo to disk accretion rate ratio, shock properties, accretion rates, and the nature of the quasi-periodic oscillations' (if observed) entire outburst is classified into four different spectral states: hard, hard-intermediate, soft-intermediate, and soft. From the time variation of intrinsic flow parameters, it appears that their evolutions in the declining phase do not retrace the path of the rising phase. Since our current model does not include magnetic fields, spectral turnover at energies beyond 500-600 keV cannot be explained.

79 ASTRONOMY AND ASTROPHYSICS↗

Discovery of multiple p-mode pulsation frequencies in the roAp star, HD 86181

ABSTRACT We report the frequency analysis of a known roAp star, HD 86181 (TIC 469246567), with new inferences from Transiting Exoplanet Survey Satellite (TESS) data. We derive the rotation frequency to be νrot = 0.48753 ± 0.00001 d−1. The pulsation frequency spectrum is rich, consisting of two doublets and one quintuplet, which we interpret to be oblique pulsation multiplets from consecutive, high-overtone dipole, quadrupole, and dipole modes. The central frequency of the quintuplet is 232.7701 d−1 (2.694 mHz). The phases of the sidelobes, the pulsation phase modulation, and a spherical harmonic decomposition all show that the quadrupole mode is distorted. Following the oblique pulsator model, we calculate the rotation inclination, i, and magnetic obliquity, β, of this star, which provide detailed information about the pulsation geometry. The i and β derived from the best fit of the pulsation amplitude and phase modulation to a theoretical model, including the magnetic field effect, slightly differ from those calculated for a pure quadrupole, indicating the contributions from ℓ = 4, 6, 8,... are small. Non-adiabatic models with different envelope convection conditions and physics configurations were considered for this star. It is shown that models with envelope convection almost fully suppressed can explain the excitation at the observed pulsation frequencies.

Shi, Fangfei (ORCID:0000000323998216)↗

Total absorption spectroscopy of the β decay of 101,102 Zr and 109 Tc

The β decay of 101,102 Zr and 109 Tc was studied using the technique of total absorption spectroscopy. The experiment was performed at the National Superconducting Cyclotron Laboratory using the Summing NaI(Tl) (SuN) detector in the first-ever application of total absorption spectroscopy with a fast beam produced via projectile fragmentation. The β -decay feeding intensity and Gamow-Teller transition strength distributions were extracted for these three decays. The extracted distributions were compared to three different quasiparticle random-phase approximation (QRPA) models based on different mean-field potentials. A comparison with calculations from one of the QRPA models was performed to learn about the ground-state shape of the parent nucleus. For 101 Zr and 102 Zr, calculations assuming a pure shape configuration (oblate or prolate) were not able to reproduce the extracted distributions. These results may indicate that some type of mixture between oblate and prolate shapes is necessary to reproduce the extracted distributions. For 109 Tc, a comparison of the extracted distributions with QRPA calculations suggests a dominant oblate configuration. The other two QRPA models are commonly used to provide β-decay properties in r-process network calculations. This work shows the importance of making comparisons between the experimental and theoretical β-decay distributions, rather than just half-lives and β-delayed neutron emission probabilities, as close to the r-process path as possible.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark photon superradiance: Electrodynamics and multimessenger signals

We study the electrodynamics of a kinetically mixed dark photon cloud that forms through superradiance around a spinning black hole, and design strategies to search for the resulting multimessenger signals. A dark photon superradiance cloud sources a rotating dark electromagnetic field which, through kinetic mixing, induces a rotating visible electromagnetic field. Standard model charged particles entering this field initiate a transient phase of particle production that populates a plasma inside the cloud and leads to a system which shares qualitative features with a pulsar magnetosphere. We study the electrodynamics of the dark photon cloud with resistive magnetohydrodynamics methods applicable to highly magnetized plasma, adapting techniques from simulations of pulsar magnetospheres. We identify turbulent magnetic field reconnection as the main source of dissipation and electromagnetic emission, and compute the peak luminosity from clouds around solar-mass black holes to be as large as 10 43 erg/s for observationally allowed dark photon parameter space. The emission is expected to have a significant x-ray component and to potentially be periodic, with period set by the dark photon mass. The luminosity is comparable to the brightest x-ray sources in the Universe, allowing for searches at distances of up to hundreds of Mpc with existing telescopes. Finally, we discuss observational strategies, including targeted electromagnetic follow-ups of solar-mass black hole mergers and targeted continuous gravitational wave searches of anomalous pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗

Effect of Enantiomeric Excess on the Phase Behavior of Antiferroelectric Liquid Crystals

Null transmission ellipsometry and resonant x-ray diffraction are employed to study the effect of enantiomeric excess (EE) on the phase behavior of antiferroelectric liquid crystal 10OTBBB1M7. Phase sequence, layer spacing, and pitch of the helical structures of the smectic-C*{sub {alpha}} and smectic-C* phases are studied as a function of temperature and EE. Upon reducing EE, a liquid-gas-type critical point of the smectic-C*{sub {alpha}} to smectic-C* transition is observed, as well as the disappearance of the smectic-C*{sub d4} and the smectic-C*{sub d3} phases. Results are analyzed in a mean-field model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ultrafast Melting of Carbon Induced by Intense Proton Beams

Laser-produced proton beams have been used to achieve ultrafast volumetric heating of carbon samples at solid density. The isochoric melting of carbon was probed by a scattering of x rays from a secondary laser-produced plasma. From the scattering signal, we have deduced the fraction of the material that was melted by the inhomogeneous heating. The results are compared to different theoretical approaches for the equation of state which suggests modifications from standard models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum computation of phase transition in the massive Schwinger model

Abstract As pointed out by Coleman, physical quantities in the Schwinger model depend on a parameter θ that determines the background electric field. There is a phase transition for θ = π only. We develop a momentum space formalism on a lattice and use it to perform a quantum computation of the critical point of this phase transition on the NISQ device IMB Q Lima. After error mitigation, our results give strong indication of the existence of a critical point at m / e ≃ 0.32, where m is the bare fermion mass and e is the coupling strength, in good agreement with the classical numerical result m / e ≃ 0.3335.

Physics↗

Mesoscale modeling to inform Bison models of accident tolerant fuel concepts

U 3 Si 2 and doped UO 2 have been under investigation in recent years as potential accident-tolerant fuel concepts. In this report, lower-length scale studies of these fuel concepts carried out during Fiscal Year 2020 are detailed. U 3 Si 2 is a potential accident-tolerant fuel that shows promise due to its high thermal conductivity and higher uranium density relative to UO 2 . However, its swelling and fission gas release behavior in light water reactor (LWR) conditions is relatively unknown. To provide mechanistic insight and determine parameters for engineering-scale fuel performance modeling of pellet-form U 3 Si 2 , phase-field simulations of the growth, interconnection, and venting of intergranular fission gas bubbles were performed. The fractional coverage of the grain boundary and the fraction of bubble area that is vented were calculated as a function of time. From the simulation data, the fractional grain boundary coverage at saturation, an important parameter needed in engineering-scale modeling of swelling and fission gas release, was determined. Multiple simulations were run to determine the uncertainty in the calculated value. The effect of model assumptions and input parameters that are not well known was evaluated. Simulation results are compared to related theoretical and computational work. Based on the simulation results, a value of 0.60 for the fractional grain boundary coverage at saturation is recommended for U 3 Si 2 fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Spontaneous Breaking of Scale Invariance in U(N) Chern-Simons Gauge Theories in Three Dimensions

I explore the existence of a massive phase in a conformally invariant U(N) Chern-Simons gauge theories in D = 3 with matter fields in the fundamental representation. These models have attracted recent attention as being dual, in the conformal phase, to theories of higher spin gravity on AdS 4. Using the 0t Hooft large N expansion, exact solutions are obtained for scalar current correlators in the massive phase where the conformal symmetry is spontaneously broken. A massless dilaton appears as a composite state, and its properties are discussed. Solutions exist for matters field that are either bosons or fermions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Chiral thermodynamics of dense hadronic matter

We discuss phases of hot and dense hadronic matter using chiral Lagrangians. A two-flavored parity doublet model constrained by the nuclear matter ground state predicts chiral symmetry restoration. The model thermodynamics is shown within the mean-field approximation. A field-theoretical constraint on possible phases from the anomaly matching is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mass and density effects in the QCD phase transition

We describe the critical temperature dependence on the pion mass and on the isospin chemical potential in remarkably good agreement with lattice data. Our framework incorporates explicit dependence on quark masses, isospin and baryonic chemical potentials for the case of two flavors through ingredients from well-known high-and low-energy theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Grizzly: A Multi-scale and Multi-Physics Tool to Model Aging of Nuclear Power Plant Components

Grizzly is a multi-scale multi-physics tool being developed at Idaho National Laboratory (INL) as part of the US Department of Energy's Light Water Reactor Sustainability program to provide improved safety assessments of systems, components, and structures in nuclear power plants subjected to age-related degradation. Its goal is to provide an improved scientific basis for decisions surrounding license renewal, which would permit operation of commercial nuclear power plants beyond 60 years. Grizzly is based on INL's MOOSE framework, which enables multi-physics simulations in a massively parallel computing environment. The reactor pressure vessel (RPV) has been chosen as the initial application for Grizzly. The primary concern for the RPV is that it becomes increasingly susceptible to fracture initiation at the sites of pre-existing flaws over time. This is because exposure to irradiation and high temperatures embrittles the steel. Assessing the safety of an RPV under transient conditions involves simulations at multiple scales. At the engineering scale, strongly coupled equations of heat conduction and solid mechanics are solved to simulate the global response of the RPV to accident conditions. Sub-models are used to represent regions with pre-existing flaws and the stress intensity factors are evaluated from the domain integrals. Crack propagation is then determined using a model that provides the temperature and neutron fluence dependent fracture toughness of the RPV steel. The ability of Grizzly to provide accurate engineering scale estimates of RPV safety beyond 60 years of operation strongly relies on the predictability of the embrittlement model. Multiple models for embrittlement have been proposed, and the physically based empirical model is implemented in Grizzly. However, the formation of late blooming phases greatly reduces the confidence on the extrapolated values obtained from embrittlement model, which is further exacerbated by the absence of real field data for validation. Experimentally validated lower length scale models can significantly improve the predictability of the embrittlement model beyond the calibrated range and the development of such models is another ongoing effort in Grizzly. The modeling activity in Grizzly spans across different length scales, and can be divided into: (i) Irradiation induced microstructure evolution at the atomistic and single crystal length scale using molecular dynamics, atomistic kinetic Monte-Carlo, phase-field and rate theory; (ii) Effect of the irradiation induced defects on the flow stress and fracture processes at the polycrystalline scale using crystal plasticity. Subsequently, the prediction from these models can be used to obtain irradiation induced ductile-to-brittle transition temperature (DBTT) shifts from specimen-scale simulations; (iii) Global thermo-mechanical response of the RPV under normal and accident scenarios, and probabilistic risk analyses at the engineering scale using the DBTT shifts obtained from the lower length scale models. The item (i) has been discussed in a previous report. In this extended abstract an overview of items (ii) and (iii) is provided. The extended abstract summarizes the current status of the modeling effort under Grizzly. At the engineering scale, reduced order model development is underway using the statistical tool RAVEN. This will allow probabilistic analysis considering a wide range of flaw topologies in a computationally efficient manner. Coupling the irradiation induced defect growth predictions to crystal plasticity is also in progress. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗