湍流反应流模拟.pptxVIP

第一章 湍流、湍流反应流、数值模拟;1.1.2 湍流流动的性质 ;图1.1 Titan IV火箭发动机形成的湍流羽流图片。 喷嘴出口直径为3m，估计的羽流高度为1500m，估计Re数为 2×108。;;图1.2 湍流射流浓度场的平面影像 a) Re=5000, b) Re=20000。 Dahm和Dimotakis的实验（1990）;图1.3 表示一个湍流射流，U1(t)表示在中心线上，轴向速度随时间变化。水平线表示平均速度，用U1 表示。 Tong和Warhaft的实验（1995）;;图1.4 表示在垂直来流的方向上，平均速度U1的变化曲线。U10 表示中心线上平均速度；x2/x1 表示离中心的相对距离。可以发现：平均速度U1有光滑的轮廓，没有细微的结构。即：湍流平均速度U1的曲线轮廓与层流射流没有什么不同。 Hussein、Capp 和 George 的结果（1994）;;[3]湍流的输运、混合性质 ;;[5] Reynolds管道流动实验中的层流和湍流;[2]染料色线前进很小的距离时，就完全与周围的水混合，减低了染料的浓度。;[6]湍流加强输运和混合的重要性; 同样，湍流流动加强了热量传递。如：固－液、液－气界面上的热和质量传递速率在湍流中被大大加强。 ;1.2 湍流反应流;1.3 数值模拟;1.3.2 数值解法特点;1.3.3 数值解基本步骤;1.3.4 数值解法种类（包括离散方法、模型、技术）;[5]直接数值模拟DNS；（技术？）;注意： 计算结果与采用的： （1）空间离散格式； （2）时间积分方法； （3）边界条件处理； （4）模型，如燃烧模型，湍流模型，等 （5）求解策略，如求解方程中的次序； ……等等，有关。可以说：不同程序有不同的“个性”。;问题： 1）FDS采用的什么数值离散方法？ 2）FDS的时间离散精度是几阶？ 3）FDS的空间离散精度是几阶？ 4）FDS采用了什么模拟技术？燃烧模型 ;1.3.5 有限分析法举例;;讨论一 控制方程组的求解策略？ 需要求解的量为： [1]密度 RHO(I,J,K)； [2]速度 U(I,J,K)，V(I,J,K)，W(I,J,K)； [3]压力 P(I,J,K)； [4]温度TEMP(I,J,K)； [5]组分浓度Yk(I,J,K) 方程组为：（1）连续性方程；（2）动量方程；（3）能量方程；（4）状态方程；（5）组分方程.;求解策略一;非均匀网格：[1]采用TANH函数，所谓拉伸网格； [2]FDS中采用了多项式和分段线性拉伸网格; CALL INIT ！计算量初始化（TIME=0时），如密度、速度、组分浓度、温度（内能） C +++ NCYC=0 ! 循环数 TIME=0 ！初始时间赋值 C +++ C +++ TIME ADVANCE C +++ MAIN_LOOP: DO; CALL TIMESTEP ! 计算时间步长DT，涉及 ！ CFL数等内容 TIME=TIME+DT NCYC=NCYC+1 IF(NCYC .EQ. LNCY) THEN STOP_FLAG=.TRUE. ENDIF C +++ IF(TIME .GE. T_MAX) THEN STOP_FLAG=.TRUE. ENDIF IF (STOP_FLAG) EXIT MAIN_LOOP;DO NRK=1，Nk ! Runge-Kutta ! Adams-Bashforth 时间显式向 ！前积分 CALL RHS ！计算控制方程的右手边 GET-RHO ！获得密度 GET-U ! 获得速度 GET-TEMP（or Internal Energy）！获得温度等 C +++ CALL BC ！给定边界条件 END DO ！时间积分结束 ENDDO MAIN_LOOP ! 主循环结束; CALL OUTP