[jianfeng]

Hi, everyone. I want to test the performance of the solver of interFoam by computing the solid rotation of Zalesak's disk. I set up a simulation as following:

1) First, I comment the Pressure-velocity PIMPLE corrector loop in the main source file 'interFoam.C' lived in the directory of interFoam solver.

// --- Pressure-velocity PIMPLE corrector loop
/*
while (pimple.loop())
{
#include "UEqn.H"

// --- Pressure corrector loop
while (pimple.correct())
{
#include "pEqn.H"
}

if (pimple.turbCorr())
{
turbulence->correct();
}
}
*/

2) Generate a intial field with funkySetField. The funkySetFieldsDict is as following:

expressions
(
circleVel
{
field U;
expression "vector(0,0,6.28)^(pos()-vector(0.0,0.0,0.0))";
}
circleP
{
field p_rgh;
expression "0";
}
diskGamma0
{
field alpha1;
expression "1";
condition "pow(pos().x-0.0,2) + pow(pos().y-0.25,2) < pow(0.15,2)";
}
diskGamma1
{
field alpha1;
expression "0";
condition "pos().x > -0.025 && pos().x < 0.025 && pos().y > 0.15 && pos().y < 0.5";
}
);

3) But the Zalesak's disk loses its sharp shape very soon, the following is the result at t=0.25 s.

[Henning86]

what is your courant number?

[ndtrong]

Originally Posted by jianfeng

Hi, everyone. I want to test the performance of the solver of interFoam by computing the solid rotation of Zalesak's disk. I set up a simulation as following:

1) First, I comment the Pressure-velocity PIMPLE corrector loop in the main source file 'interFoam.C' lived in the directory of interFoam solver.

// --- Pressure-velocity PIMPLE corrector loop
/*
while (pimple.loop())
{
#include "UEqn.H"

// --- Pressure corrector loop
while (pimple.correct())
{
#include "pEqn.H"
}

if (pimple.turbCorr())
{
turbulence->correct();
}
}
*/

2) Generate a intial field with funkySetField. The funkySetFieldsDict is as following:

expressions
(
circleVel
{
field U;
expression "vector(0,0,6.28)^(pos()-vector(0.0,0.0,0.0))";
}
circleP
{
field p_rgh;
expression "0";
}
diskGamma0
{
field alpha1;
expression "1";
condition "pow(pos().x-0.0,2) + pow(pos().y-0.25,2) < pow(0.15,2)";
}
diskGamma1
{
field alpha1;
expression "0";
condition "pos().x > -0.025 && pos().x < 0.025 && pos().y > 0.15 && pos().y < 0.5";
}
);

3) But the Zalesak's disk loses its sharp shape very soon, the following is the result at t=0.25 s.

Did you solve your problems, can you share with me your final setup for this simulations?