[naoki]

I've used refineMesh utility of OpenFOAM-1.4, and it has worked well.
However, when I tried refineMesh utility using v.1.4.1, it didn't work well.

In system/refineMeshDict, "useHexTopology" is true, and "geometricCut" is false.

The error message I got is the following:
^^^^^^^

Exec : refineMesh . . -dict
Date : Dec 05 2007
Time : 17:11:33
Host : hoge
PID : 10800
Root : /home/ohnishi/test/refineTest
Case : .
Nprocs : 1
Create time

Create polyMesh for time = 0

Mesh edge statistics:
x aligned : number:4620 minLen:0.5 maxLen:0.5
y aligned : number:4410 minLen:0.5 maxLen:0.5
z aligned : number:4620 minLen:0.5 maxLen:0.5
other : number:0 minLen:1e+15 maxLen:-1e+15

Refining according to refineMeshDict

Read 2000 cells from cellSet "constant/polyMesh/sets/c0"

Writing refined cells (2537) to cellSet "0.005/polyMesh/sets/refinedCells"

#0 Foam::error::printStack(Foam:stream&) in "/home/ohnishi/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/home/ohnishi/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 main in "/home/ohnishi/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/refineMe sh"
#4 __libc_start_main in "/lib64/libc.so.6"
#5 Foam::regIOobject::readIfModified() in "/home/ohnishi/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/refineMe sh"
^^^^

Then, now I try refineMesh with "useHexTopology false", and it looks working fine.
So, I wonder there are any bugs in somewhere concerned with "useHexTopylogy".

regards,

Naoki

[mattijs]

Can you send me a simple case where this happens?

[mattijs]

Dear Naoki Onishi,

thanks for the testcase. The bug was just in the picking up of the cells to refine. Replace 1.4.1 dynamicMesh/meshCut/meshModifiers/multiDirRefinement/multiDirRefinement.C with attached version

multiDirRefinement.C

[jaswi]

Dear Mattijs

Good Morning (guess we are in same time zone)

I am having trouble with running nozzle2D test case and problem lies where the Allrun script executes the refineMesh . nozzle2D step.

The error says:

<pre>
Exec : refineMesh . nozzleFlow2D -dict
Root : /home/openfoam/OpenFOAM/OpenFOAM-1.4.1/tutorials/lesInterFoam
Case : nozzleFlow2D
Nprocs : 1
Create time

Create polyMesh for time = 0

Mesh edge statistics:
x aligned : number:4987 minLen:4e-05 maxLen:4.0039947e-05
y aligned : number:15200 minLen:1.1640181e-06 maxLen:5.3384144e-05
z aligned : number:7600 minLen:6.2039102e-07 maxLen:0.00013971502
other : number:10088 minLen:4.0040048e-05 maxLen:4.0320766e-05

Refining according to refineMeshDict

Read 3537 cells from cellSet "constant/polyMesh/sets/c0"

Global Coordinate system:
normal : (0 0 1)
tan1 : (1 0 0)
tan2 : (0 1 0)

--> FOAM FATAL ERROR : Cannot move points: size of given point list smaller than the number of active points

From function primitiveMesh::movePoints(const pointField& newPoints, const pointField& oldPoints)
in file meshes/primitiveMesh/primitiveMesh.C at line 251.

FOAM aborting
</pre>

Forum search lead me to this post where you have posted the corrected version of dynamicMesh/meshCut/meshModifiers/multiDirRefinement/multiDirRefinement.C

I did as suggested but recompiling the library gives the following error:

<pre>
openfoam@nari:~/OpenFOAM/OpenFOAM-1.4.1/src/dynamicMesh> wmake libso
Making dependency list for source file meshCut/meshModifiers/multiDirRefinement/multiDirRefinement.C
SOURCE=meshCut/meshModifiers/multiDirRefinement/multiDirRefinement.C ; g++ -m64 -Dlinux64 -DDP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -march=opteron -O3 -DNoRepository -ftemplate-depth-40 -I/home/openfoam/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude -I/home/openfoam/OpenFOAM/OpenFOAM-1.4.1/src/meshTools/lnInclude -I/home/openfoam/OpenFOAM/OpenFOAM-1.4.1/src/triSurface/lnInclude -IlnInclude -I. -I/home/openfoam/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/multiDirRefinement.o
meshCut/meshModifiers/multiDirRefinement/multiDirRefinement.C: In member function 'void Foam::multiDirRefinement::refineHex8(Foam::polyMes h&, const Foam::labelList&, bool)':
meshCut/meshModifiers/multiDirRefinement/multiDirRefinement.C:351: error: 'class Foam::polyTopoChange' has no member named 'changeMesh'
make: *** [Make/linux64GccDPOpt/multiDirRefinement.o] Error 1

</pre>

Please , any clues what I am doing wrong

With Best Regards
Jaswinder

[jaswi]

Hallo Mattijs

Update to the above posted message:

I am running the development version.

I just tried the same tutorial in the official version and it went through. No error messages as such from the refineMesh utility The log.refineMesh is clean and case is running.

Sorry for the incorrect post. Perhaps there is some incompatibilty between the src/dynamicMesh folders of official and development version. (Just a guess, apologies if that is not the case).

Please do not look into this matter for now. I will run diff over the related classes and report the outcome.

Thanks

Best Regards
Jaswinder

[mattijs]

Thanks for letting me know. The handling of retired points is probably not consistent in 1.4.1. If you don't care about them you can usually trim them by running subsetMesh with a cellset containing all cells.