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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Consistent large-eddy simulation of a temporal mixing layer laden with evaporating drops : Part 1. Direct numerical simulation formulation and a priori analysis

This paper is organized as follows. In § 2, we present the DNS formulation. This formulation consists of a set of coupled conservation equations for the gas phase and the liquid phase (i.e. the drops). The coupling describing the interaction between the two phases is provided in the gas-phase conservation equations by source terms. In § 3, we analyse the DNS results; in particular, we evaluate the adopted form of the heat and mass fluxes. The gas-phase LES equations are developed in § 4 and assumptions that may simplify these equations are examined using the DNS database. Thereafter, the contribution of the source terms is examined by considering the budget of the LES equations. The final form of the LES equations are then presented, in which only the SGS fluxes and the FSTs need be modelled. A parallel perspective on the importance of the source terms is presented in § 5, wherein the irreversible entropy production (i.e. the dissipation) expression for two-phase flow with phase change is derived and evaluated for transitional DNS and the filtered flow fields. Models for the FSTs appearing in the LES equations are investigated in § 6, while several constantcoefficient SGS models are proposed and calibrated in § 7. The present study aims to create simple models so that the modelling approach can be validated by means of an SGS models are proposed and calibrated in § 7. The present study aims to create simple models so that the modelling approach can be validated by means of an a posteriori study.

"Okong'o, Nora A.", 'Bellan, Josette↗
DOE OSTI2025

Theory and simulations of cross-beam energy transfer between speckled laser beams

In this paper, we introduce a realistic modeling of laser beams smoothed by random phase plates to investigate cross-beam energy transfer (CBET). We compute the total field of the smoothed laser beams, which allows us to perform accurate simulations and to build a model accounting for the realistic structure of the speckles. This lets us confirm the results already discussed in a previous work given by Oudin et al. [Phys. Plasmas 29, 112112 (2022)]. When CBET is induced by a plasma flow, plane wave models can predict the exchange between smoothed laser beams. By contrast, when the exchange is caused by a wavelength shift between the laser fields of each beam, plane wave models overestimate the energy transfer, even for strongly Landau-damped ion acoustic waves. We also show that the laser angular spread, which results from its focusing, affects the resonance between the ponderomotive force induced by the laser fields and the ion acoustic wave. It increases the resonance width, which weakly depends on the Landau damping rate when there is a wavelength shift and no flow. In addition, the power transfer obtained with or without a phase plate is essentially the same in the perturbative regime and for small crossing angles.

Laser plasma interactions↗
DOE OSTI2026

Capturing the fractocohesive length scale in elastomers through a statistical mechanics-based gradient enhanced damage model

Here, this study aims to examine modeling flaw sensitivity in elastomers. The direct incorporation of polymer chain statistical mechanics considerations into a continuum stretch-based gradient-enhanced damage formulation, in turn, allows a representation of diffuse chain damage and fracture events that align with known micromechanical mechanisms. Through a series of numerical experiments, we simulate crack propagation and extract the fracture energy as an output of the model, while keeping track of the micromechanical signatures of diffuse chain damage that accommodate fracture propagation and eventually influence flaw sensitivity. Finally, by combining the fracture toughness and the work to rupture, we identify a fractocohesive length of the material, corresponding to the full width of the damage process zone. As the damage-to-fracture cascade in the proposed GED model is influenced by the introduction of a length scale associated with network imperfection and long-range load transfer, the emerging relationship of the two length scales is discussed, providing a potential link between microscopic damage mechanisms and the observed macroscopic fracture response.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Growth-controlled twinning and magnetic anisotropy in CeSb 2

Cerium diantimonide (CeSb 2 ) is a layered heavy-fermion Kondo lattice material that hosts complex magnetism and pressure-induced superconductivity. The interpretation of its in-plane anisotropy has remained unsettled due to structural twinning, which superimposes orthogonal magnetic responses. Here, in this study, we combine controlled crystal growth with magnetization and rotational magnetometry to disentangle the effects of twinning. Nearly untwinned high-quality single crystals reveal the intrinsic in-plane anisotropy: The in-plane easy axis saturates at 𝑀 easy ⁡(4 T) ≈ 1.8 𝜇 B /Ce, while the in-plane hard axis magnetization is strongly suppressed, nearly linear, and comparable to the out-of-plane response. These results resolve long-standing discrepancies in reported magnetic measurements, in which in-plane metamagnetic transition fields and saturation magnetization varied significantly across previous studies. Growth experiments demonstrate that avoiding the proposed 𝛼 𝛽 structural transition—through Sb-rich flux and slower cooling—systematically reduces twinning. However, powder x-ray diffraction and differential thermal analysis measurements show no clear evidence of a distinct 𝛽 phase. Our results establish a consistent magnetic phase diagram and provide essential constraints for crystal-electric field models, enabling a clearer understanding of the interplay between anisotropic magnetism and unconventional superconductivity in CeSb 2 .

Weber, Jan T. [Goethe University, Frankfurt (Germa↗
DOE OSTI2025

Demonstration of controlled spatial incoherence for beam smoothing

A concept for beam smoothing of high-energy laser systems is demonstrated. Multiple coherent beams are multiplexed so that their incoherent mix decreases the far-field nonuniformity after random spatial phase modulation introduced by a phase plate. Modeling shows that the incoherent mix effectively leads to the expected scaling in terms of the number of uncorrelated beams, and that spatial multiplexing can be combined with smoothing by spectral dispersion. The results of an experimental demonstration of this concept with up to six beams, leading to more than 40 modes with polarization multiplexing and angular dispersion, are consistent with expectations. The potential practical advantages of this smoothing scheme and future developments are discussed.

47 OTHER INSTRUMENTATION↗
NASA NTRS2003

High-Fidelity Micromechanics Model Enhanced for Multiphase Particulate Materials

This 3-year effort involves the development of a comprehensive micromechanics model and a related computer code, capable of accurately estimating both the average response and the local stress and strain fields in the individual phases, assuming both elastic and inelastic behavior. During the first year (fiscal year 2001) of the investigation, a version of the model called the High-Fidelity Generalized Method of Cells (HFGMC) was successfully completed for the thermo-inelastic response of continuously reinforced multiphased materials with arbitrary periodic microstructures (refs. 1 and 2). The model s excellent predictive capability for both the macroscopic response and the microlevel stress and strain fields was demonstrated through comparison with exact analytical and finite element solutions. This year, HFGMC was further extended in two technologically significant ways. The first enhancement entailed the incorporation of fiber/matrix debonding capability into the two-dimensional version of HFGMC for modeling the response of unidirectionally reinforced composites such as titanium matrix composites, which exhibit poor fiber/matrix bond. Comparison with experimental data validated the model s predictive capability. The second enhancement entailed further generalization of HFGMC to three dimensions to enable modeling the response of particulate-reinforced (discontinuous) composites in the elastic material behavior domain. Next year, the three-dimensional version will be generalized to encompass inelastic effects due to plasticity, viscoplasticity, and damage, as well as coupled electromagnetothermomechanical (including piezoelectric) effects.

Pindera, Marek-Jerzy↗
DOE OSTI2026

Site-decorated model for unconventional frustrated magnets: Ultranarrow phase crossover and two-dimensional spin reversal transition

Here, the site-decorated Ising model is introduced to advance the understanding and experimental realization of the recently discovered one-dimensional (1D) finite-temperature ultranarrow phase crossover in an external magnetic field, while mitigating the geometric complexities of traditional bond-decorated models. The unconventional frustration and physics are clarified by exactly mapping the 1D site-decorated Ising model in a magnetic field onto a zero-field bond-decorated 𝐽 1 −𝐽 2 Ising model with conventional geometrical frustration. Furthermore, although higher-dimensional Ising models in an external field remain unsolved exactly, an exact solution for a spin-reversal transition—driven by an exotic, hidden half-ice, half-fire state induced by site decoration—is derived. This transition, triggered by a slight variation in temperature or magnetic field—without changing its direction—even in the weak-field limit, offers a promising route toward energy-efficient applications such as data storage and processing. The results suggest that site decoration offers an avenue for materials and device design, particularly in systems such as mixed 𝑑−𝑓 compounds, optical lattices, and neural networks, calling for further studies with site-decorated Heisenberg models. In addition, the site-decorated model offers a rigorous test ground for artificial intelligence (AI) in science, as the analytic derivation of the present results was not only validated but also improved by a general-purpose large language model, inspiring the use of AI as scientific discoverer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2012

Progress in Turbulence Detection via GNSS Occultation Data

The increased availability of radio occultation (RO) data offers the ability to detect and study turbulence in the Earth's atmosphere. An analysis of how RO data can be used to determine the strength and location of turbulent regions is presented. This includes the derivation of a model for the power spectrum of the log-amplitude and phase fluctuations of the permittivity (or index of refraction) field. The bulk of the paper is then concerned with the estimation of the model parameters. Parameter estimators are introduced and some of their statistical properties are studied. These estimators are then applied to simulated log-amplitude RO signals. This includes the analysis of global statistics derived from a large number of realizations, as well as case studies that illustrate various specific aspects of the problem. Improvements to the basic estimation methods are discussed, and their beneficial properties are illustrated. The estimation techniques are then applied to real occultation data. Only two cases are presented, but they illustrate some of the salient features inherent in real data.

Cornman, L. B.↗
NASA NTRS1971

The propagation of solar cosmic-ray bursts.

Description of a model showing analytically the three phases of anisotropy which occur during solar cosmic-ray events in the 7.5 to 21 MeV kinetic-energy interval and reported by McCracken et al. (1971): (1) a highly anisotropic, near field-aligned, initial phase, (2) a convective phase, and (3) a late-time phase in which the anisotropy is approximately perpendicular to the mean interplanetary magnetic field. The model is based on the cosmic-ray particles being convectively transported out from the sun, undergoing anisotropic diffusion along the interplanetary magnetic-field lines, and losing energy by adiabatic deceleration or by collision processes. The event is seen simply as a pulse moving outward from the sun after a cosmic-ray burst with a negative density-gradient in front of it and a positive gradient behind.

Ng, C. K.↗
NASA NTRS1987

Beaming of Uranian kilometric radiation and implications for the planetary magnetic field

Radio waves from Uranus at 481.2 kHz detected by Voyager 2 show a strong beaming effect. When the data are analyzed according to the magnetic latitude of the observer in an offset tilted dipole field model, a sharply defined emission beam extending from 25 to 55 deg latitude is evident consistent with a source over the south magnetic pole. The lack of prominent emission from the north magnetic pole can be attributed to the large spatial offset of the magnetic dipole resulting in much lower source frequencies. For proper alignment of the beam, an advance of the dipole axis longitude by 23 deg from the standard model is required. This phase shift is interpreted as the result of higher order magnetic moments producing a skew of the field lines at low altitudes in the polar source region. This conclusion further substantiates models of the interior of Uranus which predict a complex intrinsic field with significant nondipolar contributions to the surface magnetic field.

Barbosa, D. D.↗
NASA NTRS1986

Production of flickering aurora and field-aligned electron flux by electromagnetic ion cyclotron waves

Recent observations have suggested that flickering aurora is produced by a modulation of the field-aligned component of the electron flux within an auroral arc. It is proposed that a portion of the field-aligned electrons are of ionospheric origin and that these electrons are accelerated and their flux modulated by electromagnetic ion cyclotron waves that occur below the main acceleration region on auroral arc field lines. A model of the electromagnetic ion cyclotron wave shows that the parallel phase velocity of the wave increase as the wave propagates toward the ionosphere. A test particle calculation shows that ionospheric electrons trapped or reflected by the wave are accelerated to energies of several keV and that their flux is modulated at the wave frequency. The relative amplitudes of the model wave electric fields are consistent with the observations of small-scale low-frequency ionospheric and magnetospheric electric fields near auroral arcs of approximately 10 mV/m and 100 mV/m, respectively. The large-amplitude ion cyclotron waves also produce a ponderomotive force and a self-consistent ambipolar electric field. Energy considerations show that the downward energy flux in the electromagnetic ion cyclotron wave can be several percent of the total downward auroral electron energy flux.

Temerin, M.↗
DOE OSTI2025

Light control of intramolecular nuclear dynamics by vortex electron localization

In strong-field ionization of molecules, intense light pulses are thought to have a negligible direct influence on atomic nuclei. Molecular dissociation is thus expected to be determined by the geometrical configuration of the molecular ion at the ionization instant. Contrary to this picture, we observe a counterintuitive electron-proton angular correlation and the formation of proton vortices following strong-field ionization of H 2 molecules by bicircularly polarized two-color laser fields. We explain this phenomenon by the pathway interference and localization of the residual H 2 ⁡+ electron in different angular-momentum states formed in the tail of the driving laser pulse. We validate this interpretation by combining a quantum-mechanical numerical simulation of the field-driven coupled electronic-nuclear dynamics and a semiclassical-trajectory model for the phase accumulation of the laser-driven electronic-nuclear wave packet. Our joint experimental and theoretical study reveals a general picture of vortex electron localization which can be used for controlling molecular-bond breaking with circularly polarized laser fields.

Atomic & molecular processes in external fields↗
NASA NTRS1994

Magnitude of B(sub Z) in the neutral sheet of the magnetotail

Statistical estimates of the average value of B(sub Z) in the magnetotail neutral sheet between 10 and 22 R(sub E) in X, and -10 to 10 R(sub E) in Y (GSM coordinates) are given for different phases of geomagnetic activity. With few exceptions the average value of B(sub Z) is found to be between 5 and 8 nT, but generally near 7 nT, irrespective of the type or phase of activity. This is considerably higher than is predicted in current magnetic field models of the magnetotail. The magnetic field within 2 R(sub E) of the neutral sheet is examined to determine whether B(sub Z) remains approximately constant. We do not find this to be true, B(sub Z) decreasing significantly over this distance from Delta-Z = 0. The structure of the cross-tail current and variations in current density with substorm phase appear to be more complicated than is commonly accepted.

Huang, C. Y.↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Dynamical capture of field stars by interstellar clouds

Evidence suggesting the existence of old stars in young galactic clusters, has motivated an investigation of a model in which such stars are captured from the field during the gas cloud collapse phase of star formation. A hyperbolic field star cannot be gravitationally captured by a static cloud since the star's total energy is conserved. However, if the cloud's potential is an explicit function of time, the energy of a penetrating star can be decreased. As a first step in modeling the dynamical capture of field stars, capture by freely falling, homogeneous spherical cloud is considered. This simplified model is solved, and a prediction is obtained for number of captured field stars in terms of the cloud parameters and dispersion velocities of the field stars.

Whitman, P. G.↗
NASA NTRS2009

Forecasting Lightning Threat using Cloud-resolving Model Simulations

As numerical forecasts capable of resolving individual convective clouds become more common, it is of interest to see if quantitative forecasts of lightning flash rate density are possible, based on fields computed by the numerical model. Previous observational research has shown robust relationships between observed lightning flash rates and inferred updraft and large precipitation ice fields in the mixed phase regions of storms, and that these relationships might allow simulated fields to serve as proxies for lightning flash rate density. It is shown in this paper that two simple proxy fields do indeed provide reasonable and cost-effective bases for creating time-evolving maps of predicted lightning flash rate density, judging from a series of diverse simulation case study events in North Alabama for which Lightning Mapping Array data provide ground truth. One method is based on the product of upward velocity and the mixing ratio of precipitating ice hydrometeors, modeled as graupel only, in the mixed phase region of storms at the -15\dgc\ level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domainwide statistics of the peak values of simulated flash rate proxy fields against domainwide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. A blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Weather Research and Forecast Model simulations of selected North Alabama cases show that this model can distinguish the general character and intensity of most convective events, and that the proposed methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because models tend to have more difficulty in correctly predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models, the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of cloud-allowing forecasts become available.

McCaul, E. W., Jr.↗
NASA NTRS1988

Modelling of SAR polarisation phase difference from trees

The data for polarization phase difference Delta Phi between the HH- and VV-polarized backscattered waves from tree-covered fields were obtained with an airborne synthetic aperture radar at 1.225 GHz. The mean values over tree-covered fields were derived from the images of the phase difference and were examined as a function of incident beam angle from 15 to 55 deg. A theoretical model for simulating these data, based on the electromagnetic wave scatterings from the tree trunk and its branches, both of which are assumed as very long dielectric cylinders was developed. The radius and direction of a tree branch are taken as random variables and are chosen by a Monte Carlo method to encounter the incident waves in producing the scattering events. The Monte Carlo simulated results are in good agreement with the observations within experimental uncertainty.

Mo, Tsan↗
DOE OSTI2025

Spin-flop-like transition as quantum critical point in Cs2⁢RuO4

We report thermodynamic, neutron diffraction, and inelastic neutron scattering measurements on Cs2⁢RuO4, a member of the celebrated family of frustrated magnets Cs2⁢𝑀⁢𝑋4 (𝑀 = Cu, Co, 𝑋 = Br, Cl). Unlike the previously studied members, it is based on 4⁢𝑑 transition metal ions with 𝑆=1. Mapping out the 𝐻−𝑇 magnetic phase diagram reveals an unusual continuous spin-flop-like phase transition associated with a quantum critical point within the antiferromagnetically ordered phase. A quantitative analysis of the complex magnetic excitation spectrum measured in zero field allows us to derive a model magnetic Hamiltonian for this compound. Its main feature is a frustration of magnetic anisotropy on a level that is much higher than in any of the previously studied species. This frustration naturally explains the peculiar phase transition observed.

Nabi, S [Laboratory for Solid State Physics, ETH Z↗