FLUENT中的degassing边界条件UDF

Example2-DegassingBoundaryConditionThefollowingUDFsareusedtodefinethebottomsurfaceasastandardvelocityinletforthegas(primary)phase.TheinletVOFofthedropletphaseis0andanegativesourcetermforsecondaryphasemassconservationissetforthelayerofcellsnexttotheoutlet.Thesourceremovesallsecondaryphasemassinthecellduringonetimestep.Therecoilforceduetothemasssourceisalsocalculated./*******************************************************************//*ThisUDFisanimplementationofthedegassingboundarycondition*//*******************************************************************/#includeudf.h#includesg.h#includesg_mphase.h#includeflow.h#includemem.h#includemetric.hDEFINE_SOURCE(degassing_source,cell,thread,dS,eqn){realsource;Thread*tm=THREAD_SUPER_THREAD(thread);source=-C_R(cell,thread)*C_VOF(cell,thread)/CURRENT_TIMESTEP;C_UDMI(cell,tm,0)=source;dS[eqn]=-C_R(cell,thread)/CURRENT_TIMESTEP;returnsource;}DEFINE_SOURCE(x_prim_recoil,cell,tp,dS,eqn){realsource;Thread*tm=THREAD_SUPER_THREAD(tp);Thread*ts;ts=THREAD_SUB_THREAD(tm,1);source=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_U(cell,tp);dS[eqn]=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP;returnsource;}DEFINE_SOURCE(x_sec_recoil,cell,ts,dS,eqn){realsource;Thread*tm=THREAD_SUPER_THREAD(ts);source=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_U(cell,ts);dS[eqn]=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP;returnsource;}DEFINE_SOURCE(y_prim_recoil,cell,tp,dS,eqn){realsource;Thread*tm=THREAD_SUPER_THREAD(tp);Thread*ts;ts=THREAD_SUB_THREAD(tm,1);source=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_V(cell,tp);dS[eqn]=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP;returnsource;}DEFINE_SOURCE(y_sec_recoil,cell,ts,dS,eqn){realsource;Thread*tm=THREAD_SUPER_THREAD(ts);source=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_V(cell,ts);dS[eqn]=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP;returnsource;}DEFINE_SOURCE(z_prim_recoil,cell,tp,dS,eqn){realsource;Thread*tm=THREAD_SUPER_THREAD(tp);Thread*ts;ts=THREAD_SUB_THREAD(tm,1);source=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_W(cell,tp);dS[eqn]=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP;returnsource;}DEFINE_SOURCE(z_sec_recoil,cell,ts,dS,eqn){realsource;Thread*tm=THREAD_SUPER_THREAD(ts);source=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_W(cell,ts);dS[eqn]=-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP;returnsource;}