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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Lunar gravitational field estimation and the effects of mismodeling upon lunar satellite orbit prediction

Lunar spherical harmonic gravity coefficients are estimated from simulated observations of a near-circular low altitude polar orbiter disturbed by lunar mascons. Lunar gravity sensing missions using earth-based nearside observations with and without satellite-based far-side observations are simulated and least squares maximum likelihood estimates are developed for spherical harmonic expansion fit models. Simulations and parameter estimations are performed by a modified version of the Smithsonian Astrophysical Observatory's Planetary Ephemeris Program. Two different lunar spacecraft mission phases are simulated to evaluate the estimated fit models. Results for predicting state covariances one orbit ahead are presented along with the state errors resulting from the mismodeled gravity field. The position errors from planning a lunar landing maneuver with a mismodeled gravity field are also presented. These simulations clearly demonstrate the need to include observations of satellite motion over the far side in estimating the lunar gravity field. The simulations also illustrate that the eighth degree and order expansions used in the simulated fits were unable to adequately model lunar mascons.

Davis, John H.↗
DOE OSTI2024

Wave-packet manipulation of He Rydberg states by a seeded free-electron laser

We report a two-dimensional pump-control-probe spectroscopic study of the dynamics of singly excited He Rydberg-state wave packets with a seeded extreme ultraviolet (XUV) free-electron laser (FEL) source. A pair of coherent XUV pulses, defined by their coarse time separation and relative phase, created and manipulated the wave packets. The He atoms were postionized by infrared (IR) pulses, and the ion yield was measured as a function of XUV phase and IR arrival time. We tagged and sorted the relative phase of the XUV pulse pair on a single-shot basis by fitting each FEL spectrum with a suitable function that accounts for nonidealities of the XUV pulse pairs, associated with the seeding process; more generally, the fit returns the time-dependent electric field of the FEL spectra. The experimental two-dimensional maps of ion yields, measured as a function of IR (probe) delay and of XUV (pump-control) phase, were compared with the solution of the first-order time-dependent Schrödinger equation for this field. Despite the fact that the experimental conditions imply strong excitation, beyond the approximations of first-order perturbation theory, the simulated map satisfactorily reproduces the experimental one for temporally well-separated pulses. We show that by selecting data at appropriate values of pump-control phase, we enhance or suppress the amplitude of chosen wave-packet components consisting of two or more Rydberg states. When the temporal overlap of the pulse pair cannot be neglected, the phase reconstruction is underdetermined, and we provide a simplified comparison between data and simulations. Published by the American Physical Society 2024

Dumergue, Mathieu (ORCID:0000000236984406)↗
NASA NTRS2020

Beam Propagation Through Atmospheric Turbulence Using an Altitude-Dependent Structure Profile with Non-Uniformly Distributed Phase Screens

Modeling the effects of atmospheric turbulence on optical beam propagation is a key element in the design and analysis of free-space optical communication systems. Numerical wave optics simulations provide a particularly useful technique for understanding the degradation of the optical field in the receiver plane when the analytical theory is insufficient for characterizing the atmospheric channel. Motivated by such an application, we use a split-step method modeling the turbulence along the propagation path as a series of thin random phase screens with modified von Karman refractive index statistics using the Hufnagel-Valley turbulence profile to determine the effective structure constant for each screen. In this work, we employ a space-to-ground case study to examine the irradiance and phase statistics for both uniformly and non-uniformly spaced screens along the propagation path and compare to analytical results. We find that better agreement with the analytical theory is obtained using a non-uniform spacing with the effective structure constant for each screen chosen to minimize its contribution to the scintillation in the receiver plane. We evaluate this method as a flexible alternative to other standard layered models used in astronomical imaging applications.

Chahine, Yousef K.↗
DOE OSTI2023

EUTERPE: A global gyrokinetic code for stellarator geometry

The current state of the EUTERPE code is described with emphasis on the implemented models and their numerical implementation. The code solves the multi-species electromagnetic gyrokinetic equations in the full volume of a three-dimensional domain. Noise reduction of the particle-in-cell method is achieved by using a δf-method and Fourier filters. The field equations are discretized with B-splines and the resulting system of equations is solved iteratively. For linear simulations a phase-factor transformation is applied in order to strongly reduce the necessary grid resolution. Apart from the full gyrokinetic model, other numerically less expensive hybrid models are also implemented. They are mainly tailored for comparison with fluid theory and for studying the interaction of the bulk plasma with fast particles. The code is parallelized for CPUs by particle and domain decomposition. Good scalability up to several thousand nodes is demonstrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Phase-field approach to cellular blebbing

Bulges in the plasma membrane of cells known as blebs can form spontaneously in a wide range of biological processes, but what controls their shape and stability remains incompletely understood. Here, to address this we introduce a dual phase-field model with coupled order parameters representing the cell cortex and plasma membrane that can quantitatively model blebbing in three dimensions. Simulations and sharp-interface analyses reveal that, depending on whether blebbing occurs by detachment of the plasma membrane or rupture of the actin cortex, blebs can form discontinuously through a saddle-node bifurcation or continuously with increasing cortical tension. The model predictions are in good quantitative agreement with existing experimental data for laser-induced cortex rupture.

Biological and medical sciences↗
DOE OSTI2022

Lefschetz thimble quantum Monte Carlo for spin systems

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Furthermore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Modeling of Grid-Forming and Grid-Following Inverters for Dynamic Simulation of Large-Scale Distribution Systems

Historically, distribution system planning studies mainly focused on steady state and quasi-steady state analysis, with limited attention paid to dynamic analysis. This paper develops three-phase, electromechanical models for both grid-forming and grid-following inverters, and integrates them into an open source, three-phase distribution network solver, thereby enabling dynamic simulation of large-scale, three-phase unbalanced distribution systems with high penetration of inverter-based DERs. The proposed inverter models are validated against electromagnetic simulation and field test data from the CERTS/AEP microgrid testbed, and simulated in an islanded 5252 node distribution system in the GridLAB-D simulation environment. Simulation verifies the effectiveness of the proposed inverter models for large-scale distribution system analysis. Finally, study results show that compared to traditional grid-following inverters, the high penetration of grid-forming inverters can improve the voltage and frequency stability of islanded distribution systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Marmot V2

MARMOT is a robust numerical tool for mesoscale modeling of fuel performance developed under the NEAMS Fuels technical area to predict the coevolution of microstructure and properties in fuel and cladding materials. MARMOT accomplishes this using the phase field method coupled with finite strain mechanics and heat conduction. MARMOT is based on the open source Multiphysics Object-Oriented Simulation Environment (MOOSE) and solves the coupled partial differential equations defining the physics using the finite element method. MARMOT is being developed in order to facilitate the development of improved materials models for fuel performance, but it is also being developed as a powerful tool in and of itself for the simulation of mesoscale fuel performance.

Aagesen, LarryK.↗
DOE OSTI2021

Development of High-Performance Ni-base Alloys for Gas Turbine Wheels Using a Coprecipitation Approach (Final Technical Report)

In this project, a synergistic approach involving experiments informed by modelling was used to develop γ'/γ'' coprecipitation-strengthened Ni-base superalloys for land-based gas turbine wheels. In thick sections of such turbine wheels, the slow cooling rate normally observed, as well as variations in cooling rate could potentially result in drastically inhomogeneous and coarsened microstructure. Thus, ten alloy compositions were selected based on coprecipitation and sluggish γ’ growth and coarsening during three iterations of alloy development, assisted by extensive CALPHAD modeling. Phase equilibrium calculations, informed by the characterization of early alloys, were successful in predicting suitable aging treatments to obtain bimodal coprecipitate distribution in several alloys. These aging treatments were applied to the subsequent round of alloys and the resulting microstructures were characterized in detail. Two different cooling rates were used during each of these heat treatments to simulate the range of cooling rates observed in a thick section forging. The microstructures were used to calibrate a multi-phase field model (MPF), which was used to study the formation mechanisms of various coprecipitate morphologies observed in the experiments. Following this, the most promising alloys in each round were selected and studied using high temperature tensile and creep tests. Considering that a slower cooling rate would lead to a strength debit, while a faster cooling rate might lead to a decline in creep rupture life of a full-scale wheel, tensile tests were carried out for slower cooled samples, while creep tests were carried out for faster cooled samples. The deformation microstructures were then characterized using advanced characterization techniques such as Diffraction Contrast-Scanning Transmission Electron Microscopy (DC-STEM) and Electron Channeling Contrast Imaging (ECCI). Further, a mean-field creep model was developed and calibrated using these experimental findings to predict long-term creep behavior of these coprecipitation-strengthened alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Parabolic Wave Equation Generalized to a Spatially Varying Plasma in a Vertical Magnetic Field

This report walks through the process of deriving the wave equation for a vertical wave being propagated through a spatially varying plasma in a vertical magnetic field. Additionally, this report shows how to derive the index of refraction of the plasma corresponding to three different modes of the vertical magnetic field. Finally, this paper briefly explores the next steps planned to successfully simulate and test the equations in a pre-existing phase screen scintillation code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Defect-phase-dynamics approach to statistical domain-growth problem of clock models

The growth of statistical domains in quenched Ising-like p-state clock models with p = 3 or more is investigated theoretically, reformulating the analysis of Ohta et al. (1982) in terms of a phase variable and studying the dynamics of defects introduced into the phase field when the phase variable becomes multivalued. The resulting defect/phase domain-growth equation is applied to the interpretation of Monte Carlo simulations in two dimensions (Kaski and Gunton, 1983; Grest and Srolovitz, 1984), and problems encountered in the analysis of related Potts models are discussed. In the two-dimensional case, the problem is essentially that of a purely dissipative Coulomb gas, with a sq rt t growth law complicated by vertex-pinning effects at small t.

Kawasaki, K.↗
NASA NTRS2003

Thermocapillary-Induced Phase Separation with Coalescence

Research has been undertaken on interactions of two or more deformable drops (or bubbles) in a viscous fluid and subject to a temperature, gravitational, or flow field. An asymptotic theory for nearly spherical drops shows that small deformations reduce the coalescence and phase separation rates. Boundary-integral simulations for large deformations show that bubbles experience alignment and enhanced coalescence, whereas more viscous drops may break as a result of hydrodynamic interactions. Experiments for buoyancy motion confirm these observations. Simulations of the sedimentation of many drops show clustering phenomena due to deformations, which lead to enhanced phase separation rates, and simulations of sheared emulsions show that deformations cause a reduction in the effective viscosity.

Davis, Robert H.↗
DOE OSTI2021

Numerical Modelling of the Optical Stochastic Cooling Experiment at IOTA

A proof-of-principle optical-stochastic cooling (OSC) experiment is currently in its commissioning phase at the Fermilab’s IOTA ring. In support of this experiment, we recently implemented an OSC element in the ELEGANT tracking program. The model, based on a semi-analytic description of OSC [*], supports the simulation of a large number of macroparticles (10⁴-10⁶) over many turns (10⁶). This paper showcases the simulation capabilities to investigate the beam dynamics in the presence of cooling (or self-interacting radiation field in general) and quantify the impact of various sources of error (e.g. transverse and phase jitter), guide data analysis.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Robust and ultrafast state preparation by ramping artificial gauge potentials

Abstract The implementation of static artificial magnetic fields in ultracold atomic systems has become a powerful tool, e.g. for simulating quantum-Hall physics with charge-neutral atoms. Taking an interacting bosonic flux ladder as a minimal model, we investigate protocols for adiabatic state preparation via magnetic flux ramps. Considering the fact that it is actually the artificial vector potential (in the form of Peierls phases) that can be experimentally engineered in optical lattices, rather than the magnetic field, we find that the time required for adiabatic state preparation dramatically depends on which pattern of Peierls phases is used. This can be understood intuitively by noting that different patterns of time-dependent Peierls phases that all give rise to the same magnetic field ramp, generally lead to different artificial electric fields during the ramp. As an intriguing result, we find that an optimal choice allows for preparing the ground state almost instantaneously in the non-interacting system, which can be related to the concept of counterdiabatic driving. Remarkably, we find extremely short preparation times also in the strongly-interacting regime. Our findings open new possibilities for robust state preparation in atomic quantum simulators.

Wang, Botao (ORCID:0000000282202452)↗
DOE OSTI2025

Thermodynamic Cloud Phase Classifications Using Machine Learning at NSA and ANX

Vertically resolved thermodynamic cloud phase classifications are essential for studies of atmospheric cloud and precipitation processes. The Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) THERMOCLDPHASE Value-Added Product (VAP) uses a multi-sensor approach to classify thermodynamic cloud phase by combining lidar backscatter and depolarization, radar reflectivity, Doppler velocity, spectral width, microwave radiometer-derived liquid water path, and radiosonde temperature measurements. The measured voxels are classified as ice, snow, mixed-phase, liquid (cloud water), drizzle, rain, and liq_driz (liquid+drizzle). We use this product as the ground truth to train three machine learning (ML) models to predict the thermodynamic cloud phase from multi-sensor remote sensing measurements taken at the ARM North Slope of Alaska (NSA) observatory: a random forest (RF), a multilayer perceptron (MLP), and a convolutional neural network (CNN) with a U-Net architecture. Evaluations against the outputs of the THERMOCLDPHASE VAP with one year of data show that the CNN outperforms the other two models, achieving the highest test accuracy, F1-score, and mean Intersection over Union (IOU). Analysis of ML confidence scores shows ice, rain, and snow have higher confidence scores, followed by liquid, while mixed, drizzle, and liq_driz have lower scores. Feature importance analysis reveals that the mean Doppler velocity and vertically resolved temperature are the most influential datastreams for ML thermodynamic cloud phase predictions. The ML models’ generalization capacity is further evaluated by applying them at another Arctic ARM site in Norway using data taken during the ARM Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) field campaign. Finally, we evaluate the ML models’ response to simulated instrument outages and signal degradation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Transient Coupled Chemical-Thermal-Fluid Field Simulation for Sealed Aluminum-clad Spent Nuclear Fuel Storage Canister

As the first step toward developing three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model for unsealed and vented canister storage system, a 3D CFD model coupled with bulk gas radiolysis reactions was developed first for sealed DOE standard canisters filled with inert gas and trace amount of air and water. The workflow for constructing canister-scale 3D CFD models and coupling with gas phase radiolysis reactions were established, which can be readily extended to unsealed, vented canister storage system. This interim milestone report documents the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields and hydrogen gas concentrations within sealed DOE standard canisters over long period of time. In addition, a mesh refinement test was done to show that increasing the models mesh refinement had negligible impact upon the temperature profiles.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1991

Evolution of magnetic flux ropes associated with flux transfer events and interplanetary magnetic clouds

Spacecraft observations suggest that flux transfer events and interplanetary magnetic clouds may be associated with magnetic flux ropes which are magnetic flux tubes containing helical magnetic field lines. In the magnetic flux ropes, the azimuthal magnetic field is superposed on the axial field. The time evolution of a localized magnetic flux rope is studied. A two-dimensional compressible MHD simulation code with a cylindrical symmetry is developed to study the wave modes associated with the evolution of flux ropes. It is found that in the initial phase both the fast magnetosonic wave and the Alfven wave are developed in the flux rope. After this initial phase, the Alfven wave becomes the dominant wave mode for the evolution of the magnetic flux rope and the radial expansion velocity of the flux rope is found to be negligible. Numerical results further show that even for a large initial azimuthal component of the magnetic field, the propagation velocity along the axial direction of the flux rope remains the Alfven velocity. It is also found that the localized magnetic flux rope tends to evolve into two separate magnetic ropes propagating in opposite directions. The simulation results are used to study the evolution of magnetic flux ropes associated with flux transfer events observed at the earth's dayside magnetopause and magnetic clouds in the interplanetary space.

Wei, C. Q.↗
DOE OSTI2022

Multi-frame, ultrafast, x-ray microscope for imaging shockwave dynamics

Inertial confinement fusion (ICF) holds increasing promise as a potential source of abundant, clean energy, but has been impeded by defects such as micro-voids in the ablator layer of the fuel capsules. It is critical to understand how these micro-voids interact with the laser-driven shock waves that compress the fuel pellet. At the Matter in Extreme Conditions (MEC) instrument at the Linac Coherent Light Source (LCLS), we utilized an x-ray pulse train with ns separation, an x-ray microscope, and an ultrafast x-ray imaging (UXI) detector to image shock wave interactions with micro-voids. To minimize the high- and low-frequency variations of the captured images, we incorporated principal component analysis (PCA) and image alignment for flat-field correction. After applying these techniques we generated phase and attenuation maps from a 2D hydrodynamic radiation code (xRAGE), which were used to simulate XPCI images that we qualitatively compare with experimental images, providing a one-to-one comparison for benchmarking material performance. Moreover, we implement a transport-of-intensity (TIE) based method to obtain the average projected mass density (areal density) of our experimental images, yielding insight into how defect-bearing ablator materials alter microstructural feature evolution, material compression, and shock wave propagation on ICF-relevant time scales.

47 OTHER INSTRUMENTATION↗