|
[Sponsors] |
To modify the interFoam solver for cavitating |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 17 ![]() |
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} |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 17 ![]() |
![]() |
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 17 ![]() |
Could anybody give me some hints?
|
|
![]() |
![]() |
![]() |
![]() |
#4 |
Member
Sung-Eun Kim
Join Date: Mar 2009
Posts: 76
Rep Power: 17 ![]() |
For stable solution of Gamma, you will need Sp as well in case the source term becomes negative.
|
|
![]() |
![]() |
![]() |
![]() |
#5 |
New Member
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 17 ![]() |
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. |
|
![]() |
![]() |
![]() |
![]() |
#6 |
New Member
morteza
Join Date: Mar 2009
Posts: 18
Rep Power: 17 ![]() |
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! |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
About interFoam solver | zou_mo | OpenFOAM Running, Solving & CFD | 129 | December 2, 2019 05:39 |
How to modify interFoam to keep a constant velocity field | thibault_pringuey | OpenFOAM | 0 | January 9, 2009 05: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 |