Fluid Model of Plasma and Computational Methods for Solution
|
|
- Lora Fleming
- 5 years ago
- Views:
Transcription
1 WDS'06 Proceedings of Contributed Papers, Part III, , ISBN X MATFYZPRESS Fluid Model of Plasma and Computational Methods for Solution E. Havlíčková Charles University, Faculty of Mathematics and Physics, V Holešovičkách 2, Prague, Czech Republic. Abstract. This contribution gives basic insight into the fluid modelling in plasma physics. The main aim of the paper is to offer brief conspectus of equations that form the fluid model of plasma and to point out the differences between the theoretical model in general form as derived from the Boltzmann equation and the fluid model which is commonly used in practical simulations. Second part of the contribution summarizes computational methods which are often used for solution of the described model, as they were the most widely mentioned in literature. Introduction The simulation of plasma processes can be based on various approaches. The microscopic description of plasma is for many practical purposes too detailed or computationally too complicated and in some applications we need not to know the information about the behaviour of individual particles. Fluid model describes macroscopic plasma phenomena and reveals how the statistical plasma parameters evolve in time and space. Fluid model of plasma Theoretical description of plasma Fluid model of plasma is based on partial differential equations (PDEs) which describe the macroscopic quantities such as density, flux, average velocity, pressure, temperature or heat flux. Governing PDEs can be derived from the Boltzmann equation (BE) by taking velocity moments. Zero moment of the Boltzmann equation ( + (BE)d3 v) yields continuity equation for the particle density n t + (n u) = Qρ m Density and average velocity are defined as n = + f d3 v and u = 1 + n vf d3 v. The source term on the right side of the equation (1) corresponds to the collision term of the Boltzmann equation and describes mass production and annihilation. It is defined as Q ρ = m + ( f/ t) c d3 v. Similarly, the momentum transport equation (equation of motion) (2, 3) and heat transport equation (energy equation) (4, 5) can be found as first moment (m + v (BE)d3 v) and second moment ( 1 2 m + v2 (BE)d 3 v) of the Boltzmann equation. m (n u) t + P +m (n u u) nq( E + u B) = uq ρ + Q p (2) mn u t + mn( u ) u + P nq( E + u B) = Q p (3) ( 1 t γ 1 p + 1 ) ( 1 2 mnu2 + γ 1 p u + u P + 1 ) 2 mnu2 u + L nqe u = 1 2 u2 Q ρ + u Q p + Q E (4) ( ) 1 p γ 1 t + (p u) + ( P ) u + L = Q E (5) Mean quantities are defined by the following equations. L = 1 2 m + ( v u)( v u)2 f d 3 v defines heat flux L, P = m + ( v u)( v u)f d3 v defines pressure tensor P and P ij = pδ ij defines scalar pressure p. γ in (4) and (5) is the ratio of specific heats. Source terms in (2) and (4) describe transport of (1) 180
2 momentum and heat transfer due to collisions and they are defined as Q p = m + ( v u)( f/ t) c d3 v, Q E = 1 2 m + ( v u)2 ( f/ t) c d 3 v. We can continue to derive the moments for high order terms, however the equation chain must be truncated somewhere. In many practical problems this is made in the first order by substituting the energy equation by an equation of the state, or in the second order by using algebraic expression for the heat flux as closure approximation. Therefore the equations (1), (2) and (4) establish the typical set of fluid equations which are used in many simulations. Moreover, to complete the system of governing equations and to obtain the self-consistent description of plasma, we must consider the Maxwell s equations describing the electromagnetic behaviour of plasma. The meaning of individual terms in the fluid equations and their derivation is explained in Golant et al. [1980], Chen [1984] or Braginskii [1965] in more detail. Specific application of the fluid theory leads to various forms of the model. Generally we can use one set of fluid equations for each plasma species and considering a simple two component (ion and electron) plasma one obtains the so-called two-fluid model. Magnetohydrodynamics, the most widely known plasma theory, is another modification - plasma is considered as a single fluid in the center of mass frame [Chen, 1984; Vold et al., 1991]. A choice of the model depends on the character of the simulated plasma. Classical fluid formulation Practical use of the fluid model of plasma is not based on the general form of the equations (1), (2) and (4) as they were derived from the Boltzmann equation, but it is always connected with various approximations. Therefore let s introduce the classical fluid modelling approach used most commonly in simulations [Chen et al., 2004; Bukowski et al., 1996]. Approximative form of the equations is connected especially with a simplification of the source terms which are very complicated functions of velocity. The mass balance collision term gives the rates of creation and loss of species. The productivity of the reactions is defined by the reaction rate coefficients k r corresponding to the collision of type r. The coefficients k r are the input parameters of the simulation. Assuming an electron collision with neutral, the source term for the formation of species k can be evaluated as Q ρ k = m k r l rkk r (T e )n e n n, where l rk is the number of particles of species k created or lost per collision type r [Bukowski et al., 1996]. The useful approximation of the collision term in the momentum transport equation (2) and the energy equation (4) is the Krook s approximation Q p k = m k m l l m k +m l n k ( u k u l ) ν kl and Q E k = l 3 2 k n k collision frequency [Golant, 1980; Chen et al., 2004]. Fluid equations as used in simulations are often simplified by neglecting some terms: 2m k m l (m k +m l ) 2 (T k T n ) ν kl, where ν kl is the mean - neglecting viscosity effect in the transfer equation characterizing the anisotropic part of the pressure, which yields the reduction of pressure tensor to the scalar pressure [Golant et al., 1980]. - neglecting kinetic energy contribution compared to the thermal one, which leads to the significant simplification of the energy equation [Bukowski et al., 1996; Chen et al., 2004]. - neglecting the first term in the thermal flux closure for electrons L e = 5 2 Γ e kt e 5 2 kd en e T e, which is known as Fourier s approximation [Chen et al., 2004]. The expressions for thermal flux of various species can be found in Golant et al. [1980]. The most widely mentioned approximation in the classical fluid formulation is the so-called driftdiffusion approximation [Herrebout et al., 2001; Herrebout et al., 2002; Chen et al., 2004; Boeuf, 1987; Donkó, 2001; Passchier et al., 1993]. This approximation reduces the number of partial differential equations included in model by the use of the algebraic expression for particle flux Γ k (7) instead of full equation of motion (2). The equation (7) describes the transport of species due to density gradient and the transport of charged species under the influence of the electric field. n k t + Γ k = Qρ k m (6) Γ k = ±µ k n k E (Dk n k ) (7) Coefficients of mobility µ and diffusion D are input parameters of the simulation. The expression (7) is equivalent to the equation of motion (3) after neglecting unsteady (first term) and inertial (second term) contributions. Conditions under which the drift-diffusion approximation and the approximations 181
3 mentioned above are fulfilled can be found in Golant et al. [1980] and Chen et al. [2004]. These publications also include more detailed characteristics of the individual terms in the equations and the definition of all quantities. Computational methods of solution The physical aspects of any fluid flow are governed by three fundamental principles - mass conservation, Newton s second law and energy conservation. These principles can be expressed in terms of mathematical equations, which are usually partial differential equations. This section briefly resumes the methods of computational fluid dynamics (CFD) that give the solution of PDEs describing the plasma as a fluid. The basic computational techniques are based on the construction of a discrete grid and the replacement of individual differentiated terms in PDEs by algebraic expressions connecting nodal values on a grid (Fig. 1). GOVERNING PARTIAL DIFF. EQS. DISCRETIZATION SYSTEM OF ALGEBRAIC EQUATIONS EQUATION SOLVER APPROXIMATE SOLUTION Figure 1. Overview of the computational solution procedure. The most common choices for converting the PDE to the algebraic one are finite difference, finite element and finite volume methods. In practice, time derivations in the time-dependent equation are discretized almost exclusively using the finite difference method and spatial derivatives are discretized by either the finite difference, finite element or finite volume method, typically. Finite difference method The method of finite differences (FDM) is widely used in CFD [Meeks et al., 1993; Montierth et al., 1992; Baboolal, 2002]. Finite difference representation of derivatives is based on Taylor s series expansions. The discretization process introduces an error dependent on the order of terms in the Taylor s series which are truncated. Derivation of elementary finite differences of various order accuracy is described in Wendt [1992]. To represent the derivatives by the differences, a number of choices is available, especially when the dependent variable appearing in the governing equation is a function of both coordinates and time. In this case, the finite difference approaches can be divided into implicit and explicit ones [Wendt, 1992]. The explicit approaches are relatively simple to set up and program but there are stability constraints, given by Von Neumann method, which can result in long computer running times. On the other hand the implicit approaches are stable even for larger values of the time step, however a system of algebraic equations must be solved at each time step and implicit techniques are more complicated to implement. Publications devoted to fluid modelling in plasma mention a variety of methods based on the finite difference scheme, namely the explicit Lax-Wendroff and MacCormack s methods or implicit Crank- Nicholson scheme, all described in Wendt [1992] in more detail. The summary of some techniques is given in Vold et al. [1991]. As mentioned in the first section, the drift-diffusion approximation is widely used to simplify the fluid model of plasma. Then the continuity equation is of convection-diffusion type n t + Γ x = 0 (8) Γ = v D n D n (9) x In equation (9) v D = ±µe and µ and D are transport coefficients. Standard discretization of the convection-diffusion equation requires very small grid spacing x to be stable, therefore several special discretization schemes have been developed. Scharfetter-Gummel implicit scheme given by equations (10), (11) and (12), which is the most popular one [Boeuf, 1987; Chen et al.; 2004, Passchier et al., 1993; Mareš et al., 2002], provides an optimum way to discretize the drift-diffusion equation for particle transport. It has been developed in the frame of semiconductor device simulations [Scharfetter and Gummel, 1969]. n k+1 i n k Γ k+1 Γ k+1 i i i 1 2 = 0 (10) t x 182
4 Γ i+ 1 2 = v D [ ni+1 n i exp ( vd x D )] 1 exp ( vd x D ) (11) Γ i 1 = v [ D ni n i 1 exp ( vd x D )] 2 1 exp ( vd x D ) (12) Examing the limiting cases of the scheme, for diffusion-dominated problems the limit gives classical central difference second-order scheme and for convection-dominated problems the limit produces the upwind difference form of the first-order accuracy [Passchier et al., 1993]. Finite element method First essential characteristic of finite element method (FEM) is to divide the continuum field (domain) into non-overlapping elements. The elements have either a triangular or a quadrilateral form, they can be curved and cover the whole domain. The grid formed by elements need not be structured, in opposite to FDM. The basic philosophy of FEM is that an approximate solution of the discrete problem is assumed a priori to have prescribed form, the solution has to belong to a function space u = N u j φ j (x,y,z) (13) j=1 In equation (13) u is approximative solution, u j are unknown coefficients and φ j are basis (shape) functions. The basis functions are chosen almost exclusively from low-order piecewise polynomials and are assumed to be non-zero in the smallest possible number of elements associated to the index of the basis function (Fig. 2), which is computationally advantageous. Finite element method does not look for a solution of the PDE itself, but looks for a solution of integral form of the PDE obtained from a weighted residual formulation [Wendt, 1992; Fletcher, 1991]. The unknown coefficients u j are determined by requiring that the integral of the weighted residual of the partial differential equation R over the computational domain Ω is zero W m (x,y,z)rdx dy dz = 0 m = 1,2,...,N (14) Ω Different choices for the weight function W m in (14) give rise to different methods in the class of methods of weighted residuals. In FEM the weight functions are chosen from the same family as the basis functions W m = φ m, which is the most popular choice. Equation (14) results in a system of algebraic equations for coefficients u j, or a system of ordinary differential equations for time-dependent problems. The finite element approach is applied in Charrada et al. [1996] for example. Figure 2. One-dimensional linear approximating functions [Fletcher, 1991]. Finite volume method The basic idea of the finite volume method (FVM) is to discretize the integral form of the equations instead of the differential form and it can be thought as a special case of the so-called subdomain 183
5 method [Fletcher, 1991]. The computational domain is subdivided into a set of cells that cover the whole domain. The volumes on which the integral forms of conservation laws (such as continuity equation or momentum transport equation) are applied need not coincide with the cells of the grid and they can even be overlapping. Different choices of volumes (Fig. 3) determine different formulations of the finite volume technique. nodes a b c cell-centered volume cell-vertex volume (non-overlapping) cell-vertex volume (overlapping) cell-vertex volume (on triangular cells) Figure 3. Typical choice of volumes in FVM [Wendt, 1992]. Let s consider the cell-vertex formulation (Fig. 4) and the PDE of the form of equation u t + f x + g y = 0 (15) Integrating the equation (15) over the control volume ABCD (Fig. 4) and applying Green s theorem yields the following equation d u dv + (f dy g dx) = 0 (16) dt ABCD Discretizing the equation (16) to obtain an approximate evaluation of the integral form of the governing PDE, one ends with equation (17) for each nodal point (j,k). The equation (17) then leads to the discretization scheme of the cell-vertex formulation of FVM. S ABCD du j,k dt + DA AB (f y g x) = 0 (17) Detailed analysis of steps that result in the discretization scheme is described in Fletcher [1991] and Wendt [1992]. Application of the finite volume approach can be found in Maruzewski et al. [2002] for instance. Figure 4. Two-dimensional finite volume [Fletcher, 1991]. 184
6 Discussion HAVLÍČKOVÁ: FLUID MODELLING IN PLASMA PHYSICS Both finite difference techniques and the methods based on the weighted residual formulation are widely used in simulations. The advantage of FDM is its relatively simple implementation, especially in the case of problems which don t require to transform the coordinates. By contrast, practical use of FEM and FVM in simulations is usually connected with commercial solvers such as Fluent [Freton et al., 2000; Gonzalez et al., 2002; Blais et al., 2003] or Comsol Multiphysics [Bartoš et al., 2005]. Comparison of the individual methods can be deduced from Wendt [1992] and Fletcher [1991], who demonstrate various applications of the methods by particular examples. The most important advantage of FEM and FVM is given by the possibility to use unstructured grid. Due to the unstructured form, very complex geometries can be handled with ease. This geometric flexibility is not shared by FDM [Wendt, 1992]. In addition, FVM provides a simple way of discretization without the need to introduce generalized coordinates even when the global grid is irregular [Fletcher, 1991]. The use of FVM is also supported by situations, where the conservation laws can not be represented by PDEs but only the integral forms are guaranteed (discontinuities, etc.). However, the main problem of FVM are difficulties in the definition of derivatives which can not be based on a Taylor-expansion. FVM becomes more complicated when applied to the PDE containing the second derivatives, therefore the FVM is best suitable for flow problems in primitive variables where the viscous terms are absent or are not very important [Wendt, 1992]. Accuracy of the particular methods is analyzed by Fletcher [1991] and Wendt [1992]. As a rough guide, the use of linear approximating functions in FEM generates solutions of about the same accuracy as second-order FDM and the accuracy of FEM with quadratic approximating functions is comparable with third-order FDM [Fletcher, 1991]. The accuracy of the finite volume techniques is determined by the formulation type and can depend on irregularity of the grid. Generally, finite volume approaches are first-order or second-order accurate in space [Wendt, 1992]. Conclusion Presented contribution provides a brief description of the fluid model of plasma and reviews elementary computational approaches. Only the basic aspects of the methods were described with emphasis on the comparison of the individual techniques and with a view to accentuate their essential advantages. Extended analysis of both the fluid model of plasma and the computational methods for solution is included in the referenced publications. Acknowledgments. The work is a part of the research plan MSM that is financed by the Ministry of Education of the Czech Republic and was partly supported by the Grant Agency of Charles University Prague, Grants No. 296/2004 and 220/2006. References Baboolal S., Boundary conditions and numerical fluid modelling of time-evolutionary plasma sheaths, Journal of Physics D: Applied Physics, 35, 658, Bartoš P., Jelínek P., Hrach R., Fluid Modelling of Plasma Behaviour in the Vicinity of Planar and Cylindrical Probes, WDS 05 Proceedings on Contributed Papers, Part III, 625, Matfyzpress, Blais A., Proulx P., Boulos M. I., Three-dimensional numerical modelling of a magnetically deflected dc transferred arc in argon, Journal of Physics D: Applied Physics, 36, 488, Boeuf J.-P., Numerical model of rf glow discharges, Physical Review A, 36, 2782, Braginskii S. I., Transport processes in a plasma, Reviews of Plasma Physics, 1, 205, 1965 Bukowski J. D., Graves D. B., Two-dimensional fluid model of an inductively coupled plasma with comparison to experimental spatial profiles, Journal of Applied Physics, 80, 2614, Charrada K., Zissis G., Aubes M., Two-temperature, two-dimensional fluid modelling of mercury plasma in high-pressure lamps, Journal of Physics D: Applied Physics, 29, 2432, Chen F. F., Introduction to Plasma Physics, Academia, Prague, 1984 Chen G., Raja L. L., Fluid modeling of electron heating in low-pressure, high-frequency capacitively coupled plasma discharges, Journal of Applied Physics, 96, 6073, Donkó Z., Apparent secondary-electron emission coefficient and the voltage-current characteristics of argon glow discharges, Physical Review E, 64, , Fletcher C. A. J., Computational Techniques for Fluid Dynamics, Vol. I, Springer-Verlag, Freton P., Gonzalez J. J., Gleizes A., Comparison between a two- and a three-dimensional arc plasma configuration, Journal of Physics D: Applied Physics, 33, 2442, Golant V. E., Zhilinsky A. P., Sakharov I. E., Fundamentals of Plasma Physics, Wiley, New York,
7 Gonzalez J. J., Freton P., Gleizes A., Comparisons between two- and three-dimensional models: gas injection and arc attachment, Journal of Physics D: Applied Physics, 35, 3181, Herrebout D. et al., One-dimensional fluid model for an rf methane plasma of interest in deposition diamond-like carbon layers, Journal of Applied Physics, 90, 570, Herrebout D. et al., Modeling of a capacitively coupled radio-frequency methane plasma: Comparison between a one-dimensional and a two-dimensional fluid model, Journal of Applied Physics, 92, 2290, Mareš R., Vicher M., Hrach R., PIC and fluid model of plasma sheath, Czechoslovak Journal of Physics, 52, D705, Maruzewski P. et al., Simulation of arc-electrode interaction using sheath modelling in SF6 circuit-breakers, Journal of Physics D: Applied Physics, 35, 891, Meeks E., Cappelli M. A., Two-temperature fluid model for high-pressure plasmas in contact with cooled electrodes, Journal of Applied Physics, 73, 3172, Montierth L. M., Neuman W. A., Morse R. L., Fluid model treatment of surface plasma structures, Physics of Fluids B, 4, 784, Passchier J. D. P., Goedheer W. J., A two-dimensional fluid model for an argon rf discharge, Journal of Applied Physics, 74, 3744, Scharfetter D. L., Gummel H. K., Large-Signal Analysis of Silicon Read Diode Oscilator, IEEE Transaction on Electron Devices, 16, 64, Vold E. L., Najmabadi F., Conn R. W., Fluid model equations for the tokamak plasma edge, Physics of Fluids B, 3, 3132, Wendt J. F. (editor) et al., Computational Fluid Dynamics - An Introduction, Springer-Verlag,
The Computational Simulation of the Positive Ion Propagation to Uneven Substrates
WDS' Proceedings of Contributed Papers, Part II, 5 9,. ISBN 978-8-778-85-9 MATFYZPRESS The Computational Simulation of the Positive Ion Propagation to Uneven Substrates V. Hrubý and R. Hrach Charles University,
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
More informationPlasma Modeling with COMSOL Multiphysics
Plasma Modeling with COMSOL Multiphysics Copyright 2014 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
More informationLattice Boltzmann Method
3 Lattice Boltzmann Method 3.1 Introduction The lattice Boltzmann method is a discrete computational method based upon the lattice gas automata - a simplified, fictitious molecular model. It consists of
More informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationSimulations of Electrical Arcs: Algorithms, Physical Scales, and Coupling. Henrik Nordborg HSR University of Applied Sciences Rapperswil
Simulations of Electrical Arcs: Algorithms, Physical Scales, and Coupling Henrik Nordborg HSR University of Applied Sciences Rapperswil What is an electrical arc? 2 Technical applications of arcs and industrial
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
More informationEffect of Gas Flow Rate and Gas Composition in Ar/CH 4 Inductively Coupled Plasmas
COMSOL CONFERENCE BOSTON 2011 Effect of Gas Flow Rate and Gas Composition in Ar/CH 4 Inductively Coupled Plasmas Keisoku Engineering System Co., Ltd., JAPAN Dr. Lizhu Tong October 14, 2011 1 Contents 1.
More informationThe Role of Secondary Electrons in Low Pressure RF Glow Discharge
WDS'05 Proceedings of Contributed Papers, Part II, 306 312, 2005. ISBN 80-86732-59-2 MATFYZPRESS The Role of Secondary Electrons in Low Pressure RF Glow Discharge O. Brzobohatý and D. Trunec Department
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationMultidimensional Numerical Simulation of Glow Discharge by Using the N-BEE-Time Splitting Method
Plasma Science and Technology, Vol.14, No.9, Sep. 2012 Multidimensional Numerical Simulation of Glow Discharge by Using the N-BEE-Time Splitting Method Benyssaad KRALOUA, Ali HENNAD Electrical Engineering
More informationPhysics of fusion power. Lecture 13 : Diffusion equation / transport
Physics of fusion power Lecture 13 : Diffusion equation / transport Many body problem The plasma has some 10 22 particles. No description is possible that allows for the determination of position and velocity
More informationOne dimensional hybrid Maxwell-Boltzmann model of shearth evolution
Technical collection One dimensional hybrid Maxwell-Boltzmann model of shearth evolution 27 - Conferences publications P. Sarrailh L. Garrigues G. J. M. Hagelaar J. P. Boeuf G. Sandolache S. Rowe B. Jusselin
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
More informationNumerical Methods for Partial Differential Equations: an Overview.
Numerical Methods for Partial Differential Equations: an Overview math652_spring2009@colorstate PDEs are mathematical models of physical phenomena Heat conduction Wave motion PDEs are mathematical models
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationPRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING
PRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING Second Edition MICHAEL A. LIEBERMAN ALLAN J, LICHTENBERG WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC PUBLICATION CONTENTS PREFACE xrrii PREFACE
More informationFINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES
FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES SHRUTI JAIN B.Tech III Year, Electronics and Communication IIT Roorkee Tutors: Professor G. Biswas Professor S. Chakraborty ACKNOWLEDGMENTS I would like
More informationChapter 5 MAGNETIZED PLASMAS. 5.1 Introduction. 5.2 Diamagnetic current
Chapter 5 MAGNETIZED PLASMAS 5.1 Introduction We are now in a position to study the behaviour of plasma in a magnetic field. In the first instance we will re-examine particle diffusion and mobility with
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
More informationEffect of an axial magnetic field on a DC argon arc
Vol 17 No 2, February 2008 c 2008 Chin. Phys. Soc. 1674-1056/2008/17(02)/0649-06 Chinese Physics B and IOP Publishing Ltd Effect of an axial magnetic field on a DC argon arc Li Lin-Cun( ) and Xia Wei-Dong
More informationSolution of Time-dependent Boltzmann Equation
WDS'5 Proceedings of Contributed Papers, Part III, 613 619, 5. ISBN 8-8673-59- MATFYZPRESS Solution of Time-dependent Boltzmann Equation Z. Bonaventura, D. Trunec, and D. Nečas Masaryk University, Faculty
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationContents: 1) IEC and Helicon 2) What is HIIPER? 3) Analysis of Helicon 4) Coupling of the Helicon and the IEC 5) Conclusions 6) Acknowledgments
Contents: 1) IEC and Helicon 2) What is HIIPER? 3) Analysis of Helicon 4) Coupling of the Helicon and the IEC 5) Conclusions 6) Acknowledgments IEC:! IEC at UIUC modified into a space thruster.! IEC has
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationThe fast model for ionic wind simulation
Andrey Samusenko, Yury Stishkov, Polina Zhidkova The fast model for ionic wind simulation Research and Educational Center Electrophysics Saint Petersburg State University Faculty of Physics Ionic wind
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationA Magnetohydrodynamic study af a inductive MHD generator
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan A Magnetohydrodynamic study af a inductive MHD generator Augusto Montisci, Roberto Pintus University of Cagliari, Department of Electrical
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationMHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION
MHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION Marty Goldman University of Colorado Spring 2017 Physics 5150 Issues 2 How is MHD related to 2-fluid theory Level of MHD depends
More information8.1 Drift diffusion model
8.1 Drift diffusion model Advanced theory 1 Basic Semiconductor Equations The fundamentals of semiconductor physic are well described by tools of quantum mechanic. This point of view gives us a model of
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationModelling of low-temperature plasmas: kinetic and transport mechanisms. L.L. Alves
Modelling of low-temperature plasmas: kinetic and transport mechanisms L.L. Alves llalves@tecnico.ulisboa.pt Instituto de Plasmas e Fusão Nuclear Instituto Superior Técnico, Universidade de Lisboa Lisboa,
More informationStudy of DC Cylindrical Magnetron by Langmuir Probe
WDS'2 Proceedings of Contributed Papers, Part II, 76 8, 22. ISBN 978-737825 MATFYZPRESS Study of DC Cylindrical Magnetron by Langmuir Probe A. Kolpaková, P. Kudrna, and M. Tichý Charles University Prague,
More informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationFLOW OVER A GIVEN PROFILE IN A CHANNEL WITH DYNAMICAL EFFECTS
Proceedings of Czech Japanese Seminar in Applied Mathematics 2004 August 4-7, 2004, Czech Technical University in Prague http://geraldine.fjfi.cvut.cz pp. 63 72 FLOW OVER A GIVEN PROFILE IN A CHANNEL WITH
More informationLow Temperature Plasma Technology Laboratory
Low Temperature Plasma Technology Laboratory Equilibrium theory for plasma discharges of finite length Francis F. Chen and Davide Curreli LTP-6 June, Electrical Engineering Department Los Angeles, California
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationEquilibrium and transport in Tokamaks
Equilibrium and transport in Tokamaks Jacques Blum Laboratoire J.-A. Dieudonné, Université de Nice Sophia-Antipolis Parc Valrose 06108 Nice Cedex 02, France jblum@unice.fr 08 septembre 2008 Jacques Blum
More informationNumerical modeling of a cutting torch
Journal of Physics: Conference Series OPEN ACCESS Numerical modeling of a cutting torch To cite this article: B R Mancinelli et al 2014 J. Phys.: Conf. Ser. 511 012071 View the article online for updates
More information14. Energy transport.
Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation
More informationIntroduction. Chapter Plasma: definitions
Chapter 1 Introduction 1.1 Plasma: definitions A plasma is a quasi-neutral gas of charged and neutral particles which exhibits collective behaviour. An equivalent, alternative definition: A plasma is a
More informationPHYSICS OF HOT DENSE PLASMAS
Chapter 6 PHYSICS OF HOT DENSE PLASMAS 10 26 10 24 Solar Center Electron density (e/cm 3 ) 10 22 10 20 10 18 10 16 10 14 10 12 High pressure arcs Chromosphere Discharge plasmas Solar interior Nd (nω) laserproduced
More informationSW103: Lecture 2. Magnetohydrodynamics and MHD models
SW103: Lecture 2 Magnetohydrodynamics and MHD models Scale sizes in the Solar Terrestrial System: or why we use MagnetoHydroDynamics Sun-Earth distance = 1 Astronomical Unit (AU) 200 R Sun 20,000 R E 1
More information3D numerical simulation of the electric arc motion between bus-bar electrodes
3D numerical simulation of the electric arc motion between bus-bar electrodes M. Lisnyak 1 *, M. Chnani 2, A. Gautier 2, J-M. Bauchire 1 1 GREMI, UMR 7344, CNRS/Université d Orléans, 14 Rue d'issoudun,
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationIntroduction to numerical simulation of fluid flows
Introduction to numerical simulation of fluid flows Mónica de Mier Torrecilla Technical University of Munich Winterschool April 2004, St. Petersburg (Russia) 1 Introduction The central task in natural
More informationFlood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum
Archives of Hydro-Engineering and Environmental Mechanics Vol. 56 (29), No. 3 4, pp. 121 137 IBW PAN, ISSN 1231 3726 Flood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum Dariusz
More informationPhysique des plasmas radiofréquence Pascal Chabert
Physique des plasmas radiofréquence Pascal Chabert LPP, Ecole Polytechnique pascal.chabert@lpp.polytechnique.fr Planning trois cours : Lundi 30 Janvier: Rappels de physique des plasmas froids Lundi 6 Février:
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
More informationIntroduction to Heat and Mass Transfer. Week 8
Introduction to Heat and Mass Transfer Week 8 Next Topic Transient Conduction» Analytical Method Plane Wall Radial Systems Semi-infinite Solid Multidimensional Effects Analytical Method Lumped system analysis
More informationENHANCEMENT OF CONVECTIVE HEAT TRANSFER IN INTERNAL FLOWS USING AN ELECTRICALLY-INDUCED CORONA JET
ENHANCEMENT OF CONVECTIVE HEAT TRANSFER IN INTERNAL FLOWS USING AN ELECTRICALLY-INDUCED CORONA JET Reza Baghaei Lakeh Ph.D. Candidate PRESENTATION OUTLINE Corona Discharge Corona Wind and Ion-Drag Flows
More informationModélisation de sources plasma froid magnétisé
Modélisation de sources plasma froid magnétisé Gerjan Hagelaar Groupe de Recherche Energétique, Plasma & Hors Equilibre (GREPHE) Laboratoire Plasma et Conversion d Énergie (LAPLACE) Université Paul Sabatier,
More informationFINAL REPORT. DOE Grant DE-FG03-87ER13727
FINAL REPORT DOE Grant DE-FG03-87ER13727 Dynamics of Electronegative Plasmas for Materials Processing Allan J. Lichtenberg and Michael A. Lieberman Department of Electrical Engineering and Computer Sciences
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationCFD in Heat Transfer Equipment Professor Bengt Sunden Division of Heat Transfer Department of Energy Sciences Lund University
CFD in Heat Transfer Equipment Professor Bengt Sunden Division of Heat Transfer Department of Energy Sciences Lund University email: bengt.sunden@energy.lth.se CFD? CFD = Computational Fluid Dynamics;
More informationarxiv: v1 [physics.plasm-ph] 10 Nov 2014
arxiv:1411.2464v1 [physics.plasm-ph] 10 Nov 2014 Effects of fast atoms and energy-dependent secondary electron emission yields in PIC/MCC simulations of capacitively coupled plasmas A. Derzsi 1, I. Korolov
More informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationFeature-level Compensation & Control
Feature-level Compensation & Control 2 Plasma Eray Aydil, UCSB, Mike Lieberman, UCB and David Graves UCB Workshop November 19, 2003 Berkeley, CA 3 Feature Profile Evolution Simulation Eray S. Aydil University
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationFigure 1.1: Ionization and Recombination
Chapter 1 Introduction 1.1 What is a Plasma? 1.1.1 An ionized gas A plasma is a gas in which an important fraction of the atoms is ionized, so that the electrons and ions are separately free. When does
More informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More informationDiffusion during Plasma Formation
Chapter 6 Diffusion during Plasma Formation Interesting processes occur in the plasma formation stage of the Basil discharge. This early stage has particular interest because the highest plasma densities
More informationDraft Notes ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer
Draft Notes ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer Instructor: Jayathi Y. Murthy School of Mechanical Engineering Purdue University Spring 00 c 1998 J.Y. Murthy and S.R. Mathur.
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationA Plasma Torch Model. 1. Introduction
A Plasma Torch Model B. Chinè School of Materials Science and Engineering, Costa Rica Institute of Technology, Cartago, Costa Rica P.O. Box 159-7050, Cartago, Costa Rica, bchine@itcr.ac.cr Abstract: Plasma
More informationDiscretization of Convection Diffusion type equation
Discretization of Convection Diffusion type equation 10 th Indo German Winter Academy 2011 By, Rajesh Sridhar, Indian Institute of Technology Madras Guides: Prof. Vivek V. Buwa Prof. Suman Chakraborty
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More information2 GOVERNING EQUATIONS
2 GOVERNING EQUATIONS 9 2 GOVERNING EQUATIONS For completeness we will take a brief moment to review the governing equations for a turbulent uid. We will present them both in physical space coordinates
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationDedicating Finite Volume Method to Electromagnetic Plasma Modeling: Circuit Breaker Application
Dedicating Finite Volume Method to Electromagnetic Plasma Modeling: Circuit Breaker Application Loïc Rondot, Vincent Mazauric, Yves Delannoy, Gérard Meunier To cite this version: Loïc Rondot, Vincent Mazauric,
More informationLangmuir Probes as a Diagnostic to Study Plasma Parameter Dependancies, and Ion Acoustic Wave Propogation
Langmuir Probes as a Diagnostic to Study Plasma Parameter Dependancies, and Ion Acoustic Wave Propogation Kent Lee, Dean Henze, Patrick Smith, and Janet Chao University of San Diego (Dated: May 1, 2013)
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationHelicon Plasma Thruster Experiment Controlling Cross-Field Diffusion within a Magnetic Nozzle
Helicon Plasma Thruster Experiment Controlling Cross-Field Diffusion within a Magnetic Nozzle IEPC-2013-163 Presented at the 33rd International Electric Propulsion Conference, The George Washington University
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationSIMULATION MODEL OF INDUCTION HEATING IN COMSOL MULTIPHYSICS
Acta Electrotechnica et Informatica, Vol. 15, No. 1, 2015, 29 33, DOI: 10.15546/aeei-2015-0005 29 SIMULATION MODEL OF INDUCTION HEATING IN COMSOL MULTIPHYSICS Matúš OCILKA, Dobroslav KOVÁČ Department of
More informationPIC-MCC simulations for complex plasmas
GRADUATE SUMMER INSTITUTE "Complex Plasmas August 4, 008 PIC-MCC simulations for complex plasmas Irina Schweigert Institute of Theoretical and Applied Mechanics, SB RAS, Novosibirsk Outline GRADUATE SUMMER
More informationCONTROL OF UNIFORMITY IN CAPACITIVELY COUPLED PLASMAS CONSIDERING EDGE EFFECTS*
CONTROL OF UNIFORMITY IN CAPACITIVELY COUPLED PLASMAS CONSIDERING EDGE EFFECTS* Junqing Lu and Mark J. Kushner Department of Electrical and Computer Engineering at Urbana-Champaign mjk@uiuc.edu, jqlu@uiuc.edu
More informationImplicit kinetic relaxation schemes. Application to the plasma physic
Implicit kinetic relaxation schemes. Application to the plasma physic D. Coulette 5, E. Franck 12, P. Helluy 12, C. Courtes 2, L. Navoret 2, L. Mendoza 2, F. Drui 2 ABPDE II, Lille, August 2018 1 Inria
More informationNONLINEAR ELECTROMAGNETICS MODEL OF AN ASYMMETRICALLY DRIVEN CAPACITIVE DISCHARGE
NONLINEAR ELECTROMAGNETICS MODEL OF AN ASYMMETRICALLY DRIVEN CAPACITIVE DISCHARGE M.A. Lieberman Department of Electrical Engineering and Computer Sciences University of California Berkeley, CA 94720 Collaborators:
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More informationDepartment of Aerospace Engineering and Engineering Mechanics, The University of Texas at Austin, Austin, Texas 78712, USA
1 MAGNETIZED DIRECT CURRENT MICRODISCHARGE, I: EFFECT OF THE GAS PRESSURE Dmitry Levko and Laxminarayan L. Raja Department of Aerospace Engineering and Engineering Mechanics, The University of Texas at
More informationA Finite Element Model for Numerical Analysis of Sintering
A Finite Element Model for Numerical Analysis of Sintering DANIELA CÂRSTEA High-School Group of Railways, Craiova ION CÂRSTEA Department of Computer Engineering and Communication University of Craiova
More informationHIGH ACCURACY NUMERICAL METHODS FOR THE SOLUTION OF NON-LINEAR BOUNDARY VALUE PROBLEMS
ABSTRACT Of The Thesis Entitled HIGH ACCURACY NUMERICAL METHODS FOR THE SOLUTION OF NON-LINEAR BOUNDARY VALUE PROBLEMS Submitted To The University of Delhi In Partial Fulfillment For The Award of The Degree
More informationKinetic relaxation models for reacting gas mixtures
Kinetic relaxation models for reacting gas mixtures M. Groppi Department of Mathematics and Computer Science University of Parma - ITALY Main collaborators: Giampiero Spiga, Giuseppe Stracquadanio, Univ.
More informationMHD SIMULATIONS IN PLASMA PHYSICS
MHD SIMULATIONS IN PLASMA PHYSICS P. Jelínek 1,2, M. Bárta 3 1 University of South Bohemia, Department of Physics, Jeronýmova 10, 371 15 České Budějovice 2 Charles University, Faculty of Mathematics and
More informationDRIFT-DIFFUSION SYSTEMS: VARIATIONAL PRINCIPLES AND FIXED POINT MAPS FOR STEADY STATE SEMICONDUCTOR MODELS
DRIFT-DIFFUSION SYSTEMS: VARIATIONAL PRINCIPLES AND FIXED POINT MAPS FOR STEADY STATE SEMICONDUCTOR MODELS Joseph W. Jerome Department of Mathematics Northwestern University Evanston, IL 60208 Abstract
More informationBoundary Conditions for the Child Langmuir Sheath Model
IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 28, NO. 6, DECEMBER 2000 2207 Boundary Conditions for the Child Langmuir Sheath Model Mikhail S. Benilov Abstract A collision-free space-charge sheath formed by
More informationEXAMINATION QUESTION PAPER
Faculty of Science and Technology EXAMINATION QUESTION PAPER Exam in: Fys-2009 Introduction to Plasma Physics Date: 20161202 Time: 09.00-13.00 Place: Åsgårdvegen 9 Approved aids: Karl Rottmann: Matematisk
More information