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
      冠
      冠竹