Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multi-fluid model”

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 19 records

A multi-fluid model of an H2O-dominated dusty cometary atmosphere

The solution of the simultaneous set of differential equations representing conservation of number density, momentum, and energy, together with solar radiation transfer in streams that are responsible for major photolytic processes and nucleus heating, yields a self-consistent multifluid solution for the thermal and dynamical structure of an H2O-dominated, two-phase dusty gas cometary atmosphere, whose validity is restricted to the collision-dominated region where ions and neutrals are assumed to achieve a common temperature and velocity. The model is transonic, due to the dust in the inner coma, and involves the subsonic expansion of the ions and neutrals from the nucleus until, in the outermost part of the collision-dominated coma, temperatures of 300-400 K are reached due to UV pyrolytic heating.

Marconi, M. L.↗

Chemical and Hydrodynamical Models of Cometary Comae

Multi-fluid modelling of the outflowing gases which sublimate from cometary nuclei as they approach the Sun is necessary for understanding the important physical and chemical processes occurring in this complex plasma. Coma chemistry models can be employed to interpret observational data and to ultimately determine chemical composition and structure of the nuclear ices and dust. We describe a combined chemical and hydrodynamical model [1] in which differential equations for the chemical abundances and the energy balance are solved as a function of distance from the cometary nucleus. The presence of negative ions (anions) in cometary comae is known from Giotto mass spectrometry of 1P/Halley. The anions O(-), OH(-), C(-), CH(-) and CN(-) have been detected, as well as unidentified anions with masses 22-65 and 85-110 amu [2]. Organic molecular anions such as C4H(-) and C6H(-) are known to have a significant impact on the charge balance of interstellar clouds and circumstellar envelopes and have been shown to act as catalysts for the gas-phase synthesis of larger hydrocarbon molecules in the ISM, but their importance in cometary comae has not yet been fully explored. We present details of new models for the chemistry of cometary comae that include atomic and molecular anions and calculate the impact of these anions on the coma physics and chemistry af the coma.

Charnley, Steven↗

Development of a Computationally Efficient, High Fidelity, Finite Element Based Hall Thruster Model

This report documents the development of a two dimensional finite element based numerical model for efficient characterization of the Hall thruster plasma dynamics in the framework of multi-fluid model. Effect of the ionization and the recombination has been included in the present model. Based on the experimental data, a third order polynomial in electron temperature is used to calculate the ionization rate. The neutral dynamics is included only through the neutral continuity equation in the presence of a uniform neutral flow. The electrons are modeled as magnetized and hot, whereas ions are assumed magnetized and cold. The dynamics of Hall thruster is also investigated in the presence of plasma-wall interaction. The plasma-wall interaction is a function of wall potential, which in turn is determined by the secondary electron emission and sputtering yield. The effect of secondary electron emission and sputter yield has been considered simultaneously, Simulation results are interpreted in the light of experimental observations and available numerical solutions in the literature.

Jacobson, David↗

Propagation of Interplanetary Disturbances in the Outer Heliosphere

Contents include the following: 1. We have developed a one-dimensional, spherically symmetric, multi-fluid MHD model that includes solar wind protons and electrons, pickup ions, and interstellar neutral hydrogen. This model advances the existing solar wind models for the outer heliosphere in two important ways: one is that it distinguishes solar wind protons from pickup ions, and the other is that it allows for energy transfer from pickup ions to the solar wind protons. Model results compare favorably with the Voyager 2 observations. 2. 2. Solar wind slowdown and interstellar neutral density. The solar wind in the outer heliosphere is fundamentally different from that in the inner heliosphere since the effects of interstellar neutrals become significant. 3. ICME propagation from the inner to outer heliosphere. Large coronal mass ejections (CMEs) have major effects on the structure of the solar wind and the heliosphere. The plasma and magnetic field can be compressed ahead of interplanetary CMEs. 4. During the current solar cycle (Cycle 23), several major CMEs associated with solar flares produced large transient shocks which were observed by widely-separated spacecraft such as Wind at Earth and Voyager 2 beyond 60 AU. Using data from these spacecraft, we use the multi-fluid model to investigate shock propagation and interaction in the heliosphere. Specifically, we studied the Bastille Day 2000, April 2001 and Halloween 2003 events. 5. Statistical properties of the solar wind in the outer heliosphere. In a collaboration with L.F. Burlaga of GSFC, it is shown that the basic statistical properties of the solar wind in the outer heliosphere can be well produced by our model. We studied the large-scale heliospheric magnetic field strength fluctuations as a function of distance from the Sun during the declining phase of a solar cycle, using our numerical model with observations made at 1 AU during 1995 as input. 6. Radial heliospheric magnetic field events. The heliospheric magnetic field (HMF) direction, on average, conforms well to the Parker spiral.

Wang, Chi↗

Science Enhancements by the MAVEN Participating Scientists

NASA implemented a Participating Scientist Program and released a solicitation for the Mars Atmosphere and Volatile EvolutioN mission (MAVEN) proposals on February 14, 2013. After a NASA peer review panel evaluated the proposals, NASA Headquarters selected nine on June 12, 2013. The program's intent is to enhance the science return from the mission by including new investigations that broaden and/or complement the baseline investigations, while still addressing key science goals. The selections cover a broad range of science investigations. Included are: a patching of a 3D exosphere model to an improved global ionosphere-thermosphere model to study the generation of the exosphere and calculate the escape rates; the addition of a focused study of upper atmosphere variability and waves; improvement of a multi-fluid magnetohydrodynamic model that will be adjusted according to MAVEN observations to enhance the understanding of the solar-wind plasma interaction; a global study of the state of the ionosphere; folding MAVEN measurements into the Mars International Reference Ionosphere under development; quantification of atmospheric loss by pick-up using ion cyclotron wave observations; the reconciliation of remote and in situ observations of the upper atmosphere; the application of precise orbit determination of the spacecraft to measure upper atmospheric density and in conjunction with other Mars missions improve the static gravity field model of Mars; and an integrated ion/neutral study of ionospheric flows and resultant heavy ion escape. Descriptions of each of these investigations are given showing how each adds to and fits seamlessly into MAVEN mission science design.

Thermosphere↗

Science Enhancements by the MAVEN Participating Scientists

NASA implemented a Participating Scientist Program and released a solicitation for the Mars Atmosphere and Volatile EvolutioN mission (MAVEN) proposals on February 14, 2013. After a NASA peer review panel evaluated the proposals, NASA Headquarters selected nine on June 12, 2013. The program’s intent is to enhance the science return from the mission by including new investigations that broaden and/or complement the base-line investigations, while still addressing key science goals. The selections cover a broad range of science investigations. Included are: a patching of a 3D exosphere model to an improved global ionosphere-thermosphere model to study the generation of the exosphere and calculate the escape rates; the addition of a focused study of upper atmosphere variability and waves; improvement of a multi-fluid magnetohydrodynamic model that will be adjusted according to MAVEN observations to enhance the understanding of the solar-wind plasma interaction; a global study of the state of the ionosphere; folding MAVEN measurements into the Mars International Reference Ionosphere under development; quantification of atmospheric loss by pick-up using ion cyclotron wave observations; the reconciliation of remote and in situ observations of the upper atmosphere; the application of precise orbit determination of the spacecraft to measure upper atmospheric density and in conjunction with other Mars missions improve the static gravity field model of Mars; and an integrated ion/neutral study of ionospheric flows and resultant heavy ion escape. Descriptions of each of these investigations are given showing how each adds to and fits seamlessly into MAVEN mission science design.

Mars↗

Propagation of Interplanetary Disturbances in the Outer Heliosphere

Work finished during 2002 included: (1) Finished a multi-fluid solar wind model; (2) Determined the solar wind slowdown and interstellar neutral density; (3) Studied shock propagation and evolution in the outer heliosphere; (4) Investigated statistical properties of the solar wind in the outer heliosphere.

Wang, Chi↗

Physics of the Inner Heliosphere 1-10 Rs: Plasma Diagnostics and Models

While the mechanisms responsible for the solar corona and the high-speed solar wind streams are still unknown, model computations offer means of predicting the properties of such mechanisms in light of the empirical constraints currently available. Modeling and data analysis efforts were aimed at understanding the plasma properties of the acceleration of the solar wind, its filamentary nature, and the conditions needed to account for a rapidly accelerating solar wind, reaching its terminal speed within 10 R(sub s). A sequence of models ranging from steady one-fluid descriptions of the solar wind to multi-fluid time-dependent models were developed. Plasma diagnostics evolved from the analysis of data acquired from Skylab to SOHO, and complemented by ground-based observations.

Habbal, Shadia R.↗

Multi-Fluid Block-Adaptive-Tree Solar Wind Roe-Type Upwind Scheme: Magnetospheric Composition and Dynamics During Geomagnetic Storms, Initial Results

The magnetosphere contains a significant amount of ionospheric O{+}, particularly during geomagnetically active times. The presence of ionospheric plasma in the magnetosphere has a notable impact on magnetospheric composition and processes. We present a new multifluid MHD version of the BATS-R-US model of the magnetosphere to track the fate and consequences of ionospheric outflow. The multi-fluid MHD equations are presented as are the novel techniques for overcoming the formidable challenges associated with solving them. Our new model is then applied to the May 4, 1998 and March 31, 2001 geomagnetic storms. The results are juxtaposed with traditional single- fluid MHD and multispecies MHD simulations from a previous study, thereby allowing us to assess the benefits of using a more complex model with additional physics. We find that our multi-fluid MHD model (with outflow) gives comparable results to the multi-species MHD model (with outflow), including a more strongly negative Dst, reduced CPCP, and a drastically improved magnetic field at geosynchronous orbit, as compared to single-fluid MHD with no outflow. Significant differences in composition and magnetic field are found between the multi-species and multi-fluid approach further away from the Earth. We further demonstrate the ability to explore pressure and bulk velocity differences between H{+} and O(+}, which is not possible when utilizing the other techniques considered.

Gkocer, A.↗

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗

A First-Principles Model Based on Saturation of the Electron Cyclotron Drift Instability for Electron Transport in Hydrodynamics Simulations of Hall Thruster Plasmas

We present a computational model based on the hypothesis that the higher-than-expected electron mobility across magnetic field lines observed in Hall thrusters is due to the growth and eventual saturation of the electron cyclotron drift instability. The key improvement of this model with respect to previous work is that it decouples the saturation of the instability from the correlation that exists between mobility and magnitude of the wave perturbations in linear theory, thereby enabling first-principles simulations. The model has been incorporated in the 2-D (r-z) multi-fluids code Hall2De and simulations of the H6 thruster at 300 V, 20 A are numerically stable and achieve a steady-state solution at a computational cost that is not significantly higher than previous Hall2De simulations. The computed location of the acceleration region is found to be within 10% of the length of the acceleration channel compared to that inferred by experiments. The simulations also capture well the plasma gradients along the channel centerline of this thruster. To further establish the validity of this model, we plan to conduct simulations of other thrusters and/or operating conditions.

Lopez Ortega, Alejandro↗

Slow Solar Wind: Observations and Modeling

While it is certain that the fast solar wind originates from coronal holes, where and how the slow solar wind (SSW) is formed remains an outstanding question in solar physics even in the post-SOHO era. The quest for the SSW origin forms a major objective for the planned future missions such as the Solar Orbiter and Solar Probe Plus. Nonetheless, results from spacecraft data, combined with theoretical modeling, have helped to investigate many aspects of the SSW. Fundamental physical properties of the coronal plasma have been derived from spectroscopic and imaging remote-sensing data and in situ data, and these results have provided crucial insights for a deeper understanding of the origin and acceleration of the SSW. Advanced models of the SSW in coronal streamers and other structures have been developed using 3D MHD and multi-fluid equations.

Abbo, L.↗

Cometary comae

The basis of our knowledge about the structure and composition of a comet nucleus is examined. Such knowledge is inferred from observations of the cometary atmosphere (coma). However, photodissociation, photoionization, and chemistry destroy the mother molecules evaporating from the nucleus. To extract the primary information, the chemical kinetics and the physics of the coma are modeled with a computer and the results are compared with coma observations. The physics and chemistry for a dust free coma are described taking into account energy balance, multi-fluid flow for fast atomic and molecular hydrogen and the bulk fluid, and the transition from a collision dominated inner region to the free molecular flow outer region. Special attention is paid to the molecular data requirements for the current modeld and for extended models which will include solar wind interaction and dust. Such models are an important tool in support of the Giotto mission to Halley's comet, in the analysis and interpretation of coma observations, and in the understanding of the earliest history of the solar system.

Huebner, W. F.↗

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗

Preferential Heating of Oxygen 5+ Ions by Finite-Amplitude Oblique Alfven Waves

Minor ions in the fast solar wind are known to have higher temperatures and to flow faster than protons in the interplanetary space. In this study we combine previous research on parametric instability theory and 2.5D hybrid simulations to study the onset of preferential heating of Oxygen 5+ ions by large-scale finite-amplitude Alfven waves in the collisionless fast solar wind. We consider initially non-drifting isotropic multi-species plasma, consisting of isothermal massless fluid electrons, kinetic protons and kinetic Oxygen 5+ ions. The external energy source for the plasma heating and energization are oblique monochromatic Alfven-cyclotron waves. The waves have been created by rotating the direction of initial parallel pump, which is a solution of the multi-fluid plasma dispersion relation. We consider propagation angles theta less than or equal to 30 deg. The obliquely propagating Alfven pump waves lead to strong diffusion in the ion phase space, resulting in highly anisotropic heavy ion velocity distribution functions and proton beams. We discuss the application of the model to the problems of preferential heating of minor ions in the solar corona and the fast solar wind.

Maneva, Yana G.↗

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work is limited to very simple models. In this project, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. The proposed work is based on previously funded NASA work that allowed us to develop a two-dimensional numerical method for boiling flows and to demonstrate the ability of the method to simulate film boiling. While numerical simulations of multi-fluid flows have been advanced in a major way during the last five years, or so, similar capability for flows with phase change are still in their infancy. Although the feasibility of the proposed approach has been demonstrated, it has yet to be extended and applied to fully three-dimensional simulations. Here, a fully three-dimensional, parallel, grid adaptive code will be developed. The numerical method will be used to study nucleate boiling in microgravity, with particular emphasis on two aspects of the problem: 1) Examination of the growth of bubbles at a wall nucleation site and the instabilities of rapidly growing bubbles. Particular emphasis will be put on accurately capturing the thin wall layer left behind as a bubble expands along a wall, on computing instabilities on bubble surfaces as bubbles grow, and on quantifying the effects of both these phenomena on heat transfer; and 2) Examination of the effect of shear flow on bubble growth and heat transfer.

Tryggvason, Gretar↗

Plasma Sheet Circulation Pathways

Global simulations of Earth's magnetosphere in the solar wind compute the pathways of plasma circulation through the plasma sheet. We address the pathways that supply and drain the plasma sheet, by coupling single fluid simulations with Global Ion Kinetic simulations of the outer magnetosphere and the Comprehensive Ring Current Model of the inner magnetosphere, including plasmaspheric plasmas. We find that the plasma sheet is supplied with solar wind plasmas via the magnetospheric flanks, and that this supply is most effective for northward IMF. For southward IMF, the innermost plasma sheet and ring current region are directly supplied from the flanks, with an asymmetry of single particle entry favoring the dawn flank. The central plasma sheet (near midnight) is supplied, as expected, from the lobes and polar cusps, but the near-Earth supply consists mainly of slowly moving ionospheric outflows for typical conditions. Work with the recently developed multi-fluid LFM simulation shows transport via plasma "fingers" extending Earthward from the flanks, suggestive of an interchange instability. We investigate this with solar wind ion trajectories, seeking to understand the fingering mechanisms and effects on transport rates.

Moore, Thomas E.↗