收藏本站
《吉林大学》 2004年
收藏 | 手机打开
二维码
手机客户端打开本文

若干时间相关、非线性偏微分方程的数值解法

尹丽  
【摘要】:大多数的自然现象与物理规律都可归结为与时间相关的非线性偏微分方程,诸如热传导方程、声波与弹性波方程、反应扩散与对流扩散方程、流体与气体力学方程组等发展方程。绝大多数情形这些问题的解不能用解析的公式表达出来,因此研究这些发展方程的数值求解问题备受人们关注。 在本文中我们主要构造了三类与时间相关的偏微分方程的数值计算方法,并作了相应的理论分析。本文主要内容可分为三部分,分别研究了含源项的浅水波方程组,描述浅水域表面小振幅、长波传播的“good”Boussinesq(简记为GB)方程,以及四阶非线性扩散模型Cahn-Hilliard方程的数值计算方法。 1.浅水波问题及其数值计算方法 考虑在矩形通道中孤立障碍物上方的一维流体运动问题,此时的数学模型为: h_t+(hu)_x=0, (hu)_t+(hu~2+1/2gh~2)_x=-ghB'(x),其中常数g代表重力加速度,障碍物B(x)是关于原点(x=0)对称的凸函数。h代表在障碍物上方流体的厚度,u代表流体水平方向的速度。当h(x),u(x)充分光滑时,系统(1)可以转化成完全散度形式: h_t+(hu)_x=0, u_t+(1/2u~2+gh+gB)_x=0系统(1)与系统(2)的稳态解(又称定常解)(h(x),u(x))都满足如下平衡条件: h(x)u(x)=const, 1/2u~2(x)+gh(x)+gB(x)=const, 为了便于对定常间题进行理论分析,引进如下无量纲变量: 凡一揣,M*一告“一;睽+1,Mc一会 其中ho和。。分别代表来流的流体厚度与速度,Hc代表障碍物顶点的高度. 我们证明了连续定常解的存在唯一性(见定理l)以及在典型情形下含间断定常解的 存在性(见定理2和定理3). 定理1(连续定常解的存在唯一性)对任意给定的参数0凡笋1,当McM*时, 定常问题(3)存在唯一的一条连续解曲线(h(x),诚x)),x任R.并且由B(x)=B(一x)可 知,此时流体厚度h(x),速度可x)和自由表面侧x)十B(x)是关于原点(x=0)对称的函 数. 一些实验表明:当McM*时,定常解不再为连续函数,通常含有间断.对于间断 定常解的分析有如下基本假设: (i)临界状态出现在障碍物的顶点位置,即tL。=两孚二,其中功。和。。分别代表流体 在障碍物顶点处的厚度与速度. (ii)在间断点处,定常解满足激波间断条件. 基于上述基本假设,定常解的结构分析可按上游和下游两部分进行. 首先,考虑流体在障碍物上游(一cox兰0)发生了激波间断,并记流体在间断线后 方的稳态厚度和稳态速度分别为hA和。A,则上述上游稳态解满足下列方程: hAu注一houo hA一ho c‘一。一恤了亲磊, (F卜H间斯条件), (4) +必c+刀七 uc功。=uAhA 。。=斌万砚, 一爵+h一凡,(平衡条件), =怕, (5) 屺一助 (临界条件). (6) 定理2当McM*时,存在唯一的数组(hA,。A,c:,功。,恢)〔R”满足物理条件:hAh0, 。,二。并使方程(4)一(6)成立.因此,在基本假设(i),(11)T,当Mc*时,在降 碍物上游存在含间断的定常物理解(城x),u(x)),一cox兰0. 2 其次,考虑一类典型的下游定常解结构.这里,将上游定常解按平衡条件延伸到下 游,由此在障碍物的后面获得一个新的稳态速度tLB与稳态厚度hB.由于障碍物对定常 解的影响是局部的,故流体在远离障碍物的状态将恢复为来流状态(h。,u。).这与黎曼问 题相似,我们考虑状态(hB,。B)与状态(ho,坳)通过一个激波间断和一个2一稀疏波连 接,则上述下游稳态解满足下列方程: 二2 .2 箭+“”一Ka‘一箭+“A),(平衡条件), tLBhB二从(=。AhA), (7) h刀牡刀一hxtLx hB一hx (R-H间断条件), (8) Ct.=牡B一hx 。x一2而石 (2一稀疏波). (9) 定理3当McM*时,。标h0且助十标了亲藉咐ho斌薪藉.成立,肪在 唯一的数组(hB,。B,hx,ux,cr)〔R“满足物理条件:hBhA,。A。B,hoxB, 二x二B,并使方程(7)一(g)成立.从而,在基本饭设(O,(词下,当Mc*时,在 降碍物下游存在一个由激波与稀硫波相连接的定常物理解(城x),可x)),Ox+co. 从浅水波方程组的完全散度形式(2)出发,用迎风差分格式、流矢量分裂法和基于 流失量分裂法的TVD格式来模拟各种情形下的非定常解.数值试验结果与前面的理论 分析相吻合.并且证明了迎风差分格式和流矢量分裂法的数值近似解在各个界面处满足 具有二阶精度的平衡条件(见命题1和命题2).进而证明了流矢量分裂法保持了厚度变 量h非负性的物理要求(见定理,). 从浅水波方程组的完全散度形式出发,在单元Ix,一1/2,肠+1/z]上进行积分,并对时 间用向前Euler差分格式离散后,得到如下全离散格式: 0典 (l时 畔十‘一岭 矛,夕.,,。一尸夕·价 朴二竺认王一竺- 其中叮为t。时刻在单元冈一1/2,xj十1/aj上的U一(h,司T的平均值,峨1/2为‘ 定义在单元边界xj+1/:处的数值流量F= 1。__‘甲 (h“,乏“‘+gh+g万)‘ 3 迎风差分格式中数值流量象:/2定义为: 瓜:/2一叠!二(二卜;(二1)卜看s、·(,麒毛些生))!二(二1卜二(二)}·(11) 其中 Sign(A)=TSign(A)T一’,(Sign(A))灯=Sign(a‘)占巧,A
【学位授予单位】:

知网文化
【相似文献】
中国期刊全文数据库 前20条
1 楚泽甫 ,吴集基;弗雷德霍姆方程的一种数值解法[J];河南大学学报(自然科学版);1988年04期
2 姚九存;;一类多维的线性Fredholm积分方程的数值解[J];科技创新导报;2009年29期
3 谢柳辉;;用加权余量法求非线性振动系统的定常解[J];铁道科学与工程学报;1985年03期
4 向新民,陈安乐;高维广义Kuramoto-Sivashinsky型方程全离散Fourier拟谱方法的大时间性态研究[J];高等学校计算数学学报;1998年01期
5 陈方权;;奇异积分方程的一种数值解法——离散涡旋法[J];四川师范大学学报(自然科学版);1988年S1期
6 黄明瑞;;一类含积分的非线性方程的数值解法[J];昆明理工大学学报(理工版);1986年03期
7 贾红艳;;第二类Fredholm积分方程的小波解法[J];安阳师范学院学报;2010年02期
8 阿不都热西提·阿不都外里;对于热传导方程的一种新数值解法的定式化[J];新疆大学学报(自然科学版);1996年02期
9 石东洋;裴丽芳;;细菌模型的非协调有限元分析[J];应用数学学报;2011年03期
10 邹继高;关于抛物型方程非线性边值问题的一个数值解法[J];数值计算与计算机应用;1983年01期
11 姜波,徐家旺;非线性代数方程组的数值解法比较[J];沈阳航空工业学院学报;2003年03期
12 陈文娟;;论方程根的数值解法[J];现代商贸工业;2009年13期
13 胡裕秀;周仁根;;考虑徐變影響的正交異性板的近似解[J];太原理工大学学报;1963年02期
14 胡海昌;广义变分原理在近似解中的合理应用[J];力学学报;1982年01期
15 吴微;与时间相关的二阶双曲方程的全离散有限元逼近[J];吉林大学学报(理学版);1985年02期
16 陈四海;;SLY—1.8型双头离心摇床的振动分析[J];鞍山科技大学学报;1986年02期
17 李才中,张玉平,严国政;非齐次弹性力学方程组近似解的收敛性[J];数学物理学报;1989年04期
18 郑韫瑜;障碍问题的一个非协调有限元[J];复旦学报(自然科学版);1989年01期
19 李莉;;关于罚函数法与障碍函数法的改进[J];曲阜师范大学学报(自然科学版);1991年03期
20 张近苇;;渐变折射率光纤中光波场方程的近似解与严格解[J];暨南大学学报;1992年03期
中国重要会议论文全文数据库 前10条
1 亢一澜;王娟;;时间相关界面力学性能表征与界面参数反演识别[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
2 孙国芹;尚德广;;时间相关单多轴疲劳寿命预测[A];第十五届全国疲劳与断裂学术会议摘要及论文集[C];2010年
3 丁俊;康国政;刘宇杰;郭严;;35CrMo钢单轴棘轮-疲劳交互作用实验研究[A];中国力学学会学术大会'2009论文摘要集[C];2009年
4 王娟;亢一澜;秦庆华;李晓奇;;基于实验的反演识别方法及其在黏接界面力学性能表征中的应用[A];中国力学学会学术大会'2009论文摘要集[C];2009年
5 沈建和;陈树辉;;一类三维非线性自治系统极限环分析的改进L-P法[A];第九届全国振动理论及应用学术会议论文集[C];2007年
6 邵雪娇;康国政;郭素娟;;颗粒增强金属基复合材料时间相关棘轮行为的有限元模拟[A];第十五届全国复合材料学术会议论文集(下册)[C];2008年
7 吴宠雄;徐炜;;无耗二口网络参数的直接寻优数值解法[A];1989年全国微波会议论文集(上)[C];1989年
8 许婧;傅克祥;李象远;;Stokes频移的时间相关性[A];第九届全国化学动力学会议论文摘要集[C];2005年
9 李杰超;魏德敏;;大跨索网结构风速时程数值模拟[A];第15届全国结构工程学术会议论文集(第Ⅲ册)[C];2006年
10 张建明;戎卫国;;马丁-候方程数值解法的研究[A];2007年山东省制冷空调学术年会论文集[C];2007年
中国重要报纸全文数据库 前10条
1 潘 峰;大盘又到变局时?[N];证券日报;2004年
2 四川 严中剑;车载CD机不工作检修[N];电子报;2005年
3 丁心怡;老人用“阿奇”注意防过敏[N];医药养生保健报;2007年
4 本报记者 陈刚;阳光文化主攻出版分销[N];中国图书商报;2003年
5 ;南京工业大学自主创新成果[N];中国化工报;2007年
6 ;维系生命之源[N];科技日报;2006年
7 陈锦屏;有些药物易致心肌炎[N];健康时报;2006年
8 泰克公司中国市场经理 张权;实时频谱分析仪:RF测试有力工具[N];中国电子报;2006年
9 马文方;感知Google下一代搜索痕迹[N];中国计算机报;2007年
10 记者  杨茜 通讯员  孙水观;我市向“五小”行业“开刀”[N];嘉兴日报;2006年
中国博士学位论文全文数据库 前10条
1 尹丽;若干时间相关、非线性偏微分方程的数值解法[D];吉林大学;2004年
2 金承日;某些延迟微分方程的数值方法[D];哈尔滨工业大学;2006年
3 张秀峰;超短脉冲激光调制在时间高分辨光谱仪中应用的研究[D];中国科学院研究生院(长春光学精密机械与物理研究所);2005年
4 陈守婷;半离散AKNS系统的对称与代数结构[D];上海大学;2011年
5 王秀花;非线性方程的一些数值解法及其理论分析[D];上海大学;2011年
6 曾光;椭圆型方程的数值解法与稳定性分析[D];电子科技大学;2011年
7 李焕荣;土壤水动力系统的数值解法的分析与应用研究[D];首都师范大学;2007年
8 张冉;自由流和多孔介质流耦合问题的数学和数值分析[D];山东大学;2011年
9 白冬梅;非线性Schr(?)dinger方程及其相关耦合方程的数值解法[D];南京航空航天大学;2010年
10 王泽军;二维轴对称活塞问题激波解的存在性[D];复旦大学;2004年
中国硕士学位论文全文数据库 前10条
1 斯琴;同伦摄动法在非线性方程求解中的应用[D];内蒙古师范大学;2008年
2 杜萍;延迟微分方程的再生核数值解法[D];哈尔滨工业大学;2006年
3 卢明;改进的近似解析离散化方法及弹性波波场模拟[D];清华大学;2004年
4 豆福全;非线性演化方程的几种求解方法研究[D];西北师范大学;2005年
5 丁玲;二维水流水质黎曼近似解模型性能分析及应用[D];河海大学;2003年
6 肖曼玉;大型周期块三对角线性方程组的并行算法[D];西北工业大学;2006年
7 姚金兰;一类不存在本质连通区的问题的研究[D];贵州大学;2006年
8 吴荣俏;一类可压缩的Heisenberg链方程组的数值解[D];华南师范大学;2007年
9 郑生富;有损均匀传输线的时域有限差分法研究[D];重庆大学;2007年
10 刘晓蕾;局部相平衡模型(LPEM)方程数值解法的实现与分析[D];中国地质大学(北京);2007年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978