求极坐标对流扩散方程推导 - 数学 - 计算数学 - 第3页小木虫论坛-学术科研互动平台

桩基工程及其耐久性

21楼2015-04-13 15:37:45

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

数学EPI: 6
应助: 413 (硕士)
贵宾: 0.002
金币: 7116.5
散金: 15
红花: 63
帖子: 1639
在线: 798.9小时
虫号: 2469437
注册: 2013-05-14
专业: 计算数学与科学工程计算

引用回帖:

21楼: Originally posted by yuezhuwen at 2015-04-13 15:37:45
可以告知一下这个是哪本书里的吗？...

这个不用书吧？我是自己写的。比较直截了当

22楼2015-04-14 12:46:33

yuezhuwen

铜虫 (小有名气)

feixiaolin: 屏蔽内容, 屏蔽此楼 2015-04-23 16:52:28

本帖内容被屏蔽

23楼2015-04-14 16:04:59

yuezhuwen

铜虫 (小有名气)

应助: 2 (幼儿园)
金币: 3455.6
散金: 85
红花: 2
帖子: 275
在线: 219小时
虫号: 1303441
注册: 2011-05-23
性别: GG
专业: 结构工程

引用回帖:

22楼: Originally posted by pippi6 at 2015-04-14 12:46:33
这个不用书吧？我是自己写的。比较直截了当...

亲不是知道看到的的留言了没？

桩基工程及其耐久性

24楼2015-04-17 14:30:00

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

数学EPI: 6
应助: 413 (硕士)
贵宾: 0.002
金币: 7116.5
散金: 15
红花: 63
帖子: 1639
在线: 798.9小时
虫号: 2469437
注册: 2013-05-14
专业: 计算数学与科学工程计算

引用回帖:

24楼: Originally posted by yuezhuwen at 2015-04-17 14:30:00
亲不是知道看到的的留言了没？...

对不起，因为出差没能及时回复。我看了一下，有以下几点
（1）方程的推导和速度方向（或符号）无关，在不同区域可以给定不同速度分布即可
（2）你的差分方程选择了逆风差分（如果v>0），这样容易形成不稳定性。不知道你的速度多大。而且扩散的一阶导数也是单边差分。如果，你的速度相对于扩散很大,局部Peclet number=v*dr/D >2, 那么混合差分或迎风差分(upwind)就是必要的；如果速度不大，局部Peclet number < 2，那么可以考虑中心差分。在任何场合都不应用逆风差分。而且一阶扩散导数总应该用中心差分
（3） matlab我不会
（4）你的计算结果需要解释一下，纵轴是浓度分布？不同曲线对应不同速度？

25楼2015-04-19 19:36:06

yuezhuwen

铜虫 (小有名气)

应助: 2 (幼儿园)
金币: 3455.6
散金: 85
红花: 2
帖子: 275
在线: 219小时
虫号: 1303441
注册: 2011-05-23
性别: GG
专业: 结构工程

引用回帖:

25楼: Originally posted by pippi6 at 2015-04-19 19:36:06
对不起，因为出差没能及时回复。我看了一下，有以下几点
（1）方程的推导和速度方向（或符号）无关，在不同区域可以给定不同速度分布即可
（2）你的差分方程选择了逆风差分（如果v>0），这样容易形成不 ...

对于计算结果两个曲线都是相同的速度只是有两种计算时间两种计算方法一种计算方法是沿径向分为对流扩散区-扩散区-对流扩散区三种区域，另一种方法是作为对比考虑沿径向为纯扩散区域

桩基工程及其耐久性

26楼2015-04-20 09:58:37

wuxun1989

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 168.5
帖子: 14
在线: 22.7小时
虫号: 2579505
注册: 2013-08-03
专业: 岩土与基础工程

引用回帖:

14楼: Originally posted by pippi6 at 2015-02-07 07:17:01
对微元体直接推导 diffusion equation 是一种很原始的推导方式，比较直观。

快照118.jpg

快照119.jpg
...

式（862）左边第二项的括号内是否多个r，请教用坐标变换的结果为什么不一样？

27楼2018-08-03 18:39:22