Design of a Semi-analytical Method to Describe Plasma-Wall Interactions

Gonzalez Muñoz, Jorge (2013). Design of a Semi-analytical Method to Describe Plasma-Wall Interactions. Thesis (Master thesis), E.T.S.I. Aeronáuticos (UPM).


Title: Design of a Semi-analytical Method to Describe Plasma-Wall Interactions
  • Gonzalez Muñoz, Jorge
  • Donoso Vargas, Jose Manuel
Item Type: Thesis (Master thesis)
Masters title: Ciencia, Tecnología e Ingeniería Aeroespacial
Date: 27 September 2013
Faculty: E.T.S.I. Aeronáuticos (UPM)
Department: Física Aplicada a la Ingeniería Aeronáutica [hasta 2014]
Creative Commons Licenses: None

Full text

PDF - Requires a PDF viewer, such as GSview, Xpdf or Adobe Acrobat Reader
Download (6MB) | Preview


Electric probes are objects immersed in the plasma with sharp boundaries which collect of emit charged particles. Consequently, the nearby plasma evolves under abrupt imposed and/or naturally emerging conditions. There could be localized currents, different time scales for plasma species evolution, charge separation and absorbing-emitting walls. The traditional numerical schemes based on differences often transform these disparate boundary conditions into computational singularities. This is the case of models using advection-diffusion differential equations with source-sink terms (also called Fokker-Planck equations). These equations are used in both, fluid and kinetic descriptions, to obtain the distribution functions or the density for each plasma species close to the boundaries. We present a resolution method grounded on an integral advancing scheme by using approximate Green's functions, also called short-time propagators. All the integrals, as a path integration process, are numerically calculated, what states a robust grid-free computational integral method, which is unconditionally stable for any time step. Hence, the sharp boundary conditions, as the current emission from a wall, can be treated during the short-time regime providing solutions that works as if they were known for each time step analytically. The form of the propagator (typically a multivariate Gaussian) is not unique and it can be adjusted during the advancing scheme to preserve the conserved quantities of the problem. The effects of the electric or magnetic fields can be incorporated into the iterative algorithm. The method allows smooth transitions of the evolving solutions even when abrupt discontinuities are present. In this work it is proposed a procedure to incorporate, for the very first time, the boundary conditions in the numerical integral scheme. This numerical scheme is applied to model the plasma bulk interaction with a charge-emitting electrode, dealing with fluid diffusion equations combined with Poisson equation self-consistently. It has been checked the stability of this computational method under any number of iterations, even for advancing in time electrons and ions having different time scales. This work establishes the basis to deal in future work with problems related to plasma thrusters or emissive probes in electromagnetic fields.

More information

Item ID: 38761
DC Identifier:
OAI Identifier:
Deposited by: Sr. Jorge González Muñoz
Deposited on: 14 Dec 2015 15:28
Last Modified: 17 Oct 2016 15:05
  • Open Access
  • Open Access
  • Sherpa-Romeo
    Check whether the anglo-saxon journal in which you have published an article allows you to also publish it under open access.
  • Dulcinea
    Check whether the spanish journal in which you have published an article allows you to also publish it under open access.
  • Recolecta
  • e-ciencia
  • Observatorio I+D+i UPM
  • OpenCourseWare UPM