Zero cell porosity in v5.0 => modify array c_solid_flag
Posted: Tue Oct 30, 2018 2:18 pm
Dear Users, Developers,
After some trial and error, I was able to use porosity in Code_Saturne 5.0. If you set the cell porosity to zero in some of the cells (= no fluid), then you will have to modify the array c_solid_flag. It is located in cs_glob_mesh_quantities. It should be zero where fluid exists and one where no fluid exists.
Failure to do so will most likely produce NaNs. They originate from the LU factorization (_fact_lu33 in alge/cs_sles.it).
With the hope that it will be helpful,
Cédric
EDIT: I just noticed that there is also a Fortran array isolid / isolid_0. It is updated in caltri only when iporos=3 and I was using iporos=1. I am not sure how redundant they are.
After some trial and error, I was able to use porosity in Code_Saturne 5.0. If you set the cell porosity to zero in some of the cells (= no fluid), then you will have to modify the array c_solid_flag. It is located in cs_glob_mesh_quantities. It should be zero where fluid exists and one where no fluid exists.
Failure to do so will most likely produce NaNs. They originate from the LU factorization (_fact_lu33 in alge/cs_sles.it).
With the hope that it will be helpful,
Cédric
EDIT: I just noticed that there is also a Fortran array isolid / isolid_0. It is updated in caltri only when iporos=3 and I was using iporos=1. I am not sure how redundant they are.