CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Fluent UDF and Scheme Programming (https://www.cfd-online.com/Forums/fluent-udf/)
-   -   Particle transport and deposition & Drift Flux Model (https://www.cfd-online.com/Forums/fluent-udf/73985-particle-transport-deposition-drift-flux-model.html)

 Peter023 March 22, 2010 01:01

Particle transport and deposition & Drift Flux Model

Hello all:

please, could anyone help me with the, so called, Drift Flux Model which is, among another tools for the evaluation of the particle transport and deposition onto different surfaces, a possible solution of the particle transport phenomena. Unfortunately, the original model and other improved models depicted in the literature are not clearly explained and I have a real problem with their transformation into the UDF.
the literature, for example, is :
http://www.sciencedirect.com/science...63e83283d44ef9

or:

http://www.sciencedirect.com/science...1128c844815001

If somebody has the solution or some idea or source code, please let me know :)

Pete.

 Peter023 March 22, 2010 01:07

supplementary code

As a source I'm trying to use this:

real vs = (9.81 * (1400.0 - C_R(c,t)) * dp * dp) / (18.0 * C_MU_EFF(c,t));
source = (vs + C_W(c,t)) * C_UDSI(c,t,0);
dS[eqn] = 0.0;

As a diffusivity term (DEFINE_DIFFUSIVITY):
D = (Kb * 300.0 * Cc(dp, Lambda)) / (3.0 * 3.14 * C_MU_EFF(cell,thread) * dp); // Brownian diffusivity

and deposition onto surfaces:

{
face_t f;
cell_t c;
real jay_a = 0.0, jay_b = 0.0;
real vs = 0.0; // Settling velocity
real v_dd = 0.0; // Dep. velocity, downward horiz. surface
real D = 0.0; // Brownian diffusivity

{

vs = (9.81 * (1400.0 - C_R(c,t)) * dp * dp) / (18.0 * C_MU_EFF(c,t)); // Settling vel.

if (selector)
{
v_dd = vs / (exp(vs * integral / frict_vel) - 1.0); // Deposition vel.
jay_b = sign * (v_dd * C_UDSI(c,t,0));
}else
{
D = (Kb * 300.0 * Cc(dp, Lambda)) / (3.0 * 3.14 * C_MU_EFF(c,t) * dp); // Brownian diffusivity
jay_a = sign * ((D + C_MU_T(c,t)) * C_UDSI_G(c,t,0)[2]) + vs * C_UDSI(c,t,0);
}
}