[vmbraga]

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.

[rezaamininiaki]

Hi,
did you find how can you solve this problem?
if yes, may i ask you please shear it with me?
thanks

[vmbraga]

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.

[rezaamininiaki]

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?

[rezaamininiaki]

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??

[vmbraga]

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.

[rezaamininiaki]

really thanks

[rezaamininiaki]

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

[vmbraga]

Well, it depends. VOF1 is the primary phase, which you defined in Fluent.