Tianxu
Subscriber
Here is my parrallel inlet UDF:n#include udf.hnDEFINE_PROFILE(inlet_velocity,th,i)n{nnface_t f;nreal po;nreal pcr;nreal To;nreal gamma;nreal R;nreal p_cur;npo = 1013250; ngamma = 1.29;npcr = po * pow(2/(gamma + 1), gamma/(gamma-1)); nTo = 300;  nR = 368.9; nn#if !RP_HOSTn    n    begin_f_loop(f,th) n    {nn        p_cur = F_P(f,th); nn        if(p_cur >= po)n        F_PROFILE(f,th,i) = 0.0000001;nn        else if(p_cur <= po && p_cur > pcr)n        F_PROFILE(f,th,i) = sqrt( 2*gamma/(gamma-1) *R*To * (1 - pow(p_cur/po, gamma-1/gamma)) );nn        else if(p_cur <= pcr)n        F_PROFILE(f,th,i) = sqrt(2*gamma/(gamma+1) *R*To); n        n    }n    end_f_loop(f,th)n#endifnn}n