CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Bugs

Cannot refineMesh with useHexTopology true ver141

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   December 5, 2007, 05:03
Default I've used refineMesh utility o
  #1
New Member
 
Naoki Onishi
Join Date: Mar 2009
Location: Tokyo, Japan
Posts: 4
Rep Power: 7
naoki is on a distinguished road
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
naoki is offline   Reply With Quote

Old   December 6, 2007, 05:47
Default Can you send me a simple case
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
Can you send me a simple case where this happens?
mattijs is offline   Reply With Quote

Old   December 7, 2007, 06:53
Default Dear Naoki Onishi, thanks f
  #3
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
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
mattijs is offline   Reply With Quote

Old   June 11, 2008, 04:26
Default Dear Mattijs Good Morning (
  #4
Senior Member
 
Join Date: Mar 2009
Posts: 248
Rep Power: 8
jaswi is on a distinguished road
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 is offline   Reply With Quote

Old   June 11, 2008, 04:50
Default Hallo Mattijs Update to the
  #5
Senior Member
 
Join Date: Mar 2009
Posts: 248
Rep Power: 8
jaswi is on a distinguished road
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
jaswi is offline   Reply With Quote

Old   June 11, 2008, 05:07
Default Thanks for letting me know. Th
  #6
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
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.
mattijs is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
BlockMesh cellSet refineMesh mattijs OpenFOAM Mesh Utilities 34 August 19, 2014 05:54
RefineMesh tetrahedral markc OpenFOAM Mesh Utilities 3 December 7, 2010 10:16
Is there a simple way to make refineMesh write the cellLevel chtrapp OpenFOAM Pre-Processing 0 February 16, 2009 11:25
RefineMesh warning mgz1985 OpenFOAM Native Meshers: blockMesh 1 August 29, 2008 09:45
Using refineMesh matteo_gautero OpenFOAM Mesh Utilities 0 February 11, 2008 10:07


All times are GMT -4. The time now is 01:21.