[求助] 请教gallongallon 高手，局部单向化出口的边界条件

niuaiming

请教高手gallongallong：这是出口边界条件的介绍，该怎样用UDF表示呢？
按照计算流体力学和数值传热学的方法，考虑出口边界上网格节点的参数值对于网格内邻近节点上的参数值无影响，即局部单向化坐标假定。故只给定流场压力的标定值。
    出口设在无回流处，流线方向与出口界面相垂直，除流场大气压力外，各变量在出口界面上沿流动方向的梯度为0。

waitingforu

这个好做，就是与出口面相临的内接点的值 赋值给 出口面对应点的值，就行了。
用到C0,T0。 如F_T(f,t) = C_T(c0,t0);

waitingforu

F_T（）要改成F_PROFILE()

niuaiming

好，谢谢楼上热心帮助，请看以下程序要表示成局部单向化边界条件，该怎样改写呢？请指教，谢谢！！
/****************************************************/
/*outlet.c*/
/****************************************************/
#include"udf.h"
DEFINE_PROFILE(outlet_u,thread,position)
{
cell_t c;
face_t f;
Thread *tc;
real U;
begin_f_loop(f,thread)
{
c=F_C0(f,thread);/*获得边界处相邻网格C0的cell index*/
tc=THREAD_T0(thread);/*获得C0的cell thread*/
U=C_U(c,tc);/*获得C0的速度分量U*/
F_PROFILE(f,thread,position)=U;
}
end_f_loop(f,thread)
}
DEFINE_PROFILE(outlet_k,thread,position)
{
cell_t c;
face_t f;
Thread *tc;
real K;
begin_f_loop(f,thread)
{
c=F_C0(f,thread);/*获得边界处相邻网格C0的cell index*/
tc=THREAD_T0(thread);/*获得C0的cell thread*/
K=C_K(c,tc);/*获得C0的湍动能K*/
F_PROFILE(f,thread,position)=K;
}
end_f_loop(f,thread)
}
DEFINE_PROFILE(outlet_d,thread,position)
{
cell_t c;
face_t f;
Thread *tc;
real D;
begin_f_loop(f,thread)
{
c=F_C0(f,thread);/*获得边界处相邻网格C0的cell index*/
tc=THREAD_T0(thread);/*获得C0的cell thread*/
D=C_D(c,tc);/*获得C0的耗散系数D*/
F_PROFILE(f,thread,position)=D;
}
end_f_loop(f,thread)
}

gallongallon

我给你的这个UDF就是这个意思，你没有看懂

niuaiming

那为什么初始化时，会出现错误
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()
[br][br][以下内容由 niuaiming 在 2006年11月16日 06:22pm 时添加] [br]
谢谢楼上二位帮助，谢谢！