质量源项是压力关于时间的偏导

jwj

如果质量源项是压力关于时间的偏导该怎么表示，我写了一个编译没错但一是运行FLUENT就没了这是怎么回事，请各位高手指点指点
#include "udf.h"
#define a 24.114
#define b 2.62
DEFINE_SOURCE(my_source,c,t,ds,eqn)
{
                real source;
        real time=RP_Get_Real("physical-time-step");
        real time_step=RP_Get_Integer("time-step");
    real P=C_P(c,t)/1000000;     /*兆帕*/
        real P1=C_P_M1(c,t)/1000000; /*兆帕*/
    real D=C_R(c,t);
        if(time_step>1)
         source=(a*b/((1+b*P)*(1+b*P))*(P1-P)*D/time)*1.39;  /* 公式源项单位 */
        return source;
}

swimming02

#include "udf.h"
#define a 24.114
#define b 2.62
DEFINE_SOURCE(my_source,c,t,dS,eqn)
{
real source;
real time=RP_Get_Real("physical-time-step");
real P=C_P(c,t)/1000000;     /*兆帕*/
real P1=C_P_M1(c,t)/1000000; /*兆帕*/
real D=C_R(c,t);
int time_step=RP_Get_Integer("time-step");

if(time_step>1)
source=(a*b/((1+b*P)*(1+b*P))*(P1-P)*D/time)*1.39;  /* 公式源项单位 */
dS[eqn]=0.0;
return source;
}

cumtsafe

学习了，谢谢

hdlishuhua

学习了 谢谢 我想问一下 你这个源项从FLUENT中的哪个口导入/

易阳

学习了，但有个问题，我用k = F_U(f,thread);
        ki = F_U_M1(f,thread);
    dq = 0.0125*(k-ki)/0.00005;
来表示，速度对时间的导数，，为什么有问题呢？，，编译可以通过，，但不能初始化
#include "udf.h"

#define H 0.625
#define W 0.221
#define fic0 0.034
#define R 0.326
#define le 0.63
#define li 0.63
#define Ud 102
#define mu 0.00002
#define a 3.2
#define m 1.75
#define Vp 2.5
#define as 340
#define Lc 1.11
#define Den 1
#define Ac 0.25
#define I 0.03
#define c1 0.02
#define c2 0.017

DEFINE_PROFILE(unsteady_pressure, thread, position)
{
  face_t f;
  real p,k,ki;
  real q,dq,dj,b,lc,mB,ent,La1,u1,u2,Z,X;

  begin_f_loop(f, thread)
  {
    real t = RP_Get_Real("flow-time");
    p = F_P(f,thread);
    k = F_U(f,thread);
        ki = F_U_M1(f,thread);
   
    q = 0.0125*k;
    b = 2*as*(sqrt(Ac*Lc/Vp));
    lc = li+le+1/a;
    mB = (1-m)*b-1;
    ent = -q*q;
    La1 = (Den*R*R*R*Ac*b)/(I*Ud);
    dq = 0.0125*(k-ki)/0.00005;
        Z = lc*dq/H+(p-fic0)/H-1+0.5*(q/W-1)*(q/W-1)*(q/W-1)+(le*Ud*La1*ent*q)/(b*H);
        dj = (1-(q/W-1)*(q/W-1)-1/2+Z/(2*q/W-2)-(mu*W)/(3*a*H)-(2*Ud*La1*(m-1)*W*ent)/(3*b*H))*((3*a*H)/((1-mB*a)*W));
    u1 = lc*dq+p+c1*(q-2*W);
        u2 = H*(Z*dq/(q-W)-1.5*(q/W-1-Z)*dj-1.5*(q/W-1)*(q/W-1)*dq/W+(3*le*Ud*La1*q*q*dq)/(b*H))+c1*dq;   
        X = -c2*(p-u1)+(q-2*W)/lc+u2;
        F_PROFILE(f, thread, position) = 0.00005*X+p;
  }
  end_f_loop(f, thread)

}

lushenghui79

楼主，问题解决了没有，我也遇到类似问题，求指导。。。