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

decomposePar looses patchField entries

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

Like Tree1Likes
  • 1 Post By letzel

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 15, 2013, 19:02
Default decomposePar looses patchField entries
  #1
Member
 
Marcus Letzel
Join Date: Sep 2012
Location: Bremen
Posts: 35
Rep Power: 13
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: Bremen
Posts: 35
Rep Power: 13
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 Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Help with if statement CHARLES OpenFOAM Programming & Development 17 August 22, 2021 04:14
Elementwise multiplication operator johndeas OpenFOAM Running, Solving & CFD 3 March 9, 2019 14:03
adding compressible option to ptot immortality OpenFOAM Programming & Development 13 June 15, 2015 16:00
decomposePar gives errors of_user_ OpenFOAM 1 July 4, 2011 06:27
A stupid question luckyluke OpenFOAM Running, Solving & CFD 14 August 13, 2007 05:25


All times are GMT -4. The time now is 22:50.