UDF救助

lihua_0629

二维喷管入口总温，总压随时间变化如下:
T(t)=300+150000*t        （0≤t≤0.02）
T(t)=3300

P(t)=200000+30000000*t  （0≤t≤0.02）
P(t)=800000

求DEFINE-PROLIFE     UFD  程序。

lihua_0629

各位高手，没人懂这个吗？浪子兄呢？

浪子轻狂

没大看懂你这个分段函数的意思：
T(t)=300+150000*t        （0≤t≤0.02）
T(t)=3300
t＞0.02，T就一直=3300，还是什么意思？

znon

如果没猜错的话，下面的程序应该有帮助。
不过，我想说的是，这个程序，直接看帮助就可以找到，我只是简单修改了一下，建议提问前自己多做些功课。
#include "udf.h"

#define TA 300.
#define TB 150000.
#define TC 3300.

DEFINE_PROFILE(T, t, nv)
{
  face_t f ;

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

  if (flow_time <= 0.2 )
    {

      begin_f_loop(f,t)
       {
         F_PROFILE(f,t,nv) = TA + TB. * flow_time ;
       }
      end_f_loop(f,t)
     }
  else
    {
   
      begin_f_loop(f,t)
        {
          F_PROFILE(f,t,nv) = TC ;
        }
      end_f_loop(f,t)
     }
}

浪子轻狂

znon兄，你这个udf，把if……else语句写在begin……end这一个loop中不是更简单吗？

znon

语言上应该是更简单，我是把UDF帮助中的例子简单改了一下，改完后，我也想把时间判断写在loop中，后来想是不是这样效率会高一些。
要是判断写在里面，循环内每次都需要判断一下；这样写，判断一下时间后，直接计算。
不知道我理解的对不对，请多指导。

浪子轻狂

一个是内部判断，一个是外部判断。过程都是：判断时间+遍历循环，效率不会有什么差异吧。当然，没做过测试，个人习惯，我都是写里面，有空多交流！

znon

我也没有测试，以后多交流。

invader

浪子的提法 要 更好一些，