Segmentation fault

[Di LI]

Hello everyone,

I am doing a test with an external reduced electric field file. The electron density is calculated self-consistently

However,  a problem occurs as follows:

ZDPlasKin INFO: set autosave in QTplaskin format ON ...

          TIME

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:

#0  0xffffffff

#1  0xffffffff

#2  0xffffffff

#3  0xffffffff

#4  0xffffffff

#5  0xffffffff

#6  0xffffffff

#7  0xffffffff

#8  0xffffffff

#9  0xffffffff

#10  0xffffffff

#11  0xffffffff

#12  0xffffffff

#13  0xffffffff

#14  0xffffffff

#15  0xffffffff

#16  0xffffffff

#17  0xffffffff

#18  0xffffffff

#19  0xffffffff

#20  0xffffffff

#21  0xffffffff

I cannot find the problem and hope you can help me.

Thanks

[Prof. Maikel Y. Ballester]

Hello,

As Francisco C Parra-Rojas mentioned before, I also think this is related to the number of species used in Bolsig section, which has a maximum of 32. I recommend checking this.

Best,

Maikel

---

Prof. M.Y. Ballester Furones
Professor de Física
Universidade Federal de Juiz de Fora
www.ufjf.br/lfmt
Eu apoio a Cientista Empoderada
https://cientistaempoderada.com.br/

--
You received this message because you are subscribed to the Google Groups "ZDPlasKin" group.
To unsubscribe from this group and stop receiving emails from it, send an email to zdplaskin+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/zdplaskin/59982335-51fc-4824-be2f-21f290ac4decn%40googlegroups.com.

[Kevin Veer]

Dear all,

Based on the ZDPlasKin change log of version 2.0 this is not the case:

CHANGES in 2.0

- Latest BOLSIG+ ver. 03/2016 solver is integrated.
- Maximum number of species and individual electron cross sections loading
  to BOLSIG+ are increased to 200 and 4000, respectively.
- Run-time generated log file with the list of all loaded electron processes
  is available after each run. This might be used for checking links between
  reactions and electron cross sections.

And I don't think one would easily reach this limit.

The problem can be ill defined initial conditions (density, Temperature maybe)

or in your own code if it uses arrays.

Kind regards,

Kevin

Op dinsdag 13 april 2021 om 15:12:21 UTC+2 schreef Maikel:

[Di LI]

I have tried my best but still cannot find what is the problem. 

I am not sure if the external file should be written in special format or there is an incompatible issue of X64 version of MinGW.

[Di LI]

Hello,

My initial conditions are: 

double precision, parameter  :: gas_pressure     =  1.0d0,   &  ! pressure, atm

                                 gas_temperature = 353.0d0, &  ! temperature, K

                                 n2_ratio       = 0.5d0 ,  & ! N2 ratio

                                 o2_ratio       = 0.5d0 , &  ! o2 ratio  

 o2_density     = o2_ratio * 101325.0d0 * gas_pressure &

                                                 / gas_temperature &

                                                 / 1.38d-17, &  ! o2 density, cm-3

                                 n2_density     = n2_ratio * 101325.0d0 * gas_pressure &

                                                 / gas_temperature &

                                                 / 1.38d-17, &  ! n2 density, cm-3

My data_in file is:

Time_s Field_Td   

0 73.564079

0.0000001 62.2465284

0.0000002 67.9053037

0.0000003 67.9053037

0.0000004 62.2465284

0.0000005 62.2465284

0.0000006 67.9053037

0.0000007 73.564079

0.0000008 62.2465284

0.0000009 73.564079

0.000001 73.564079

0.0000011 79.2228543

0.0000012 90.540405

0.0000013 73.564079

0.0000014 73.564079

0.0000015 79.2228543

0.0000016 73.564079

0.0000017 79.2228543

0.0000018 84.8816296

0.0000019 101.857956

0.000002 50.9289778

0.0000021 84.8816296

.....

The maximum value is 204 Td

I think there is no problem with this part. Btw, I have 195 lines of data in data_in file, but the processor said 'found 97 points in data_in.dat", sot sure what's wrong here.

I still cannot find any solutions for this problem...

在2021年4月13日星期二 UTC+8 下午9:18:19<Kevin Veer> 写道：

[Di LI]

Hello,

Thank you for your reply, I've tried to reduce the number of species in Bolsig section to 26, but the problem still occurred. 

The 2.0 version can load more species, so I think the problem should be located somewhere else.

[Kevin Veer]

Then I would look into that code: reading the data_in.dat.

Which shoudl have nothing to do with ZDPK.

SIGSEGV error can occur when you try to access an index of an array but with the index that you are requesting being out of bounds.

You report that your file should be larger than what the code reads, so maybe you are out of bound due to the expected length?

Op woensdag 14 april 2021 om 07:46:50 UTC+2 schreef disherl...@gmail.com:

[Di LI]

Thank you for your advice, I tried to divide the value of reduced field by 10, then the code started to run successfully. 

It could be some values exceeded the boundary under the real reduced field.

But it doesn't make sense as the reduced field is calculated by the voltage file and also located within the typical range of DBD plasma.

I am now confused ......

[Di LI]

zdp also reported an error:

.At line 141 of file dbd0D.F90 (unit = 5, file = 'fort.5')

Fortran runtime error: End of file

When the reduced field is divided by 10...

在2021年4月14日星期三 UTC+8 下午5:16:53<Kevin Veer> 写道：

[Nicolas Minesi]

Hello,

Your field varies by approximatelly 5 Td by ns. This is a lot I think.

Nicolas

[Di LI]

Hello,

Thank you for your reply, the data_in is directly from the oscilloscope, I tried to simplify the data to avoid the repeating and spike points adjacently.  But the problem still occurred.

I assume it is due to some values exceeding the bounds, but I cannot figure it out.

[Kevin Veer]

Hi,

Sorry my last comment was not specifically about the value of E/N, anything between 1 and 1000 Td should be "fine", however for a filamentary plasma, I typically observe 200 Td maximum.

Large gradients can give problems.

My comment about being out of bounds was simply about arrays.

You mentioned something like your look up table (file) having 125 rows, but the program only reports 95 orso?

I was then wondering if you get an error because you might have hard coded to look for the 96th row+, while you only have 95?

Anyway assuming that that works:

- Such files must be interpolated (best to do that with your fortran code itself) such that you can take smaller time steps than the resolution of the file ofcourse (i recommend trying 1 ns, but try to see if it better with as little as 1 ps)

- I have never had much succes with imposing fixed E/N values, while electron densities are freely calculated.

- If the error is due to plasma conditions: look at and report n_E, E/N, T_E which you can all get as output from your model.

(E/N goes to BS+ to calculate the EEDF which returns T_E (average of the EEDF), and the EEDF in principle calculates k for e-impact reactions, determining for example ionization, thus determining your electron density. If the problem is with your imposed plasma conditions, it is probably there that you will find odd behaviour)

Op donderdag 15 april 2021 om 04:03:41 UTC+2 schreef disherl...@gmail.com:

[Di LI]

Hello, 

Thank you for your help.

I still cannot solve the Segmentation fault. I show you the reported results of Te, Ne,  and EN with the data_in file.

It seems zdplaskin can calculate 3 time steps at the beginning as the result shown. Then the process went wrong and showed Segmentation fault.

From the result, I didn't find any odd results or some values out of bounds.

By the way, now the 241 lines can all be found by zdplaskin solver.

[Di LI]

The data_in file :

Time_s      Field_Td    Electrons_cm-3

0         5.66 2.24E+10

2.00E-07 11.32 4.49E+09

4.00E-07 16.98 3.74E+09

6.00E-07 28.29 2.92E+09

8.00E-07 29.21 3.69E+09

1.00E-06 28.29 6.69E+09

1.20E-06 30.29 1.01E+10

1.40E-06 32.29 3.81E+09

1.60E-06 39.61 1.60E+09

1.80E-06 45.27 1.68E+09

2.00E-06 50.93 2.12E+09

0.0000022 45.27 1.54E+09

0.0000024 56.59 1.35E+09

0.0000026 56.59 2.02E+09

2.80E-06 56.59 1.46E+09

3.00E-06 56.59 1.35E+09

0.0000032 62.25 1.94E+09

0.0000034 67.91 1.59E+09

0.0000036 56.59 2.80E+09

3.80E-06 84.88 1.87E+09

4.00E-06 79.22 6.41E+08

4.20E-06 84.88 8.23E+08

0.0000044 90.54 2.31E+09

0.0000046 90.54 8.41E+08

0.0000048 107.52 1.06E+09

5.00E-06 96.2 1.45E+09

5.20E-06 101.86 7.48E+08

0.0000054 107.52 7.68E+08

0.0000056 113.18 6.73E+08

0.0000058 107.52 8.27E+08

6.00E-06 96.2 2.31E+09

6.20E-06 113.18 1.57E+09

6.40E-06 118.83 4.27E+08

0.0000066 113.18 1.29E+09

0.0000068 124.49 1.12E+09

0.000007 124.49 1.58E+09

7.20E-06 135.81 4.21E+08

7.40E-06 130.15 6.34E+08

7.60E-06 135.81 6.08E+08

0.0000078 135.81 5.14E+08

0.000008 147.13 7.34E+08

0.0000082 141.47 5.83E+08

8.40E-06 147.13 5.61E+08

8.60E-06 152.79 7.06E+08

0.0000088 147.13 6.47E+08

0.000009 158.45 4.41E+08

0.0000092 158.45 6.81E+08

9.40E-06 152.79 7.48E+08

9.60E-06 158.45 5.61E+08

9.80E-06 175.42 5.43E+08

0.00001     169.76 5.24E+08

0.0000102 169.76 4.86E+08

0.0000104 181.08 6.66E+08

1.06E-05 169.76 4.49E+08

1.08E-05 169.76 1.76E+09

0.000011 175.42 4.71E+08

0.0000112 192.4 3.96E+08

0.0000114 181.08 4.56E+08

1.16E-05 186.74 4.76E+08

1.18E-05 186.74 3.40E+08

1.20E-05 181.08 3.86E+08

0.0000122 186.74 3.74E+08

0.0000124 198.06 2.89E+08

0.0000126 198.06 3.21E+08

1.28E-05 203.72 3.12E+08

1.30E-05 203.72 2.80E+08

0.0000132 192.4 6.60E+08

0.0000134 192.4 1.32E+09

0.0000136 209.37 2.12E+08

1.38E-05 215.03 1.77E+08

1.40E-05 209.37 3.03E+08

1.42E-05 215.03 2.66E+08

0.0000144 215.03 2.36E+08

0.0000146 215.03 2.36E+08

0.0000148 220.69 6.04E+08

1.50E-05 220.69 2.59E+08

1.52E-05 215.03 3.84E+08

1.54E-05 226.35 2.80E+08

0.0000156 226.35 1.96E+08

0.0000158 226.35 3.09E+08

0.000016 237.67 1.87E+08

1.62E-05 237.67 1.60E+08

1.64E-05 232.01 2.46E+08

0.0000166 232.01 2.46E+08

0.0000168 237.67 2.14E+08

0.000017 232.01 1.92E+08

1.72E-05 243.33 2.35E+08

1.74E-05 232.01 1.29E+09

1.76E-05 237.67 3.74E+08

0.0000178 243.33 5.22E+07

0.000018 254.64 1.25E+08

0.0000182 243.33 1.83E+08

1.84E-05 237.67 2.14E+08

1.86E-05 243.33 3.13E+08

0.0000188 248.99 1.78E+08

0.000019 248.99 1.53E+08

0.0000192 248.99 1.02E+08

1.94E-05 248.99 1.02E+08

1.96E-05 248.99 1.27E+08

1.98E-05 254.64 7.48E+07

0.00002     254.64 9.97E+07

0.0000202 254.64 1.50E+08

0.0000204 254.64 1.25E+08

0.0000206 254.64 2.99E+08

0.0000208 260.3 9.76E+07

0.000021 260.3 9.76E+07

0.0000212 260.3 7.32E+07

0.0000214 254.64 9.97E+07

0.0000216 260.3 1.22E+08

0.0000218 265.96 4.54E+08

0.000022 254.64 4.99E+07

0.0000222 260.3 9.76E+07

0.0000224 265.96 7.16E+07

0.0000226 265.96 7.16E+07

0.0000228 260.3 7.32E+07

0.000023 260.3 7.32E+07

0.0000232 260.3 7.32E+07

0.0000234 260.3 7.32E+07

0.0000236 260.3 4.88E+07

0.0000238 265.96 4.77E+07

0.000024 265.96 4.77E+07

0.0000242 260.3 4.88E+07

0.0000244 260.3 2.44E+07

0.0000246 265.96 4.77E+07

0.0000248 260.3 2.44E+07

0.000025 265.96 2.39E+07

0.0000252 260.3 2.44E+07

0.0000254 265.96 2.39E+07

0.0000256 260.3 0.00E+00

0.0000258 254.64 0.00E+00

0.000026 254.64 2.49E+07

0.0000262 260.3 0.00E+00

0.0000264 254.64 0.00E+00

0.0000266 260.3 2.44E+07

0.0000268 265.96 2.39E+07

0.000027 260.3 2.44E+07

0.0000272 254.64 0.00E+00

0.0000274 260.3 2.44E+07

0.0000276 260.3 0.00E+00

0.0000278 254.64 0.00E+00

0.000028 254.64 2.49E+07

0.0000282 254.64 2.49E+07

0.0000284 254.64 2.49E+07

0.0000286 248.99 2.55E+07

0.0000288 254.64 2.49E+07

0.000029 254.64 2.49E+07

0.0000292 248.99 2.55E+07

0.0000294 248.99 2.55E+07

0.0000296 248.99 2.55E+07

0.0000298 243.33 5.22E+07

0.00003     243.33 2.61E+07

0.0000302 243.33 7.83E+07

0.0000304 243.33 7.83E+07

0.0000306 243.33 5.22E+07

0.0000308 237.67 5.34E+07

0.000031 237.67 8.01E+07

0.0000312 237.67 8.01E+07

0.0000314 232.01 8.21E+07

0.0000316 237.67 5.34E+07

0.0000318 226.35 8.41E+07

0.000032 237.67 8.01E+07

0.0000322 232.01 8.21E+07

0.0000324 226.35 8.41E+07

0.0000326 232.01 1.09E+08

0.0000328 226.35 1.12E+08

0.000033 220.69 1.15E+08

0.0000332 220.69 1.44E+08

0.0000334 226.35 1.12E+08

0.0000336 215.03 1.48E+08

0.0000338 220.69 1.15E+08

0.000034 215.03 1.18E+08

0.0000342 209.37 1.52E+08

0.0000344 215.03 1.18E+08

0.0000346 209.37 1.21E+08

0.0000348 203.72 1.56E+08

0.000035 209.37 1.21E+08

0.0000352 198.06 1.92E+08

0.0000354 192.4 1.65E+08

0.0000356 192.4 1.98E+08

0.0000358 198.06 1.92E+08

0.000036 192.4 1.98E+08

0.0000362 192.4 1.98E+08

0.0000364 192.4 1.65E+08

0.0000366 192.4 1.98E+08

0.0000368 181.08 2.10E+08

0.000037 175.42 2.17E+08

0.0000372 181.08 2.10E+08

0.0000374 181.08 2.45E+08

0.0000376 181.08 2.10E+08

0.0000378 169.76 2.24E+08

0.000038 169.76 2.24E+08

0.0000382 164.1 2.71E+08

0.0000384 158.45 2.80E+08

0.0000386 152.79 3.32E+08

0.0000388 158.45 2.80E+08

0.000039 152.79 2.91E+08

0.0000392 147.13 3.02E+08

0.0000394 152.79 2.91E+08

0.0000396 141.47 3.14E+08

0.0000398 147.13 3.45E+08

0.00004     135.81 3.27E+08

0.0000402 141.47 403913647.4

0.0000404 135.81 327244853.2

0.0000406 130.15 390254731.8

0.0000408 124.49 458992781.1

0.000041 124.49 407993583.2

0.0000412 113.18 504892059.2

0.0000414 118.83 480849580.2

0.0000416 118.83 534277311.4

0.0000418 107.52 531465325.5

0.000042 118.83 480849580.2

0.0000422 101.86 623323529.9

0.0000424 107.52 472413622.7

0.0000426 107.52 531465325.5

0.0000428 90.54 631115074.1

0.000043 84.88 673189412.3

0.0000432 90.54 701238971.2

0.0000434 84.88 747988235.9

0.0000436 90.54 631115074.1

0.0000438 79.22 801415967.1

0.000044 73.56 863063349.1

0.0000442 73.56 776757014.2

0.0000444 67.91 934985294.9

0.0000446 67.91 934985294.9

0.0000448 62.25 917985562.3

0.000045 56.59 1121982354

0.0000452 50.93 1371311766

0.0000454 50.93 1121982354

0.0000456 45.27 1402477942

0.0000458 45.27 1542725737

0.000046 50.93 2617958826

0.0000462 33.95 2243964708

0.0000464 39.61 1763115128

0.0000466 33.95 2991952944

0.0000468 22.64 2804955885

0.000047 22.64 5329416181

0.0000472 28.29 2243964708

0.0000474 16.98 10845829421

0.0000476 11.32 8414867654

0.0000478 11.32 5048920592

0.000048 5.66 13463788246

[Di LI]

[Nicolas Minesi]

Hello,

Following a parallel conversation on the forum, what happens if you just set the reduced field to 50 Td? Or 100 Td?

If this work, you're simply not providing correctly E/N to the solver.

Best,
Nicolas