CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Compute Zalesak's disk

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

Like Tree1Likes
  • 1 Post By jianfeng

LinkBack Thread Tools Display Modes
Old   May 21, 2013, 08:23
Default Compute Zalesak's disk
New Member
jianfeng zou
Join Date: Feb 2013
Posts: 4
Rep Power: 6
jianfeng is on a distinguished road
Sponsored Links
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())

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

field U;
expression "vector(0,0,6.28)^(pos()-vector(0.0,0.0,0.0))";
field p_rgh;
expression "0";
field alpha1;
expression "1";
condition "pow(pos().x-0.0,2) + pow(pos().y-0.25,2) < pow(0.15,2)";
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
File Type: jpg visit0000.jpg (24.4 KB, 60 views)
dduque likes this.
jianfeng is offline   Reply With Quote
Sponsored Links

Old   May 23, 2013, 10:19
Join Date: Nov 2012
Posts: 71
Rep Power: 6
Henning86 is on a distinguished road
what is your courant number?
Henning86 is offline   Reply With Quote


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
actuator disk Atit CFX 9 July 26, 2016 09:39
ActuatorDiskExplicitForce in OF2.1. Help be_inspired OpenFOAM Programming & Development 8 July 3, 2014 10:54
vof rotating disk 2d bene FLUENT 1 October 24, 2010 13:18
Actuator Disk Model sujialin CFX 0 April 18, 2010 16:29

Sponsored Links

All times are GMT -4. The time now is 00:43.