Fluent UDF中根据zone的名字获取ID

Fluent UDF中根据zone的名字获取ID

Fluent UDF中经常用到thread*类型的指针,一般可以通过如下函数来获取。

Lookup_Thread(Domain*domain, int id)

其中id是边界面的ID或者cell区域的ID,如下图中outlet边界的ID是2。

 

 

 

 

 

 

 

 

 

 

 

遗憾的是,网格载入Fluent后,其ID都是无法事先确定的。很多朋友只有在UDF源码开头用如下代码根据事后查到的ID手动定义,每次网格一变化又得重新在源码里面修改,重新编译,十分影响通用性。

#define OUTLET_ID 2 //每次不同网格需要根据情况修改后重新编译
Domain* domain=Get_Domain(1); //非多相流或多相流混合物的domain永远是1
Thread*tf=Lookup_Thread(domain, OUTLET_ID);

有没有一种办法能够一次性源码编译达到适应于所有网格呢?答案当然是肯定的,你可以通过zone名字来获取其ID号,然后画网格的时候只要取相同名字即可。实现该功能的函数源代码如下(插件VC++ UDF Studio 2022R1学术版上编译通过)

#include "udf.h"
#include "SuperUdfExtension.h" //VC++ UDF Studio自带的扩展库头文件,具体参考该软件中的编程手册
#pragma comment(lib, "SuperUdfExtension.lib")  //VC++ UDF Studio自带的扩展库的lib文件 
int GetZoneIdByName(CString zoneName)  //适用于所有Fluent版本
{
    int returnID=-1;
    Domain*domain=Get_Domain(1);
    CString strCurrentFluentVersion;
    strCurrentFluentVersion.Format("%d.%d", RampantReleaseMajor, RampantReleaseMinor);//格式化当前Fluent版本为字符串形式
    double fCurrentFluentVersion = atof(strCurrentFluentVersion.GetBuffer()); //当前Fluent版本转为double类型

    if(fCurrentFluentVersion<=19.2)  // 对于Fluent6.3-19.2,只能调用VC++ UDF Studio扩展库
    {
        SuperUdf_Initialize(AfxGetInstanceHandle()); //调用VC++ UDF Studio扩展库中任何函数之前必须调用此初始化函数,具体参考该软件中的编程手册
#if !RP_NODE
        returnID=SuperUdf_GetZoneIdByName(zoneName.GetBuffer()); //调用VC++ UDF Studio扩展库中的SuperUdf_GetZoneIdByName函数,具体参考该软件中的编程手册
#endif
        host_to_node_int_1(returnID);
    }
    else // 对于Fluent version >=19.3,有直接UDF函数可以实现
    {        
        Thread*tf;
        thread_loop_f(tf, domain) //对所有面的thread进行循环查找
        {
            if(0==zoneName.CompareNoCase(THREAD_NAME(tf))) //对比名字是否相同
            {
                returnID=THREAD_ID(tf);
                break;
            }
        } 
        if(-1==returnID)  //如果面的thread中无法找到匹配名字
        {
            Thread*tc;
            thread_loop_c(tc, domain) //对所有网格的thread进行循环查找
            {
                if(0==zoneName.CompareNoCase(THREAD_NAME(tc))) //对比名字是否相同
                {
                    returnID=THREAD_ID(tc);
                    break;
                }
            }
        }
    }

    return returnID;
}


DEFINE_EXECUTE_ON_LOADING(get_id, libudf)
{
    int theID=GetZoneIdByName("inlet"); //根据边界名字获取其ID,如果返回-1表示找不到
    Message("the zone id of inlet is %d\n",theID);
}

以上源代码实现了通过zone名字来获取其ID号的功能。对于Fluent19.3或更高版本,可以利用THREAD_NAME与需要的名字进行对比,匹配情况下用THREAD_ID获得其ID。但对于Fluent19.2或更低版本, THREAD_NAME不起作用(可能Fluent的bug),我们只能依赖于插件VC++ UDF Studio (https://vcudfstudio.github.io)中的拓展库函数SuperUdf_GetZoneIdByName来实现。是时候抛弃傻傻的#define ID了。

 

posted @ 2022-02-12 15:14  SuperUDF  阅读(1706)  评论(0)    收藏  举报