CFD考试题答案参考

1、.计算流体动力学及其应用 习题1、扩散 (质量、动量、能量)现象的方程和边界条件为：u2u0 t T0 x Ltx2初始条件： u(x,0)=0 (除x=0.5外)；u(0.5,0)=0.1边界条件： u(0,t)=1 ，u(L ,t)=1 (除 t=0外)；u(x,0)=0设 =1,x=0.1, L =1, T=2, I =10；请在 t=0.005 ， t=0.01两种情况下各进行 10个时间步的数值计算，并找出某些规律。提示：采用时间前差空间中差格式。看两种情况是否稳定？t n t u0u1u2 u3u4u5u6 u7u8u9 u100000000.1000000.005110.0111

2、0.015110.02110.025110.0311.2、已知：d 2 u2u 4x22x40x1dx2边界条件： x=0,du1 ； x=1,u= -1dx利用 FDM 和FVM ，求各节点的 u值，并与精确解进行比较。提示：参考 P6-8 Neumann 边界条件3、用一般方法 (前向差分、后向差分、中心差分)导出一阶微分u具有三 (六)阶精度的差分格式。p20x4、采用 Von Neumann 方法判断下式的稳定性。ui n 1ui n 1(ui 1n2ui nui 1n )O( t 2 , x2 )2tx2P37-39 可采用反证法证明其无条件不稳定。.5、已知：2uf ( x, y)

3、精确解为：u2x2 y 2应用第一类边界条件，有：u1, u2 ,u3,u6 ,u90u48, u88f58,f820采用有限体积法求 2 2单位网格的数值解。 p1126、利用 CFD 软件 (如 FLUENT 或CFX 或STAR-CD 等)计算一个实例。要求：写出实例的具体内容、边界条件、初始条件并打印出模型、网格及计算结果。1. 采用时间前差，空间中差uin 1uind (uin 12uinuin 1)t0.005tU0U1U2U3U4U5U6U7U8U9U100000000.1000000.00510000.0500.0500010.0110.500.02500.0500.02500

4、.510.01510.50.262500.037500.037500.26250.510.0210.63120.250.1500.037500.150.250.631210.02510.6250.39060.1250.093800.09380.1250.39060.62510.0310.69530.3750.24220.06250.09380.06250.24220.3750.695310.003510.68750.46880.21880.1680.06250.1680.21880.46880.687510.0410.73440.45310.31840.14060.1680.14060.318

5、40.45310.734410.04510.72660.52640.29690.24320.14060.24320.29690.52640.726610.0510.76320.51170.38480.21880.24320.21880.38480.51170.76321t=0.01tU0U1U2U3U4U5U6U7U8U9U1000.000.000.000.000.000.100.000.000.000.000.000.011.000.000.000.000.10-0.100.100.000.000.001.000.021.001.000.000.10-0.200.30-0.200.100.0

6、01.001.00.0.031.000.001.10-0.300.60-0.700.60-0.301.100.001.000.041.002.10-1.402.00-1.601.90-1.602.00-1.402.101.000.051.00-2.505.50-5.005.50-5.105.50-5.005.50-2.501.000.061.009.00-13.0016.00-15.6016.10-15.6016.00-13.009.001.000.071.00-21.0038.00-44.6047.70-47.3047.70-44.6038.00-21.001.000.081.0060.00

7、-103.60130.30-139.60142.70-139.60130.30-103.6060.001.000.091.00-162.60293.90-373.50412.60-421.90412.60-373.50293.90-162.601.000.101.00457.50-830.001080.00-1208.001247.10-1208.001080.00-830.00457.501.00从时间步长为 t=0.005 情况看，速度逐步衰减，是收敛的；时间步长为t=0.01 情况下，速度存在突变，是发散的。由此可以看出选择时间步长对计算结果收敛情况起决定作用。分析：采用 Von Neu

8、mann方法，分别计算 T=0.005s和 T=0.01s时的扩散因子 d，可发现 T=0.005s时 0 d1/ 2 ，计算结果稳定，而 T=0.01s时 d1/ 2 ，计算结果不稳定。1. FDM不采用虚拟节点时：u2u311/ 3节点 2： u32u2u12u2f 2(1/ 3)2节点 3： u42u3u22u3f 3(1/ 3)2代入 u41解得： u30.3516采用在节点 1 给出 du 的二阶精度公式3u14u2 u31d x2 x可解得： u329FVM节点 1： u2u11 f2xx2节电 2： u3u2u2u12u2 x f 2xxx节点 3： u4u3u3u22u3 x

9、f 3xxx代入 u41, x1/ 3, f238/ 9, f332 / 9解得： u32 / 9精确解 u2x2x.对比可得 FVM 得出的解等于精确解，即 Neumann边界条件在 FVM 中精确实现。而 FDM 的解有一定误差，因为在 FDM 中难以实现 Neumann边界条件，但可采用虚拟节点法，二阶精度法等来改进 FDM 以得到精确解。2. 前向差分ui 1uix(u( x)2(2u( x)33u( x)44u)i2!x2 ) i3!(x3 )i4!(x4 )ixui 2ui2 x( u) i(2 x) 2 (2u2 )i(2 x)3(3u3 ) i(2 x) 4 (4u4 )ix2

10、!x3!x4!xui 3ui3 x( u)i(3 x)2(2u2 )i(3 x)3(3u3 )i(3 x)4(4 u4 )ix2!x3!x4!x令(u)iauibui 1cui 2dui 3(x3)xx列方程解得： a11, b3,c3 , d1623后向差分ui 1ui(x)( u)i(x) 2 (2u2)i(x)3(3u3 )i(x) 4 (4u4 )ix2!x3!x4!xui 2ui(2x)(u )i(2x)2(2u2 ) i(2x)3(3u3 ) i(2x)4(4 u4 )ix2!x3!x4!xui 3ui( 3 x)(u)i( 3 x) 2(2u( 3 x)3(3u( 3 x)4(4

11、 ux2!x2 ) i3!x3 ) i4!x4 )i令 (u)iauibui 1cui 2dui 3(x3)xx11 ,b3 , d1解得： a3,c623中心差分ui 3ui( 3 x)( u ) i( 3 x)2(2 u2 )i( 3 x)3 (3u3 )i( 3 x)4(4u4 )i( 3 x)5(5u5 )i( 3 x) 6 (6u6 )ix2!x3!x4!x5!x6!xui 2ui( 2 x)( u)i( 2 x)2(2u2 )i( 2 x)3 (3u3 )i( 2 x)4(4u4 )i( 2 x)5(5u5 )i( 2 x)6(6u6 )ix2!x3!x4!x5!x6!xui1ui

12、(x)(u )i(x)2(2u2 )i(x)3(3u3 )i(x)4(4u4 )i(x) 5(5u5 )i(x)6(6u6 ) ix2!x3!x4!x5!x6!xui1ui(x)(u )i(x)2(2 u2 )i(x)3(3u3 ) i(x) 4(4u4 )i(x)5(5u5 )i( x)6(6 u6 )ix2!x3!x4!x5!x6!x.ui2ui(2x)(u )i(2x) 2(2u2 )i(2 x)3(3u3 )i(2 x) 4(4u4 )i(2 x)5(5u5 )i(2 x)6(6u6 )ix2!x3!x4!x5!x6!xui3ui(3x)(u )i(3x) 2(2u2 )i(3 x)3

13、(3u3 ) i(3 x)4(4u4 )i(3 x) 5(5u5 )i(3 x)6(6u6 )i令x2!x3!x4!x5!x6!x( u5u8 ) s5,8( u9u8 ) s9,8(u7u8 ) s7,8f 8 A8( u )iaui 3bui 2cui 1duieui 1fui 2gui 3( x6 )xxa0.0167, b0.15, c0.75解得：0, e0.75, f0.15, g 0.0167d4.uin 1uin 1(uin12uinuin1 )( t 2 , x2 ) （1）2tx2uini假定任一个时间步 n 的数值解uninn 1n 1nnn代入得 ii( i 12ii

14、1) （2）2tx2将 in 分解为傅立叶级数nNnIije.（ 3）ijN将（3）式代入（ 2）式得n 1n1(nnn)IieI (i 1)2 eIieI ( i 1)2tex2令 d2tx2n 1n1dn2 e I) .（ 4）(eIn令误差放大因子 fn(n 1) 代入（ 4）得fn 11cos) .（ 5）2d (1f nfn1cos) .（6）2d (1f n 1.这里采用反证法，假设其有条件稳定则由（ 5）（6）可推出 d0, fn 0 再由（5）（6）得fn 111而该方程与假设（ fn 1 f n , fnf n 1 ）矛盾fnfn 1f n固假设不成立，即（ 1）式无条件不稳定。5.节点 5 的差分方程为( u2u5) s2,5( u6u5 ) s6,5(u8u5 ) s8,5(u4 u5 )s4,5 f