Error in using createBaffles
 

Dear Foam Fathers,

Good afternonn!

OpenFOAM-1.5
OpenSUSE-11.0
gcc-4.3.1

Procedure to reproduce the error messages is:
Click to download the case: Attachment 140
1. blockMesh
2. faceSet
3. Add another entry in file boundary, name this new patch "plate", as followed,
4. createBaffles plate1 plate
Then:

Thanks for any help.
Regards,

Daniel

Why do you have fixedValue bc on a cyclic?
The number of the startface is not 321432.

so apart from a few mistakes in the supplied case it works for me using 1.5.x

So I would say if it is a bug and not a mistake on your side, this has been fixed already.

Sorry for that I have uploaded a wrong case.

Actually, I'm sure I did not make any mistakes there (I mean fixedValues for cyclic bc, and startSurface number as you mentioned).

So I am compiling version-1.5.x, and see tomorrow if it works or not.

Thank you Niklas!

Shame! I made a mistake again, I spent too much time sitting before computer.
It works now. Not a bug, sorry!
And thank you, Niklas!

Dear developers,

Hi, I want to report to you again, though not knowing where the bug is, but in my case now, I could run createBaffles using OpenFOAM-1.5-dev, but using OpenFOAM-1.6 ended with a segv. So I guess this is a bug in OpenFOAM-1.6.

So now I will turn on the OpenFOAM-1.5-dev when I need to createBaffles, and then turn back 1.6 for my pisoFoam case, and this works correctly.

I don't know whether you are interested in this.:confused:

Did you try 1.6.x? If so can you post some more information about the failure?

Thanks, Mattijs

My pleasure, I will build 1.6.x when my simulation is finished.

BTW, I really love the idea in version 1.6, that the turbulence modeling structure is rearranged, this is perfect in my opinion, for RANS and LES are going nearer and nearer recent years.

Best regards,

FOAM FATAL ERROR:
Cannot find faceZone f0
Valid zones are
0
(
)

the 1.6.x version takes a faceZone instead of a faceSet as in 1.5. For its usage see e.g. Allrun script in the lagrangian/porousExplicitSourceReactingParcelFoam/filter tutorial.

It works, thanks a lot.

One more question, is there a way to view the number of each face in paraview, I found faceSet or cellSet utilities are difficult to use.
In my case, in order to select some faces (with different angles), I need to use an inclined box first, and then test each number manually, which is very inconvenient.

Hi, mattijs!

Could you help and have a look at this, please.
Attachment 1416
It still goes wrong. The process to reproduce the error is
1. Download the A6 case;
2. blockMesh;
3. use the boundary_AfterPlateAdded to replace the boundary file;
4. (No cellSet and faceSet and setsToZones are needed, since I have already done it, the files are already exist);
5. createBaffles f0 plate -overwrite;

Thank you.

The boundary_AfterPlateAdded has an extra patch 'plate' but it has nFaces != 0 so is not a drop-in replacement. Run checkMesh to check.

(why not add plate to your blockMeshDict)

2) you'll probably want to unset FOAM_SIGFPE because the patchFields created on patch plate have no value (since they cannot be mapped from existing patch faces). If they happen to have an illegal value you will get a sigfpe.

Thank you, mattijs.

The case is an external flow past a bridge deck (twin box), but under each box (near the center of the whole deck) there is a guide vane. The guide vanes are too small to use blockMesh to handle in a convenient way, so I decided to treat it as a zero thickness plate. So first, I disregard the guide vanes, and build the whole mesh with blockMesh (quite complicated), then I want to use createBaffles to create the guide vanes with zero thicknesses.

That's why I did not add plate to my original blockMeshDict.

In an old post, I remember the way to createBaffles is to modify the boundary file after faceSet is done, (1) modify the patches number, (2) add another patch entry but leave the nFaces empty meanwhile modify the startFace.

You know, I managed to do it successfully in a 2D mesh which is exactly the same as in this case A6 except the spanwise grid number is unity, but in this 3D case I met the segv, and not knowing what was going wrong here.

Howto? Just comment out "unset FOAM_SIGFPE=" in etc/bashrc?

I have tried many possibilities, but to no avail. I don't know why the utility works in 2D case, but fails in 3D.

So sorry, I did not quite understand your answer, I hope my description are clear and enough.

Many thanks.

It works on your case.

You can add 'plate' to your blockMeshDict as any other patch but just leave the list of faces empty. Then when you run blockMesh it will create a patch 'plate' with zero faces in it.

FOAM_SIGFPE is a environment variable. See your shell on how to unset it.

O, great, it works now, thank you very much.

Dear lakeat,

I am using version-1.6.
First I generate fluent mesh which has an interior face int_WALL.
Then I use fluent3DMeshToFoam,

in faceZones file has this int_WALL.

so I use
createBaffles int_WALL WALL command.
but it display:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Converting faces on zone int_WALL1 into baffles.

Found 1121 faces on zone int_WALL1



Cannot find patch WALL1
Valid patches are
6
(
INLET
OUTLET
LEFT
RIGHT
TOP
GROUND
)


From function findPatchID(const polyMesh&, const word&)
in file createBaffles.C at line 115.

FOAM exiting
.

Then I add
WALL
{
type wall;
nFaces 0;
startFace 127074;
}
to boundary file,
WALL
{
type wall;
nFaces 0;
startFace 127074;
}

but still have error:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0



Expected a ')' or a '}' while reading PtrList, found on line 56 the word 'WALL1'

file: /home/panda60/OpenFOAM/run/exp/constant/polyMesh/boundary at line 56.

From function Istream::readEndList(const char*)
in file db/IOstreams/IOstreams/Istream.C at line 152.

FOAM exiting

.
How could I solve this problem ?
Thank you.

You'll need to adapt the number at the top of the boundary file. It probably still says 6 but you've got 7 entries in there now.

Dear ,Mattijs

Thank you very much.
Now already OK.

but I have nearly 30 internal walls, it is indeed a large work.
for 1,000,000 mesh, I do this again and again and wait for 2 hours.:)

If we have command to convert these 30 faceZones to internal walls for only one time ?

No, there is no command to do multiple faceZones.

Dear Mattijs,
If I want to create an internal patch for sampling, what boundary condition should I set for this patch ?

this patch belongs to internal domain, is not boundary. how should I set for it, then it can be an internal patch, only used for sampling, but doesn't influence the flow ?

If I should use createPatch ? then boundary condition ?

Thank you very much.