CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   coriolis force in MRFSimpleFoam (

renyun0511 December 23, 2009 02:49

coriolis force in MRFSimpleFoam
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


waynezw0618 December 23, 2009 03:40

Just remove :

but i don`t think it is a topic for study

sivakumar October 28, 2013 08:27

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....


All times are GMT -4. The time now is 21:07.