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 289 records · Page 16

Energy buildup in coronal magnetic flux tubes

A time-dependent two-dimensional MHD simulation is used to study the response of the magnetic field in coronal loops to photospheric motion. From an initially uniform field, circular sections of the ends of the loop are slowly rotated to represent the photospheric motion. The evolution of the field and flow is characterized by three phases: (1) a phase of negligible kinetic energy where the current and field are predominantly parallel; (2) a phase where the field twist increases, the axial field at and near the axis increases, and the axial field decreases in two cylindrical regions away from the axis; and (3) a phase in which a significant portion of the field makes several rotations at large radii, with a corresponding reducton in the axial field to a few percent of the initial value.

Steinolfson, R. S.↗

Evidence of a Plasmoid-Looptop Interaction and Magnetic Inflows During a Solar Flare/Coronal Mass Ejection Eruptive Event

Observational evidence is presented for the merging of a downward-propagating plasmoid with a looptop kernel during an occulted limb event on 2007 January 25. RHESSI light curves in the 9-18 keV energy range, as well as that of the 245 MHz channel of the Learmonth Solar Observatory, show enhanced nonthermal emission in the corona at the time of the merging suggesting that additional particle acceleration took place. This was attributed to a secondary episode of reconnection in the current sheet that formed between the two merging sources. RHESSI images were used to establish a mean downward velocity of the plasmoid of 12 km/s. Complementary observations from the SECCHI suite of instruments on board STEREO-B showed that this process occurred during the acceleration phase of the associated coronal mass ejection (CME). From wavelet-enhanced EUV Imager, image evidence of inflowing magnetic field lines prior to the CME eruption is also presented. The derived inflow velocity was found to be 1.5 km/s. This combination of observations supports a recent numerical simulation of plasmoid formation, propagation, and subsequent particle acceleration due to the tearing mode instability during current sheet formation.

Milligan, Ryan O.↗

Global Images of Trapped Ring Current Ions During Main Phase of 17 March 2015 Geomagnetic Storm as Observed by TWINS

A unique view of the trapped particles in the inner magnetosphere provided by energetic neutral atom (ENA) imaging is used to observe the dynamics of the spatial structure and the pitch angle anisotropy on a global scale during the last 6 h of the main phase of a large geomagnetic storm (minimum SYM-H 230 nT) that began on 17 March 2015. Ion flux and pressure anisotropy obtained from Two Wide-angle Imaging Neutral-atom Spectrometers (TWINS) ENA images are shown. The ion flux shows two peaks, an inner one at approximately radii 34 RE in the dusk-to-midnight sector and an outer peak at radii 89 RE prior to midnight. The inner peak is relatively stationary during the entire period with some intensification during the final steep decline in SYM-H to its minimum. The outer peak shows the significant temporal variation brightening and dimming and finally disappearing at the end of the main phase. The pressure anisotropy shows the expected perpendicular pitch angles inside of L 6 but shows parallel pitch angles at greater L values. This is interpreted as consistent with pitch angle-dependent drift as modeled in the Tsy05 magnetic field and Comprehensive Inner Magnetosphere-Ionosphere simulations. The TWINS results are compared directly with Radiation Belt Storm Probes Ion Composition Experiment (RBSPICE)-A measurements. Using 15 min snapshots of flux and pressure anisotropy from TWINS along the path of RBSPICE-A during the 6 h focused upon in this study, the essential features displayed in the TWINS global images are supported.

Perez, J. D.↗

A quantitative approach to the topology of large-scale structure

A quantitative measure of the topology of large-scale structure: the genus of density contours in a smoothed density distribution, is described and applied. For random phase (Gaussian) density fields, the mean genus per unit volume exhibits a universal dependence on threshold density, with a normalizing factor that can be calculated from the power spectrum. If large-scale structure formed from the gravitational instability of small-amplitude density fluctuations, the topology observed today on suitable scales should follow the topology in the initial conditions. The technique is illustrated by applying it to simulations of galaxy clustering in a flat universe dominated by cold dark matter. The technique is also applied to a volume-limited sample of the CfA redshift survey and to a model in which galaxies reside on the surfaces of polyhedral 'bubbles'. The topology of the evolved mass distribution and 'biased' galaxy distribution in the cold dark matter models closely matches the topology of the density fluctuations in the initial conditions. The topology of the observational sample is consistent with the random phase, cold dark matter model.

Gott, J. Richard, III↗

Unmanned Aircraft Systems (UAS) Integration in the National Airspace System (NAS) Project

Phase 1 of the UAS-NAS (Unmanned Aircraft Systems-National Airspace System) project focused on MOPS (Minimum Operational Performance Standards) development for large UAS transitioning through Class D, E, and G airspace. Phase 2 activities are currently focused on extended UAS operations in Class D, E, and G airspace, as well as sensors and architectures that will enable DAA (Detect and Avoid) equipment to be installed on a wider range of UAS. Encounters with non-cooperative intruders in low altitude airspace under 10,000 feet will be explored for UAS with a low Size, Weight, and Power (SWaP) radar. New sensors for detecting non-cooperative intruder aircraft will have a more limited detection range and field of regard; therefore, the analysis of non-cooperative encounter geometries may help in developing the requirements for an onboard low SWaP radar (and electro-optical sensor). The current simulation, led by HSI (Human Systems Integration Division at NASA Ames Research Center) as part of the DAA subproject, will support efforts toward developing a modified DAA Well Clear (DWC) definition that would be more appropriate for UAS equipped with limited surveillance and aircraft performance capabilities compared to the Phase 1 DAA system. The findings will inform Phase 2 display requirements for alternative surveillance use cases.

Collision Avoidance↗

Perpendicular heating of electrons by upper hybrid waves generated by a ring distribution

Satellite observations of electron conical distributions with enhanced fluxes just outside the loss cone suggest that telectrons have been heated perpendicularly to the magnetic field in the mid-altitude polar magnetosphere. To understand electron conical distributions, plasma simulations are conducted to examine an upper hybrid wave instability of a ring electron distribution perpendicular to the magnetic field in a cold electron background. The simulations indicate that both the cold and ring distributions are heated perpendicularly during the saturation stage. From the plasma data, a ring distribution can be identified as a trapped distribution function with an enhancement near 90-deg pitch angle in the phase space density plot. It is suggested that the ring distribution might provide an additional free energy source for generating upper hybrid waves associated with electron conical events.

Lin, C. S.↗

Gamma radiation from pulsar magnetospheric gaps

We investigate the production of gamma rays in two pulsar emission models: the 'polar cap' model and the 'outer cap' model. For the former, we have performed detailed simulations of energetic electrons flowing in the vacuum dipole open field line region. In the outer gap case, we generate light curves for various magnetosphere geometries. Using data from radio and optical observations, we construct models for specific viewing angles appropriate to the Crab and Vela pulsars. Phase-resolved spectra are also computed in the polar cap case and provide signatures for testing the models. The calculations have been extended to include millisecond pulsars, and we have been able to predict fluxes and spectra for populations of recycled pulsars, which are compared to COS B data for globular cluster populations.

Chiang, James↗

Anisotropic turbulence in the solar wind

Solar wind turbulence has been viewed traditionally as composed of parallel propagating ('slab' fluctuations) or otherwise as isotropic turbulence. A variety of recent investigations, reviewed here, indicate that the spectrum may contain a significant admixture of two dimensional fluctuations, having variations mainly perpendicular to the local magnetic field. These indications come from simulations, from the theory of nearly incompressible MHD, from cosmic ray transport studies and from transport theory for solar wind turbulence, as well as from interpretations of direct observations. Thus, solar wind turbulence may be more like bundles of spaghetti than like parallel phase fronts.

Matthaeus, W. H.↗

The ZBOT-NC Experiment - Effects of Non-Condensable Gases on Propellant Tank Pressurization and Pressure Control

Introduction Integral to all phases of NASA’s projected planetary expeditions is affordable and reliable cryogenic fluid storage for use in propellant or life support systems. It is greatly advantageous to develop innovative vent-less pressure control designs based on cooling/mixing of the bulk tank fluid to allow storage of the cryogenic fluid with zero or reduced boil-off. The presence of noncondensable gases, can interfere with the condensation at the interface impacting tank pressure control during subcooled jet mixing especially in microgravity. The Zero-Boil-Off Tank (ZBOT) Experiments are a series of small-scale experiments aboard the International Space Station (ISS) that use a transparent volatile simulant fluid in a transparent sealed tank to delineate various fundamental fluid flow, heat and mass transport, and phase change phenomena associated with storage tank pressurization and pressure control in microgravity. The ZBOT-1 experiment was performed on the ISS 2017-2018 timeframe and collected data to validate a state-of-the-art CFD model for tank pressurization and pressure control for a pure system. The ZBOT-NC Experiment is the second experiment in the series to be performed on the International Space Station (ISS) in 2025. Its goal is to investigate the effects of noncondensable gases on interfacial evaporation and condensation during self -pressurization and jet-mixing pressure control in microgravity for a two-component system. Materials & Methods In this work, we will describe the detailed features of the ZBOT-NC experimental hardware and Diagnostics that includes nonintrusive Quantum Dot Thermometry (QDT) for whole field temperature measurement. All microgravity pressurization and pressure control tests will be performed for both the pure and the two-components systems with Xenon and Neon as the two noncondensable gases spanning small and large molecular weights and sizes. The two-phase CFD model that is developed as part of the project will be also presented and discussed. Results Ground-based pressurization and jet mixing experiments and CFD simulation results are compared to each other to validate both the fidelity of the CFD model predictions and the accuracy of the QDT measurement. Model simulations for noncondensable gas effects will also be compared against large Cryogenic LH2-GHe experiments to indicate the noncondensable gas effects on tank pressure control during jet mixing in 1G. Finally, CFD results will be presented to predict the effects of the noncondensable gas during subcooled jet mixing during the ZBOT-NC in advance of the microgravity experiment in 2015.

Evaporation Condensation↗

The Zero-Boil-Off Tank Experiment: the Effects of Non-Condensable Gases on Pressurization and Pressure Control of Propellant Tanks in Microgravity

Introduction Integral to all phases of NASA’s projected planetary expeditions is affordable and reliable cryogenic fluid storage for use in propellant or life support systems. It is greatly advantageous to develop innovative vent-less pressure control designs based on cooling/mixing of the bulk tank fluid to allow storage of the cryogenic fluid with zero or reduced boil-off. The presence of noncondensable gases, can interfere with the condensation at the interface impacting tank pressure control during subcooled jet mixing especially in microgravity. The Zero-Boil-Off Tank (ZBOT) Experiments are a series of small-scale experiments aboard the International Space Station (ISS) that use a transparent volatile simulant fluid in a transparent sealed tank to delineate various fundamental fluid flow, heat and mass transport, and phase change phenomena associated with storage tank pressurization and pressure control in microgravity. The ZBOT-1 experiment was performed on the ISS 2017-2018 timeframe and collected data to validate a state-of-the-art CFD model for tank pressurization and pressure control for a pure system. The ZBOT-NC Experiment is the second experiment in the series to be performed on the International Space Station (ISS) in 2025. Its goal is to investigate the effects of noncondensable gases on interfacial evaporation and condensation during self -pressurization and jet-mixing pressure control in microgravity for a two-component system. Materials & Methods In this work, we will describe the detailed features of the ZBOT-NC experimental hardware and Diagnostics that includes nonintrusive Quantum Dot Thermometry (QDT) for whole field temperature measurement. All microgravity pressurization and pressure control tests will be performed for both the pure and the two-components systems with Xenon and Neon as the two noncondensable gases spanning small and large molecular weights and sizes. The two-phase CFD model that is developed as part of the project will be also presented and discussed. Results Ground-based pressurization and jet mixing experiments and CFD simulation results are compared to each other to validate both the fidelity of the CFD model predictions and the accuracy of the QDT measurement. Model simulations for noncondensable gas effects will also be compared against large Cryogenic LH2-GHe experiments to indicate the noncondensable gas effects on tank pressure control during jet mixing in 1G. Finally, CFD results will be presented to predict the effects of the noncondensable gas during subcooled jet mixing during the ZBOT-NC in advance of the microgravity experiment in 2015.

Evaporation Condensation↗

Use of Polarimetric Radar Measurements to Constrain Simulated Convective Cell Evolution: A Pilot Study with Lagrangian Tracking

To probe the potential value of a radar-driven field campaign to constrain simulation of isolated convection subject to a strong aerosol perturbation, convective cells observed by the operational KHGX weather radar in the vicinity of Houston, Texas, are examined individually and statistically. Cells observed in a single case study of onshore flow conditions during July 2013 are first examined and compared with cells in a regional model simulation. Observed and simulated cells are objectively identified and tracked from observed or calculated positive specific differential phase (K(sub DP)) above the melting level, which is related to the presence of supercooled liquid water. Several observed and simulated cells are subjectively selected for further examination. Below the melting level, we compare sequential cross sections of retrieved and simulated raindrop size distribution parameters. Above the melting level, we examine time series of KDP and radar differential reflectivity (Z(sub DR)) statistics from observations and calculated from simulated supercooled rain properties, alongside simulated vertical wind and supercooled rain mixing ratio statistics. Results indicate that the operational weather radar measurements offer multiple constraints on the properties of simulated convective cells, with substantial value added from derived K(sub DP) and retrieved rain properties. The value of collocated three-dimensional lightning mapping array measurements, which are relatively rare in the continental US, supports the choice of Houston as a suitable location for future field studies to improve the simulation and understanding of convective updraft physics. However, rapid evolution of cells between routine volume scans motivates consideration of adaptive scan strategies or radar imaging technologies to amend operational weather radar capabilities. A 3-year climatology of isolated cell tracks, prepared using a more efficient algorithm, yields additional relevant information. Isolated cells are found within the KHGX domain on roughly 40 % of days year-round, with greatest concentration in the northwest quadrant, but roughly 5-fold more cells occur during June through September. During this enhanced occurrence period, the cells initiate following a strong diurnal cycle that peaks in the early afternoon, typically follow a south-to-north flow, and dissipate within 1 h, consistent with the case study examples. Statistics indicate that ∼ 150 isolated cells initiate and dissipate within 70 km of the KHGX radar during the enhanced occurrence period annually, and roughly 10 times as many within 200 km, suitable for multi-instrument Lagrangian observation strategies. In addition to ancillary meteorological and aerosol measurements, robust vertical wind speed retrievals would add substantial value to a radar-driven field campaign.

Fridlind, Ann M.↗

SPS phase control system performance via analytical simulation

A solar power satellite transmission system which incorporates automatic beam forming, steering, and phase control is discussed. The phase control concept centers around the notation of an active retrodirective phased array as a means of pointing the beam to the appropriate spot on Earth. The transmitting antenna (spacetenna) directs the high power beam so that it focuses on the ground-based receiving antenna (rectenna). A combination of analysis and computerized simulation was conducted to determine the far field performance of the reference distribution system, and the beam forming and microwave power generating systems.

Lindsey, W. C.↗

Development of an Atomization Methodology for Spray Combustion

In liquid rocket propulsion, the knowledge and the understanding of liquid-gas interfacial phenomena are very important. This is important for predicting the onset of cavitation occurring in swirl injection elements used in STME, as well as atomization processes in shear-induced injectors (co-axial) and impinging injector elements. From the fact that all the physical processes including droplet size distribution, droplet dispersion, mixing and combustion are controlled by atomization processes, it is expected that the successful incorporation of the volume of fraction (VOF) will greatly enhance the analytical capability of predicting spray combustion processes in liquid-fueled engines. In this paper, a methodology is developed to define and track interfaces between two fluids in non-orthogonal, body-fitted grids using a single fractional volume of fluid (VOF) variable to describe the distribution of the liquid phase in a gas-liquid flow field. This method was implemented in a mature CFD code MAST (Multiphase All-Speed Transient) utilizing the general PISO-C algorithm. For the preliminary study on the analysis of spray combustion and tracking of the interface between two phases, we will report on the progress of the simulation of the instability on the liquid column; the surface wave instability and the droplet breakup from the liquid surface.

Seung, S. P.↗

Computational Modeling of Meteor-Generated Ground Pressure Signatures

We present a thorough validation of a computational approach to predict infrasonic signatures of centimeter-sized meteoroids. We assume that the energy deposition along the meteor trail is dominated by atmospheric drag and simulate the steady, inviscid flow of air in thermochemical equilibrium to compute the meteoroid's near-body pressure signature. This signature is then propagated through a stratified and windy atmosphere to the ground using a methodology adapted from aircraft sonic-boom analysis. An assessment of the numerical accuracy of the near field and the far field solver is presented. The results show that when the source of the signature is the cylindrical Mach-cone, the simulations closely match the observations. The prediction of the shock rise-time, the zero-peak amplitude of the waveform, and the duration of the positive pressure phase are consistently within 10% of the measurements. Uncertainty in the shape of the meteoroid results in a poorer prediction of the trailing part of the waveform. Overall, our results independently verify energy deposition estimates deduced from optical observations.

Computational↗

Stormtime transport of ring current and radiation belt ions

This is an investigation of stormtime particle transport that leads to formation of the ring current. Our method is to trace the guiding-center motion of representative ions (having selected first adiabatic invariants mu) in response to model substorm-associated impulses in the convection electric field. We compare our simulation results qualitatively with existing analytically tractable idealizations of particle transport (direct convective access and radial diffusion) in order to assess the limits of validity of these approximations. For mu approximately less than 10 MeV/G (E approximately less than 10 keV at L equivalent to 3) the ion drift period on the final (ring-current) drift shell of interest (L equivalent to 3) exceeds the duration of the main phase of our model storm, and we find that the transport of ions to this drift shell is appropriately idealized as direct convective access, typically from open drift paths. Ion transport to a final closed drift path from an open (plasma-sheet) drift trajectory is possible for those portions of that drift path that lie outside the mean stormtime separatrix between closed and open drift trajectories, For mu approximately 10-25 MeV/G (110 keV approximately less than E approximately less than 280 keV at L equivalent to 3) the drift period at L equivalent to 3 is comparable to the postulated 3-hr duration of the storm, and the mode of transport is transitional between direct convective access and transport that resembles radial diffusion. (This particle population is transitional between the ring current and radiation belt). For mu approximately greater than 25 MeV/G (radiation-belt ions having E approximately greater than 280 keV at L equivalent to 3) the ion drift period is considerably shorter than the main phase of a typical storm, and ions gain access to the ring-current region essentially via radial diffusion. By computing the mean and mean-square cumulative changes in 1/L among (in this case) 12 representative ions equally spaced in drift time around the steady-state drift shell of interest (L equivalent to 3), we have estimated (from both our forward and our time-reversed simulations) the time-integrated radial-diffusion coefficients D(sup sim)(sub LL) for particles having selected values of mu approximately greater than 15 MeV/G. The results agree surprisingly well with the predictions (D(sup ql)(sub LL)) of quasilinear radial diffusion theory, despite the rather brief duration (approximately 3 hrs) of our model storm and despite the extreme variability (with frequency) of the spectral-density function that characterizes the applied electric field during our model storm. As expected, the values of D(sup sim)(sub LL) deduced (respectively) from our forward and time-reversed simulations agree even better with each other and with D(sup sim)(sub LL) when the impulse amplitudes which characterize the individual substorms of our model storm are systematically reduced.

Chen, Margaret W.↗

Acoustic Optimization for Anti-Phase Asymmetric Rotor

This investigation seeks to optimize the implementation of anti-phase alternating trailing edge (TE) patterns for rotor noise suppression. The design objective is to maximize reduction of noise perceived by the community while maintaining the aerodynamic thrust. Computations using a three-dimensional Unsteady-Reynolds-Averaged-Navier-Stokes (URANS) with k-w Shear Stress Transport (SST) turbulence model and Ffowcs-Williams and Hawkings (FW-H) formula are used to obtain aerodynamic thrust and far-field noise level. A parametric acoustic study of 13 configurations of KDE rotor with variable alternating trailing edge period, alternating trailing edge length, and trailing edge deflection angle is conducted. The best design candidate for the KDE rotor has a four-period TE waveform which results in a reduction in far-field noise level of 2.1 dB in the hover condition and a reduction of 1.1 dB in the forward flight condition at 9.7 m/s. A further parametric acoustic study is conducted for a different rotor manufactured by APC. Six APC rotor design candidates are simulated. The best design candidate 4H for the APC rotor results in a reduction in far-field noise level of 4.0 dB in the hover condition and a reduction of 1.3 dB in the forward flight condition at 9.7 m/s. A series of acoustic experiments in the Penn State University (PSU) anechoic chamber have been conducted. In the forward flight condition at 9.7 m/s, the APC anti-phase 4H rotor offers clear evidence of noise suppression capability across a wide range of the azimuthal angle. In the broadband frequency range of 2000-4000 Hz, the APC anti-phase 4H rotor produces as much as 6 dB noise reduction. The experimental results appear to confirm the noise suppression capability of the proposed anti-phase rotor design concepts.

Xiong, Juntao↗

On Modeling Air/Space-Borne Radar Returns in the Melting Layer

The bright band is the enhanced radar echo associated with the melting of hydrometeors in stratiform rain where the melting process usually occurs below 0 C isotherm over a distance of about 500m. To simulate this radar signature, a scattering model of melting snow is proposed in which the fractional water content is prescribed as a function of the radius of a spherical mixed- phase particle consisting of air, ice and water. The model is based on the observation that melting starts at the surface of the particle and then gradually develops towards the center. To compute the scattering parameters of a non-uniform melting particle, the particle is modeled as a sphere represented by a collection of 64(exp 3) cubic cells of identical size where the probability of water at any cell is prescribed as a function of the radius. The internal field of the particle, used for deriving the effective dielectric constant, is computed by the Conjugate Gradient and Fast Fourier Transform (CGFFT) numerical methods. To make computations of the scattering parameters more efficient, a multi-layer stratified-sphere scattering model is introduced after demonstrating that the scattering parameters of the non-uniformly melting particle can be accurately reproduced by the stratified sphere. In conjunction with a melting layer model that describes the melting fractions and fall velocities of hydrometeors as a function of the distance from the 0 C isotherm, the stratified-sphere model is used to simulate the radar bright band profiles. These simulated profiles are shown to compare well with measurements from the Precipitation Radar (PR) aboard the Tropical Rainfall Measuring Mission (TRMM) satellite and a dual-wavelength airborne radar. The results suggest that the proposed model of a melting snow particle may be useful in studying the characteristics of the bright-band in particular and mixed- phase hydrometeors in general.

Liao, Liang↗

Computational Aerothermodynamics in Aeroassist Applications

Aeroassisted planetary entry uses atmospheric drag to decelerate spacecraft from super-orbital to orbital or suborbital velocities. Numerical simulation of flow fields surrounding these spacecraft during hypersonic atmospheric entry is required to define aerothermal loads. The severe compression in the shock layer in front of the vehicle and subsequent, rapid expansion into the wake are characterized by high temperature, thermo-chemical nonequilibrium processes. Implicit algorithms required for efficient, stable computation of the governing equations involving disparate time scales of convection, diffusion, chemical reactions, and thermal relaxation are discussed. Robust point-implicit strategies are utilized in the initialization phase; less robust but more efficient line-implicit strategies are applied in the endgame. Applications to ballutes (balloon-like decelerators) in the atmospheres of Venus, Mars, Titan, Saturn, and Neptune and a Mars Sample Return Orbiter (MSRO) are featured. Examples are discussed where time-accurate simulation is required to achieve a steady-state solution.

Gnoffo, Peter A.↗