I am using the PIC method with the version of 19.2.0(newest)
First, I calculate the case on my local computer for a while
Then, I transfer the case to the cluster and restart it, but there is an error occurred
Dear Dr. Musser
It does works and thanks for your suggestion.
However, I have got another problem with chemical reaction.
1, When I enable the species equation of gas (species_gas(0)= true) without any chemical equations, in my mind, the species equations should not be solved, but the residual of the gas species still get some values rather than zero
If you enable the species equations, then they are solved. If there are no chemical reactions or phase changes are defined, (i.e., usr_rates.f and usr_rates_des.f have not been modified), then you are solving a basic scalar transport equation.
By default, MFIX initializes the domain to the first species and sets all other species to zero. You can change this by modifying the gas composition in the initial conditions. In you case, all species except the first have a value of zero everywhere. The result is they have no-error (as there is no value to transport). However, the non-zero first species has some error given that it is non-zero, but the amount should be rather small. If it isn’t, it could indicate that your mesh quality is poor.
Check to see if there are errors in your LOG file. If you ran in parallel with MPI, there could be multiple LOG files. The route shown in the back trace, check_data_30.f assesses the value of field variables and stops the simulation if unphysical values are detected.
Many thanks!
I have refined the mesh but still the problem
The problem occurred only five timesteps without reactions. I guess the reason is not due to the unbalance of the mass fraction. Some other properties lead to this but I am not sure.
Besides, the value of the fraction is one, why it says it beyonds the bound