The use and application of FEMLAB . S.H.Lee and J.K.Lee. Plasma Application Modeling Lab. Department of Electronic and Electrical Engineering Pohang University of Science and Technology. 24. Apr. 2006. Plasma Application Modeling, POSTECH. What is FEMLAB?.
S.H.Lee and J.K.Lee
Plasma Application Modeling Lab.
Department of Electronic and Electrical Engineering
Pohang University of Science and Technology
24. Apr. 2006
What is FEMLAB?
*Reference: Manual of FEMLAB Software
1. Application of Chemical engineering Module
2. Application of Electromagnetics Module
3. Application of the Structural Mechanics Module
Classical PDE modes
PDE modes ( General, Coefficient, Weak)
2. Add used eqs. by using ‘add’ button.
1. Select eqs.
3. Multi-eqs. are displayed here.
FEMLAB Modeling Flow
In FEMLAB, use solid modeling or boundary modeling to create objects
in 1D, 2D, and 3D.
Physics and Mesh Menus
Solve and Postprocessing Menus
Magnetic Field of a Helmholtz Coil
1. Choose 3D, Electromagnetic Module, Quasi-statics mode in
2. After Application Mode Properties in Model Navigator is clicked,
the potential and Default element type are set to magnetic and
vector, respectively. Gauge fixing is off.
3. In the Options and setting menu, select the constant dialog box.
Define constant value (J0=1) in the constant dialog box.
Procedure of Simulation (2)
4. In the Geometry Modeling menu, open Work Plane Settings dialog
box, and default work plane is selected in x-y plane.
5. In the 2D plane, set axes and grid for drawing our simulation
geometry easily as follows,
6. Draw two rectangles by using Draw menu, then select these
rectangles . Click Revolve menu to revolve them in 3D.
In the 3D, add a sphere with radius of 1 and center of zero position.
It determines a calculation area.
Addition of a sphere with radius of 1 and center of zero position.
7. In the Physics Settings menu, select boundary conditions, and use
default for boundary conditions.
Select the Subdomain Settings, then fill in conductivity and external
current density in the Subdomain Settings dialog box.
Procedure of Simulation (4)
8. Element growth rate is set to 1.8 in Mesh Parameters dialog box
in Mesh Generation menu, and initialize it.
Result of a Helmholtz Coil
9. By using Postprocessing and Visualization menu, optimize your results.
Heated Rod in Cross Flow
1. Select 2D Fluid Dynamic, Incompressible Navier-Stokes, steady-state
analysis in the Model Navigator.
2. By using Draw menu, rectangle and half circle.
3. In the Subdomain Settings of Physics settings, enter v(t0)=0.5 in init tab.
Subdomain settings (physics tab)
Subdomain settings (init tab)
4. In the Boundary Settings dialog box, all boundaries are set to
Slip/Symmetry. Boundaries of 7 and 8 are no-slip.
Boundary Settings and Mesh Generation
5. Generate Mesh, and click Solve button.
Result of Velocity Flow
6. Add the Convection and Conduction mode in the Model Navigator.
7. In the Subdomain Settings, enter T(t0)=23 in the init tab of subdomain
of 1, 2.
Solving Convection and Conduction Eq.
8. In the Boundary Settings dialog box, all boundary conditions are thermal
insulation. 2 and 5 have the following boundary conditions.
9. In the Solver Manager, click Solver for tab, and select convection and
conduction. Click a Solve button.
Temperature Result of Heated Rod in Cross Flow
Steady-State 2D Axisymmetric Heat Transfer with Conduction
Boundary condition variations - General Heat Transfer
Electrostatic Potential Between Two Cylinder
This 3D model computes the potential field in vacuum around two cylinders, one with a potential of +1 V and the other with a potential of -1 V.
Porous Reactor with Injection Needle
Inlet species A
Inlet species C
Inlet species B
A + B C
Thin Layer Diffusion
D: diffusion coefficient(5e-5)
R: reaction rate(0)
Electromagnetic module(II) – Copper Plate
Electromagnetic module – Copper Plate
The plot shows the electric potential in copper plate.
The arrows show the current density.
The hole in the middle of geometry affects the potential and the current leading to a higher current density above and below the hole.
2D Steady-State Heat Transfer with Convection
Heat Transfer - 2D Steady-State Heat Transfer with Convection
556 elements is used as mesh.
2D symmetric Transient Heat Transfer
Heat Transfer - 2D symmetric Transient Heat Transfer
Semiconductor Diode Model
: symmetric boundary conditions
neumann boundary conditions
Semiconductor Diode Model
: Navier-Stokes equation
It is clear and intuitive that the magnitude of the velocity vector decreases as the cross-sectional area for the flow increases.