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

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

small test case
 

Hello Mattijs

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

Thanks and Best Regards

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

It works
 

Hello Mattijs

Thanks, it works