Growth submodel based on cell pressure #194
Replies: 1 comment 1 reply
-
I can only suggest a strategy to troubleshoot the problem: save or print every variable in every cell and make sure that everything is indeed constant first. How do you recover the pressure field from the growth kernel function? |
Beta Was this translation helpful? Give feedback.
1 reply
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
-
I wrote a populationBalanceSubModel based on grid pressure 'p_' , a constant density 'rho_' and a constant pressure 'pV_'. At the end of the Kg function
return Cg_.value()sqrt(2.0/(3.0rho_mag(pV_-p_[celli])))(pV_-p_[celli])
So the Kg function needs to be provided an additional celli parameter. I have added this parameter to the Kg function of the rest of the submodel, added the celli parameter to the Kg function of the parent class of growthModel, and added the celli argument to all calls to the Kg function. The compilation was successful, the pureGrowth example under pbeFoam was used to verify my increased growth model and I get correct result. However, when I used case with more cells and set the constant pressure field and 0 velocity, I found that the particle growth rate was not the same on each cell. Does anyone know why this is so? Or I can provide more detailed information.
Beta Was this translation helpful? Give feedback.
All reactions