龙格库塔算法

开心果 · 发表于 2005-7-8 18:02:55

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

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

x

请问，龙格库塔算法的显式算法应该就是一步一步算下来，不用再用直接法，迭代法等解方程组了把，确定一下，就开始编程了！！请高人不吝赐教！

开心果 · 发表于 2005-7-8 18:06:09

我记得anyone说过是这样的

anyone · 发表于 2005-7-9 06:41:57

楼下要买程序的是你的马甲换是你的老师? 不会这么巧吧, HEHE.
关于龙格库塔显式算法的问题, 其实我也有些疑问. 今天有时间和高手聊了一会, 澄清了一些看法.
搞低速不可压这么久了没听说谁用显式算法的. 纠其原因, 大盖是因为不可压对质量守恒要求太高, 用显式很难满足. 用龙格库塔做隐式每个时间步计算量太大. 可以考虑用CRANK-NICOLSON.时间是2阶精度.
非牛顿流有很多种. 如果只是粘性的变化, 可以试试还按牛顿流来解, 只是每次迭代时显式把粘性改了.
当然最好是先查查文献, 看别人是怎么做的, 可以少走很多弯路.

开心果 · 发表于 2005-7-10 19:30:00

哪个楼下？什么马甲？我不穿马甲的！哈哈
我现在又考虑还用simple，就是动量方程不能按照牛顿流体的线性关系变成N－S方程，直接带着应力求，每一步的应力根据本构关系用速度求出。

anyone · 发表于 2005-7-12 11:20:20

压力要不要分出来?
如果分出来的话, 就是我上次说的显式改粘性的方法. 这是目前很多人用的.

开心果 · 发表于 2005-7-12 11:29:12

压力分出来？我不明白什么意思阿！显式改粘度？粘度更新的时候用什么公式我不知道阿，是不是跟温度，压力，剪切速率都有关系？

anyone · 发表于 2005-7-12 11:41:25

应力中有一部分是压力.
粘度更新的时候用什么公式要根据流体介质属性来定. 可以跟温度，压力，剪切速率有关系. 用显式算粘度并不困难, 不过可能要迭代才能收敛.
你是要解决具体问题,还是要搞研究?

开心果 · 发表于 2005-7-12 12:06:14

ai,断消息太难用了，老是好不容易写好了，腾一下全没了，不知刚才断消息收到没有，我是搞具体问题

开心果 · 发表于 2005-7-12 12:12:40

粘度和T，P，剪切速率的关系要靠经验公式还是要把你的具体问题规到一种特定的非牛大的分类里面？我真佩服你，anyone，我自己都烦自己了。

anyone · 发表于 2005-7-12 12:23:23

收到了. 经验公式最好. 不过可以先用简单模型试试. 把塑性材料的大意先做出来.
我要下了. 下次再聊.

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

龙格库塔算法

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