CFD中文网

我没太看过相关的讨论。目前这种只有周期性边界才不会发生这种问题。你说的这个问题再往深了我就没关注过了，遗憾，老铁整的挺深啊:146:

@李东岳 尖角部位，加一个竖直的翼墙，这样平分后角度好点，但是附近网格不再是水平的，这样对interFoam计算有影响吗？

@李东岳 明白了，谢谢李老师耐心细致的解答。

@Jacobian 我尝试了您的方法，其实不需要设置rho.water、rho.air。我设置了水的密度随高度按照双曲函数变化，初始时刻以及边界密度是没有问题的。但是运算开始后出现问题，密度分布立即按照thermalphysicalProperties设置变化，也就是说按照目前的方法无法控制密度与其他物理场的联系。实际上我的密度场是可以随着流场变化的，因此我觉得应该添加一个浓度输运方程，然后重新定义热物性，使密度随之变化。

@bestucan
感谢您的回复！“不同方向的流动计算差异就大”这个解释很直接明了，我突然就有点明白这个网格线出现的原因了:xinxin:

兄弟，这个是树木简化为多孔介质吧，这个比较简单，不要用UDF就能实现，直接采用多孔介质模型就可以了，就是需要转化一下，多孔介质模型所填的阻力系数=2*LAD*Cd。另外，兄弟，我也是在做这个，我想在多孔介质区域加一个附加项，考虑树木对例子的吸附，会吗？欢迎加我讨论QQ1225220118

对简单的长方体进行网格划分，再检查网格质量时发现在长方体的周围有很多有问题的face产生，也不知原因，很疑惑，求大神指点，万分感谢！

0T5OC4QVS4OWJO998JSUAU7.png

顺便测试一下侧边栏回复

fluent初始化显示mass diffusivity：invalid diffusivity please check material。用的组分输运涡耗散模型

请问问题解决了吗？我的不管怎么加密也是这个趋势

@东岳 东岳老师，我知道您是化工方向，我当你的小迷弟已经快三年了，我是说您在网上总结应该是生化环材方向的吧哈哈哈哈啊哈哈哈哈啊哈哈哈哈🌚

缺这个包：binutils

或者路径没配置对，可以用echo $PATH看看有没有/usr/bin这一项

@星星星星晴 OpenFOAM中的DPM粒子追踪功能，在运动网格(ALE)上的支持并行计算吗？有没有什么限制？

@yuhx 请问前辈，这样的错误一般跟什么设置错误有关

@aiweimo 也行，有链接不

@hurricane007 我去，忘了。。。这个我当时打算自建一个免费的计算中心。后来发现一直需要给客户装机，电力都放在那面了，这个自己用的一直没搞起来。

我距离百亿市场只是因为电路容量不够 :136:

@李东岳 那理解了，谢谢李老师

你的图就是节点对齐。

有文献里是中心六面体网格，边界层为四面体问题网格。

你可以试试hex-core，在做体网格的时候，可以成圣hex-core

各位前辈好！最近在关于喷油器动网格问题上陷入了困境，希望能够得到大家的帮助。
先简单说一下喷油器的工作原理。图一为三维模型的剖面截图。
5eff8661692327fa42f84dfe24b6f95.png
5deae10e01c9e42e75ad3500a12d353.png
cd0e8a4b2e3be20ae6dee2b9c90b6dc.png
上面那个是球阀，下面那个是针阀，球阀、针阀最初均只受到液体力、弹簧力共同作用静止不动，然后给球阀施加一个电磁力，球阀开始运动，进而改变了流场的压力，从而使得下面的针阀也随之运动。研究的内容主要是通过向fluent中导入模拟两个阀门运动的udf文件，如下。球阀udf```
#include "udf.h"
#include "dynamesh_tools.h"

static real v_prev = 0.0;//initial velocity
static real loc_prev = 0.0;//spring deformation

DEFINE_CG_MOTION(ballball_six_three_one_parallel3, dt, vel, omega, time, dtime)
{
real ELECTROMAGNETIC_FORCE_MAGNITUDE = 155.0;
real SPRING_FORCE_MAGNITUDE = 100.0;
real SPRING_STIFFNESS = 20000.0;
real MASS = 0.00184;
real DISPLACEMENT_LIMIT = 0.0001;
real f1 = 0.0;
real f2, f_s, dv;
real NV_VEC(A);

#if !RP_HOST
Thread *t;
face_t f;
NV_S(vel, =, 0.0);
NV_S(omega, =, 0.0);
if (!Data_Valid_P())return;

t = DT_THREAD(dt); begin_f_loop(f, t) { if (PRINCIPAL_FACE_P(f, t)) { F_AREA(A, f, t); f1 = f1 + F_P(f, t) * A[1]; } } end_f_loop(f, t)

#endif

#if RP_NODE
f1 = PRF_GRSUM1(f1);
#endif

node_to_host_real_1(f1);

#if !RP_NODE
if (time <= 0.0001)
{
if (loc_prev < DISPLACEMENT_LIMIT)
{
f_s = SPRING_FORCE_MAGNITUDE + SPRING_STIFFNESS * loc_prev;//spring force
f2 = f1 + ELECTROMAGNETIC_FORCE_MAGNITUDE - f_s;
dv = dtime * f2 / MASS;
v_prev = v_prev + dv;
loc_prev = loc_prev + v_prev * dtime;
}
else
{
v_prev = 0.0;
}
}
if (time > 0.0001)
{
if (loc_prev > 0)
{
f_s = SPRING_FORCE_MAGNITUDE + SPRING_STIFFNESS * loc_prev;//spring force
f2 = - f_s;
dv = dtime * f2 / MASS;
v_prev = v_prev + dv;
loc_prev = loc_prev + v_prev * dtime;
}
else
{
v_prev = 0;
}
}
Message("\n\ntime=%f,y_vel=%f,f1=%f,f2=%f,f_s=%f,loc_prev=%f,dv=%f\n", time, v_prev, f1, f2, f_s, loc_prev, dv);
#endif
host_to_node_real_1(v_prev);
vel[1] = v_prev;
}

针阀udf：

#include "udf.h"
#include "dynamesh_tools.h"

static real v_prev = 0.0;//initial velocity
static real loc_prev = 0.0;//spring deformation

DEFINE_CG_MOTION(needle_six_three_one_parallel2, dt, vel, omega, time, dtime)
{
real SPRING_FORCE_MAGNITUDE = 30.0;
real SPRING_STIFFNESS = 20000.0;
real DISPLACEMENT_LIMIT = 0.00018;
real MASS = 0.01;
real f1 = 0.0;
real f2, f_s, dv;
real NV_VEC(A);

#if !RP_HOST
Thread *t;
face_t f;
NV_S(vel, =, 0.0);
NV_S(omega, =, 0.0);
if (!Data_Valid_P())return;

t = DT_THREAD(dt); begin_f_loop(f, t) { if (PRINCIPAL_FACE_P(f, t)) { F_AREA(A, f, t); f1 = f1 + F_P(f, t) * A[1]; } } end_f_loop(f, t)

#endif
#if RP_NODE
f1 = PRF_GRSUM1(f1);
#endif

node_to_host_real_1(f1);

#if !RP_NODE
f_s = SPRING_FORCE_MAGNITUDE + SPRING_STIFFNESS * loc_prev;
f2 = f1 - f_s;
dv = dtime * f2 / MASS;
if (f2 <= 0.0)
{
if (loc_prev <= 0.0)
{
v_prev = 0.0;
}
else
{
v_prev = v_prev + dv;
loc_prev = loc_prev + v_prev * dtime;
}
}
else
{
if (loc_prev < DISPLACEMENT_LIMIT)
{
v_prev = v_prev + dv;
loc_prev = loc_prev + v_prev * dtime;
}
else
{
v_prev = 0.0;
}
}
Message("\n\nn_time=%f,n_y_vel=%f,n_f1=%f,n_f2=%f,n_f_s=%f,n_loc_prev=%f,n_dv=%f\n",time,v_prev,f1,f2,f_s,loc_prev,dv);
#endif
host_to_node_real_1(v_prev);
vel[1] = v_prev;
}

但是导入计算后计算到一半fluent总是卡在一半，没有报错， 也没有退出，就是卡在一个界面无法继续算下去，就是上图 的一直卡在0.000077s，就算不下去了，等了快半周了还是在 那个位置，不知道是代码的原因还是什么原因。 希望能够得到前辈们的解答。万分感谢！！！！