Boundary conditions for PIC simulation

Dear Jeff and all,

I will come back to this topic instead of creating a new one. This is also in combination with also another topic I am dealing with, which involves a complex particle-level code:

https://mfix.netl.doe.gov/forum/t/implementation-of-a-particle-level-model/

Regarding the BC setting for this specific case (ignoring for the moment the particle-level). I set up the case as indicated previously, where the two black boxes covers part of the tube and then I have a mass inflow of fresh biomass close to the inclined wall.

The problem I am observing is that the case runs smoothly but then I am getting an error, T_g < TMIN after around 0.2 s. I am monitoring the minimum T in the system and it seems ok until I reach this problem and it jumps from 293 K to TMIN (I don’t know why). In fact, it also occurs when removing the chemical reactions out of the system (setting A = 0 in all of them). So it is not a matter of the reactions. What could I be doing wrong? Is it possible that the mesh is creating this disturbance? Why is the minimum T jumping to TMIN and not any other value?

I am aware that some parcels are seeded out of the domain and that the SW for PIC is very big, but in view of the case running for some time and then finding this error I believe this is unrelated. Any piece of advice on how to improve the model is highly appreciated.

Please see attached the .stl and the .mfx files that run the model.

Thank you for your kind help,
Eduardo

geometria_coarser.stl (1.9 MB)
test_pequeno.mfx (37.1 KB)