[renyun0511]

Hi, all,

Recently, I try to simulate the transient inner flow of centrifugal pump by PimpleDyMFoam on OpenFOAM-2.1.0. But it always ended with the floating point exception error after some steps. I looked into the internal Field of my case by Ensight software and I found that there was a wrong pressure distribution in this centrifugal pump, although the directions of both velocity at inlet and velocity at outlet are right and the direction of impeller rotation is also right. I don't know why it was this? Can someone give me some hints? I will be very appreciated it!

In order to explain my questions better, I give the pressure distributions of this pump with steadystate and transient simulation at the first output time step, respectively. Here are:

Attachment 21641,Attachment 21642

The boundary file:
boundary.txt
The fvSchemes file:
fvSchemes.txt
The fvSolution file:
fvSolution.txt
The controlDict file:
controlDict.txt

The initial BCs:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "3500";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [0 1 -1 0 0 0 0];
internalField   nonuniform List<vector> 
1196068
(
(0.716798 0.840057 -7.4744)
(0.638485 0.861783 -7.48279)
(0.681562 0.835391 -7.58056)
...  ...
    wall
    {
        type            fixedValue;
        value           nonuniform List<vector>
88020
(
(0 0 0)
(0 0 0)
(0 0 0)
...  ...
    wallimp4
    {
        type            movingWallVelocity;
        value           uniform (0 0 0);
    }
    wallimp3
    {
        type            movingWallVelocity;
        value           uniform (0 0 0);
    }
    wallimp2
    {
        type            movingWallVelocity;
        value           uniform (0 0 0);
    }
    wallimp1
    {
        type            movingWallVelocity;
        value           uniform (0 0 0);
    }
    outlet
    {
        type            zeroGradient;
    }
    inlet
    {
        type            fixedValue;
        value           uniform (0 -2.708 0);
    }
}

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "3500";
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [0 2 -2 0 0 0 0];
internalField   nonuniform List<scalar>
1196068
(
0.38443
0.248261
0.21847
...  ...
    wall
    {
        type            zeroGradient;
    }
    wallhgb
    {
        type            zeroGradient;
    }
    wallqgb
    {
        type            zeroGradient;
    }
    interafcevol1
    {
        type            cyclicAMI;
        value           nonuniform List<scalar>
1446
(
5.50122
5.60699
5.11278
...  ...
    interfacevol
    {
        type            cyclicAMI;
        value           nonuniform List<scalar>
22196
(
-92.0152
-109.644
-123.262
...  ...
    interfacesuc
    {
        type            cyclicAMI;
        value           nonuniform List<scalar>
2368
(
-675.839
-676.411
-675.763
...  ...
    interfaceimp1
    {
        type            cyclicAMI;
        value           nonuniform List<scalar>
1120
(
-676.74
-675.488
-676.987
...  ...
    interfaceimp
    {
        type            cyclicAMI;
        value           nonuniform List<scalar>
8680
(
-105.911
-104.788
-102.254
...  ...
    interfaceckyc
    {
        type            cyclicAMI;
        value           nonuniform List<scalar>
748
(
-6.02945
-7.33232
-6.99204
...  ...
    outlet
    {
        type            fixedValue;
        value           uniform 0;
    }
    inlet
    {
        type            zeroGradient;
    }
}

Thank you very much!


Yun Ren