资源描述：
《重庆大学高等流体作业流体作业.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、研究生课程考核试卷科目：高等流体力学教师：何川姓名：苗闪闪学号：专业：动力工程及工程热物理类别：学术考生成绩：卷面成绩平时成绩课程综合成绩阅卷评语：阅卷教师（签名）重庆大学研究生院制四、设的气体以的速度以零攻角的速度定常绕流长度为的大平板，试用数值解讨论边界层内的流动规律。解：设以匀速U运动的流体沿切线方向绕过一静止平板，由于流体的粘性作用，板面上的流体流速被降为零，板壁附近极薄区域的流体沿板的法向存在很大的速度梯度，该区域被称为流动边界层。，为大雷诺数绕流。沿板面取x坐标，板的发现方向取y坐标舍边界层的厚度为δ=δ(x)，因厚度很薄，既有。根据问题的特性，忽略质量力，
2、边界层内的速度矢量,压强，流体力学基本方程具体化为：C.E.(4.1)M.E.(4.2)(4.3)下面，应数量级比较的方法分析方程中各项的数量级关系。首先，选特征参数U、L，奖方程组中个变物无量纲化，并在无量纲化过程中注意使用各无量纲变量的取值范围为(0，1)或取1的数量级（记作(~1)）。即令：注：大雷诺数流动时，压力采用流动造成的动压作为参照具有相当的数量级关系。将各无量纲变量代入式（4.1），并比较其各项的数量级关系，有：将各无量纲变量代入式（4.2），并比较其各项的数量级关系，有：1将各无量纲变量代入式（4.3），并比较其各项的数量级关系，有：1根据物理依据和工
3、程事实，在连续型方程中，认为A=B；在x方向的大量方程中，去掉数量级较低的第一项，保留第二项，即边界层中的粘性作用，认为该项于压力梯度项具有相同的数量级，则有；分析y方向的动量方程，对流项与粘性扩散项具有相同的数量级，而压力梯度项却具有高得多的数量级。得到简化的近似方程组：C.E.(4.4)M.E.(4.5)(4.6)由式（4.6）可以认为，在任一过断流面上，边界层内各点的压力与其外边界上的流势压力pe相等，即，而边界层外势流区满足这里，即有边界层方程即可表示为C.E.(4.4)M.E.(4.7)对于平板绕流，Ve=U=c，则平板然刘的边界层发成可以简化为C.E.(4.
4、4)M.E.(4.8)其定界边界条件为：y=0:u=0,v=0(4.9a)y→∞:u=U(4.9b)下面采用无量纲相似性解法。i引入流函数，将未知函数的数量由两个化为一个令是连续性方程自动满足，动量方程变形为(4.10)下面用但参数群论发证明上述方程所对应的边界层流存在相似性解：令，带入上述方程：即与元方程对比，当时，与原方程有相同的形式，因此，该边界层流动存在相似性解。ii作相似性变换引入无量纲相似性变换则有式（4.10）可变换为常微分方程(4.11)边界条件相应变为：f(0)=0，f’(0)=0，f’(∞)=1(4.12)iii将初边值问题变换为初值问题令f’’(0
5、)=A，这里，A为非零常数。引入新变量则有：则式（4.11）变为(4.13)相应的边界条件变为：f”(0)=A，F”(0)=1定界条件相应变换为：F(0)=0，F’(0)=0，F”(0)=1(4.14)，iv龙格库塔解方程（1）标准四阶龙格-库塔法解微分方程组的原理对于下面的一阶常微分方程组初值问题：，其中，当时取步长，逐步做递进式运算其中分别为：（2）对三阶非线性常微分方程求解根据龙格-库塔法求解常微分方程的原理，要先将本问题中的原三阶非线性常微分方程化为一阶常微分方程组的形式然后再求解。对由方程，我们可设(1)我们再设则（1）式可变为，再由初值条件，，可推得t=0时
6、，，，于是对于该问题的一阶常微分方程组，当时，取步长为h则有：，，则原方程化简为设置一定的迭代停滞条件，这样逐步迭代就可以算出，从而得到A。同理得到A后就可以对进行化简得到：进一步用龙格-库塔法迭代求解。（3）编写matlab程序h=0.005;%步长h为0.005F=zeros(3,1000);F(3,1)=1;K=zeros(3,4);forn=1:1000,K(1,1)=h*F(2,n);K(2,1)=h*F(3,n);K(3,1)=h*(-F(1,n)*F(3,n));K(1,2)=h*(F(2,n)+K(2,1)/2);K(2,2)=h*(F(3,n)+K(3
7、,1)/2);K(3,2)=h*(-(F(1,n)+K(1,1)/2)*(F(3,n)+K(3,1)/2));K(1,3)=h*(F(2,n)+K(2,2)/2);K(2,3)=h*(F(3,n)+K(3,2)/2);K(3,3)=h*(-(F(1,n)+K(1,2)/2)*(F(3,n)+K(3,2)/2));K(1,4)=h*(F(2,n)+K(2,3));K(2,4)=h*(F(3,n)+K(3,3));K(3,4)=h*(-(F(1,n)+K(1,3))*(F(3,n)+K(3,3)));F(:,n+1)=F(:,n)+(K(:,
