|
[Sponsors] |
Can not get correct result in inducer simulation by pimpleDyMFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 19, 2017, 01:58 |
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 |
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 |
|
February 19, 2017, 03:43 |
|
#2 |
New Member
TeiGyou
Join Date: Oct 2016
Location: Tokyo,Japan
Posts: 20
Rep Power: 10 |
I show you pics to let you know my uncorrect answer easily
|
|
|
|
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 |