help~~UDF on solid pressure

Yurong

Dear all:

I would like to revise the particle pressure in gas-solid two phase flow with Euler-Euler model. Could you please help me to look at the following udf codes about solid pressure? I don't think it is correct coz no right results come out.

I am very grateful for your help in advance.

#include "udf.h"

#define rou_s 1000

DEFINE_PROPERTY(solid_pressure,c,t) { real Ps,alfa_s,theta_s,e_ss,g0_ss;

alfa_s = C_VOF(c,t); theta_s = PROP_granular_t; e_ss = PROP_restitution_coeff; g0_ss = PROP_mpg_radial_f;

Ps = alfa_s*rou_s*theta_s + 2.*rou_s*(1.+e_ss)*alfa_s*alfa_s*g0_ss*theta_s ; return Ps; }

Yurong