setFields after snappyHexMesh

mo.houssami · May 11, 2015, 03:41

Dear foamers, I have been trying to do the following

blockMesh
surfaceFeatureExtract
cp system/snappyHexMeshDict.stl system/snappyHexMeshDict
snappyHexMesh -overwrite //no error
topoSet
createPatch -overwrite
extrudeMesh // going from 3D to 2D
checkMesh // mesh ok!
setFields
decomposePar

I noticed that even though setFields doesn't show any error message, I find that there's nothing changed in 0/

I also get this error message after decomposePar:
--> FOAM FATAL IO ERROR:
size 2369124 is not equal to the given value of 82122

file: /scratch/mhoussami/job/foam1/PostFMG/2Dmesh/0/ccx from line 18 to line 2377264.

From function Field<Type>::Field(const word& keyword, const dictionary&, const label)
in file /home/mhoussami/OpenFOAM/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude/Field.C at line 236.

FOAM exiting

I guess it's because after doing extrudeMesh, the number of cells is not the same anymore.

Why doesn't setFields do anything? how can I decomposePar?

Thank you

nimasam · May 12, 2015, 12:31

remove Cxx and Cyy and so on

mo.houssami · May 12, 2015, 12:47

Hi NimaSam,

Thank you for your quick answer, actually I have done that before, but I wasn't sure it was the right way.

I had to remove cc*, pointLevel.gz and cellLevel.gz in 0/ and then decomposePar runs with no error message.

The only problem remaining is that I still don't see the effect of setFields. The log file from setFields shows no problem at all, but clearly no region values are added in 0/

Any experience with that?
Thanks

Mo

nimasam · May 12, 2015, 16:42

Play with setFields, order of commands is correct

mo.houssami · May 13, 2015, 11:44

Yes, my problem was in one of the coordinates (zmin zmax) which apparently was too small.

Now it works, thanks!

May 11, 2015, 03:41	setFields after snappyHexMesh	#1
mo.houssami New Member mo houssami Join Date: May 2011 Posts: 14 Rep Power: 14	Dear foamers, I have been trying to do the following blockMesh surfaceFeatureExtract cp system/snappyHexMeshDict.stl system/snappyHexMeshDict snappyHexMesh -overwrite //no error topoSet createPatch -overwrite extrudeMesh // going from 3D to 2D checkMesh // mesh ok! setFields decomposePar I noticed that even though setFields doesn't show any error message, I find that there's nothing changed in 0/ I also get this error message after decomposePar: --> FOAM FATAL IO ERROR: size 2369124 is not equal to the given value of 82122 file: /scratch/mhoussami/job/foam1/PostFMG/2Dmesh/0/ccx from line 18 to line 2377264. From function Field<Type>::Field(const word& keyword, const dictionary&, const label) in file /home/mhoussami/OpenFOAM/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude/Field.C at line 236. FOAM exiting I guess it's because after doing extrudeMesh, the number of cells is not the same anymore. Why doesn't setFields do anything? how can I decomposePar? Thank you

May 12, 2015, 12:31		#2
nimasam Senior Member Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,266 Blog Entries: 1 Rep Power: 24	remove Cxx and Cyy and so on __________________ My Personal Website (http://nimasamkhaniani.ir/) Telegram channel (https://t.me/cfd_foam)

May 12, 2015, 16:42		#4
nimasam Senior Member Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,266 Blog Entries: 1 Rep Power: 24	Play with setFields, order of commands is correct __________________ My Personal Website (http://nimasamkhaniani.ir/) Telegram channel (https://t.me/cfd_foam)

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[snappyHexMesh] snappyhexmesh doesn't creat mesh in parallel issue?	klausb	OpenFOAM Meshing & Mesh Conversion	1	March 7, 2015 11:55
[snappyHexMesh] Error in SnappyHexMesh	gooya_kabir	OpenFOAM Meshing & Mesh Conversion	2	October 23, 2013 04:41
Strange Results With snappyHexMesh	calebamiles	OpenFOAM Running, Solving & CFD	0	August 14, 2011 16:02
[snappyHexMesh] stitchMesh and snappyHexMesh	gdbaldw	OpenFOAM Meshing & Mesh Conversion	0	December 23, 2009 02:09
[snappyHexMesh] SnappyHexMesh not generate mesh first time	mavimo	OpenFOAM Meshing & Mesh Conversion	4	August 26, 2008 07:08

May 12, 2015, 12:47		#3
mo.houssami New Member mo houssami Join Date: May 2011 Posts: 14 Rep Power: 14	Hi NimaSam, Thank you for your quick answer, actually I have done that before, but I wasn't sure it was the right way. I had to remove cc*, pointLevel.gz and cellLevel.gz in 0/ and then decomposePar runs with no error message. The only problem remaining is that I still don't see the effect of setFields. The log file from setFields shows no problem at all, but clearly no region values are added in 0/ Any experience with that? Thanks Mo

May 13, 2015, 11:44		#5
mo.houssami New Member mo houssami Join Date: May 2011 Posts: 14 Rep Power: 14	Yes, my problem was in one of the coordinates (zmin zmax) which apparently was too small. Now it works, thanks!