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

discretizing void fraction transport equation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 9, 2013, 11:40
Default discretizing void fraction transport equation
  #1
New Member
 
Join Date: Sep 2010
Location: Windsor, Ontario
Posts: 17
Rep Power: 6
mehrdad_kbg is on a distinguished road
Hi everybody,

I am having a really hard time for discretizing the transport equation for void fraction in OpenFOAM. The equation is as following:

d(alpha)/dt+U.Δ(alpha)=Δ.J

and J is a function of alpha. I appreciate it if you help me to write it in OpenFOAM format.

Thanks.
mehrdad_kbg is offline   Reply With Quote

Old   September 10, 2013, 02:23
Default
  #2
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
What did you try up until now? Where are you stuck exactly? Have a look at $FOAM_SOLVERS/basic/scalarTransportFoam
Bernhard is offline   Reply With Quote

Old   September 10, 2013, 13:25
Default
  #3
New Member
 
Join Date: Sep 2010
Location: Windsor, Ontario
Posts: 17
Rep Power: 6
mehrdad_kbg is on a distinguished road
Thanks for your response. I describe for you the whole story. The equation that I have put in my first post d(alpha)/dt+U.Δ(alpha)=Δ.J is the governing equation for the transport of particles in a suspension due to normal stresses (Δ is the divergence operator here). The normal stress is the J term in the right hand side.

The process that I have done to reach the J equation is:

SP=-etha0*gamma*ethaN*Q

Here, etha0 and ethaN are viscosity, gamma is the magnitude of the strain rate and Q is a tensor. In OpenFOAM format it is :

volTensorField SP = -1.0*etha0*sqrt(magSqr(symm(fvc::grad(U))))*ethan*Q

To find J, i have done:

volVectorField divSP = fvc::div(SP);

volVectorField J= (-2.0*aa*aa/9.0/etha0)*fphi*divSP;

As you see J is a vector which is calculated based on a divergence process on a tensor field. I do not know if the divergence calculation is correct or not. Having the J vector field, I have written the transport equation of void fraction:


fvScalarMatrix alphaEqn
(
fvm::ddt(alpha)
+ fvm::div(phic, alpha, scheme)
- fvc::div(J)
);


As you see another divergence is applied to J. So now I do not know
first: if the two divergence process that I have done are written correctly in OpenFOAM format or not,

second: is my transport equation which is the discretized form of the equation presented in my first post correct or not?

I appreciate if any body can help me.
mehrdad_kbg is offline   Reply With Quote

Old   September 11, 2013, 08:27
Default
  #4
New Member
 
Join Date: Sep 2010
Location: Windsor, Ontario
Posts: 17
Rep Power: 6
mehrdad_kbg is on a distinguished road
any answer?!
mehrdad_kbg 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
Transport equation in CFX D.ko CFX 0 August 9, 2013 07:57
Sub-iterations for transport equation in drift-flux Siemen OpenFOAM Programming & Development 2 July 5, 2013 06:06
how to make a scalar transport equation to be conserved? immortality OpenFOAM Programming & Development 3 May 5, 2013 08:11
Implement species transport equation Tobi OpenFOAM Programming & Development 0 June 2, 2012 13:26
Source Term on Scalar Transport Equation - crashing simulation alessio.nz OpenFOAM 19 February 8, 2011 08:38


All times are GMT -4. The time now is 17:36.