原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
基于matlab及龙格库塔法求解布拉修斯方程
Runge—Kutta法求解布拉修斯解
摘要
薄剪切层方程主要有三种解法，即相似解，非相似条件下对偏微分方程组的数值解和近似解。布拉修斯解是布拉修斯于1908年求出的，它是零攻角沿平板流动的相似解。本文用四阶Runge—Kutta法求解高阶微分方程的方法，并用matlab编程实现，求得了与实际层流边界层相符合的数值解。
关键词：布拉修斯解，相似解，Runge—Kutta法，数值解。
1 布拉修斯近似解方程
二维定常不可压缩层流边界层的方程为：
（1）
（2）
边界条件为
将式（1）和式（2）进行法沃克纳—斯坎变换（简称F—S变换），将边界层方程无量纲化，即设
（3）
（4）
得出F—S变换后的动量方程
（5）
其中k为流动类型指标，横曲率项t为
（6）
m是量纲一的压力梯度参数，定义为
（7）
其边界条件变为
对于二维平面实壁流动（）可以忽略横曲率项t的轴对称流动，式（5）成为
（8）
根据相似解的定义，方程（8）中的函数f若式相似的，则它应只与η有关而与x无关，即对x的偏导数应为零。于是方程（8）应成为
（9）
若fw为常数，则方程（9）的边界条件为
；
2 布拉修斯解
布拉修斯于1908年求出了零攻角沿平板流动的解。这时
因而方程（9）成为
（10）
此即布拉修斯方程。对于实壁，，边界条件成为
；
3 Runge—Kutta法求解
Runge—Kutta通过将高阶微分方程化为一阶线性方程组，从而解出高阶方程的数值解。在方程（10）中令
（11）
于是方程（10）变为
（12）
当区步长为h，有四阶Runge—Kutta的形式如下
（13）
使用matlab软件取步长为0.2，迭代100步视作η→无穷大。迭代到第40步左右就收敛了，迭代结果如下（本文附录有全程序源代码）
表格 1平板层流边界层方程的数值解
f 0 0 0 0.33206 0.2 0.006641 0.066408 0.33199 0.4 0.02656 0.13277 0.33147 0.6 0.059736 0.19894 0.33008 0.8 0.10611 0.26471 0.32739 1 0.16557 0.32978 0.32301 1.2 0.23795 0.39378 0.31659 1.4 0.32299 0.45627 0.30787 1.6 0.42033 0.51676 0.29667 1.8 0.52952 0.57476 0.28293 2 0.65003 0.62977 0.26675 2.2 0.7812 0.68132 0.24835 2.4 0.9223 0.72899 0.22809 2.6 1.0725 0.77246 0.20646 2.8 1.231 0.81151 0.18401 3 1.3968 0.84605 0.16136 3.2 1.5691 0.87609 0.13913 3.4 1.747 0.90177 0.11788 3.6 1.9295 0.92333 0.098087 3.8 2.116 0.94112 0.080126 4 2
文档评论（0）