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

coriolis force in MRFSimpleFoam

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

LinkBack Thread Tools Display Modes
Old   December 23, 2009, 02:49
Smile coriolis force in MRFSimpleFoam
Join Date: Jun 2009
Posts: 75
Rep Power: 10
renyun0511 is on a distinguished road
hi all,
i want to study the effect of coriolis force in pump by MRFSimpleFoam,when i gedit MRFSimpleFoam.C,it occurred:

for (runTime++; !runTime.end(); runTime++)
Info<< "Time = " << runTime.timeName() << nl << endl;
# include "readSIMPLEControls.H"
// Pressure-velocity SIMPLE corrector
// Momentum predictor
tmp<fvVectorMatrix> UEqn
fvm::div(phi, U)
+ turbulence->divDevReff(U)
solve(UEqn() == -fvc::grad(p));
volScalarField rAU = 1.0/UEqn().A();
U = rAU*UEqn().H();
phi = fvc::interpolate(U, "interpolate(HbyA)") & mesh.Sf();
adjustPhi(phi, U, p);

and here is my problem:
if i want to delete the coriolis force from the MRFSimpleFoam,what can i do?only delete the red line above?right? or did i need any other operations?for example,the coriolis force described in MRFZones.C also:

void Foam::MRFZone::addCoriolis(fvVectorMatrix& UEqn) const
if (cellZoneID_ == -1)
const labelList& cells = mesh_.cellZones()[cellZoneID_];
const scalarField& V = mesh_.V();
vectorField& Usource = UEqn.source();
const vectorField& U = UEqn.psi();
const vector& Omega = Omega_.value();
forAll(cells, i)
Usource[cells[i]] -= V[cells[i]]*(Omega ^ U[cells[i]]);

should i delete the code above ?
thanks in advance

renyun0511 is offline   Reply With Quote

Old   December 23, 2009, 03:40
Senior Member
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 308
Rep Power: 11
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to waynezw0618

Just remove :

but i don`t think it is a topic for study
waynezw0618 is offline   Reply With Quote

Old   October 28, 2013, 08:27
Senior Member
sivakumar selvaraju
Join Date: Mar 2009
Location: Cape Town - South Africa
Posts: 186
Rep Power: 10
sivakumar is on a distinguished road
Send a message via Skype™ to sivakumar
Dear Renyun0511 and Wayne,
I am working on axial flow fan, I want to see the effect of Coriolis force in the fan.

I am using MRFSimpleFoam, how can I switch of Coriolis force effect.

Do I need to write a new solver or can I command out some of the line in the existing code?

Please help me in this regard....

sivakumar 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
Urgently Need the code of Lift force and VM force Kai Yan Main CFD Forum 0 July 16, 2008 07:07
DEFINE_CG_MOTION and pressure force Teo Fumagalli FLUENT 0 April 11, 2008 10:25
MRFZonesC questions what is the mesh_V and why only Coriolis force no centrifugal force waynezw0618 OpenFOAM Running, Solving & CFD 49 April 8, 2008 04:23
pressure force and viscous force CFD FLUENT 0 March 7, 2006 02:03
Coriolis force Saturn FLUENT 5 September 30, 2004 15:56

All times are GMT -4. The time now is 05:53.