[PLNTRF]

Dear all,

I am trying to simulate free gas expansion with non-Maxwellian initial distribution function f(p, r, t=0) with dsmcFoam. I assume that cloud of particles is injected once and then expand.

I set initial conditions for distributions rhoN, rhoM, ect. in 0.orig folder. I have issues with dsmcInitialise procedure, because dsmcInitialiseDictionary requires temperature of particles. I consider temperature as a moment of distribution function f(p, r, t). What kind of distribution function determines temperature in dsmcInitialiseDictionary? How can I set particle's temperature that corresponds to my initial distribution function? Is dsmcFoam solver suitable for that kind of problem?

Thank you.