24小时热门版块排行榜

>论坛更新日志 (3725)
>虫友互识 (614)
>文献求助 (214)
>休闲灌水 (213)
>导师招生 (198)
>考博 (129)
>论文投稿 (118)
>博后之家 (90)
>硕博家园 (81)
>招聘信息布告栏 (71)
>论文道贺祈福 (68)
>基金申请 (57)
>教师之家 (49)
>考研 (47)
>找工作 (45)
>公派出国 (21)

返回列表

左岸流年

新虫 (正式写手)

应助: 2 (幼儿园)
金币: 187.9
散金: 213
红花: 2
帖子: 368
在线: 98.8小时
虫号: 2118709
注册: 2012-11-10
专业: 传热传质学

[求助] 关于udf 变量梯度的问题已有1人参与

请问 {
  Alloc_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
  Scalar_Reconstruction(pDomain, SV_VOF,-1,SV_VOF_RG,NULL);
  Scalar_Derivatives(pDomain,SV_VOF,-1,SV_VOF_G,SV_VOF_RG,Vof_Deriv_Accumulate);
}

{
   Alloc_Storage_Vars(domain, SV_T_RG, SV_T_G,  SV_NULL);
   T_derivatives(domain);

这里面的变量什么意思  该如何修改
/**************************************************************/
/* UDF for specifying an interfacail area density          */
/**************************************************************/

DEFINE_ADJUST(area_density, domain)
{
  Thread *t;
  Thread **pt;
  cell_t c;
  Domain *pDomain = DOMAIN_SUB_DOMAIN(domain,P_PHASE);
  real voidx, voidy, voidz=0;


  {
  Alloc_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
  Scalar_Reconstruction(pDomain, SV_VOF,-1,SV_VOF_RG,NULL);
  Scalar_Derivatives(pDomain,SV_VOF,-1,SV_VOF_G,SV_VOF_RG,Vof_Deriv_Accumulate);
}

{
   Alloc_Storage_Vars(domain, SV_T_RG, SV_T_G,  SV_NULL);
   T_derivatives(domain);
   Free_Storage_Vars(domain, SV_T_RG, SV_NULL);
}

   mp_thread_loop_c (t,domain,pt)
if (FLUID_THREAD_P(t))
   {
      Thread *tp = pt[P_PHASE];

      begin_c_loop (c,t)
      {

#if RP_3D
C_UDMI(c,t,0) = (C_VOF_G(c,tp)[0]*C_T_G(c,t)[0]+C_VOF_G(c,tp)[1]*C_T_G(c,t)[1]+C_VOF_G(c,tp)[2]*C_T_G(c,t)[2]);
#endif

#if RP_2D
C_UDMI(c,t,0) = (C_VOF_G(c,tp)[0]*C_T_G(c,t)[0]+C_VOF_G(c,tp)[1]*C_T_G(c,t)[1]);
#endif

      }
      end_c_loop (c,t)
      }

Free_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
Free_Storage_Vars(domain, SV_T_G, SV_NULL);
}

回复此楼