Temprature parameter in input.deck
naokikatsura opened this issue · 2 comments
Dear All,
I am currently using EPOCH2D version 4.19.2 to simulate the expansion process of high-density (n0 = 10^23 m^-3) and high-temperature (T0 = 1 keV) plasma in an external divergent magnetic field. Before conducting simulations involving the magnetic field, I performed two test calculations of plasma expansion without the magnetic field using identical parameters, except for the temperature settings. The temperature in the first simulation is set by:
number_density = if(sqrt(xx+yy) lt r0, n0*(1-sqrt(xx+yy)/r0), 0.0)
temp_ev = if(sqrt(xx+yy) lt r0, T0, 0.0)
In the second simulation, it is set by:
number_density = if(sqrt(xx+yy) lt r0, n0*(1-sqrt(xx+yy)/r0), 0.0)
temprature_ev = if(sqrt(xx+yy) lt r0, T0, 0.0).
The spatial distribution of electron number density shows differing characteristics between the two simulations (please see the attached figures).
The distribution set by temp_ev demonstrates spreading over time, whereas in the second simulation, the electrons remain at their initial positions. I attach the input.deck of this simulation just in case.
input_deck.txt
Could anyone explain why these differences occur?
Thank you for your assistance.
Regards,
Naoki Katsura
It's hard to tell without your full second input deck (the one which shows no expansion), but I suspect the issue you are seeing is due to the spelling of temprature_ev
being incorrect. Try temperature_ev
.
Oh, I didn't notice that 'temperature' was mistyped as 'temperature.' That's so embarrassing!
Thank you for your comment.