Can not get correct result in inducer simulation by pimpleDyMFoam
 

Hello, everyone

Recently, I am confusing with a problem. Last week, I have spent much time but still could not solve it.
I want to simulate a running of inducer. I have tested it at 6000 rpm and 600 rpm. Even at 6000 rpm, the edge of inducer is just running at 20m/s. But I got a very strange result which above 1000m/s in 0~0.0006s. Pressure is also very strange. I calculated it by head coefficient. The data I got is 100 times than my calculating result. Now, I will show you all the file I tested. I am using the file in pimpleDyMFoam propeller.


FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField
{


INLET
{
type zeroGradient;
}

OUTLET
{
type fixedValue;
value uniform 0;
}

INTERFACE_11
{
type cyclicAMI;
}

INTERFACE_12
{
type cyclicAMI;
}

INTERFACE_21
{
type cyclicAMI;
}

INTERFACE_22
{
type cyclicAMI;
}

WALL_OUTER_ROTATESIDE
{
type zeroGradient;
}

WALL_OUTER_INLETSIDE
{
type zeroGradient;
}

WALL_OUTER_OUTLETSIDE
{
type zeroGradient;
}

WALL_ROTATE_CORE_FIX
{
type zeroGradient;
}

WALL_ROTATE_CORE
{
type zeroGradient;
}
}



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

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

internalField uniform (0 0 0);

boundaryField
{


WALL_ROTATE_CORE
{
type movingWallVelocity;
value uniform (0 0 0);
}

INLET
{
type pressureInletVelocity;
value uniform (0 0 0);
}

OUTLET
{
type inletOutlet;
inletValue uniform (0 0 0);
}

WALL_ROTATE_CORE_FIX
{
type noSlip;
}

WALL_OUTER_OUTLETSIDE
{
type noSlip;
}

WALL_OUTER_INLETSIDE
{
type noSlip;
}

WALL_OUTER_ROTATESIDE
{
type noSlip;
}

INTERFACE_11
{
type cyclicAMI;
}

INTERFACE_12
{
type cyclicAMI;
}

INTERFACE_21
{
type cyclicAMI;
}

INTERFACE_22
{
type cyclicAMI;
}

}



FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application pimpleDyMFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 1;

deltaT 0.00001;
writeControl adjustableRunTime;
writeInterval 0.0001;

////- For testing with moveDynamicMesh

///deltaT 1e-5;
//writeControl timeStep;
//writeInterval 1;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression on;

timeFormat general;

timePrecision 6;

runTimeModifiable true;

adjustTimeStep yes;

maxCo 1;

functions
{

}


FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh solidBodyMotionFvMesh;

motionSolverLibs ( "libfvMotionSolvers.so" );

solidBodyMotionFvMeshCoeffs
{
cellZone FLUID_ROTATE_CORE_TRI;

solidBodyMotionFunction rotatingMotion;
rotatingMotionCoeffs
{
origin (0 0 0);
axis (0 1 0);
omega 628; // rad/s
}
}


FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
pcorr
{
solver GAMG;
tolerance 1e-2;
relTol 0;
smoother DICGaussSeidel;
cacheAgglomeration no;
maxIter 50;
}

p
{
$pcorr;
tolerance 1e-5;
relTol 0.01;
}

pFinal
{
$p;
tolerance 1e-6;
relTol 0;
}

"(U|k|epsilon)"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-6;
relTol 0.1;
}

"(U|k|epsilon)Final"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-6;
relTol 0;
}
}

PIMPLE
{
correctPhi no;
nOuterCorrectors 3;
nCorrectors 2;
nNonOrthogonalCorrectors 1;
}

relaxationFactors
{
"(U|k|epsilon).*" 1;
}

cache
{
grad(U);
}

In order to decrease the simulation time, I was using laminar model to test it. Could you anyone tell what is wrong with my case? Thank you very much.

ZX

I show you pics to let you know my uncorrect answer easily