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

Problem running PimpleDyMFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 26, 2015, 17:07
Default Problem running PimpleDyMFoam
  #1
New Member
 
Join Date: Apr 2015
Posts: 10
Rep Power: 3
DOliveira is on a distinguished road
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?
DOliveira is offline   Reply With Quote

Old   August 18, 2015, 10:44
Default Use Propeller Tutorial
  #2
Member
 
Werner
Join Date: Jul 2015
Location: West Lafayette, USA
Posts: 34
Rep Power: 2
WernerW is on a distinguished road
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
WernerW is offline   Reply With Quote

Reply

Tags
pimpledymfoam

Thread Tools
Display Modes

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem running movingCylinders case in parallel with foam-extend-3.1 mhkenergy OpenFOAM Running, Solving & CFD 4 April 12, 2016 03:43
Problem with parallel run of my solver based on pimpleDyMFoam o.kotsur OpenFOAM Running, Solving & CFD 0 October 6, 2013 03:44
Problem while running in Highperformance computing environment Phanipavan STAR-CD 1 September 11, 2013 06:42
problem with running in parallel dhruv OpenFOAM 3 November 25, 2011 06:06
Kubuntu uses dash breaks All scripts in tutorials platopus OpenFOAM Bugs 8 April 15, 2008 07:52


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