The Free Molecular Flow Interface
In the free molecular flow regime, the mean free path of a gas is larger than the dimensions of the object being considered. Equivalently the Knudsen number is much greater than 1. Gas molecules therefore interact with surfaces more frequently than they interact with one another and the flow of gas is determined by collisions with the surfaces in the system.
The Free Molecular Flow (fmf) interface (), found under the Rarefied Flow branch () is intended for molecular flows interacting with objects that are moving slowly with respect to the speed of the molecules, such as vacuum systems. Diffuse reflection from all surfaces is assumed (this is reasonable in the majority of practical situations) with molecules from all directions effectively adsorbed onto the surface and subsequently reemitted according to Knudsen’s law (that is, with an intensity that varies as the cosine of the angle of emission to the normal to the surface).
When this physics interface is added, these default nodes are also added to the Model Builder: Molecular Flow, Wall, and Initial Values. Then, from the Physics toolbar, add other nodes that implement, for example, boundary conditions and volume forces. You can also right-click Free Molecular Flow to add physics from the context menu.
Settings
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics interface. Refer to such physics interface variables in expressions using the pattern <name>.<variable_name>. In order to distinguish between variables belonging to different physics interfaces, the name string must be unique. Only letters, numbers, and underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is fmf.
Integration Settings
Select an Integration methodHemicube (the default) or Direct area integration. These methods define how the incoming molecular flux, G, is computed at each element.
Hemicube
Hemicube is the default. The more sophisticated and general hemicube method uses a z-buffered projection on the sides of a hemicube (with generalizations to 2D and 1D) to account for shadowing effects. Think of it as rendering digital images of the geometry in five different directions (in 3D; in 2D only three directions are needed), and counting the pixels in each mesh element to evaluate its view factor.
Its accuracy can be influenced by setting the integration resolution of the virtual snapshots. The number of z-buffer pixels on each side of the 3D hemicube equals the specified resolution squared. Thus the time required to evaluate the flux increases quadratically with resolution. In 2D, the number of z-buffer pixels is proportional to the resolution property, and thus the time is, as well.
For an axisymmetric geometry, Gm and Famb must be evaluated in a corresponding 3D geometry obtained by revolving the 2D boundaries about the axis. COMSOL creates this virtual 3D geometry by revolving the 2D boundary mesh into a 3D mesh. You can control the resolution in the azimuthal direction by setting the number of azimuthal sectors, which is the same as the number of elements to a full revolution. Try to balance this number against the mesh resolution in the rz-plane.
Direct Area Integration
COMSOL Multiphysics evaluates the integrals directly, without considering which face elements are obstructed by others. This means that shadowing effects (that is, surface elements being obstructed in nonconvex cases) are not taken into account. Elements facing away from each other are, however, excluded from the integrals.
You can control the accuracy by specifying an integration order. Sharp angles and small gaps between surfaces can require a higher integration order for accuracy but also more time to evaluate the flux.
For Hemicube select an Integration resolution256 is the default.
For Direct area integration select an Integration order4 is the default.
For 2D axisymmetric components, also select the Integration sectors. To compute the integrals in an axisymmetric component COMSOL Multiphysics must construct a virtual geometry in 3 dimensions. The number of sectors in this geometry is determined by this setting.
Compute
In order to solve a free molecular flow it is always necessary to solve for the incident molecular flux, G for each species. The Pressure, Number density, and Heat flux on the surfaces can be optionally computed by performing additional integrations. The Pressure and Number density are computed by default.
Dependent Variables
Enter the Number of species. The default is 1, but it is possible to add more. It is often convenient to rename the species to something more descriptive for cases when multiple species exist, for example, H2 or SiH4.
The Incident molecular fluxes G (SI unit: 1/(m2s)) dependent variable is always solved for. It represents the incoming flux at the surface (the outgoing flux is available as fmf.J). Other dependent variables might not be solved for, depending on the settings in the Compute section (see above).
The Pressure p (SI unit: Pa) represents the total normal force acting on the surface per unit area (including contributions from both incoming and outgoing molecules).
The Number density n (SI unit: 1/(m3)), gives the total number density in the vicinity of the surface.
The Outward heat flux Q (SI unit: W/m2) is the heat flux out of the Free Molecular Flow domain as a result of heat adsorbed by the wall as it interacts with the gas.
For each of the dependent variables, the name can be changed in the corresponding field, but the name of fields and dependent variables must be unique within a model.
Discretization
Select the element order for the dependent variable. The default is constant, meaning the molecular flux, pressure and number density will be constant within each mesh element.
In the COMSOL Multiphysics Reference Manual see Table 2-4 for links to common sections and Table 2-5 to common feature nodes. You can also search for information: press F1 to open the Help window or Ctrl+F1 to open the Documentation window.
The following models are found at this Application Library path: Molecular_Flow_Module/Industrial_Applications/
Adsorption and Desorption of Water in a Load Lock Vacuum System: application name: water_adsorption_desorption
Differential Pumping: application name: differential_pumping
Evaporator: application name: evaporator
Molecular Flow in an Ion-Implant Vacuum System: application name: ion_implanter
Molecular Flow Through an S-Bend: application name: s_bend_benchmark
Molecular Flow Through a Microcapillary: application name: vacuum_capillary
Molecular Flow Through an RF Coupler: application name: rf_coupler
Outgassing Pipes: application name: outgassing_pipes
Rotating Plate in a Unidirectional Molecular Flow: application name: rotating_plate