基于matlab的龍格库塔法求解布拉修斯方程.docVIP

Runge—Kutta法求解布拉修斯解 摘要 薄剪切层方程主要有三种解法，即相似解，非相似条件下对偏微分方程组的数值解和近似解。布拉修斯解是布拉修斯于1908年求出的，它是零攻角沿平板流动的相似解。本文用四阶Runge—Kutta法求解高阶微分方程的方法，并用matlab编程实现，求得了与实际层流边界层相符合的数值解。 关键词：布拉修斯解，相似解，Runge—Kutta法，数值解。 1 布拉修斯近似解方程 二维定常不可压缩层流边界层的方程为： （1） （2） 边界条件为 将式（1）和式（2）进行法沃克纳—斯坎变换（简称F—S变换），将边界层方程无量纲化，即设 （3） （4） 得出F—S变换后的动量方程