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 613 records · Page 34

Lattice simulations of inflation

The scalar field theory of cosmological inflation constitutes nowadays one of the preferred scenarios for the physics of the early universe. In this paper we aim at studying the inflationary universe making use of a numerical lattice simulation. Various lattice codes have been written in the last decades and have been extensively used for understating the reheating phase of the universe, but they have never been used to study the inflationary phase itself far from the end of inflation (i.e. about 50 e-folds before the end of inflation). In this paper we use a lattice simulation to reproduce the well-known results of some simple models of single-field inflation, particularly for the scalar field perturbation. The main model that we consider is the standard slow-roll inflation with an harmonic potential for the inflaton field. We explore the technical aspects that need to be accounted for in order to reproduce with precision the nearly scale invariant power spectrum of inflaton perturbations. We also consider the case of a step potential, and show that the simulation is able to correctly reproduce the oscillatory features in the power spectrum of this model. Even if a lattice simulation is not needed in these cases, that are well within the regime of validity of linear perturbation theory, this sets the basis to future work on using lattice simulations to study more complicated models of inflation.

79 ASTRONOMY AND ASTROPHYSICS↗

Design and characterization of a W‐band wide‐beam horn for radio holography

Abstract A W‐band wide‐beam horn suitable for prime‐focus feed is described in this paper. It was designed to provide better illuminations for outer panels of larger reflector antennas and to improve performance of radio holography systems. Much effort has been spent on precise characterization of the amplitude and phase patterns of the wide beam horn which has a direct impact on the measurement accuracy of the radio holography. A far‐field measurement setup with the correction of phase errors due to cable flexing has been developed. Good beam pattern agreement has been achieved between simulation and measurement.

Zhang, Xiaoling↗

Studies of the charging of a thin dust layer in a plasma

The present study is designed to extend Wilson's model to include refinements in the actual operational design of the simulation code itself and to include additional physics in order to assess the resulting electric field structure and grain charges. Most of the initial phases of the work have focussed on the fact that the code, though in principle reliable, suffers from severe time constraints, requiring as must as 10-50 hours of CPU time for typical runs on a VAX computer. It is likely that by staging the code into steps in which early steps are carried out crudely and later steps successively fine tuned the total CPU time can be cut significantly. Preliminary results indicate the possibility to cut CPU time to as much as 1/3 or 1/5 the former time. The original Wilson code has been carefully examined, tested and documented. Shortcuts and staging systems have been introduced and tested which should make the code produce more accurate results with a finer grid system than would have been feasible heretofor. Continued testing and implementation of these refinements will be carried out at the author's home institution in the near future. In addition to the effort at reducing CPU time to a manageable amount, work is in progress to include a dust-size distribution instead of assuming all dust grains have the same size.

Peterson, Lennart R.↗

Global Simulation of Proton Precipitation Due to Field Line Curvature During Substorms

The low latitude boundary of the proton aurora (known as the Isotropy Boundary or IB) marks an important boundary between empty and full downgoing loss cones. There is significant evidence that the IB maps to a region in the magnetosphere where the ion gyroradius becomes comparable to the local field line curvature. However, the location of the IB in the magnetosphere remains in question. In this paper, we show simulated proton precipitation derived from the Field Line Curvature (FLC) model of proton scattering and a global magnetohydrodynamic simulation during two substorms. The simulated proton precipitation drifts equatorward during the growth phase, intensifies at onset and reproduces the azimuthal splitting published in previous studies. In the simulation, the pre-onset IB maps to 7-8 RE for the substorms presented and the azimuthal splitting is caused by the development of the substorm current wedge. The simulation also demonstrates that the central plasma sheet temperature can significantly influence when and where the azimuthal splitting takes place.

Gilson, M. L.↗

Optimization of Lyot Coronagrah focal plane masks for improved low-order wavefront sensing

The WFIRST Coronagraph uses a Low-Order Wavefront Sensing and Control (LOWFS/C) subsystem to stabilize the stellar wavefront. This is true for the Hybrid Lyot Coronagraph (HLC), the Shaped Pupil Coronagraph (SPC) spectrograph, and the SPC wide field of view modes. For the SPC spectrograph, we have found via modeling and simulations that a focal plane mask (FPM) with an elliptical phase dimple can reduce the LOWFS errors significantly compared to the conventional, circular phase dimple. For HLC designs, we concurrently optimized the FPM parameters and deformable mirror shapes using our open-source software toolbox, the Fast Linearized Coronagraph Optimizer (FALCO). This process resulted in new FPMs with non-uniform central-area phase modulation. We have found in numerical simulations that such FPMs still provide sufficient LOWFS performance. In this paper, we will present the LOWFS performance for several new HLC and SPC FPM designs.

Riggs, A.J. Eldorado↗

Reduced gravity boiling and condensing experiments simulated with the COBRA/TRAC computer code

A series of reduced-gravity two-phase flow experiments has been conducted with a boiler/condenser apparatus in the NASA KC-135 aircraft in order to obtain basic thermal-hydraulic data applicable to analytical design tools. Several test points from the KC-135 tests were selected for simulation by means of the COBRA/TRAC two-fluid, three-field thermal-hydraulic computer code; the points were chosen for a 25-90 percent void-fraction range. The possible causes for the lack of agreement noted between simulations and experiments are explored, with attention to the physical characteristics of two-phase flow in one-G and near-zero-G conditions.

Cuta, Judith M.↗

A Three-Dimensional MHD Simulation of the Solar Wind for a Tilted-Dipole Magnetic Field on the Sun

Using a three-dimensional MHD model, we simulate the global steady-state structure of the solar corona and solar wind for a dipole magnetic field on the Sun inclined by 30 degrees to the solar rotation axis. This represents the solar conditions typical for a declining phase of solar cycle. The computations can extend from the coronal base out to 100-AU and at large heliospheric distances includes the effects of interstellar neutral hydrogen and their interaction with solar wind protons. The simulations can model the formation of corotating interaction regions and the heliospheric current sheet. The simulations are also capable of describing very strong rarefaction regions that include embedded sub-Alfvenic regions that form on the trailing edge of a fast flows.

Goldstein, Melvyn L.↗

Bipropellant rocket exhaust plume analysis on the Galileo spacecraft

This paper describes efforts to quantify the contaminant flow field produced by 10 N thrust bipropellant rocket engines used on the Galileo spacecraft. The prediction of the composition of the rocket exhaust by conventional techniques is found to be inadequate to explain experimental observations of contaminant deposition on moderately cold (200 K) surfaces. It is hypothesized that low volatility contaminants are formed by chemical reactions which occur on the surfaces. The flow field calculations performed using the direct simulation Monte Carlo method give the expected result that the use of line-of-sight plume shields may have very little effect on the flux of vapor phase contaminant species to a surface, especially if the plume shields are located so close to the engine that the interaction of the plume with the shield is in the transition flow regime. It is shown that significant variations in the exhaust plume composition caused by nonequilibrium effects in the flow field lead to very low concentrations of species which have high molecular weights in the more rarefied regions of the flow field. Recommendations for the design of spacecraft plume shields and further work are made.

Guernsey, C. S.↗

Imaging Field–Driven Melting of a Molecular Solid at the Atomic Scale

Solid–liquid phase transitions are basic physical processes, but atomically resolved microscopy has yet to capture their full dynamics. A new technique is developed for controlling the melting and freezing of self–assembled molecular structures on a graphene field–effect transistor (FET) that allows phase–transition behavior to be imaged using atomically resolved scanning tunneling microscopy. This is achieved by applying electric fields to 2,3,5,6–tetrafluoro–7,7,8,8–tetracyanoquinodimethane–decorated FETs to induce reversible transitions between molecular solid and liquid phases at the FET surface. Nonequilibrium melting dynamics are visualized by rapidly heating the graphene substrate with an electrical current and imaging the resulting evolution toward new 2D equilibrium states. An analytical model is developed that explains observed mixed–state phases based on spectroscopic measurement of solid and liquid molecular energy levels. The observed nonequilibrium melting dynamics are consistent with Monte Carlo simulations.

36 MATERIALS SCIENCE↗

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗

Hybrid simulations of the saturated Farley-Buneman instability in the ionosphere

Numerical simulations of the Farley-Buneman instability in 2-1/2 dimensions using particle ions and fluid electrons show the growth, saturation and nonlinear behavior of two-stream waves. This hybrid technique models the saturated state of the instability for a much longer period of time than the pure particle codes that preceded it. While focusing principally on modeling the topside E region equatorial electrojet, many of these results apply to the auroral two-stream instability as well. The following features are seen in all our hybrid simulations: (1) wave growth at an angle offset from the electron drift direction where the angle depends on the strength of the driving electric field, (2) nonlinear coupling to waves traveling perpendicular to the propagation direction of the principal two-stream waves, (3) a saturated wave phase velocity at or above the sound speed but well below the velocity predicted by linear theory and (4) phase velocities which remain almost constant as a simulated radar sweeps from a horizontal direction to nearly vertical. The nonlinear electron motion dominates the behavior of these waves. Further, these simulations indicate that ion kinetic effects are not essential for the saturation of the instability and that electron temperature effects have a minor impact on the final saturated state.

Oppenheim, Meers↗

Machine-learning force-field models for dynamical simulations of metallic magnets

We review recent advances in machine-learning (ML) force-field methods for Landau–Lifshitz–Gilbert simulations of itinerant electron magnets, focusing on their scalability and transferability. Built on the principle of locality, a deep neural-network model is developed to efficiently and accurately predict electron-mediated forces governing spin dynamics. Symmetry-aware descriptors constructed through a group-theoretical approach ensure rigorous incorporation of both lattice and spin-rotation symmetries. The framework is demonstrated using the prototypical s-d exchange model widely employed in spintronics. ML-enabled large-scale simulations reveal novel nonequilibrium phenomena, including anomalous coarsening of tetrahedral spin order on the triangular lattice and the freezing of phase-separation dynamics in lightly hole-doped, strong-coupling square-lattice systems. These results establish ML force-field frameworks as scalable, accurate, and versatile tools for modeling nonequilibrium spin dynamics in itinerant magnets.

Artificial neural networks↗

Formation of Fourier phase shifts in the solar Ni I 6768 A line

A formalism is developed to understand better how Doppler shifts of spectrum lines as inferred from phase shifts in the Fourier transforms of line profiles are related to the underlying velocity structures which they are intended to measure. With a standard model atmosphere and a simplified, quasi-LTE treatment of line formation, the formalism is applied to the Ni I 6768 A line, which has been selected for use with a network of imaging interferometers under development by the Global Oscillations Network Group for research in helioseismology. Fourier phase shifts are found to be a remarkably linear measure of velocity even in the presence of gradients and unresolved lateral variations in the assumed velocity field. An assumed outward increase in amplitude of a model oscillatory velocity is noticeably reflected in the center-to-limb behavior of the simulated velocity measure, and a sample model of solar granulation is found to have a strong influence on the formation of the Fourier phase.

Jones, Harrison P.↗

PyJMAK: An Open-Source Python Toolkit for Modeling Solid-State Metallurgical Phase Transformations

Accurate prediction of metallurgical phase transformations is an essential basis for autonomous optimization and rapid part qualification. Several methods can be used to estimate the evolution of phase fractions such as JMAK kinetics-based models, phase-field models, thermodynamic models, and data-driven machine learning models. Thermodynamic and phase-field-based methodologies solve multiphysics equations requiring numerous calibration parameters and significant computational resources. As a result, the computation domain is limited to a point or on order of micron-meters. The data-driven models rely on large datasets from experiments and simulations. While the JMAK model only provides information about phase fraction evolution, it can predict this evolution in near real-time using thermal history and thermodynamic data without restriction on the domain. JMAK models have been popularly used by researchers to model phase transformations occuring during additive manufacturing or over arbitrary temperature profiles. Commercial proprietary software such as Abaqus and Ansys or closed-source in-house implementations offer the ability to model JMAK based kinetics to predict phase transformation. However, these software packages are not open-source or freely available for use and development in conjunction with manufacturing machines, sensors, and machine learning algorithms. In addition, the use of the model is restricted by a license token. In contrast, given temperature profiles at multiple points in the domain, this Python-based PyJMAK model can compute phase evolution in parallel due to its stand-alone modular, voxel-based structure, and it can be executed on high-performance computing resources without any license restrictions.

Prabhune, Bhagya [Oak Ridge National Laboratory (O↗

Implementation of real‐time TDDFT for periodic systems in the open‐source PySCF software package

Abstract We present a new implementation of real‐time time‐dependent density functional theory (RT‐TDDFT) for calculating excited‐state dynamics of periodic systems in the open‐source Python‐based PySCF software package. Our implementation uses Gaussian basis functions in a velocity gauge formalism and can be applied to periodic surfaces, condensed‐phase, and molecular systems. As representative benchmark applications, we present optical absorption calculations of various molecular and bulk systems and a real‐time simulation of field‐induced dynamics of a (ZnO) 4 molecular cluster on a periodic graphene sheet. We present representative calculations on optical response of solids to infinitesimal external fields as well as real‐time charge‐transfer dynamics induced by strong pulsed laser fields. Due to the widespread use of the Python language, our RT‐TDDFT implementation can be easily modified and provides a new capability in the PySCF code for real‐time excited‐state calculations of chemical and material systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gas Bubble Evolution in Polycrystalline UMo Fuels Under Elastic-Plastic Deformation: A Phase-Field Model With Crystal-Plasticity

In monolithic UMo fuels, the interaction between the Al cladding and large gas bubble volumetric swelling causes both elastic-plastic and creep deformation. In this work, a phase-field model of gas bubble evolution in polycrystalline UMo under elastic-plastic deformation was developed for studying the dynamic interaction between evolving gas bubble/voids and deformation. A crystal plasticity model, which assumes that the plastic strain rate is proportional to resolved shear stresses of dislocation slip systems on their slip planes, was used to describe plastic deformation in polycrystalline UMo. Xe diffusion and gas bubble evolution are driven by the minimization of chemical and deformation energies in the phase-field model, while evolving gas bubble structure was used to update the mechanical properties in the crystal plasticity model. With the developed model, we simulated the effect of gas bubble structures (different volume fractions and internal gas pressures) on stress-strain curves and the effect of local stresses on gas bubble evolution. The results show that 1) the effective Young’s modulus and yield stress decrease with the increase of gas bubble volume fraction; 2) the hardening coefficient increases with the increase of gas bubble volume fraction, especially for gas bubbles with higher internal pressure; and 3) the pressure dependence of Xe thermodynamic and kinetic properties in addition to the local stress state determine gas bubble growth or shrinkage. The simulated results can serve as a guide to improve material property models for macroscale fuel performance modeling.

36 MATERIALS SCIENCE↗

Gravity field determination for Mars Observer

Mars Observer will be the first near circular, low altitude and short-periodic orbiter of Mars. From previous Mars orbiters, such as Mariner 9 and Vikings 1 and 2, a large gravitational oblateness and in general a gravity field ten times stronger than the earth's field have been determined. Because of these, the gravity field will dominate the evolution of the spacecraft's orbital motion and is the principal error source for reconstructing the spacecraft's motion. Thus a gravity calibration (GC) period has been established prior to the mapping phase in order to refine the gravitational field model. This paper describes the plan being made by Navigation to implement the GC strategy. By analyzing simulated Doppler data acquired during the seven day GC period, we give (1) the expected improvement in the gravity field model and (2) the resulting improvement in reconstructed spacecraft orbital motion. In addition, we have developed a strategy for converging the GC gravity model.

Esposito, Pasquale B.↗

Analysis of the Defect Structure of B2 Feal Alloys

The Bozzolo, Ferrante and Smith (BFS) method for alloys is applied to the study of the defect structure of B2 FeAI alloys. First-principles Linear Muffin Tin Orbital calculations are used to determine the input parameters to the BFS method used in this work. The calculations successfully determine the phase field of the B2 structure, as well as the dependence with composition of the lattice parameter. Finally, the method is used to perform 'static' simulations where instead of determining the ground state configuration of the alloy with a certain concentration of vacancies, a large number of candidate ordered structures are studied and compared, in order to determine not only the lowest energy configurations but other possible metastable states as well. The results provide a description of the defect structure consistent with available experimental data. The simplicity of the BFS method also allows for a simple explanation of some of the essential features found in the concentration dependence of the heat of formation, lattice parameter and the defect structure.

Bozzolo, Guillermo↗