[renyun0511]

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"
p.storePrevIter();
// Pressure-velocity SIMPLE corrector
{
// Momentum predictor
tmp<fvVectorMatrix> UEqn
(
fvm::div(phi, U)
+ turbulence->divDevReff(U)
);
mrfZones.addCoriolis(UEqn());
UEqn().relax();
solve(UEqn() == -fvc::grad(p));
p.boundaryField().updateCoeffs();
volScalarField rAU = 1.0/UEqn().A();
U = rAU*UEqn().H();
UEqn.clear();
phi = fvc::interpolate(U, "interpolate(HbyA)") & mesh.Sf();
mrfZones.relativeFlux(phi);
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)
{
return;
}
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

jennifer

[waynezw0618]

Just remove ：
mrfZones.addCoriolis(UEqn());

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

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

Thanks,
Sivakumar