CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (http://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   surfaceToPatch bug? (http://www.cfd-online.com/Forums/openfoam-bugs/69592-surfacetopatch-bug.html)

bruce October 28, 2009 07:07

surfaceToPatch bug?
 
Hello

I have strange bug in surfaceToPatch utility.

I have polyMesh with patch called "wall".

I have created faceSet of wall called "wall" even though it is not necessary.

I have stl file called "in.stl"

I have to create new polyPatch called "in" using stl file where it matches tolerence of 0.001.

So i use surfaceToPatch utility :

I use OpenFOAM v1.6.x

Here is the bug report:

[bruce@cfd4 patchCheck]$ surfaceToPatch in.stl -faceSet wall
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6-f802ff2d6c5a
Exec : surfaceToPatch in.stl -faceSet wall
Date : Oct 28 2009
Time : 16:34:46
Host : cfd4
PID : 12067
Case : /SNAPPY/01.full/patchCheck
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Reading surface from "in.stl" ...
Repatching only the faces in faceSet wall according to nearest surface triangle ...
All boundary faces further away than 0.001 of mesh bounding box (2729571.5 6066501 -50) (2766569.25 6103499 3230.00509) will keep their patch label ...
Before patching:
patch size
upperWall 25452
lowerWall 25452
wall 9234

Surface region 0 becomes boundary patch 0 with name stl_file
Created patch stl_file at 3


Could not find (2.73002e+06 6.08084e+06 919.998) in octree.#0 Foam::error:rintStack(Foam::Ostream&) in "/volume/cfd4_shared/2009-ct-PI-004//OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/volume/cfd4_shared/2009-ct-PI-004//OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::octreeDataFaceList::getSampleType(Foam:ctree <Foam:ctreeDataFaceList> const&, Foam::Vector<double> const&) const in "/volume/cfd4_shared/2009-ct-PI-004//OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libdynamicMesh.so"
#3 Foam::treeNode<Foam::octreeDataFaceList>::setSubNo deType(int, Foam:ctree<Foam::octreeDataFaceList>&, Foam::octreeDataFaceList const&) in "/volume/cfd4_shared/2009-ct-PI-004//OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libdynamicMesh.so"
#4 Foam::octree<Foam::octreeDataFaceList>::octree(Foa m::treeBoundBox const&, Foam::octreeDataFaceList const&, int, double, double) in "/volume/cfd4_shared/2009-ct-PI-004//OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libdynamicMesh.so"
#5 Foam::boundaryMesh::getNearest(Foam::primitiveMesh const&, Foam::Vector<double> const&) const in "/volume/cfd4_shared/2009-ct-PI-004//OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libdynamicMesh.so"
#6 main in "/volume/cfd4_shared/2009-ct-PI-004/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/surfaceToPatch"
#7 __libc_start_main in "/lib64/libc.so.6"
#8 __gxx_personality_v0 in "/volume/cfd4_shared/2009-ct-PI-004/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/surfaceToPatch"


From function octreeDataFaceList::getSampleType(octree<octreeDat aFaceList>&, const point&)
in file boundaryMesh/octreeDataFaceList.C at line 153.

FOAM aborting

Aborted

Can anyone help me to fix the problem?

Kind Regards

mattijs November 2, 2009 07:59

Do you have a small testcase you can post or send (m.janssens, domain opencfd.co.uk)

bruce November 6, 2009 07:50

small test case
 
1 Attachment(s)
Hello Mattijs

I have posted a very small test case that reflects my problem.

Thanks and Best Regards

mattijs November 11, 2009 13:35

1 Attachment(s)
I've pushed some changes to 1.6.x. If I use attached controlDict (moved the minX? patch to coincide with the inlet.stl) it works. Please let me know if there any more problems with it.

Thanks,

Mattijs

bruce November 12, 2009 09:23

It works
 
Hello Mattijs

Thanks, it works


All times are GMT -4. The time now is 08:13.