UDF for Contact Angle
 

Hi,

I want to write a UDF for defining Contact angle boundary condition for the case of 2 phase problem solved with unsteady VOF method. My requirement is to give input to the case file with the updated contact angle at evry new time step.

I am including the code i have written...i appreciate if any one can find mistake in the code......

#include "udf.h"


DEFINE_PROFILE(contact_angle,c, mixture_thread)

{

real theta;

float current_time = CURRENT_TIME;

if (current_time = 0)

theta = 140;

if (0 < current_time <= 0.00005)

theta = 130;

if (0.00005 < current_time <= 0.0001)

theta = 125;


else

theta = 135;

}

sorry, it doesn`t work T-T