Zero cell porosity in v5.0 => modify array c_solid_flag

Questions and remarks about Code_Saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Post Reply
Posts: 20
Joined: Fri Apr 08, 2016 2:19 pm

Zero cell porosity in v5.0 => modify array c_solid_flag

Post by C.FLAG. » 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/

With the hope that it will be helpful,

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.

Post Reply