[DOliveira]

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?

[WernerW]

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

[HenningW]

I know this is a old thread, but nonetheless:
-Your U boundaries, as far as I can see, only have an inflow, but no face that allows the flow to exit.
-The way you approach this looks like a deforming mesh to me (but knowing your boundary and dynamicMeshDict would help. This is something I recommend not to do. Your deformation will be too much. Use a slidingMesh, if possible, instead. It allows rotation and does not require any changes in your cellMotionU