CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新
    1. 主页
    2. 冠竹
    冠
    • 资料
    • 关注 1
    • 粉丝 0
    • 主题 43
    • 帖子 111
    • 群组 0

    冠竹

    @冠竹

    38
    资料浏览
    111
    帖子
    0
    粉丝
    1
    关注
    注册时间 最后登录

    冠竹 取消关注 关注

    冠竹 发布的最新帖子

    • RE: fluent中在动量方程中添加弹性项以及添加5个用户自定义标量函数UDS后

      @李东岳谢谢老师的建议,我担心学习of的时间不太够

      发布在 Fluent
      冠
      冠竹
    • RE: fluent中在动量方程中添加弹性项以及添加5个用户自定义标量函数UDS后

      @bestucan 谢谢老师,我还想问一下我uds方程里对流项没有密度,我直接把udf帮助文档里面对流项的代码中密度的地方改为1可以吧?
      6dd6d91a-8e8a-46f6-b499-d4ce49fde849-68216d0572e0cd03f7d2bae473bb343.png

      DEFINE_UDS_FLUX(No_Cxx_FLUX,f,t,i)                 //Cxx对流项    
      {
          	cell_t c0, c1 = -1;
          	Thread *t0, *t1 = NULL;
          	real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;
          	c0 = F_C0(f,t);
          	t0 = F_C0_THREAD(f,t);
          	F_AREA(A, f, t);
          
         	if (BOUNDARY_FACE_THREAD_P(t)) 
          	{
             		real dens;
             
             		if (NNULLP(THREAD_STORAGE(t,SV_DENSITY)))
               		dens = F_R(f,t);  
             		else
               		dens = C_R(c0,t0); 
             		NV_DS(psi_vec, =, F_U(f,t), F_V(f,t), F_W(f,t), *, 1);
            		flux = NV_DOT(psi_vec, A); 
      	}
          	else
          	{
           		c1 = F_C1(f,t);   
           		t1 = F_C1_THREAD(f,t);
           		NV_DS(psi_vec, =, C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),*,1);
           		NV_DS(psi_vec, +=, C_U(c1,t1),C_V(c1,t1),C_W(c1,t1),*,1);
           		flux = NV_DOT(psi_vec, A)/2.0;
          	}
      	return flux;
      } 
      
      
      发布在 Fluent
      冠
      冠竹
    • RE: fluent中在动量方程中添加弹性项以及添加5个用户自定义标量函数UDS后

      不会是因为这个警告吧?
      502eeef3-d4b5-40d8-bb39-e9fa15299eba-dc8a5ccdc03793d837070a5fdbcebea.png

      发布在 Fluent
      冠
      冠竹
    • fluent中在动量方程中添加弹性项以及添加5个用户自定义标量函数UDS后

      我想在fluent中使用黏弹性流体(FENE-P模型),需要使用udf在动量方程中添加弹性项源项,并且使用uds添加6个标量方程(FENE-P模型的本构方程Cxx,Cxy,Cxz,Cyy,Cyz,Czz),与N-S方程耦合求解,编译成功后可以计算,但是uds的值变化很小(初始化时uds的值给的是1,计算结束后基本还是1,知识小数点后几位有变化),几乎体现不出弹性项的作用。求助各位大佬
      (1)会不会是多个uds方程需要对i进行定义
      在帮助手册找到这个(如果要求解多个标量,则可以在UDF中使用条件IF语句为每个i定义不同的通量函数。i=0与uds-0(要求解的第一个标量方程)相关联。)但是没找到具体怎么操作
      FENE-P 模型的本构方程是
      1db2aecf-bb28-4e18-83b2-b78012d300a0-8de57d7561f6f0ce3fa3533604fefd3.png
      直角坐标下为
      2c357f13-0b20-4832-8817-2d63da569bd0-4ad7186b8fd9973a10024ee26d56f59.png
      49efe7aa-242a-40a1-99f0-77492f023026-e3f29d15bfaa922b3ad458e6135be3c.png
      2f7eef5f-0dc2-481c-ae2e-fbf59aaf90f4-2675ed868f2111033ffc1d7af482a07.png
      添加弹性项后的X方向上的动量方程
      7da0ce9c-10e5-467d-bbc7-a5e298f58678-8985c92d1edb7e6c3b4912f0ccb09ae.png
      添加弹性项后的y方向上的动量方程
      3c27a4d7-a3de-4795-97bf-85f89b1ca171-7f1fc03c84026987030e60dd8546d30.png
      添加弹性项后的z方向上的动量方程
      3ba2dfa4-8db9-4ae7-ac16-e72c39f52553-4506fbbd7bf7138e779729d1836a718.png
      具体的udf代码

      #include "mem.h"
      #define miup (0.00736)
      #define lamda (2.47024)
      #define L (40)
      
      
      DEFINE_UDS_UNSTEADY(No_Cxx_unsteady,c,t,i,apu,su)        //cxx非稳态项
       {
       real physical_dt, vol, rho, phi_old;
       physical_dt = RP_Get_Real("physical-time-step");
       vol = C_VOLUME(c,t);
       
       *apu = -vol / physical_dt;/*implicit part*/
       phi_old = C_STORAGE_R(c,t,SV_UDSI_M1(0));
       *su = vol*phi_old/physical_dt;/*explicit part*/
       }
      
      
      DEFINE_UDS_UNSTEADY(No_Cxy_unsteady,c,t,i,apu,su)         //cxy非稳态项
       {
       real physical_dt, vol, rho, phi_old;
       physical_dt = RP_Get_Real("physical-time-step");
       vol = C_VOLUME(c,t);
       
       *apu = -vol / physical_dt;/*implicit part*/
       phi_old = C_STORAGE_R(c,t,SV_UDSI_M1(1));
       *su = vol*phi_old/physical_dt;/*explicit part*/
       }
      
       DEFINE_UDS_UNSTEADY(No_Cxz_unsteady,c,t,i,apu,su)       //cxz非稳态项
       {
       real physical_dt, vol, rho, phi_old;
       physical_dt = RP_Get_Real("physical-time-step");
       vol = C_VOLUME(c,t);
       
       *apu = -vol / physical_dt;/*implicit part*/
       phi_old = C_STORAGE_R(c,t,SV_UDSI_M1(2));
       *su = vol*phi_old/physical_dt;/*explicit part*/
       }
      
       
      
       DEFINE_UDS_UNSTEADY(No_Cyy_unsteady,c,t,i,apu,su)       //cyy非稳态项
       {
       real physical_dt, vol, rho, phi_old;
       physical_dt = RP_Get_Real("physical-time-step");
       vol = C_VOLUME(c,t);
       
       *apu = -vol / physical_dt;/*implicit part*/
       phi_old = C_STORAGE_R(c,t,SV_UDSI_M1(3));
       *su = vol*phi_old/physical_dt;/*explicit part*/
       }
      
       DEFINE_UDS_UNSTEADY(No_Cyz_unsteady,c,t,i,apu,su)       //cyz非稳态项
       {
       real physical_dt, vol, rho, phi_old;
       physical_dt = RP_Get_Real("physical-time-step");
       vol = C_VOLUME(c,t);
       
       *apu = -vol / physical_dt;/*implicit part*/
       phi_old = C_STORAGE_R(c,t,SV_UDSI_M1(4));
       *su = vol*phi_old/physical_dt;/*explicit part*/
       }
      
       
      
       DEFINE_UDS_UNSTEADY(No_Czz_unsteady,c,t,i,apu,su)       //czz非稳态项
       {
           	real physical_dt, vol, rho, phi_old;
           	physical_dt = RP_Get_Real("physical-time-step");
           	vol = C_VOLUME(c,t);
       
           	*apu = -vol / physical_dt;/*implicit part*/
           	phi_old = C_STORAGE_R(c,t,SV_UDSI_M1(5));
           	*su = vol*phi_old/physical_dt;/*explicit part*/
       }
      
      DEFINE_UDS_FLUX(No_Cxx_FLUX,f,t,i)                 //Cxx对流项    
      {
          	cell_t c0, c1 = -1;
          	Thread *t0, *t1 = NULL;
          	real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;
          	c0 = F_C0(f,t);
          	t0 = F_C0_THREAD(f,t);
          	F_AREA(A, f, t);
          
         	if (BOUNDARY_FACE_THREAD_P(t)) 
          	{
             		real dens;
             
             		if (NNULLP(THREAD_STORAGE(t,SV_DENSITY)))
               		dens = F_R(f,t);  
             		else
               		dens = C_R(c0,t0); 
             		NV_DS(psi_vec, =, F_U(f,t), F_V(f,t), F_W(f,t), *, 1);
            		flux = NV_DOT(psi_vec, A); 
      	}
          	else
          	{
           		c1 = F_C1(f,t);   
           		t1 = F_C1_THREAD(f,t);
           		NV_DS(psi_vec, =, C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),*,1);
           		NV_DS(psi_vec, +=, C_U(c1,t1),C_V(c1,t1),C_W(c1,t1),*,1);
           		flux = NV_DOT(psi_vec, A)/2.0;
          	}
      	return flux;
      } 
      
      
      
      DEFINE_UDS_FLUX(No_Cxy_FLUX,f,t,i)                 //Cxy对流项    
      {
          	cell_t c0, c1 = -1;
          	Thread *t0, *t1 = NULL;
          	real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;
          	c0 = F_C0(f,t);
          	t0 = F_C0_THREAD(f,t);
          	F_AREA(A, f, t);
          
         	if (BOUNDARY_FACE_THREAD_P(t)) 
          	{
             		real dens;
             
             		if (NNULLP(THREAD_STORAGE(t,SV_DENSITY)))
               		dens = F_R(f,t);  
             		else
               		dens = C_R(c0,t0); 
             		NV_DS(psi_vec, =, F_U(f,t), F_V(f,t), F_W(f,t), *, 1);
            		flux = NV_DOT(psi_vec, A); 
      	}
          	else
          	{
           		c1 = F_C1(f,t);   
           		t1 = F_C1_THREAD(f,t);
           		NV_DS(psi_vec, =, C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),*,1);
           		NV_DS(psi_vec, +=, C_U(c1,t1),C_V(c1,t1),C_W(c1,t1),*,1);
           		flux = NV_DOT(psi_vec, A)/2.0;
          	}
      	return flux;
      } 
      
      
      DEFINE_UDS_FLUX(No_Cxz_FLUX,f,t,i)                 //Cxz对流项    
      {
          	cell_t c0, c1 = -1;
          	Thread *t0, *t1 = NULL;
          	real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;
          	c0 = F_C0(f,t);
          	t0 = F_C0_THREAD(f,t);
          	F_AREA(A, f, t);
          
         	if (BOUNDARY_FACE_THREAD_P(t)) 
          	{
             		real dens;
             
             		if (NNULLP(THREAD_STORAGE(t,SV_DENSITY)))
               		dens = F_R(f,t);  
             		else
               		dens = C_R(c0,t0); 
             		NV_DS(psi_vec, =, F_U(f,t), F_V(f,t), F_W(f,t), *, 1);
            		flux = NV_DOT(psi_vec, A); 
      	}
          	else
          	{
           		c1 = F_C1(f,t);   
           		t1 = F_C1_THREAD(f,t);
           		NV_DS(psi_vec, =, C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),*,1);
           		NV_DS(psi_vec, +=, C_U(c1,t1),C_V(c1,t1),C_W(c1,t1),*,1);
           		flux = NV_DOT(psi_vec, A)/2.0;
          	}
      	return flux;
      } 
      
      
      
      DEFINE_UDS_FLUX(No_Cyy_FLUX,f,t,i)                 //Cyy对流项    
      {
          	cell_t c0, c1 = -1;
          	Thread *t0, *t1 = NULL;
          	real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;
          	c0 = F_C0(f,t);
          	t0 = F_C0_THREAD(f,t);
          	F_AREA(A, f, t);
          
         	if (BOUNDARY_FACE_THREAD_P(t)) 
          	{
             		real dens;
             
             		if (NNULLP(THREAD_STORAGE(t,SV_DENSITY)))
               		dens = F_R(f,t);  
             		else
               		dens = C_R(c0,t0); 
             		NV_DS(psi_vec, =, F_U(f,t), F_V(f,t), F_W(f,t), *, 1);
            		flux = NV_DOT(psi_vec, A); 
      	}
          	else
          	{
           		c1 = F_C1(f,t);   
           		t1 = F_C1_THREAD(f,t);
           		NV_DS(psi_vec, =, C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),*,1);
           		NV_DS(psi_vec, +=, C_U(c1,t1),C_V(c1,t1),C_W(c1,t1),*,1);
           		flux = NV_DOT(psi_vec, A)/2.0;
          	}
      	return flux;
      } 
      
      DEFINE_UDS_FLUX(No_Cyz_FLUX,f,t,i)                 //Cyz对流项    
      {
          	cell_t c0, c1 = -1;
          	Thread *t0, *t1 = NULL;
          	real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;
          	c0 = F_C0(f,t);
          	t0 = F_C0_THREAD(f,t);
          	F_AREA(A, f, t);
          
         	if (BOUNDARY_FACE_THREAD_P(t)) 
          	{
             		real dens;
             
             		if (NNULLP(THREAD_STORAGE(t,SV_DENSITY)))
               		dens = F_R(f,t);  
             		else
               		dens = C_R(c0,t0); 
             		NV_DS(psi_vec, =, F_U(f,t), F_V(f,t), F_W(f,t), *, 1);
            		flux = NV_DOT(psi_vec, A); 
      	}
          	else
          	{
           		c1 = F_C1(f,t);   
           		t1 = F_C1_THREAD(f,t);
           		NV_DS(psi_vec, =, C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),*,1);
           		NV_DS(psi_vec, +=, C_U(c1,t1),C_V(c1,t1),C_W(c1,t1),*,1);
           		flux = NV_DOT(psi_vec, A)/2.0;
          	}
      	return flux;
      } 
      
      
      
      DEFINE_UDS_FLUX(No_Czz_FLUX,f,t,i)                 //Czz对流项    
      {
          	cell_t c0, c1 = -1;
          	Thread *t0, *t1 = NULL;
          	real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;
          	c0 = F_C0(f,t);
          	t0 = F_C0_THREAD(f,t);
          	F_AREA(A, f, t);
          
         	if (BOUNDARY_FACE_THREAD_P(t)) 
          	{
             		real dens;
             
             		if (NNULLP(THREAD_STORAGE(t,SV_DENSITY)))
               		dens = F_R(f,t);  
             		else
               		dens = C_R(c0,t0); 
             		NV_DS(psi_vec, =, F_U(f,t), F_V(f,t), F_W(f,t), *, 1);
            		flux = NV_DOT(psi_vec, A); 
      	}
          	else
          	{
           		c1 = F_C1(f,t);   
           		t1 = F_C1_THREAD(f,t);
           		NV_DS(psi_vec, =, C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),*,1);
           		NV_DS(psi_vec, +=, C_U(c1,t1),C_V(c1,t1),C_W(c1,t1),*,1);
           		flux = NV_DOT(psi_vec, A)/2.0;
          	}
      	return flux;
      } 
      
      DEFINE_SOURCE(No_Cxx_source, c, t, dS, eqn)             //Cxx源项
      {  
      	real source;
          real fR;
          fR=(L*L-3)/(L*L-C_UDSI(c,t,0)-C_UDSI(c,t,3)-C_UDSI(c,t,5));
      	source=2*(C_UDSI(c,t,0)*C_DUDX(c,t)+C_UDSI(c,t,1)*C_DUDY(c,t)+C_UDSI(c,t,2)*C_DUDZ(c,t))-(fR*C_UDSI(c,t,0)-1)/lamda;
      	dS[eqn]=0.0;
      
      	return source;
      }
      
      DEFINE_SOURCE(No_Cxy_source, c, t, dS, eqn)             //Cxy源项
      {  
      	real source;
          real fR;
          fR=(L*L-3)/(L*L-C_UDSI(c,t,0)-C_UDSI(c,t,3)-C_UDSI(c,t,5));
      	source=(C_UDSI(c,t,0)*C_DVDX(c,t)+C_UDSI(c,t,1)*C_DVDY(c,t)+C_UDSI(c,t,2)*C_DVDZ(c,t)+\
          C_UDSI(c,t,1)*C_DUDX(c,t)+C_UDSI(c,t,3)*C_DUDY(c,t)+C_UDSI(c,t,4)*C_DUDZ(c,t))-\
          (fR*C_UDSI(c,t,1))/lamda;
      	dS[eqn]=0.0;
      
      	return source;
      }
      
      
      DEFINE_SOURCE(No_Cxz_source, c, t, dS, eqn)             //Cxz源项
      {  
      	real source;
          real fR;
          fR=(L*L-3)/(L*L-C_UDSI(c,t,0)-C_UDSI(c,t,3)-C_UDSI(c,t,5));
      	source=(C_UDSI(c,t,0)*C_DWDX(c,t)+C_UDSI(c,t,1)*C_DWDY(c,t)+C_UDSI(c,t,2)*C_DWDZ(c,t)+\
          C_UDSI(c,t,2)*C_DUDX(c,t)+C_UDSI(c,t,4)*C_DUDY(c,t)+C_UDSI(c,t,5)*C_DUDZ(c,t))-\
          (fR*C_UDSI(c,t,2))/lamda;
      	dS[eqn]=0.0;
      
      	return source;
      }
      
      
      
      
      DEFINE_SOURCE(No_Cyy_source, c, t, dS, eqn)             //Cyy源项
      {  
      	real source;
          real fR;
          fR=(L*L-3)/(L*L-C_UDSI(c,t,0)-C_UDSI(c,t,3)-C_UDSI(c,t,5));
      	source=2*(C_UDSI(c,t,1)*C_DVDX(c,t)+C_UDSI(c,t,3)*C_DVDY(c,t)+C_UDSI(c,t,4)*C_DVDZ(c,t))-(fR*C_UDSI(c,t,3)-1)/lamda;
      	dS[eqn]=0.0;
      
      	return source;
      }
      
      
      DEFINE_SOURCE(No_Cyz_source, c, t, dS, eqn)             //Cyz源项
      {  
      	real source;
          real fR;
          fR=(L*L-3)/(L*L-C_UDSI(c,t,0)-C_UDSI(c,t,3)-C_UDSI(c,t,5));
      	source=(C_UDSI(c,t,1)*C_DWDX(c,t)+C_UDSI(c,t,3)*C_DWDY(c,t)+C_UDSI(c,t,4)*C_DWDZ(c,t)+\
          C_UDSI(c,t,2)*C_DVDX(c,t)+C_UDSI(c,t,4)*C_DVDY(c,t)+C_UDSI(c,t,5)*C_DVDZ(c,t))-\
          (fR*C_UDSI(c,t,4))/lamda;
      	dS[eqn]=0.0;
      
      	return source;
      }
      
      
      
      
      
      
      
      DEFINE_SOURCE(No_Czz_source, c, t, dS, eqn)             //Czz源项
      {  
      	real source;
          real fR;
          fR=(L*L-3)/(L*L-C_UDSI(c,t,0)-C_UDSI(c,t,3)-C_UDSI(c,t,5));
      	source=2*(C_UDSI(c,t,2)*C_DWDX(c,t)+C_UDSI(c,t,4)*C_DWDY(c,t)+C_UDSI(c,t,5)*C_DWDZ(c,t))-(fR*C_UDSI(c,t,5)-1)/lamda;
      	dS[eqn]=0.0;
      
      	return source;
      }
      
      
      DEFINE_SOURCE(x_momentum, c, t, dS, eqn)             //x方向动量源项
      {  
      	real source;
          real cxx = C_UDSI(c,t,0);
          real cxy = C_UDSI(c,t,1);
          real cxz = C_UDSI(c,t,2);
          real cyx = C_UDSI(c,t,1);
          real cyy = C_UDSI(c,t,3);
          real cyz = C_UDSI(c,t,4);
          real czx = C_UDSI(c,t,2);
          real czy = C_UDSI(c,t,4);
          real czz = C_UDSI(c,t,5);
          real NV_VEC(gcxx);
          real NV_VEC(gcxy);
          real NV_VEC(gcxz);
          real NV_VEC(gcyx);
          real NV_VEC(gcyy);
          real NV_VEC(gcyz);
          real NV_VEC(gczx);
          real NV_VEC(gczy);
          real NV_VEC(gczz);
          NV_DS(gcxx, =,C_UDSI_G(c,t,0)[0],C_UDSI_G(c,t,0)[1],C_UDSI_G(c,t,0)[2],*,1);
          NV_DS(gcxy, =,C_UDSI_G(c,t,1)[0],C_UDSI_G(c,t,1)[1],C_UDSI_G(c,t,1)[2],*,1);
          NV_DS(gcxz, =,C_UDSI_G(c,t,2)[0],C_UDSI_G(c,t,2)[1],C_UDSI_G(c,t,2)[2],*,1);
          NV_DS(gcyx, =,C_UDSI_G(c,t,1)[0],C_UDSI_G(c,t,1)[1],C_UDSI_G(c,t,1)[2],*,1);
          NV_DS(gcyy, =,C_UDSI_G(c,t,3)[0],C_UDSI_G(c,t,3)[1],C_UDSI_G(c,t,3)[2],*,1);
          NV_DS(gcyz, =,C_UDSI_G(c,t,4)[0],C_UDSI_G(c,t,4)[1],C_UDSI_G(c,t,4)[2],*,1);
          NV_DS(gczx, =,C_UDSI_G(c,t,2)[0],C_UDSI_G(c,t,2)[1],C_UDSI_G(c,t,2)[2],*,1);
          NV_DS(gczy, =,C_UDSI_G(c,t,4)[0],C_UDSI_G(c,t,4)[1],C_UDSI_G(c,t,4)[2],*,1);
          NV_DS(gczz, =,C_UDSI_G(c,t,5)[0],C_UDSI_G(c,t,5)[1],C_UDSI_G(c,t,5)[2],*,1);
          source = miup*(L*L-3)/lamda/pow(L*L-cxx-cyy-czz,2)*((gcxx[0]+gcxy[1]+gcxz[2])*(L*L-cxx-cyy-czz)+\
          cxx*(gcxx[0]+gcyy[0]+gczz[0])+cxy*(gcxx[1]+gcyy[1]+gczz[1])+cxz*(gcxx[2]+gcyy[2]+gczz[2]));
      	return source;
      }
      
      
      DEFINE_SOURCE(y_momentum, c, t, dS, eqn)             //y方向动量源项
      {  
      	real source;
          real cxx = C_UDSI(c,t,0);
          real cxy = C_UDSI(c,t,1);
          real cxz = C_UDSI(c,t,2);
          real cyx = C_UDSI(c,t,1);
          real cyy = C_UDSI(c,t,3);
          real cyz = C_UDSI(c,t,4);
          real czx = C_UDSI(c,t,2);
          real czy = C_UDSI(c,t,4);
          real czz = C_UDSI(c,t,5);
          real NV_VEC(gcxx);
          real NV_VEC(gcxy);
          real NV_VEC(gcxz);
          real NV_VEC(gcyx);
          real NV_VEC(gcyy);
          real NV_VEC(gcyz);
          real NV_VEC(gczx);
          real NV_VEC(gczy);
          real NV_VEC(gczz);
          NV_DS(gcxx, =,C_UDSI_G(c,t,0)[0],C_UDSI_G(c,t,0)[1],C_UDSI_G(c,t,0)[2],*,1);
          NV_DS(gcxy, =,C_UDSI_G(c,t,1)[0],C_UDSI_G(c,t,1)[1],C_UDSI_G(c,t,1)[2],*,1);
          NV_DS(gcxz, =,C_UDSI_G(c,t,2)[0],C_UDSI_G(c,t,2)[1],C_UDSI_G(c,t,2)[2],*,1);
          NV_DS(gcyx, =,C_UDSI_G(c,t,1)[0],C_UDSI_G(c,t,1)[1],C_UDSI_G(c,t,1)[2],*,1);
          NV_DS(gcyy, =,C_UDSI_G(c,t,3)[0],C_UDSI_G(c,t,3)[1],C_UDSI_G(c,t,3)[2],*,1);
          NV_DS(gcyz, =,C_UDSI_G(c,t,4)[0],C_UDSI_G(c,t,4)[1],C_UDSI_G(c,t,4)[2],*,1);
          NV_DS(gczx, =,C_UDSI_G(c,t,2)[0],C_UDSI_G(c,t,2)[1],C_UDSI_G(c,t,2)[2],*,1);
          NV_DS(gczy, =,C_UDSI_G(c,t,4)[0],C_UDSI_G(c,t,4)[1],C_UDSI_G(c,t,4)[2],*,1);
          NV_DS(gczz, =,C_UDSI_G(c,t,5)[0],C_UDSI_G(c,t,5)[1],C_UDSI_G(c,t,5)[2],*,1);
          source = miup*(L*L-3)/lamda/pow(L*L-cxx-cyy-czz,2)*((gcyx[0]+gcyy[1]+gcyz[2])*(L*L-cxx-cyy-czz)+\
          cyx*(gcxx[0]+gcyy[0]+gczz[0])+cyy*(gcxx[1]+gcyy[1]+gczz[1])+cyz*(gcxx[2]+gcyy[2]+gczz[2]));
      	return source;
      }
      
      DEFINE_SOURCE(z_momentum, c, t, dS, eqn)             //z方向动量源项
      {  
      	real source;
          real cxx = C_UDSI(c,t,0);
          real cxy = C_UDSI(c,t,1);
          real cxz = C_UDSI(c,t,2);
          real cyx = C_UDSI(c,t,1);
          real cyy = C_UDSI(c,t,3);
          real cyz = C_UDSI(c,t,4);
          real czx = C_UDSI(c,t,2);
          real czy = C_UDSI(c,t,4);
          real czz = C_UDSI(c,t,5);
          real NV_VEC(gcxx);
          real NV_VEC(gcxy);
          real NV_VEC(gcxz);
          real NV_VEC(gcyx);
          real NV_VEC(gcyy);
          real NV_VEC(gcyz);
          real NV_VEC(gczx);
          real NV_VEC(gczy);
          real NV_VEC(gczz);
          NV_DS(gcxx, =,C_UDSI_G(c,t,0)[0],C_UDSI_G(c,t,0)[1],C_UDSI_G(c,t,0)[2],*,1);
          NV_DS(gcxy, =,C_UDSI_G(c,t,1)[0],C_UDSI_G(c,t,1)[1],C_UDSI_G(c,t,1)[2],*,1);
          NV_DS(gcxz, =,C_UDSI_G(c,t,2)[0],C_UDSI_G(c,t,2)[1],C_UDSI_G(c,t,2)[2],*,1);
          NV_DS(gcyx, =,C_UDSI_G(c,t,1)[0],C_UDSI_G(c,t,1)[1],C_UDSI_G(c,t,1)[2],*,1);
          NV_DS(gcyy, =,C_UDSI_G(c,t,3)[0],C_UDSI_G(c,t,3)[1],C_UDSI_G(c,t,3)[2],*,1);
          NV_DS(gcyz, =,C_UDSI_G(c,t,4)[0],C_UDSI_G(c,t,4)[1],C_UDSI_G(c,t,4)[2],*,1);
          NV_DS(gczx, =,C_UDSI_G(c,t,2)[0],C_UDSI_G(c,t,2)[1],C_UDSI_G(c,t,2)[2],*,1);
          NV_DS(gczy, =,C_UDSI_G(c,t,4)[0],C_UDSI_G(c,t,4)[1],C_UDSI_G(c,t,4)[2],*,1);
          NV_DS(gczz, =,C_UDSI_G(c,t,5)[0],C_UDSI_G(c,t,5)[1],C_UDSI_G(c,t,5)[2],*,1);
          source = miup*(L*L-3)/lamda/pow(L*L-cxx-cyy-czz,2)*((gczx[0]+gczy[1]+gczz[2])*(L*L-cxx-cyy-czz)+\
          czx*(gcxx[0]+gcyy[0]+gczz[0])+czy*(gcxx[1]+gcyy[1]+gczz[1])+czz*(gcxx[2]+gcyy[2]+gczz[2]));
      	return source;
      }
      code_text
      
      发布在 Fluent
      冠
      冠竹
    • OpenFOAMv2206中dictionaryConstructorTable程序问题

      各位大佬,我想在OpenFOAMv2206中植入我在OpenFOAM9中修改的模型。但是在编译过程中出现了错误:

      9cb8d1c3-cfce-407c-a09c-cdecbf702dd3-图片.png

      我的代码在OF9中是可以运行的,但是在OpenFOAMv2206中就报错了。看其他的代码感觉OpenFOAM9和OpenFOAMv2206没有太大的区别,这个错误是为什么呢?

      Foam::autoPtr<Foam::viscoelasticLogLaw> Foam::viscoelasticLogLaw::New
      (
          const word& name,
          const volVectorField& U,
          const surfaceScalarField& phi,
          const dictionary& dict
      )
      {
          const word typeName(dict.get<word>("type"));
      
          Info<< "Selecting viscoelasticLog model " << typeName << endl;
      
          dictionaryConstructorTable::iterator cstrIter =
              dictionaryConstructorTablePtr_->find(typeName);
      
          if (cstrIter == dictionaryConstructorTablePtr_->end())
          {
              FatalErrorIn
              (
                  "viscoelasticLogLaw::New(const word& name, const volVectorField&, "
                  "const surfaceScalarField&)"
              )   << "Unknown viscoelasticLogLaw type " << typeName
                  << endl << endl
                  << "Valid viscoelasticLogLaw types are :" << endl
                  << dictionaryConstructorTablePtr_->sortedToc()
                  << exit(FatalError);
          }
      
          return autoPtr<viscoelasticLogLaw>(cstrIter()(name, U, phi, dict));
      }
      
      发布在 OpenFOAM
      冠
      冠竹
    • RE: 槽道流Umean设置疑惑

      @vbcwl 老哥我这么理解对不对,就是速度设置成U+,也就是U/Ut,那这个按照自带的算例Ret=395的话,是不是就是U=Ubar/ut=0.1335/0.079呢

      发布在 OpenFOAM
      冠
      冠竹
    • 槽道流Umean设置疑惑

      各位大佬,我想复现一下Ret=395的一个槽道流,但是有一个疑惑的点是平均速度应该如何设置呢?Ret和传统的Re之间有换算关系么

      发布在 OpenFOAM
      冠
      冠竹
    • DNS求取yplus

      各位大佬,在DNS槽道流case中想求取yPlus进行数据分析,但OpenFOAM里边的yPlus工具只适用于LES和Rans,有适用于DNS的yPlus工具吗

      发布在 OpenFOAM
      冠
      冠竹
    • RE: 在编译求解器时引用弹性项

      @guangzhou 这个不影响吧,我用rheoTool的目的是他的粘弹性流体库。他的求解器对我来说限制有点多,所以我就直接把他当库文件用自己新写了一个DNS求解器。如果觉得他的不适合的话就引用库自己改一个就好了啊

      发布在 OpenFOAM
      冠
      冠竹
    • RE: 在编译求解器时引用弹性项

      @guangzhou 这个求解器在高维森贝格数下会产生稳定性问题,所以一般采用对数重构法去解决这个问题。论文可以参考张红娜老师在17年的一篇论文《Numerical Simulation of Heat Transfer Process of Viscoelastic Fluid Flow at High Weissenberg Number by Log-Conformation Reformulation》,或者可以看看rheoTool工具箱相关介绍

      发布在 OpenFOAM
      冠
      冠竹