|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Join Date: Oct 2016
Posts: 3
Rep Power: 9 ![]() |
Hi all,
i want to implement to mass flow rate of a secondary phase (eulerian model) through a given face by using a DEFINE_UDS_FLUX, in Fluent. Basically I want to use a code similar to the one in the Fluent manual (attached below), but I want the velocity for the flow to be the average between the velocity of the primary phase and the velocity of the secondary phase. I do not understand how to define the threads for each phase. Can you help me? /************************************************** ********************/ /* UDF that implements a simplified advective term in the */ /* scalar transport equation */ /************************************************** ********************/ #include "udf.h" DEFINE_UDS_FLUX(my_uds_flux,f,t,i) { cell_t c0, c1 = -1; Thread *t0, *t1 = NULL; real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0; c0 = F_C0(f,t); t0 = F_C0_THREAD(f,t); F_AREA(A, f, t); /* If face lies at domain boundary, use face values; */ /* If face lies IN the domain, use average of adjacent cells. */ if (BOUNDARY_FACE_THREAD_P(t)) /*Most face values will be available*/ { real dens; /* Depending on its BC, density may not be set on face thread*/ if (NNULLP(THREAD_STORAGE(t,SV_DENSITY))) dens = F_R(f,t); /* Set dens to face value if available */ else dens = C_R(c0,t0); /* else, set dens to cell value */ NV_DS(psi_vec, =, F_U(f,t), F_V(f,t), F_W(f,t), *, dens); flux = NV_DOT(psi_vec, A); /* flux through Face */ } else { c1 = F_C1(f,t); /* Get cell on other side of face */ t1 = F_C1_THREAD(f,t); NV_DS(psi_vec, =, C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),*,C_R(c0,t0)); NV_DS(psi_vec, +=, C_U(c1,t1),C_V(c1,t1),C_W(c1,t1),*,C_R(c1,t1)); flux = NV_DOT(psi_vec, A)/2.0; /* Average flux through face */ } /* ANSYS Fluent will multiply the returned value by phi_f (the scalar’s value at the face) to get the ‘‘complete’’ advective term. */ return flux; } [/COLOR] Last edited by Waid; April 3, 2019 at 11:27. |
|
![]() |
![]() |
![]() |
Tags |
eulerian model, flux, multiphase model, uuds |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
kklOmega omega boundary condition at the wall - OF3.0 | Artur | OpenFOAM Running, Solving & CFD | 1 | June 17, 2016 15:42 |
Turbulence model choice | dave13 | CFX | 3 | December 22, 2015 07:12 |
problem with solving lagrange reaction cloud | Polli | OpenFOAM Running, Solving & CFD | 0 | April 30, 2014 07:53 |
manualInjection model in sprayFoam | Mentalo | OpenFOAM Running, Solving & CFD | 1 | April 2, 2014 09:29 |
Problems bout CFD model of biomass gasification, Downdraft gasifier | wanglong | FLUENT | 2 | November 25, 2009 23:27 |