CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

twoPhaseEulerFoam & Source term

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

LinkBack Thread Tools Display Modes
Old   September 10, 2013, 06:26
Default twoPhaseEulerFoam & Source term
Senior Member
Aurelien Thinat
Join Date: Jul 2010
Posts: 165
Rep Power: 9
Aurelien Thinat is on a distinguished road
Hello everyone,

I'm willing to include a mass source term in twoPhaseEulerFoam solver :

S_Mass = - rho1 * alpha1 / dT, with dT the time step.

It should extract the mass of the phase1 (dispersed phase) in a zone I define with a function :
F = 1 if I activate the source term,
0 elsewhere.

That's for the theory. I'm stuck with the OpenFOAM code now.

Since the equation is in volumic form, I wrote :

alphaEqn :
- fvc::ddt(alpha1)
- fvm::laplacian

pEqn :
fvm::laplacian(Dp,p) == fvc::div(phiHByA) + F*alpha1/(runTime.deltaT().value)

Does it seems correct to you ?

Thank you.
Aurelien Thinat is offline   Reply With Quote


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
transient source term strakakl OpenFOAM 38 November 19, 2013 02:18
Problem of SOURCE term gradient in UDS wind Fluent UDF and Scheme Programming 5 June 21, 2013 05:39
build problem swak4Foam OF 2.2.0 mcathela OpenFOAM Installation 14 April 23, 2013 13:59
Compiling gmshFoam with OpenFOAM-1.5 BlGene Open Source Meshers: Gmsh, Netgen, CGNS, ... 10 August 6, 2009 04:26
UDFs for Scalar Eqn - Fluid/Solid HT Greg Perkins FLUENT 0 October 11, 2000 03:43

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