平板附面层方程的数值解.pdfVIP

平板附面层方程的数值解 陈楚平 (中山大学应用力学与工程系，广州，510275) 一耍奉文从一组积分方程出发，用积分的近似公式把它化为一组非线性方程组．进而用单纯形 法求解，得到平扳的阻力分布． 关■词平板附面层方程，数值解，单纯形法 O前言 不可压缩粘性流体绕流平板的层流附面层是普朗特提出的一个经典的问题，勃拉休斯求 解这个问题时，把在y=6上，“=U的边界条件移至无穷处．本文从一组与哥鲁别娃积 分关系式相当的式子出发，利用端点导数表示积分的近似公式，得到一个非线性方程组，用 单纯形数值方法求解这个非线性方程组，从而得到平板所受的阻力分布，其结果与勃拉休斯 解比较是一致的． 1基本方程 平板附面层方程为 l“磊栅面刮萨 u’㈤ f“罢+专=V雾 J 1娶+娑：o． n’㈣ I面+面_o· L 边界条件为 “=柚=0， 当y=0，工≥O时； (3) “=【，， 当y=研工J时． (4) 式中：耐JJ。U分别为附面层厚度和沿平板方向的来流速度． 根据连续方程(2)，引入流函数y∽y)使 。：掣，u：一掣． ”亩瑚一言。 则动量方程(2)变为 ㈣… ay孤西孤av2曲3 等意一警窘=V争． 边界条件式(3)，(4)则变为 ry2詈乩当y。慵呲 l孚：u， 当y：占(工)时． (6) 设 妒∽y)=￡晒(工)，(叶) (7) 式中：tl=y，占(算) 则有 『-詈卅㈣争=等M守=等肼 { 7 (8) I晏￥：一孚咿”∽)，娑：【席(，m)一形-伪)) L孤av 占…”h。～……7 将式(B)代入方程式(5)，则有 P’W艿2． (9) 万”2咄 因为式(9)左边是竹的函数，右边是工的函数，所以A应为常数，并且有 堡：n兰． (10) dx U ，。l+丑扩”=0． (11a) 而边界条件式(6)化为 I，(0)=，