參考文獻
[1] Sandham, N. D and Reynolds, W. C. , “Three dimensional simulations of large eddies in the compressible mixing layer,’’ J. Fluid Mech, Vol.224, 1991, pp.
133-158.
[2] Magi, V. and Abraham, J. , “Exploring velocity and density ratio effects in a mixing layer using DNS, ” IJCFD, Vol. 8,1997, pp. 147-151.
[3] Brown, G.L. and Roshko, A., “On density effect and large structure in turbulent mixing layers,” J. Fluid Mech, Vol.64, pp.775-816, (1974).
[4] Sandham, N. D and Reynolds, W. C. , “Growth of oblique waves in the mixing layers at high Mach number,’’ Turbulent Shear Flows 7, 1989, pp. 1-6.
[5] Grinstein, F.F. and Guirguis, R.H., “Effective viscosity in the simulation of
spatially evolving shear flows with monotonic FCT models,” J.Comput.Phys, Vol.101, p.p.165-175, (1992)
[6] Choi, H, and Lasheras, J. C. , “Three dimensional instability of a plane free shear layer : An experimental study of the formation and evolution of streamwise votices,” J. Fluid Mech, Vol. 189, 1988, pp. 53-88.
[7] Bogdanoff, D. W. , “Compressibility effects in turbulent shear layer, ” AIAA J, Vol.
21, 1983, pp. 926-927.
[8] Roshko, A. and Papamoschou, D. , “The compressible turbulent shear：An experimental study,” J. Fluid Mech, Vol. 197, 1988, pp. 453-477.
[9] Samimy, M. , Reeder, M. F. and Elliott, G. S. , “Compresibility effects on large structures in the free shear flows,” Phys. Fluids A, Vol. 4,1992, pp. 1251-1258.
[10] Grinstein, F. F. , “Open boundary conditions in the simulation of subsonic turbulent shear flows, ” J. Comput. Phys, Vol. 114, 1994, pp. 43-55.
[11] Tsai, Y. P. and Christiansen, W. H. , “Two dimensional numerical simulation of shear layer optics,” AIAA J, Vol. 28, 1990, pp. 2092-2097.
[12] Soh, W. Y. , “Numerical simulation of free shear flows and far field sound pressure directivity,” International Journal For Numeriacl Methods In Fluids, Vol.
18, 1994, pp. 337-359.
[13] Ansari, A.,”Self similarity and mixing characteristics of turbulent mixing layers starting from laminar initial conditions,” Phys. Fluid, Vol.9, p.p.1714-1728, (1997)
[14] Liu, X.D., Osher, S. and Chan, T., “Weighted Essentially Non-Oscillatory Schemes,” J.Comput. Phys, Vol.115, p.p.200-212, (1994)
[15] Jiang, G.S. and Shu, C.W., “Efficient Implementation of Weighted ENO Scheme,” J.Comput.Phys, Vol.126, p.p.202-228, (1996)
[16] Grinstein, F.F., Oran, E.S. and Boris, J.P., “Numerical simulations of asymmetric mixing in planar shear flows,” J.Fluid Mech, Vol.165, p.p.201-220, (1986)
[17] 韋立興，”The Physics of Supersonic Mixing Layer：A Numerical Inverstigation”, 中華大學機械與航太工程研究所碩士班論文,2002 年 7 月。
[18] Messersmith, N.L., Goebel, S.G., Frantz, W.H., Krammer, E.A., Renie, J.P., Dutton, J.C. and Krier, H., “Experimental and analytical investigations of supersonic mixing layers,” AIAA paper, (1988)
[19] Papamoschou, D. and Rochko, A., “Observations of Supersonic Free Shear Layers,”
AIAA paper 86-0162, (1986)
[20] Wu, K. C. and Lu, P. J, “Numerical investigation on the structure of a confined supersonic mixing layer,” Phys. Fluids A, Vol.3,1991, pp.3063-3079.
圖4-1 次音速流體剪切混合層之瞬間等密度線圖 密度比值λρ=0.64,速度比值λu=0.55(U1=0.5M)
圖4-2 次音速流體剪切混合層之瞬間等渦漩線圖 密度比值λρ=0.64,速度比值λu=0.55(U1=0.5M)
圖4-3 次音速流體剪切混合層之瞬間等密度線圖 密度比值λρ=0.64,速度比值λu=0.55(U1=0.8M)
圖4-4 次音速流體剪切混合層之瞬間等渦漩線圖 密度比值λρ=0.64,速度比值λu=0.55(U1=0.8M)
圖4-5 超音速流體剪切混合層之瞬間等渦漩線圖 (密度比值λρ=1.55,速度比值λu=0.43,時間Time=0.005s, M1=1.4)
圖4-6 超音速流體剪切混合層之瞬間等渦漩線圖 (密度比值λρ=0.5,速度比值λu=0.35,時間Time=0.002s, M1=1.6)
圖4-7 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=0.64,速度比值λu=0.35,時間Time=0.0019s, M1=1.8)
圖4-8 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=1.55,速度比值λu=0.57,時間Time=0.004s, M1=1.9)
圖4-9 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=0.6,速度比值λu=0.16,時間Time=0.001s, M1=2.35)
圖4-10 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=0.37,速度比值λu=0.24,時間Time=0.0008s, M1=3.01)
圖4-11 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=0.64,速度比值λu=0.35,時間Time=0.001s,M1=1.80)
圖4-12 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=0.64,速度比值λu=0.35,時間Time=0.0013s,M1=1.80)
圖4-13 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=0.64,速度比值λu=0.35,時間Time=0.0019s,M1=1.80)
圖4-14 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=0.64,速度比值λu=0.35,時間Time=0.0022s,M1=1.80)
圖4-15 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=1.55,速度比值λu=0.57,時間Time=0.002s,M1=1.91)
圖4-16 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=1.55,速度比值λu=0.57,時間Time=0.003s,M1=1.91)
圖4-17 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=1.55,速度比值λu=0.57,時間Time=0.0035s,M1=1.91)
圖4-18 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=1.55,速度比值λu=0.57,時間Time=0.004s,M1=1.91)
圖4-19 超音速流體剪切混合層之瞬間等密度線圖 (密度比值λρ=1.55,速度比值λu=0.43,時間Time=0.005s, M1=1.4)
圖4-20 超音速流體剪切混合層之瞬間等密度線圖 (密度比值λρ=0.5,速度比值λu=0.35,時間Time=0.002s, M1=1.6)
圖4-21 超音速流體剪切混合層之瞬間等密度線圖
(密度比值λρ=0.64,速度比值λu=0.35,時間Time=0.0019s, M1=1.8)
圖4-22 超音速流體剪切混合層之瞬間等密度線圖
(密度比值λρ=1.55,速度比值λu=0.57,時間Time=0.004s, M1=1.91)
圖4-23 超音速流體剪切混合層之瞬間等密度線圖 (密度比值λρ=0.6,速度比值λu=0.16,時間Time=0.001s, M1=2.35)
圖4-24 超音速流體剪切混合層之瞬間等密度線圖
(密度比值λρ=0.37,速度比值λu=0.24,時間Time=0.0008s, M1=3.01)
圖4-25 次音速流體剪切混合層之瞬間等密度線圖
(密度比值λρ=1.55,速度比值λu=0.80,時間Time=0.004s, M1=0.7)
圖4-26 超音速流體剪切混合層之瞬間等密度線圖
(密度比值λρ=1.55,速度比值λu=0.43,時間Time=0.004s, M1=1.4)
圖4-27 超音速流體剪切混合層之瞬間等密度線圖
(密度比值λρ=1.55,速度比值λu=0.57,時間Time=0.004s, M1=1.91)
圖4-28 次音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=1.55,速度比值λu=0.80,時間Time=0.004s, M1=0.7)
圖4-29 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=1.55,速度比值λu=0.43,時間Time=0.004s, M1=1.4)
圖4-30 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λ =1.55,速度比值λ =0.57,時間Time=0.004s, M =1.91)
圖4-31 次音速流體剪切混合層之瞬間等渦漩線圖 (密度比值λρ=1,速度比值λu=0.35,時間Time=0.002s, M1=0.6)
圖4-32 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=0.5,速度比值λu=0.35,時間Time=0.002s, M1=1.6)
圖4-33 超音速流體剪切混合層之瞬間等渦漩線圖
(密度比值λρ=0.64,速度比值λu=0.35,時間Time=0.002s, M1=1.8)
圖4-34 次音速流體剪切混合層之瞬間等密度線圖 (密度比值λρ=1,速度比值λu=0.35,時間Time=0.002s, M1=0.6)
圖4-35 超音速流體剪切混合層之瞬間等密度線圖
(密度比值λρ=0.5,速度比值λu=0.35,時間Time=0.002s, M1=1.6)
圖4-36 超音速流體剪切混合層之瞬間等密度線圖
(密度比值λρ=0.64,速度比值λu=0.35,時間Time=0.002s, M1=1.8)
Sarlar and Lakshmanan
0.0 0.2 0.4 0.6 0.8 1.0
Langley curve Elliott and Samimy Hall,Dimotalis &Rosemann U2/U1=0.55, =0.64 U2/U1=0.54, =0.64 U2/U1=0.57, =1.55 U2/U1=0.35, =0.64 U2/U1=0.24, =0.37
0.0 0.5 1.0 1.5 2.0
M c
圖4-37 無因次的剪切混合層成長率與對流馬赫數之關係圖
圖4-38 超音速流體對流馬赫數MC=0.51 之平均流線方向速度輪廓之自我相似說明圖 (密度比值λρ=0.64,速度比值λu=0.35, M1=1.8)
圖4-39 超音速流體對流馬赫數MC=0.45 之平均流線方向速度輪廓之自我相似說明圖 (密度比值λρ=1.55,速度比值λu=0.57, M1=1.91)
圖4-40 超音速流體對流馬赫數MC=0.86 之平均流線方向速度輪廓之自我相似說明圖 (密度比值λρ=0.37,速度比值λu=0.24, M1=2.35)
圖4-41 超音速流體對流馬赫數MC=0.86 之平均流線方向速度輪廓之自我相似說明圖 (密度比值λρ=0.6,速度比值λu=0.16, M1=3.01)
(A) 密度
(B) 壓力
圖4-42 超音速情況下瞬時之密度與壓力圖
(M1=1.91 馬赫，M2=1.36 馬赫，MC=0.45 馬赫，時間Time=0.004s)
(A) 流線方向速度擾動輪廓圖
(B) 跨流線方向速度擾動輪廓圖
(C) 雷諾主剪應力擾動輪廓圖
圖4-43 超音速流體(密度比值λρ=1.55,速度比值λu=0.43,時間Time=0.005s, M1=1.4)
(A) 流線方向速度擾動輪廓圖
(B) 跨流線方向速度擾動輪廓圖
(C) 雷諾主剪應力擾動輪廓圖
圖4-44 超音速流體(密度比值λρ=0.5,速度比值λu=0.53,時間Time=0.002s, M1=1.6)
(A) 流線方向速度擾動輪廓圖
(B) 跨流線方向速度擾動輪廓圖
(C) 雷諾主剪應力擾動輪廓圖
圖4-45 超音速流體(密度比值λρ=0.64,速度比值λu=35,時間Time=0.0019s, M1=1.81)
(A) 流線方向速度擾動輪廓圖
(B) 跨流線方向速度擾動輪廓圖
(C) 雷諾主剪應力擾動輪廓圖
圖4-46 超音速流體(密度比值λρ=1.55,速度比值λu=0.57,時間Time=0.004s, M1=1.91)
表4- 1 混合剪切流之計算條件
M1 U2/U1
2/ 1 McCase 1 0.5 0.64 0.55 0.1
Case 2 0.6 1.0 0.35 0.2
Case 3 0.7 1.55 0.80 0.08
Case 4 0.8 0.64 0.55 0.16
Case 5 1.4 1.55 0.43 0.44
Case 6 1.6 0.5 0.35 0.42
Case 7 1.8 0.64 0.35 0.51
Case 8 1.91 1.55 0.57 0.45
Case 9 2.35 0.6 0.16 0.86
Case 10 3.01 0.37 0.24 0.86
附錄一
網格測試
343*43
523*83
743*103