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

Can not get correct result in inducer simulation by pimpleDyMFoam

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 19, 2017, 01:58
Default Can not get correct result in inducer simulation by pimpleDyMFoam
  #1
New Member
 
TeiGyou
Join Date: Oct 2016
Location: Tokyo,Japan
Posts: 20
Rep Power: 10
zxzx is on a distinguished road
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
zxzx is offline   Reply With Quote

Old   February 19, 2017, 03:43
Default
  #2
New Member
 
TeiGyou
Join Date: Oct 2016
Location: Tokyo,Japan
Posts: 20
Rep Power: 10
zxzx is on a distinguished road
I show you pics to let you know my uncorrect answer easily
Attached Images
File Type: png Screenshot from 2017-02-20 18-41-40.png (168.1 KB, 12 views)
File Type: png Screenshot from 2017-02-20 18-40-59.png (158.6 KB, 12 views)
zxzx is offline   Reply With Quote

Reply


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
PimpleDymFOAM simulation problem of Savonius rotor Sunzida OpenFOAM Running, Solving & CFD 0 March 14, 2016 02:45
Confusion abt expt / non-dim simulation result comparison quarkz Main CFD Forum 8 November 1, 2015 09:03
Simulation Problem Where Suddenly Finished With No Result Produced PeiSan Fidelity CFD 4 August 21, 2014 08:06
how to make sure the simulation result is correct? sham81 CFX 3 March 22, 2004 17:41
Result is not correct Li FLUENT 1 December 30, 2000 01:21


All times are GMT -4. The time now is 15:34.