Search NASA⌕ Search

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.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Primordial black holes and magnetic fields in conformal neutrino mass models

Sufficiently strong and long-lasting first-order phase transitions can produce primordial black holes (PBHs) that contribute substantially to the dark matter abundance of the Universe, and can produce large-scale primordial magnetic fields. Here, we study these mechanisms in a generic class of conformal U(1)' models that also explain active neutrino oscillation data via the type-I seesaw mechanism. We find that phase transitions that occur at seesaw scales between 10 4 GeV and 10 11 GeV produce gravitational wave signals (from the dynamics of the phase transition and from the decay of cosmic string loops) at LISA/ET that can be correlated with microlensing signals of PBHs at the Roman Space Telescope, while scales near 10 11 GeV can be correlated with Hawking evaporation signals at future gamma-ray telescopes. LISA can probe the entire range of PBH masses between 1 × 10 -16 M ⊙ and 8 × 10 -11 M ⊙ if PBHs fully account for the dark matter abundance. For Z' masses between 40 TeV and 10 4 TeV, and 10 TeV right-handed neutrinos, helical magnetic fields can be produced with magnitudes ≳ 0.5 pG and coherence lengths ≳ 0.008 Mpc, above current blazar lower bounds.

cosmological phase transitions↗
DOE OSTI2025

Strong Kitaev Interaction in BaCo 2⁢ (AsO 4 ) 2

The inelastic neutron scattering results and their analysis unequivocally point to a dominant Kitaev interaction in the honeycomb-lattice cobaltate BaCo 2 ⁢(AsO 4 ) 2 . Our anisotropic-exchange model closely describes all available neutron scattering data in the material’s field-polarized phase. Furthermore, the density-matrix renormalization group results for our model are in close accord with the unusual double-zigzag magnetic order and the low in-plane saturation field of BaCo 2 ⁢(AsO 4 ) 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1991

Preliminary calibration plan for the Advanced Particles and Field Observatory (APAFO) magnetometer experiment

Prelaunch and postlaunch calibration plans for the APAFO magnetometer experiment are presented. A study of tradeoffs between boom length and spacecraft field is described; the results are summarized. The prelaunch plan includes: calibration of the Vector Fluxgate Magnetometer (VFM), Star Sensors, and Scalar Helium Magnetometer (SHM); optical bench integration; and acquisition of basic spacecraft field data. Postlaunch calibration has two phases. In phase one, SHM data are used to calibrate the VFM, total vector magnetic field data are used to calibrate a physical model of the spacecraft field, and both calibrations are refined by iteration. In phase two, corrected vector data are transformed into geocentric coordinates, previously undetected spacecraft fields are isolated, and initial geomagnetic field models are computed. Provided the SHM is accurate to the required 1.0 nT and can be used to calibrate the VFM to the required 3.0- nT accuracy, the tradeoff study indicates that a 12 m boom and a spacecraft field model uncertainty of 5 percent together allow the 1.0 nT spacecraft field error requirement to be met.

Voorhies, C. V.↗
NASA NTRS2012

Broadband Spectroscopy Using Two Suzaku Observations of the HMXB GX 301-2

We present the analysis of two Suzaku observations of GX 301-2 at two orbital phases after the periastron passage. Variations in the column density of the line-of-sight absorber are observed, consistent with accretion from a clumpy wind. In addition to a CRSF, multiple fluorescence emission lines were detected in both observations. The variations in the pulse profiles and the CRSF throughout the pulse phase have a signature of a magnetic dipole field. Using a simple dipole model we calculated the expected magnetic field values for different pulse phases and were able to extract a set of geometrical angles, loosely constraining the dipole geometry in the neutron star. From the variation of the CRSF width and energy, we found a geometrical solution for the dipole, making the inclination consistent with previously published values.

Suchy, Slawomir↗
NASA NTRS1991

Ordering and phase transitions in random-field Ising systems

An exact analysis of the Ising model with infinite-range interactions in a random field and a local mean-field theory in three dimensions is carried out leading to a phase diagram with several coexistence surfaces and lines of critical points. The results show that the phase diagram depends crucially on whether the distribution of random fields is symmetric or not. Thus, Ising-like phase transitions in a porous medium (the asymmetric case) are in a different universality class from the conventional random-field model (symmetric case).

Maritan, Amos↗
NASA NTRS1989

Novel local symmetries and chiral-symmetry-broken phases in S = 1/2 triangular-lattice Heisenberg model

Using a nonmean-field approach the triangular-lattice S = 1/2 Heisenberg antiferromagnet with nearest- and next-nearest-neighbor couplings is shown undergo an Ising-type phase transition into a chiral-symmetry-broken phase (Kalmeyer-Laughlin-like state) at small T. Removal of next-nearest-neighbor coupling introduces a local Z2 symmetry, thereby suppressing any finite-T chiral order.

Baskaran, G.↗
NASA NTRS2023

Breaking Wave Field Statistics With A Multi-Layer Model

The statistics of breaking wave fields are characterised within a novel multi-layer framework, which generalises the single-layer Saint-Venant system into a multi-layer and non-hydrostatic formulation of the Navier–Stokes equations. We simulate an ensemble of phase-resolved surface wave fields in physical space, where strong nonlinearities, including directional wave breaking and the subsequent highly rotational flow motion, are modelled, without surface overturning. We extract the kinematics of wave breaking by identifying breaking fronts and their speed, for freely evolving wave fields initialised with typical wind wave spectra. The Λ(c) distribution, defined as the length of breaking fronts (per unit area) moving with speed c to c + dc following Phillips (J. Fluid Mech., vol. 156, 1985, pp. 505–531), is reported for a broad range of conditions. We recover the Λ(c) ∝ c−6 scaling without wind forcing for sufficiently steep wave fields. A scaling of Λ(c) based solely on the root-mean-square slope and peak wave phase speed is shown to describe the modelled breaking distributions well. The modelled breaking distributions are in good agreement with field measurements and the proposed scaling can be applied successfully to the observational data sets. The present work paves the way for simulations of the turbulent upper ocean directly coupled to a realistic breaking wave dynamics, including Langmuir turbulence, and other sub-mesoscale processes.

Jiarong Wu↗
DOE OSTI2025

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

The influence of cooling rates on strain phase diagrams and domain structures of ferroelectric thin films: A case study of PbTiO 3

Strain engineering has been established as an effective approach to control phase equilibria, domain configurations, and functional properties of ferroelectric thin films. Temperature-strain phase diagrams have been used as powerful tools for providing insights into strain engineering. However, almost all existing phase diagrams established using the phase-field approach assume quenching conditions without considering actual cooling rates during the post-deposition annealing process of ferroelectric thin films. Within this work, we systematically investigate the influence of cooling rates on domain structures and the strain-phase diagram of ferroelectric thin films using phase-field simulations, taking PbTiO 3 thin films as a model system. We found that both the position of phase boundaries in the strain phase diagrams and the domain morphology are significantly influenced by the cooling rates. It is revealed that while the paraelectric-ferroelectric phase boundary remains invariant, the phase boundaries between single-phase and multi-phase regions tend to shift toward the corresponding multi-phase region as the cool rate reduces. Slow cooling generally leads to more ordered domain structures with increased domain size. Using the obtained equilibrium domain structures, we calculated effective thermal conductivities and found significant variations that can be tuned by the cooling rates. In conclusion, this work reveals an underexplored yet critical impact of cooling rates on phase equilibria and domain structures in ferroelectric thin films, which may inspire further fine-tuning of domains and domain walls in low-dimensional ferroelectrics for multifunctional applications.

Domain structures↗
NASA NTRS1993

CDAW 9 analysis of magnetospheric events on May 3, 1986 - Event C

An intense geomagnetic substorm event on May 3, 1986, occurring toward the end of a strong storm period, is studied. The auroral electrojet indices and global imaging data from both the Northern and Southern Hemispheres clearly revealed the growth phase and expansion phase development for a substorm with an onset at 0111 UT. An ideally located constellation of four spacecraft allowed detailed observation of the substorm growth phase in the near-tail region. A realistic time-evolving magnetic field model provided a global representation of the field configuration throughout the growth and early expansion phase of the substorm. Evidence of a narrowly localized substorm onset region in the near-earth tail is found. This region spread rapidly eastward and poleward after the 0111 UT onset. The results are consistent with a model of late growth phase formation of a magnetic neutral line. This reconnection region caused plasma sheet current diversion before the substorm onset and eventually led to cross-tail current disruption at the time of the substorm onset.

Baker, D. N.↗
NASA NTRS1977

Assessment of analytical techniques for predicting solid propellant exhaust plumes

The calculation of solid propellant exhaust plume flow fields is addressed. Two major areas covered are: (1) the applicability of empirical data currently available to define particle drag coefficients, heat transfer coefficients, mean particle size and particle size distributions, and (2) thermochemical modeling of the gaseous phase of the flow field. Comparisons of experimentally measured and analytically predicted data are made. The experimental data were obtained for subscale solid propellant motors with aluminum loadings of 2, 10 and 15%. Analytical predictions were made using a fully coupled two-phase numerical solution. Data comparisons will be presented for radial distributions at plume axial stations of 5, 12, 16 and 20 diameters.

Tevepaugh, J. A.↗
DOE OSTI2025

Flux effects in precipitation under irradiation simulation of Fe-Cr alloys

Radiation-enhanced precipitation of Cr-rich α’ in irradiated Fe-Cr alloys, which results in hardening and embrittlement, depends on the irradiating particle and the displacement per atom (dpa) rate. Here, we utilize a Cahn-Hilliard phase-field based approach, that includes simple models for nucleation, irradiating particle and rate dependent radiation-enhanced diffusion and cascade mixing to simulate α’ evolution under neutrons, heavy ions, and electron irradiations. Different irradiating particles manifest very different cascade mixing efficiencies. The model was calibrated using neutron data. For cascade inducing neutron/heavy-ion dpa rates at 300 °C between 10-8 and 10-6 dpa/s the model predicts approximately constant number density, decreasing radius, decreasing α’ Cr composition, and lower α’ volume fraction. The model then predicts a dramatic transition to no α’ formation above approximately 10-5 dpa/s, while electron irradiation, with weak mixing, had little effect at dpa rates up to 10-3 dpa/s. These model predictions are consistent with experiments. We explain the results in terms of the flux dependence of the radiation enhanced diffusion, cascade mixing, and their ratio, which all vary significantly in relevant flux ranges for neutron and cascade inducing ion irradiations. These results show that both cascade mixing and radiation enhanced diffusion must be accounted for when attempting to emulate neutron-irradiation effects using accelerated ion irradiations.

Ke, Jia-Hong↗
NASA NTRS1984

Numerical simulation of reconnection in an emerging magnetic flux region

Numerical solutions in two dimensions are presented for the resistive MHD equations of an initial boundary value problem, simulating reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. This numerical model displays four main phases, which are interpreted in terms of (1) a slowly evolving quasi-steady phase during which most of the magnetic flux emerges, with approximate equilibrium between magnetic and pressure forces; (2) an impulsive phase in which either the reconnection or continuing emergence of the first phase comes to disrupt the equilibrium, and extensive acceleration occurs as the high pressure region within the emerging region drives the fluid upwards and outwards; (3) a second quasi-steady phase; and (4) a potential-static phase in which continuing reconnection, ohmic dissipation, and fluid transport through boundaries depletes the system of all currents and flows.

Forbes, T. G.↗
DOE OSTI2025

Unit-density SU(3) Fermi-Hubbard model with spin-flavor imbalance

The advent of ultracold alkaline-earth atoms in optical lattices has established a platform for investigating correlated quantum matter with SU⁡(𝑁) symmetry, offering highly tunable model parameters that allow experiments to access phenomena that are unavailable in conventional materials. Understanding the ground-state physics of SU⁡(𝑁) Fermi-Hubbard models away from the Heisenberg limit and from the spin-flavor balanced setting is important, as examining the flavor imbalance reveals new physics in Fermi-Hubbard models and shows how SU⁡(𝑁) phases react to practical experimental imperfections in optical lattices. Here, in this study, mean-field phase diagrams are presented for the unit-density SU(3) Fermi-Hubbard model at two sets of flavor densities, ($\frac{1}{3}$ − 𝛿, $\frac{1}{3}$ + 𝛿, $\frac{1}{3}$) and ($\frac{1}{4}$ − 𝛿, $\frac{1}{4}$ +𝛿, $\frac{1}{2}$), with the flavor imbalance introduced as 𝛿. Novel phases are identified at moderate interaction strengths for both densities, and their robustness is investigated in the presence of flavor imbalance. Furthermore, we provide microscopic explanations of the phases found and their stability. Analysis of thermal ensembles of random mean-field solutions indicate that, at temperatures accessible in state-of-the-art cold atom experiments, some spin orders are hard for conventional scattering or local observable measurements to detect, but can be more accessible with quantum gas microscopy in optical lattice experiments. This work also shows that nesting and Mottness, intertwined in the usual SU(2) Hubbard model in stark contrast to generic materials, can be tuned in the SU(3) model and play distinct roles. The resulting phase diagrams not only deepen our understanding of SU⁡(𝑁) Fermi-Hubbard models but also inform future experimental search for new phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1983

Spinning mode acoustic radiation from the flight inlet

A mathematical model was developed for spinning mode acoustic radiation from a thick wall duct without flow. This model is based on a series of experiments (with and without flow). A nearly pure azimuthal spinning mode was isolated and then reflection coefficients and far field pressure (amplitude and phase) were measured. In our model the governing boundary value problem for the Helmholtz equation is first converted into an integral equation for the unknown acoustic pressure over a disk, S1, near the mouth of the duct and over the exterior surface, S2, of the duct. Assuming a pure azimuthal mode excitation, the azimuthal dependence is integrated out which yields an integral equation over the generator C1 of S1 and the generator C2 of S2. The sound pressure on C1 was approximated by a truncated modal expansion of the interior acoustic pressure. Piecewise linear spline approximation on C2 was used.

Moss, W. F.↗
NASA NTRS2006

Overview of the Field Phase of the NASA Tropical Cloud Systems and Processes (TCSP)Experiment

The Tropical Cloud Systems and Processes experiment is sponsored by the National Aeronautics and Space Administration (NASA) to investigate characteristics of tropical cyclone genesis, rapid intensification and rainfall using a three-pronged approach that emphasizes satellite information, suborbital observations and numerical model simulations. Research goals include demonstration and assessment of new technology, improvements to numerical model parameterizations, and advancements in data assimilation techniques. The field phase of the experiment was based in Costa Rica during July 2005. A fully instrumented NASA ER-2 high altitude airplane was deployed with Doppler radar, passive microwave instrumentation, lightning and electric field sensors and an airborne simulator of visible and infrared satellite sensors. Other assets brought to TCSP were a low flying uninhabited aerial vehicle, and a surface-based radiosonde network. In partnership with the Intensity Forecasting Experiment of the National Oceanic and Atmospheric Administration (NOAA) Hurricane Research Division, two NOAA P-3 aircraft instrumented with radar, passive microwave, microphysical, and dropsonde instrumentation were also deployed to Costa Rica. The field phase of TCSP was conducted in Costa Rica to take advantage of the geographically compact tropical cyclone genesis region of the Eastern Pacific Ocean near Central America. However, the unusual 2005 hurricane season provided numerous opportunities to sample tropical cyclone development and intensification in the Caribbean Sea and Gulf of Mexico as well. Development of Hurricane Dennis and Tropical Storm Gert were each investigated over several days in addition to Hurricane Emily as it was close to Saffir-Simpson Category 5 intensity. An overview of the characteristics of these storms along with the pregenesis environment of Tropical Storm Eugene in the Eastern Pacific will be presented.

Hood, Robbie E.↗
NASA NTRS2012

Wavefront Sensing for WFIRST with a Linear Optical Model

In this paper we develop methods to use a linear optical model to capture the field dependence of wavefront aberrations in a nonlinear optimization-based phase retrieval algorithm for image-based wavefront sensing. The linear optical model is generated from a ray trace model of the system and allows the system state to be described in terms of mechanical alignment parameters rather than wavefront coefficients. This approach allows joint optimization over images taken at different field points and does not require separate convergence of phase retrieval at individual field points. Because the algorithm exploits field diversity, multiple defocused images per field point are not required for robustness. Furthermore, because it is possible to simultaneously fit images of many stars over the field, it is not necessary to use a fixed defocus to achieve adequate signal-to-noise ratio despite having images with high dynamic range. This allows high performance wavefront sensing using in-focus science data. We applied this technique in a simulation model based on the Wide Field Infrared Survey Telescope (WFIRST) Intermediate Design Reference Mission (IDRM) imager using a linear optical model with 25 field points. We demonstrate sub-thousandth-wave wavefront sensing accuracy in the presence of noise and moderate undersampling for both monochromatic and polychromatic images using 25 high-SNR target stars. Using these high-quality wavefront sensing results, we are able to generate upsampled point-spread functions (PSFs) and use them to determine PSF ellipticity to high accuracy in order to reduce the systematic impact of aberrations on the accuracy of galactic ellipticity determination for weak-lensing science.

Jurling, Alden S.↗