Fluent UDF（八）DEFINE_ADJUST

DEFINE_ADJUST 是一个通用宏，用于调整或修改未作为参数传递的 Ansys Fluent 变量。例如，可以使用 DEFINE_ADJUST 来修改流动变量（例如，速度、压力）和计算积分。您还可以使用它在整个计算域内对标量进行积分，并根据结果调整边界条件。使用 DEFINE_ADJUST 定义的函数在每次迭代时执行，并在每次迭代开始时，在求解传输方程之前被调用。

宏形式	参数类型	返回值
DEFINE_ADJUST(name , d)	symbol name, Domain *d	void

DEFINE_ADJUST宏有两个参数：name和d：

• name：用户定义的 UDF名称： UDF的名称。
• Domain *d：指向其应用调整函数的域指针，由 Fluent 求解器传入。

DEFINE_ADJUST 定义的 UDF 编译完成后需要在 “User-Defined Function Hooks” 对话框进行挂载。

请注意，可以将多个 adjust 函数挂接到模型中。

示例：湍流耗散率求和

如下定义了一个名为 my_adjust 的 UDF 函数，该函数在每次迭代开始时被调用。它的作用是计算流场中所有单元的湍流耗散总和，并将结果打印到控制台中。代码通过遍历每个单元线程中的所有单元，累加每个单元的湍流耗散乘以其体积的值，最终得到整个流场的湍流耗散体积积分。

#include "udf.h"

// 定义一个通用的调整宏，用于在每次迭代中执行操作
DEFINE_ADJUST(my_adjust, d)
{
    Thread *t;          // 定义一个线程指针，用于遍历所有单元线程
    real sum_diss = 0.; // 初始化湍流耗散的总和
    cell_t c;           // 定义一个单元变量，用于遍历每个单元

    // 遍历所有单元线程
    thread_loop_c(t, d)
    {
        // 遍历线程中的每个单元
        begin_c_loop(c, t)
        {
            // 计算当前单元的湍流耗散并累加到总和中
            sum_diss += C_D(c, t) * C_VOLUME(c, t);
        }
        end_c_loop(c, t)
    }

    // 在控制台中打印湍流耗散的体积积分结果
    printf("湍流耗散的体积积分: %g\n", sum_diss);
}

示例：获取面心坐标

如下 UDF 获取指定面（由 FACE_ID 标识）的质心坐标，并在每次迭代时打印这些坐标信息

#include "udf.h"

// 定义 ADJUST 宏，每次迭代执行
DEFINE_ADJUST(face_centroid, domain)
{
    real FC[ND_ND];  // 定义一个数组用于存储坐标
    face_t face;  // 定义一个 face_t 类型的变量，用于存储面索引
    int FACE_ID = 16;  // 定义面 ID，在 Fluent 中查看，可以用宏常量定义：#define FACE_ID 16
    real fx, fy;  // 定义变量存储坐标值

    // 查找在domain中ID为FACE_ID的面线程face_thread
    Thread *face_thread;
    face_thread = Lookup_Thread(domain, FACE_ID);//后续数据获取宏介绍

    // 检查face_thread是否为空
    if (face_thread == NULL)
    {
        Message("Error: face_thread is NULL. Check FACE_ID.\n");
        return;
    }

    // 遍历 face_thread 中的所有面face
    begin_f_loop(face, face_thread)
    {
        F_CENTROID(FC, face, face_thread);  // 获取面 face 的质心坐标，并存储在 FC 数组中
        fx = FC[0];  // 获取 x 坐标
        fy = FC[1];  // 获取 y 坐标

        // 打印坐标信息
        Message("Centroid Coordinates: x = %f, y = %f\n", fx, fy);
    }
    end_f_loop(face, face_thread)  // 结束面循环
}

演示案例网格如下，几何尺寸 $20\times 40\text{mm}$ ，为直观展现上述 UDF 功能，案例网格尺寸设置为$5mm$ ，现演示如何使用面心坐标 (F_CENTROID)宏获取入口面面心坐标。

本案例网格比较少，可以直接得出面网格面心坐标为 $(0,2.5),(0,7.5),(0,12.5),(0,17.5)$ 此处单位为mm

导入网格后，可在边界条件信息栏查看面ID

将上述inlet ID=16 赋值给UDF中的FACE_ID，并将UDF文件保存为centroid.c，读入fluent进行编译。

编译完成后，ADJUST宏需要挂载使用

简单设置边界条件，inlet类型为速度入口，outlet为压力出口，其余默认。混合初始化完成后设置迭代次数为1

点击 Calculate 进行迭代计算，在控制台面板可查看输出结果与预想结果吻合。