利用动网格模拟筒内刚体在气体作用下运动的问题

swj662002

问题描述：在一个圆柱形筒内，有一个类似于活塞的刚体，初始时刻位于筒的底端。燃气和蒸汽的混合气体从intlet进入。随着时间的推移，刚体在气体的作用下向上运动。
物理模型：如图
这里我有几个问题想请教：
1 如果建立二维模型计算该问题，动网格更新方法才采用动态分层法。可否用四边形网格划分流体区域？划分网格时，是要在活塞上部、下部都画分，还是只在活塞的下部划分网格？
2 如何处理活塞与圆柱筒壁之间的关系？是否要在两者之间划interface?
3 对于气体作用活塞运动，如何自定义UDF？可否提供一个相关的UDF例子？
我刚开始接触动网格，有很多基本的问题不懂。请知道的和做过的高手帮帮忙，在这里谢谢了！

gallongallon

当然都要划分
可以用interface，也可以不用
一个标准的CG-MOTION的UDF就可以满足要求吧，最多做活塞上下两个面压力计算来计算推动力

swj662002

万分感谢 gallongallon 的帮助
如果我用Dynamic Layering 动网格更新方法 可否用四边形划分网格呢？

youzhizhe

用四边形网格完全可以的

tongliang

能否传一个关于这个问题的UDF看看
我正好也碰见这个问题

swj662002

这是一个DEFINE CG MOTION 的例子
Example
Consider the following example where the linear velocity is computed from a simple force
/************************************************************
* 1-degree of freedom equation of motion (x-direction)
* compiled UDF
************************************************************/
#include "udf.h"
static real v_prev = 0.0;
DEFINE_CG_MOTION(piston,dt,vel,omega,time,dtime)
{
Thread *t;
face_t f;
real NV_VEC(A);
real force, dv;
/* reset velocities */
NV_S(vel, =, 0.0);
NV_S(omega, =, 0.0);
if (!Data_Valid_P())
return;
/* get the thread pointer for which this motion is defined */
t = DT_THREAD(dt);
/* compute pressure force on body by looping through all faces */
force = 0.0;
begin_f_loop(f,t)
{
F_AREA(A,f,t);
force += F_P(f,t) * NV_MAG(A);
}
end_f_loop(f,t)
/* compute change in velocity, i.e., dv = F * dt / mass
velocity update using explicit Euler formula */
dv = dtime * force / 50.0;
v_prev += dv;
Message ("time = %f, x_vel = %f, force = %f\n", time, v_prev,
force);
/* set x-component of velocity */
vel[0] = v_prev;
}
详细见附件

swj662002

tongliang

请教一下用interface怎么做这个算例子啊
急求