CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

how to write header file for coalescence model in twoPhaseEulerFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 20, 2013, 14:33
Default how to write header file for coalescence model in twoPhaseEulerFoam
  #1
New Member
 
Join Date: Sep 2010
Location: Windsor, Ontario
Posts: 17
Rep Power: 6
mehrdad_kbg is on a distinguished road
Hi everyone,
I am quite new to openfoam and trying to simulate the two phase flow using twoPhaseEulerFoam. I am adding coalescence model (attached here from an online thesis) to the solver but I do not have the header file for that. I appreciate if you help me to develop the header file.
Thanks,



namespace Foam
{
defineTypeNameAndDebug(HibikiIshiiC, 0);

addToRunTimeSelectionTable
(
coalesenceModel,
HibikiIshiiC,
dictionary
);
}

Foam::HibikiIshiiC::HibikiIshiiC
(
const dictionary& interfaceDict,
const phaseModel& phasea,
const phaseModel& phaseb
)
:
coalescenceModel(interfaceDict, phasea, phaseb)
{}

Foam::HibikiIshiiC::HibikiIshiiC()
{}

Foam::tmp<Foam::volScalarField>Foam::HibikiIshiiC: :C
(
const volScalarField& alpha,
const volScalarField& IAC,
const volScalarField& epsilon,
const volScalarField& Ur

) const

{

scalar GammaC = 0.005;
scalar Kc = 1.29;
scalar alphaMax = 0.52;

volScalarField dBuff = 6.0*alpha/IAC;

dimensionedScalar plusDS
(
"plusDS",
dimensionSet(0, 1, 0, 0, 0, 0, 0),
scalar (1.0e-6)
);
volScalarField sourceRCa = - GammaC*pow(alpha,2.0)*pow(epsilon,1.0/3.0);
volScalarField sourceRCb = pow(dBuff+plusDS,11.0/3.0)*(alphaMax-alpha);
volScalaField sourceRCc = exp(-Kc*pow(phaseb-.rho(),1.0/2.0)
*pow(dBuff+plusDS,5.0/6.0)*pow(epsilon,1.0/3.0)/pow(phaseb-.sig(),1.0/2.0));
volScalarField sourceRC = sourceRCa/sourceRCb*sourceRCc;

return sourceRCa*scalar(1.0);
}
mehrdad_kbg is offline   Reply With Quote

Old   May 20, 2013, 21:23
Default
  #2
Senior Member
 
kmooney's Avatar
 
Kyle Mooney
Join Date: Jul 2009
Location: Amherst, MA USA - San Diego, CA USA
Posts: 285
Rep Power: 9
kmooney is on a distinguished road
Honestly you might be better off just contacting the author of the thesis you got the source code from. If they are willing to publish the C file I can't imagine they wouldn't send you the header if you asked.
kmooney 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
funkySetFields compilation error tayo OpenFOAM 39 December 3, 2012 06:18
problem on installing swak4Foam navid2 OpenFOAM Installation 2 May 30, 2012 04:32
error message with modeling a cube with a hold at the center hsingtzu OpenFOAM Native Meshers: blockMesh 2 March 14, 2012 10:56
Problem installing on Ubuntu 9.10 -> 'Cannot open : No such file or directory' mfiandor OpenFOAM Installation 2 January 25, 2010 10:50
Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00


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