Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase field simulation”

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 487 records · Page 27

Spatially Separated Electron and Proton Beams in a Simulated Solar Coronal Jet

Magnetic reconnection is widely accepted to be a major contributor to non-thermal particle acceleration in the solar atmosphere. In this paper we investigate particle acceleration during the impulsive phase of a coronal jet, which involves bursty reconnection at a magnetic null point. A test-particle approach is employed, using electromagnetic fields from a magnetohydrodynamic simulation of such a jet. Protons and electrons are found to be accelerated non-thermally both downwards toward the domain’s lower boundary and the solar photosphere, and outwards along the axis of the coronal jet and into the heliosphere. A key finding is that a circular ribbon of particle deposition on the photosphere is predicted, with the protons and electrons concentrated in different parts of the ribbon. Furthermore, the outgoing protons and electrons form two spatially separated beams parallel to the axis of the jet, signatures that may be observable inin-situobservations of the heliosphere.

Ross Pallister↗

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work is limited to very simple models. In this project, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. The proposed work is based on previously funded NASA work that allowed us to develop a two-dimensional numerical method for boiling flows and to demonstrate the ability of the method to simulate film boiling. While numerical simulations of multi-fluid flows have been advanced in a major way during the last five years, or so, similar capability for flows with phase change are still in their infancy. Although the feasibility of the proposed approach has been demonstrated, it has yet to be extended and applied to fully three-dimensional simulations. Here, a fully three-dimensional, parallel, grid adaptive code will be developed. The numerical method will be used to study nucleate boiling in microgravity, with particular emphasis on two aspects of the problem: 1) Examination of the growth of bubbles at a wall nucleation site and the instabilities of rapidly growing bubbles. Particular emphasis will be put on accurately capturing the thin wall layer left behind as a bubble expands along a wall, on computing instabilities on bubble surfaces as bubbles grow, and on quantifying the effects of both these phenomena on heat transfer; and 2) Examination of the effect of shear flow on bubble growth and heat transfer.

Tryggvason, Gretar↗

Scheduling Improvements Following the Phase 1 Field Evaluation of the ATD-2 Integrated Arrival, Departure, and Surface Concept

NASA is conducting the Airspace Technology Demonstration-2 to evaluate an Integrated Arrival, Departure, and Surface (IADS) traffic management system that extends traffic sequencing for the entire life-cycle of a flight from departure gate to arrival gate within multi-airport, metroplex environments. After development and testing in human-in-the-loop simulations, the IADS system was deployed to Charlotte Douglas International Airport for a three-year field evaluation. From the initial IADS concept development through the end of the Phase 1 field evaluation many lessons were learned with regards to the IADS scheduler. In this paper we describe how data from the Phase 1 field evaluation helped identify scheduler improvements and guided the implementation of refinements. The improvements in the IADS scheduler described in this paper are incorporated into the IADS Phase 2 scheduler enabling strategic Surface Metering Programs and will be evaluated during the field evaluation.

Coupe, William J.↗

Phase field modeling of hydrogen release in nitrile-butadiene rubber composites after high-pressure hydrogen exposure

Experimental data showed that a simple diffusion model doesn’t suffice to explain the hydrogen (H 2 ) release kinetics in nitrile butadiene rubber (NBR) composites (one with no filler or plasticizer, one with plasticizer only, one with fillers and plasticizer, and one with fillers only) after exposed to hydrogen gas at a pressure of 27.6 MPa. In this work, we developed a phase field model that considers the effect of H 2 gas bubble evolution on H 2 release kinetics. With the model, we simulated the effect of thermodynamic and kinetics properties on the nucleation, growth, and shrinkage of gas bubbles and H 2 release kinetics during and after decompression. Here, the results demonstrated that 1) the model built upon certain thermodynamic and kinetics properties of a given material can well describe its H 2 release kinetics measured in experiments, and 2) the model can predict the gas bubble evolution, which is associated with the material property degradation and failures, with accurate thermodynamic and kinetics properties or H 2 release data from experiments.

08 HYDROGEN↗

Shape-driven, emergent behavior in active particle mixtures

Abstract Active particle systems can vary greatly from one-component systems of spheres to mixtures of particle shapes at different composition ratios. We investigate computationally the combined effect of anisotropy and stoichiometry on the collective behavior of two-dimensional active colloidal mixtures of polygons. We uncover three emergent phenomena not yet reported in active Brownian particle systems. First, we find that mixtures containing hexagons exhibit micro-phase separation with large grains of hexagonal symmetry. We quantify a measurable, implicit ‘steric attraction’ between the active particles as a result of shape anisotropy and activity. This calculation provides further evidence that implicit interactions in active systems, even without explicit attraction, can lead to an effective preferential attraction between particles. Next, we report stable fluid clusters in mixtures containing one triangle or square component. We attribute the fluidization of the dense cluster to the interplay of cluster destabilizing particles, which introduce grain boundaries and slip planes into the system, causing solid-like clusters to break up into fluid clusters. Third, we show that fluid clusters can coexist with solid clusters within a sparse gas of particles in a steady state of three coexisting phases. Our results highlight the potential for a wide variety of behavior to be accessible to active matter systems and establish a route to control active colloidal systems through simple parameter designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Scheduling Improvements Following the Phase 1 Field Evaluation of the ATD-2 Integrated Arrival, Departure, and Surface Concept

NASA is conducting the Airspace Technology Demonstration-2 (ATD-2) to evaluate an Integrated Arrival, Departure, and Surface (IADS) traffic management system that extends traffic sequencing for the entire life-cycle of a flight from departure gate to arrival gate within multi-airport, metroplex environments. After development and testing in human-in-the-loop simulations, the IADS system was deployed to Charlotte Douglas International Airport for a three-year field evaluation. From the initial IADS concept development through the end of the Phase 1 field evaluation many lessons were learned with regards to the IADS scheduler. In this paper we describe how data from the Phase 1 field evaluation helped identify scheduler improvements and guided the implementation of refinements. The improvements in the IADS scheduler described in this paper are incorporated into the IADS Phase 2 scheduler enabling strategic Surface Metering Programs and will be evaluated during the field evaluation.

Airspace Technology Demonstration-2↗

Accelerate microstructure evolution simulation using graph neural networks with adaptive spatiotemporal resolution

Abstract Surrogate models driven by sizeable datasets and scientific machine-learning methods have emerged as an attractive microstructure simulation tool with the potential to deliver predictive microstructure evolution dynamics with huge savings in computational costs. Taking 2D and 3D grain growth simulations as an example, we present a completely overhauled computational framework based on graph neural networks with not only excellent agreement to both the ground truth phase-field methods and theoretical predictions, but enhanced accuracy and efficiency compared to previous works based on convolutional neural networks. These improvements can be attributed to the graph representation, both improved predictive power and a more flexible data structure amenable to adaptive mesh refinement. As the simulated microstructures coarsen, our method can adaptively adopt remeshed grids and larger timesteps to achieve further speedup. The data-to-model pipeline with training procedures together with the source codes are provided.

36 MATERIALS SCIENCE↗

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗

Some characteristics of electric field momentum coupling with the neutral atmosphere

A three-dimensional model has been developed to describe momentum coupling between high-latitude electric fields, neutral winds, temperature, and composition. The Hall drag is found to be the main source for atmospheric winds and the small divergence component of winds is due to the Pedersen drag and the Hall drag. Adiabatic heat transfer is responsible for the back pressure which damps the divergence field and for the reversal in circulation of the divergence field at higher altitudes. Back pressure causes a decrease in total wind velocity of about 10% at exospheric heights and by a factor of about 2 at 120 km. The wind field with the pressure feedback may be simulated by neglecting pressure variations and the Coriolis force. Density variations of Ar, N2, O, and He, induced by the momentum source, are in phase above 120 km and out of phase with the temperature amplitude above 150 km. The electrostatic field momentum source is ineffective for directly inducing density and temperature variations in the upper thermosphere.

Mayr, H. G.↗

PT symmetry, pattern formation, and finite-density QCD

A longstanding issue in the study of quantum chromodynamics (QCD) is its behavior at nonzero baryon density, which has implications for many areas of physics. The path integral has a complex integrand when the quark chemical potential is nonzero and therefore has a sign problem, but it also has a generalized ${\mathscr{PT}}$ symmetry. We review some new approaches to ${\mathscr{PT}}$-symmetric field theories, including both analytical techniques and methods for lattice simulation. We show that ${\mathscr{PT}}$-symmetric field theories with more than one field generally have a much richer phase structure than their Hermitian counterparts, including stable phases with patterning behavior. The case of a ${\mathscr{PT}}$-symmetric extension of a φ 4 model is explained in detail. The relevance of these results to finite density QCD is explained, and we show that a simple model of finite density QCD exhibits a patterned phase in its critical region.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Global MHD Simulations of the Substorm Current Wedge and Dipolarization

This paper presents results from global MHD simulations showing the evolution of the plasma and field in the near-Earth tail during the substorm phases. The late growth phase is characterized by pronounced thinning of the plasma sheet and stretching of the field in the region between approximately -6 R(sub E) to -30 R(sub E). A pre-existing X-line moves tailward to beyond -50 R(sub E). Close to onset, a new X-line forms near -18 R(sub E) in the midnight sector. Earthward flows emanating from this X-line dipolarize the near-Earth field, leading to a reduction of the cross-tail current in the midnight sector, but not elsewhere. The magnetic shear between the dipolarized field near midnight and the stretched field elsewhere is equivalent to currents flowing through the ionosphere in a region I sense, and so forming the current wedge. Later in the expansion phase, the dipolarization spreads in local time at a rate of about 0.3 hours MLT per minute. A strong electric field and a rapid increase of the plasma pressure is associated with the dipolarization. Near midnight the dipolarization appears to occur at all distances between 6.6 and 13 R(sub E) at the same time within the resolution (+/- 2 min) of our model. However, the model results indicate that dipolarization starts before ground onset in the pre-midnight sector and propagates both earthward and eastward.THus, dipolarization may be much more complex than simple earthward/tailward and/or azimuthal expansion.

Raeder, J.↗

Convective Systems Over the South China Sea: Cloud-Resolving Model Simulations

The South China Sea Monsoon Experiment (SCSMEX) was conducted in May-June 1998. One of its major objectives is to better understand the key physical processes for the onset and evolution of the summer monsoon over Southeast Asia and southern China. Multiple observation platforms (e.g., upper-air soundings, Doppler radar, ships, wind profilers, radiometers, etc.) during SCSMEX provided a first attempt at investigating the detailed characteristics of convective storms and air pattern changes associated with monsoons over the South China Sea region. SCSMEX also provided rainfall estimates which allows for comparisons with those obtained from the Tropical Rainfall Measuring Mission (TRMM), a low earth orbit satellite designed to measure rainfall from space. The Goddard Cumulus Ensemble (GCE) model (with 1-km grid size) is used to understand and quantify the precipitation processes associated with the summer monsoon over the South China Sea. This is the first (loud-resolving model used to simulate precipitation processes in this particular region. The GCE-model results captured many of the observed precipitation characteristics because it used a fine grid size. For example, the temporal variation of the simulated rainfall compares quite well to the sounding-estimated rainfall variation. The time and domain-averaged temperature (heating/cooling) and water vapor (drying/ moistening) budgets are in good agreement with observations. The GCE-model-simulated rainfall amount also agrees well with TRMM rainfall data. The results show there is more evaporation from the ocean surface prior to the onset of the monsoon than after the on-~et of monsoon when rainfall increases. Forcing due to net radiation (solar heating minus longwave cooling) is responsible for about 25% of the precipitation in SCSMEX The transfer of heat from the ocean into the atmosphere does not contribute significantly to the rainfall in SCSMEX. Model sensitivity tests indicated that total rain production is reduced 17-18% in runs neglecting the ice phase. The SCSMEX results are compared to other GCE-model-simulated weather systems that developed during other field campaigns (i.e., west Pacific warm pool region, eastern Atlantic region and central USA). Large-scale forcing vie temperature and water vapor tendency, is the major energy source for net condensation in the tropical cases. The effects of large-scale cooling exceed that of large-scale moistening in the west pacific warm pool region and eastern Atlantic region. For SCSMEX, however, the effects of large-scale moistening predominate. Net radiation and sensible and latent hc,it fluxes play a much more important role in the central USA.

Tao, Wei-Kuo↗

Light control of surface–bulk coupling by terahertz vibrational coherence in a topological insulator

The demand for disorder-tolerant quantum logic and spin electronics can be met by generating and controlling dissipationless spin currents protected by topology. Dirac fermions with helical spin-locking surface transport offer a way of achieving such a goal. Yet, surface-bulk coupling can lead to strong Dirac electron scattering with bulk carriers and phonons as well as impurities, assisted by such dissipative channel, which results in “topological breakdown”. Here, we demonstrate that coherent lattice vibrations periodically driven by a single-cycle terahertz (THz) pulse can significantly suppress such dissipative channel in topological insulators. This is achieved by reducing the phase space in the bulk available for Dirac fermion scattering into during coherent lattice oscillations in Bi 2 Se 3 . This light-induced suppression manifests as a remarkable transition exclusively in surface transport, absent for bulk, above the THz electric fields for driving coherent phonons, which prolongs the surface transport lifetime. These results, together with simulations, identify the critical role of spin–orbit coupling for the “phase space contraction” mechanism that suppresses the surface-bulk coupling. Imposing vibrational quantum coherence into topological states of matter may become a universal light control principle for reinforcing the symmetry-protected helical transport.

36 MATERIALS SCIENCE↗

Simulations of Arctic Mixed-Phase Boundary Layer Clouds: Advances in Understanding and Outstanding Questions

An overview of mixed-phase boundary layer cloud simulations emphasizes what detailed studies reveal about the level of understanding of microphysical processes, using analogous liquid-phase boundary layer clouds as a reference for the dynamical conditions. We focus primarily on modeling results from three major field campaign case studies, spanning a range of liquid water path, aerosol loading, cloud temperatures, and active processes (including drizzle, aggregation, and riming). Simulations with detailed microphysics are able to reproduce basic aspects of all three case studies, including continuous ice formation within the context of a well-mixed liquid-cloud-topped layer. However, simulations constrained with measurements of ice nuclei generally underpredict the observed amount of ice present, indicating a possible role for ice multiplication that is poorly established. Assumed ice properties may also bear a significant influence on the water vapor budget and reflectivity properties, but remain imprecisely quantified on a case study basis.

Arctic clouds; Mixed-phase clouds; Mixed-phase mic↗

Quantitative Simulations of Siloxane Adsorption in Metal–Organic Frameworks

We present a transferable force field (FF) for simulating the bulk properties of linear and cyclic siloxanes and the adsorption of these species in metal–organic frameworks (MOFs). Unlike previous FFs for siloxanes, our FF accurately reproduces the vapor–liquid equilibria of each species in the bulk phase. The quality of our FF combined with the Universal Force Field using standard Lorentz–Berthelot combining rules for MOF atoms was assessed in a wide range of MOFs without open metal sites, showing good agreement with dispersion-corrected density functional theory calculations. Predictions with this FF show good agreement with the limited experimental data for siloxane adsorption in MOFs that is available. As an example of using the FF to predict adsorption properties in MOFs, we present simulations examining entropy effects in binary linear and cyclic siloxane mixture coadsorption in the large-pore MOF with structure code FOTNIN.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unraveling fundamental mechanisms of silicon nitride crystallization in microelectronics manufacturing

This research project investigates the fundamental mechanisms of silicon nitride (SiN) crystallization, aiming to enhance the understanding of this critical material in microelectronics manufacturing. Through a collaborative effort between Sandia National Laboratories, the University of Tennessee, and the University of Florida, we developed a comprehensive framework that integrates experimental techniques, atomistic modeling, meso-scale simulations, and an integrated multi-scale model to capture this physical phenomenon on multiple time and length scales . The project developed a new machine learning based atomistic potential and utilized advanced phase field modeling to capture the complexities of polycrystalline growth and the influence of mechanical stresses on crystallization dynamics. By employing a grain tracker algorithm, the meso-scale model effectively identified and tracked individual crystal grains, enabling the simulation of anisotropic growth behaviors reflective of SiN’s physical properties. The integration of atomistic simulations with meso-scale modeling created a powerful multi-scale framework that validated atomistic inputs and enhanced predictive accuracy for crystallization dynamics at larger scales, validated experimentally. This adaptable modeling capability not only accelerates development times by informing manufacturing processes but also serves as a valuable starting point for understanding crystallization physics in similar materials. The insights gained from this research unlock new opportunities for the development of advanced materials tailored for future microelectronics and photonics applications. Overall, this project represents a significant advancement in understanding of fundamental physics of SiN and establishes a foundation for future research in material science, bridging the gap between atomic-level phenomena and macroscopic material behaviors for practical applications.

36 MATERIALS SCIENCE↗

Pressure-Induced Phase Transformations

The study of phase transitions in solids under high pressure conditions is a very active and vigorous research field. In recent decades, thanks to the development of experimental techniques and computer simulations, a plethora of important discoveries has been made under high-pressure conditions. Many of the achievements accomplished in recent years a affect various research fields, from solid-state physics to chemistry, materials science, and geophysics. They not only contribute to a deeper understanding of solid-sold phase transitions but also to a better understanding of melting under compression.

36 MATERIALS SCIENCE↗

Phase instability-coupled fracture behavior in garnet LLZO solid electrolytes: a machine learning-enabled atomistic study

Fracture in the garnet-type solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) poses a critical threat to both the performance and safety of solid-state batteries. To unravel the coupled chemomechanical processes that govern fracture in LLZO under external loading, we carry out large-scale molecular-dynamics simulations with a validated machine-learning force field. Our results demonstrate that triaxial stresses at crack flanks trigger a localized cubic-to-tetragonal phase transformation, which is accompanied by Li-ion rearrangement. The emergent tetragonal domains feature lattice contraction normal to the fracture plane, imposing coherent misfit strains that provide an additional driving force for further crack propagation. Crucially, introducing Li deficiencies stabilizes the cubic phase, postponing the phase transition and thereby delaying fracture initiation. These findings highlight the role of intrinsic phase instability in dictating LLZO's fracture resistance and its critical connection to local Li concentration. This chemomechanical coupling points toward targeted strategies to enhance the mechanical robustness of garnet electrolytes, including tuning Li content, ensuring dopant homogeneity, and refining processing protocols.

Materials science↗