[ishan_ae]

Hello.


When I run decomposePar on my simpleFoam case it decomposes only the internalField data and not the boundaryField data.

I want to run a simpleFoam case with U internalFields values generated from potentialFoam.

The geometry/mesh/patch names I am using is same for both the solvers.

But I had to simplify the BC for potentialFoam run.

Accordingly, I mapped the U field from potentialFoam runto U/0 in simpleFoam folder.

I want to use this mapped field as values to internalFields to simpleFoam run.

Since I am dealing with large mesh, I copied that data into a separate foamFile callled initialU and its format looks like this:

initialU //name of the field
List<vector>
9528301
(
(x1 y1 z1)
(x2 y2 z2)
....
);

This way I can easily update the BC in U/0 without loading that data every time.

By the way, this foamFile is located in 0 folder.

The U/0 folder looks like this:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


 #include "initialU"

dimensions [0 1 -1 0 0 0 0];

internalField nonuniform $initialU;


boundaryField
{
    //my boundary conditions

}

After I run decomposePar on this setup, only the internalField values are decomposed and not the boundaryField data.


The decomposed U file from processor0 looks like this:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   nonuniform List<vector> 
531740
(

(1.92541 -0.998319 0.87524)
(1.81862 -1.05479 0.762567)
(1.78588 -1.01315 0.756257)
...
)
;
    }
}

The last two curly brackets are suspicious to me. Not sure what they are closing.



In the original mapped file, i.e. before I did all these shenanigans, the U/0 looked like this:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   nonuniform List<vector> 

9528301
(
(137.464 -2.17855 -0.202383)
(137.466 -2.17814 -0.20957)
(140.515 -2.36609 -0.160376)
..............
);


boundaryField
{
    //my boundary conditions 

}

Am i complicating things here ? Can someone please suggest what's going wrong ?