Falkner-Skan方程的近似解析解.pdf

第25卷第4期 计算力学学报 V01．25，No．4 Chinese of Mechanics 2008 Journal August 2008年8月 Computational Falkner—Skan方程的近似解析解 郑连存“， 温安国1， 张欣欣2 (1．北京科技大学数学力学系，北京100083f2．北京科技大学热能工程系，北京100083) Crocco变量变换将无穷区间的边界值问题转为初值问题并利用Pad6逼近技巧确定初值，给出了一种有效的解析 分解方法。进一步．本文设计了一种数值解法，将本文得到的近似解析解及数值结果与早期研究者Hartree等人 的结果进行了比较，证明了本文提出的解法的有效性和可靠性． 关键词：边界层；Falkner—Skan方程；Adomian拆分法；Crocco变换 中图分类号：0175．1 文献标识码：A 解析分解方法可以应用于求解许多非线性问题，具 1 Falkner-Skan方程 有广泛的理论意义和应用价值。 绕流楔型物体的Falkner-Skan方程是一个典 2 方程(1．1)的求解方法及过程 型的层流边界层相似解问题，1931年首先由V． 2．1 M．Falkner和S．W．Skan提出。Falkner-Skan方Adomian拆分法 程一直受到学术界的极大重视，很多流体力学工作 Adomian拆分法[卜叼又称逆算符法，是美国数 者都对其进行过研究。Falkner—skan方程表述为学物理学家G．Adomian提出的求解线性、非线性 (1．1) 厂+∥+fl(1一／2)=0 数学物理方程近似解析解的一个数学方法，其特点 边界条件为 是适用范围广，计算过程简单，收敛速度快，它不需 ，(o)=厂(o)=0，厂(∞)一1(1．2)任何近似条件，就可给出方程的高精度逼近解析解 式中卢为压力梯度，卢0表示顺压梯度；卢一0表 甚至精确解。对处理强非线性问题既不需要借助线 示压力梯度为零，卢0表示逆压梯度。 性、摄动、迭代或简化模型方程等途径，也不需要数 当楔形半顶夹角为叫2(o卢≤2)，如图1值方法。分解法的基本精神主要为：一是把一个真 所示，此时的物体前驻点附近的流动具有相似解。 解分解为若干个解分量之和，设法分别求出各阶解 本文目的在于寻求方程(1．1)在边界条件(1．2) 分量，然后让这些解分量之和以任意所需的高精度 下的近似解析解并探求解的物理意义。本文提供的 逼近真值；二是把整个方程恰当地分解为若干部 分。主要按照算符分解为线性、非线性、确定及随机 性各部分，原则上可以任意地分解，但要讲究技巧， 如所选的确定项线性算符是可逆的，从而易于求得 该线性算符相应方程的部分解，然后利用已知初值 或边值条件，从中设法找出方程中的其余部分解与