# Einzel_Lensx - Comsol

Πολεοδομικά Έργα

15 Νοε 2013 (πριν από 4 χρόνια και 5 μήνες)

642 εμφανίσεις

Einzel Lens

rights reserved.

Introduction

An
Einzel

lens is an electrostatic device used for focusing
charged particle beams.

The beam focusing depends on the initial particle energy, the
voltage on the
Einzel

lens and the initial beam collimation
(initial radius and transverse velocity of the charged particles).

The problem is almost axisymmetric: the electrostatic field is
axisymmetric but the particle velocity transverse to the beam
axis is random, which spoils the symmetry. Thus, the model is
solved in 3D.

The model uses COMSOL
Multiphysics

with the Particle
Tracing Module.

3 identical cylinders arranged on same axis

Outer cylinders are grounded

Middle cylinder has non
-
zero voltage

This particular model uses cylinders, but other geometries are
possible.

Geometry

𝑉
=
0

𝑉
=
0

𝑉
=
𝑉
0

0

Focusing By Fringe Fields

Since the focusing is accomplished via the fringe fields, then it is important to model
this region accurately. This includes an accurate model of the cylinder edges (filleted)
and a dense mesh since the spatial variation of the field is highest near the edges.

Notes on COMSOL Model

The dimensions are listed in the Parameters table.

The geometry is constructed in a plane, which is then rotated
180
o

and mirrored.

Domains and boundaries are grouped in Selections (under
Definitions) to make the rest of the model clear.

The mesh is also constructed in the plane and then swept
180
o

in both directions.

The electrostatics problem is solved first. The particle
trajectories are computed in a separate study.

Initial Set Up

The first part is a 3D
e
lectrostatics calculation.

Parameters

Add a table of parameters under
Global Definitions.

Geometry

Work Plane

Draw a 2D cross
-
section of the geometry in a Work Plane. Later, the plane will
be revolved and mirrored to get the 3D geometry. The Work Plane is
𝑦
=
0
.

Geometry

Rectangle With Filleted Corners

Draw a Rectangle in the Work Plane, and
then Fillet all 4 corners.

Geometry

Array

Make copies of the filleted rectangle using a
1
×
3

array.

Geometry

Vacuum Chamber Rectangle

the vacuum chamber, including an
inner box using “Layers”.

Geometry

Work Plane

Add a point that will eventually
trace out a circle for the beam inlet.

The finished 2D geometry.

Geometry

Revolve

Revolve all objects in
the Work Plane by
180
o

Geometry

Mirror

Mirror all of the revolved objects about the
xz
-
plane, whose normal is
0
,
1
,
0
.
Keep the original set of objects as well!

Selections

Add Selections to clearly group domains and boundaries for later use.

Right click on Definitions and choose “Selections > Explicit”.

Selections

Cylinder 1

1)
Click on “Wireframe Rendering” in the Graphics window.

2)
Click on “Go to ZX View” in the Graphics window.

3)
Click on “Select Box” in the Graphics window and drag the box around the cylinder
on the left. Right
-
click to add the domain numbers to the selection list.

4)
Rename “Explicit 1” to “Cylinder 1”.

Step 1

Step
2

Step
3

Selections

Cylinder 2

1)
Right
-
click on Definitions and choose
“Selections > Explicit”.

2)
Click on “Select Box” in the Graphics window and drag the box around the cylinder
in the middle. Right
-
click to add the domain numbers to the selection list.

3)
Rename “Explicit 2” to “Cylinder 2”.

Step
2

Selections

Cylinder 3

1)
Right
-
click on Definitions and choose
“Selections > Explicit”.

2)
Click on “Select Box” in the Graphics window and drag the box around the cylinder
on the right. Right
-
click to add the domain numbers to the selection list.

3)
Rename “Explicit 3” to “Cylinder 3”.

Step
2

Selections

Vacuum Domains

1)
Right
-
click on Definitions and choose
“Selections >
Complement”.

2)
Click on “Add” (the “+” sign) and choose all 3 cylinders.

3)
Rename “Complement 1” to “Vacuum domain”.

Selections

Boundaries

Define Selections to group the boundaries of the cylinders individually and
collectively.

Right
-
click on “Cylinder 1” and choose “Duplicate”.

Change “Output Entities” from “Selected domains” to “Adjacent boundaries”.

Rename “Cylinder 1.2” to “Cylinder 1

boundaries.

Repeat for the other 2 cylinders.

Selections

All cylinder boundaries

Right
-
click on Definitions and choose “Selections > Union”.

Under “Geometric Entity Level”, change Level from Domain to Boundary.

Press the Add button and select all of the cylinder boundaries defined previously.

Rename “Union 1” to “All cylinder boundaries”.

Electrostatics

Right
-
click on Electrostatics and choose
“Sort by Space Dimension”.

Under “Domain Selection”, change “All
domains” to “Vacuum domains”. (The
potential inside each cylinder is equal to
the value on the surface, so it does not
need to be computed.)

Electrostatics

Vacuum Properties

We could define a material called “Vacuum”, but only the permittivity is needed so
just set the value in “Charge Conservation 1”.

Click on “Charge Conservation 1” and under “Electric Field”, change “Relative
permittivity” from “From material” to “User defined”. Leave the value as 1.

Electrostatics

Boundary Conditions

The default boundary condition is “Zero Charge”, which means


𝐷
=


𝜖𝛻𝑉
=
0
, that is, the surface charge density is zero.

The “Zero Charge” boundary condition is fine for the ends of the vacuum chamber.
For the other surfaces, the potential will be fixed.

Right
-
click on “Boundaries” and choose Ground. Select all of the surfaces on the
sides of the vacuum chamber. (2, 3, 6, 8, 11, 13, 59, 61, 63, 111, 112, 113)

Rename “Ground 1” to “Ground

vacuum chamber walls”.

Electrostatics

Outer Cylinders

Right
-
click on “Boundaries” and choose “Ground”.

Under “Boundary Selection”, change “Manual” to “Cylinder 1
-

boundaries”.

Rename “Ground 2” to “Ground

Cylinder 1”.

Repeat for Cylinder 3.

Electrostatics

Middle Cylinder

Right
-
click on “Boundaries” and choose “Electric Potential”.

Under “Boundary Selection”, change “Manual” to “Cylinder 2
-

boundaries”.

Change “Electric potential” to “V0”, defined in the Parameters table above.

Rename “Electric Potential 1”
to “Electric Potential
-

Cylinder 2
”.

Mesh

Edge

The particle trajectories will be near the cylindrical axis and so construct a dense
mesh there.

Right
-
click on “Mesh 1” and choose “More Operations > Edge”.

Add the edges on the axis to the Selection list. (122, 125, 128)

Right
-
click on “Edge 1” and choose “Distribution”.

Change “Number of elements” from 5 to 120.

Press “Build Selected”.

Mesh

Free Triangular

Add boundaries 85, 87 and 89 to the Selection list.

Right
-
click on “Free Triangular 1” and choose “Size”.

In “Size 1”, clear boundaries 85 and 89 from the Selection
list, so that only boundary 87 remains.

Under “Element Size”, change “Predefined” from
“Normal” to “Extremely fine”.

Select “Custom” and change “Maximum element size” to
0.01.

Press “Build Selected”.

Right
-
click on “Mesh 1” and choose “More Operations > Free Triangular”.

Mesh

Swept

Right
-
click on “Mesh 1” and choose “Swept”.

Add the domains, source faces and destination
faces as shown on the left.

Right
-
click on “Swept 1” and choose “Distribution”.

Change “Number of elements” from 5 to 10.

Press “Build Selected”.

Mesh

Swept

Right
-
click on “Swept 1” and choose “Duplicate”.

In the newly created “Swept 2”, change the
domains to 1,3,5, which are the other half of the
vacuum chamber.

Press “Build All”.

Right click on “Study 1” and choose Compute!

Data Sets

Cut Plane: y=0

Create 3 new data sets to plot a cut
-
away view of the equipotential surfaces near
the
Einzel

lens

Right
-
click on “Data Sets” and choose “Cut Plane”

Under “Plane Data”, change Plane to “
xz
-
planes”

Rename “Cut Plane 1” to “Cut Plane: y=0”

Data Sets

Solution 1

behind lenses only

Under “Data Sets”, right
-
click on “Solution 1” and choose
Duplicate

Rename “Solution 2” to “Solution 1

behind lenses only”

Right
-
click on
“Solution 1

behind lenses only
” and choose
Selections

Change “Geometric entity level” to “Domain” and add
domain 3, which is one of the half cylinder sections near the
Einzel

lens

Data Sets

Solution
3

cylinder surfaces

Right
-
click on “Data Sets” and choose “Solution”

Rename “Solution 3” to “
Solution 3

cylinder
surfaces”

Right
-
click on
“Solution 3

cylinder surfaces
” and choose Selection

Change “Geometric entity level” to Boundary

Change Selection to “All cylinder boundaries”

Add a “View” for the upcoming plot so that the perspective and zoom level can
remain fixed, even after more model steps are added

Right
-
click on “Show” and check “Advanced Results Options”, which has the effect
of adding a new node called “Views” under Results

Right
-
click on “Results > Views” and choose “View 3D”.

In “View 3D 3”, uncheck “Show grid” and “Show axis orientation”.

Plots

Equipotential Surfaces

Rename “Electric Potential” to “Equipotential surfaces near
Einzel

lens”

Under “Plot Settings”, change the View to “View 3D 3” and uncheck the box “Plot
data set edges”

Right
-
click on
“Equipotential surfaces near
Einzel

lens
” and choose Contour. Also
Isosurface

and Surface.

The plot settings for each are shown on the next slide.

Plot Settings

Equipotential Surfaces

Contour

Isosurface

Surface

Equipotential Surfaces

Rotate and zoom
-
in to obtain a plot similar to the one below

2D Plot of Fringe Field

Plot the fringe field in the cut plane y=0.

Right
-
click on Results and choose “2D Plot Group”.

Rename “2D Plot Group 2” to “Fringe Field”.

Change “Data set” to “Cut Plane: y=0”.

Right
-
click on “Fringe Field” and choose “Surface”. Also add “Contour” and
“Arrow Surface”.

The plot settings are shown on the next slide.

Plot Settings for Fringe Field

Surface

Contour

Arrow Surface

Fringe Field

Zoom
-
in to obtain a plot of the fringe field similar to the one
below.

To add “Charged Particle Tracing” to the model:

1)
Right
-
click on “Model 1” and choose “Add Physics”.

2)
Choose “AC/DC > Charged Particle Tracing”.

3)
Choose “Time Dependent” and uncheck “Electrostatics” below.

Particle Tracing

Vacuum Domains

Once “Charged Particle Tracing” has been added to the model, right
-
click on it and
choose “Sort By Space Dimension”.

Change “Domain Selection” to “Vacuum domains”.

Particle Tracing

Electric Force

By default, there are no forces on the particles. Right
-
click on “Charged Particle
Tracing > Domains” and choose “Electric Force”.

The electric force only needs to be computed near the
Einzel

lens. (It is negligible
elsewhere.) Add only domains 3 and 4 to the selection list.

The electric field is computed in “Electrostatics”, so change “Electric field” from
“User defined” to “Electric field (
es
/ccn1)”.

The gradient of the electrostatic potential is not continuous across mesh element
boundaries, so then neither is the electric field. Check the box beside “Use
piecewise polynomial recovery on field” to apply smoothing.

Particle Tracing

Coulomb Repulsion

The particle density in the beam is very low so the Coulomb repulsion
between particles is negligible.

If the particle density were high enough so that the Coulomb repulsion
between particles is significant, then this effect could be added by right
-
clicking on Domains and choosing “Particle
-
Particle Interaction”.

“Particle
-
Particle Interaction” is not added to the model.

Particle Tracing

Random Function

Define a random function in the model under “Global Definitions”. This will be
used in the initial particle velocity components transverse to the beam axis.

Change the “Distribution” of random numbers to Normal.

The number of arguments is 1, which is the seed for the random number
generator.

Initial Conditions of the Particles

The initial conditions are the position and the velocity.

The initial position is on the small circle at the end of
the vacuum chamber. Right
-
click on Boundaries and
choose Inlet.

Select boundaries 57, 58, 84 and 86.

Add the initial velocity. The x
-

and y
-
components are
small and random using seeds “x/

and “y/
”, respectively.

With “Initial position = Mesh based” and “Refinement
factor = 1”, there will be 1 particle released on each
mesh element on the boundary.

Electrons

Click on “Global > Particle Properties 1” and verify that the particles are electrons.

The charge number is
-
1 for the electron, that is, the charge on the particle is
(charge number) x (charge on proton).

The mass is set to “
me_const
”, where “_
const
” signifies a built
-
in physical
constant.

To see a list of all built
-
in physical constants, click on “Help > Documentation”,
search with “_
const
” and click on “Physical Constants”.

Particle Tracing

Time Dependent Study

Click on “Study 2 > Step 1: Time Dependent”
and specify the times for the particle tracing
calculation.

Note that “T” is the time it would take a
particle to travel from one end of the vacuum
chamber to the other in the absence of any
external force. The end time of the study is
5% higher in case the particles slow down
near the
Einzel

lens. (“T” is defined in
“Global Definitions > Parameters”.)

Verify that only “Charged Particle Tracing” will
be computed, not “Electrostatics”.

To use the electrostatics solution from the
previous study, check the box beside “Values
of variables not solved for”, change “Method”
to “Solution” and change “Study” to “Study 1,
Stationary”.

Right
-
click on “Study 2” and choose Compute!

Create a New Plot View

To create a new plot view, right
-
click on “Results > Views” and choose “View 3D”.

Clear all of the boxes.

The angle and zoom level for the following plots can be set to this particular
view, “View 3D 5”. If we want to change the view while doing something else in
the Model above, the view for the plots using “View 3D 5” will not be changed.

Plot

Particle Trajectories

Change the plot settings for “Particle Trajectories (
cpt
)” to those below.

The default is to plot the particle position using a point at the final time. Change
this to a line to trace out the trajectory.

Change the color so that it is the ratio of the particle kinetic energy, “
cpt.Ep
”, to the
initial kinetic energy, “E0”.

Particle Trajectories and Energies

Modify the perspective to obtain a plot similar to the one below.

Focusing depends on the initial position and velocity as well as the charge
-
to
-
mass
ratio and the voltage on the
Einzel

lens.

The initial speed of the electrons is
~
0
.
3𝑐
. The relativistic kinetic energy is
𝛾

𝑐
2
, where
𝛾
=
1
1

𝑣
2
𝑐
2

1
+
1
2
𝑣
2
𝑐
2
+
3
8
𝑣
4
𝑐
4
+

The first term in this expansion is

𝑐
2
, the rest mass
energy (which does not change in this example!), and the second term is
1
2

𝑣
2
, the
Galilean kinetic energy. The lowest order correction to the kinetic energy is
~
𝑣
4
𝑐
4
, so let’s
plot it!

Right
-
click on “Particle Trajectories (
cpt
)” and choose Duplicate.

Rename the new plot group and change “Color Expression” as shown below.

Relativistic effects < 1%, even at 0.3c!