Molecular Flow Module
Molecular Flow Module
For Modeling Low Pressure Gas Flow in Vacuum Systems
In an ion implanter, the average number density of outgassing molecules along the beam path is used as a figure of merit to evaluate the design. It must be computed as a function of wafer angle, with rotation about one axis.
Accurate Modeling of Low Pressure, Low Velocity Gas Flows
The Molecular Flow Module is designed to offer previously unavailable simulation capabilities for the accurate modeling of low pressure, low velocity gas flows in complex geometries. It is ideal for the simulation of vacuum systems including those used in semiconductor processing, particle accelerators and mass spectrometers. Small channel applications (e.g. shale gas exploration and flow in nanoporous materials) can also be addressed.
The Molecular Flow Module uses a fast angular coefficient method to simulate steady-state free molecular flows. You can model isothermal and nonisothermal molecular flows, and automatically calculate the heat flux contribution from the gas molecules. The discrete velocity method is also included in the module for the simulation of transitional flows.
Two Methods for Modeling Free Molecular and Transitional Flows
The Molecular Flow Module offers two alternatives to these methods, allowing you to solve for low velocity and low pressure flows in a manageable and accurate fashion. Two specific physics interfaces, configured to receive model inputs via the graphical user interface (GUI) to fully specify a set of equations, are available:
Free Molecular Flows
The Free Molecular Flow interface uses the angular coefficient method to model flows with Knudsen numbers that are greater than ten. This physics interface avoids solving the physics in the volumes of the modeled geometries, and requires meshing only of the surfaces. Completely diffuse scattering (total accommodation) and emission are assumed at all surfaces in the geometry, and flow is computed by integrating the flux arriving at a surface from all other surfaces in its line-of-sight. This means that the dependent variables exist only on the surfaces of the geometry, and the solution process is much faster than the DSMC method. Furthermore, it is not subject to statistical scatter. Number densities are reconstructed using a method included in the Free Molecular Flow interface.
The Transitional Flow interface solves the Boltzmann BGK equation by employing a modified form of the Lattice Boltzmann/Discrete Velocity method to solve transitional flows. Unlike the DSMC method, the solutions are not subject to statistical noise. Diffuse reflection of gas molecules is also assumed at all surfaces, with molecules from all directions effectively adsorbed onto the surface and subsequently re-emitted according to Knudsen’s law. In this interface, the model geometry is meshed to discretize the physical space, and a velocity quadrature is chosen, which provides a set of dependent variables that represent a mesh in velocity space. Both the mesh and the quadrature can be independently adjusted to ensure the problem is resolved in both physical and velocity space.
Optimized Methods for Fast and Accurate Simulations
Gases at low pressures cannot be modeled using conventional computational fluid dynamics tools. That is due to the fact that kinetic effects become important as the mean free path of the gas molecules becomes comparable to the length scale of the flow. Flow regimes are categorized quantitatively via the Knudsen number (Kn), which represents the ratio of the molecular mean free path to the flow geometry size for gases:
|Flow type||Knudsen Number|
|Free molecular flow||Kn>10|
While the Microfluidics Module is used for modeling slip and continuum flows, the Molecular Flow Module is designed for accurately simulating flows in the free molecular flow and transitional flow regimes. Historically, flows in this regime have been modeled by the direct simulation Monte Carlo (DSMC) method. This computes the trajectories of large numbers of randomized particles through the system, but introduces statistical noise to the modeling process. For low velocity flows, such as those encountered in vacuum systems, the noise introduced by DSMC renders the simulations unfeasible. COMSOL uses alternative approaches: employing a discrete velocity method for transitional flows (using a Lattice Boltzmann velocity quadrature) and the angular coefficient method for molecular flows.
Differentially pumped vacuum systems use a small orifice or tube to connect two parts of a vacuum system that are at very different pressures. Such systems are necessary when processes run at higher pressures and are monitored by detectors that require UHV for operation. In this model, gas flow through a narrow tube and into a high vacuum chamber ...
Molecular Flow in an Ion Implant Vacuum System
This example shows how to model an ion implantation system using the Molecular Flow interface available in the Microfluidics Module. In ion implantation, outgassing molecules interact with the ion beam to produce undesirable species. The average number density of outgassing molecules along the beam path is used as a figure of merit to evaluate the ...
Molecular Flow Through a Microcapillary
Computing molecular flows in arbitrary geometries produces complex integral equations that are very difficult to compute analytically. Analytic solutions are, therefore, only available for simple geometries. One of the earliest problems solved was that of gas flow through tubes of arbitrary length, which was first treated correctly by Clausing. ...
Rotating Plate in a Unidirectional Molecular Flow
This model computes the particle flux, number density and pressure on the surface of a plate that rotates in a highly directional molecular flow. The results obtained are compared with those from other, approximate, techniques for computing molecular flows.
Adsorption and Desorption of Water in a Load Lock Vacuum System
This model shows how to simulate the time-dependent adsorption and desorption of water in a vacuum system at low pressures. The water is introduced into the system when a gate valve to a load lock is opened and the subsequent migration and pumping of the water is modeled.
This benchmark model computes the pressure in a system of outgassing pipes with a high aspect ratio. The results are compared with a 1D simulation and a Monte-Carlo simulation of the same system from the literature.
Molecular Flow in an RF Coupler
This model computes the transmission probability through an RF coupler using both the angular coefficient method available in the Free Molecular Flow interface and a Monte Carlo method using the Mathematical Particle Tracing interface. The computed transmission probability determined by the two methods is in excellent agreement with less than a 1% ...