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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

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↗
DOE OSTI2024

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↗
NASA NTRS1985

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.↗
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 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↗
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.↗
NASA NTRS1983

The structure of the heliospheric current sheet - 1978-1982

Continuing the study near solar minimum in 1976-1977 (Hoeksema et al., 1982), the configuration of the heliospheric magnetic field for the period 1976-1982 is calculated using a potential field model. Particular attention is given to the structure during the rising phase, maximum, and early decline of sunspot cycle 21, from 1978 to 1982. Four warps in the current sheet (the boundary between interplanetary magnetic field toward and away from the sun) are seen early in this interval; these give rise to a four-sector structure in the interplanetary magnetic field observed at earth. It is noted that the location of the current sheet changes slowly and extends to a heliographic latitude of approximately 50 deg. The strength of the polar field correction throughout this period is determined and included in the model calculations. The lower latitude magnetic fields become much stronger as the polar fields weaken and reverse polarity near maximum, decreasing the effect of the polar field correction.

Hoeksema, J. T.↗
DOE OSTI2026

Spatially resolved polarization swings in the supermassive binary black hole candidate OJ 287 with first Event Horizon Telescope observations

We present the first Event Horizon Telescope 1.3 mm observations of the supermassive binary black hole candidate OJ 287. The observations achieved an unprecedented angular resolution of 18 μas and reveal significant structural and polarization variability over just five days, marking the shortest timescale on which such changes have been directly imaged in this source. The inner jet exhibits a twisted ridgeline structure, with features displaying apparent superluminal motions up to about 22 c. The linear polarization maps reveal three main polarized features whose electric-vector position angles (EVPAs) change substantially over the time span of our observations, including a component with a radial polarization consistent with being produced by a recollimation shock. Most notably, we directly resolved two innermost jet components whose EVPAs rotate in opposite directions. The faster component, moving at 2.4 ± 0.9 μas/day (17.4 ± 6.5 c), exhibits counterclockwise EVPA swings of roughly 3.7° per day, while the slower component, with a proper motion of 1.4 ± 0.3 μas/day (10.2 ± 2.2 c), rotates clockwise at approximately 2.5° per day. Previous studies inferred helical magnetic fields in AGN jets from time-resolved or integrated polarization variability but lacked the angular resolution to directly image this effect. Our results provide spatially resolved evidence that a helical magnetic field threads the jet’s collimation and acceleration zone, ruling out models based on the superposition of unresolved components. Our analysis suggests that propagating shocks interact with a Kelvin–Helmholtz plasma instability, illuminating different phases of the helical magnetic field and producing the observed polarization spatial and temporal variability. Moreover, our model naturally accounts for the more rapid polarization rotation observed in the faster moving component. Our model predicts even more rapid swings in polarization, which could be tested with future observations featuring a more densely sampled time coverage.

OJ 287↗
NASA NTRS2014

Multi-Scale Computational Modeling of Two-Phased Metal Using GMC Method

A multi-scale computational model for determining plastic behavior in two-phased CMSX-4 Ni-based superalloys is developed on a finite element analysis (FEA) framework employing crystal plasticity constitutive model that can capture the microstructural scale stress field. The generalized method of cells (GMC) micromechanics model is used for homogenizing the local field quantities. At first, GMC as stand-alone is validated by analyzing a repeating unit cell (RUC) as a two-phased sample with 72.9% volume fraction of gamma'-precipitate in the gamma-matrix phase and comparing the results with those predicted by finite element analysis (FEA) models incorporating the same crystal plasticity constitutive model. The global stress-strain behavior and the local field quantity distributions predicted by GMC demonstrated good agreement with FEA. High computational saving, at the expense of some accuracy in the components of local tensor field quantities, was obtained with GMC. Finally, the capability of the developed multi-scale model linking FEA and GMC to solve real life sized structures is demonstrated by analyzing an engine disc component and determining the microstructural scale details of the field quantities.

nickel alloys↗
NASA NTRS2002

A model of hollow cathode plasma chemistry

We have developed a new model of hollow cathode plasma chemistry based on the observation that xenon ion mobility is diffusion limited due to resonant charge exchange reactions. The model shows that vapor phase barium atoms are ionized almost immediately and electric fields accelerate the ions upstream from the emission zone. We have also applied the model to the orifice region, where the resultant ion generation profile correlates with previously reported orifice erosion.

hollow cathode↗
DOE OSTI2024

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

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

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1980

Magnetic variations during substorms

Extensive observations of magnetic field variations during auroral substorms have been used to construct phenomenological models of the substorm process. In one model, three substorm phases, growth, expansion and recovery correspond respectively to energy storage, energy release and return to equilibrium. New data and more detailed studies show the simple model required modification. Consideration is given to a number of recent studies of substorm magnetic field variations which illustrate the limitations of the simple model. These include the existence of several substorm expansions following a single growth phase; the frequent occurrence of multiple onsets of a single substorm; plasma sheet thinning prior to each onset; correlations between northward fluctuations of the IMF and expansion onset. These observations are incorporated into a revised phenomenological model.

Mcpherron, R. L.↗
NASA NTRS2008

Forecasting Lightning Threat using Cloud-Resolving Model Simulations

Two new approaches are proposed and developed for making time and space dependent, quantitative short-term forecasts of lightning threat, and a blend of these approaches is devised that capitalizes on the strengths of each. The new methods are distinctive in that they are based entirely on the ice-phase hydrometeor fields generated by regional cloud-resolving numerical simulations, such as those produced by the WRF model. These methods are justified by established observational evidence linking aspects of the precipitating ice hydrometeor fields to total flash rates. The methods are straightforward and easy to implement, and offer an effective near-term alternative to the incorporation of complex and costly cloud electrification schemes into numerical models. One method is based on upward fluxes of precipitating ice hydrometeors in the mixed phase region at the-15 C 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 domain-wide statistics of the peak values of simulated flash rate proxy fields against domain-wide 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. Our blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Exploratory tests for selected North Alabama cases show that, because WRF can distinguish the general character of most convective events, our methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because the models tend to have more difficulty in 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 forecasts become available.

McCaul, Eugene W., Jr.↗
NASA NTRS1977

A design handbook for phase change thermal control and energy storage devices

Comprehensive survey is given of the thermal aspects of phase change material devices. Fundamental mechanisms of heat transfer within the phase change device are discussed. Performance in zero-g and one-g fields are examined as it relates to such a device. Computer models for phase change materials, with metal fillers, undergoing conductive and convective processes are detailed. Using these models, extensive parametric data are presented for a hypothetical configuration with a rectangular phase change housing, using straight fins as the filler, and paraffin as the phase change material. These data are generated over a range of realistic sizes, material properties, and thermal boundary conditions. A number of illustrative examples are given to demonstrate use of the parametric data. Also, a complete listing of phase change material property data are reproduced herein as an aid to the reader.

Humphries, W. R.↗
NASA NTRS1983

The structure of the heliospheric current sheet, 1978 - 1982

The structure of the heliospheric magnetic field changes substantially during the 11 year sunspot cycle. Its configuration for the period 1976 through 1982 using a potential field model was calculated. The structure during the rising phase, maximum, and early decline of sunspot cycle 21, from 1978 to 1982 is considered.

Hoeksema, J. T.↗