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

Eulerian multiphase tracer transport - value of scalar?

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

LinkBack Thread Tools Display Modes
Old   January 14, 2015, 10:44
Default Eulerian multiphase tracer transport - value of scalar?
New Member
Dominik Schmidt
Join Date: Mar 2014
Posts: 11
Rep Power: 5
dschmidt is on a distinguished road
Dear Foamers,

I believe the answer to my question is quite simple, but I couldn't find a clear definition for it:

I'm adding a scalar transport equation to an incompressible multiphase solver, where the scalar should only be transported by/in the liquid phase,
representing e.g. a solution of salt in water.

\frac{\partial {\alpha_L Tr}}{\partial {t}} + \nabla \cdot (\alpha_L U_L Tr) = \nabla \cdot (\alpha_L \Gamma) \nabla Tr)

fvm::ddt(alpha1, Tracer)    
+ fvm::div(fvc::interpolate(alpha1)*phi1, Tracer)
- fvm::laplacian(fvc::interpolate(alpha1)*DT, Tracer)
Assuming this part is correct (please correct me, if it isn't),
my questions is about the "value" of the transported scalar:

1.) "Tracer" = mass(Tracer) per cell-volume (= concentration in each cell)?

... so to get the concentration in the liquid phase in each cell one has to divide "Tracer"-value by alpha1 (volumetric phase fraction of liquid phase)?


2.) the value for "Tracer" already represents the concentration in the liquid phase?


3.) anything else?

In short: Does the Tracer-value represent a concentration based on the cell volume or based on the liquid volume in each cell?
Or is a completely different approach necessary to get the concentration of e.g. a salt in the water phase for the evaluation of homogeneity inside the domain.

Thanks for your help!


Reading helped me to understand that the second option should be correct.

Last edited by dschmidt; January 19, 2015 at 04:10.
dschmidt is offline   Reply With Quote

Old   February 2, 2015, 00:21
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,911
Rep Power: 28
alberto will become famous soon enoughalberto will become famous soon enough
Yes, your equation is correct and refers to a tracer into the liquid phase (note that alpha_L is in all the terms).
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status:
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
alberto 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
problems concerning mass conservativity in bubbleFoam with custom scalar transport cutter OpenFOAM Programming & Development 3 February 10, 2015 05:25
Diverging solution in transonicMRFDyMFoam tsalter OpenFOAM Running, Solving & CFD 30 July 7, 2014 06:20
compressible flow in turbocharger riesotto OpenFOAM 50 May 26, 2014 01:47
whats the cause of error? immortality OpenFOAM Running, Solving & CFD 11 April 22, 2014 12:32
is internalField(U) equivalent to zeroGradient? immortality OpenFOAM Running, Solving & CFD 7 March 29, 2013 02:27

All times are GMT -4. The time now is 19:51.