请问关于UDS源项和密度UDF的问题，

ayeelan

遇到的问题是，VOF模型，两相，空气和自定义液体，密度单独可以运行；UDS0单独不可以运行，两个同时设置的时候，不能运行；提示如下的错误：
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()
请问版上可能指点一下？
我的UDS和UDF如下：
DEFINE_SOURCE(uds_source_W,cell,thread,dS,eqn)
{
           int phase_domain_index=0;
        Thread *mixture_thread;
        Thread *foam_thread=THREAD_SUB_THREAD(mixture_thread,phase_domain_index);
        real aF=C_VOF(cell,foam_thread);
        real R_F=C_R(cell,foam_thread);
        real T=C_T(cell,thread);
        real Xw=C_UDSI(cell,thread,0);

        real Aw=25.859;
        real Ew=4.412e4;
        real R=8.314472;
        real Rbaliq=741.8;
        real Rpu=1100.;
        real Lliq=0.005;

        real source_W;
        source_W=aF*R_F*Aw*exp(-Ew/(R*T))*(1-Xw)*(1/(1+Lliq*Rpu/Rbaliq));
        dS[eqn]=-aF*R_F*Aw*exp(-Ew/(R*T))*(1/(1+Lliq*Rpu/Rbaliq));
//        dS[eqn]=0;

        return source_W;
}

DEFINE_PROPERTY(cell_density_test,cell,thread)
{

        real T=C_T(cell,thread);
        real P=C_P(cell,thread)+101325;
        real Xw=C_UDSI(cell,thread,0);
        real R=8.314472;
        real Mh2o=18.;
        real xw0=0.02;
        real density_liq22=(Xw*R*T*xw0*1000/(P*Mh2o));
       
        return density_liq22;
}

[ 本帖最后由 ayeelan 于 2010-12-28 23:08 编辑 ]

ayeelan

DEFINE_SOURCE(source01,cell,thread,dS,eqn)
{
           int phase_domain_index=0;
        Thread *foam_thread=THREAD_SUB_THREAD(thread,phase_domain_index);

        real R_F=C_R(cell,foam_thread);
        real aF=C_VOF(cell,foam_thread);

      ······
}

VOF中，thread是混合相吧，在C_R(cell,foam_thread)中使用foam_thread的话就会错误，请问这个主相foam_thread的使用有什么错误？

zzulxl1977

楼主做的课题与我做的工作比较接近，如果方便希望能向您请教

没晓晓

楼主解决了吗？我遇到了同样的问题，求帮助