24小时热门版块排行榜    

查看: 1004  |  回复: 4

sdmengjia

金虫 (小有名气)

[求助] 求教一个一维非线性偏微分方程的离散方法 已有2人参与

有一个简单的一维非线性对流输运方程,如下图所示。其中p代表对流物质的密度,Vp代表物质的对流速度,是一个变量,但和x和t之间没有简单的函数对应关系。
现在需要把这个方程离散化。我从简单的方法入手选择了有限差分法中的迎风算法,但是发现当对流速度Vp为变量时,无论用低阶还是高阶的迎风算法还是其他的一些算法都发现存在对流物质总量不守恒的问题。而当Vp为恒定值时则不存在这个问题,请问怎样用合适的方法离散这个方程?
求教一个一维非线性偏微分方程的离散方法
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

m06z511

木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
sdmengjia: 金币+5, ★★★很有帮助 2015-10-13 12:38:32
建议:
1:首先确定方程是正确的。要求该方程物理意义明确,方程式是否需要其他的约束条件,方程式在什么范围内有意义等。
2:在物理意义确定,方程完整正确的条件下,再想办法进行离散,差分求解。
3:只要能进行差分求解,不管是迎风算法还是龙格库塔法等,其实没啥大的本质的区别。
4:关于偏微分方程,有专门的书籍可供参考,在数学上已经讨论很清楚了。我们需要的就是搞清楚前提条件,明白物理上的含义。
2楼2015-10-12 16:23:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sdmengjia

金虫 (小有名气)

引用回帖:
2楼: Originally posted by m06z511 at 2015-10-12 16:23:52
建议:
1:首先确定方程是正确的。要求该方程物理意义明确,方程式是否需要其他的约束条件,方程式在什么范围内有意义等。
2:在物理意义确定,方程完整正确的条件下,再想办法进行离散,差分求解。
3:只要能进行 ...

可以确定方程式正确且意义明确,没有其他约束条件。
Vp的值需要查一个数据曲线来确定。或者可以将它假定为一个随机值。
以我对这个公式的理解,Vp的不同只会影响物质输送的快慢,不会影响物质的总量。不知我这样理解是否正确?
能否赐教类似的书籍吗?
3楼2015-10-12 19:46:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shenyxtata

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
sdmengjia: 金币+5, ★★★很有帮助 2015-10-13 12:38:46
用有限体积法吧,就没有守恒的问题了。
4楼2015-10-12 20:06:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

m06z511

木虫 (职业作家)

引用回帖:
3楼: Originally posted by sdmengjia at 2015-10-12 19:46:08
可以确定方程式正确且意义明确,没有其他约束条件。
Vp的值需要查一个数据曲线来确定。或者可以将它假定为一个随机值。
以我对这个公式的理解,Vp的不同只会影响物质输送的快慢,不会影响物质的总量。不知我这样 ...

Vp的不同只会影响物质输送的快慢,不会影响物质的总量。这样的理解是正确的。本身就是输运方程。随机值也是定值。你的问题可以求解。书籍吗?请参考流体力学。
5楼2015-10-13 08:28:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sdmengjia 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 娱乐 +3 Tide man 2026-08-03 3/150 2026-08-04 09:58 by gltch
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 jRl3mE6ddZGq 2026-08-03 10/500 2026-08-04 09:26 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +9 锦衣卫寒战 2026-07-28 13/650 2026-08-04 07:21 by Ermito
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见