I'm modelling urban airflow using buoyantsimplefoam. We are using actual terrain instead of flat terrain. I created mesh using the following command
mpirun -np 64 snappyHexMesh -parallel -overwrite
Log snappyhexmesh ends with
Writing mesh to time constant
Wrote mesh in = 350.41 s.
Mesh snapped in = 2558.65 s.
Checking final mesh ...
Checking faces in error :
non-orthogonality > 60 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-15 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.02 : 0
faces on cells with determinant < 0.001 : 0
Finished meshing without any errors
Finished meshing in = 5887.69 s.
Finalising parallel run
After reconstructParMesh and while checking mesh
I get one quality failed
Overall domain bounding box (0 0 -5.739853e-14) (5500 5500 1000)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (1.53412e-16 -3.11233e-15 1.939486e-13) OK.
Max cell openness = 4.438246e-16 OK.
Max aspect ratio = 9.576651 OK.
Minimum face area = 0.007753621. Maximum face area = 421.7486. Face area magnitudes OK.
Min volume = 0.05158672. Max volume = 6093.245. Total volume = 3.007628e+10. Cell volumes OK.
Mesh non-orthogonality Max: 60.00908 average: 3.920135
Non-orthogonality check OK.
Face pyramids OK.
***Max skewness = 7.219928, 10 highly skew faces detected which may impair the quality of the results
<<Writing 10 skew faces to set skewFaces
Coupled point location match (average 0) OK.
Failed 1 mesh checks.
I'm not sure why mesh fails after reconstructing
When I run the simulation with single core the simulation starts without error
But while decompose again and run parallel my decomposition fails with errors like
--> FOAM FATAL IO ERROR:
error in IOstream "/home/users/nus/ijaz01/PinnacleCaseStudy1_buoyantsimplefoam/20230920/case2_reducedmesh/simulation1/processor16/constant/polyMesh/faces" for operation Foam::Ostream& Foam::operator<<(Foam::Ostream&, int64_t)
file: /home/users/nus/ijaz01/PinnacleCaseStudy1_buoyantsimplefoam/20230920/case2_reducedmesh/simulation1/processor16/constant/polyMesh/faces at line 2867957.
 --> FOAM FATAL ERROR:
 Cannot find file "points" in directory "polyMesh" in times "0" down to constant
my decomposePar dict is as below
n ( 2 2 2 );
n ( 1 1 1 );
roots ( );
I'm using openfoam 1912, is there any issue with decomposition and reconstructPar that causes error while trying to run the simulation in parallel?