《计算流体力学》结课作业解读
流体力学答案解析

流体力学答案流体力学课后答案 分析答案 解答BP1.1.1 根据阿佛迦德罗定律,在标准状态下(T = 273°K ,p = 1.013×105 Pa )一摩尔空气(28.96ɡ)含有6.022×10 23个分子。
在地球表面上70 km 高空测量得空气密度为8.75×10 -5㎏/m 3。
试估算此处 10 3μm 3体积的空气中,含多少分子数n (一般认为n <106 时,连续介质假设不再成立)答: n = 1.82×10 3提示:计算每个空气分子的质量和103μm 3体积空气的质量 解: 每个空气分子的质量为 g 1081.410022.6g 96.282323-⨯=⨯=m 设70 km 处103μm 3体积空气的质量为M g 1075.8)m 1010)(kg/m 1075.8(20318335---⨯=⨯⨯=M323201082.1g1081.4g 1075.8⨯=⨯⨯==--m M n 说明在离地面70 km 高空的稀薄大气中连续介质假设不再成立。
BP1.3.1 两无限大平行平板,保持两板的间距δ= 0.2 mm 。
板间充满锭子油,粘度为μ= 0.01Pa ⋅s ,密度为ρ= 800 kg / m 3。
若下板固定,上板以u = 0.5 m / s 的速度滑移,设油内沿板垂直方向y 的速度u (y)为线性分布,试求: (1) 锭子油运动的粘度υ; (2) 上下板的粘性切应力τ1、τ2 。
答: υ= 1.25×10 – 5 m 2/s, τ1=τ2 = 25N/m 2。
提示:用牛顿粘性定侓求解,速度梯度取平均值。
解:(1 ) /s m 1025.1kg/m800/sm kg 0.0125-3⨯===ρμν (2)沿垂直方向(y 轴)速度梯度保持常数,δμμττ/21u dydu==== (0.01Ns /m 2)(0.5m/s)/(0.2×10-3m)=25N/m 2BP1.3.2 20℃的水在两固定的平行平板间作定常层流流动。
《计算流体力学》作业答案

计算流体力学作业答案问题1:什么是计算流体力学?计算流体力学(Computational Fluid Dynamics,简称CFD)是研究流体力学问题的一种方法,它使用数值方法对流体流动进行数值模拟和计算。
主要包括求解流体运动的方程组,通过空间离散和时间积分等计算方法,得到流体在给定条件下的运动和相应的物理量。
问题2:CFD的应用领域有哪些?CFD的应用领域非常广泛,包括但不限于以下几个方面:1.汽车工业:CFD可以用于汽车流场的模拟和优化,包括空气动力学性能和燃烧过程等。
2.航空航天工业:CFD可以用于飞机、火箭等流体动力学性能的预测和优化,包括机身、机翼的设计和改进等。
3.能源领域:CFD可以用于燃烧、热交换等能源领域的流体力学问题求解和优化。
4.管道流动:CFD可以用于石油、化工等行业的管道流动模拟和流体输送优化。
5.空气净化:CFD可以用于大气污染物的传输和分布模拟,以及空气净化设备的设计和改进。
6.生物医药:CFD可以用于生物流体输送和生物反应过程的模拟和分析,包括血液流动、药物输送等。
问题3:CFD的数值方法有哪些?CFD的数值方法一般包括以下几种:1.有限差分法(Finite Difference Method,FDM):将模拟区域划分为网格,并在网格上离散化流体运动的方程组,利用有限差分近似求解。
2.有限体积法(Finite Volume Method,FVM):将模拟区域划分为有限体积单元,通过对流体流量和通量的控制方程进行离散化,求解离散化方程组。
3.有限元法(Finite Element Method,FEM):将模拟区域划分为有限元网格,通过对流体运动方程进行弱形式的变分推导,将流动问题转化为求解线性方程组。
4.谱方法(Spectral Method):采用谱方法可以对流体运动方程进行高精度的空间离散,通常基于傅里叶变换或者基函数展开的方式进行求解。
5.计算网格方法(Meshless Methods):不依赖网格的数值方法,主要包括粒子方法(Particle Methods)、网格自适应方法(Gridless Method)等。
流体力学研究生就业(3篇)

流体力学研究生就业(3篇)以下是网友分享的关于流体力学研究生就业的资料3篇,希望对您有所帮助,就爱阅读感谢您的支持。
篇一2013~2014学年第1学期12级研究生《计算流体力学》结课作业适用专业:供热供燃气通风及空调工程一、试论述有限差分法和有限容积法的基本思想、技术要点和求解步骤。
(不少于4千字)。
(30分)二、试说明有限差格式中耗散与色散的含义(15分)三、用分离变量法求解热方程(15分)∂u∂2u= 0≤x≤1 ∂t∂x2边界条件:u(t,0)=0u(t,1)=0初始条件:u(0,x)=sin(2πx) 四、试确定如下关于∂u∂y的一个差分近似格式的截断误差(均匀网格)。
(25分)∂u-3ui,j+4ui,j+1-ui,j+2 ≈∂y2∆y五、写出下列波方程的迎风差分格式,并分析其稳定性条件。
(25分)ut+c(ux+uy)=0 0≤x≤1篇二《流体力学》科目研究生考试大纲一、流体力学课程的性质与任务-中国在职研究生招生网官网―流体力学‖是工程类本科专业(化工、机械、能源及相关专业)的一门必修学科基础课。
它是研究流体受力及其宏观运动规律的一门学科,既有基础学科的性质,又有鲜明的应用学科特点,而工程流体力学则更侧重于后者。
热能与动力工程专业的流体力学课程的任务是使学生掌握流体力学的基本原理、基本方法及其在工程实际问题中的应用,从而为分析研究流体流动规律及其相关传递过程,以及设计开发及应用新型高效热能与动力工程装备奠定必备的基础。
二、课程教学基本要求学生应达到如下要求:掌握流体的力学基础及运动学的基本概念,能够根据流体力学的基本原理及分析方法解决静止流体、流体流动的部分工程实际问题。
三、理论内容和基本要求(一)绪论1.流体力学学习的相关问题2.流体的连续介质模型3.流体的主要物理性质4.牛顿流体和非牛顿流体基本要求:1.理解连续介质与理想流体的概念和在流体力学研究中的意义;2.理解流体的主要物理力学性质,重点掌握流体粘滞性、牛顿内摩擦定律及其适用条件;3.掌握物理量的基本量纲、基本单位及导出量的单位。
计算流体力学大作业

计算液体力学基础及应用课程期末作业-----程序调试最终版学号:134212059 姓名:徐影ContentsCFD模型示意图一、拟一维喷管理论解求解二、拟一维喷管的CFD求解三、理论值与CFD解的对比CFD模型示意图两圆弧直径为10米,喉部直径为0.59米,长为3米clear all;I=imread('xuying.png'); imshow(I)一、拟一维喷管理论解求解喷管内马赫数的变化公依赖于面积比A/A0,所以可以将Ma作为x的函数1.2.采用隐函数绘图给出理论的马赫数解gamma=1.4;h0=59/100;% 取学生学号后两位数的十分之一作喉部直径syms x Ma A_x y;% xz为x坐标,Ma为马赫数A_x=((10.59-2*sqrt(25-(x-1.5)^2))/0.59)^2;% A_x为面积系数figure('Color',[1 1 1]);set(gcf,'position',[0,0,1.5*468,468]);plot_Ma=A_x^2-(2/(gamma+1)+(gamma-1)/(gamma+1)*y^2)^((gamma+1)/(gamma-1))/y^2;subplot(1,2,1);gca=ezplot(plot_Ma,[0,3]);xlabel('x');ylabel('马赫数');title('采用隐函数求解的马赫数结果');grid on; % 得到两条曲线,由递增规律选取上升曲线段,从该曲线上得到一系列点的坐标为[x0,Ma0]load tk.mat;x_0=tk(:,1);Ma_0=tk(:,2);% 这里load的数据采用某算法从上面出的图取点拟合得到,用到polyval和polyfit函数subplot(1,2,2);plot(x_0,Ma_0);xlabel('x');ylabel('马赫数');title('马赫数的理论解');grid on;求出马赫数后,压力、密度、温度的变化都是Ma的函数,求出理论值并绘图1.2.3.p_0=(1+(gamma-1)/2*Ma_0.^2).^(-gamma/(gamma-1));rho_0=(1+(gamma-1)/2*Ma_0.^2).^(-1/(gamma-1));t_0=(1+(gamma-1)/2*Ma_0.^2).^-1;figure('Color',[1 1 1]);set(gcf,'position',[0,0,1.5*468,1.5*468]);subplot(3,1,1);plot(x_0,p_0);title('压力比理论值');xlabel('x');ylabel('p');grid on; subplot(3,1,2);plot(x_0,rho_0);title('密度比理论值');xlabel('x');ylabel('rho');grid on; subplot(3,1,3);plot(x_0,t_0);title('温度比理论值');xlabel('x');ylabel('T');grid on;二、拟一维喷管的CFD求解clear all;L=3;N=31;dx=L/(N-1);x=linspace(0,L,N);C=0.5;n=2000;student_num=59;A=((10+student_num/100-2*((25-((x-1.5).^2))).^0.5)/(student_num/100)).^2;%面积比A/A_0与x坐标的关系第一步,密度比、温度比、速度比的初始条件设定1.2.3.Rou=1-0.3146*x;rhobi=zeros(1,n);T=1-0.2314*x;V=(0.1+1.09*x).*sqrt(T);P_rou_t=zeros(size(Rou));P_v_t=zeros(size(Rou));P_T_t=zeros(size(Rou));P_rou_t_2=zeros(size(Rou));P_v_t_2=zeros(size(Rou));P_T_t_2=zeros(size(Rou));第二步,预估步第三步,并求Δt,求rou, V, T的预测量1.2.3.第四步,修正步第五步,求平均时间导数1.2.3.最后,得到t+Delta t时刻流动参数的修正值为1.2.3.第七步,边界条件处理for j=1:ntemp=Rou(16);% 第二步,预估步for i=2:30P_rou_t(i)=-V(i)*((Rou(i+1)-Rou(i))/dx)-Rou(i)*((V(i+1)-V(i))/dx)-Rou(i)*V(i)*((log(A(i+1))-log(A(i)))/dx);P_v_t(i)=-V(i)*((V(i+1)-V(i))/dx)-((T(i+1)-T(i))/dx+((Rou(i+1)-Rou(i))/dx)*T(i)/Rou(i))*1/1.4;P_T_t(i)=-V(i)*((T(i+1)-T(i))/dx)-0.4*T(i)*(((V(i+1)-V(i))/dx)+V(i)*((log(A(i+1))-log(A(i)))/dx));end% 第三步,并求Δt,求rou, V, T的预测量dt=C*(dx./(V(2:30)+sqrt(T(2:30))));dt=min(dt);Rou1(2:30)=Rou(2:30)+P_rou_t(2:30).*dt;V1(2:30)=V(2:30)+P_v_t(2:30).*dt;T1(2:30)=T(2:30)+P_T_t(2:30).*dt;V1(1)=V(1);T1(1)=T(1);Rou1(1)=Rou(1);% 第四步,修正步%for i=2:30P_rou_t_2(i)=-V1(i)*((Rou1(i)-Rou1(i-1))/dx)-Rou1(i)*((V1(i)-V1(i-1))/dx)-Rou1(i)*V1(i)*((log(A(i))-log(A(i-1)))/dx); P_v_t_2(i)=-V1(i)*((V1(i)-V1(i-1))/dx)-((T1(i)-T1(i-1))/dx+((Rou1(i)-Rou1(i-1))/dx)*T1(i)/Rou1(i))*1/1.4;P_T_t_2(i)=-V1(i)*((T1(i)-T1(i-1))/dx)-0.4*T1(i)*(((V1(i)-V1(i-1))/dx)+V1(i)*((log(A(i))-log(A(i-1)))/dx));end% 第五步,求平均时间导数P_rou_av=(P_rou_t+P_rou_t_2)/2;P_v_av=(P_v_t+P_v_t_2)/2;P_T_av=(P_T_t+P_T_t_2)/2;% 最后,得到t+Delta t时刻流动参数的修正值为Rou(2:30)=Rou(2:30)+P_rou_av(2:30).*dt;T(2:30)=T(2:30)+P_T_av(2:30).*dt;V(2:30)=V(2:30)+P_v_av(2:30).*dt;P(2:30)=Rou(2:30).*T(2:30);% 第七步,边界条件处理V(1)=2*V(2)-V(3);V(31)=2*V(30)-V(29);Rou(31)=2*Rou(30)-Rou(29);T(31)=2*T(30)-T(29);p=Rou.*T;Ma=V./sqrt(T);rhobi(j)=abs((temp-Rou(16))/temp); % 计算后一次时间步与前一时间步之间的密度比的变化情况,以此检验CFD过程收敛性质end最终结果的绘图figure('Color',[1 1 1]);set(gcf,'position',[0,0,1.2*468,1.5*468]);subplot(3,1,1);plot(1:n,rhobi);xlabel('x');ylabel('Ma');title('相对密度比');grid on;% 密度比收敛情况绘图subplot(3,1,2);plot(x,Ma);title('喷管内马赫数分布');xlabel('x');ylabel('Ma');grid on;% 马赫数CFD值绘图subplot(3,1,3);plot(x,p);title('喷管内压力分布');xlabel('x');ylabel('p');grid on; % 压力分布CFD值绘图shu=[x;A;Ma;V;T;p;Rou];显示各参量最终计算结果fprintf('%6s\t%12s\t%12s\t%12s\t%12s\t%12s\t%12s\r\n','x','A/A_0','Ma','v/v_0','T/T_0','p/p_0','rho')% 依次显示坐标点、形状参数、马赫数、速度、温度、压力的结果fprintf('%6.1f\t%12.4f\t%12.4f\t%12.4f\t%12.4f\t%12.4f\t%12.4f\r\n',shu)x A/A_0 Ma v/v_0 T/T_0 p/p_0 rho0.0 3.1709 0.1859 0.1859 1.0000 1.0000 1.00000.1 2.8156 0.2124 0.2121 0.9975 0.9915 0.99390.2 2.5056 0.2389 0.2383 0.9956 0.9847 0.98900.3 2.2361 0.2711 0.2700 0.9922 0.9728 0.98050.4 2.0030 0.3056 0.3038 0.9885 0.9602 0.97140.5 1.8022 0.3451 0.3422 0.9834 0.9433 0.95910.6 1.6303 0.3882 0.3838 0.9775 0.9234 0.94470.7 1.4844 0.4364 0.4298 0.9700 0.8989 0.92670.8 1.3617 0.4891 0.4794 0.9611 0.8701 0.90540.9 1.2600 0.5469 0.5331 0.9502 0.8362 0.88001.0 1.1771 0.6096 0.5903 0.9374 0.7974 0.85071.1 1.1116 0.6776 0.6508 0.9224 0.7536 0.81701.2 1.0620 0.7507 0.7142 0.9051 0.7053 0.77921.3 1.0273 0.8289 0.7800 0.8855 0.6532 0.73761.4 1.0068 0.9119 0.8475 0.8636 0.5982 0.69271.5 1.0000 0.9998 0.9160 0.8394 0.5416 0.64521.6 1.0068 1.0921 0.9849 0.8132 0.4847 0.59601.7 1.0273 1.1887 1.0534 0.7853 0.4288 0.54611.8 1.0620 1.2893 1.1210 0.7559 0.3753 0.49641.9 1.1116 1.3934 1.1869 0.7255 0.3250 0.44802.0 1.1771 1.5009 1.2507 0.6943 0.2788 0.40152.1 1.2600 1.6113 1.3119 0.6629 0.2371 0.35762.2 1.3617 1.7245 1.3705 0.6315 0.2001 0.31682.3 1.4844 1.8398 1.4258 0.6006 0.1678 0.27952.4 1.6303 1.9576 1.4782 0.5702 0.1400 0.24552.5 1.8022 2.0764 1.5269 0.5408 0.1163 0.21512.6 2.0030 2.1983 1.5732 0.5122 0.0962 0.1879。
计算流体力学实验报告

一、实验目的1. 了解计算流体力学的基本原理和方法;2. 掌握计算流体力学软件的使用方法;3. 通过实验验证计算流体力学在工程中的应用。
二、实验原理计算流体力学(Computational Fluid Dynamics,简称CFD)是一种利用数值方法求解流体运动和传热问题的学科。
其基本原理是利用数值方法将连续的物理问题离散化,将其转化为求解偏微分方程组的问题。
在计算流体力学中,常用的数值方法有有限差分法、有限元法和有限体积法。
本实验采用有限体积法进行流体运动的数值模拟。
有限体积法将计算区域划分为若干个控制体,在每个控制体上应用守恒定律,将连续的偏微分方程转化为离散的代数方程组。
通过求解这些代数方程组,可以得到流体在各个控制体内的速度、压力和温度等参数。
三、实验内容1. 实验一:二维不可压缩流体的稳态流动模拟(1)实验目的:通过模拟二维不可压缩流体的稳态流动,验证计算流体力学在流体运动模拟中的应用。
(2)实验步骤:① 建立二维流场模型,包括进口、出口、壁面和障碍物等;② 划分计算区域,选择合适的网格划分方法;③ 设置边界条件和初始条件;④ 选择合适的数值方法和湍流模型;⑤ 运行计算流体力学软件,得到流场参数;⑥ 分析结果,绘制流线图、速度矢量图等。
(3)实验结果与分析:通过模拟二维不可压缩流体的稳态流动,得到流场参数,并绘制流线图、速度矢量图等。
根据实验结果,可以分析流场特征,验证计算流体力学在流体运动模拟中的应用。
2. 实验二:三维不可压缩流体的瞬态流动模拟(1)实验目的:通过模拟三维不可压缩流体的瞬态流动,验证计算流体力学在流体运动模拟中的应用。
(2)实验步骤:① 建立三维流场模型,包括进口、出口、壁面和障碍物等;② 划分计算区域,选择合适的网格划分方法;③ 设置边界条件和初始条件;④ 选择合适的数值方法和湍流模型;⑤ 运行计算流体力学软件,得到流场参数;⑥ 分析结果,绘制流线图、速度矢量图等。
计算流体力学课程大作业

《计算流体力学》课程大作业——基于涡量-流函数法的不可压缩方腔驱动流问题数值模拟张伊哲 航博1011、 引言和综述2、 问题的提出,怎样使用涡量-流函数方法建立差分格式3、 程序说明4、 计算结果和讨论5、 结论1引言虽然不可压缩流动的控制方程从形式上看更为简单,但实际上,目前不可压缩流动的数值方法远远不如可压缩流动的数值方法成熟。
考虑不可压缩流动的N-S 方程:01()P t νρ∇⋅=⎧⎪∂⎨+∇⋅=-∇+∆⎪∂⎩U UUU f U (1.1)其中ν是运动粘性系数,认为是常数。
将方程组写成无量纲的形式:01()Re P t∇⋅=⎧⎪∂⎨+∇⋅=-∇+∆⎪∂⎩U UUU f U (1.2) 其中Re 是雷诺数。
从数学角度看,不可压缩流动的控制方程中不含有密度对时间的偏导数项,方程表现出椭圆-抛物组合型的特点;从物理意义上看,在不可压缩流动中,压力这一物理量的波动具有无穷大的传播速度,它瞬间传遍全场,以使不可压缩条件在任何时间、任何位置满足,这就是椭圆型方程的物理意义。
这就造成不可压缩的N-S 方程不能使用比较成熟的发展型...偏微分方程的数值求解理论和方法。
如果将动量方程和连续性方程完全耦合求解,即使使用显示的离散格式,也将会得到一个刚性很强的、庞大的稀疏线性方程组,计算量巨大,更重要的问题是不易收敛。
因此,实际应用中,通常都必须将连续方程和动量方程在一定程度上解耦。
目前,求解不可压缩流动的方法主要有涡量-流函数法,SIMPLE 法及其衍生的改进方法,有限元法,谱方法等,这些方法各有优缺点。
其中涡量-流函数法是解决二维不可压缩流动的有效方法。
作者本学期学习了研究生计算流体课程,为了熟悉计算流体的基本方法,选择使用涡量-流函数法计算不可压缩方腔驱动流问题,并且对于不同雷诺数下的解进行比较和分析,得出一些结论。
本文接下来的内容安排为:第2节提出不可压缩方腔驱动流问题,并分析该问题怎样使用涡量-流函数方法建立差分格式、选择边界条件。
计算流体力学大作业

1 提出问题[问题描述]Sod 激波管问题是典型的一类Riemann 问题。
如图所示,一管道左侧为高温高压气体,右侧为低温低压气体,中间用薄膜隔开。
t=0 时刻,突然撤去薄膜,试分析其他的运动。
Sod 模型问题:在一维激波管的左侧初始分布为:0 ,1 ,1111===u p ρ,右侧分布为:0 ,1.0 ,125.0222===u p ρ,两种状态之间有一隔膜位于5.0=x 处。
隔膜突然去掉,试给出在14.0=t 时刻Euler 方程的准确解,并给出在区间10≤≤x 这一时刻u p , ,ρ的分布图。
2 一维Euler 方程组分析可知,一维激波管流体流动符合一维Euler 方程,具体方程如下: 矢量方程:0U ft x∂∂+=∂∂ (0.1)分量方程:连续性方程、动量方程和能量方程分别是:222,,p u ρ()()()()2000u tx u u pt x x u E p E tx ρρρρ∂⎧∂+=⎪∂∂⎪⎪∂∂∂⎪++=⎨∂∂∂⎪⎪∂+⎡⎤∂⎣⎦+=⎪∂∂⎪⎩ (0.2)其中 22v u E c T ρ⎛⎫=+ ⎪⎝⎭对于完全气体,在量纲为一的形式下,状态方程为:()2p T Ma ργ∞=(0.3)在量纲为一的定义下,定容热容v c 为:()211v c Ma γγ∞=- (0.4)联立(1.2),(1.3),(1.4)消去温度T 和定容比热v c ,得到气体压力公式为:()2112p E u γρ⎛⎫=-- ⎪⎝⎭(0.5)上式中γ为气体常数,对于理想气体4.1=γ。
3 Euler 方程组的离散3.1 Jacibian 矩阵特征值的分裂Jacibian 矩阵A 的三个特征值分别是123;;u u c u c λλλ==+=-,依据如下算法将其分裂成正负特征值:()12222k k k λλελ±±+=(0.6)3.2 流通矢量的分裂这里对流通矢量的分裂选用Steger-Warming 分裂法,分裂后的流通矢量为()()()()()()()12312322232121212122f u u c u c u u c u c w γλλλργλλλγλλγλ⎛⎫⎪-++ ⎪=-+-++ ⎪ ⎪ ⎪-+-+++ ⎪⎝⎭+++++++++++(0.7)()()()()()()()12312322232121212122f u u c u c u u c u c w γλλλργλλλγλλγλ⎛⎫⎪-++ ⎪=-+-++ ⎪ ⎪ ⎪-+-+++ ⎪⎝⎭-----------(0.8)其中:()()()223321c w γλλγ±±±-+=- c 为量纲为一的声速:22Tc Ma ∞=(0.9)联立(1.3),(1.9)式,消去来流马赫数得:ργp c =3.3 一阶迎风显示格式离散Euler 方程组 10n n i i x i x i U U f f t xδδ+-++--++=∆∆ (0.10)得到()()n+1nj j 11U =U j j j j t f f f f x++---+∆⎡⎤--+-⎣⎦∆ 算法如下:① 已知初始时刻t=0的速度、压力及密度分布000,,j j j u P ρ,则可得到特征值分裂值0k λ±,从而求出流通矢量0j f ±;② 应用一阶迎风显示格式可以计算出1t t =∆时刻的组合变量1j U ,从而得到1t t =∆时刻的速度、压力及密度分布111,,j j j u P ρ;③ 利用1t t =∆时刻的速度、压力及密度分布111,,j j j u P ρ可得特征值分裂值1k λ±,从而求出流通矢量1j f ±;④ 按照步骤2的方法即可得到2t t =∆时刻的速度、压力及密度分布222,,j j j u P ρ;⑤ 循环以上过程即可得到()1t n t =+∆时刻的速度、压力及密度分布n+1n+1n+1,,j j j u P ρ。
计算流体力学课后题作业

课后习题第一章1.计算流体动力学的基本任务是什么计算流体动力学是通过计算机数值计算和图像显示,对包含有流体流动和热传导等相关物理现象的系统所做的分析。
2.什么叫控制方程?常用的控制方程有哪几个?各用在什么场合?流体流动要受物理守恒定律的支配,基本的守恒定律包括:质量守恒定律、动量守恒定律、能量守恒定律。
如果流动包含有不同组分的混合或相互作用,系统还要遵守组分守恒定律。
如果流动处于湍流状态,系统还要遵守附加的湍流输运方程。
控制方程是这些守恒定律的数学描述。
常用的控制方程有质量守恒方程、动量守恒方程、能量守恒方程、组分质量守恒方程。
质量守恒方程和动量守恒方程任何流动问题都必须满足,能量守恒定律是包含有热交换的流动系统必须满足的基本定律。
组分质量守恒方程,在一个特定的系统中,可能存在质的交换,或者存在多种化学组分,每种组分都需要遵守组分质量守恒定律。
4.研究控制方程通用形式的意义何在?请分析控制方程通用形式中各项的意义。
建立控制方程通用形式是为了便于对各控制方程进行分析,并用同一程序对各控制方程进行求解。
各项依次为瞬态项、对流项、扩散项、源项。
6.CFD商用软件与用户自行设计的CFD程序相比,各有何优势?常用的商用CFD软件有哪些?特点如何?由于CFD的复杂性及计算机软硬件条件的多样性,用户各自的应用程序往往缺乏通用性。
CFD商用软件的特点是功能比较全面、适用性强。
具有比较易用的前后处理系统和其他CAD及CFD软件的接口能力,便于用户快速完成造型、网格划分等工作。
具有比较完备的容错机制和操作界面,稳定性高。
可在多种计算机、多种操作系统,包括并行环境下运行。
常用的商用CFD软件有PHOENICS、CFX、SRAR-CD、FIDAP、FLUENT。
PHOENICS除了通用CFD软件应该拥有的功能外,PHOENICS软件有自己独特的功能:开放性、CAD接口、运动物体功能、多种模型选择、双重算法选择、多模块选择。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2012~2013学年第1学期12级研究生《计算流体力学》结课作业适用专业:供热供燃气通风及空调工程一、结合某一具体学科,阐述纯理论方法、实验方法及数值方法在科学研究中的各自优缺点,在此基础上论述数值模拟方法的发展前景。
(不少于4千字)。
流体力学是力学的一个重要分支, 是研究流体(液体和气体)的力学运动规律及其应用的学科, 主要研究在各种力的作用下,流体本身的静止状态和运动状态特征,以及流体和相邻固体界面有相对运动时的相互作用和流动规律。
在人们的生活和生产活动中随时随地都可遇到流体,流体力学与人类的日常生活和生产事业密切相关。
按其研究内容的侧重点不同,分为理论流体力学和工程流体力学。
其中理论流体力学主要采用严密的数学推理方法,力求准确性和严密性,工程流体力学侧重于解决工程实际中出现的问题,而不追求数学上的严密性。
当然由于流体力学研究的复杂性,在一定程度上,两种方法都必须借助于实验研究,得出经验或半经验的公式。
在实际工程的诸多领域流体力学都起着十分重要的作用。
如气象、水利的研究,船舶、飞行器、叶轮机械和核电站的设计及其运行,可燃气体或炸药的爆炸,都广泛地用到流体力学知识。
许多现代科学技术所关心的问题既受流体力学的指导,同时也促进了流体力学自身的不断发展。
1950年后,计算机的发展给予流体力学以极大的推动作用。
目前,解决流体力学问题的方法主要有实验方法、理论分析方法和数值方法三种。
实验方法同物理学、化学等学科一样,流体力学的研究离不开实验,尤其是对新的流体运动现象的研究。
实验能显示运动特点及其主要趋势,有助于形成概念,检验理论的正确性。
二百年来流体力学发展史中每一项重大进展都离不开实验。
流体力学实验研究方法有实物实验、比拟研究和模型研究三类:实物实验是用仪器实测原型系统的流动参数,适用于较小的原型;比拟实验是利用电场和磁场来模拟流场,实施起来限制条件较多;模型研究是实验流体力学最常用的研究方法。
实验研究的一般过程是:在相似理论的指导下建立实验模型,用流体测量技术测量流动参数,处理和分析实验数据。
建立实验模型要求模型与原型满足相似理论,即满足两个流场相似。
流体力学中两个流场相似要求: 几何相似、运动相似、动力相似、边界条件、初始条件相似。
两个流场动力相似,则两个流场所有的动力相似准则应分别相等。
但要做到两个独立的动力相似准则同时分别与原型的同名准则相等是不可能的,所以只能部分相似,即近似模型实验。
模拟实验在流体力学中占有重要地位。
根据模型实验所得的数据可以用像换算单位制那样的简单算法求出原型的数据。
实验方法有诸多优点:实验方法可靠性高,能反映工程中的实际流动规律,发现新现象,检验理论结果等;工程实际中,由于控制方程多为非线性方程,大多问题无法得到理论解析结果,而必须借助于实验的方法,尤其是对于目前尚未有合适数学模型的复杂湍流流动、某些非牛顿流体的流动、多相流等问题,实验测试则是唯一的研究方法。
但实验方法受到模型尺寸、流动扰动、人身安全和测量精度的限制,有时可能通过实验无法得到结果;另外实验中还会遇到经费投入不足,人力、物力的巨大耗费及周期长等诸多困难。
理论分析方法理论分析(理论研究方法)是根据流体运动的普遍规律如质量守恒、动量守恒、能量守恒等,利用数学分析的手段,研究流体的运动,解释已知的现象,预测可能发生的结果。
理论分析的一般过程是:建立力学模型,用物理学基本定律推导流体力学数学方程,用数学方法求解方程,检验和解释求解结果。
理论研究方法的关键在于提出理论模型,并能运用数学方法求出理论结果,达到揭示液体运动规律的目的。
流体力学中最常用的基本模型有:连续介质、牛顿流体、不可压缩流体、理想流体、平面流动等。
对这样的理论模型,根据机械运动的普遍规律,用数学语言将质量守恒、动量守恒、能量守恒等定律表达出来,从而得到连续性方程、动量方程和能量方程。
此外,还要加上某些联系流动参量的关系式(例如状态方程),或者其他方程,构成流体力学基本方程组。
将原来的具体流动问题转化为数学问题,在相应的边界条件和初始条件下求解。
求出方程组的解后,可结合具体流动,解释这些解的物理含义和流动机理。
理论分析优点:能揭示流动的内在规律,具有普遍适用性,成本最低,结果最理想,影响因素表达清楚。
但理论分析方法局限于非常简单的问题。
数值方法数值研究的一般过程是:对流体力学数学方程作简化和数值离散化,编制程序作数值计算,将计算结果与实验结果比较。
在流体力学理论研究和工程应用中,描述流体运动的数学方程是非线性偏微分方程组,只对极少数的简化模型可以通过数学方法,获得理论分析解,多数情况下,只能通过数值计算的途径进行求解。
这里说的“数值计算”,是指利用高速电子计算机,对描述流体力学具体问题的偏微分方程初边值问题进行离散化计算,从而获得流动区域中离散点上的流体物理量的求解方法。
这种通过数值计算获得流动区域中离散点的数值解的方法,通常称为流体力学数值解法,也可称之为计算流体力学。
随着高速电子计算机的发展与普及,数值方法越来越受到重视,已成为流体力学理论研究和工程应用的重要手段。
计算流体力学是以计算机为工具、以流体力学的基本方程(纳维-斯托克斯方程)为理论依据,采用离散化的数值方法对流体力学问题进行数值模拟和分析的流体力学分支学科。
常用的方法有:有限差分法、有限元法、有限体积法、边界元法、谱分析法等。
计算的内容包括:飞机、汽车、河道、桥梁、涡轮机等流场计算;湍流、流动稳定性、非线性流动等数值模拟。
大型工程计算软件已成为研究工程流动问题的有力武器。
数值方法的优点是能计算理论分析方法无法求解的数学方程,比实验方法省时省钱,但毕竟是一种近似解方法,适用范围受数学模型的正确性和计算机的性能所限制。
三种方法各有优缺点,我们应取长补短,互为补充。
流体力学的研究不仅需要深厚的理论基础,而且需要很强的动手能力。
学习流体力学应注重理论与实践结合,理论分析、实验研究和数值计算并重。
数值模拟方法的发展前景任何流体运动的动力学特征都是由质量守恒、动量守恒和能量守恒定律所确定的,这些基本定律可以由流体流动的控制方程组来描述。
而这些控制方程大多是一些极其复杂的偏微分方程,在考虑粘性作用时更是如此,如果不靠计算机,就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。
20世纪30~40年代,对于复杂而又特别重要的流体力学问题,曾组织过人力用几个月甚至几年的时间做数值计算,比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。
数学的发展,计算机的不断进步,以及流体力学各种计算方法的发明,使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性,此时,数值方法显现出了极大地优越性。
这又促进了流体力学计算方法的发展,并形成了“计算流体力学”。
计算流体力学(Computational Fluid Dynamics,简称CFD)是21世纪流体力学领域的重要技术之一,使用数值方法在计算机中对流体力学的控制方程进行求解,从而可预测流场的流动。
流体力学的运动方程是极其复杂的非线性偏微分方程,具有各种不同的类型,而且往往还是混合型的。
计算流体力学在很大程度上就是针对不同性质的偏微分方程采用和发展相应的数值解方法。
经过40年来的发展,计算流体力学己经成为一种有力的数值实验与设计手段,在许多工业领域如航天航空、汽车、船舶等部门解决了大量的工程设计实际问题,其中在航天航空领域所取得的成绩尤为显著。
现在人们已经可以利用计算流体力学方法来设计飞机的外形,确定其气动载荷,从而有效地提高了设计效率,减少了风洞试验次数,大大地降低了设计成本。
此外,计算流体力学也己经大量应用于大气、生态环境、车辆工程、船舶工程、传热以及工业中的化学反应等各个领域,显示了计算流体力学强大的生命力。
随着计算机技术、网络技术、计算方法和后处理技术的迅速发展,利用计算流体力学解决流动问题的能力越来越高,现在许多复杂的流动问题可以通过数值计算手段进行分析并给出相应的结果。
计算流体力学也己经发展成为以数值手段求解流体力学物理模型、分析其流动机理为主线。
目前计算流体力学主要向二个方向发展:一方面是研究流动非定常稳定性以及湍流流动机理,开展高精度、高分辩率的计算方法和并行算法等的流动机理与算法研究;另一方面是将计算流体力学直接应用于模拟各种实际流动,解决工业生产中的各种问题。
随着科技的进步,许多关于数值模拟方法的商业软件产生了。
计算流体力学商业软件最早出现于上世纪八十年代初,目前已经在工业领域和学术研究领域发挥着积极的作用。
这些软件的使用减少了计算流体力学研究和开发人员的工作量,降低了对研究人员计算机知识的要求,从而使研究者可以把精力集中在对计算流体力学本质问题的研究和技术开发上。
流体力学和其他学科一样,是通过理论分析和实验研究两种手段发展起来的。
很早就已有理论流体力学和实验流体力学两大分支。
解决流体力学问题时、实验方法、理论分析方法和数值方法是相辅相成的。
实验需要理论指导,才能从分散的、表面上无联系的现象和实验数据中得出规律性的结论;而理论分析和数值模拟也要依靠实验方法给出物理图形或数据,以建立流动的力学模型和数学模式;最后,还须依靠实验来检验这些模型的完善程度。
此外,实际流动往往异常复杂,理论分析和数值计算会遇到巨大的数学和计算方面的困难,得不到具体结果,只能通过实验方法进行研究。
任何一种方法都有它自身的局限性,而数值方法在弥补理论分析的不足之处有着极大的优越性。
理论分析是用数学方法求出问题的定量结果,鉴于流体力学控制方程组的特点,能用这种方法求出结果的问题毕竟是少数,计算流体力学正是为弥补分析方法的不足而发展起来的。
数值模拟方法的优点显而易见,在解决工程问题中已得到了很好的证明,特别是近些年关于数值模拟方法的软件日趋成熟,使得数值模拟方法的应用更为方便准确。
数值模拟已经成为人类改造世界的第三种手段,今后将会得到广泛应用与长足发展。
二、概述有限差分方法的基本思想、技术要点及应用步骤等(不少于5千字)。
物理学和其他学科领域的许多问题在被分析研究之后, 往往可以归结为常微分方程或偏微分方程的求解问题。
一般说来,处理一个特定的物理问题,除了需要知道它满足的数学方程外,还应当同时知道这个问题的定解条件,然后才能设计出行之有效的计算方法来求解。
有限差分法以变量离散取值后对应的函数值来近似微分方程中独立变量的连续取值。
在有限差分方法中,我们放弃了微分方程中独立变量可以取连续值的特征,而关注独立变量离散取值后对应的函数值。
但是从原则上说,这种方法仍然可以达到任意满意的计算精度。
因为方程的连续数值解可以通过减小独立变量离散取值的间格,或者通过离散点上的函数值插值计算来近似得到。