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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Gamma-Ray Emission in Dissipative Pulsar Magnetospheres: from Theory to Fermi Observations

We compute the patterns of gamma-ray emission due to curvature radiation in dissipative pulsar magnetospheres. Our ultimate goal is to construct macrophysical models that are able to reproduce the observed gamma-ray light curve phenomenology recently published in the Second Fermi Pulsar Catalog. We apply specific forms of Ohm's law on the open field lines using a broad range for the macroscopic conductivity values that result in solutions ranging, from near-vacuum to near-force-free. Using these solutions, we generate model gamma-ray light curves by calculating realistic trajectories and Lorentz factors of radiating particles under the influence of both the accelerating electric fields and curvature radiation reaction. We further constrain our models using the observed dependence of the phase lags between the radio and gamma-ray emission on the gamma-ray peak separation. We perform a statistical comparison of our model radio-lag versus peak-separation diagram and the one obtained for the Fermi standard pulsars. We find that for models of uniform conductivity over the entire open magnetic field line region, agreement with observations favors higher values of this parameter. We find, however, significant improvement in fitting the data with models that employ a hybrid form of conductivity, specifically, infinite conductivity interior to the light cylinder and high but finite conductivity on the outside. In these models the gamma-ray emission is produced in regions near the equatorial current sheet but modulated by the local physical properties. These models have radio lags near the observed values and statistically best reproduce the observed light curve phenomenology. Additionally, they also produce GeV photon cut-off energies.

stars↗
DOE OSTI2024

Machine learning for the identification of phase transitions in interacting agent-based systems: A Desai-Zwanzig example

Deriving closed-form analytical expressions for reduced-order models, and judiciously choosing the closures leading to them, has long been the strategy of choice for studying phase- and noise-induced transitions for agent-based models (ABMs). In this paper, we propose a data-driven framework that pinpoints phase transitions for an ABM—the Desai-Zwanzig model—in its mean-field limit, using a smaller number of variables than traditional closed-form models. To this end, we use the manifold learning algorithm Diffusion Maps to identify a parsimonious set of data-driven latent variables, and we show that they are in one-to-one correspondence with the expected theoretical order parameter of the ABM. We then utilize a deep learning framework to obtain a conformal reparametrization of the data-driven coordinates that facilitates, in our example, the identification of a single parameter-dependent ordinary differential equation (ODE) in these coordinates. Additionally, we identify this ODE through a residual neural network inspired by a numerical integration scheme (forward Euler). We then use the identified ODE—enabled through an odd symmetry transformation—to construct the bifurcation diagram exhibiting the phase transition.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2004

Time Changes of the European Gravity Field from GRACE: A Comparison with Ground Measurements from Superconducting Gravimeters and with Hydrology Model Predictions

We investigate the time-variable gravity changes in Europe retrieved from the initial GRACE monthly solutions spanning a 18 month duration from April 2002 to October 2003. Gravity anomaly maps are retrieved in Central Europe from the monthly satellite solutions we compare the fields according to various truncation levels (typically between degree 10 and 20) of the initial fields (expressed in spherical harmonics to degree 120). For these different degrees, an empirical orthogonal function (EOF) decomposition of the time-variable gravity field leads us to its main spatial and temporal characteristics. We show that the dominant signal is found to be annual with an amplitude and a phase both in agreement with predictions in Europe modeled using snow and soil-moisture variations from recent hydrology models. We compare these GRACE gravity field changes to surface gravity observations from 6 superconducting gravimeters of the GGP (Global Geodynamics Project) European sub-network, with a special attention to loading corrections. Initial results suggest that all 3 data sets (GRACE, hydrology and GGP) are responding to annual changes in near-surface water in Europe of a few microGal (at length scales of approx.1000 km) that show a high value in winter and a summer minimum. We also point out that the GRACE gravity field evolution seems to indicate that there is a trend in gravity between summer 2002 and summer 2003 which can be related to the 2003 heatwave in Europe and its hydrological consequences (drought). Despite the limited time span of our analysis and the uncertainties in retrieving a regional solution from the network of gravimeters, the calibration and validation aspects of the GRACE data processing based on the annual hydrology cycle in Europe are in progress.

Hinderer, J.↗
DOE OSTI2025

XY-like incommensurate magnetic order in Ce 2 ⁢SnS 5

We report the synthesis of single crystals of Ce 2 ⁢SnS 5 through a two-stage chemical vapor transport method. The Ce 2 ⁢SnS 5 system is a member of the orthorhombic Pbam (No. 55) space group and realizes a distorted trigonal tricapped prism (TTP) crystal field around each cerium site. We characterized the sample through orientation-dependent magnetization and heat capacity measurements to probe the magnetic anisotropy in the system characteristic of XY-like anisotropic Heisenberg model behavior. Ce 2 ⁢SnS 5 furthermore enters a zero-field ordered phase under 𝑇 𝑁 =2.4 K; powder neutron diffraction measurements reveal incommensurate magnetic order near 𝑇 𝑁 . Furthermore, the system then locks into a commensurate, two-𝑞 magnetic structure below approximately 1.2 K. This commensurate structure belongs to the Shubnikov group 𝑃⁢𝑏′⁢𝑎′⁢𝑚′ ⁡(MSG 55.359) and realizes the propagation vectors $\overrightarrow{𝑞}$ = (1/3, 0, 0) and $\overrightarrow{𝑞}$ = (0, 0, 0).

Antiferromagnetism↗
DOE OSTI2024

Defect Production and Microstructural Feature Impact for Radiation Damage in Additively Manufactured 316 Stainless Steel

This milestone presents multi-scale modeling research results for additively manufactured 316 stainless steel. A combination of phase field, cluster dynamics, molecular dynamics, and density functional theory with machine learning is used, allowing for predictions of radiation-driven microstructural evolution in additively manufactured 316 stainless steel over a range of temperatures, damage rates, neutron spectra, and microstructures, and supporting the development of combined ion and neutron irradiation for material qualification. Informed by ion irradiation and neutron irradiation results across the Advanced Materials and Manufacturing Technologies program, we investigate the unique aspects of radiation-driven microstructure evolution in additively manufactured 316 stainless steel. In particular, we focus on understanding the impact of carbon concentration (varying, for example, between the 316L and 316H standards) on void formation; radiation-induced segregation at dislocation cells and grain boundaries; and the evolution of dislocation loops and network dislocation populations. We find that the unique characteristics of the additively manufactured microstructures must be accounted for in understanding the evolution of dislocation populations under thermal and irradiation conditions, such as the variation in sink strengths arising due to the variation in dislocation density. We also find that the radiation-induced segregation of Cr and Ni to grain boundaries and cell walls differs due to the differences in their defect sink biases. We also find that increasing the Ni content can slow vacancy diffusion, which may provide a mechanism for the observed reduction in transient swelling rate for austenitic Fe-Cr-Ni alloys with increasing Ni content. In addition, ion irradiations have shown that increasing carbon content in 316 SS results in a larger population of smaller voids, suggesting reduced vacancy diffusion. Our results show that the carbon content of additively manufactured 316 SS has a significant impact on the migration rate of defect clusters. The presence of carbon atoms results in carbon-vacancy trapping, significantly reducing the diffusion rate of vacancies. Carbon atoms may also be trapped near the surface of a void, which may reduce void growth by trapping vacancies that diffuse toward the void.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1972

Magnetism, dimensional changes, and magnetic transitions in hydrated cesium manganese chloride

Dimensional changes (strain) along the three principal crystal axes of the antiferromagnet CsMnCl3-2H2O are studied as a function of magnetic field and temperature in the antiferromagnetic, spin flopped, and paramagnetic phases. Changes in dimensions through the phase transitions between the magnetic states are examined. By applying the molecular field model and utilizing all available information, magnetic properties of CsMnCl3-2H2O are determined. The possible usefulness of this material in a magnetic refrigeration cycle is evaluated.

Woollam, J. A.↗
NASA NTRS1981

Semi-actuator disk theory for compressor choke flutter

A mathematical anaysis predict the unsteady aerodynamic utilizing semi actuator theory environment for a cascade of airfoils harmonically oscillating in choked flow was developed. A normal shock is located in the blade passage, its position depending on the time dependent geometry, and pressure perturbations of the system. In addition to shock dynamics, the model includes the effect of compressibility, interblade phase lag, and an unsteady flow field upstream and downstream of the cascade. Calculated unsteady aerodynamics were compared with isolated airfoil wind tunnel data, and choke flutter onset boundaries were compared with data from testing of an F100 high pressure compressor stage.

Micklow, J.↗
DOE OSTI2024

The Impact of the 8–10 March 2012 Geomagnetic Storm on Inner Zone Protons as Measured by Van Allen Probes

The Relativistic Electron Proton Telescope (REPT) instrument on the Van Allen Probes observed a double-peaked inner zone proton population throughout the 7 year lifetime of the mission. M. Hudson et al. (2023) showed that a strong SEP event accompanied by a CME-shock in early March 2012 provided the Solar Energetic Proton (SEP) source for the higher L trapped proton population, which then diffused radially inward to be observed by REPT at L = 2. The study followed trajectories of SEP protons launched isotropically from a sphere at 7 Re for 2.5 hr in fields calculated by the LFM-RCM global MHD model, which includes electric fields needed to model the transport and trapping of the protons by the shock, and then a radial diffusion simulation was run for 2 years using the result from the test-particle simulation as the initial condition. The simulation result was compared with REPT measurement in November 2013 and showed reasonable agreement. However, the simulation overestimated the Phase Space Density by a factor of four due to lack of field line curvature scattering during the storm in the model. In this study, a test-particle simulation is performed for 2 days following the injection and trapping of protons in March 2012 using TS05 fields to simulate the field line curvature scattering of the trapped SEP due to the buildup of the ring current during the geomagnetic storm. The resulting sample distribution was then weighted using the flux at the end of the two-hour MHD-test particle simulation. A radial diffusion simulation is then run using the initial profile that included the loss effect, with improved comparison with REPT measurements after 2 years.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Phased Array Characterization of Slat Noise Radiation from a High-Lift Common Research Model

A test campaign was conducted in the open-jet test section of the NASA Langley 14- by 22-Foot Subsonic Tunnel on a 10%-scale, semispan version of a High-Lift Common Research Model (CRM-HL) incorporating a leading-edge slat, trailing-edge flap and removable high-fidelity main landing gear. Far field acoustic measurements were obtained on baseline and acoustically treated model configurations using a traversing 97-microphone phased array that viewed the pressure side of the airframe. A modified version of the DAMAS deconvolution method incorporating corrections for shear layer acoustic wave decorrelation was employed to determine the locations and strengths of relevant noise sources along the span of the slat and in the vicinity of the flap edges and landing gear. The main goal of the test campaign was to use the array data to evaluate the effectiveness of various slat noise reduction concepts. It was found that slat-gap fillers produced significant broadband noise reduction for all Mach numbers and airfoil angles of attack that were investigated. Mach number scaling revealed an approximate velocity to the sixth power relationship for the slat noise at higher frequencies. Directivity measurements obtained from integration of DAMAS pressure-squared values over defined geometric zones around the model showed that baseline slat noise radiation approximated a dipole as a function of polar emission angle, becoming approximately omnidirectional with the introduction of the slat-gap filler. The test campaign clearly demonstrated the noise reduction benefits that can be obtained by applying appropriate treatments to leading edge slats on commercial transport-class aircraft.

William M Humphreys, Jr.↗
NASA NTRS1989

Generation of macroscopic magnetic-field-aligned electric fields by the convection surge ion acceleration mechanism

The 'convection surge' model for ion acceleration, designed by Mauk (1986) to explain the observed ion distributions and the field-aligned character of middle magnetospheric ion distributions during the expansion phase of a substorm, was extended to include the self-consistent generation of magnetic-field-aligned electric fields. Results from the modified model show that the convection surge mechanism leads to the generation of dynamical macroscopic magnetic field-aligned electric fields that begin their strongest developments very near the magnetic equator and then propagate to higher latitudes. Potential drops as high as 1 to 10 kV might be expected, depending on the mass species of the ions and on the electron temperatures. It is speculated that the convection surge mechanism could be a key player in the transient field-aligned electromagnetic processes observed to operate within the middle magnetosphere.

Mauk, B. H.↗
DOE OSTI2026

Anomalies of coset non-invertible symmetries

Anomalies of global symmetries provide important information on the quantum dynamics. We show the dynamical constraints can be organized into three classes: genuine anomalies, fractional topological responses, and integer responses that can be realized in symmetry-protected topological (SPT) phases. Coset symmetry can be present in many physical systems including quantum spin liquids, and the coset symmetry can be a non-invertible symmetry. We introduce twists in coset symmetries, which modify the fusion rules and the generalized Frobenius-Schur indicators. We call such coset symmetries twisted coset symmetries, and they are labeled by the quadruple (G, K, ω D+1 , α D ) in D spacetime dimensions where G is a group and K ⊂ G is a discrete subgroup, ω D+1 is a (D + 1)-cocycle for group G, and α D is a D-cochain for group K. We present several examples with twisted coset symmetries using lattice models and field theory, including both gapped and gapless systems (such as gapless symmetry-protected topological phases). We investigate the anomalies of general twisted coset symmetry, which presents obstructions to realizing the coset symmetry in (gapped) symmetry-protected topological phases. We show that finite coset symmetry G/K becomes anomalous when G cannot be expressed as the bicrossed product G = H ⋈ K , and such anomalous coset symmetry leads to symmetry-enforced gaplessness in generic spacetime dimensions. We illustrate examples of anomalous coset symmetries with A 5 / \Z 2 symmetry, with realizations in lattice models.

Hsin, Po-Shen [King’s College London (United Kingd↗
DOE OSTI2026

Microstructure Validation of Graph Theory Model-Derived Cooling Rates in the Wire Arc Additive Manufacturing of ER70S-6 Steel

Wire arc additive manufacturing (WAAM) enables high-rate fabrication of large metallic components, but spatial variations in thermal history can lead to microstructural heterogeneity that requires efficient process models to evaluate. This study evaluates whether cooling rates extracted from a graph theory model (GTM)-based thermal simulation are consistent with the microstructural evolution observed in an ER70S-6 WAAM wall. Thermal histories from the model were analyzed at selected build heights, and cooling rates were extracted from the final thermal excursion through the austenite phase field. Microstructures at corresponding locations were characterized using electron backscatter diffraction (EBSD) to quantify grain size distributions, and pearlite interlamellar spacing was used as an additional indicator of cooling behavior. The modeled cooling rates were highest near the substrate and generally decreased with build height, consistent with the observed reduction in the fine grain fraction and the progressive shift in the grain size distribution as build height increased. Pearlite spacing trends also supported the modeled cooling rate variation. These results indicate that GTM-derived thermal histories can be post-processed into metallurgically meaningful cooling rate estimates for WAAM steel builds and linked to dataset specific empirical grain size distribution relationships for process–thermal history–microstructure assessment.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Phase V of the single-breath washout test

A downward-deflecting phase V is often seen following the terminal rise (phase IV) in single-breath washout tests. To investigate the nature of phase V, experiments using simultaneous washouts of N2 and tracer boluses of Ar were performed on eight normal nonsmoking subjects aged 27-41 who exhibited a phase V. Phase V is found to occur in all subjects at expiratory flow rates between 0.1 and 2.0 l/sec shortly after expiration became flow limited. Volumes of both phases IV and V increase with increasing flow rate. The difference between the exhaled volumes at which flow became limited and phase V appeared is shown to be approximately equal to the anatomic dead space. Results support a model of lung emptying in a gravitational field in which flow limitation occurs first in the lower lung regions and then progresses toward the upper regions, causing phase IV. A decrease in the amount of flow from the upper relative to the lower regions after all regions have become flow limited then causes phase V.

Nichol, G. M.↗
NASA NTRS2015

Next-Generation Ion Thruster Design Tool

Computational tools that accurately predict the performance of electric propulsion devices are highly desirable and beneficial to NASA and the broader electric propulsion community. The current state of the art in electric propulsion modeling relies heavily on empirical data and numerous computational "knobs." In Phase I of this project, Tech-X Corporation developed the most detailed ion engine discharge chamber model that currently exists. This kinetic model simulates all particles in the discharge chamber along with a physically correct simulation of the electric fields. In addition, kinetic erosion models are included for modeling the ion-impingement effects on thruster component erosion. In Phase II, Tech-X developed a user-friendly computer program for NASA and other governmental and industry customers. Tech-X has implemented a number of advanced numerical routines to bring the computational time down to a commercially acceptable level. NASA now has a highly sophisticated, user-friendly ion engine discharge chamber modeling tool.

Stolz, Peter↗
DOE OSTI2025

Data-Efficient Methods for Determining Flory–Huggins χ Parameters in Multicomponent Polymer Formulations

Polymer formulations are essential in diverse applications including personal care products, coatings, paints, adhesives, and plastic materials. Designing these formulations requires navigating large, complex design spaces, where phase and self-assembly behavior critically impact performance. The Flory–Huggins χ parameter, which quantifies segmental miscibility, is widely used to parametrize the excess free energy of mixing in formulation models. In this work, we introduce two data-efficient, top-down methods for estimating χ parameters using the Random Phase Approximation (RPA): (i) Boundary Nonlinear Regression (Boundary-NLR), which fits theoretical spinodal boundaries to experimental phase boundaries, and (ii) Surrogate Model Inverse Parameter Estimation (SMIPE), which uses a Gaussian Process Classifier to fit sparse phase maps via a surrogate model. Both methods allow rapid parametrization of polymer field-theoretic models without the need for additional experiments. We evaluate these approaches on data sets involving polymer–solvent–nonsolvent ternary mixtures and block copolymer–solvent systems, demonstrating their robustness to experimental noise and their relevance for real-world formulation design.

copolymers↗
NASA NTRS1993

Three dimensional modeling of cirrus during the 1991 FIRE IFO 2: Detailed process study

A three-dimensional model of cirrus cloud formation and evolution, including microphysical, dynamical, and radiative processes, was used to simulate cirrus observed in the FIRE Phase 2 Cirrus field program (13 Nov. - 7 Dec. 1991). Sulfate aerosols, solution drops, ice crystals, and water vapor are all treated as interactive elements in the model. Ice crystal size distributions are fully resolved based on calculations of homogeneous freezing of solution drops, growth by water vapor deposition, evaporation, aggregation, and vertical transport. Visible and infrared radiative fluxes, and radiative heating rates are calculated using the two-stream algorithm described by Toon et al. Wind velocities, diffusion coefficients, and temperatures were taken from the MAPS analyses and the MM4 mesoscale model simulations. Within the model, moisture is transported and converted to liquid or vapor by the microphysical processes. The simulated cloud bulk and microphysical properties are shown in detail for the Nov. 26 and Dec. 5 case studies. Comparisons with lidar, radar, and in situ data are used to determine how well the simulations reproduced the observed cirrus. The roles played by various processes in the model are described in detail. The potential modes of nucleation are evaluated, and the importance of small-scale variations in temperature and humidity are discussed. The importance of competing ice crystal growth mechanisms (water vapor deposition and aggregation) are evaluated based on model simulations. Finally, the importance of ice crystal shape for crystal growth and vertical transport of ice are discussed.

Jensen, Eric J.↗
DOE OSTI2024

CFT and lattice correlators near an RG domain wall between minimal models

Conformal interfaces separating two conformal field theories (CFTs) provide maps between different CFTs, and naturally exist in nature as domain walls between different phases. One particularly interesting construction of a conformal interface is the renormalization group (RG) domain wall between CFTs. For a given Virasoro minimal model \mathcal{M}_{k+3,k+2} ℳ k + 3 , k + 2 , an RG domain wall can be generated by a specific deformation which triggers an RG flow towards its adjacent Virasoro minimal model \mathcal{M}_{k+2,k+1} ℳ k + 2 , k + 1 with the deformation turned on over part of the space. An algebraic construction of this domain wall was proposed by Gaiotto in [J. High Energy Phys. 12, 103 (2012)]. In this paper, we will provide a study of this RG domain wall for the minimal case k=2 k = 2 , which can be thought of as a nonperturbative check of the construction. In this case the wall is separating the Tricritical Ising Model (TIM) CFT and the Ising Model (IM) CFT. We will check the analytical results of correlation functions from the RG brane construction with the numerical density matrix renormalization group (DMRG) calculation using a lattice model proposed in [arXiv.1206.1332, Science 344, 280 (2014)], and find a perfect agreement. We comment on possible experimental realizations of this RG domain wall.

Cogburn, Cameron V.↗
NASA NTRS2019

Translating Phased Array Measurements of a Low-Noise Top-Mounted Propulsion Installation for a Supersonic Airliner

A novel translating-phased-array technique was applied during a model-scale validation test of low-noise propulsion concepts for supersonic transport aircraft. The test involved particle image velocimetry, far field and translating phased array acoustic measurements. This paper covers the translating phased array portion. Data were acquired for a series of nozzles with different chevron designs, both uninstalled and installed on a representative aircraft planform. The far field data confirmed the noise reduction potential of top-mounted installation, while the phased array results gave insights into how the reduction occurred and highlighted areas for improvement. This was the inaugural test for the new translating phased array, which measures source maps as seen by observers at different polar angles. Propagation effects through the free jet were documented for the different viewing angles. The source maps obtained were used to understand important aspects of the far field noise measurements.

Beamforming↗