找回密码
 注册
查看: 839|回复: 0

【求助】用UDF添加电动势,为什么出现这种结果

[复制链接]
发表于 2016-7-24 04:11:34 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 CranberryCFD 于 2016-7-24 04:16 编辑

很简单一个模型,矩形截面的长直管道流动,壁面具有电势,且已知电势分布的泊松方程,想用UDF添加一个user define scalar

[math]\nabla^2\psi=\frac{-\rho_e}{\epsilon\epsilon_0}[/math]
大致思想是这样:因为是泊松方程,等式左边二阶导项就是扩散项,右边是个关于电势的函数,看做源项。
UDF里定义了扩散系数和源项。FLUENT里编译没有问题,非定常项和对流项都是none,边界条件是在壁面定义Specific value然后恒等于我设定的壁面电势,入口和出口定义为Specific flux,恒等于0。
Governing equation本身没有问题,我在matlab里差分得到的分布应该是这样(下图为管道横截面:壁面=壁面电势,bulk flow电势为0)

Matlab plotting

Matlab plotting

但是Fluent算出来后电势却处处恒等于壁面电势。

有遇到一样问题的朋友吗?在这个问题上卡了好几天了,急急TwT

下面是我的UDF,改了后缀名
psi_UDF.doc (3.02 KB, 下载次数: 14)
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表