Phase stability in the three-dimensional open-source code for the chiral mean-field model
Not Available
Cruz-Camacho, Nikolas (ORCID:0009000478700039)↗
SEARCH · Search NASA
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.
Interactions Enhance Ramp Reversal Memory in Locally Phase Separated Materials
The ramp-reversal memory (RRM) effect in metal–insulator transition metal oxides (TMOs), a non-volatile resistance change induced by repeated temperature cycling, has attracted considerable interest in neuromorphic computing and non-volatile memory devices. Our previous defect motion model successfully explained RRM in vanadium dioxide (VO 2 ), capturing observed critical temperature shifts and memory accumulation throughout the sample. However, this approach lacked interactions between metallic and insulating domains. Here, we extend our model by combining a correlated Random Field Ising Model with defect diffusion-segregation, enabling accurate hysteresis modeling while predicting the relationship between RRM and domain interactions. Our simulations demonstrate that the maximum RRM occurs when the turnaround temperature approaches the inflection point. This peak in RRM vs. turnaround temperature is consistent with prior transport measurements, as well as our own optical measurements reported here. Significantly, we find that increasing nearest-neighbor interactions enhances the maximum memory effect, thus providing a clear mechanism for optimizing RRM performance. Since our model employs minimal assumptions, we predict that RRM should be a widespread phenomenon in materials exhibiting patterned phase coexistence of electronic domains. This work not only advances fundamental understanding of memory behavior in TMOs but also establishes a much-needed theoretical framework for optimizing device applications.
36 MATERIALS SCIENCE↗
Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids
Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.
Aluminum alloy↗
Highly tunable magnetocrystalline anisotropy energy in Fe3+-doped BaTiO3
Magnetic dopants in ferroelectric oxide host materials provide a platform for electric field control of isolated spins, facilitated by tuning of the magnetocrystalline anisotropy energy (MCAE). We present first-principles calculations of the MCAE experienced by isolated Fe3+ dopants in the tetragonal, orthorhombic, and rhombohedral phases of the prototypical ferroelectric BaTiO3. We identify an order-of-magnitude decrease in the MCAE in the rhombohedral phase relative to the tetragonal and orthorhombic phases. We explain this dramatic decrease, as well as the formation of a spin-easy plane in the tetragonal phase and spin-easy axes in the orthorhombic and rhombohedral phases, using crystal field theory arguments. Building a superposition model from crystal field theory, we show how a set of simple criteria based on crystalline environment can be used to estimate the MCAE. We suggest this as a route to rapidly screen candidate ferroelectric hosts and magnetic dopants that possess phases with spin-easy axes and maximal MCAE tunability.
Barker, Bradford A↗
Navier-Stokes simulations of unsteady transonic flow phenomena
Numerical simulations of two classes of unsteady flows are obtained via the Navier-Stokes equations: a blast-wave/target interaction problem class and a transonic cavity flow problem class. The method developed for the viscous blast-wave/target interaction problem assumes a laminar, perfect gas implemented in a structured finite-volume framework. The approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the full viscous terms are computed by central differencing. Comparisons of unsteady numerical, analytical, and experimental results are made in two- and three-dimensions for Couette flows, a starting shock-tunnel, and a shock-tube blockage study. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing of the predominantly inviscid flows. Viscous effects were increasingly significant at large post-interaction times. While the blast-wave/target interaction problem benefits from high-resolution methods applied to the Euler terms, the transonic cavity flow problem requires the use of an efficient scheme implemented in a geometrically flexible overset mesh environment. Hence, the Reynolds averaged Navier-Stokes equations implemented in a diagonal form are applied to the cavity flow class of problems. Comparisons between numerical and experimental results are made in two-dimensions for free shear layers and both rectangular and quieted cavities, and in three-dimensions for Stratospheric Observatory For Infrared Astronomy (SOFIA) geometries. The acoustic behavior of the rectangular and three-dimensional cavity flows compare well with experiment in terms of frequency, magnitude, and quieting trends. However, there is a more rapid decrease in computed acoustic energy with frequency than observed experimentally owing to numerical dissipation. In addition, optical phase distortion due to the time-varying density field is modelled using geometrical constructs. The computed optical distortion trends compare with the experimentally inferred result, but underpredicts the fluctuating phase difference magnitude.
Atwood, C. A.↗
Historical Total Ozone Radiative Forcing Derived from CMIP6 Simulations
Radiative forcing (RF) time series for total ozone from 1850 up to present-day are calculated based on historical simulations of ozone from 10 climate models contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). In addition, RF is calculated for ozone fields prepared as input for CMIP6 models without chemistry schemes and from a chemical transport model simulation. A radiative kernel for ozone is constructed and used to derive the RF. The ozone RF in 2010 (200530 -2014) relative to 1850 is 0.35 W m-2 [0.08 to 0.61] (5-95% uncertainty range) based on models with both tropospheric and stratospheric chemistry. One of these models has a negative present-day total ozone RF. Excluding this model, the present-day ozone RF increases to 0.39 W m-2 [0.27 to 0.51] (533 - 95% uncertainty range). The rest of the models have RF close to or stronger than the RF time series assessed by the Intergovernmental Panel on Climate Change in the fifth assessment report with the primary driver likely being the new precursor emissions used in CMIP6. The rapid adjustments beyond stratospheric temperature are estimated to be weak and thus the RF is a good measure of effective radiative forcing.
ozone↗
Chemistry of a polluted cloudy boundary layer
A one-dimensional photochemical model for cloud-topped boundary layers has been developed to include descriptions of gas- and aqueous-phase chemistry and the radiation field in and below the cloud. The model is applied to the accumulation of pollutants during a wintertime episode with low stratus over Bakersfield, CA. The mechanisms of sulfate production and the balance between the concentrations of acids and bases are examined. It is shown that most of the sulfate production may be explained by the Fe(III)-catalyzed autoxidation of S(IV). Another source of sulfate is the oxidation of SO2 by OH in both the gas and the aqueous phase. It is shown that the sulfate production in the model is controlled by the availability of NH3. It is suggested that this explains the balance observed between total concentration of acids and bases.
Jacob, Daniel J.↗
Robust Spin-Moiré-Superlattice-Driven Gap Opening in EuAg 4 Sb 2 under in-Plane Magnetic Field
Moiré superlattices introduce new length and energy scales, enabling discoveries, such as unconventional superconductivity, in 2D systems. This concept has recently been extended to bulk materials with multiple-q spin textures, opening exciting opportunities for spin moiré physics. A notable example is EuAg 4 S b2 , where a spin moiré superlattice (SMS), manifested as a double-q spin modulation, induces a superzone gap opening. In this work, we investigate the tunability of this gap by an in-plane magnetic field in EuAg 4 Sb 2 by using neutron scattering, magnetization, and transport measurements. We reveal a cycloidal ground state and multiple spin-reoriented phases induced by the in-plane field, highlighting the critical role of in-plane components in driving magnetic transitions. Moreover, we demonstrate that a robust gap opening persists in the double-q phase, regardless of in-plane field orientation. Model calculations attribute this robustness to the stability of the SMS under tilted fields. Here, our results establish EuAg 4 Sb 2 as a tunable platform for exploring the spin-texture-driven superzone gap opening in electronic states.
36 MATERIALS SCIENCE↗
Modeling Macro- and Micro-Scale Turbulent Mixing and Chemistry in Engine Exhaust Plumes
Simulation of turbulent mixing and chemical processes in the near-field plume and plume-vortex regimes has been successfully carried out recently using a reduced gas phase kinetics mechanism which substantially decreased the computational cost. A detailed mechanism including gas phase HOx, NOx, and SOx chemistry between the aircraft exhaust and the ambient air in near-field aircraft plumes is compiled. A reduced mechanism capturing the major chemical pathways is developed. Predictions by the reduced mechanism are found to be in good agreement with those by the detailed mechanism. With the reduced chemistry, the computer CPU time is saved by a factor of more than 3.5 for the near-field plume modeling. Distributions of major chemical species are obtained and analyzed. The computed sensitivities of major species with respect to reaction step are deduced for identification of the dominant gas phase kinetic reaction pathways in the jet plume. Both the near field plume and the plume-vortex regimes were investigated using advanced mixing models. In the near field, a stand-alone mixing model was used to investigate the impact of turbulent mixing on the micro- and macro-scale mixing processes using a reduced reaction kinetics model. The plume-vortex regime was simulated using a large-eddy simulation model. Vortex plume behind Boeing 737 and 747 aircraft was simulated along with relevant kinetics. Many features of the computed flow field show reasonable agreement with data. The entrainment of the engine plumes into the wing tip vortices and also the partial detrainment of the plume were numerically captured. The impact of fluid mechanics on the chemical processes was also studied. Results show that there are significant differences between spatial and temporal simulations especially in the predicted SO3 concentrations. This has important implications for the prediction of sulfuric acid aerosols in the wake and may partly explain the discrepancy between past numerical studies (that employed parabolic or temporal approximations) and the measured data. Finally to address the major uncertainty in the near-field plume modeling related to the plume processing of sulfur compounds and advanced model was developed to evaluate its impact on the chemical processes in the near wake. A comprehensive aerosol model is developed and it is coupled with chemical kinetics and the axisymmetric turbulent jet flow models. The integrated model is used to simulate microphysical processes in the near-field jet plume, including sulfuric acid and water binary homogeneous nucleation, coagulation, non-equilibrium heteromolecular condensation, and sulfur-induced soot activation. The formation and evolution of aerosols are computed and analyzed. The computed results show that a large number of ultra-fine (0.3--0.6 nm in radius) volatile HSO4 - HO embryos are generated in the near-field plume. These embryos further grow in size by self coagulation and condensation. Soot particles can be activated by both heterogeneous nucleation and scavenging of H2SO4-H2O aerosols. These activated soot particles can serve as water condensation nuclei for contrail formation. Conditions under which ice contrails can form behind aircrafts are studied. The sensitivities of the threshold temperature for contrail formation with respect to aircraft propulsion efficiency, relative humidity, and ambient pressure are evaluated. The computed aerosol properties for different extent of fuel sulfur conversion to S(VI) (SO3 and H2SO4) in engine are examined and the results are found to be sensitive to this conversion fraction.
Menon, Suresh↗
Topological excitonic insulator with tunable momentum order
Correlated topological materials often maintain a delicate balance among physical symmetries. Many topological orders are symmetry protected, whereas most correlated phenomena arise from spontaneous symmetry breaking. Cases where symmetry breaking induces a non-trivial topological phase are rare. Here we demonstrate the presence of two such phases in Ta 2 Pd 3 Te 5 , where Coulomb interactions form excitons that condense below 100 K, one with zero and the other with finite momentum. We observed a full spectral bulk gap, which stems from exciton condensation. This topological excitonic insulator state spontaneously breaks mirror symmetries but involves a weak structural coupling. Scanning tunnelling microscopy shows gapless boundary modes in the bulk insulating phase. Their magnetic field response, together with theoretical modelling, indicates a topological origin. These observations establish Ta 2 Pd 3 Te 5 as a topological excitonic insulator in a three-dimensional crystal. Finally, our results manifest a unique sequence of topological exciton condensations in a bulk crystal, offering exciting opportunities to study critical behaviour and excitations.
36 MATERIALS SCIENCE↗
Gravity field determination for Mars Observer
Mars Observer will be the first near circular, low altitude and short-periodic orbiter of Mars. From previous Mars orbiters, such as Mariner 9 and Vikings 1 and 2, a large gravitational oblateness and in general a gravity field ten times stronger than the earth's field have been determined. Because of these, the gravity field will dominate the evolution of the spacecraft's orbital motion and is the principal error source for reconstructing the spacecraft's motion. Thus a gravity calibration (GC) period has been established prior to the mapping phase in order to refine the gravitational field model. This paper describes the plan being made by Navigation to implement the GC strategy. By analyzing simulated Doppler data acquired during the seven day GC period, we give (1) the expected improvement in the gravity field model and (2) the resulting improvement in reconstructed spacecraft orbital motion. In addition, we have developed a strategy for converging the GC gravity model.
Esposito, Pasquale B.↗
Current Status on the use of Parallel Computing in Turbulent Reacting Flow Computations Involving Sprays, Monte Carlo PDF and Unstructured Grids
The state of the art in multidimensional combustor modeling as evidenced by the level of sophistication employed in terms of modeling and numerical accuracy considerations, is also dictated by the available computer memory and turnaround times afforded by present-day computers. With the aim of advancing the current multi-dimensional computational tools used in the design of advanced technology combustors, a solution procedure is developed that combines the novelty of the coupled CFD/spray/scalar Monte Carlo PDF (Probability Density Function) computations on unstructured grids with the ability to run on parallel architectures. In this approach, the mean gas-phase velocity and turbulence fields are determined from a standard turbulence model, the joint composition of species and enthalpy from the solution of a modeled PDF transport equation, and a Lagrangian-based dilute spray model is used for the liquid-phase representation. The gas-turbine combustor flows are often characterized by a complex interaction between various physical processes associated with the interaction between the liquid and gas phases, droplet vaporization, turbulent mixing, heat release associated with chemical kinetics, radiative heat transfer associated with highly absorbing and radiating species, among others. The rate controlling processes often interact with each other at various disparate time 1 and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and liquid phase evaporation in many practical combustion devices.
Raju, M. S.↗
Superalloy composition modeling
Development of a predictive method for determination of the gamma/gamma prime phase fields, i.e., gamma prime volume fraction as a function of the multicomponent composition, is described. The cluster variation method used for binary alloys in which the precipitated phase is coherent with the matrix phase is extended for application to the multicomponent coherent gamma/gamma prime nickel-based superalloys. It is shown that the cluster variation method can accurately describe the equilibrium (incoherent) gamma/gamma prime phase fields in the binary Ni-Al phase diagram. The gamma/gamma prime phase field for the Ni-Cr-Al ternary phase diagram is computed as a function of temperature. A reasonable fit results between the calculated and the experimental diagrams. The modeling of the six-component Ni-Cr-Al-Co-Mo-Ti base superalloy and the effect of Ni substitution of Co are discussed.
Barefoot, J.↗
Hybird state of the tail mangetic configuration during steady convection events
Previous observations have shown that during periods of steady magnetospheric convection (SMC) a large amount of magnetic flux crosses the plasma sheet (corresponding to approximately 10 deg wide auroral oval at the nightside) and that the magnetic configuration in the midtail is relaxed (the curent sheet is thick and contains enhanced B(sub Z). These signatures are typical for the substorm recovery phase. Using near-geostationary magnetic field data, magnetic field modeling and a noval diagostic technique (isotropic boundary algorithm), we show that in the near-Earth tail the magnetic confirguration is very stretched during the SMC events. This stretching is caused by an intense, thin westward current. Because of the srongly depressed B(sub Z), there is a large radial gradient in the near-tail magetic field. These signatures have been peviously associated only with the substorm growth phase. Our results indicate that during the SMC periods the magnetic configuration is very peculiar, with co-existing thin near-Earth current sheet and thick midtail plasma sheet. The deep local minimum of the equatorial B(sub Z) that devleops at R approximately 12 R(sub E) is consistent with steady, adiabatic, Earthward convection in the midtail. These results impose contraints on the existing substorm theories, and call for an explanation of how such a stressed configuration can persist for such a long time without tail current disruptions that occur at the end of a substorm growth phase.
Sergeev, V. A.↗
Generation of Large-Volume High-Pressure Plasma by Spatio-Temporal Control of Space Charge
Due to the fundamental limitations of scaling, very little progress has been made towards achieving a large-volume dense non-equilibrium plasma near atmospheric pressures. Commercially available state-of-the art glow-like atmospheric plasma devices for industrial applications are either narrow tubes or wide slits with narrow openings. Often multiple sources are put together in various configurations to process larger surfaces. The traditional methods of exciting electrodes create spatially fixed electric fields. As a result, the space charge at atmospheric pressure tends to be spatially confined resulting in non-uniformity which eventually leads to instability as the discharge is scaled. Theoretical work done under this project showed that it is possible to generate a spatially rotating electric field by exciting a set of electrodes with phase staggered sinusoidal waveforms. The modeling and simulations were done using plasma fluid models. It was shown that such a field can produce a uniform plasma. At the conclusion of the project, experimental proof of concept with an eight-electrode system in various gases (Air, Helium and Argon) was demonstrated. Power measurements and spectral investigation show that the concept can be used to generate a uniform stable plasma. This plasma source has the potential of opening new applications of nonthermal plasma including combustion of carbon-free fuel. The current limitation of the proposed method in scaling to higher volume and pressure is the need for multiple high voltage amplifiers.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Gallery of soft modes: Theory and experiment at a ferromagnetic quantum phase transition
We examine the low-energy excitations in the vicinity of the quantum critical point in LiHoF4, a physical realization of the Transverse Field Ising Model, focusing on the long-range fluctuations which soften to zero energy at the ferromagnetic quantum phase transition. Microwave spectroscopy in tunable loop-gap resonator structures identifies and characterizes the electronuclear soft mode and higher-energy electronuclear states as a function of frequency and magnetic fields applied transverse and parallel to the Ising axis. These are understood in the context of a theoretical model of a soft electronuclear mode that interacts with soft photons as well as soft phonons. We identify competing infrared divergences at the quantum critical point, coming from the photons and the electronuclear soft mode. It is an incomplete cancellation of these divergences that leads to the muted but distinct signatures observed in the experiments. The application of a longitudinal magnetic field gaps the soft mode. As a result, measurements well away from the quantum critical point reveal a set of "Walker'' modes associated with ferromagnetic domain dynamics.
Dipolar interaction↗
A study of backscattering and emission from closely packed inhomogeneous media
The effects of close spacing between small scattering spheres were examined by keeping the distance-dependent terms in the expressions for the transverse scattered fields. The phase matrix was then derived from these fields and was used in the radiative transfer formulation to model scattering and emission from a densely populated, inhomogeneous layer. Computed results were compared with those obtained when the phase matrix was specialized to the far-field condition. It was found that the use of the far-zone condition tended to underestimate both the level of the copolarized backscattering and the cross-polarized backscattering. In emission computations, the use of the far-zone condition overestimated the level of the brightness temperature. These effects decreased with a decrease in the volume fraction or an increase in the exploring frequency, as expected. An improvement on the snow parameter (density and crystal size) estimation was shown to be possible when this new phase matrix was used.
Fung, A. K.↗
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↗