c - 如何在我的问题中访问单元格的压力值并使用 UDF 在 Ansys fluent 中将其设置为零？

Question

我正在尝试定义一个函数以在每次迭代时遍历每个单元格并检查该单元格的压力，如果压力为负，则将该单元格的压力设置为零。我该怎么做呢？

这个函数正确吗？

#include "udf.h"

DEFINE_ADJUST(pressure_zero, d)
{
    cell_t c;
    Thread *c_thread;
    begin_c_loop(c, c_thread)
        {
            if (C_P(c, c_thread) < 0)
            {
                   C_P(c, c_thread) = 0; 
            }
        }
    }
    end_c_loop(c, c_thread)
}