Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Fusion-Enabled Pluto Orbiter and Lander

Direct Fusion Drive (DFD) is a unique fusion engine concept based on the Princeton Field-Reversed Configuration (PFRC) fusion reactor conceived by Dr. Sam Cohen of the Princeton Plasma Physics Laboratory. DFD would enable the Pluto orbiter and lander context mission and more broadly enable true "rapid transit" to outer-planet and near interstellar space. The truly game-changing levels of thrust and power in a modestly sized package could integrate with our current launch infrastructure while radically expanding the science capability of these missions. Our Phase I was our first funded work on the DFD, with previous work at PSS occurring only under internal R&D. We established the feasibility of our Pluto mission trajectories using straight-line and planar models, including a departure spiral from Earth and insertion at Pluto. We developed our first thrust and specific impulse model using the results of the UEDGE multi-fluid code. Our specific power model was improved. During this Phase II effort, we continued our efforts to increase the fidelity of the designs for the RF, magnet, and shielding subsystems. Dedicated thrust augmentation experiments were run on the PFRC experiment, using a supersonic gas puffing valve. For the first time, we analyzed the design of a closed-loop operation mode and estimated the hardware that would be required for a dual-mode engine. In an exciting new development, we have invented a new thermophotovoltaic thermal conversion method that has the potential to have efficiencies of a Brayton or Stirling system. Our report presents details of these analyses. Our roadmap to bringing DFD to flight predicts that with sufficient support, a first flight unit could be built by 2040. We anticipate that three machine generations are required before this point: a ~1 T PFRC-3 machine hitting new plasma temperature and density levels, a ~5 T PFRC-4 machine with first demonstration of D-3He fusion, and a flight. (The current experiment, PFRC-2, is limited to about 0.1 T). In order to achieve a flight in 2035-2040, the TRL of the supporting systems must be increased in parallel, including low mass radiators, cryogenic propellant storage, and large (>100 kW) thermal conversion systems. Fortunately, many of these systems are dual-use and are required for other technologies including fission systems, so DFD would contribute to and benefit from those programs. This NIAC support and the results of our work have led to multiple follow-on contracts. We won two NASA STTRs focused on DFD subsystems, one on the RF system and one on the superconducting magnets, and our magnet STTR is now midway through a Phase II. We will be receiving a superconducting test magnet for experiments at PPPL this summer. In addition, we won an ARPA-E OPEN grant; this is the first time that OPEN has included fusion technologies, and we are part of a cohort of three alternative fusion companies now supported directly by DOE. We are very optimistic that if we are able to meet our experimental milestones in the next 12-18 months, we will be competitive for additional DOE grants to build PFRC-3.

Propulsion↗
NASA NTRS2002

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↗
DOE OSTI2025

Ponderomotive electron physics captured in a single-fluid extended magnetohydrodynamics model

The ponderomotive force, arising from the interaction between electromagnetic waves and plasma, plays a critical role in laser fusion, astrophysical plasmas, and laser diagnostics. Traditionally, modeling this force requires multi-fluid or particle-in-cell simulations due to its strong coupling to electron-scale dynamics. In this work, we demonstrate that a one-fluid, two-temperature extended magnetohydrodynamics (XMHD) model—augmented with a generalized Ohm's law (GOL) including electron inertia—can accurately reproduce key ponderomotive effects. We derive the ponderomotive force within this framework using a phasor-based approach and then validate its nonlinear manifestations through direct numerical simulations in the PERSEUS code, where steepening and density modulation phenomena typically associated with kinetic-scale models are reproduced. These results and prior work establish XMHD as a robust and efficient alternative for modeling nonlinear laser–plasma dynamics, bridging the gap between ideal MHD and fully kinetic approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

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

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

Self-Consistent Implementation of a Zero-Equation Transport Model Into a Predictive Model for a Hall Effect Thruster

The performance of an axisymmetric multi-fluid Hall thruster code that incorporates a self-consistent, data-driven closure model for the anomalous electron transport is investigated. Five different operating conditions of the H9 magnetically shielded Hall thruster are simulated with the Jet Propulsion Laboratory’s Hall2De. In order to capture the inherent uncertainty associated with the closure model, O(100) simulations are run for each condition, each of using a coefficient set sampled randomly from a probability distribution. The results of these simulations provide probabilistic predictions of thruster performance quantities including thrust, and discharge current, as well as several component efficiencies and centerline plasma properties. The model is found to yield converged solutions at all conditions, with large 10 kHzrange oscillations and performance trends with voltage and flow rate similar to experiment. The model under-predicts the thrust by 15-25% and over-predicts the discharge current by 20% on average compared to experiments at the same discharge voltage and mass flow rate. This performance discrepancy is due to lower beam utilization, mass utilization, and divergence efficiency than experiment, resulting from high Hall parameters in the acceleration region, which lead to a protracted ion acceleration region. The physical processes underlying this result are discussed in the context of future data-driven modeling efforts.

Jorns, Benjamin A.↗
NASA NTRS2010

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↗
NASA NTRS2017

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↗
NASA NTRS2016

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

Multiscale Hybrid Modeling of the Impact Response of the Earth's Magnetotail to Ionospheric O+ Outflow

Ionospheric outflow plays an important role in coupling the ionosphere with the solar wind-magnetosphere system. Previous multi-fluid MHD studies explored the global influence of oxygen ions of ionospheric origin (O+) on magnetospheric dynamics. A detailed exploration of the interaction of ionospheric ions with the magnetotail requires kinetic treatment for ions. We perform a self-consistent investigation of these processes with a three-dimensional space-time adaptive hybrid code, HYPERS, powered by an intelligent Event-driven Multi-Agent Planning System (EMAPS). By comparing simulations with and without outflow we conclude that oxygen ions, flowing from the ionosphere through the lobes into the tail or directly entering the inner magnetosphere, are able to significantly modify the magnetotail configuration and induce X-points and current sheet structures at locations where magnetic reconnection does not occur in a simulation without outflow, potentially very close to the Earth. This finding may have implications for interpreting substorms and magnetotail reconnection events observed for southward magnetic field simultaneously with significant contents of oxygen ions of ionospheric origin.

multiscale↗
NASA NTRS1985

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

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling a divertor with mid-leg pumping for high-power H-mode scenarios in DIII-D considering E x B drift flows

Edge-plasma simulations of a baffled, long-legged divertor in DIII-D, performed using the multi-fluid code UEDGE, indicate that the position of the detachment front is constrained to the location of the pump duct along the low-field side (LFS) baffle. Simulations including magnetic and E x B drifts were performed for 12.5 MW deuterium plasmas including intrinsic carbon and seeded neon to assess the optimal location of the LFS divertor pump to create a stable detachment front between the target and the X-point. The radiation front position in the simulations, taken to be indicative of the detachment front, can be controlled between the pump and X-point in the favorable magnetic field direction for H-mode access by moving the pump duct location upstream of the target along the LFS baffle. In the unfavorable magnetic field direction, the radial E θ x B drift flows are directed towards the pumping surface, efficiently removing the injected deuterium gas and limiting the sensitivity of the radiation front location to the gas injection rate. The role of pumping rate and drift direction on the pumping efficiency are also found to affect the divertor plasma conditions and detachment front location in UEDGE simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2014

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↗
NASA NTRS2016

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

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↗
NASA NTRS2008

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