# Problem running PimpleDyMFoam

 Register Blogs Members List Search Today's Posts Mark Forums Read

 July 26, 2015, 17:07 Problem running PimpleDyMFoam #1 New Member   Join Date: Apr 2015 Posts: 10 Rep Power: 3 Hi everyone! I'm trying to model a running turbine in OpenFoam. I'm using a library that I compiled based on the angularOscillatingVelocity. I simply removed the oscillating part and changed the equation scalar angle = angle0 + amplitude*sin(omega*t.value()) to scalar angle =angle= + omega*t.value() My cellMotionU dictionary is: boundaryField { default { type empty; } inlet { type fixedValue; value uniform (0 0 0); } outlet { type fixedValue; value uniform (0 0 0); } down { type fixedValue; value uniform (0 0 0); } up { type fixedValue; value uniform (0 0 0); } front { type fixedValue; value uniform (0 0 0); } back { type fixedValue; value uniform (0 0 0); } turbine_patch28539 { type rotatingWallVelocity; axis (1 0 0); origin (0 0 0); angle0 0; omega 20.94 value uniform (1 0 0); } and my U dictionary is: boundaryField { inlet { type fixedValue; value uniform (1 0 0); } outlet { type fixedValue; value uniform (0 0 0); } down { type fixedValue; value uniform (0 0 0); } up { type fixedValue; value uniform (0 0 0); } front { type fixedValue; value uniform (0 0 0); } back { type fixedValue; value uniform (0 0 0); } turbine_patch28539 { type fixedValue; value uniform (0 0 0); My fvSolution is: solvers { pcorr { solver GAMG; tolerance 1e-06; relTol 0.01; smoother DICGaussSeidel; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; maxIter 50; } p { \$pcorr; tolerance 1e-6; relTol 0.01; } pFinal { \$p; tolerance 1e-6; relTol 0; } cellMotionU { solver GAMG; smoother DICGaussSeidel; tolerance 1e-07; reTol 0; cacheAgglomeration no; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; maxIter 50; } U { solver GAMG; smoother DILUGaussSeidel; agglomerator faceAreaPair; nCellsInCoarsestLevel 10; cacheAgglomeration true; tolerance 1e-6; relTol 0.1; } UFinal { \$U tolerance 1e-6; relTol 0; } /* "(U|k|omega)" { solver PBiCG; preconditioner DILU; tolerance 1e-06; relTol 0.1; } "(U|k|omega)Final" { \$U; tolerance 1e-06; relTol 0; } */ k { solver PBiCG; preconditioner DILU; relTol 0.1; } } PIMPLE { correctPhi yes; nOuterCorrectors 2; nCorrectors 3; nNonOrthogonalCorrectors 0; pRefCell 0; pRefValue 0; } relaxationFactors { fields { } equations { "U.*" 1; } } When I run moveMesh everything is ok, but when I run mpirun -np 8 pimpleDyMFoam -parallel I have a problem when it gets to the part that needs to solve the U equation Solving for Ux, Initial residual = -nan, Final residual =nan, No Iterations 1001 I've tryed to changed the solvers but that always happens... can anyone help me?

 August 18, 2015, 10:44 Use Propeller Tutorial #2 Member   Werner Join Date: Jul 2015 Location: West Lafayette, USA Posts: 34 Rep Power: 2 Hi ! , I'm sorry to tell you I can't help you very much with my shallow knowledge of OpenFoam, since I haven't experienced such problem before, But I have a couple questions that may guide you. I'm simulating an axial turbine as well and I just imitated the propeller case of the tutorials. According to that case you must define a Cylinder that contains your turbine and with which the sliding mesh rotates. Do you have such Cylinder somewhere ? If not, since I don't see it in the boundaries, how are you defining the sliding mesh? Regards, Werner

 Tags pimpledymfoam

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post mhkenergy OpenFOAM Running, Solving & CFD 4 April 12, 2016 03:43 o.kotsur OpenFOAM Running, Solving & CFD 0 October 6, 2013 03:44 Phanipavan STAR-CD 1 September 11, 2013 06:42 dhruv OpenFOAM 3 November 25, 2011 06:06 platopus OpenFOAM Bugs 8 April 15, 2008 07:52

All times are GMT -4. The time now is 03:45.