COMSOL 6.4 - Plasma Pulsed Arc

PDF

Plasma Pulsed Arc

Thermal plasmas have nowadays a large range of industrial applications including cutting, welding, spraying, waste destruction, and surface treatment. Thermal plasmas are assumed to be under partial to complete local thermodynamic equilibrium (LTE) conditions. Under LTE, the plasma can be considered a conductive fluid mixture and therefore, be modeled using the magnetohydrodynamics (MHD) equations. This model shows how to use the Equilibrium Discharges, In-Plane Currents interface (available in 2D and 2D axisymmetric) interface to simulate the plasma generated in a pulsed arc.

This application requires the Plasma Module and AC/DC Module.

Model Definition

This model is based on the work presented in Ref. 1. In Ref. 1 The authors develop a complex model that includes the description of the weld pool under the action of a pulsed arc. In this work it is only simulated the plasma and the transfer of heat and currents in the metals neglecting the weld pool. This model starts by opening a model of a DC arc which solution are going to be used as initial conditions to the time-dependent problem.

The applied current consists of a pulse with a frequency of 1 Hz, a peak current of 160 A, a floor current of 80 A, and a duty cycle of 0.5. The current source is set at the cathode and the bottom plate is grounded. In the 5 mm gap between the electrodes an argon plasma arc is created that heats the metal electrodes and surrounding gas. A shielding flow is added along the cathode.

The temperature-dependent physical properties of argon are loaded from the material library under Equilibrium Discharge. The temperature range of the physical properties span from 500 K to 25,000 K. A minimum electric conductivity of 1 S/m is used for numerical stability reasons. Another important aspect to keep in mind is that the model used is not valid to describe the plasma sheath region since in this regions there is charge separation and deviations from equilibrium. From the practical point of view, having a fine resolution in the plasma-electrode region causes numerical instabilities (and does not bring a better description of the physics). To make the model more stable use a mesh coarse enough so that the plasma sheath is averaged out.

Results and Discussion

Figure 1 and Figure 2 show the fluid velocity and temperature for instants in the current peak and floor of 160 A and 80 A at 5 s. As expected, at the peak value the fluid velocity and temperature are higher. In the results shown the temperature of the subtract still has not reached a steady state that is reached at about 50 s. The final time of computation was set to 5s for practical reasons. The weld pool (not solved here) also has significant changes during the pulse different phases as shown in Ref. 1.

Figure 1: Plot of the velocity magnitude of the fluid at the current peak of 160 A (left) and current floor of 80 A (right).

Figure 2: Plot of the plasma temperature at the current peak of 160 A (left) and current floor of 80 A (right).

1. A. Traidia, Multiphysics Modeling and Numerical Simulation of GTA Weld Pools, École Polytechnique, PhD Thesis, Paris, France, 2011.

Plasma_Module/Equilibrium_Discharges/plasma_pulsed_arc

Modeling Instructions

This model uses an existing model from a DC arc and sets an AC excitation instead. The solution of the DC arc model are used as initial conditions to the time-dependent model.

Application Libraries

1

From the menu, choose .

2

In the window, select > > in the tree.

3

button in the toolbar.

Add a parameter to define the floor of the current pulse.

Global Definitions

1

In the window, under click .

2

In the window for , locate the section.

3

In the table, enter the following settings:


Name	Expression	Value	Description
delta	80[A]	80 A

Add a square wave function to define the pulsed excitation.

Component 1 (comp1)

In the window, expand the node.

Waveform 1 (wv1)

1

In the window, expand the > node.

2

Right-click and choose > .

3

In the window for , type current in the text field.

4

Locate the section. From the list, choose .

5

In the text field, type 0.05.

6

In the T text field, type 1.

7

In the φ text field, type 180.

8

In the A text field, type delta/2.

Duplicate the existent Normal Current Density feature and set a pulsed current. The Normal Current Density features are disabled at the study level when needed.

Magnetic and Electric Fields (mef)

In the window, expand the > node.

Normal Current Density 2

1

In the window, expand the > > node.

2

Right-click > > > and choose .

3

In the window for , locate the section.

4

In the Jn text field, type -(I0+delta/2+current(t[1/s]))/(pi*(1.6[mm])^2).

Disable the pulsed excitation for the stationary study. This way, if this study needs to be solved later the correct current feature is used.

Step 1: Stationary

1

In the window, expand the node, then click .

2

In the window for , locate the section.

3

Select the checkbox.

4

In the tree, select > > > .

5

Electric conductivity, Magnetic Flux, Temperature, Temperature 3D, Velocity

1

In the window, under , Ctrl-click to select , , , , and .

2

Right-click and choose .

In the window for , type DC Arc in the text field.

1

In the toolbar, click

to open the window.

2

Go to the window.

3

Find the subsection. In the tree, select > .

4

Click the button in the window toolbar.

5

In the toolbar, click

to close the window.

Step 1: Time Dependent

Prepare the time-dependent study. Here, a few changes are necessary: disable the DC excitation, use previous solutions as initial conditions, use a fully coupled solver.

1

In the window for , locate the section.

2

Select the checkbox.

3

In the tree, select > > > .

4

5

Locate the section. In the text field, type range(0,0.1,5).

6

In the toolbar, click

.

Electric Field (mef), Magnetic Flux Density (mef), Magnetic Flux Density, Revolved Geometry (mef), Pressure (spf), Temperature (ht), Velocity (spf)

1

In the window, under , Ctrl-click to select , , , , , and .

2

Right-click and choose .

In the window for , type Pulsed Arc in the text field.

Temperature Pulsed

1

In the toolbar, click

.

2

In the window for , type Temperature Pulsed in the text field.

3

Locate the section. From the list, choose .

1

Right-click and choose .

2

In the window for , locate the section.

3

In the text field, type T.

4

Locate the section. From the list, choose .

Step 1: Time Dependent

1

In the window, under click .

2

In the window for , click to expand the section.

3

Find the subsection. From the list, choose .

4

From the list, choose .

5

From the list, choose .

Solver Configurations

In the window, expand the > node.

Solution 2 (sol2)

1

In the window, expand the > > node.

2

Right-click and choose .

3

In the window for , click to expand the section.

4

In the text field, type 0.8.

5

From the list, choose .

Step 1: Time Dependent

1

In the window, under click .

2

In the window for , click to expand the section.

3

Select the checkbox.

4

In the table, enter the following settings:


Plot group	Plot window
Temperature Pulsed	Graphics

5

From the list, choose .

6

In the toolbar, click

.

Magnetic Flux Density (mef)

Prepare plots to show the velocity and the temperature for the peak and floor values of the current side by side.

1

In the window, click .

2

In the window for , click to expand the section.

3

From the list, choose .

4

Click to expand the section. From the list, choose .

1

In the window, expand the node, then click .

2

In the window for , locate the section.

3

From the list, choose .

4

From the list, choose .

5

Locate the section. From the list, choose .

1

Right-click and choose .

2

In the window for , locate the section.

3

From the list, choose .

4

Click to expand the section. From the list, choose .

5

In the toolbar, click

.

6

button in the toolbar.

1

In the window, click .

2

In the window for , locate the section.

3

Select the checkbox.

Temperature Pulsed

1

In the window, click .

2

In the window for , locate the section.

3

From the list, choose .

4

Locate the section. Select the checkbox.

5

Locate the section. From the list, choose .

1

In the window, click .

2

In the window for , locate the section.

3

From the list, choose .

4

From the list, choose .

5

Locate the section. From the list, choose .

1

Right-click > > > and choose .

2

In the window for , locate the section.

3

From the list, choose .

4

Locate the section. From the list, choose .

5

In the toolbar, click

.

Temperature Pulsed

Plot the temporal evolution of the applied current and the temperature of the bottom plate.

Current vs. Time

1

In the toolbar, click

.

2

In the window for , type Current vs. Time in the text field.

3

Locate the section. From the list, choose .

4

Click to expand the section. From the list, choose .

1

Right-click and choose .

2

Select Point 4 only.

3

In the window for , locate the section.

4

In the text field, type -mef.mi1.ncd2.nJ*(pi*(1.6[mm])^2).

5

In the toolbar, click

.

Temperature vs. Time

1

In the toolbar, click

.

2

In the window for , type Temperature vs. Time in the text field.

3

Locate the section. From the list, choose .

4

Locate the section. From the list, choose .

1

Right-click and choose .

2

Select Point 1 only.

3

In the window for , locate the section.

4

In the text field, type T.

5

In the toolbar, click

.