西安交大陶老师SIMPLER教学程序的一些疑问

yka8150 · 发表于 2005-3-19 13:39:28

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

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

x

在SIMPLER教学程序的example-8中（就是冲击旋转圆盘的那道题），有一些疑问。
在行号为502的上一行，
CON(I,J)=RHOM*RSWM**2/RMN(J)**3
给源项赋值的分母怎么是3次方呢？按书上的解释应该是1次方才对。
那位高手给解释一下，谢谢！！

xuyu · 发表于 2005-3-19 14:00:18

你看看下面的文章，有一篇我们也在讨论这个程序呢，我们可以一起讨论

xuyu · 发表于 2005-3-19 14:10:24

对了，是哪本书上的什么位置提到了？

yka8150 · 发表于 2005-3-19 16:01:33

我说的是陶文铨的《数值传热学》，书中P196页有提到这一内容，就是说动量方程的源项设置，不过和这道题还有点不一样。这道题的方程得根据完整的圆柱坐标系N-S方程整理一下。

yka8150 · 发表于 2005-3-19 16:09:51

还有，你说的下面的帖子是指哪一个呀？
下面是由一个相关的帖子，但那是讨论simpler算法本身的，而我是要问教学程序中自带的例子example-8中的问题。

xuyu · 发表于 2005-3-19 16:47:14

我没有找到你说的那一页上有这个内容，可能我们看的版本不一样，请你告诉我是哪一章哪一节，目前我也在看单一计算不可压流的速度和压力分布的问题，参考的就是圆盘的那个例子，我是说我们看的方向差不多

yka8150 · 发表于 2005-3-19 17:03:55

陶文铨《数值传热学》2001年5月版的
第6章求解椭圆型流动与换热问题的原始变量法
6.1动量方程的源项及流场求解中的关键问题
6.1.1动量方程中的源项
其中有个表，3种坐标系中U方程和V方程的源项。
对于这个例子，着重看圆柱坐标系的。书中给的是假定旋转方向速度为0的情况。
但这个例子还有些不同，因为考虑了旋转方向的速度，源项得从N-S方程中去分离出来。

xuyu · 发表于 2005-3-19 17:17:54

我用的是旧版的没找到，不过我印象里，V方程的源项应该是 -miu*v/r**2
那么你是怎么对它进行线性化的呢
考虑切向速度w谁不是要加上这样一项（记不太清楚了）rho*w**2/r

yka8150 · 发表于 2005-3-19 17:28:46

那你在看一下我在顶楼的那一行程序，分母是r**3，不是r

xuyu · 发表于 2005-3-19 17:35:57

RHOM 和 RSWM 分别代表什么意思？我还没看这

		自动登录	找回密码
密码			注册

西安交大陶老师SIMPLER教学程序的一些疑问

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