Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

DBD plasma actuator model

Please login with a confirmed email address before reporting spam

Hello everyone,
I am trying to study gas discharge at atmospheric pressure using a dielectric barrier discharge plasma actuator (2D geometry) consisting of two electrodes, one exposed to air gas at atmospheric pressure and the other covered by a dielectric material. Air plasma kinetic model was modelled by air cross-sections from file, with N2, O2, N2_1p, O2_1p as species. Volumetric and surface reactions were also considered. N2 species has the largest mass fraction and it was chosen to fulfil the mass constraint. Initial mass fraction of O2 species is equal to 0.22. Initial number density of O2 ions is zero. Electric potential BC is applied on the exposed electrode (Vrf = 5e3[V]*sin(2*pi*1e3*t[1/s])), Ground BC on the other electrode.
I set a time-dependent study step from 0 to 10 ms, but I get the error:

Repeated error test failures. May have reached a singularity.
Time : 1.985019325484047e-005
Undefined value found.
Undefined value found in the stiffness matrix.
There are 2 equations giving NaN/Inf in the matrix rows for the variable mod1.wN2_1p.
at coordinates: (0.00100937,0.00120891), (0.00102875,0.0012092), ...
and similarly for the degrees of freedom, NaN/Inf in the matrix columns.
Last time step is not converged.

The same model was successful solved by applying an electric potential with same freq but lower amplitude (1000V), so I think that the problem is the too much high amplitude of the electric potential Vrf.
I tried to increase the Vrf’s amplitude gradually by using previous solutions (obtained with lower amplitude) as initial values of dependent variables (both variables solved for and variables not solved for), but I tried the same error with amplitudes > 1000 V.

Can someone help me?

Chiara D.P.

3 Replies Last Post Jan 29, 2017, 6:53 a.m. EST
Luke Gritter Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 10 years ago Nov 24, 2014, 10:04 a.m. EST
Chiara,

When the error message includes the comment "May have reached a singularity", it is often the case that a dependent variable is tending toward infinity or negative infinity on one of the node points. The software should save the solution at the time of failure, so you can plot the dependent variable that is referenced in the error message (wN2_1p) at this time and check whether a singularity appears anywhere in the solution. If not, check the other dependent variables. If you see a singularity, try refining the mesh in that location and running the model again.

Convergence errors are not uncommon for this type of analysis, and an overly coarse mesh is often the cause. Generally, higher plasma currents require finer meshes. It also appears that you may be using the linear equation formulation. If so, you may want to try the log formulation, which is often more stable.

--
Luke Gritter
AltaSim Technologies
Chiara, When the error message includes the comment "May have reached a singularity", it is often the case that a dependent variable is tending toward infinity or negative infinity on one of the node points. The software should save the solution at the time of failure, so you can plot the dependent variable that is referenced in the error message (wN2_1p) at this time and check whether a singularity appears anywhere in the solution. If not, check the other dependent variables. If you see a singularity, try refining the mesh in that location and running the model again. Convergence errors are not uncommon for this type of analysis, and an overly coarse mesh is often the cause. Generally, higher plasma currents require finer meshes. It also appears that you may be using the linear equation formulation. If so, you may want to try the log formulation, which is often more stable. -- Luke Gritter AltaSim Technologies

Please login with a confirmed email address before reporting spam

Posted: 10 years ago Nov 27, 2014, 8:31 a.m. EST
Hi Luke,
I solved my model by using extremely fine custom meshing.
I thank you so much for your help!

Chiara
Hi Luke, I solved my model by using extremely fine custom meshing. I thank you so much for your help! Chiara

Please login with a confirmed email address before reporting spam

Posted: 8 years ago Jan 29, 2017, 6:53 a.m. EST
I'm trying to simulate a DBD plasma actuator using the DC discharge module.A sinusoidal signal of 1kV and 10 kHz is applied to the electrode,the simulation is carried for argon first and i have a problem of convergence if any one can please help me.
I'm trying to simulate a DBD plasma actuator using the DC discharge module.A sinusoidal signal of 1kV and 10 kHz is applied to the electrode,the simulation is carried for argon first and i have a problem of convergence if any one can please help me.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.