# Compute Zalesak's disk

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

May 21, 2013, 08:23
Compute Zalesak's disk
#1
New Member

jianfeng zou
Join Date: Feb 2013
Posts: 4
Rep Power: 5
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.
Attached Images
 visit0000.jpg (24.4 KB, 55 views)

 May 23, 2013, 10:19 #2 Member   Join Date: Nov 2012 Posts: 67 Rep Power: 5 what is your courant number?

 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 Atit CFX 9 July 26, 2016 09:39 be_inspired OpenFOAM Programming & Development 8 July 3, 2014 10:54 bene FLUENT 1 October 24, 2010 13:18 sujialin CFX 0 April 18, 2010 16:29

All times are GMT -4. The time now is 01:04.

 Contact Us - CFD Online - Top