CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   coriolis force in MRFSimpleFoam (http://www.cfd-online.com/Forums/openfoam-solving/71279-coriolis-force-mrfsimplefoam.html)

renyun0511 December 23, 2009 01: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"
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 December 23, 2009 02:40


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

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

Thanks,
Sivakumar


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