CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Case decomposition of a scene with a small object

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 26, 2016, 10:59
Default Case decomposition of a scene with a small object
  #1
New Member
 
Vladislav Galas
Join Date: Jun 2016
Posts: 5
Rep Power: 2
gallon is on a distinguished road
Hello CFD, hope this is the right forum to ask the question.

The goal is to place an stl model and refine the mesh.
Everything goes fine in serial but since I want to refine up to 5-6 levels it takes too long.

So here's my script:
Code:
blockMesh
decomposePar [-force]
At this decomposePar knows nothing about WALL10 (the surface read from .stl by sHM).
OK, then I run
Code:
surfaceFeautureExtract
mpirun -np 2 snappyHexMesh -overwrite -parallel
mpirun -np 2 sonicFoam -parallel
The surface is read and refined in parallel OK. If a given processors's decomposed domain does not contain a part of the surface, it is still mentioned in processorX/constant/polyMesh/boundary as an empty boundary:
Code:
[v.galas@node1 tmp]$ cat processor0/constant/polyMesh/boundary
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0.1                                |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

4
(
    OUTL2
    {
        type            patch;
        nFaces          22816;
        startFace       418376;
    }
    INLE1
    {
        type            patch;
        nFaces          0;
        startFace       441192;
    }
    WALL10
    {
        type            wall;
        inGroups        1(wall);
        nFaces          0;
        startFace       441192;
    }
    procBoundary0to1
    {
        type            processor;
        inGroups        1(processor);
        nFaces          1360;
        startFace       441192;
        matchTolerance  0.0001;
        transform       unknown;
        myProcNo        0;
        neighbProcNo    1;
    }
)

// ************************************************************************* //
Code:
mpirun -np 2 sonicFoam -parallel
But now the solver run at the processors without WALL10 cannot find the BC for it.
I have tried mapFields and redistributePar but to no avail.

p.s. sHM (serial) -> decomposePar -> sonicFoam (parallel) works ok.
Any ideas?
gallon 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
Odd lagrangian particle behavior at small size phandy OpenFOAM Running, Solving & CFD 1 November 5, 2015 16:54
Simple channel case using cyclicAMI will not converge cbcoutinho OpenFOAM Running, Solving & CFD 3 August 4, 2015 12:28
Changing the grid on the same set-up Katya FLUENT 12 October 8, 2009 16:31
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 12:24


All times are GMT -4. The time now is 10:11.