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

decomposePar looses patchField entries

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 15, 2013, 19:02
Default decomposePar looses patchField entries
  #1
Member
 
Marcus Letzel
Join Date: Sep 2012
Location: Aurich
Posts: 35
Rep Power: 4
letzel is on a distinguished road
Dear Foamers,

my following workflow does not work, I hope it just requires a simple correction:

Code:
blockMesh > blockMesh.log 2>&1
decomposePar -force > decomposePar1.log 2>&1
surfaceFeatureExtract > surfaceFeatureExtract.log 2>&1
foamJob -parallel snappyHexMesh
foamJob -parallel checkMesh -latestTime
cp -pv 0/* 0.2 > cp.log 2>&1
setFields > setFields.log 2>&1
decomposePar -fields -latestTime -constant > decomposePar2.log 2>&1
foamJob -parallel pimpleFoam
snappy performs two steps, dt is 0.1 s, hence I copy the initial fields from 0 to latestTime=0.2.

pimpleFoam complains:
Code:
Create mesh for time = 0.2

Reading field p

[0] 
[0] 
[0] --> FOAM FATAL IO ERROR: 
[0] Cannot find patchField entry for alexanderplatz_alexanderplatz
[0] 
[0] file: /home/meineruser/OpenFOAM/meineruser-2.2.2/run/[...]/processor0/0.2/p.boundaryField from line 26 to line 41.
[0] 
[0]     From function GeometricField<Type, PatchField, GeoMesh>::GeometricBoundaryField::readField(const DimensionedField<Type, GeoMesh>&, const dictionary&)
[0]     in file /opt/openfoam222/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 206.
[0] 
FOAM parallel run exiting
p in processor0/0.2 simply reads
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.2                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0.2";
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    defaultName
    {
        type            zeroGradient;
    }
    procBoundary0to1
    {
        type            processor;
        value           uniform 0;
    }
    procBoundary0to2
    {
        type            processor;
        value           uniform 0;
    }
    procBoundary0to3
    {
        type            processor;
        value           uniform 0;
    }
}


// ************************************************************************* //
But p in the 0.2 subdirectory of the case directory still has all the patches specified:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

libs (
         "libgroovyBC.so"
     ); 

boundaryField
{
    alexanderplatz_patch_ausgang_db
    {
        type            groovyBC;
        value           uniform 0;
        valueExpression "0";
        gradientExpression "0";
        fractionExpression "inVel<0 ? 1 : 0";
        timelines (
          {
            name        inVel;
            outOfBounds clamp;
            fileName    "$FOAM_CASE/input_data/patch_ausgang_db.txt";
          }
                  );
    }

    "(alexanderplatz_patch_ausgang_a.*|alexanderplatz_patch_ausgang_b.*|alexanderplatz_patch_ausgang_l.*|alexanderplatz_patch_ausgang_u.*|alexanderplatz_patch_t.*)"
    {
        type            fixedValue;
    value        uniform 0;
    }

    "(alexanderplatz_alexanderplatz|treppen_treppen|saeulen_saeulen|defaultName)"
    {
        type            zeroGradient;
    }
}

// ************************************************************************* //
Looking forward to hear any suggestions to revise my workflow.

Cheers,
Marcus
letzel is offline   Reply With Quote

Old   November 18, 2013, 05:06
Default
  #2
Member
 
Marcus Letzel
Join Date: Sep 2012
Location: Aurich
Posts: 35
Rep Power: 4
letzel is on a distinguished road
Forgot to add that a slightly modified workflow including a reconstructParMesh step works as follows:
Code:
blockMesh > blockMesh.log 2>&1
decomposePar -force > decomposePar1.log 2>&1
surfaceFeatureExtract > surfaceFeatureExtract.log 2>&1
foamJob -parallel snappyHexMesh
foamJob -parallel checkMesh -latestTime
reconstructParMesh -latestTime -mergeTol 1e-06
cp -pv 0/* 0.2 > cp.log 2>&1
setFields > setFields.log 2>&1
decomposePar -latestTime -force > decomposePar2.log 2>&1
foamJob -parallel pimpleFoam
Is it possible to avoid this extra step?
letzel 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
adding compressible option to ptot immortality OpenFOAM Programming & Development 13 June 15, 2015 15:00
Help with if statement CHARLES OpenFOAM Programming & Development 15 September 18, 2013 18:29
Elementwise multiplication operator johndeas OpenFOAM Running, Solving & CFD 2 May 14, 2012 09:45
decomposePar gives errors of_user_ OpenFOAM 1 July 4, 2011 05:27
A stupid question luckyluke OpenFOAM Running, Solving & CFD 14 August 13, 2007 04:25


All times are GMT -4. The time now is 00:45.