找回密码
 注册
查看: 256|回复: 0

代码第三行有错误,有人能看出来吗?

[复制链接]
发表于 2023-12-23 21:42:26 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x


#include "udf,h"
#define PI 3.1415926535
DEFINE_DPM_LAW(xie,tp,ci)
{
  Tracked_Particle *tp;/*定位粒子*/
  int ci=1;/*连续相与离散相耦合*/
  Cell_t c=TP_CELL(tp);/*粒子所在的单元格*/
  Thread *t=TP_CELL_THREAD(tp);/*定义单元格所在链表*/
  real kr;/*扩散系数*/
  real kd;/*动力系数*/
  real C1=5.0e-12.0;/*与扩散有关的经验常数*/
  real tpa=TP_T(tp);/*粒子当前的温度*/
  real tc=C_T(c,t);/*单元格的温度*/
  real dp=TP_DIAM(tp);/*定义粒子的直径*/
  kd=C1*pow(((tpa+ts)/2),0.75)/dp;/*计算扩散系数*/
  kr=3.3e+19.0*exp(-167400.0/(UNIVERSAL_GAS_CONSTANT*tpa));/*计算动力系数*/
  real pox;/*定义氧气的分压*/
  pox=C_R(c,t)*UNIVERSAL_GAS_CONSTANT*tpa*TP_OXID_SPECIES_INDEX(tp)/32.0;/*计算氧气的环境分压*/
  real dm,dt;/*定义粒子的质量变化率与时间步长*/
  dt=TP_DT(tp);
  dm=-PI*dp*dp*pox*kr*kd/(kd+kr);/*焦炭因为燃烧消耗的质量*/
  TP_MASS(tp)=TP_MASS0(tp)+dm*(TP_TIME(tp)-TP_TIME0(tp));/*当前的粒子质量等于初始粒子的质量加上质量变化率与时间的乘积*/
  C_DPMS_BURNOUT(c,t)=dm;

显示第三行有问题,不知道哪里错了。
我正在做有关焦炭燃烧的模拟。请问有没有人有fluent动力扩散模型的代码呢



您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表