The Boltzmann Equation, Two-Term Approximation Interface
In a previous blog post, we introduced readers to different kinds of electron energy distribution functions (EEDFs) and their importance in plasma modeling. Today, we focus our attention on the Boltzmann Equation, Two-Term Approximation interface, demonstrating its use with an example from our Application Library.
Editor’s note: The original version of this blog post was published on April 8, 2015. It has since been updated to reflect new functionality in version 6.0 of the COMSOL Multiphysics® software.
An Introduction to the Boltzmann Equation, Two-Term Approximation Interface
In a plasma model, the electron energy distribution function as well as the transport properties of the electrons (e.g., electron mobility) are needed. For the simplest cases, a Maxwellian EEDF and a constant value for the electron mobility can be used. The other transport properties are then computed in COMSOL Multiphysics using Einstein’s relation. However, in some cases, it might be advantageous to use an EEDF obtained from the solution of the Boltzmann equation and define the electron transport properties as a function of the mean electron energy. But how do we obtain this data?
The answer: The Boltzmann Equation, Two-Term Approximation interface in COMSOL Multiphysics. Examples illustrating the use of this interface are available in the Application Library, one of which is the Argon Boltzmann Analysis model. To compute the Boltzmann equation in the two-term approximation, parameters such as the ionization degrees of the plasma are required. These parameters are not known a priori. Thus, the procedure is an iterative process.
The process begins with making an initial guess for the parameters and solving the Boltzmann equation. Then, the EEDF and transport properties are imported, if needed, to your plasma model. Finally, the plasma model is computed and the Boltzmann equation is re-solved with the new parameters from the plasma model. You can continue to repeat these steps until convergence is reached.
Let’s walk through the steps of creating, exporting, and importing the data into the plasma model.
Electron Energy Distribution Function and Electron Transport Properties
Create the Data from the Boltzmann Equation, Two-Term Approximation Interface
The first step is to create the data by solving the Boltzmann equation in the two-term approximation. The figure below shows a screenshot from the Boltzmann Equation, Two-Term Approximation interface used for this step. You need to define a constant Maximum energy value for the electron energy. In our example, this is set to
Emax = 100 V. In addition, you need a Mean Energies study to compute the EEDF for a range of mean electron energies.
The figure below shows computed EEDFs in an argon plasma for several mean electron energies. The plasma features a gas temperature of 400 K, an electron density of 1018 1/m³, an ionization degree of 10-6, and a mole fraction of excited argon of 0.01 %. It can be observed that the EEDF is a function of the electron energy.
Computed EEDFs in an argon plasma.
The next figure depicts the corresponding reduced electron transport properties computed with the Boltzmann Equation, Two-Term Approximation interface. The data is a function of the mean electron energy.
Reduced electron transport properties in an argon plasma.
Exporting the Data from the Boltzmann Equation, Two-Term Approximation Interface
The EEDF has to be imported to the plasma model as a spreadsheet consisting of three rows. The first row (x data) has to be the electron energy (eV), while the second row (y data) has to be the mean electron energy (eV). Meanwhile, the third row must contain the value for the distribution function (
eV^(-3/2)). Finally, you need to export a 2D plot that looks like the following image.
A 2D plot of the computed EEDF. Here, the x-axis indicates the electron energy and the y-axis represents the mean electron energy. The color illustrates the value of the distribution function.
Note on units: In the Plasma Module, the electron energy and electron mean energy appear with units of V but internally are treated as eV. As such, in this context, V should be read as eV. The extra dimension that represents the electron energy and on which the EEDF is solved appears with units of meters but is treated internally as eV.
To export the EEDF in the desired format, use the Parametric Extrusion dataset. The Parametric Extrusion extends the dataset with the parameter (the mean energies in this case). Right-click Datasets and choose Parametric Extrusion 1D.
The Boltzmann Equation, Two-Term Approximation interface is a 0D interface and uses an extra dimension to represent the electron energy on a 1D axis. As the extra dimension is normalized to the Maximum energy value, it has to be manually scaled with the Maximum energy value before exporting the data. To do this, we can use a Transformation dataset. With a Transformation dataset, you can scale, rotate, and move datasets. So right-click Datasets and choose Transformation 2D. Set the x-scale to the maximum energy (
Emax in this case).
Then, right-click the Transformation 2D data and choose Add Data to Export. Enter
be.f in the Expression window. Select a filename and click Export.
The transport properties can be exported in a straightforward manner from the respective 1D plot. Right-click the Global node in the 1D plot and choose Add Plot Data to Export. Select a file name and click Export.
Importing the Data to the Plasma Model
To import the EEDF to the plasma model, create an Interpolation function. Select File as the data source and enter 2 in the Number of arguments field. Click Browse and import the file.
Once you have done this, you can choose the interpolation function as an EEDF in the Electron Energy Distribution Function settings in the main node of the plasma model. This is illustrated in the screenshot below.
The functions for the transport properties, such as electron mobility and diffusivity, can also be imported as an interpolation function. Here, however, the number of arguments is 1. In the plasma model node, you can use this interpolation function by entering
int2(plas.ebar). In this case,
int2 is the name of the function,
plas is the tag of the interface, and
ebar is the mean electron energy.
- COMSOL Now
- Fluid & Heat
- Structural & Acoustics
- Today in Science