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 325 records · Page 18

A generic multi-flex-body dynamics, controls simulation tool for space station

An order (n) multiflex body Space Station simulation tool is introduced. The flex multibody modeling is generic enough to model all phases of Space Station from build up through to Assembly Complete configuration and beyond. Multibody subsystems such as the Mobile Servicing System (MSS) undergoing a prescribed translation and rotation are also allowed. The software includes aerodynamic, gravity gradient, and magnetic field models. User defined controllers can be discrete or continuous. Extensive preprocessing of 'body by body' NASTRAN flex data is built in. A significant aspect, too, is the integrated controls design capability which includes model reduction and analytic linearization.

London, Ken W.↗

Reversible electron heating vs. wave-particle interactions in quasi-perpendicular shocks

The energy necessary to explain the electron heating in quasi-perpendicular collisionless shocks can be derived either from the electron acceleration in the d.c. cross shock electric potential, or by the interactions between the electrons and the waves existing in the shock. A Monte Carlo simulation has been performed to study the electron distribution function evolution through the shock structure, with and without particle diffusion on waves. This simulation has allowed us to clarify the relative importance of the two possible energy sources; in particular it has been shown that the electron parallel temperature is determined by the d.c. 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.

Veltri, P.↗

The Effects of Neutral Inertia on Ionospheric Currents in the High-Latitude Thermosphere Following a Geomagnetic Storm

Results of an experimental and theoretical investigation into the effects of the time dependent neutral wind flywheel on high-latitude ionospheric electrodynamics are presented. The results extend our previous work which used the National Center for Atmospheric Research Thermosphere/Ionosphere General Circulation Model (NCAR TIGCM) to theoretically simulate flywheel effects in the aftermath of a geomagnetic storm. The previous results indicated that the neutral circulation, set up by ion-neutral momentum coupling in the main phase of a geomagnetic storm, is maintained for several hours after the main phase has ended and may dominate height-integrated Hall currents and field-aligned currents for up to 4-5 hours. We extend the work of Deng et al. to include comparisons between the calculated time-dependent ionospheric Hall current system in the storm-time recovery period and that measured by instruments on board the Dynamics Explorer 2 (DE 2) satellite. Also, comparisons are made between calculated field-aligned currents and those derived from DE 2 magnetometer measurements. These calculations also allow us to calculate the power transfer rate (sometimes called the Poynting flux) between the magnetosphere and ionosphere. The following conclusions have been drawn: (1) Neutral winds can contribute significantly to the horizontal ionospheric current system in the period immediately following the main phase of a geomagnetic storm, especially over the magnetic polar cap and in regions of ion drift shear. (2) Neutral winds drive Hall currents that flow in the opposite direction to those driven by ion drifts. (3) The overall morphology of the calculated field-aligned current system agrees with previously published observations for the interplanetary magnetic field (IMF) B(sub Z) southward conditions, although the region I and region 2 currents are smeared by the TI(ICM model grid resolution. (4) Neutral winds can make significant contributions to the field-aligned current system when B(sub Z) northward conditions prevail following the main phase of a storm, but can account for only a fraction of the observed currents. (5) DE 2 measurements provide a demonstration of "local" (satellite-altitude) flywheel effects. (6) On the assumption that the magnetosphere acts as an insulator, we calculate neutral-wind-induced polarization electric fields of approx. 20-30 kV in the period immediately following the geomagnetic storm.

Deng, W.↗

Contribution of Meteorological Downscaling to Skill and Precision of Seasonal Drought Forecasts

Research in meteorological prediction on sub-seasonal to seasonal (S2S) timescales has seen growth in recent years. Concurrent with this, demand for seasonal drought forecasting has risen. While there is obvious synergy between these fields, S2S meteorological forecasting has typically focused on low resolution global models, while the development of drought can be sensitive to the local expression of weather anomalies and their interaction with local surface properties and processes. This suggests that downscaling might play an important role in the application of meteorological S2S forecasts to skillful forecasting of drought. Here, we apply the Generalized Analog Regression Downscaling (GARD) algorithm to downscale meteorological hindcasts from the NASA Goddard Earth Observing System (GEOS) global S2S forecast system. Downscaled meteorological fields are then applied to drive offline simulations with the Catchment Land Surface Model (CLSM) to forecast United States Drought Monitor (USDM) style drought indicators derived from simulated surface hydrology variables. We compare the representation of drought in these downscaled hindcasts to hindcasts that are not downscaled, using the North American Land Data Assimilation System Phase 2 (NLDAS-2) dataset as an observational reference. We find that downscaling using GARD improves hindcasts of temperature and temperature anomalies, but the results for precipitation are mixed and generally small. Overall, GARD downscaling led to improved hindcast skill for total drought across the Contiguous United States (CONUS), and improvements were greatest for extreme (D3) and exceptional (D4) drought categories.

Ryan A Zamora↗

Effects of Subsurface Sampling & Processing on Martian Simulant Containing Varying Quantities of Water

The presence of water-ice in the Martian subsurface is a subject of much debate and excited speculation. Recent results from the gammaray spectrometer (GRS) on board NASA's Mars Odyssey spacecraft indicate the presence of large amounts of hydrogen in regions of predicted ice stability. The combination of chemistry, low gravitational field (3.71 m/s(exp 2)) and a surface pressure of about 6.36 mbar at the mean radius, place limits on the stability of H2O on the surface, however, results from the GRS indicate that the hydrogen rich phase may be present at a depth as shallow as one meter in some locations on Mars. The potential for water on Mars leads directly to the speculation that life may once have existed there, since liquid water is the unifying factor for environments known to support life on Earth. Lubricant-free drilling has been considered as a means of obtaining water-rich subsurface samples on Mars, and two recent white papers sponsored by the Mars Program have attempted to identify the problems associated with this goal. The two major issues identified were: the engineering challenges of drilling into a water-soil mixture where phase changes may occur; and the potential to compromise the integrity of in-situ scientific analysis due to contamination, volatilization, and mineralogical or chemical changes as a result of processing. This study is a first attempt to simulate lubricantfree drilling into JSC Mars-1 simulant containing up to 50% water by weight. The goal is to address the following: 1) Does sample processing cause reactions or changes in mineralogy which will compromise the interpretation of scientific measurements conducted on the surface? 2) Does the presence of water-ice in the sample complicate (1)? 3) Do lubricant-free drilling and processing leave trace contaminants which may compromise our understanding of sample composition? 4) How does the torque/power required for drilling change as a function of water content and does this lead to unexpected thermal effects?

Menard, J.↗

Magnetic fine structure of solar coronal loops

A numerical simulation of the effect of a random photospheric flow on the magnetic structure of a coronal loop is presented. An initially uniform field embedded in a perfectly conducting plasma is assumed, extending between two flat parallel plates representing the solar photosphere at both ends of the loop. The field is perturbed by a sequence of randomly phased, sinusoidal flow patterns applied at one of the boundary plates, and the corresponding sequence of force-free fields is determined. It is found that the electric currents generated by these flows develop a fine structure on a scale significantly smaller than the wavelength of the velocity patterns. This suggests that magnetic energy is transferred to smaller scale via a cascade process.

Van Ballegooijen, A. A.↗

Radiative transfer theory for polarimetric remote sensing of pine forest

The radiative transfer theory is applied to interpret polarimetric radar backscatter from pine forest with clustered vegetation structures. The scattering function of each cluster is calculated by incorporating the phase interference of scattered fields from each component. The resulting phase matrix is used in the radiative transfer equations to evaluate the polarimetric backscattering coefficients from random medium layers embedded with vegetation clusters. Upon including multiscale structures (trunks, primary and secondary branches, and needles), polarimetric radar responses from pine forest for different frequencies and looking angles are interpreted and simulated. Preliminary results are shown to be in good agreement with the measured backscattering coefficients at the Landes maritime pine forest during the MAESTRO-1 experiment.

Hsu, C. C.↗

CFD simulation of coaxial injectors

The development of improved performance models for the Space Shuttle Main Engine (SSME) is an important, ongoing program at NASA MSFC. These models allow prediction of overall system performance, as well as analysis of run-time anomalies which might adversely affect engine performance or safety. Due to the complexity of the flow fields associated with the SSME, NASA has increasingly turned to Computational Fluid Dynamics (CFD) techniques as modeling tools. An important component of the SSME system is the fuel preburner, which consists of a cylindrical chamber with a plate containing 264 coaxial injector elements at one end. A fuel rich mixture of gaseous hydrogen and liquid oxygen is injected and combusted in the chamber. This process preheats the hydrogen fuel before it enters the main combustion chamber, powers the hydrogen turbo-pump, and provides a heat dump for nozzle cooling. Issues of interest include the temperature and pressure fields at the turbine inlet and the thermal compatibility between the preburner chamber and injector plate. Performance anomalies can occur due to incomplete combustion, blocked injector ports, etc. The performance model should include the capability to simulate the effects of these anomalies. The current approach to the numerical simulation of the SSME fuel preburner flow field is to use a global model based on the MSFC sponsored FNDS code. This code does not have the capabilities of modeling several aspects of the problem such as detailed modeling of the coaxial injectors. Therefore, an effort has been initiated to develop a detailed simulation of the preburner coaxial injectors and provide gas phase boundary conditions just downstream of the injector face as input to the FDNS code. This simulation should include three-dimensional geometric effects such as proximity of injectors to baffles and chamber walls and interaction between injectors. This report describes an investigation into the numerical simulation of GH2/LOX coaxial injectors. The following sections will discuss the physical aspects of injectors, the CFD code employed, and preliminary results of a simulation of a single coaxial injector for which experimental data is available. It is hoped that this work will lay the foundation for the development of a unique and useful tool to support the SSME program.

Landrum, D. Brian↗

Simulated rarefied entry of the Galileo probe into the atmosphere of Jupiter

Flow properties and aerodynamics are computed with a direct simulation Monte Carlo (DSMC) method for rarefied entry of the Galileo Probe into the atmosphere of Jupiter. Accurate predictions of vehicle drag coefficients are needed in order to assess atmospheric properties from the onboard Atmosphere Structure Experiment where highly-sensitive accelerometers will measure the drag force to within 10-6 barr during the initial entry phase at high altitudes. The corresponding flow rarefraction extends from the free molecule limit to the near continuum transition regime (Re less than 1000). Simulation results indicate that C(sub D) varies from 2.1 at the free molecule limit down to 1.6 at Re(infinity) = 1,000. Temperatures, densities, and internal energies throughout the flow field were also computed at each altitude ranging from 735 km to 353 km above the 1 barr level in the Jovian atmosphere. Surface heating and temperatures of the probe were computed directly in the DSMC code by assuming radiative equilibrium. Material response was re-asssessed accurately during entry by accounting for conductivity, heat capacity, and pyrolysis which led to surface material mass efflux several times that of the freestream mass influx. The simulation also accounted for the quantum nature of the rotational energy mode of the dominant atmospheric species H2 through partial internal excitation in the freestream gas.

Hass, Brian L.↗

Global-Scale Consequences of Magnetic-Helicity Injection and Condensation on the Sun

In the recent paper of Antiochos, a new concept for the injection of magnetic helicity into the solar corona by small-scale convective motions and its condensation onto polarity inversion lines (PILs) has been developed. We investigate this concept through global simulations of the Sun's photospheric and coronal magnetic fields and compare the results with the hemispheric pattern of solar filaments. Assuming that the vorticity of the cells is predominately counter-clockwise/clockwise in the northern/southern hemisphere, the convective motions inject negative/positive helicity into each hemisphere. The simulations show that: (i) On a north-south orientated PIL, both differential rotation and convective motions inject the same sign of helicity which matches that required to reproduce the hemispheric pattern of filaments. (ii) On a high latitude east-west orientated polar crown or sub-polar crown PIL, the vorticity of the cells has to be approximately 2-3 times greater than the local differential rotation gradient in order to overcome the incorrect sign of helicity injection from differential rotation. (iii) In the declining phase of the cycle, as a bipole interacts with the polar field, in some cases helicity condensation can reverse the effect of differential rotation along the East-West lead arm, but not in all cases. The results show that this newly developed concept of magnetic helicity injection and condensation is a viable method to explain the hemispheric pattern of filaments in conjunction with the mechanisms used in Yeates et al. (2008). Future observational studies should focus on determining the vorticity component within convective motions to determine, both its magnitude and latitudinal variation relative to the differential rotation gradient on the Sun.

magnetic fields↗

Chromospheric-coronal coupling during solar flares: Current systems and particle acceleration

Two-dimensional (three velocity) electrostatic particle simulations are used to investigate the particle heating and acceleration associated with the impulsive phase of a solar flare. A crossfield current in the high corona (which is presumably driven by reconnection processes) is used to initiate the flare. Due to the differential motion of the electrons and ions, currents, and associated quasi-static electric fields are generated with the primary current and balancing return current being on adjacent field lines. These currents extend from the corona down into the chromosphere. Electrons can be accelerated to energies exceeding 100 keV on short time scales via the quasi-static fields and wave-particle interactions. The spectra of these electrons has a broken power-law distribution which hardens in time. The spatially separate primary and return currents are closed by the cross-field acceleration of the ambient ions into the primary current regions. These ions are then accelerated upwards into the corona by the same quasi-static electric field accelerating the electrons downwards. This acceleration can account for the broadened stationary and weak blue shifted component seen in soft x ray line emissions and enhancements in heavy ion abundances seen in the solar wind in associations with solar flares.

Winglee, Robert M.↗

The Diurnal Cycle in TOGA-COARE: Regional Scale Model Simulations

The diurnal variation of precipitation processes over the tropics is a well-known phenomenon and has been studied using surface rainfall data, radar reflectivity data, and satellite-derived cloudiness and precipitation. Recently, Sui (1997) analyzed observations from TOGA COARE in the tropical western Pacific ocean to study the relevant mechanisms producing diurnal variation of precipitation. They found that the diurnal SST cycle is important for afternoon showers in the undisturbed periods and diurnal radiative processes for nocturnal rainfall. Takayabu (1996) found a quasi-2-day cycle in precipitation during TOGA COARE and they suggested that inertia-gravity waves may be associated with this 2-day cycle. Chen and Houze (1997), however, suggested that the quasi-2-day oscillation is mainly a function of the time required by the lower-tropospheric moisture field to recover from the drying caused by deep convection. In this study, the Penn State/NCAR Mesoscale Model (MM5) with improved physics (i.e., cloud microphysics, radiation, land-soil-vegetation-surface processes, and TOGA COARE flux scheme) and a multiple level nesting technique will be used to simulate two TOGA COARE convective espisodes, one convectively suppressed phase (mid to late January 1993) and one convectively active phase (mid to late December 1992). We will examine precipitation processes over the open ocean and over land by comparing MM5 simulated rainfall and out-going longwave radiation. By examining the OLR and precipitation, we can determine if there is a temporal lag between the maximum precipitation and the coldest upwelling longwave radiation (the time-lag between stratiform-cirrus and convective towers). The boundary layer response by (or recovery from) precipitation processes will also be shown by examining the PBL thermodynamic structure and the sensible and latent heat fluxes. A preliminary MM5 simulation showed a clear diurnal variation in rainfall over both land and open ocean for the convectively active phase. The results also indicated that a quasi-2-day cycle in precipitation was simulated only over part of the ocean.

Tao, W.-K.↗

Contribution to Solar Brightness of Small-Size Magnetic Elements

Irradiance variability is mostly driven by surface magnetism, each magnetic feature contributing in a different manner. The contribution of small-size magnetic elements observed ubiquitously at high-spatial resolution on the solar photosphere, is still debated, as such features are mostly unresolved on full-disk images employed to model irradiance variability. Understanding the contribution to solar brightness of small-size magnetic elements, especially in quiet regions, will help understanding irradiance variability on the decadal and longer temporal scales, which, in turn, are fundamental to understand the role of the Sun on the Earth’s climate. We present a preliminary study of the brightness of small-size magnetic elements using high-spatial resolution spectro-polarimetric observations acquired with the National Science Foundations Daniel K. Inouye Solar Telescope (DKIST) during Commissioning Phase 1. We focus on the Intensity contrast vs magnetic field relation, which is a fundamental observable in irradiance studies. Previous studies mostly focused on the FeI 630.1/630.2 nm spectral range and so to contextualize our results we focus here on this spectral region as well. By comparing DKIST results with MURaM magneto-hydrodynamic simulations, and with previous results published in the literature, we conclude that our results are in agreement with studies conducted on HINODE observations. We plan to extend the analysis to UV and IR ViSP observations and VBI blue continuum in order to better constrain the quiet sun contribution to spectral irradiance variations.

SMD↗

Modeling and numerical simulation of microwave pulse propagation in air breakdown environment

It is shown that ionization occurs wherever the field intensity of the pulse exceeds the local breakdown threshold field of the background air. The produced plasma then attenuates the pulse and gives rise to a tail erosion phenomenon that plays the primary role in limiting the energy transfer of the pulse from source to destination. A theoretical model describing the propagation of an intense microwave pulse in an air breakdown environment is developed that includes the possible focusing effect introduced by either using phase array antennas or the other arrangements. The self-consistent description of the propagation process is provided by a set of two modal equations. These include a continuity equation (Poynting's equation) for the energy density of the pulse and a rate equation of the electron density. A forward wave approximation is used to simplify Poynting's equation, and a semiempirical formula is used for the ionization frequency, nu(sub i). This frequency provides the coupling between the two modal equations, and is used to express the electron rate equation explicitly. In terms of the relevant parameters of the atmosphere, these two equations are normalized for numerical analysis of pulse propagation in the atmosphere. The dependencies of the propagation characteristics of the pulse on intensity, frequency, width, and shape of the pulse are determined. The numerical simulations lead to a useful empirical relation p(exp 3)w = alpha = constant, where p and w are the incident power and width of the pulse and alpha depends on the percentage of pulse energy transferred from the source point to a destined position. The density distribution of the pulse's self-generated plasma is also evaluated. The results also show that for ionization caused by a single unfocused microwave pulse transmitted upwards from the ground, the maximum electron density produced at, for example, 50 km altitude is limited by the tail erosion effect to below 10(exp 6) cm(exp -3). Repetitive pulse and focused beam approaches are also examined. Both approaches can increase the maximum electron density by no more than an order of magnitude. A scheme using two obliquely propagating pulses intersecting at the destined height, e.g. 50 km, is considered. It is shown that the electron density generated at the lowest intersecting position can easily reach a value of 6.6 x 10(exp 8) cm(exp -3), which is considered to be high enough for artificial ionospheric mirror (AIM) application.

Kim, J.↗

CFD Modeling & Validation of Single Bubble Growth During Cavitation in Zero Gravity

Cryogenic Fluid Management (CFM) in zero gravity represents a vital area of research as NASA pursues future long-duration space missions. Success of these missions relies heavily on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. In the Zero Boil-Off Tank (ZBOT-1) Experiment in 2018, unexpected and intense phase change occurred at the screened Liquid Acquisition Device (LAD) during tank pressure control that was later attributed to cavitation. The cavitation ultimately resulted in pump failure. Since cryogen cavitation during storage and transfer of propellants in microgravity has been a source of serious concern for the aerospace community, a series of drop tower experiments with liquid methane were also performed at the University of Bremen to carefully study cavitation from a single bubble nucleus in an engineered micron-sized wall cavity. As the test cell is depressurized during the 10-second drop, cavitation causes bubble growth from the nucleus. Here, in a collaboration between the Case Western Reserve University and University of Bremen, we present a two-phase Computational Fluid Dynamics (CFD) model of bubble growth from a single wall cavity during depressurization in microgravity. Both 2D axisymmetric and 3D parametric simulation case studies are performed to show the evolution of the bubble growth. The model uses a Volume of Fluid (VOF) approach with and without adaptive mesh refinement and adaptive time stepping to capture the growth dynamics using an explicit formulation with algebraic reconstruction of the interface. CFD results show good agreement with corresponding experiment. Based on the velocity, temperature, and volume fraction fields generated by the validation case studies, detailed discussion and analysis of fluid flow and heat transfer during bubble growth process are presented to delineate the intricacies of the cavitation process. Strengths and shortcomings of the present numerical approach in capturing this phenomenon are also highlighted.

Two Phase Flow↗

Variability of the Reconnection Guide Field in Solar Flares

Solar flares may be the best-known examples of the explosive conversion of magnetic energy into bulk motion, plasma heating, and particle acceleration via magnetic re-connection. The energy source for all flares is the highly sheared magnetic field of a filament channel above a polarity inversion line (PIL). During the flare, this shear field becomes the so-called reconnection guide field (i.e., the non-reconnecting component), which has been shown to play a major role in determining key properties of the recon-nection including the efficiency of particle acceleration. We present new high-resolution, three-dimensional, magnetohydrodynamics simulations that reveal the detailed evolution of the magnetic shear/guide field throughout an eruptive flare. The magnetic shear evolves in three distinct phases: shear first builds up in a narrow region about the PIL, then expands outward to form a thin vertical current sheet, and finally is transferred by flare reconnection into an arcade of sheared flare loops and an erupting flux rope. We demonstrate how the guide field may be inferred from observations of the sheared flare loops. Our results indicate that initially the guide field is larger by about a factor of 5 than the reconnecting component, but it weakens by more than an order of magnitude over the course of the flare. Instantaneously, the guide field also varies spatially over a similar range along the three-dimensional current sheet. We discuss the implications of the remarkable variability of the guide field for the timing and localization of efficient particle acceleration in flares.

Solar magnetic reconnection↗

Immiscible materials and alloys

Only immiscible materials which involve a fluid phase have a potential advantage when considered for space processing. Processing of such systems at low gravity can lead to a material which is very homogeneous and contains a finely dispersed mixture of phases. Materials with such structures may exhibit potentially useful properties for such applications as superconductors, dispersion-strengthened materials, superplastic materials, permanent magnets, etc. These conclusions have resulted from review of past efforts in the field of low-g processing of materials containing a liquid phase miscibility gap and from present efforts dealing with such systems. These studies have concentrated on precipitation of liquid droplets in a host liquid during cooling through the miscibility gap. This effort has been to study the agglomeration of the droplets both experimentally and through computer simulation in order to understand the mechansism involved and to deduce the effect of gravity on the droplet distribution and resulting structure of the solidified material. The agglomeration mechanisms treated in detail are those due to diffusional growth and collision processes caused by Stokes on convection current-induced droplet migration.

Markworth, A. J.↗

Noncircular Cross Sections Could Enhance Mixing in Sprays

A computational study has shown that by injecting drops in jets of gas having square, elliptical, triangular, or other noncircular injection cross sections, it should be possible to increase (relative to comparable situations having circular cross section) the entrainment and dispersion of liquid drops. This finding has practical significance for a variety of applications in which it is desirable to increase dispersion of drops. For example, in chemical-process sprays, increased dispersion leads to increases in chemical- reaction rates; in diesel engines, increasing the dispersion of drops of sprayed fuel reduces the production of soot; and in household and paint sprays, increasing the dispersion of drops makes it possible to cover larger surfaces. It has been known for some years that single-phase fluid jets that enter flow fields through noncircular inlets entrain more fluid than do comparable jets entering through circular inlets. The computational study reported here was directed in part toward determining whether and how this superior mixing characteristic of noncircular single phase jets translates to a similar benefit in cases of two-phase jets (that is, sprays). The study involved direct numerical simulations of single- and two-phase free jets with circular, elliptical, rectangular, square, and triangular inlet cross sections. The two-phase jets consisted of gas laden with liquid drops randomly injected at the inlets. To address the more interesting case of evaporating drops, the carrier gas in the jets was specified to be initially unvitiated by the vapor of the liquid chemical species and the initial temperature of the drops was chosen to be smaller than that of the gas. The mathematical model used in the study was constructed from the conservation equations for the two-phase flow and included complete couplings of mass, momentum, and energy based on thermodynamically self-consistent specification of the enthalpy, internal energy, and latent heat of vaporization of the vapor.

Bellan, Josette↗