|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Vitor Braga
Join Date: Feb 2019
Posts: 6
Rep Power: 7 ![]() |
Hello! So, as far as I know, when the Mixture model is selected for a two-phase flow, Fluent calculates the viscosity of the mixture with the equation:
mu_mixture = alpha_1 * mu_1 + alpha_2 * mu_2. I would like to change the viscosity calculation to more complex models, such as the one proposed by Dukler (1964). Does anyone know how to proceed? What kind of UDF do I need? Where to implement it? I really appreciate any help I can get! Thank you very much. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
reza
Join Date: Aug 2019
Posts: 14
Rep Power: 7 ![]() |
Hi,
did you find how can you solve this problem? if yes, may i ask you please shear it with me? thanks |
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Vitor Braga
Join Date: Feb 2019
Posts: 6
Rep Power: 7 ![]() |
Yes, I found a solution. You don't need to change the way Fluent calculates the mixture viscosity. You can just change by UDF the viscosity of each phase to represent the viscosity model that you want. Thus, each phase will have the same value of viscosity (calculated based on the specific model that you want to implement), hence when Fluent calculates:
mu_mixture = alpha_1 * mu_1 + alpha_2 * mu_2 the value it obtains is the value of your model, since mu_1 = mu_2 (and alpha_2 = 1 - alpha_1). I'm sorry if I didn't make myself clear... English is not my native language. If you need any clarifications, please, feel free to ask, I'll try to improve my explanation. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
New Member
reza
Join Date: Aug 2019
Posts: 14
Rep Power: 7 ![]() |
Hi,
Really thank you, another question: you used same udf for each phase? may i ask you please if you have your udf, share with me please? |
|
![]() |
![]() |
![]() |
![]() |
#5 |
New Member
reza
Join Date: Aug 2019
Posts: 14
Rep Power: 7 ![]() |
did you use DEFINE_PROFILE for your work? my main question is that can we write UDF as we can write for single-phase after we can implement this UDF for each phase??
|
|
![]() |
![]() |
![]() |
![]() |
#6 |
New Member
Vitor Braga
Join Date: Feb 2019
Posts: 6
Rep Power: 7 ![]() |
Yes, I used the same UDF.
This is the UDF: #include"udf.h" #include "math.h" #define mu_phase1 0.004452 #define mu_phase2 1.208e-05 DEFINE_PROPERTY(beatie_whalley,c,t_phase) { real mu_lam; real vof1,vof2; Thread *tc = THREAD_SUPER_THREAD(t_phase); Thread **pt = THREAD_SUB_THREADS(tc); Thread *tp = pt[0]; vof1 = C_VOF(c,tp); vof2 = 1-C_VOF(c,tp); mu_lam = vof2 * mu_phase2 + mu_phase1 * vof1 * (1 + 2.5 * vof2); return mu_lam; } Yes, you use the UDF for single phase in each phase for the viscosity. |
|
![]() |
![]() |
![]() |
![]() |
#7 |
New Member
reza
Join Date: Aug 2019
Posts: 14
Rep Power: 7 ![]() |
really thanks
|
|
![]() |
![]() |
![]() |
![]() |
#8 |
New Member
reza
Join Date: Aug 2019
Posts: 14
Rep Power: 7 ![]() |
Hi, I have a question, I would be grateful if you could help me.
for my problem: VOF_fluid = 0.8 VOF_solid = 0.2 base on your code: vof1 = 0.8 vof2 = 0.2 ???? or vof1 = 0.2 vof2 = 0.8 ??? or if I give my UDF to each phase, automatically VOF is equal to volume fraction of that phase? really thanks for your great help |
|
![]() |
![]() |
![]() |
![]() |
#9 |
New Member
Vitor Braga
Join Date: Feb 2019
Posts: 6
Rep Power: 7 ![]() |
Well, it depends. VOF1 is the primary phase, which you defined in Fluent.
|
|
![]() |
![]() |
![]() |
Tags |
mixture model, multiphase flow, two-phase flow, viscosity model |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
What viscosity to use for solid phase in VOF model | redwanamit034 | FLUENT | 0 | June 14, 2017 17:09 |
Using a new implemented viscosity model with simpleFoam | TemC | OpenFOAM Running, Solving & CFD | 6 | March 8, 2017 03:07 |
Enabling Open Channel Flow Sub-Model in Mixture model | cod213 | FLUENT | 0 | January 10, 2017 13:40 |
Overflow Error in Multiphase Modelling with Two Continuous Fluids | ashtonJ | CFX | 6 | August 11, 2014 14:32 |
species transport model or mixture model? | achaokaoyan | Main CFD Forum | 0 | July 10, 2010 10:52 |