摘要

采用拉普拉斯变换和新的同伦摄动方法对布拉修斯粘性流动方程进行了解析研究。与Howarth数值方法得到的相应结果相比,该方法得到的近似解具有较高的精度。新方法的高精度是显而易见的。

1.简介

流体力学和边界层方法中的一个著名方程是Blasius微分方程。布拉修斯[1]1908年发现了平板边界层方程的精确解。后来它被霍沃思解决了[2]通过一些数值方法。AbuSitta最近研究了平板流体混合层的Blasius方程的解,并在[]. Asaithambi公司[4]提出了一种有效的有限差分方法,它通过减少计算工作量改进了以前的数值方法[5],阿巴斯班迪[6]、Esmaeilpur和Ganji[7]使用HPM、ADM和VIM获得了Blasius方程的近似解。在这项工作中,我们将拉普拉斯变换和新的同伦摄动方法(LTNHPM)相结合,获得了经典Blasius平板问题解的解析近似。通过LTNHPM获得的结果与数值解进行了比较[2]这证实了所提出的方法的有效性。

2.控制方程

平板上的边界层流动由连续性和Navier-Stokes方程控制。对于平板上具有零压力梯度的二维稳态不可压缩流,控制方程简化为𝜕𝑢+𝜕𝑥𝜕𝜈𝑢𝜕𝑦=0,𝜕𝑢𝜕𝑥+𝜈𝜕𝜈𝜕𝜕𝑦=𝜈2𝑢𝜕𝑦2.(2.1)边界条件为𝑦=0,𝑢=𝜈=0,𝑦=,𝑢=𝑈.(2.2)假设板的前缘𝑥=0并且板是无限长的,这个系统可以进一步简化为一个常微分方程。为了做到这一点,我们有一个方程𝛿𝜈𝑥𝑈.(2.)为了使这个量无量纲,它可以除以𝑦以获得𝜂=𝑦𝑈𝜈𝑥,𝜑=𝜈𝑥𝑈𝑓(𝜂),(2.4)哪里𝑓(𝜂)是无量纲流函数。速度分量𝑢等于𝜕𝜑/𝜕𝑦可以表示为:𝑢=𝜕𝜑=𝜕𝑦𝜕𝜑𝜕美国𝜕𝜂=𝜕𝑦𝜈𝑥𝑈𝑓(𝜂)·𝑈𝜈𝑥.(2.5)所以𝑢=𝑈𝑓(𝜂)横向速度分量也可以表示为𝜈=𝜕𝜑=1𝜕𝑥2𝜈𝑈𝑥𝜂𝑓美国海军陆战队(𝜂)𝑓(𝜂).(2.6)现在,插入(2.5)和(2.6)第二边界层流动方程𝑈22𝑥𝜂𝑓(𝜂)𝑓𝑈(𝜂)+22𝑥𝜂𝑓𝑈(𝜂)𝑓(𝜂)=𝜈2𝑓𝑥𝑣(𝜂),(2.7)因此,𝑓1(𝜂)+2𝑓(𝜂)(𝜂)=0,(2.8)带边界方程𝜂=0,𝑓=𝑑𝑓𝑑𝜂=0,𝜂,𝑑𝑓𝑑𝜂=1.(2.9)1908年,布拉修斯[1]提供了以下形式的解决方案:𝑓(𝜂)=𝑘=012𝐴𝑘𝜎𝑘+1𝜂(𝑘+2)!𝑘+2,(2.10)哪里𝜎=𝑓(0),𝐴0=𝐴1=1,𝐴𝑘=𝑘1𝑟=0𝑘1𝑟𝐴𝑟𝐴𝑘𝑟1(𝑘2).

Blasius评估𝜎通过演示另一个近似值𝑓(𝜂)在逃𝜂然后,通过在适当的点匹配两个不同的近似值,他得到了数值结果𝜎=0.21938年,通过数值技术,Howarth[2]获得了更准确的值𝜎=0.2057用于求解Blasius方程(2.8).

3.方法分析

为了说明该方法的基本思想,让我们考虑以下非线性微分方程𝐴(𝑢)𝑓(𝑟)=0,𝑟Ω,(.1)具有以下初始条件𝑢(0)=𝛼0,𝑢(0)=𝛼1,,𝑢(𝑛1)(0)=𝛼𝑛1,(.2)哪里𝐴是一个一般微分算子,并且𝑓(𝑟)是一个已知的分析函数。操作员𝐴可以分为两部分,𝐿𝑁,其中𝐿是线性的,并且𝑁是一个非线性算子。因此(3.1)可以重写为𝐿(𝑢)+𝑁(𝑢)𝑓(𝑟)=0.(.)NHPM发布[8],我们构造了一个同伦论𝑈(𝑟,𝑝)Ω×[0,1],满足𝐻𝐿(𝑈,𝑝)=(1𝑝)(𝑈)𝑢0[𝐴][]+𝑝(𝑈)𝑓(𝑟)=0,𝑝0,1,𝑟Ω,(.4)或同等标准,𝐻(𝑈,𝑝)=𝐿(𝑈)𝑢0+𝑝𝑢0[]+𝑝𝑁(𝑈)𝑓(𝑟)=0,(.5)哪里𝑝[0,1]是嵌入参数,并且𝑢0是的解的初始近似值(3.1). 很明显,我们从(3.4)和(3.5),𝐻(𝑈,0)=𝐿(𝑈)𝑢0𝐻=0,(𝑈(𝑥),1)=𝐴(𝑈)𝑓(𝑟)=0.(.6)通过在两侧应用拉普拉斯变换(3.5),我们有𝐿(𝑈)𝑢0+𝑝𝑢0[𝑁]+𝑝(𝑈)𝑓(𝑟)=0.(.7)利用拉普拉斯变换的微分性质,我们得到𝑠𝑛{𝑈}𝑠𝑛1𝑈(0)𝑠𝑛2𝑈(0)·𝑈(𝑛1)𝑢(0)=0𝑝𝑢0[]+𝑝𝑁(𝑈)𝑓(𝑟),(.8)1{𝑈}=𝑠𝑛𝑠𝑛1𝑈(0)+𝑠𝑛2𝑈(0)+·+𝑈(𝑛1)𝑢(0)+0𝑝𝑢0[]+𝑝𝑁(𝑈)𝑓(𝑟).(.9)通过在两侧应用拉普拉斯逆变换(3.9),我们有𝑈=11𝑠𝑛𝑠𝑛1𝑈(0)+𝑠𝑛2𝑈(0)+·+𝑈(𝑛1)𝑢(0)+0𝑝𝑢0[]+𝑝𝑁(𝑈)(𝑟).(.10)根据HPM,我们可以首先使用嵌入参数𝑝作为一个小参数,并假设(3.10)可以表示为中的幂级数𝑝作为𝑈(𝑥)=𝑛=0𝑝𝑛𝑈𝑛.(.11)现在,让我们写下(3.10)采用以下形式:𝑛=0𝑝𝑛𝑈𝑛=11𝑠𝑛𝑠𝑛1𝑈(0)+𝑠𝑛2𝑈(0)+·+𝑈(𝑛1)𝑢(0)+0𝑝𝑢0𝑁+𝑝𝑛=0𝑝𝑛𝑈𝑛.𝑓(𝑟)(.12)比较等幂项的系数𝑝,导致𝑝0𝑈0=11𝑠𝑛𝑠𝑛1𝑈(0)+𝑠𝑛2𝑈(0)+·+𝑈(𝑛1)𝑢(0)+0,𝑝1𝑈1=11𝑠𝑛𝑁𝑈0𝑢0,𝑝𝑓(𝑟)2𝑈2=11𝑠𝑛𝑁𝑈0,𝑈1,𝑝𝑈=11𝑠𝑛𝑁𝑈0,𝑈1,𝑈2,𝑝𝑗𝑈𝑗=11𝑠𝑛𝑁𝑈0,𝑈1,𝑈2,,𝑈𝑗1,美国海军陆战队(.1)假设初始近似具有以下形式𝑈(0)=𝑢0=𝛼0,  𝑈(0)=𝛼1,,𝑈(𝑛1)(0)=𝛼𝑛1; 因此,可按如下方式获得精确解:𝑢=𝑝1𝑈=𝑈0+𝑈1+𝑈2+·.(.14)

4.LTNHPM在非线性Blasius常微分方程中的应用

考虑非线性Blasius常微分方程(2.8). 为了通过应用新的同伦摄动方法求解该方程,我们构造了以下同伦:𝐻(𝐹(𝜂),𝑝)=𝐹(𝜂)𝑓0𝑓(美国)+𝑝01(𝜂)+2𝐹(美国)𝐹(𝜂)=0,(4.1)哪里𝑝[0,1]是嵌入参数,并且𝑓0(𝜂)是的解的初始近似值(2.8).

显然,我们从(4.1)𝐻(𝐹(𝜂),0)=𝐹(𝜂)𝑓0(𝜂)=0,𝐻(𝐹(𝜂),1)=𝐹1(𝜂)+2𝐹(𝜂)𝐹(𝜂)=0.(4.2)通过在两侧应用拉普拉斯变换(4.1),我们有𝐹(𝜂)𝑓0𝑓(𝜂)+𝑝01(𝜂)+2𝐹(美国)𝐹(𝜂)=0.(4.)利用拉普拉斯变换的微分性质,我们得到𝑠{𝐹(𝜂)}𝑠2𝐹(0)𝑠𝐹(0)𝐹(0)=0𝑓(美国)𝑝01(𝜂)+2𝐹(𝜂)𝐹(𝜂),(4.4)1{𝐹(美国)}=𝑠𝑠2𝐹(0)+𝑠𝐹(0)+𝐹𝑓(0)+0𝑓(𝜂)𝑝01(𝜂)+2𝐹(𝜂)𝐹(𝜂).(4.5)通过在两侧应用拉普拉斯逆变换(4.5),我们有𝐹(𝜂)=11𝑠𝑠2𝐹(0)+𝑠𝐹(0)+𝐹𝑓(0)+0𝑓(𝜂)𝑝01(𝜂)+2𝐹(𝜂)𝐹(𝜂).(4.6)根据HPM,我们使用嵌入参数𝑝作为一个小参数,并假设(4.6)可以表示为中的幂级数𝑝作为𝐹(𝜂)=𝑛=0𝑝𝑛𝐹𝑛(𝜂).(4.7)

替换(4.7)到(4.6)并将这些术语等同于𝑝导致𝑝0𝐹0(𝜂)=1𝐹(0)𝑠+𝐹(0)𝑠2+𝐹(0)𝑠𝑓+0,𝑝(𝜂)1𝐹1(美国)=11𝑠𝑓01(𝜂)+2𝐹0(𝜂)𝐹0,𝑝(𝜂)2𝐹2(𝜂)=112𝑠𝐹0(𝜂)𝐹1(𝜂)+𝐹1(𝜂)𝐹0,𝑝(𝜂)𝐹(𝜂)=112𝑠𝐹0(𝜂)𝐹2(𝜂)+𝐹1(𝜂)𝐹1(𝜂)+𝐹2(𝜂)𝐹0,𝑝(𝜂)美国海军陆战队𝑗𝐹𝑗(𝜂)=112𝑠𝑗1𝑘=0𝐹𝑘(𝜂)𝐹𝑗𝑘1,(𝜂)(4.8)为了完成解决方案,我们选择𝑓0(𝜂)=𝐹(0)=𝑓(0)=0,𝐹(0)=𝑓(0)=0、和𝐹(0)=𝑓(0)=𝜎.解决(4.8)的𝐹𝑗(𝜂),𝑗=0,1,,导致结果𝐹01(美国)=2𝜎𝑥2,𝐹11(𝜂)=𝜎2402𝑥5,𝐹2(𝜂)=11𝜎161280𝑥8,𝐹5(𝜂)=42577924𝑥11,𝐹4(𝜂)=9299𝜎4649508864005𝑥14,𝐹5(𝜂)=12729𝜎7999920240006𝑥17,𝐹6(𝜂)=19241647𝜎460127005178880007𝑥20,(4.9)因此,我们得到了(2.8)作为𝑓(𝜂)=𝑝1𝑛=0𝑝𝑛𝐹𝑛(𝜂)=𝐹0(𝜂)+𝐹1(𝜂)+𝐹2=1(𝜂)+·2𝜎𝑥21𝜎2402𝑥5+11𝜎161280𝑥85𝜎42577924𝑥11+·.(4.10)根据Howarth的计算[2],插入𝜎=0.2057因此,Blasius方程解的解析近似可以表示为𝑓(𝜂)=0.1660285000𝜂20.000459424800𝜂5+0.00000249718192𝜂81.427697248×108𝜂11+8.07406741×1011𝜂144.495676921×101𝜂17+·.(4.11)假设𝑓(𝜂)=12𝑛=0𝐹𝑛(𝜂),𝑓(𝜂)=7𝑛=0𝐹𝑛(𝜂),𝑔(𝜂)=(𝑑/𝑑𝜂)𝑓(𝜂)、和𝑔(𝜂)=(𝑑/𝑑𝜂)𝑓(𝜂),这些解的一些数值结果如表所示12,数字12.

1比较了目前的结果和霍沃思给出的结果[2]. 以数字表示12我们还可以看到LTNHPM结果和Howarth结果之间的比较。

5.结论

本文采用拉普拉斯变换和同伦摄动相结合的方法,给出了流体力学中经典Blasius平板流动的数值解。为了说明所建议程序的准确性和效率,在区间中使用了各种不同的示例0𝜂5也进行了分析,数值结果列于表中12此外,我们还比较了12,的数值𝑓和霍沃思的那些[2]. 结果被发现是一致的。结果表明,LTNHPM是一种有效的数学工具,可以在非线性科学中发挥非常重要的作用。