A Unified Hyperbolic Model for Coupled Fluid Dynamics and Electromagnetics in Aerospace Applications (original) (raw)
Related papers
A Magneto-Fluid-Dynamic Model and Computational Solving Methodologies for Aerospace Applications
Fluid Dynamics, Computational Modeling and Applications, 2012
Computational plasma physics is concerned primarily with the study of the evolution of plasma by means of computer simulation. The main task of this computational branch is to develop methods able to obtain a better understanding of plasma physics. Therefore a close contact to theoretical plasma physics and numerical methods is necessary. Ideally computational plasma physics acts as a pathfinder to guide scientific and technical development and to connect experiment and theory. To build a valid computer simulation program means to devise a model which is sufficiently detailed to reproduce faithfully the most important physical effects, with a computational effort sustainable by modern computers in reasonable time . Computational models have played an important role in the development of plasma physics since the beginning of the computer age. Advances in our understanding of many plasma phenomena like magnetohydrodynamic instabilities, micro-instabilities, transport, wave propagation, etc. have gone hand-in-hand with the increased computational power available to researchers. Several trends are evident in how computer modelling is carried out: the models are becoming increasingly complex, for example, by coupling separate computer codes together. This allows for more realistic modelling of the plasma. Presently several efforts are carried out in different countries to develop plasma numerical tools for several applications such as fusion, electric propulsion, active control over hypersonic vehicles: these efforts lead to a growing experience in CMFD field (see Park et al. The chapter presented was carried out in the context of a research activity motivated by renewed interest in investigating the influence that electromagnetic fields can exert on the thermal and pressure loads imposed on a body invested by a high energetic flow. In this regard, spacecraft thermal protection and the opportunity to use active control surfaces during planetary (re)entry represent the driving engineering applications. The contents of the study should be considered, to a certain extent, a systematic re-examination of past work complemented with somewhat innovative ideas. So, in this chapter, methodologies for plasma modelling have been developed and then implemented and tested into a numerical code EMC3NS, developed in the frame of this
Detailed Modeling of Plasmas for Computational Aerodynamics
2016
As demonstrated by Parent, B., et al., (“Electron and Ion Transport Equations in Computational Weakly-Ionized Plasmadynamics,” Journal of Computational Physics, Vol. 259, 2014, pp. 51–69), the computational efficiency of the drift-diffusion plasma model can be increased significantly by recasting the equations such that the potential is obtained from Ohm’s law rather than Gauss’s law and by adding source terms to the ion transport equations to ensure that Gauss’s law is satisfied. Not only did doing so reduce the stiffness of the system, leading to faster convergence, but it also resulted in a higher resolution of the converged solution. The combined gains in convergence acceleration and resolution amounted to a hundredfold increase in computational efficiency when simulating nonneutral plasmas with significant quasi-neutral regions. In this paper, it is shown that such a recast of the drift-diffusion model has yet another advantage: its lack of stiffness permits the electron and ion transport equations to be integrated in coupled form along with the Favre-averaged Navier–Stokes equations. Test cases relevant to plasma aerodynamics (including nonneutral sheaths, magnetic field effects, and negative ions) demonstrate that the proposed coupled system of equations can be converged in essentially the same number of iterations as that describing nonionized flows while not sacrificing the generality of the drift-diffusion model.
A flux-limited numerical method for the MHD equations to simulate propulsive plasma flows
31st Plasmadynamics and Lasers Conference
For numerical simulations to be effective tools in plasma propulsion research, a higher order accurate solver that captures MHD shocks monotonically and works reliably for strong magnetic fields is needed. For this purpose, a characteristics based scheme for the MHD equations, with flux limiters to improve spatial accuracy, has been developed. In this method, the symmetric form of the MHD equations, accounting for waves propagating in all directions, is solved. The required eigensystem of axisymmetric MHD equations, with appropriate normalization, is presented. This scheme was validated against unsteady (Riemann problem) and force-free equilibrium (Taylor state) test cases. The capability of this method to simulate resistive plasma flows is demonstrated using a simple example.
A flux-limited numerical method for solving the MHD equations to simulate propulsive plasma flows
International Journal for Numerical Methods in Engineering, 2002
For numerical simulations to be e ective tools in plasma propulsion research, a high-order accurate solver that captures MHD shocks monotonically and works reliably for strong magnetic ÿelds is needed. For this purpose, a characteristics-based scheme for the MHD equations, with ux limiters to improve spatial accuracy, has been developed. In this method, the symmetric form of the MHD equations, accounting for waves propagating in all directions, are solved. The required eigensystem of axisymmetric MHD equations, with appropriate normalization, is presented. This scheme was validated with unsteady (Riemann problem) and force-free equilibrium (Taylor state) test cases, as well as with measured current density patterns in a magnetoplasmadynamic thruster. Copyright ? 2001 John Wiley & Sons, Ltd.
A solution procedure for the Navier-Stokes equations coupled with the full Maxwell equations is described. The approach implements a strongly conservative fluid formulation in which the Lorentz force and Ohmic heating terms are recast as convective terms. This removes explicit sources from the Navier-Stokes equations, which have previously introduced severe sti ness and demanded a very delicate numerical treatment. The coupling with the full Maxwell equations enables the displacement current to be incorporated directly. To demonstrate the e ectiveness of this technique, a fully explicit finite volume approximate Riemann solver is used to obtain numerical solutions to the Brio and Wu electromagnetic shock problem. Comparisons with the analytical solution show good agreement, and the implementation requires less than an hour of computational time on a single processor machine. Simulations using large and small conductivities confirm that the formulation captures both wave and di usion limits of the magnetic field.
Fully-implicit finite volume method for the ideal two-fluid plasma model
Computer Physics Communications, 2018
We present a novel numerical model that simulates ideal two-fluid plasmas coupled to the full set of Maxwell's equations with application to space and laboratory plasmas. We use a fully-implicit finite volume method for unstructured meshes, that uses an advection upstream splitting method (i.e., AUSM +up) for all speeds to discretize the numerical fluxes of the fluids. In addition, we discretize Maxwell's equations with a modified-Rusanov scheme. The electromagnetic numerical dissipation is scaled using the scales of the fluid-electromagnetics coupled problem that are found to be very different from those of the uncoupled problem. Our numerical scheme guarantees that the elliptical constraints of Maxwell's equations are satisfied by using hyperbolic divergence cleaning. We validate the performance and accuracy of our model by simulating the following conventional cases: a circularly polarized wave, a Brio-Wu type shock tube, and a two-fluid plasma reconnection with the GEM challenge set up. Our model reveals the complexity of the two-fluid model compared to magnetohydrodynamics (MHD) models, as the inclusion of charge separation, the displacement current and the electron dynamics present are ignored by the MHD simplifications. The two-fluid model shows the presence of electromagnetic and plasma waves and the effect that they have in even the simplest cases. We also compare our model to other available two-fluid models and find our results to be in good agreement.
Numerical Modeling of Plasma – Flow Interaction
Lecture Notes in Computer Science, 2006
In the frame of the internal project PUMA (Plasma Used to Master Aerodynamics), ONERA is conducting fundamental studies of plasma-flow interactions. In this paper, the ionic wind created by corona discharges is studied in the case of a subsonic flow over a flat plate. The proposed mechanism of the ionic wind proposed is the addition of momentum by collisions between charged and neutral particles. In order to evaluate the effect of plasma on aerodynamics, a kinetic modeling of the discharge is coupled with a Fluid Dynamics code.
Modeling and Simulating Flowing Plasmas and Related Phenomena
Space Science Reviews, 2008
Simulation has become a valuable tool that compliments more traditional methods used to understand solar system plasmas and their interactions with planets, moons and comets. The three popular simulation approaches to studying these interactions are presented. Each approach provides valuable insight to these interactions. To date no one approach is capable of simulating the whole interaction region from the collisionless to the collisional regimes. All three approaches are therefore needed. Each approach has several implicit physical assumptions as well as several numerical assumptions depending on the scheme used. The magnetohydrodynamic (MHD), test-particle/Monte-Carlo and hybrid models used in simulating flowing plasmas are described. Special consideration is given to the implicit assumptions underlying each model. Some of the more common numerical methods used to implement each model, the implications of these numerical methods and the resulting limitations of each simulation approach are also discussed.
40th AIAA Plasmadynamics and Lasers Conference, 2009
A Navier-Stokes solver for non-equilibrium aerothermodynamics is coupled with a fully implicit solver of the Maxwell equations. Transport phenomena are modelled using a high accurate Chapman-Enskog method that accounts for the effects of concentration, pressure and temperature gradients and provides anisotropic diffusion velocity, heat flux and viscosity in the presence of a magnetic field. Gas neutrality is not enforced so that the electric charge density can be locally different from zero. Numerical tests are carried out on a flatfaced cylinder to test the capabilities of the implicit Maxwell solver and of the coupled Maxwell/Navier-Stokes system.
An implicit method for electromagnetic plasma simulation in two dimensions
Journal of Computational Physics, 1982
A new method for modeling low-frequency plasma phenomena is presented. The method uses an implicit formulation of the Vlasov-Maxwell equations to relax restrictions on the time-step and mesh spacing so that larger values which correspond to the frequencies and wavelengths of interest can be used. As a result, the range of length and time scales accessible to plasma simulation is increased by orders of magnitude. The algorithm, as embodied in a new code VENUS for electromagnetic plasmas in two dimensions, is described, its stability and accuracy analyzed through linear and nonlinear analysis, and its properties, including suppression of the finite grid instability, illustrated through its application to the Weibel instability.