24小时热门版块排行榜    

查看: 2430  |  回复: 11

shenjianbiao

新虫 (小有名气)

[求助] 哪位大神给指导下如果求解一个偏微分方程组

各位大神帮忙看下这个偏微分方程,文献里给出的方法是Lax-Wendroff方法,但是具体实施起来该怎么弄呢?倒腾了好多天了,结果都不行。哪位大侠知道的话请不吝赐教,越具体越好,不胜感激。

我现在的做法是把方程整理为图2里面的形式用两步Lax-Wendroff法求解由h和u组成的偏微分方程组,但是结果发生震荡,不知道是我方法的问题还是为啥
哪位大神给指导下如果求解一个偏微分方程组
差分方程01.jpg


哪位大神给指导下如果求解一个偏微分方程组-1
差分方程02.jpg

文献的出处为:Bubble coalescence in pure liquids[ Last edited by shenjianbiao on 2013-10-13 at 21:37 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 1982_Bubble_coalescence_in_pure_liquids.pdf
  • 2013-10-13 21:37:06, 441.69 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

SEEDEN

银虫 (正式写手)

★ ★ ★ ★ ★
shenjianbiao(sweety代发): 金币+5, 辛苦应答,代楼主发点币币,以资鼓励。 2013-12-05 16:39:11
方程本身不复杂,但我很怀疑你的边界条件是否合适,你确定你提出的边界条件不多不少?
原文中好像没有R=0时的边界条件,你是怎么得到的?
如果你的边界条件没问题,在离散的时候,一般在柱坐标或者球坐标下,r=0时的边界条件需要特殊处理一下
8楼2013-10-15 19:45:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
5楼: Originally posted by shenjianbiao at 2013-10-14 10:19:39
还是不太明白,请大侠再说详细点,俺是菜鸟...

Lax-Wendroff方法的步长有两个。
6楼2013-10-14 11:01:13
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

调整步长:某一个、or 两个同时。
2楼2013-10-13 21:57:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

shenjianbiao

新虫 (小有名气)

引用回帖:
2楼: Originally posted by feixiaolin at 2013-10-13 21:57:36
调整步长:某一个、or 两个同时。

步长我已经用了0.0000001试过了,还是不行。不知道您说的“两个同时”是什么意思
3楼2013-10-14 09:15:31
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
3楼: Originally posted by shenjianbiao at 2013-10-14 09:15:31
步长我已经用了0.0000001试过了,还是不行。不知道您说的“两个同时”是什么意思...

步长有两个。
4楼2013-10-14 09:38:12
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

shenjianbiao

新虫 (小有名气)

引用回帖:
4楼: Originally posted by feixiaolin at 2013-10-14 09:38:12
步长有两个。...

还是不太明白,请大侠再说详细点,俺是菜鸟
5楼2013-10-14 10:19:39
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

shenjianbiao

新虫 (小有名气)

引用回帖:
6楼: Originally posted by feixiaolin at 2013-10-14 11:01:13
Lax-Wendroff方法的步长有两个。...

还是不太明白,我是热能专业的,数值方法方面确实欠缺太多了,不过还是多谢您了
7楼2013-10-15 18:03:54
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

shenjianbiao

新虫 (小有名气)

引用回帖:
8楼: Originally posted by SEEDEN at 2013-10-15 19:45:36
方程本身不复杂,但我很怀疑你的边界条件是否合适,你确定你提出的边界条件不多不少?
原文中好像没有R=0时的边界条件,你是怎么得到的?
如果你的边界条件没问题,在离散的时候,一般在柱坐标或者球坐标下,r=0时 ...

您好,首先很感谢您的回复,R=0处的边界条件是我依据 对称条件 得到的,因为数值计算的时候要用到,所以我就写出来了,原文里确实没有给;

r=0时的边界条件,我的处理过程为附图所示,不知道对不对,忘大神指导。

还有希望大神告诉我下用Lax-Wendroff格式求解我的问题的具体过程,麻烦您了
哪位大神给指导下如果求解一个偏微分方程组-2
差分方程03.jpg

9楼2013-10-17 11:08:55
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

SEEDEN

银虫 (正式写手)

引用回帖:
9楼: Originally posted by shenjianbiao at 2013-10-17 11:08:55
您好,首先很感谢您的回复,R=0处的边界条件是我依据 对称条件 得到的,因为数值计算的时候要用到,所以我就写出来了,原文里确实没有给;

r=0时的边界条件,我的处理过程为附图所示,不知道对不对,忘大神指导 ...

首先我不是大神,呵呵
差分方式很多种,哪怕最简单的向前向后中心差分都可以解决你的问题,未必要L-W格式。就算是L-W格式,偏微分方程数值解的书也有很多介绍,好好看看,再算一下看看
我认为你的在r=0时的边界条件处理是合适的,在接近0的地方是要处理一下,我也看过类似的处理方法
个人认为你最关键的问题还是边界条件的建立,跟问题切合的最合适合理的边界条件到底是什么。要根据物理情景好好考虑你的边界条件。我不做气泡方面的东西,也没法给建议。只是觉得你多找相关文献,看看别人对气泡是如何处理的。毕竟气泡除了热平衡,好像还有力平衡(如张力)等,不容易分析的。祝好运!

» 本帖已获得的红花(最新10朵)

10楼2013-10-22 14:21:41
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shenjianbiao 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见