Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field modeling”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Detecting Quantum Critical Points of Correlated Systems by Quantum Convolutional Neural Network Using Data from Variational Quantum Eigensolver

Machine learning has been applied to a wide variety of models, from classical statistical mechanics to quantum strongly correlated systems, for classifying phase transitions. The recently proposed quantum convolutional neural network (QCNN) provides a new framework for using quantum circuits instead of classical neural networks as the backbone of classification methods. We present the results from training the QCNN by the wavefunctions of the variational quantum eigensolver for the one-dimensional transverse field Ising model (TFIM). We demonstrate that the QCNN identifies wavefunctions corresponding to the paramagnetic and ferromagnetic phases of the TFIM with reasonable accuracy. The QCNN can be trained to predict the corresponding ‘phase’ of wavefunctions around the putative quantum critical point even though it is trained by wavefunctions far away. The paper provides a basis for exploiting the QCNN to identify the quantum critical point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Correction to “Frank–Kasper Phases of Diblock Copolymer Melts Studied with the DPD Model: SCF Results”

In a recent study, (1) we performed the self-consistent field (SCF) calculations of a dissipative particle dynamics (DPD) model of conformationally asymmetric diblock copolymer A-B melts to study the stability of various Frank–Kasper (FK) phases formed by such A-B melts. While all the equations in ref (1) are correct, we did not calculate the dimensionless (mean-field) Helmholtz free energies per chain βf c = βu c,χ + βu c,κ – s c /k B and the dimensionless stresses dβf c /dθ of ordered phases accurately. For example, the dimensionless (nonbonded) internal energy per chain due to the system compressibility βu c,κ = (N/2κ)∑ q βû 0 (|q|)(ϕ̂ A (q) + ϕ̂ B (q) – δ q,0 )(ϕ̂ A (−q) + ϕ̂ B (−q) – δ– q,0 ) should be calculated in the reciprocal space as given (note that βû 0 (|q|) is obtained analytically, and that since ϕ A (r) is real 𝜙̂ A (−𝐪)= $\overline{𝜙̂ A}$(𝐪), i.e., the complex conjugate of ϕ̂ A (q)); this corresponds to numerically evaluating βu c,κ = (N/2κ)∫dr(ϕ A (r) + ϕ B (r) – 1)ψ(r) with ψ(r) ≡ (∫dr′βu 0 (|r–r′|)(ϕ A (r′) + ϕ B (r′) – 1)/V in the real space via the (composite) trapezoidal rule. While the trapezoidal rule is less accurate than the Romberg integration (2) in general, the opposite occurs when the integration domain is periodic; this was pointed out by Matsen in ref (3) and explained in detail in refs (4and5). With our ignorance of this exceptional convergence property of the trapezoidal rule for periodic functions, we calculated all convolutions having the form of ψ(r) in the reciprocal space via the fast Fourier transforms (6) but those having the form of βu c,κ in the real space via Romberg integration; the latter leads to inaccurate results and unfortunately incorrect phase diagrams reported in ref (1), which are corrected here.

He, Juntong↗
DOE OSTI2011

Holography, heavy-quark free energy, and the QCD phase diagram

We use gauge/string duality to investigate the free energy of two static color sources (a heavy-quark-antiquark pair) in a Yang-Mills theory in strongly interacting matter, varying temperature and chemical potential. The dual space geometry is anti-de Sitter with a charged black hole to describe finite temperature and density in the boundary theory, and we also include a background warp factor to generate confinement. The resulting deconfinement line in the {mu}-T plane is similar to the one obtained by lattice and effective models of QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1994

A three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawn and dusk sides of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H+ fluxes at tens of keV, which are always over-estimated. A newly-invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm-time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, j(sub o)(1+Ay(exp n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗
NASA NTRS1995

Three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawnside and duskside of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H(+) fluxes at tens of keV, which are always overestimated. A newly invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, J(sub o)(1 + Ay(sup n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗
NASA NTRS1996

Photon Number-Phase Uncertainty Relation in the Evolution of the Field in a Kerr-Like Medium

A model of a single-mode field, initially prepared in a coherent state, coupled to a two-level atom surrounded by a nonlinear Kerr-like medium contained inside a very good quality cavity is considered. We derive the photon number-phase uncertainty relation in the evolution of the field for a weak and strong nonlinear coupling respectively, within the Hermitian phase operator formalism of Pegg and Barnett, and discuss the effects of nonlinear coupling of the Kerr-like medium on photon number-phase uncertainty relation of the field.

Fan, An-Fu↗
DOE OSTI2016

Changes in geophysical properties caused by fluid injection into porous rocks: analytical models: Geophysical changes in porous rocks

Analytical models are provided that describe how the elastic compliance, electrical conductivity, and fluid-flow permeability of rocks depend on stress and fluid pressure. In order to explain published laboratory data on how seismic velocities and electrical conductivity vary in sandstones and granites, the models require a population of cracks to be present in a possibly porous host phase. The central objective is to obtain a consistent mean-field analytical model that shows how each modeled rock property depends on the nature of the crack population. We describe the crack populations by a crack density, a probability distribution for the crack apertures and radii, and the averaged orientation of the cracks. The possibly anisotropic nature of the elasticity, conductivity, and permeability tensors is allowed for; however, only the isotropic limit is used when comparing to laboratory data. For the transport properties of conductivity and permeability, the percolation effect of the crack population linking up to form a connected path across a sample is modeled. But, this effect is important only in crystalline rock where the host phase has very small conductivity and permeability. In general, the importance of the crack population to the transport properties increases as the host phase becomes less conductive and less permeable.

58 GEOSCIENCES↗
DOE OSTI2025

Ferroelectric phase transition in group-IV monochalcogenides from an equivariant machine learned force field

Group-IV monochalcogenides are a class of layered ferroelectric semiconductors that have demonstrated spontaneous intrinsic polarization above room temperature. Here, in this study, we use the multi-atomic cluster expansion (MACE) machine learning architecture to train and test a force field capable of modeling the structural properties and second-order ferroelectric-to-paraelectric phase transition in a Group-IV monochalcogenide, GeSe. The model captures the double-well potential energy surface associated with the onset of macroscopic polarization in bulk GeSe within 12.5 meV/atom, as well as near-equilibrium properties like the phonon dispersion. The development of this quantitatively accurate force field enables long-time molecular dynamics simulations, which predict the critical temperature of the ferroelectric-to-paraelectric phase transition in bulk GeSe to be T c = 600 K. This study demonstrates the capabilities of equivariant force-fields to accurately describe phenomena associated with structural symmetry breaking.

ferroelectricity↗
DOE OSTI2016

Observable signatures of a classical transition

Eternal inflation arising from a potential landscape predicts that our universe is one realization of many possible cosmological histories. One way to access different cosmological histories is via the nucleation of bubble universes from a metastable false vacuum. Another way to sample different cosmological histories is via classical transitions, the creation of pocket universes through the collision between bubbles. Using relativistic numerical simulations, we examine the possibility of observationally determining if our observable universe resulted from a classical transition. We find that classical transitions produce spatially infinite, approximately open Friedman-Robertson-Walker universes. The leading set of observables in the aftermath of a classical transition are negative spatial curvature and a contribution to the Cosmic Microwave Background temperature quadrupole. The level of curvature and magnitude of the quadrupole are dependent on the position of the observer, and we determine the possible range of observables for two classes of single-scalar field models. For the first class, where the inflationary phase has a lower energy than the vacuum preceding the classical transition, the magnitude of the observed quadrupole generally falls to zero with distance from the collision while the spatial curvature grows to a constant. For the second class, where the inflationary phase has a higher energy than the vacuum preceding the classical transition, the magnitude of the observed quadrupole generically falls to zero with distance from the collision while the spatial curvature grows without bound. We find that the magnitude of the quadrupole and curvature grow with increasing centre of mass energy of the collision, and explore variations of the parameters in the scalar field lagrangian.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Hundness versus Mottness in a three-band Hubbard–Hund model: On the origin of strong correlations in Hund metals☆

Hund metals are multi-orbital systems with moderate Coulomb interaction, U, among charges and sizeable Hund’s rule coupling, J(<U), that aligns the spins in different orbitals. They show strong correlation effects, like very low Fermi-liquid coherence scales and intriguing incoherent transport regimes, resulting in bad metallic behavior. But to what extent are these strong correlations governed by Mottness, i.e. the blocking of charge fluctuations close to a Mott insulator transition (MIT) induced by U, or by Hundness, a new route towards strong correlations induced by J? To answer this question, we study the full phase diagram of a degenerate three-band Hubbard–Hund model on a Bethe lattice at zero temperature using single-site dynamical mean-field theory and the numerical renormalization group as efficient real-frequency multi-band impurity solver. Hund metal behavior occurs in this minimal model for a filling close to n{sub d}=2, moderate U and sizeable J, the “Hund-metal regime”. In particular, strong correlations manifest themselves there by an unusually low quasiparticle weight. Generalizing previous results on this model, we show that “spin–orbital separation” (SOS) is a generic Hund’s-coupling-induced feature in the whole metallic regime of the phase diagram for 1<n{sub d}<3 and sizeable J. There orbital screening always occurs at much higher energies than spin screening below which Fermi-liquid behavior sets in. The low quasiparticle weight can then be directly explained in terms of the Hund’s-coupling-reduced Fermi-liquid scale. We carefully analyze the effect of J (Hundness), and the effect of the MIT at n{sub d}=2 and n{sub d}=3 (Mottness) on the energy scales and the nature of SOS. In the Hund-metal regime, far from any MIT, Hundness is shown to be the key player to induce strong correlations by constraining the spin rather than the charge dynamics. There, physical properties are governed by a broad incoherent energy regime of SOS where intriguing Hund metal physics occurs: large, almost unscreened spins are coupled to screened orbital degrees of freedom. With increasing proximity to an MIT correlations are further enhanced and the Fermi-liquid scale is further reduced. However, in the Hund-metal regime, this effect of Mottness is minor. In contrast, very close to the MIT at n{sub d}=2, the incoherent spin–orbital separation regime is strongly downscaled and becomes negligibly small, whereas Mottness – the localization of charges – becomes dominant in inducing strong correlations. Close to the MIT at n{sub d}=3, the SOS regime widens up because the orbital degrees of freedom get blocked by the formation of an S=(3/2) impurity spin, but its nature changes: the orbital and spin dynamics get decoupled. Our results confirm Hundness as a distinct mechanism towards strong correlations in the normal state of Hund metals, leading to various interesting implications for the nature of electronic transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

Investigation of the Alfven Wave Turbulence on the Dynamics of SEPs

Solar energetic particles (SEPs) are high-energy particles originating from the Sun that pose potential risks to space missions, especially those outside Earth's protective magnetosphere. Understanding their behavior within the heliosphere and Earth's magnetosphere is vital for the safety and functionality of space exploration. The study focuses on how the pitch angle scattering, the angle change between a particle's velocity and the magnetic field due to interactions with magnetic irregularities, impacts the decay phase of SEP events. By solving the focused transport equation, the study models the transport of SEPs along magnetic field lines, extending these lines up to 5 AU (Astronomical Units) to understand their behavior even beyond the typical proximity to the Sun. Incorporating the effects of pitch angle scattering at such vast distances is crucial for understanding the decay dynamics of a SEP event. This modeling of SEPs is integrated with simulations of other vital space phenomena: the solar wind, interplanetary magnetic field, and Alfven wave turbulence. This presentation covers the modeling approach used in this research and the role of pitch angle scattering at various heliocentric distances on the dynamics of the SEP events' decay phase.

solar energetic particles↗
DOE OSTI2011

Disordered spinor Bose-Hubbard model

We study the zero-temperature phase diagram of the disordered spin-1 Bose-Hubbard model in a two-dimensional square lattice. To this aim, we use a mean-field Gutzwiller ansatz and a probabilistic mean-field perturbation theory. The spin interaction induces two different regimes, corresponding to a ferromagnetic and antiferromagnetic order. In the ferromagnetic case, the introduction of disorder reproduces analogous features of the disordered scalar Bose-Hubbard model, consisting in the formation of a Bose glass phase between Mott insulator lobes. In the antiferromagnetic regime, the phase diagram differs more from the scalar case. Disorder in the chemical potential can lead to the disappearance of Mott insulator lobes with an odd-integer filling factor and, for sufficiently strong spin coupling, to Bose glass of singlets between even-filling Mott insulator lobes. Disorder in the spinor coupling parameter results in the appearance of a Bose glass phase only between the n and the n+1 lobes for n odd. Disorder in the scalar Hubbard interaction inhibits Mott insulator regions for occupation larger than a critical value.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm

The origin of the Finite Grid Instability (FGI) is studied by resolving the dynamics in the 1D electrostatic Particle-In-Cell (PIC) model in the spectral domain at the single particle level and at the collective motion level. The spectral fidelity of the PIC model is contrasted with the underlying physical system or the gridless model. The systematic spectral phase and amplitude errors from the charge deposition and field interpolation are quantified for common particle shapes used in the PIC models. Lastly, it is shown through such analysis and in simulations that the lack of spectral fidelity relative to the physical system due to the existence of aliased spatial modes is the major cause of the FGI in the PIC model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Extension of Ostwald Ripening Theory

The objective is to develop models based on the mean field approximation of Ostwald ripening to describe the growth of second phase droplets or crystallites. The models will include time variations in nucleation rate, control of saturation through addition of solute, precipitating agents, changes in temperature, and various surface kinetic effects. Numerical integration schemes have been developed and tested against the asymptotic solution of Liftshitz, Slyozov and Wagner (LSW). A second attractor (in addition to the LSW distribution) has been found and, contrary to the LSW theory, the final distribution is dependent on the initial distribution. A series of microgravity experiments is being planned to test this and other results from this work.

Baird, J.↗
DOE OSTI2015

β-decay of magic nuclei: Beyond mean-field description

Nuclear β-decay plays an important role not only in nuclear physics but also in astrophysics. The widely used self-consistent Random Phase Approximation (RPA) models tend to overestimate the half-lives of magic nuclei. To overcome this problem, we go beyond the mean-field description and include the effects of particle-vibration coupling (PVC) on top of the RPA model. The β-decay half-lives of {sup 34}Si, {sup 68}Ni, {sup 78}Ni, and {sup 132}Sn are studied within this approach in the case of the Skyrme interaction SkM*. It is found that the low-lying Gamow-Teller (GT) strength is shifted downwards with the inclusion of the PVC effect, and as a consequence, the half-lives are reduced due to the increase of the phase space available for β-decay, which leads to a good agreement between theoretical and experimental lifetimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Responses of estuarine circulation and salinity to the loss of intertidal flats – A modeling study

Intertidal flats in estuaries are coastal wetlands that provide critical marine habitats to support wide ranges of marine species. Over the last century many estuarine systems have experienced significant loss of intertidal flats due to anthropogenic impacts. This paper presents a modeling study conducted to investigate the responses of estuarine hydrodynamics to the loss of intertidal flats caused by anthropogenic actions in Whidbey Basin of Puget Sound on the northwest coast of North America. Changes in salinity intrusion limits in the estuaries, salinity stratification, and circulation in intertidal flats and estuaries were evaluated by comparing model results under the existing baseline condition and the no-flat condition. Model results showed that loss of intertidal flats results in an increase in salinity intrusion, stronger mixing, and a phase shift in salinity and velocity fields in the bay front areas. Model results also indicated that loss of intertidal flats enhances two-layer circulation, especially the bottom water intrusion. Loss of intertidal flats increases the mean salinity but reduces the salinity range in the subtidal flats over a tidal cycle because of increased mixing. Salinity intrusion limits extend upstream in all three major rivers discharging into Whidbey Basin when no intertidal flats are present. Changes in salinity intrusion and estuarine circulation patterns due to loss of intertidal flats affect the nearshore habitat and water quality in estuaries and potentially increase risk of coastal hazards, such as storm surge and coastal flooding. Furthermore, model results suggested the importance of including intertidal flats and the wetting-and-drying process in hydrodynamic simulations when intertidal flats are present in the model domain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

Development of integrated core disruptive accident analysis code for FBR - ASTERIA-FBR

The evaluation of consequence at the severe accident is the most important as a safety licensing issue for the reactor core of liquid metal cooled fast breeder reactor (LMFBR), since the LMFBR core is not in an optimum condition from the viewpoint of reactivity. This characteristics might induce a super-prompt criticality due to the core geometry change during the core disruptive accident (CDA). The previous CDA analysis codes have been modeled in plural phases dependent on the mechanism driving a super-prompt criticality. Then, the following event is calculated by connecting different codes. This scheme, however, should introduce uncertainty and/or arbitrary to calculation results. To resolve the issues and obtain the consistent calculation results without arbitrary, JNES is developing the ASTERIA-FBR code for the purpose of providing the cross-check analysis code, which is another required scheme to confirm the validity of the evaluation results prepared by applicants, in the safety licensing procedure of the planned high performance core of Monju. ASTERIA-FBR consists of the three major calculation modules, CONCORD, dynamic-GMVP, and FEMAXI-FBR. CONCORD is a three-dimensional thermal-hydraulics calculation module with multi-phase, multi-component, and multi-velocity field model. Dynamic-GMVP is a space-time neutronics calculation module. FEMAXI-FBR calculates the fuel pellet deformation behavior and fuel pin failure behavior. This paper describes the needs of ASTERIA-FBR development, major module outlines, and the model validation status. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRSDate not supplied

The Effect of Stochastic Acceleration on the Dynamics of Solar Energetic Particles in the Heliosphere

Solar Energetic Particles (SEPs) are crucial in space weather phenomena, impacting space-based technologies and human activities in space. Understanding the dynamics of SEPs, including their acceleration and transport mechanisms, is essential for improving predictive models and mitigating adverse effects. This study investigates the impact of stochastic acceleration (SA) on the dynamics of SEPs within the heliosphere. By integrating stochastic acceleration into a SEP transport model, we assess its impact relative to diffusive shock acceleration (DSA). The investigation utilizes the Adaptive Mesh Particle Simulator (AMPS) within the Space Weather Modeling Framework (SWMF) to simulate SEP dynamics, incorporating solar wind parameters and turbulence from the Alfven Wave Solar Model (AWSoM). To do this, we solve the focused transport equation to model the propagation of SEPs along magnetic field lines, extending these lines to 5 AU to incorporate the potential effects of pitch angle scattering beyond 1 AU. This aspect could notably impact the decay phase of a SEP event. Modeling the transport of SEPs is coupled with the simulation of solar wind dynamics, the interplanetary magnetic field, and the characteristics of Alfven wave turbulence. This presentation covers the modeling approach used in this research, the characterization of the impact of the stochastic acceleration, and the role of pitch angle scattering at various heliocentric distances on the dynamics of the SEP events' decay phase.

SEP↗