[交流] 用FLUETN UDS解电场 电荷方程 已有2人参与

用uds解两方程：电场方程与电荷密度方程 http://muchongimg.xmcimg.com/data/b ... _1453125188_614.png http://muchongimg.xmcimg.com/data/b ... _1453125188_698.png 结果与实际相差甚远，求大神指点，程序如下： #include "math.h" #include "udf.h"   #define DIELECTRIC_CONSTANT 7.0e-10 /*电荷密度方程*/ DEFINE_UDS_FLUX(current_density,f,t,i) { real NV_VEC(current); real NV_VEC(A); real NV_VEC(ion_vel); real NV_VEC(C_EE); real flux; cell_t c0; Thread *t0; c0=F_C0(f,t); t0=THREAD_T0(t); C_UDMI(c0,t0,0)=-C_UDSI_G(c0,t0,0)[0]; C_UDMI(c0,t0,1)=-C_UDSI_G(c0,t0,0)[1]; C_UDMI(c0,t0,2)=-C_UDSI_G(c0,t0,0)[2]; NV_D(ion_vel,=,C_U(c0,t0),C_V(c0,t0),C_W(c0,t0)); NV_D(C_EE,=,C_UDMI(c0,t0,0),C_UDMI(c0,t0,1),C_UDMI(c0,t0,2)); C_UDMI(c0,t0,3)=NV_MAG(C_EE); NV_V_VS(current,=,ion_vel,+,C_EE,*,2.3e-5); F_AREA(A,f,t); flux=NV_DOT(current,A); return flux; } /*电势方程源项*/ DEFINE_SOURCE(voltage_source,c,t,dS,eqn) { real vol; vol=C_UDSI(c,t,1)/DIELECTRIC_CONSTANT; dS[eqn]=0; return vol; } QQ截图20160118214105.png QQ截图20160118214116.png

回复此楼