[zhandos_baizhuma]

Dear all!
I am trying to simulate Eulerian multi phase flow around the airfoil and I need to add source terms to airfoil wall.
see Source terms and my udf file in attachments.
As example I took continuity equation
In source terms:
αd the droplet volume fraction, ρd the droplet density and A the face area vector and vd the droplet impingement vector.
I am new in UDF writing, so please can you give me any feedback,
Best
Zhandos

[AlexanderZ]

your code

#include "udf.h"
#include "mem.h"
#include "metric.h"

#define d=1000;

DEFINE_SOURCE (continuity,c,t,ds,eqn)
{ 
	real NV_VEC(v);
	real NV_VEC(A);
	F_Area(A,f,t);
	source = -NV_DOT(v,A)*C_VOF(c,t)*d;
	ds[eqn]=0;
	return source;
}

1. v is not defined
2. what is the unit of source you want to apply?
3. to use F_Area(A,f,t); macro you need face, but DEFINE_SOURCE makes a loop over cells (not faces) -> f is not defined
4. to define f inside DEFINE_SOURCE macro is challenging (from my point of view)

[zhandos_baizhuma]

Alexander Z, thank you for reply.
Actually I need to add 4 source terms, such as in picture

v - droplet velocity vector, I thought that Fluent will pass value from inlet

if face area vector is not applicable, maybe use x[ND_ND]?

[AlexanderZ]

what is the unit of source you want to apply?

[zhandos_baizhuma]

kilogram per second

[AlexanderZ]

DEFINE_SOURCE macro applying volumetric source W/m3 or kg/m3

you can get more information in Ansys Fluent customization manual

[zhandos_baizhuma]

I am still trying to solve Eulerian Muldtiphase with UDF
I made some improvement to my UDF. Please can you give some feedback on it?

#include "udf.h"

DEFINE_SOURCE(phase_cont,c,t,dS,eqn)
{
real con, source;
real NV_VEC(v), NV_VEC(A);
int ro=1000;
NV_D(v, =, C_U(c,t), C_V(c,t), C_W(c,t));
F_AREA(A,c,t);
con = NV_DOT(v,A);
source = -con*ro*C_VOF(c,t);
dS[eqn] = 0;
return source;
}

[AlexanderZ]

code seems to be correct, what problems do you have with it?

[zhandos_baizhuma]

It loads normally. But just after initialization there is a message in TUI

999999: mpt_accept: error: accept failed: No error

999999: mpt_accept: error: accept failed: No error
MPI Application rank 0 exited before MPI_Finalize() with status 2
The fl process could not be started.

I can not figure out why it happens

[AlexanderZ]

initialize first, then load UDF