24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2015  |  回复: 1

DustinA

金虫 (初入文坛)

[交流] fluent中的编译型udf怎么生成?怎么用?为什么直接compile会出现不可思议的错误?

将udf直接compile,报错,如下所示fluent中的编译型udf怎么生成?怎么用?为什么直接compile会出现不可思议的错误?
我的udf源码如下:

# include "udf.h"
#include "stdio.h"
#include "stdlib.h"
#include "math.h"
#include "string.h"

#define M 50
#define K 70
#define PULSE 1
#define X_MIN 0
#define Y_MIN 0
#define Y_MAX 1.4
#define VESSEL 0.06
#define DX 0.02
#define DY 0.02
#define NN 10 //矩阵分块后每块的维数

DEFINE_SOURCE(heat_source,cell,thread,dS,eqn)
{
        real time;
        real source0[NN][NN];
        real source;
        int i,k;
        int ln;
        int ii,kk,temp_ii,temp_kk;
        int file;
        char str[25]={'\0'};
        char filename[25]={'\0'};
        real x[ND_ND];
       
        time=RP_Get_Real("flow-time";
        if(time<=PULSE)
        {
                C_CENTROID(x,cell,thread);
                ii=floor((x[0]-X_MIN)/DX);
                kk=floor((x[1]-Y_MIN)/DY);//找到当前cell的编号
                temp_ii=floor(ii/10.0)+1;
                temp_kk=floor(kk/10.0)+1;
                //此时需要读取第x个文件,x=(temp_ii-1)*int(M/10)+temp_kk
                file=(temp_kk-1)*(M/NN)+temp_ii;
                //convert the int to char
                for(ln=0,i=file;i;i/=10,ln++);
                for(ln--,i=file;i;str[ln--]=i%10+'0',i/=10);
                strcpy(filename,"S2D_20(";
                strcat(filename,str);
                strcat(filename,".txt";
                /*读入热源文件*/
                FILE *fd;
                fd = fopen(filename,"r";
                for(k=0;k<NN;k++)
                        for(i=0;i<NN;i++)
                                fscanf(fd, "%lf", &source0[k]);
                fclose(fd);/*此时,将外部数据读入到数组source0中*/
                ii=ii%NN; kk=kk%NN;
                source=source0[kk][ii];//当前cell的热源等于对应编号下的外部数据
                dS[eqn]=0;//热源的导数为0
        }
        else
                source=0.0;
        return source;
}

请大侠指点啊!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

DustinA

金虫 (初入文坛)

沙发自己坐,坐等大神指点!
2楼2015-02-02 16:33:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 DustinA 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见