还是有关UDF的

zhangyonghai · 发表于 2003-9-25 21:35:40

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？注册

x

[这个贴子最后由zhangyonghai在 2003/09/26 09:42am 第 1 次编辑]

就是计算中的一些变量对x，y方向上面的积分，是怎么进行的？现在谢谢

cjt1234 · 发表于 2003-9-27 09:23:47

我是用的点积分，即该方向上的各个值相加，然后除以该方向上的值的个数。

badapple100 · 发表于 2003-9-27 22:39:00

下面引用由cjt1234在 2003/09/27 09:23am 发表的内容：
我是用的点积分，即该方向上的各个值相加，然后除以该方向上的值的个数。

晕，你这哪是积分呀，明明是求平均嘛。

zhangyonghai · 发表于 2003-9-29 20:31:53

大虾们，我该如何来实现这些呢〉

cjt1234 · 发表于 2003-9-30 22:02:31

我不是什么大虾，也是初学，：）你用我上面说的不是已经可以求出点平均了，求积分还有困难吗？给平均值乘上总面积，我就是这样做的。当然这对于网格不均匀的情况可能有点误差，但很小。其实FLUENT本身提供了一个求单元面积的函数，好像是F_AREA，但这个函数不是很好用，好像其相邻两个单元用该函数求的的值的差的绝对值是该单元面积的两倍（不一定对，你可以测试一下），用他应该可以实现面积分，但我觉着比较麻烦，你有精力可以试一下，欢迎共同探讨。

efang3721 · 发表于 2003-10-2 11:32:28

[这个贴子最后由efang3721在 2003/10/02 11:37am 第 2 次编辑]

我写一段“投影面加权积分”的程序，不知道是否合乎你的要求？
  face_t             /*face identifier*/
  Thread *t          /*pointer to face thread*/
  real vz[ND_ND], A[ND_ND], area_p, var;
  real total_area_p, total_var;
  real int_var;
  NV_D(vz, =, 0, 0, 1);/*the vector of xOy's plane*/

  total_var=0;
  totoal_area_p;
  begin_f_loop(f,t)
{
   F_AREA(A, f, t);    /* get area vector */
   var=F_var(f, t)    /* F_var general macro, e.g. F_T(f,t) for temperature*/
   area_p = NV_DOT(vz, A);    /* the area projection on xOy's plane */
   total_var+=var*area_p;    /* integral */
   total_area_p+=area_p;    /* integral of area */
}end_f_loop(f,t);
  int_var=totoa_var/total_area_p;  /* area-weight-integral */

另外具我所知F_AREA求的是面和向量，它的模的大小为面积的值。可以试试。

zhangyonghai · 发表于 2003-10-8 08:59:38

还是有些茫然，：（

dccat · 发表于 2010-1-8 16:52:02

这个程序好强大，用到了，谢谢

		自动登录	找回密码
密码			注册

还是有关UDF的

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。