FLUENT UDF1：当边界上的单元达到判断条件时，设置该单元上的入口大小

本案例，当质量流率入口上的单元达到温度时，设置该单元的入口大小为12。也可以设置其它判断条件。

运行代码前，一定要先打开UDM的存储单元。否则会报错。

#include "udf.h"
#define IGNITION_TEMPERATURE 2000.0     // 点火温度，单位：K
#define HIGH_TEMP 3693.0               // 高温空气温度，单位：K

// 定义按需函数
DEFINE_ON_DEMAND(set_udmx)
{
    static int thread_IDs[] = { 22 };  // 定义多个线程ID
    int num_threads = sizeof(thread_IDs) / sizeof(thread_IDs[0]);
    Domain* d = Get_Domain(1);  // 获取域指针
    Thread* tw;
    Thread* t0;  // 相邻单元线程
    face_t face;  // 面
    cell_t c0;  // 单元格

    // 初始化所有单元的用户定义内存量为0.0
    thread_loop_c(t0, d)
    {
        begin_c_loop_all(c0, t0)
        {
            C_UDMI(c0, t0, 0) = 0.0;
        }
        end_c_loop_all(c0, t0)
    }

    // 遍历所有指定的线程ID
    for (int i = 0; i < num_threads; i++)
    {
        tw = Lookup_Thread(d, thread_IDs[i]);  // 查找特定线程

        // 将特定面相邻单元的用户定义内存量设置为2.0
        begin_f_loop(face, tw)
        {
            t0 = THREAD_T0(tw);  // 获取面对应的相邻单元线程
            c0 = F_C0(face, tw);  // 获取面对应的相邻单元
            C_UDMI(c0, t0, 0) = 2.0;  // 将相邻单元的用户定义内存量设置为2.0
        }
        end_f_loop(face, tw)
    }
}

// 定义调节函数以更新温度并注入气体
DEFINE_ADJUST(update_temperature, domain)
{
    Thread* t;
    cell_t c;

    // 遍历域中的所有单元线程
    thread_loop_c(t, domain)
    {
        // 遍历单元线程中的所有单元
        begin_c_loop_all(c, t)
        {
            // 检查 C_UDMI(c, t, 0) 是否为 2.0
            if (C_UDMI(c, t, 0) == 2.0)
            {
                real temp = C_T(c, t);

                // 检查单元温度是否超过点火条件
                if (temp >= IGNITION_TEMPERATURE)
                {
                    // 将单元温度设置为高温气体温度
                    C_T(c, t) = HIGH_TEMP;
                }
            }
        }
        end_c_loop_all(c, t)
    }
}

// 设置边界条件的质量流率
DEFINE_PROFILE(mass_flow_rate_profile, t, i)
{
    face_t f;
    Thread* c_thread;
    cell_t c;

    // 遍历边界上的所有面
    begin_f_loop(f, t)
    {
        c = F_C0(f, t);  // 获取相邻单元
        c_thread = THREAD_T0(t);  // 获取相邻单元线程

        // 检查单元是否被标记为入口边界单元
        if (C_UDMI(c, c_thread, 0) == 2.0)
        {
            // 获取单元的温度
            real temp = C_T(c, c_thread);

            // 检查温度是否达到点火条件
            if (temp >= IGNITION_TEMPERATURE)
            {
                F_PROFILE(f, t, i) = 12;  // 设置高温燃气的质量流率
            }
            else
            {
                F_PROFILE(f, t, i) = 0.0;   // 如果未达到点火条件，则保持原有流率（例如，0）
            }
        }
    }
    end_f_loop(f, t)
}