请高手解决动态网格的局限

fangdian

我想模拟焊接熔池的模拟，用到fluent，想用动态网格细化。
   我的基本思路是：
    加载移动高斯热源，热源移动到哪，哪个地方就动态的网格细化。

但是我看了一些fluent中的动态网格的例子，都是wall运动之类的，都是网格运动，不是网格细化的，或者针对某个边界，比如壁面等。

  难道fluent不能实现这个功能吗？
望各位高手指点迷津，最终我是想实现三维的，目前只要能实现二维的，就行。只要有个简单的例子，能说明问题，我绝对给分，谢谢。

水若无痕

fangdian 发表于 2011-12-17 15:55
我看过几个函数的介绍：
  1) void DEFINE_CG_MOTION (UDFname,Dynamic_Thread * dt,real vel[ ], real ome ...

既然是温度 加热，请考虑温度梯度网格自适应。

fangdian

我看过几个函数的介绍：
  1) void DEFINE_CG_MOTION (UDFname,Dynamic_Thread * dt,real vel[ ], real omega[ ], real time,real dtime)。
此函数接口用于控制刚体的运动,用户把刚体质心运动速度和角速度分别赋值给vel和omega, FLUENT根据它们的值来自动计算出边界下一步的位置,从而实现动边界的控制; 刚体质心的位置可以在函数接口界面对话框中定义。Dynamic Zones中的dwall就是要控制的动边界,Motion UDF/Profile中的就是UDFname,从中可看出它已被制定成用于控制dwall,理论上 FLUEN T可以通过这种方式实现无穷多个动边界的控制; C.G.Location用于设定初始位置的质心,C.G.Orientation用于设定刚体的初始角度。一般适用于刚体本身不变形的运动。
2) void DEFINE_GEOM(char name,Domain * d,Dynamic_Thread * dt,real * position)。
此函数接口用于控制变形体的边界运动, position就是运动边界上某网格节点的位置值,用户可以通过对其赋值达到控制效果, position [0]对应边界节点的x坐标, position [1]对应y 坐标, position [2]对应z坐标; FLUENT自动遍历所有的边界节点,因此适用于有规律的可以用函数描述的运动边界。
3) void DEFINE_ GRID_MOTION(name,d,dt,time,dtime)。
此函数接口也用于控制形体的边界运动。主要用于更加复杂的控制,用户需要自己利用 FLUENT提供的其他函数来遍历运动边界上的节点,并对其位置进行控制,因此 UDF编程比前面两种复杂得多。它甚至可以事先生成好边界数据,在计算中把数据读入,完成复杂形体控制。

   从上面的三个函数看，好像都是针对边界操作。

fangdian

谢谢水若无痕 ，只是能不能具体点呢？在fluent中的哪个命令或者有没有具体的简单的例子？

hgdlhp2002

可以交流一下！QQ：2284170287

fangdian

我看了一下网格自适应，好像只是在某个时刻菜单操作，不能计算过程中动态自适应？？

fangdian

全力呼叫水若无痕，呼叫。。。呵呵

dh06

不知道这个可以达到你想要的效果不？边界层随体运动：http://forum.simwe.com/thread-958359-1-1.html

水若无痕

fangdian 发表于 2011-12-20 19:38
我看了一下网格自适应，好像只是在某个时刻菜单操作，不能计算过程中动态自适应？？ ...

恩 我说的自适应 就是菜单下的。

好像不能实现你说的情况。

fangdian

dh06 发表于 2011-12-21 11:43
不知道这个可以达到你想要的效果不？边界层随体运动：http://forum.simwe.com/thread-958359-1-1.html ...

这个是边界层的，我这个好像。。。。。不过谢谢啊。

fangdian

水若无痕 发表于 2011-12-21 14:03
恩 我说的自适应 就是菜单下的。

好像不能实现你说的情况。

谢谢，那我这个有什么办法吗？

fangdian

有没有自适应宏之类的，或者自适应命令之类的？

fangdian

我按照水若无痕的思路，试了一下，同时勾选了Dynamic mesh,设置如下：

发现似乎有作用，因为运行过程中有：
Adapting mesh (Adapt Gradient of temperature)...提示
  部分内容如下：
...........
>
Flow time = 0.18s, time step = 18
12 more time steps

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   74 solution is converged
    74 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 2.9247e-07 1.3860e-08 1.6563e-06  0:00:01   20
    75 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 1.6453e-06 1.3355e-08 1.3746e-06  0:00:01   19
    76 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 1.6579e-06 1.3194e-08 1.1408e-06  0:00:01   18
!   77 solution is converged
    77 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 5.9588e-07 1.3102e-08 9.4701e-07  0:00:04   17
(ani-monitor-update 1 #t #f)()


>
Flow time = 0.19s, time step = 19
11 more time steps

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   77 solution is converged
    77 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 5.9588e-07 1.3102e-08 9.4701e-07  0:00:04   20

Adapting mesh (Adapt Gradient of temperature)...
done.!   78 solution is converged
    78 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 2.8240e-07 1.3015e-08 7.8619e-07  0:00:07   19
(ani-monitor-update 1 #t #f)()


>
Flow time = 0.2s, time step = 20
10 more time steps
...................
全部的运行过程在附件dymesh.doc中,只是我看不了网格的变化，设置动画也不行。
   我怎么看实际的网格变化呢？比方说动画变化图。
  除此之外，还有一个问题是：
   我按照温度梯度自适应，那个值是怎么确定的呢？目前只是按照一定的计算结果大概的确定一个。
  谢谢高人指点。
希望此问题早日结束。

  

bitren

可以先在热源周围细化网格，然后让细化的网格随热源运动！！！

fangdian

bitren 发表于 2011-12-21 22:53
可以先在热源周围细化网格，然后让细化的网格随热源运动！！！

这样好像不行吧？fluent没法使得网格移动，只是使得边界移动？？或者能说具体点吗？

fangdian

今天又试了一下，我将模型全部画成六面体网格，采用勾选了Dynamic mesh的方式，以及使用温度梯度自适应，看到效果了。如下：
  
可以看到，在热源附近网格自动加密。
  只是问题又来了：怎么动态显示加密过程？那个动画对网格的设置好像有些问题，我的设置过程如下：

  这个该怎么弄呢？

fangdian

再来一张局部放大图：

  以及中间执行的时候的部分信息。

......

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   69 solution is converged
    69 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 6.1212e-07 1.6942e-08 1.5262e-08  0:00:28   20

Adapting mesh (Adapt Gradient of temperature)...
done.!   70 solution is converged
    70 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 4.8798e-07 1.6942e-08 1.5260e-08  0:00:40   19
Flow time = 0.49999998s, time step = 50
10 more time steps

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   70 solution is converged
    70 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 4.8798e-07 1.6942e-08 1.5260e-08  0:00:42   20
!   71 solution is converged
    71 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 4.7971e-07 1.6942e-08 1.5260e-08  0:00:36   19
Flow time = 0.50999998s, time step = 51
9 more time steps

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   71 solution is converged
    71 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 4.7971e-07 1.6942e-08 1.5260e-08  0:00:38   20
!   72 solution is converged
    72 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 4.3107e-07 1.6942e-08 1.5260e-08  0:00:36   19
Flow time = 0.51999998s, time step = 52
8 more time steps

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   72 solution is converged
    72 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 4.3107e-07 1.6942e-08 1.5260e-08  0:00:38   20
!   73 solution is converged
    73 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 4.2484e-07 1.6942e-08 1.5260e-08  0:00:37   19
Flow time = 0.52999998s, time step = 53
7 more time steps

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   73 solution is converged
    73 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 4.2484e-07 1.6942e-08 1.5260e-08  0:00:39   20
!   74 solution is converged
    74 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 4.2615e-07 1.6942e-08 1.5260e-08  0:00:37   19
Flow time = 0.53999998s, time step = 54
6 more time steps

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   74 solution is converged
    74 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 4.2615e-07 1.6942e-08 1.5260e-08  0:00:39   20
!   75 solution is converged
    75 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 5.9921e-07 1.6942e-08 1.5260e-08  0:00:33   19
Flow time = 0.54999998s, time step = 55
5 more time steps

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   75 solution is converged
    75 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 5.9921e-07 1.6942e-08 1.5260e-08  0:00:35   20
!   76 solution is converged
    76 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 7.0496e-07 1.6942e-08 1.5260e-08  0:00:30   19
Flow time = 0.55999998s, time step = 56
4 more time steps

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   76 solution is converged
    76 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 7.0496e-07 1.6942e-08 1.5260e-08  0:00:32   20
    77 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 1.0910e-06 1.6942e-08 1.5260e-08  0:00:28   19
!   78 solution is converged
    78 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 6.0005e-07 1.6942e-08 1.5260e-08  0:00:29   18
Flow time = 0.56999998s, time step = 57
3 more time steps

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   78 solution is converged
    78 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 6.0005e-07 1.6942e-08 1.5260e-08  0:00:32   20
!   79 solution is converged
    79 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 1.1544e-07 1.6942e-08 1.5260e-08  0:00:32   19
Flow time = 0.57999998s, time step = 58
2 more time steps

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   79 solution is converged
    79 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 1.1544e-07 1.6942e-08 1.5260e-08  0:00:33   20
!   80 solution is converged
    80 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 6.7934e-07 1.6942e-08 1.5260e-08  0:00:29   19
Flow time = 0.58999998s, time step = 59
1 more time step

Updating solution at time level N... done.
  iter continuity x-velocity y-velocity z-velocity     energy          k    epsilon     time/iter
!   80 solution is converged
    80 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 6.7934e-07 1.6942e-08 1.5260e-08  0:00:31   20

Adapting mesh (Adapt Gradient of temperature)...
done.!   81 solution is converged
    81 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 5.4153e-07 1.6942e-08 1.5260e-08  0:00:42   19
Flow time = 0.59999998s, time step = 60
.......

   这个结果作为和大家分享吧。希望不要拍砖。。。

fangdian

再补一张中间运行效果图：

  采用fluent输出了一个MPEG，发现太大，200M，就不看了，输出MPEG的视图也不好，变成了：

真是麻烦，是不是fluent动画功能太差，为何视角突然不是设定的那样呢？
  
要是知道怎么生成小动画就好。

fangdian

本来是想做个动画给大家瞧瞧，可是fluent本身带的动画效果太差，输出MPEG文件不仅大而且视图角度不对，不是按照自己设定的角度，现在有一堆的case和data，每个case对应很多的data，不知道怎么导入到Tecplot中做成动画，郁闷。。。。。看来这个得告一段落？？再过两天的话还没有人回的话我就结贴得了。

水若无痕

fangdian 发表于 2011-12-23 11:48
本来是想做个动画给大家瞧瞧，可是fluent本身带的动画效果太差，输出MPEG文件不仅大而且视图角度不对，不是 ...

批量导入tecplot。很容易。