CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Equation in Burns.C of Turbulent dispersion model of Burns et al.

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

Like Tree4Likes
  • 4 Post By geth03

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 22, 2020, 21:05
Post Equation in Burns.C of Turbulent dispersion model of Burns et al.
  #1
New Member
 
Chaoran Guan
Join Date: Aug 2019
Posts: 5
Rep Power: 7
Chaoran is on a distinguished road
Hello, Foamers. I was using the reactingTwoPhaseEulerFoam of OF6.0 Version.
A turbulent dispersion model called Burns is used which located in applications/solvers/multiphase/reactingEulerFoam/interfacialModels/turbulentDispersionModels/Burns. And the reference is

Burns, A. D., Frank, T., Hamill, I., & Shi, J. M. (2004, May).The Favre averaged drag model for turbulent dispersion in Eulerian multi-phase flows.In 5th international conference on multiphase flow, ICMF (Vol. 4, pp. 1-17). https://www.researchgate.net/publica...ti-Phase_Flows

The equation used in Burns.C from the reference is
\vec M_{\alpha}^{TD}=-\vec M_{\beta}^{TD}=-\frac{3}{4}C_D\frac{\overline r_{\beta}\rho_{\alpha}}{d_\beta}|\vec{U_\beta}-\vec{U_\alpha}|\frac{\nu_{t\alpha}}{\sigma_{r\alpha}}(\frac{1}{\overline r_{\alpha}}+\frac{1}{\overline r_{\beta}})\nabla \overline r_{\alpha}
And the correspponding source code in Burns.C from OF version 6.0 is
Code:
        0.75
       *drag.CdRe()
       *pair_.continuous().nu()
       *continuousTurbulence().nut()
       /(
            sigma_
           *sqr(pair_.dispersed().d())
        )
       *pair_.continuous().rho()
       *pair_.dispersed()
       *(
           1.0/max(pair_.dispersed(), residualAlpha_)
         + 1.0/max(pair_.continuous(), residualAlpha_)
        );
I tried to translate it into the mathematical form as
\frac{3}{4}C_D\frac{\nu_\alpha\rho_\alpha \overline r_\beta}{d^2_\beta}\frac{\nu_{t\alpha}}{\sigma_{r\alpha}}(\frac{1}{\overline r_{\alpha}}+\frac{1}{\overline r_{\beta}})
The two equations seem not to be the same. The interphase velocity term |\vec{U_\beta}-\vec{U_\alpha}| is missed in the source code and there is an additional term \frac{\nu_\alpha}{d_{\beta}} in it.
Is there any relationship between these two terms?
Chaoran is offline   Reply With Quote

Old   November 24, 2020, 03:43
Default
  #2
Senior Member
 
Join Date: Dec 2019
Location: Cologne, Germany
Posts: 365
Rep Power: 8
geth03 is on a distinguished road
drag.CdRe() in the code means C_d x Re (Reynolds number),
it cancels out the additional terms (kin. viscosity / bubble diameter)
and leaves your missing term (velocity difference).

so in the end it is implemented correctly although not in a way that
makes sense at first sight.
geth03 is offline   Reply With Quote

Reply

Tags
multiphase flows, multiphase modelling

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Domain Reference Pressure and mass flow inlet boundary AdidaKK CFX 75 August 20, 2018 05:37
How to use the turbulent dispersion in FLUENT? marialhm FLUENT 2 October 31, 2017 20:55
Confusion on the FAD (Burns) turb-dispersion model spy_hzdr Main CFD Forum 0 September 29, 2017 04:27
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 05:36
A new turbulent model without empirical coefficients s9a9m92001 OpenFOAM Programming & Development 6 October 9, 2013 23:11


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