Fluent UDF（六）时间变量宏

CFD 模拟过程中，有时需要访问和操控与时间相关的变量，例如在处理瞬态问题时，为了提高计算效率或满足特定的时间步要求，需要根据模拟时间动态调整时间步长。Fluent 提供了以下两种方式访问时间相关的变量：

• 使用求解器宏直接访问

• 使用 RP 变量宏间接访问。

求解器宏直接从求解器中提取当前的时间相关变量，如当前流动时间、当前物理时间步长等。下表列出了 Fluent 中用于访问时间相关变量的求解器宏。

宏名称	返回值
CURRENT_TIME	当前流动时间
CURRENT_TIMESTEP	当前物理时间步长
PREVIOUS_TIME	前一个流动时间
PREVIOUS_2_TIME	前两个时间步的流动时间
PREVIOUS_TIMESTEP	上一个物理时间步长
N_TIME	时间步数
N_ITER	迭代次数

注意：

• 使用 PREVIOUS_TIME 或 PREVIOUS_2_TIME 宏时， UDF 源代码中需包含 unsteady.h 头文件。
• N_ITER 宏只能在编译型 UDF (Compiled) 中使用。

某些时间相关变量，如当前物理流动时间，可以直接使用求解器宏（CURRENT_TIME）访问，或间接通过 RP 变量宏访问。

直接访问：

real current_time;
current_time = CURRENT_TIME;

间接访问：

real current_time;
current_time = RP_Get_Real("flow-time");

直接访问和间接访问宏对照表：

求解器宏	RP 变量宏
CURRENT_TIME	RP_Get_Real("flow-time")
CURRENT_TIMESTEP	RP_Get_Real("physical-time-step")
N_TIME	RP_Get_Integer("time-step")

在使用 UDF 时，不应在单元循环 (c_loop) 或面循环 (f_loop) 宏中使用 RP_GET_... 函数来访问 Scheme 变量。求解器是执行 CFD 模拟的核心计算引擎，而 Cortex 是处理用户界面和脚本的外部程序模块。两者之间的通信开销较大。在循环内部频繁调用这些函数会导致大量的通信操作，从而显著降低计算效率。为了提高 UDF 的执行效率，应该在进入循环之前获取所需的 Scheme 变量值，并在循环内部使用已经获取的值。

示例：时间相关入口速度分布

入口速度按照下列函数关系给出：
$$V_x=V_0+A\sin (\omega t);V_0=20m/s,A=5m/s,\omega =10rad/s$$
UDF 代码如下：

#include "udf.h"

// 定义宏常量
#define OMEGA 10.0  // 角速度
#define V0 20.0     // 基础速度
#define V_AMP 5.0   // 速度振幅

// 定义入口速度的 UDF
DEFINE_PROFILE(inlet_x_velocity, ft, var)
{
    real flow_time;  // 用于存储当前时间
    face_t f;        // 面索引类型

    // 获取当前流动时间
    flow_time = CURRENT_TIME;  // Fluent 提供的宏，用于获取当前时间

    // 遍历入口区域的每个面
    begin_f_loop(f, ft)  // 使用 Fluent 的面循环宏，遍历边界线程中的所有面
    {
        // 计算对应时间的速度
        real velocity = V0 + V_AMP * sin(OMEGA * flow_time);

        // 将速度值施加到面上
        F_PROFILE(f, ft, var) = velocity;
    }
    end_f_loop(f, ft)  // 结束面循环
}

示例： 动态时间步长

下面的 UDF 使用 DEFINE_DELTAT 宏在模拟中更改时间步长。使用 CURRENT_TIME 获取当前模拟时间（赋值给变量 flow_time）。在计算的前 0.5 秒内，设置时间步长为 0.1。对于剩余的模拟时间，设置时间步长为 0.2。

#include "udf.h"  

// 使用 DEFINE_DELTAT 宏定义 UDF，该宏用于更改时间步长值
DEFINE_DELTAT(mydeltat, d)
{
    real time_step;           // 定义变量来存储时间步长
    real flow_time = CURRENT_TIME; // 使用 CURRENT_TIME 宏获取当前流动时间，并赋值给 flow_time

    // 判断当前流动时间是否小于 0.5 秒
    if (flow_time < 0.5)
        time_step = 0.1;      // 如果是，则设置时间步长为 0.1 秒
    else
        time_step = 0.2;      // 否则，设置时间步长为 0.2 秒

    return time_step;         // 返回时间步长值给求解器
}

通过理解和应用时间变量相关宏，可以在 CFD 模拟中灵活地处理时间步长等时间变量，从而优化模拟过程，提高计算效率。本节介绍了这些宏的基本概念、使用方法以及注意事项，并通过示例代码展示了实际应用。