CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

coriolis force in MRFSimpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 23, 2009, 02:49
Smile coriolis force in MRFSimpleFoam
  #1
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 7
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"
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
renyun0511 is offline   Reply With Quote

Old   December 23, 2009, 03:40
Default
  #2
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 307
Rep Power: 8
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to waynezw0618

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

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

Old   October 28, 2013, 09:27
Default
  #3
Senior Member
 
sivakumar selvaraju
Join Date: Mar 2009
Location: Cape Town - South Africa
Posts: 186
Rep Power: 7
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....

Thanks,
Sivakumar
sivakumar is offline   Reply With Quote

Reply

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 08:07
DEFINE_CG_MOTION and pressure force Teo Fumagalli FLUENT 0 April 11, 2008 11: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 05:23
pressure force and viscous force CFD FLUENT 0 March 7, 2006 02:03
Coriolis force Saturn FLUENT 5 September 30, 2004 16:56


All times are GMT -4. The time now is 18:13.