CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

To modify the interFoam solver for cavitating

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 10, 2007, 09:02
Default The equation I have used is bo
  #1
New Member
 
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 8
zjucfd is on a distinguished road
The equation I have used is borrowed from the paper by Wei Shyy.

Jiongyang Wu, Yogen Utturkar, Wei Shyy. ASSESSMENT OF MODELING STRATEGIES FOR CAVITATING FLOW AROUND A HYDROFOIL, CAV2003.

/attach{eq.jpg}
zjucfd is offline   Reply With Quote

Old   July 10, 2007, 09:05
Default http://www.cfd-online.com/Open
  #2
New Member
 
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 8
zjucfd is on a distinguished road

zjucfd is offline   Reply With Quote

Old   July 11, 2007, 12:26
Default Could anybody give me some hin
  #3
New Member
 
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 8
zjucfd is on a distinguished road
Could anybody give me some hints?
zjucfd is offline   Reply With Quote

Old   July 18, 2007, 09:55
Default For stable solution of Gamma,
  #4
sek
Member
 
Sung-Eun Kim
Join Date: Mar 2009
Posts: 76
Rep Power: 8
sek is on a distinguished road
For stable solution of Gamma, you will need Sp as well in case the source term becomes negative.
sek is offline   Reply With Quote

Old   July 18, 2007, 10:43
Default I have also tried this form:
  #5
New Member
 
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 8
zjucfd is on a distinguished road
I have also tried this form:


volScalarField Sp = Cdest*rho2*min(pvap-pvap,pd-pvap)/(rho1*(rho1*Uref*Uref*0.5)*tref)

volScalarField Su = Cprod*rho2*gamma*gamma*(scalar(1)-gamma)/(rho1*tref);

MULES::explicitSolve01(gamma, phi, phiGamma, Sp, Su);


The results are still bad.
zjucfd is offline   Reply With Quote

Old   August 26, 2007, 08:03
Default Hi, I are attempted to implem
  #6
New Member
 
morteza
Join Date: Mar 2009
Posts: 18
Rep Power: 8
morteza is on a distinguished road
Hi,
I are attempted to implement SINGHAL cavitation model into lesInterFoam in order to model a supercavitation condition over a wedge.
So, I modified the gamma.H according to the SINGHAL model in the solver.
here is the gamma.H code:

for (int gCorr=0; gCorr<nGammaCorr; gCorr++)
{

# include "rhoPhiGamma.H"

fvScalarMatrix gammaEqn
(
fvm::ddt(gamma)
+ fvm::div(phi, gamma)
+ fvm::div
(
-fvc::flux(-phir, scalar(1) - gamma, gammarScheme),
gamma,
gammarScheme
)
==
mdote
- fvm::Sp(mdote+mdotc,gamma)
);

gammaEqn.solve();

rhoPhi = gammaEqn.flux()*(rho1 - rho2) + phi*rho2;
}

-------------------------------------------------------------------------------- -------
where rhoPhiGamma.H is:

volScalarField Kinetic=turbulence->k();

forAll (mesh.C(),i)
{
if ( pd[i] < psat.value() )
{
mdote[i] = Ce.value()*rho1.value()*rho2.value()/0.02*sqrt( float( 2*Kinetic[i]*(psat.value()-pd[i])/(3*rho2.value())));
mdotc[i] = 0.0;
}
else
{
mdotc[i] = Cc.value()*rho2.value()*rho2.value()/0.02*sqrt( float( 2*Kinetic[i]*(pd[i]-psat.value())/(3*rho2.value())));
mdote[i] = 0.0;
}
}

Could anybody give me some help?
thanks!
morteza 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
About interFoam solver zou_mo OpenFOAM Running, Solving & CFD 127 May 25, 2011 16:30
How to modify interFoam to keep a constant velocity field thibault_pringuey OpenFOAM 0 January 9, 2009 06:59
About interfoam solver qiu OpenFOAM Running, Solving & CFD 0 May 6, 2007 22:48
How to modify variables in Fluent solver yizhou FLUENT 1 June 19, 2006 20:00
Need documentation for interFOAM solver mer OpenFOAM Running, Solving & CFD 5 May 31, 2006 12:22


All times are GMT -4. The time now is 16:43.