CN112989674B - Simulation analysis method and simulation analysis device for arc plasma characteristics - Google Patents
Simulation analysis method and simulation analysis device for arc plasma characteristics Download PDFInfo
- Publication number
- CN112989674B CN112989674B CN202110421377.4A CN202110421377A CN112989674B CN 112989674 B CN112989674 B CN 112989674B CN 202110421377 A CN202110421377 A CN 202110421377A CN 112989674 B CN112989674 B CN 112989674B
- Authority
- CN
- China
- Prior art keywords
- setting
- current
- arc
- geometric model
- arc plasma
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
The invention discloses a simulation analysis method and a simulation analysis device for arc plasma characteristics, and relates to the simulation analysis of arc plasma in the metallurgical industry. The specific implementation mode comprises the following steps: establishing two-dimensional axisymmetric spatial dimensions, and selecting direct-current coupling discharge; establishing a geometric model of arc plasma discharge; setting material properties of the electrode part and the arc part respectively; respectively arranging a magnetic field part, a fluid heat transfer part, a laminar flow part, a multi-physical field part and a current part; carrying out mesh division on the geometric model; and configuring time units, time step lengths, time steps and tolerance of a solver, and performing solving calculation to generate temperature distribution maps of the arc plasma corresponding to different current values and axial current density distribution maps at different current values. The method can simulate the characteristic distribution of the arc plasma under different current values, realizes the research on the influence of the current on the characteristics of the arc plasma, and has the advantages of convenient operation, high efficiency and accuracy.
Description
Technical Field
The invention relates to simulation analysis of arc plasma in the metallurgical industry, in particular to simulation of arc plasma characteristics by using Comsol Multiphysics.
Background
The arc plasma has the characteristics of high temperature, concentrated energy and the like, and is widely applied to the industrial fields of metallurgy, chemical industry, waste treatment and the like. The formation process of arc plasma involves the interaction between complex electric field, magnetic field, flow field and temperature field, and is a complex process with multiple physical field changes.
At present, a plurality of restrictive factors exist in the practical measurement process, for example, the central temperature of the arc plasma can reach about 20000K, which makes the temperature measurement of the arc plasma have certain difficulties. When the test means is limited, the real-time monitoring of the characteristic parameters of the arc plasma is limited, and the characteristics of the arc plasma under different current values are difficult to study.
Therefore, in view of the above disadvantages, it is desirable to provide an accurate and effective numerical model for studying arc plasma.
Disclosure of Invention
The invention aims to solve the technical problem that real-time monitoring of arc plasma characteristic parameters is limited, and the characteristics of analysis arc plasma at different current values are difficult to study. Aiming at the defects in the prior art, a simulation analysis method and a simulation analysis device for arc plasma characteristics are provided.
In order to solve the above technical problem, according to an aspect of an embodiment of the present invention, there is provided a simulation analysis method of arc plasma characteristics, including the steps of:
(1) establishing two-dimensional axisymmetric spatial dimensions, and selecting direct-current coupling discharge;
(2) establishing a geometric model of arc plasma discharge; wherein the geometric model comprises an electrode portion and an arc portion, and the computational domain of the geometric model comprises a geometric region and an arc computational domain;
(3) setting material properties of the electrode part and the arc part, respectively;
(4) providing the magnetic field portion, the fluid heat transfer portion, the laminar flow portion, the multi-physical field portion, and the current portion, respectively;
wherein the content of the first and second substances,
setting the magnetic field part to have relative magnetic permeability, electric conductivity and relative dielectric constant from materials, setting an initial value of magnetic vector potential and a divergence condition variable scaling value, and adding vector magnetic potential gauge for repair;
setting the fluid heat transfer part as the density, constant pressure heat capacity and specific heat rate of fluid from material in geometric region, and setting the heat insulation part and the initial temperature value of the working medium gas;
in the laminar flow section, setting the fluid property to compressible flow; adding a gravity node to calculate the gravity borne by the arc calculation domain; adding a volume force node and Lorentz force to calculate the electromagnetic force of the geometric model; setting working medium gas inlet conditions, including setting inlet speed boundary conditions and initial normal inflow speed values; setting working medium gas outlet conditions, including setting outlet boundary pressure conditions, initial pressure, hydrostatic pressure compensation approximation, normal flow and backflow inhibition;
in the multi-physical field part, a balanced discharge heat source is arranged to act on the geometric model, the electric arc part is provided with a current electromagnetic coupling part and a fluid heat transfer part, the heat source components comprise enthalpy transfer, Joule heat and volume net radiation loss, and the total volume radiation system is arranged to be from a material; setting a static current density component and an induced current density component to act on the geometric model; setting Lorentz force to act on an arc part of the geometric model; setting the laminar flow part as flow coupling; setting the whole geometric area as temperature coupling; setting a balance discharge boundary heat source in the geometric model, and setting the surface work function of an anode, the surface work function of a cathode, the effective Richardson number, the effective work function and the plasma ionization potential;
arranging a current terminal and a ground at the current part, and setting at least two current values for the current terminal;
(5) carrying out mesh division on the geometric model, wherein the number of mesh units is 9380;
(6) configuring a time unit, a time step length, a time step and tolerance of a solver, and performing solution calculation to obtain a simulation result of each current value;
(7) and processing each simulation result to generate a temperature distribution map and an axial current density distribution map of the arc plasma corresponding to different current values.
Optionally, step (3) comprises:
in Comsol multiprophy,
selecting working medium material for the electrode part, and setting the relative magnetic permeability of the working medium material to be 1 and the electric conductivity to be 3 multiplied by 105S/m and a relative dielectric constant of 1; and
and setting working medium gas for the fluid part, and setting the density, constant-pressure heat capacity, heat conductivity coefficient, dynamic viscosity and electric conductivity of the working medium gas in an interpolation method mode.
Optionally, step (4) further comprises:
in Comsol multiprophy,
setting an initial value of a magnetic vector potential of the magnetic field part to 0 and a divergence condition variable scaling value to 1A/m;
setting argon gas in the working medium gas of the fluid heat transfer part, wherein the initial temperature value of argon gas ionization is 6000K;
setting the initial normal inflow speed value of a working medium gas inlet to be 20m/s and the initial pressure of a working medium gas outlet to be 0Pa in the laminar flow part;
the surface work function of the anode is set to be 4.15V, the surface work function of the cathode is set to be 4.15V, and the effective Richcson number is set to be 120A/(m)2·K2) The effective work function was set to 2.6V and the plasma ionization potential was set to 15.7V.
Optionally, step (5) comprises:
in Commol Multiphysics, performing free triangular mesh division on the geometric model in a self-defined form; the unit size of the electrode part is selected to be conventional, the electrode boundary is set as a boundary layer grid, the number of boundary layers is set to be 8, the boundary layer tension factor is set to be 1.2, the first layer thickness is selected automatically, and the thickness adjusting factor is set to be 1; the cell size of the arc portion is selected to be refined.
Optionally, step (6) comprises:
in the Commol Multiphysics, a time unit is selected to be ms, a time step is set to be 0.1ms, a time step is set to be 20ms, tolerance is set, and a solver is used for solving simulation data corresponding to each current value based on the time unit, the time step and the tolerance to obtain a simulation result under each current value.
Optionally, step (7) comprises:
in Comsol multiprophy,
selecting a two-dimensional drawing group from the temperature distribution diagram of the arc plasma, adopting a surface treatment mode, selecting temperature according to an expression, and representing and displaying temperature fields of the arc plasma with different current values;
selecting the simulation result of each current value of the last time step from the axial current density distribution diagram, and selecting a one-dimensional drawing group to represent and display the axial current density distribution of the arc plasma with different current values; wherein the horizontal axis is a distance of the vertical coordinate of the geometric model, and the vertical axis is a current density value.
In order to solve the above technical problem, according to still another aspect of an embodiment of the present invention, there is provided a simulation analysis apparatus for arc plasma characteristics, including:
the creating module is used for creating two-dimensional axisymmetric space dimensionality and selecting direct-current coupling discharge;
the modeling module is used for establishing a geometric model of arc plasma discharge; wherein the geometric model comprises an electrode portion and an arc portion, and the computational domain of the geometric model comprises a geometric region and an arc computational domain;
a first setting module for setting material properties of the electrode part and the arc part, respectively;
a second setting module for setting the magnetic field part, the fluid heat transfer part, the laminar flow part, the multi-physical field part, and the current part, respectively;
wherein the content of the first and second substances,
setting the magnetic field part to have relative magnetic permeability, electric conductivity and relative dielectric constant from materials, setting an initial value of magnetic vector potential and a divergence condition variable scaling value, and adding vector magnetic potential gauge for repair;
setting the fluid heat transfer part as the density, constant pressure heat capacity and specific heat rate of fluid from material in geometric region, and setting the heat insulation part and the initial temperature value of the working medium gas;
in the laminar flow section, setting the fluid property to compressible flow; adding a gravity node to calculate the gravity borne by the arc calculation domain; adding a volume force node and Lorentz force to calculate the electromagnetic force of the geometric model; setting working medium gas inlet conditions, including setting inlet speed boundary conditions and initial normal inflow speed values; setting working medium gas outlet conditions, including setting outlet boundary pressure conditions, initial pressure, hydrostatic pressure compensation approximation, normal flow and backflow inhibition;
in the multi-physical field part, a balanced discharge heat source is arranged to act on the geometric model, the electric arc part is provided with a current electromagnetic coupling part and a fluid heat transfer part, the heat source components comprise enthalpy transfer, Joule heat and volume net radiation loss, and the total volume radiation system is arranged to be from a material; setting a static current density component and an induced current density component to act on the geometric model; setting Lorentz force to act on an arc part of the geometric model; setting the laminar flow part as flow coupling; setting the whole geometric area as temperature coupling; setting a balance discharge boundary heat source in the geometric model, and setting the surface work function of an anode, the surface work function of a cathode, the effective Richardson number, the effective work function and the plasma ionization potential;
arranging a current terminal and a ground at the current part, and setting at least two current values for the current terminal;
the division module is used for carrying out grid division on the geometric model, and the number of grid units is 9380;
the solving module is used for configuring time units, time step lengths, time steps and tolerance of a solver, and carrying out solving calculation to obtain a simulation result of each current value;
and the generating module is used for processing each simulation result and generating a temperature distribution map and an axial current density distribution map of the arc plasma corresponding to different current values.
Optionally, the first setting module is further configured to:
in Comsol multiprophy,
selecting working medium material for the electrode part, and setting the relative magnetic permeability of the working medium material to be 1 and the electric conductivity to be 3 multiplied by 105S/m and a relative dielectric constant of 1; and
and setting working medium gas for the fluid part, and setting the density, constant-pressure heat capacity, heat conductivity coefficient, dynamic viscosity and electric conductivity of the working medium gas in an interpolation method mode.
Optionally, the second setting module is further configured to:
in Comsol multiprophy,
setting an initial value of a magnetic vector potential of the magnetic field part to 0 and a divergence condition variable scaling value to 1A/m;
setting argon gas in the working medium gas of the fluid heat transfer part, wherein the initial temperature value of argon gas ionization is 6000K;
setting the initial normal inflow speed value of a working medium gas inlet to be 20m/s and the initial pressure of a working medium gas outlet to be 0Pa in the laminar flow part;
the surface work function of the anode is set to be 4.15V, the surface work function of the cathode is set to be 4.15V, and the effective Richcson number is set to be 120A/(m)2·K2) The effective work function was set to 2.6V and the plasma ionization potential was set to 15.7V.
Optionally, the dividing module is further configured to:
in Commol Multiphysics, performing free triangular mesh division on the geometric model in a self-defined form; the unit size of the electrode part is selected to be conventional, the electrode boundary is set as a boundary layer grid, the number of boundary layers is set to be 8, the boundary layer tension factor is set to be 1.2, the first layer thickness is selected automatically, and the thickness adjusting factor is set to be 1; the cell size of the arc portion is selected to be refined.
Optionally, the solving module is further configured to:
in the Commol Multiphysics, a time unit is selected to be ms, a time step is set to be 0.1ms, a time step is set to be 20ms, tolerance is set, and a solver is used for solving simulation data corresponding to each current value based on the time unit, the time step and the tolerance to obtain a simulation result under the condition of each current value.
Optionally, the generating module is further configured to:
in Comsol multiprophy,
selecting a two-dimensional drawing group from the temperature distribution diagram of the arc plasma, adopting a surface treatment mode, selecting temperature according to an expression, and representing and displaying temperature fields of the arc plasma with different current values;
selecting the simulation result of each current value of the last time step from the axial current density distribution diagram, and selecting a one-dimensional drawing group to represent and display the axial current density distribution of the arc plasma with different current values; wherein the horizontal axis is a distance of the vertical coordinate of the geometric model, and the vertical axis is a current density value.
To achieve the above object, according to still another aspect of an embodiment of the present invention, there is provided an electronic device simulating characteristics of an arc plasma.
The electronic equipment for simulating the arc plasma characteristics of the embodiment of the invention comprises: one or more processors; a storage device for storing one or more programs which, when executed by the one or more processors, cause the one or more processors to implement a method for simulation analysis of arc plasma characteristics in accordance with an embodiment of the present invention.
To achieve the above object, according to still another aspect of embodiments of the present invention, there is provided a computer-readable storage medium.
A computer-readable storage medium of an embodiment of the present invention has stored thereon a computer program that, when executed by a processor, implements a method of simulation analysis of arc plasma characteristics of an embodiment of the present invention.
The simulation analysis method and the simulation analysis device for the arc plasma characteristics have the following beneficial effects: the method can simulate the characteristic distribution of the arc plasma under different current values based on Comsol Multiphysics, realizes the research on the influence of the current on the characteristics of the arc plasma, can efficiently and accurately obtain the characteristic distribution of the arc plasma, is convenient to operate, greatly reduces the test time and cost, and provides theoretical guidance and technical support for the application of the arc plasma.
Drawings
FIG. 1 is a schematic diagram of the main steps of a simulation analysis method of arc plasma characteristics according to an embodiment of the present invention;
FIG. 2 is a schematic flow chart diagram of a simulation analysis method of arc plasma characteristics according to a reference embodiment of the present invention;
FIG. 3 is a schematic view of a geometric model of a reference embodiment of the present invention;
FIG. 4 is a graph of the temperature profile of an arc plasma corresponding to 200A of one referenced embodiment of the present invention;
FIG. 5 is a temperature profile of an arc plasma in accordance with one referenced embodiment of the present invention at 150A;
FIG. 6 is a temperature profile of an arc plasma according to an embodiment of the present invention, generally designated 100A;
FIG. 7 is an axial current density profile of a reference embodiment of the present invention;
FIG. 8 is a schematic diagram of the main blocks of an apparatus for simulation analysis of arc plasma characteristics according to an embodiment of the present invention;
FIG. 9 is an exemplary system architecture diagram in which embodiments of the present invention may be employed;
fig. 10 is a schematic block diagram of a computer system suitable for use in implementing a terminal device or server according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
Because the arc plasma forming process involves complex multi-physical field coupling, real-time monitoring of the arc plasma characteristic parameters is limited, and the characteristics of the arc plasma at different current values are difficult to study.
Comsol multiprophy is a piece of numerical simulation software based on finite element methods. When the simulation research of the arc plasma is realized by using Commol Multiphysics, the characteristic distribution of the arc plasma under different current values can be accurately obtained, so that the influence of the current on the characteristics of the arc plasma is researched.
As shown in fig. 1, a simulation analysis method for arc plasma characteristics according to an embodiment of the present invention mainly includes the following steps:
step S101, establishing two-dimensional axisymmetric spatial dimension, and selecting direct current coupling discharge.
It is possible to establish two-dimensional axisymmetric spatial dimensions in Comsol Multiphysics and select and set suitable physical fields based on the properties of the arc plasma, which can simulate real physical conditions by means of current, magnetic, fluid heat transfer, laminar and multi-physical field sections. In addition, transient studies may be selected.
And S102, establishing a geometric model of arc plasma discharge.
The geometric model is used to simulate an arc plasma discharge. The geometric model includes an electrode part and an arc part, and the calculation domain of the geometric model includes a geometric area and an arc calculation domain, and at the same time, the geometric model is a complex formed by the electrode part and the arc part, such as shown in fig. 3.
And step S103, setting the material properties of the electrode part and the arc part respectively.
The computational domain of the geometric model comprises two parts, and thus its material properties are set for the two part computational domain separately.
In the embodiment of the present invention, step S103 may be implemented in the following manner: in Comsol Multiphysics, a working substance material is selected for an electrode part, and the relative permeability, the conductivity and the relative dielectric constant of the working substance material are set; and setting working medium gas for the fluid part, and setting the density, constant-pressure heat capacity, heat conductivity coefficient, dynamic viscosity and electric conductivity of the working medium gas in an interpolation method mode.
Interpolation, also known as "interpolation", is a method in which a function f (x) is used to insert the function values of a number of points into a certain interval, a suitable specific function is created, the known values are taken at these points, and the values of the specific function are used as the approximate values of the function f (x) at other points in the interval, and this method is called interpolation. As a preferred embodiment, the electrode part can be made of graphite material, the relative magnetic permeability is set to be 1, and the electric conductivity is set to be 3 multiplied by 105S/m, the relative dielectric constant is set to 1; the working gas can be set to argon, whose density, constant pressure heat capacity, thermal conductivity, dynamic viscosity and electrical conductivity are temperature dependent and therefore are set by interpolation.
Step S104, respectively setting a magnetic field part, a fluid heat transfer part, a laminar flow part, a multi-physical field part and a current part.
In order to study the influence of the current magnitude on the arc plasma characteristics, at least two current values can be set, so that the arc distribution characteristics under different current values can be studied.
In particular, the amount of the solvent to be used,
setting the magnetic field part as relative magnetic conductivity, electric conductivity and relative dielectric constant from materials, setting the initial value of magnetic vector potential and divergence condition variable scaling value, and adding vector magnetic potential gauge for repair;
setting the fluid heat transfer part as the density, constant pressure heat capacity and specific heat rate of fluid from material in geometric area, and setting the initial temperature values of the heat insulation part and the working medium gas;
in the laminar flow section, the fluid properties are set to compressible flow; adding a gravity node to calculate the gravity borne by the arc calculation domain; adding a volume force node and Lorentz force to calculate the electromagnetic force of the geometric model; setting working medium gas inlet conditions, including setting inlet speed boundary conditions and initial normal inflow speed values; setting working medium gas outlet conditions, including setting outlet boundary pressure conditions, initial pressure, hydrostatic pressure compensation approximation, normal flow and backflow inhibition;
in the multi-physical field part, a balanced discharge heat source is arranged to act on a geometric model, the electric arc part is provided with a current electromagnetic coupling part and a fluid heat transfer part, the heat source component comprises enthalpy transfer, Joule heat and volume net radiation loss, and the total volume radiation system is arranged to come from a material; setting a static current density component and an induced current density component to act on the geometric model; setting Lorentz force to act on an arc part of the geometric model; setting the laminar flow part as flow coupling; setting the whole geometric area as temperature coupling; setting a balance discharge boundary heat source in a geometric model, and setting the surface work function of an anode, the surface work function of a cathode, an effective Richardson constant (Richardson constant), an effective work function and plasma ionization potential;
a current terminal and a ground are provided in the current section, and at least two current values are set for the current terminal
Wherein the control equations of the various parts are as follows:
current part:
magnetic field part:
a fluid heat transfer portion:
laminar flow part:
lorentz force in the multiphysics part:
in the formula (I), the compound is shown in the specification,
is the current density vector, unit A/m2；
Is a volume current source, unit A/m3；
Is the induced electromotive force vector, unit V/m; v is the potential, in units of V;
is the vector of the electric flux density, unit C/m2；
Is the density of the externally generated current in units of A/m2(ii) a t is time;
is the magnetic vector intensity, in units of A/m;
is the magnetic induction vector, unit T;
is the magnetic vector potential, in Wb/m; σ is the conductivity, in S/m;
is density, unit kg/m3；
Is the constant pressure heat capacity, unit J/(kg. K);
is the velocity vector, in m/s;
is the heat flux in W/m2(ii) a K is the thermal conductivity, in W/(m.K); t is temperature in K, Q,
、
Is a heat source term, unit W/m3；
Is the gravitational acceleration in m/s;
is pressure, in Pa;
is the vector of the volume force (Lorentz force), in N/m3；
Is the stress term.
In this embodiment of the present invention, step S104 may further include the following specific setting contents: in Comsol Multiphysics, the initial value of the magnetic vector potential of the magnetic field portion is set to 0, and the divergence condition variable scaling value is set to 1A/m; setting argon gas as working medium gas of the fluid heat transfer part, wherein the initial temperature value of argon gas ionization is 6000K; in the laminar flow part, setting the initial normal inflow speed value of a working medium gas inlet to be 20m/s and the initial pressure of a working medium gas outlet to be 0 Pa; the surface work function of the anode is set to be 4.15V, the surface work function of the cathode is set to be 4.15V, and the effective Richcson number is set to be 120A/(m)2·K2) The effective work function was set to 2.6V and the plasma ionization potential was set to 15.7V.
Also in fig. 3, FG is set as an anode and BCDE is set as a cathode. Vector magnetomotive gauge repairs are added to calculate the magnetic field and improve the convergence of the calculations. The initial temperature value here is an initial value that is set, and the initial temperature value changes as calculated. The volume force is the non-contact force acting on all fluid elements across the space.
And step S105, carrying out mesh division on the geometric model.
The grid division can influence the calculation precision and the convergence process of the subsequent solution, and the number of grid units divided in the step is 9380.
In the embodiment of the present invention, step S105 may be implemented in the following manner: in Comsol multiprophy, the geometric model is free triangulated in a custom form.
In mesh partitioning in Comsol multiprophy, the following parameters may be set: the unit size of the electrode part is selected to be conventional, the electrode boundary is set as boundary layer grids, the number of boundary layers is set to be 8, the boundary layer tension factor is set to be 1.2, the thickness of the first layer is selected automatically, and the thickness adjusting factor is set to be 1; the cell size of the arc portion is selected to be fine. Wherein the electrode boundary is a boundary of the electrode portion.
And S106, configuring time units, time step lengths, time steps and tolerance of a solver, and performing solving calculation to obtain a simulation result of each current value.
The solver provided by Comsol Multiphysics can be configured based on requirements, such as time units, time step sizes, time steps, and tolerances. The result of solving and calculating the multiple physical fields includes a temperature field, a speed field, an electric field and a magnetic field, that is, the simulation result may include the contents of the temperature field, the speed field, the electric field and the magnetic field. The simulation environment is set through the steps S101 to S105, and in the same simulation environment, along with the progress of the experiment, the arc distribution characteristics generated by different current values are also different, that is, different current values can obtain different simulation data, the simulation data at each current value is separately solved and calculated, and the obtained calculation result is the corresponding simulation result. It should be noted that the control equation based solver is capable of performing calculations on its own logic and automatically on parameters in the control equation by default if not actively set.
In the embodiment of the present invention, step S106 may be implemented in the following manner: in the Commol Multiphysics, a time unit is selected to be ms, a time step is set to be 0.1ms, a time step is set to be 20ms, tolerance is set, and a solver solves simulation data corresponding to each current value based on the time unit, the time step and the tolerance to obtain a simulation result under the condition of each current value.
Wherein the time step comprises a step size, a calculation start time and a calculation end time. In Comsol multiprophy, tolerances can be selected for user control and relative tolerances actively set and solved.
And S107, processing and analyzing each simulation result to generate a temperature distribution map and an axial current density distribution map of the arc plasma corresponding to different current values.
The Commol Multiphysics can process and analyze the simulation result and display the simulation result in a form of a graph. Wherein the temperature profile of the arc plasma (e.g., fig. 4, 5, 6) is used to show the effect of each current value on the arc temperature profile, the axial current density profile is used to show the axial current density profiles at different current values, and the axial current density profiles at different current values can be plotted on the same graph (e.g., fig. 7).
In the embodiment of the present invention, step S107 may be implemented in the following manner: in Commol Multiphysics, a two-dimensional drawing group is selected for the temperature distribution diagram of the arc plasma, a surface treatment mode is adopted, the temperature is selected by an expression, and the temperature fields of the arc plasma with different current values are represented and displayed; and selecting simulation results of all current values of the last time step from the axial current density distribution maps at different current values, and selecting a one-dimensional drawing group to represent and display the axial current density distribution of the arc plasma at different current values.
In the axial current density distribution diagram at different current values, the horizontal axis may be a distance from the vertical coordinate of the geometric model (i.e., a distance from the origin), and the vertical axis may be a current density value.
The present invention will be further described by way of examples in order to facilitate a more complete, accurate and thorough understanding of the concepts and solutions of the present invention and to facilitate its implementation by those skilled in the art, but the scope of the present invention is not limited to these examples.
Example one
As shown in fig. 2, in the embodiment of the present invention, a simulation analysis method for arc plasma characteristics is implemented based on Comsol Multiphysics, and the method may be implemented with reference to the following process:
1. selecting spatial dimensions and physical fields
Establishing two-dimensional axisymmetric spatial dimensions in Commol Multiphysics, selecting direct current coupling discharge based on the properties of arc plasma, setting a current part, a magnetic field part, a fluid heat transfer part, a laminar flow part and a multi-physical field part, and selecting transient state research.
2. Establishing arc plasma geometric model
The computational domain of the geometric model includes an electrode portion and an arc portion, and in the geometric model shown in fig. 3, ABCDE is the electrode portion and EHGFBCD is the argon portion. In Comsol Multiphysics, BF was set to 10mm, FG was set to 9.4mm, HG was set to 14mm, BC was set to 0.3mm, DE was set to 3mm, AE was set to 0.75mm, AB was set to 4mm, and an electrode part and an arc part were set to form a united body.
3. Setting material properties according to different calculation domains
Since the computation domain comprises two parts, its material properties are set for the two part computation domain separately. In Comsol Multiphysics, the electrode part is made of graphite material, the relative permeability is set to 1, and the electrical conductivity is set to 3 × 105S/m, the relative dielectric constant is set to 1; the fluid part and the working medium gas are argon, and the density, constant-pressure heat capacity, heat conductivity coefficient, dynamic viscosity and electric conductivity of the argon are related to the temperature, so that the argon is arranged by an interpolation method.
4. Each part is arranged
That is, in Comsol Multiphysics, a magnetic field section, a fluid heat transfer section, a laminar flow section, a multi-physical field section, and a current section are separately provided.
Wherein the content of the first and second substances,
for the current part, AE is set as the current terminal, current value I0Setting FG as a variable, setting FG as a ground, and setting at least two current values for a current terminal;
setting the relative magnetic conductivity, the electric conductivity and the relative dielectric constant of the magnetic field part as materials, setting an initial value of magnetic vector potential and a divergence condition variable scaling value, and adding a vector magnetic potential gauge for repairing;
in the fluid heat transfer part, the density, constant-pressure heat capacity and specific heat rate of fluid in a geometric area are set to be from materials, BCDE and FG are set to be thermally insulated, and the thermal insulation part and the initial temperature value of working medium gas are set, namely one initial temperature value giving argon ionization conditions is 6000K;
in the laminar flow part, the fluid property is set to be compressible flow; setting EH as a working medium gas inlet, and setting inlet speed boundary conditions and an initial normal inflow speed value to be 20 m/s; setting HG as an argon outlet, setting outlet boundary pressure conditions, initial pressure, hydrostatic pressure compensation approximation, normal flow and backflow inhibition; adding a volume force node and Lorentz force to calculate the electromagnetic force of the geometric model; adding a gravity node to calculate the gravity borne by the arc calculation domain, selecting a gravity option, and inputting the gravity acceleration in the z direction "-g _ const m/s2", to add a force of gravity to the z direction, the magnitude and direction of the force of gravity must be constant after a given.
In the multi-physical field part, electrode properties are set, specifically, FG is set as an anode, and BCDE is set as a cathode. Setting a balanced discharge heat source to act on the geometric model, setting current as an electromagnetic coupling part in the arc part, setting fluid heat transfer as a heat transfer part, and setting heat source components including enthalpy transfer, joule heat and volume net radiation loss, wherein the total volume radiation system is set to come from a material; setting a static current density component and an induced current density component to act on the geometric model; setting Lorentz force to act on an arc part of the geometric model; setting the laminar flow part as flow coupling; setting the whole geometric area as temperature coupling; and setting a balance discharge boundary heat source in the geometric model, and setting the surface work function of the anode, the surface work function of the cathode, the effective Richardson constant, the effective work function and the plasma ionization potential.
The control equations for each section are as follows:
current part:
magnetic field part:
a fluid heat transfer portion:
laminar flow part:
lorentz force in the multiphysics part:
in the formula (I), the compound is shown in the specification,
is the current density vector, unit A/m2；
Is a volume current source, unit A/m3；
Is the induced electromotive force vector, unit V/m; v is the potential, in units of V;
is the vector of the electric flux density, unit C/m2；
Is the density of the externally generated current in units of A/m2(ii) a t is time;
is the magnetic vector intensity, in units of A/m;
is the magnetic induction vector, unit T;
is the magnetic vector potential, in Wb/m; σ is the conductivity, in S/m;
is density, unit kg/m3；
Is the constant pressure heat capacity, unit J/(kg. K);
is the velocity vector, in m/s;
is the heat flux in W/m2(ii) a K is the thermal conductivity, in W/(m.K); t is temperature in K, Q,
、
Is a heat source term, unit W/m3；
Is the gravitational acceleration in m/s;
is pressure, in Pa;
is the vector of the volume force (Lorentz force), in N/m3；
Is the stress term.
5. Meshing arc plasma models
In Comsol Multiphysics, an arc plasma model (i.e., a geometric model) is freely triangulated, and the number of grid cells after the freely triangulated division is 9380. The grid division mode is customized by a user, and before grid division, the calculation area is divided into two parts, wherein the unit size of the electrode part is selected as conventional, and the unit size of the arc part is selected as refined; setting electrode boundaries as boundary layer meshes; the number of boundary layers is 8, the boundary layer tension factor is 1.2, the thickness of the first layer is automatically selected, and the thickness adjusting factor is set to be 1.
6. Configure solver and compute
In Comsol multiprophy, the time unit is set to ms, the time step is set to 0.1ms, and the time step is set to 20 ms. The tolerance selection user control, with the relative tolerance set to 0.01, is solved.
It should be noted that after one computation is completed, the settings can be checked and modified if not converged.
7. Result processing and analysis
Selecting a two-dimensional drawing group for the temperature distribution diagram of the arc plasma, selecting 20ms of time, selecting a surface treatment mode, selecting temperature by an expression, and representing and displaying the temperature field of the arc plasma under different current values;
and selecting a one-dimensional drawing group for the axial current density distribution diagram, analyzing the axial current density distribution of the arc plasma with different current values, selecting the time as the last one, selecting BF for the data of the horizontal axis, and selecting a current density mode for the data of the vertical axis.
Example two
In order to study the influence of the current magnitude on the arc plasma characteristics, at least two current values can be set, so that the arc distribution characteristics under different current values can be studied.
In the previous example, except for the current values set to 200A, 150A and 100A, respectively, the rest settings were kept unchanged, and simulation analysis was performed according to steps S201 to S207, and finally temperature profiles of the arc plasma corresponding to different current values and axial current density profiles at different current values were generated, respectively.
Fig. 4, 5, and 6 are temperature profiles of arc plasmas for 200A, 150A, and 100A, respectively. It should be noted that, when the temperature distribution diagram of the arc plasma corresponding to the current value (i.e. fig. 4, 5 and 6) is shown, different colors may be used to represent different temperature levels, for example, the colors such as purple, red, orange, yellow, light green, light blue, blue and dark blue and their gradation colors represent the highest temperature to the lowest temperature, and fig. 4, 5 and 6 show that the colors mainly displayed by diffusing outward with the lower end of the electrode as the center are purple, red, orange, yellow, light green, light blue, blue and dark blue in sequence, wherein the colors except for the dark blue portion constitute "bell jar" and "bell jar".
As can be seen from fig. 4, 5 and 6, the arc temperature distribution as a whole exhibits a "bell-jar" shape, which is consistent with the measurements under the prior art conditions, with the arc having the highest axial temperature and gradually decreasing along the axial temperature. The temperature near the anode is reduced seriously, and the phenomenon of joule heat is obviously reduced. As the current value is continuously decreased, the maximum temperature of the arc plasma is also continuously decreased. The current provides energy for the formation of arc plasma, and the influence of the current on the arc is more remarkable.
Based on fig. 4, 5 and 6, the maximum temperature of the arc was 26800K, 24400K, 20200K, respectively, at currents 200A, 150A, 100A. The maximum temperature of the arc is reduced along with the reduction of the current value, but under different current values, the axial temperature distribution of the arc tends to be consistent, the maximum value appears near the cathode, the temperature gradually decreases in the arc column area, and the steep drop trend appears near the anode. The temperature near the anode is reduced seriously, and the phenomenon of joule heat is obviously reduced. The current provides energy for the formation of arc plasma, and the influence of the current on the arc is more remarkable. (it can be seen from FIG. 7 that the maximum value of the current density at the cathode gradually decreases as the current value decreases.)
Fig. 7 shows the axial current density distribution in the case of different current values, and it can be seen from fig. 7 that the maximum value of the current density at the cathode gradually decreases as the current value decreases. The horizontal axis of fig. 7 represents points on the axis of symmetry BF, and the vertical axis represents current density values corresponding to these points.
According to the two embodiments, the simulation of the arc plasma is realized by utilizing the advantage of multi-physical field coupling of Commol Multiphysics, the characteristic change of the arc plasma under different currents can be conveniently and quickly researched, and meanwhile, the characteristics of the arc plasma under different current values are calculated.
In addition, as shown in fig. 8, an embodiment of the present invention further provides a simulation analysis apparatus for arc plasma characteristics, which includes a creating module 801, a modeling module 802, a first setting module 803, a second setting module 804, a dividing module 805, a solving module 806, and a generating module 807.
Wherein the content of the first and second substances,
a creating module 801, configured to create a two-dimensional axisymmetric spatial dimension, and select dc coupling discharge;
a modeling module 802 for establishing a geometric model of arc plasma discharge; wherein the geometric model comprises an electrode portion and an arc portion, and the computational domain of the geometric model comprises a geometric region and an arc computational domain;
a first setting module 803 for setting material properties of the electrode part and the arc part, respectively;
a second setting module 804 for setting the magnetic field portion, the fluid heat transfer portion, the laminar flow portion, the multi-physical field portion, and the current portion, respectively;
wherein the content of the first and second substances,
setting the magnetic field part to have relative magnetic permeability, electric conductivity and relative dielectric constant from materials, setting an initial value of magnetic vector potential and a divergence condition variable scaling value, and adding vector magnetic potential gauge for repair;
setting the fluid heat transfer part as the density, constant pressure heat capacity and specific heat rate of fluid from material in geometric region, and setting the heat insulation part and the initial temperature value of the working medium gas;
in the laminar flow section, setting the fluid property to compressible flow; adding a gravity node to calculate the gravity borne by the arc calculation domain; adding a volume force node and Lorentz force to calculate the electromagnetic force of the geometric model; setting working medium gas inlet conditions, including setting inlet speed boundary conditions and initial normal inflow speed values; setting working medium gas outlet conditions, including setting outlet boundary pressure conditions, initial pressure, hydrostatic pressure compensation approximation, normal flow and backflow inhibition;
in the multi-physical field part, a balanced discharge heat source is arranged to act on the geometric model, the electric arc part is provided with a current electromagnetic coupling part and a fluid heat transfer part, the heat source components comprise enthalpy transfer, Joule heat and volume net radiation loss, and the total volume radiation system is arranged to be from a material; setting a static current density component and an induced current density component to act on the geometric model; setting Lorentz force to act on an arc part of the geometric model; setting the laminar flow part as flow coupling; setting the whole geometric area as temperature coupling; setting a balance discharge boundary heat source in the geometric model, and setting the surface work function of an anode, the surface work function of a cathode, the effective Richardson number, the effective work function and the plasma ionization potential;
arranging a current terminal and a ground at the current part, and setting at least two current values for the current terminal;
a division module 805, configured to perform mesh division on the geometric model, where the number of mesh units is 9380;
a solving module 806, configured to configure a time unit, a time step, and a tolerance of the solver, and perform solving calculation to obtain a simulation result of each current value;
a generating module 807 for processing each simulation result to generate a temperature distribution map and an axial current density distribution map of the arc plasma corresponding to different current values.
In this embodiment of the present invention, the first setting module 803 may further be configured to:
in Comsol multiprophy,
selecting working medium material for the electrode part, and setting the relative magnetic permeability of the working medium material to be 1 and the electric conductivity to be 3 multiplied by 105S/m and a relative dielectric constant of 1; and
and setting working medium gas for the fluid part, and setting the density, constant-pressure heat capacity, heat conductivity coefficient, dynamic viscosity and electric conductivity of the working medium gas in an interpolation method mode.
In this embodiment of the present invention, the second setting module 804 may be further configured to:
in Comsol multiprophy,
setting an initial value of a magnetic vector potential of the magnetic field part to 0 and a divergence condition variable scaling value to 1A/m;
setting argon gas in the working medium gas of the fluid heat transfer part, wherein the initial temperature value of argon gas ionization is 6000K;
setting the initial normal inflow speed value of a working medium gas inlet to be 20m/s and the initial pressure of a working medium gas outlet to be 0Pa in the laminar flow part;
the surface work function of the anode is set to be 4.15V, the surface work function of the cathode is set to be 4.15V, and the effective Richcson number is set to be 120A/(m)2·K2) The effective work function was set to 2.6V and the plasma ionization potential was set to 15.7V.
In this embodiment of the present invention, the dividing module 805 may further be configured to:
in Commol Multiphysics, performing free triangular mesh division on the geometric model in a self-defined form; the unit size of the electrode part is selected to be conventional, the electrode boundary is set as a boundary layer grid, the number of boundary layers is set to be 8, the boundary layer tension factor is set to be 1.2, the first layer thickness is selected automatically, and the thickness adjusting factor is set to be 1; the cell size of the arc portion is selected to be refined.
In this embodiment of the present invention, the solving module 806 may further be configured to:
in the Commol Multiphysics, a time unit is selected to be ms, a time step is set to be 0.1ms, a time step is set to be 20ms, tolerance is set, and a solver is used for solving simulation data corresponding to each current value based on the time unit, the time step and the tolerance to obtain a simulation result under the condition of each current value.
In this embodiment of the present invention, the generating module 807 may further be configured to:
in Comsol multiprophy,
selecting a two-dimensional drawing group from the temperature distribution diagram of the arc plasma, adopting a surface treatment mode, selecting temperature according to an expression, and representing and displaying temperature fields of the arc plasma with different current values;
selecting the simulation result of each current value of the last time step from the axial current density distribution diagram, and selecting a one-dimensional drawing group to represent and display the axial current density distribution of the arc plasma with different current values; wherein the horizontal axis is a distance of the vertical coordinate of the geometric model, and the vertical axis is a current density value.
Fig. 9 illustrates an exemplary system architecture 900 of a method for simulation analysis of arc plasma characteristics or a device for simulation analysis of arc plasma characteristics to which embodiments of the present invention may be applied.
As shown in fig. 9, the system architecture 900 may include end devices 901, 902, 903, a network 904, and a server 905. Network 904 is the medium used to provide communication links between terminal devices 901, 902, 903 and server 905. Network 904 may include various connection types, such as wired, wireless communication links, or fiber optic cables, to name a few.
A user may use the terminal devices 901, 902, 903 to interact with a server 905 over a network 904 to receive or send messages and the like. Various communication client applications can be installed on the terminal devices 901, 902, 903.
The terminal devices 901, 902, 903 may be various electronic devices having a display screen and supporting web browsing, including but not limited to smart phones, tablet computers, laptop portable computers, desktop computers, and the like.
The server 905 may be a server that provides various services, such as a background management server for providing support. The background management server can analyze and process the received data such as the product information inquiry request and feed back the processing result to the terminal equipment.
It should be noted that the simulation analysis method for arc plasma characteristics provided by the embodiment of the present invention is generally executed by the server 905, and accordingly, the simulation analysis device for arc plasma characteristics is generally disposed in the server 905.
It should be understood that the number of terminal devices, networks, and servers in fig. 9 is merely illustrative. There may be any number of terminal devices, networks, and servers, as desired for implementation.
Referring now to FIG. 10, a block diagram of a computer system 1000 suitable for use with a terminal device implementing an embodiment of the invention is shown. The terminal device shown in fig. 10 is only an example, and should not bring any limitation to the functions and the scope of use of the embodiments of the present invention.
As shown in fig. 10, the computer system 1000 includes a Central Processing Unit (CPU) 1001 that can perform various appropriate actions and processes according to a program stored in a Read Only Memory (ROM) 1002 or a program loaded from a storage section 1008 into a Random Access Memory (RAM) 1003. In the RAM 1003, various programs and data necessary for the operation of the system 1000 are also stored. The CPU 1001, ROM 1002, and RAM 1003 are connected to each other via a bus 1004. An input/output (I/O) interface 1005 is also connected to bus 1004.
The following components are connected to the I/O interface 1005: an input section 1006 including a keyboard, a mouse, and the like; an output section 1007 including a display such as a Cathode Ray Tube (CRT), a Liquid Crystal Display (LCD), and the like, and a speaker; a storage portion 1008 including a hard disk and the like; and a communication section 1009 including a network interface card such as a LAN card, a modem, or the like. The communication section 1009 performs communication processing via a network such as the internet. The driver 1010 is also connected to the I/O interface 1005 as necessary. A removable medium 1011 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, or the like is mounted on the drive 1010 as necessary, so that a computer program read out therefrom is mounted into the storage section 1008 as necessary.
In particular, according to the embodiments of the present disclosure, the processes described above with reference to the flowcharts may be implemented as computer software programs. For example, embodiments of the present disclosure include a computer program product comprising a computer program embodied on a computer readable medium, the computer program comprising program code for performing the method illustrated in the flow chart. In such an embodiment, the computer program may be downloaded and installed from a network through the communication part 1009 and/or installed from the removable medium 1011. The computer program executes the above-described functions defined in the system of the present invention when executed by the Central Processing Unit (CPU) 1001.
It should be noted that the computer readable medium shown in the present invention can be a computer readable signal medium or a computer readable storage medium or any combination of the two. A computer readable storage medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the foregoing. More specific examples of the computer readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the present invention, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. In the present invention, however, a computer readable signal medium may include a propagated data signal with computer readable program code embodied therein, for example, in baseband or as part of a carrier wave. Such a propagated data signal may take many forms, including, but not limited to, electro-magnetic, optical, or any suitable combination thereof. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to: wireless, wire, fiber optic cable, RF, etc., or any suitable combination of the foregoing.
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams or flowchart illustration, and combinations of blocks in the block diagrams or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
The modules described in the embodiments of the present invention may be implemented by software or hardware. The described modules may also be provided in a processor, which may be described as: a processor includes a creation module, a modeling module, a first setting module, a second setting module, a partitioning module, a solving module, and a generating module. The names of these modules do not in some cases form a limitation on the module itself, for example, a partitioning module may also be described as a "module that meshes the geometric model".
As another aspect, the present invention also provides a computer-readable medium that may be contained in the apparatus described in the above embodiments; or may be separate and not incorporated into the device. The computer readable medium carries one or more programs which, when executed by a device, cause the device to comprise: step S101-step S107.
In summary, the simulation analysis method and the simulation analysis apparatus for arc plasma according to the embodiments of the present invention have at least the following advantages:
the method can simulate the characteristic distribution of the arc plasma under different current values based on Comsol Multiphysics, realizes the research on the influence of the current on the characteristics of the arc plasma, can efficiently and accurately obtain the characteristic distribution of the arc plasma, is convenient to operate, greatly reduces the test time and cost, and provides theoretical guidance and technical support for the application of the arc plasma.
The invention has not been described in detail and is in part known to those of skill in the art.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A simulation analysis method for arc plasma characteristics is characterized by comprising the following steps:
(1) establishing two-dimensional axisymmetric spatial dimensions, and selecting direct-current coupling discharge;
(2) establishing a geometric model of arc plasma discharge; wherein the geometric model comprises an electrode portion and an arc portion, and the computational domain of the geometric model comprises a geometric region and an arc computational domain;
(3) setting material properties of the electrode part and the arc part, respectively;
(4) respectively arranging a magnetic field part, a fluid heat transfer part, a laminar flow part, a multi-physical field part and a current part;
wherein the content of the first and second substances,
setting the magnetic field part to have relative magnetic permeability, electric conductivity and relative dielectric constant from materials, setting an initial value of magnetic vector potential and a divergence condition variable scaling value, and adding vector magnetic potential gauge for repair;
setting the fluid heat transfer part as the density, constant pressure heat capacity and specific heat rate of fluid from material in geometric region, and setting the initial temperature values of the thermal insulation part and the working medium gas;
in the laminar flow section, setting the fluid property to compressible flow; adding a gravity node to calculate the gravity borne by the arc calculation domain; adding a volume force node and Lorentz force to calculate the electromagnetic force of the geometric model; setting working medium gas inlet conditions, including setting inlet speed boundary conditions and initial normal inflow speed values; setting working medium gas outlet conditions, including setting outlet boundary pressure conditions, initial pressure, hydrostatic pressure compensation approximation, normal flow and backflow inhibition;
in the multi-physical field part, a balanced discharge heat source is arranged to act on the geometric model, the electric arc part is provided with a current electromagnetic coupling part and a fluid heat transfer part, the heat source components comprise enthalpy transfer, Joule heat and volume net radiation loss, and the total volume radiation system is arranged to be from a material; setting a static current density component and an induced current density component to act on the geometric model; setting Lorentz force to act on an arc part of the geometric model; setting the laminar flow part as flow coupling; setting the whole geometric area as temperature coupling; setting a balance discharge boundary heat source in the geometric model, and setting the surface work function of an anode, the surface work function of a cathode, the effective Richardson number, the effective work function and the plasma ionization potential;
arranging a current terminal and a ground at the current part, and setting at least two current values for the current terminal;
(5) carrying out mesh division on the geometric model, wherein the number of mesh units is 9380;
(6) configuring a time unit, a time step length, a time step and tolerance of a solver, and performing solution calculation to obtain a simulation result of each current value;
(7) and processing each simulation result to generate a temperature distribution map and an axial current density distribution map of the arc plasma corresponding to different current values.
2. The method of claim 1, wherein step (3) comprises:
in Comsol multiprophy,
selecting working medium material for the electrode part, and setting the relative magnetic permeability of the working medium material to be 1 and the electric conductivity to be 3 multiplied by 105S/m and a relative dielectric constant of 1; and
and setting working medium gas for the electric arc part, and setting the density, constant-pressure heat capacity, heat conductivity coefficient, dynamic viscosity and electric conductivity of the working medium gas in an interpolation method mode.
3. The method of claim 1, wherein step (4) further comprises:
in Comsol multiprophy,
setting an initial value of a magnetic vector potential of the magnetic field part to 0 and a divergence condition variable scaling value to 1A/m;
setting argon gas in the working medium gas of the fluid heat transfer part, wherein the initial temperature value of argon gas ionization is 6000K;
setting the initial normal inflow speed value of a working medium gas inlet to be 20m/s and the initial pressure of a working medium gas outlet to be 0Pa in the laminar flow part;
the surface work function of the anode is set to be 4.15V, the surface work function of the cathode is set to be 4.15V, and the effective Richcson number is set to be 120A/(m)2·K2) The effective work function was set to 2.6V and the plasma ionization potential was set to 15.7V.
4. The method of claim 1, wherein step (5) comprises:
in Commol Multiphysics, performing free triangular mesh division on the geometric model in a self-defined form; the unit size of the electrode part is selected to be conventional, the electrode boundary is set as a boundary layer grid, the number of boundary layers is set to be 8, the boundary layer tension factor is set to be 1.2, the first layer thickness is selected automatically, and the thickness adjusting factor is set to be 1; the cell size of the arc portion is selected to be refined.
5. The method of claim 1, wherein step (6) comprises:
in the Commol Multiphysics, a time unit is selected to be ms, a time step is set to be 0.1ms, a time step is set to be 20ms, tolerance is set, and a solver is used for solving simulation data corresponding to each current value based on the time unit, the time step and the tolerance to obtain a simulation result under each current value.
6. The method of claim 1, wherein step (7) comprises:
in Comsol multiprophy,
selecting a two-dimensional drawing group from the temperature distribution diagram of the arc plasma, adopting a surface treatment mode, selecting temperature according to an expression, and representing and displaying temperature fields of the arc plasma with different current values;
selecting the simulation result of each current value of the last time step from the axial current density distribution diagram, and selecting a one-dimensional drawing group to represent and display the axial current density distribution of the arc plasma with different current values; wherein the horizontal axis is a distance of the vertical coordinate of the geometric model, and the vertical axis is a current density value.
7. An apparatus for simulation analysis of arc plasma characteristics, comprising:
the creating module is used for creating two-dimensional axisymmetric space dimensionality and selecting direct-current coupling discharge;
the modeling module is used for establishing a geometric model of arc plasma discharge; wherein the geometric model comprises an electrode portion and an arc portion, and the computational domain of the geometric model comprises a geometric region and an arc computational domain;
a first setting module for setting material properties of the electrode part and the arc part, respectively;
the second setting module is used for respectively setting a magnetic field part, a fluid heat transfer part, a laminar flow part, a multi-physical field part and a current part;
wherein the content of the first and second substances,
setting the magnetic field part to have relative magnetic permeability, electric conductivity and relative dielectric constant from materials, setting an initial value of magnetic vector potential and a divergence condition variable scaling value, and adding vector magnetic potential gauge for repair;
setting the fluid heat transfer part as the density, constant pressure heat capacity and specific heat rate of fluid from material in geometric region, and setting the initial temperature values of the thermal insulation part and the working medium gas;
in the laminar flow section, setting the fluid property to compressible flow; adding a gravity node to calculate the gravity borne by the arc calculation domain; adding a volume force node and Lorentz force to calculate the electromagnetic force of the geometric model; setting working medium gas inlet conditions, including setting inlet speed boundary conditions and initial normal inflow speed values; setting working medium gas outlet conditions, including setting outlet boundary pressure conditions, initial pressure, hydrostatic pressure compensation approximation, normal flow and backflow inhibition;
in the multi-physical field part, a balanced discharge heat source is arranged to act on the geometric model, the electric arc part is provided with a current electromagnetic coupling part and a fluid heat transfer part, the heat source components comprise enthalpy transfer, Joule heat and volume net radiation loss, and the total volume radiation system is arranged to be from a material; setting a static current density component and an induced current density component to act on the geometric model; setting Lorentz force to act on an arc part of the geometric model; setting the laminar flow part as flow coupling; setting the whole geometric area as temperature coupling; setting a balance discharge boundary heat source in the geometric model, and setting the surface work function of an anode, the surface work function of a cathode, the effective Richardson number, the effective work function and the plasma ionization potential;
arranging a current terminal and a ground at the current part, and setting at least two current values for the current terminal;
the division module is used for carrying out grid division on the geometric model, and the number of grid units is 9380;
the solving module is used for configuring time units, time step lengths, time steps and tolerance of a solver, and carrying out solving calculation to obtain a simulation result of each current value;
and the generating module is used for processing each simulation result and generating a temperature distribution map and an axial current density distribution map of the arc plasma corresponding to different current values.
8. The apparatus of claim 7, wherein the generating module is further configured to:
in Comsol multiprophy,
selecting a two-dimensional drawing group from the temperature distribution diagram of the arc plasma, adopting a surface treatment mode, selecting temperature according to an expression, and representing and displaying temperature fields of the arc plasma with different current values;
selecting the simulation result of each current value of the last time step from the axial current density distribution diagram, and selecting a one-dimensional drawing group to represent and display the axial current density distribution of the arc plasma with different current values; wherein the horizontal axis is a distance of the vertical coordinate of the geometric model, and the vertical axis is a current density value.
9. An electronic device for simulating an arc plasma, comprising:
one or more processors;
a storage device for storing one or more programs,
when executed by the one or more processors, cause the one or more processors to implement the method of any one of claims 1-7.
10. A computer-readable medium, on which a computer program is stored, which, when being executed by a processor, carries out the method according to any one of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110421377.4A CN112989674B (en) | 2021-04-20 | 2021-04-20 | Simulation analysis method and simulation analysis device for arc plasma characteristics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110421377.4A CN112989674B (en) | 2021-04-20 | 2021-04-20 | Simulation analysis method and simulation analysis device for arc plasma characteristics |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112989674A CN112989674A (en) | 2021-06-18 |
| CN112989674B true CN112989674B (en) | 2021-08-20 |
Family
ID=76341155
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110421377.4A Active CN112989674B (en) | 2021-04-20 | 2021-04-20 | Simulation analysis method and simulation analysis device for arc plasma characteristics |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112989674B (en) |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| ATE150163T1 (en) * | 1988-11-17 | 1997-03-15 | Tetronics Res & Dev Co Ltd | METHOD FOR MELTING SUBSTANCES |
| CN109255170B (en) * | 2018-08-28 | 2022-03-25 | 电子科技大学 | Arc voltage gradient modeling method based on magnetohydrodynamics simulation |
| CN112507583A (en) * | 2020-11-24 | 2021-03-16 | 深圳技术大学 | SF (sulfur hexafluoride)6Finite element simulation analysis method for arc plasma multi-physical field in opening process of breaker contact |
| CN112528527B (en) * | 2021-02-09 | 2021-08-24 | 北京科技大学 | Device and method for simulating and analyzing arc plasma |
-
2021
- 2021-04-20 CN CN202110421377.4A patent/CN112989674B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN112989674A (en) | 2021-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shakeri et al. | A finite volume spectral element method for solving magnetohydrodynamic (MHD) equations | |
| Dehghan et al. | Meshless local boundary integral equation (LBIE) method for the unsteady magnetohydrodynamic (MHD) flow in rectangular and circular pipes | |
| Piau et al. | Models of solar surface dynamics: impact on eigenfrequencies and radius | |
| Wu et al. | An efficient finite volume method for electric field–space charge coupled problems | |
| Giesecke et al. | Electromagnetic induction in non-uniform domains | |
| CN112528527B (en) | Device and method for simulating and analyzing arc plasma | |
| Zhang et al. | The variational multiscale element free Galerkin method for MHD flows at high Hartmann numbers | |
| CN105160117A (en) | Simulation method for plasma generator | |
| Poulos | Model for the operation of an emissive cathode in a large magnetized-plasma | |
| Makinde | Effect of arbitrary magnetic Reynolds number on MHD flows in convergent‐divergent channels | |
| Jabbari et al. | Turbulent reconnection of magnetic bipoles in stratified turbulence | |
| Sur | Decaying turbulence and magnetic fields in galaxy clusters | |
| Terentev et al. | Frequency dependencies of the characteristics of an inductively coupled radiofrequency discharge at reduced pressure | |
| CN112989674B (en) | Simulation analysis method and simulation analysis device for arc plasma characteristics | |
| Girault et al. | An algorithm for the computation of multiple Hopf bifurcation points based on Padé approximants | |
| JP2009174028A (en) | Microwave plasma analysis program | |
| Tsivilskiy et al. | An experimentally verified three-dimensional non-stationary fluid model of unloaded atmospheric pressure inductively coupled plasmas | |
| Prasath et al. | Dynamo action by turbulence in absolute equilibrium | |
| Chowdhury et al. | Influence of surface oxidation on electric potential measurements in MHD liquid metal flows | |
| Louai et al. | Implementation of an efficient element-free Galerkin method for electromagnetic computation | |
| Singh et al. | Magnetohydrodynamic mixed convection flow of liquid metals in a vertical channel: A stability analysis | |
| Arslan et al. | Exact and FDM solutions of 1D MHD flow between parallel electrically conducting and slipping plates | |
| Kawczynski et al. | Characterization of the lid-driven cavity magnetohydrodynamic flow at finite magnetic Reynolds numbers using far-field magnetic boundary conditions | |
| Icenhour et al. | Multi-physics object oriented simulation environment (moose) | |
| Guo et al. | Radial basis function interpolation supplemented lattice Boltzmann method for electroosmotic flows in microchannel |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |