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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS2008

Energy and Mass Transport of Magnetospheric Plasmas during the November 2003 Magnetic Storm

Intensive energy and mass transport from the solar wind across the magnetosphere boundary is a trigger of magnetic storms. The storm on 20-21 November 2003 was elicited by a high-speed solar wind and strong southward component of interplanetary magnetic field. This storm attained a minimum Dst of -422 nT. During the storm, some of the solar wind particles enter the magnetosphere and eventually become part of the ring current. At the same time, the fierce solar wind powers strong outflow of H+ and O+ from the ionosphere, as well as from the plasmasphere. We examine the contribution of plasmas from the solar wind, ionosphere and plasmasphere to the storm-time ring current. Our simulation shows, for this particular storm, ionospheric O+ and solar wind ions are the major sources of the ring current particles. The polar wind and plasmaspheric H+ have only minor impacts. In the storm main phase, the strong penetration of solar wind electric field pushes ions from the geosynchronous orbit to L shells of 2 and below. Ring current is greatly intensified during the earthward transport and produces a large magnetic depression in the surface field. When the convection subsides, the deep penetrating ions experience strong charge exchange loss, causing rapid decay of the ring current and fast initial storm recovery. Our simulation reproduces very well the storm development indicated by the Dst index.

Fok, Mei-Chging↗

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

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

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

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

NASA NTRS1988

Numerical simulation of cloud and precipitation structure during GALE IOP-2

A regional scale model, LAMPS (Limited Area Mesoscale Prediction System), is used to investigate cloud and precipitation structure that accompanied a short wave system during a portion of GALE IOP-2. A comparison of satellite imagery and model fields indicates that much of the large mesoscale organization of condensation has been captured by the simulation. In addition to reproducing a realistic phasing of two baroclinic zones associated with a split cold front, a reasonable simulation of the gross mesoscale cloud distribution has been achieved.

Robertson, F. R.↗

NASA NTRS1990

Long-term compression effects on elastomeric O-ring behavior

This paper discusses the results of testing performed on elastomeric seal materials that had been under compression for extended periods of time. Elastomeric seals used in the Space Shuttle redesigned solid rocket motors can experience compression times of up to six months. These seals must be capable of sealing internal motor pressure upon ignition. The tests described herein were performed in order to verify that the seals, which had experienced long-term compression could seal throughout motor operation. Testing was divided into two phases: (1) dynamic high pressure testing, and (2) resiliency testing. Dynamic testing was performed using specialized test fixtures that allowed simulation of the field joint movements during initial motor operation along with high pressure gas. Resiliency testing was performed using specialized test fixtures that also simulated field joint movements and also had the ability to measure the sealing force of the O-ring. Results from all testing indicated that the current elastomeric seals used in the redesigned solid rocket motors will seal during motor operations in the currently defined launch environments.

Clinton, R. G.↗

NASA NTRS1990

Electron heating in quasi-perpendicular shocks - A Monte Carlo simulation

To study the problem of electron heating in quasi-perpendicular shocks, under the combined effects of 'reversible' motion, in the shock electric potential and magnetic field, and wave-particle interactions a diffusion equation was derived, in the drift (adiabatic) approximation and it was solved by using a Monte Carlo method. The results show that most of the observations can be explained within this framework. The simulation has also definitively shown that the electron parallel temperature is determined by the dc electromagnetic field and not by any wave particle induced heating. Wave-particle interactions are effective in smoothing out the large gradients in phase space produced by the 'reversible' motion of the electrons, thus producing a 'cooling' of the electrons. Some constraints on the wave-particle interaction process may be obtained from a detailed comparison between the simulation and observations. In particular, it appears that the adiabatic approximation must be violated in order to explain the observed evolution of the perpendicular temperature.

Veltri, Pierluigi↗

NASA NTRS2018

Segregation and Phase Transformations Along Superlattice Intrinsic Stacking Faults in Ni-Based Superalloys

In this study, local chemical and structural changes along superlattice intrinsic stacking faults combine to represent an atomic-scale phase transformation. In order to elicit stacking fault shear, creep tests of two different single crystal Ni-base superalloys, ME501 and CMSX-4, were performed near 750 degrees Centigrade using stresses of 552 megapascals and 750 megapascals, respectively. Through high resolution Scanning Transmission Electron Microscopy (STEM) and state-of-the-art energy dispersive x-ray spectroscopy, ordered compositional changes were measured along SISFs (Superlattice Intrinsic Stacking Faults) in both alloys. For both instances, the elemental segregation and local crystal structure present along the SISFs are consistent with a nanoscale to D019 (ordering of a hexagonal close-packed crystal) phase transformation. Another notable observation is prominent Cr- and Co-rich Cottrell atmospheres and new evidence of more complex reordering processes responsible for the formation of these faults. These findings are further supported using density functional theory calculations and High Angle Annular Dark Field (HAADF) STEM image simulations.

STEM-EDS↗

NASA NTRS2014

Forced Reconnection in the Near Magnetotail: Onset and Energy Conversion in PIC and MHD Simulations

Using two-dimensional particle-in-cell (PIC) together with magnetohydrodynamic (MHD) Q1 simulations of magnetotail dynamics, we investigate the evolution toward onset of reconnection and the subsequent energy transfer and conversion. In either case, reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. The boundary deformation leads to the formation of thin embedded current sheets, which are bifurcated in the near tail, converging to a single sheet farther out in the MHD simulations. The thin current sheets in the PIC simulation are carried by electrons and are associated with a strong perpendicular electrostatic field, which may provide a connection to parallel potentials and auroral arcs and an ionospheric signal even prior to the onset of reconnection. The PIC simulation very well satisfies integral entropy conservation (intrinsic to ideal MHD) during this phase, supporting ideal ballooning stability. Eventually, the current intensification leads to the onset of reconnection, the formation and ejection of a plasmoid, and a collapse of the inner tail. The earthward flow shows the characteristics of a dipolarization front: enhancement of Bz, associated with a thin vertical electron current sheet in the PIC simulation. Both MHD and PIC simulations show a dominance of energy conversion from incoming Poynting flux to outgoing enthalpy flux, resulting in heating of the inner tail. Localized Joule dissipation plays only a minor role.

two-dimensional↗

NASA NTRS2002

Contrail Formation in Aircraft Wakes Using Large-Eddy Simulations

In this work we analyze the issue of the formation of condensation trails ("contrails") in the near-field of an aircraft wake. The basic configuration consists in an exhaust engine jet interacting with a wing-tip training vortex. The procedure adopted relies on a mixed Eulerian/Lagrangian two-phase flow approach; a simple micro-physics model for ice growth has been used to couple ice and vapor phases. Large eddy simulations have carried out at a realistic flight Reynolds number to evaluate the effects of turbulent mixing and wake vortex dynamics on ice-growth characteristics and vapor thermodynamic properties.

Paoli, R.↗

NASA NTRSDate not supplied

Sidelobe Reduction and Mode-Purity Enhancement of Vortex Beams from a Programmable, Rectangular Phased Array Antenna

Tapering the magnitude of the electric field excitation of a programmable, rectangular phased array has been employed to achieve sidelobe reduction and enhancement of mode purity in vortex beams, which carry non-zero orbital angular momentum. Far-field radiation patterns pertaining to a commercial-off-the-shelf, 256-element antenna were generated via electromagnetic simulation. The patterns were decomposed into weighted sums of Laguerre-Gaussian modes. The extracted coupling coefficients were used to gauge mode purity. In comparison to a uniform excitation of the phased array elements, our studies indicate that employing a tapering technique can result in a 10.7 dB reduction of sidelobe power and a 7.22 dB improvement in worst-case, mode-to-mode cross-talk.

orbital angular momentum, phased array, beamformin↗

NASA NTRS1991

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

NASA NTRS2012

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

NASA NTRS2019

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↗

NASA NTRS1988

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

NASA NTRS2007

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