New Accumulators Boost Particle and Ray Tracing Functionality

With the release of COMSOL Multiphysics version 5.0, the Particle Tracing Module now includes a series of features calledAccumulators, which can be used to couple the results of a particle tracing simulation to other physics interfaces. The accumulated variables may represent any physical quantity and can be defined either within domains or on boundaries, making them extremely flexible. Here, I will explain the different types of accumulators and their applications in particle tracing and ray optics models.

What Are Accumulators?

In particle tracing and ray tracing simulations, we often need to use the particle or ray properties to change a variable that is defined on a set of domains or boundaries. For example, solid particles in a fluid might exert a significant force on the surrounding fluid, and they may also erode the surfaces they hit.

In previous blog posts, I’ve discussed two other cases in greater detail:divergence of an electron beam due to self-potentialandthermal deformation of lenses in a high-powered laser system. Each of these phenomena can be modeled using Accumulators or the specialized features that are derived from them.

An Accumulator is a physics feature that communicates information from particles or rays to the underlying finite element mesh. For each Accumulator feature in a model, a corresponding dependent variable, called anaccumulated variable, is declared. These accumulated variables can be defined either within a set of domains or on a set of boundaries, and they can represent any physical quantity, making them extremely flexible.

The Accumulator features can be added to any of the physics interfaces of theParticle Tracing Module. They can also be used in theGeometrical Opticsinterface, available with theRay Optics Module, and theRay Acousticsinterface, available with theAcoustics Module.

Depending on the physics interface, more specialized versions of the Accumulator may be available for computing specific types of physical quantities. For example, theParticle Tracing for Fluid Flowinterface includes a dedicatedErosionboundary condition that includes several built-in models for computing the rate of erosive wear on a surface.

The Accumulators can be divided into three broad categories, which function in the following ways:

1. Accumulators on boundaries increment a variable defined on a boundary element whenever a particle hits it.
2. Accumulators on domains project information from each particle to the mesh elements the particle passes through.
3. Nonlocal accumulators communicate information from a particle’s current position to the location where it was originally released.

We will now investigate each of these varieties in greater detail.

Accumulators on Boundaries

When particles or rays strike a surface, they can affect that surface in a wide variety of ways. For example, a laser can cause a boundary to heat up, sediment particles can erode their surroundings, and sputtering can occur when high-velocity ions strike a wafer in a process chamber. All of these effects require the same basic modeling procedure; we define a variable on the boundary and change its value when particles or rays interact with the boundary.

To begin, let’s consider a simple case in which we want to count the number of times a boundary is hit. We first define a variable, calledrpd, for example, which can have a distinct value in every boundary mesh element. Initially, this variable is set to zero in all elements. Every time a particle hits a mesh element on this boundary, we would like to increment the variable on that element by 1.

The values of the accumulated variable on the boundary elements (illustrated as triangles) after one collision are shown below:

To implement this in COMSOL Multiphysics, we first set up the particle tracing model, then add a “Wall” node to the boundary for which we want to count collisions. In this case, let’s specify that particles are reflected at this surface by selecting the Bounce wall condition. We then add the Accumulator node as a subnode to this Wall.

The settings shown in the following screenshot cause the accumulated variable (calledrpb) to be incremented by 1 (the expression in the Source edit field) every time a particle hits the wall.

I have created an animation that demonstrates how the number of collisions with each boundary element is counted over the course of the study. Check it out:

By changing the expression in the Source edit field, it is possible to increment the accumulated variable using any combination of variables that exist on the particle and on the boundary. For example, the accumulated variable may increase by a different amount based on the velocity or mass of incoming particles. The dependent variable need not be dimensionless. In fact, it can represent any physical quantity.

In addition to the generic Accumulator subnode — which can represent anything — dedicated accumulator-based features are available in the different physics interfaces, including the following:

• In theCharged Particle Tracingphysics interface:
  • Etch(Use this to model physical sputtering of a surface by energetic ions.)
  • Current Density
  • Heat Source
  • Surface Charge Density
• In theParticle Tracing for Fluid Flowphysics interface:
  • Erosion(For computing the total mass removed from the surface or the rate of erosive wear.)
  • Mass Deposition
  • Boundary Load
  • Mass Flux
• In theGeometrical Opticsphysics interface:
  • Deposited Ray Power(For computing a boundary heat source using the power of incident rays.)

Accumulators on Domains

We may also want to transfer information from particles to all of the mesh elements they pass through, not just the boundary elements they touch. We can do so by adding an Accumulator node to the physics interface directly, instead of adding it as a subnode to a Wall or other boundary condition.

For example, we can use an Accumulator to reconstruct the number density of particles within a domain. This technique is used in a benchmark model of freemolecular flow through an s-bendin which theFree Molecular Flowinterface is used to compute the number density of molecules in a rarefied gas.

Here is the geometry of the s-bend:

The settings window for the Accumulator is shown below.

The expression in the Source edit field is a bit more complicated than in the previous case. The source termRis defined as

(1)

whereJ_{\textrm{in}}(SI unit:1/(m^2 s)) is the molecular flux at the inlet,L(SI unit:m) is the length of the inlet, andN_{p}(dimensionless) is the number of model particles.

Physically, we can interpretRas the number of real molecules per unit time, per unit length in the out-of-plane direction, that are represented by each model particle. Because this source term acts on the time derivative of the accumulated variable, each particle leaves behind a “trail” in the mesh elements it passes through, which contributes to the number density in those elements.

I have created a second animation in which the number density of molecules is computed using the Accumulator (bottom) and the result is compared to the result of theFree Molecular Flowinterface (top). Here it is:

We do see some noise in the particle tracing solution because each particle can only make a uniform contribution to the mesh element it is currently in. However, when the number of particles is large compared to the number of mesh elements, it is still possible to obtain an accurate solution.

In addition to the generic Accumulator node, which can represent anything, dedicated accumulator-based features are available in the different physics interfaces, including the following:

• In theCharged Particle Tracingphysics interface:
  • Particle-Field Interaction computes the charge density of particles, which can then be used as a source term to compute the self-potential of a beam of ions or electrons. It is also possible to compute the current density, which can create a significant magnetic field if the beam is relativistic.
• In theParticle Tracing for Fluid Flowphysics interface:
  • Fluid-Particle Interaction computes the body load exerted by particles on the surrounding fluid.
• In theGeometrical Opticsphysics interface:
  • Deposited Ray Power generates a heat source term based on the amount of power absorbed by the medium if rays propagate through an absorbing medium.

Nonlocal Accumulators

The third variety of Accumulator is a bit more advanced than the previous two. ANonlocal Accumulatoris used to communicate information from a particle’s current position to the initial position from which it was released. The Nonlocal Accumulator can be added to an “Inlet” node, causing it to declare an accumulated variable on the mesh elements on the Inlet boundary.

The Nonlocal Accumulator can be used in some advanced models of surface-to-surface radiation. In many cases, theSurface-to-Surface Radiationphysics interface (available with theHeat Transfer Module) can be used to efficiently and accurately model radiative heat transfer. However, theSurface-to-Surface Radiationinterface relies on the assumption that all surfaces reflect radiation diffusely. That is, the direction of reflected radiation is completely independent of the direction of incident radiation. It cannot be used, for example, if some of the radiation undergoes specular reflection at smooth, polished, metallic surfaces.

One approach to modeling radiative heat transfer with a combination of specular and diffuse radiation is to use theMathematical Particle Tracinginterface, as demonstrated in the example ofmixed diffuse and specular reflection between two parallel plates.

The incident heat flux on each plate is computed by releasing particles from the plate surface, querying the temperature of each surface the particles hit, and communicating this information back to the point at which the particles are initially released. The below image shows the temperature distribution between the two plates, where the top plate is heated by an external Gaussian source.

Conclusions and Next Steps

We have seen that Accumulators can be used to model interactions between particles or rays and any field that is defined on the surrounding domains of boundaries. The accumulated variables can represent any physical quantity. The Accumulator is the basic building block that allows for sophisticated one-way or two-way coupling between a particle- or ray-based physics interface and any of the other products in theCOMSOL product suite.

The Accumulators and related physics features have too many settings and applications to discuss in detail in a single blog post. To learn more about the many options available, please refer to the User’s Guide for the Particle Tracing Module (for particle tracing physics interfaces), the Ray Optics Module (for theGeometrical Opticsinterface), or the Acoustics Module (for theRay Acousticsinterface).

If you are interested in learning more about any of these products, pleasecontact us.