关键词:虚拟单元边界法;剪切层;超声速流动 Abstract An improved ghost-cell immersed boundary method proposed in this paper, coupled with a high order finite difference solver, is applied to simulate the supersonic compressive flows around the complex obstacles.The main improvement of this algorithm is the treatment of the solid boundary that both ghost points inside the solid domain and forcing points inside the fluid domain due to the extension of the boundary are chosen to reconstruct the flow information considering the effect of solid wall on fluid.This brings refined boundary with discrete points and strengthens the wall conditions, which plays the role of local mesh refinement.The fluid points are limited in a certain source space as the interpolating points of the inverse distance algorithm, which effectively avoids the fact that the interpolating points are too few to possibly lead to coincide with the forcing points.Two problems of two dimensional shock reflection (Ma=2.81) and three dimensional flow around the smooth sphere (Ma=1.2) demonstrate the significant improvement of the numerical accuracy compared to the general ghost cell method. The results reveal the instability mechanism of the free shear layer as a result of the interaction between the shear layer, the compression wave system and the wake. The thickness and Reynolds fluctuation of the shear layer experience three regimes of linear growth, large amplitude oscillation and small amplitude fluctuation, resulting in an exponential growth of wrinkling factor.The turbulent structure near the shear layer shows obvious anisotropy because the streamwise Reynolds normal stress is dominant and a spatial hysteresis exists in the effect of the tail shock on Reynolds normal stresses in different directions.
KadochB, ReimannT, SchneiderK, et al.Comparison of a spectral method with volume penalization and a finite volume method with body fitted grids for turbulent flows . Computers & Fluids, 2016, 133: 140-150 [本文引用: 1]
[4]
SotiropoulosF, YangX. Immersed boundary methods for simulating fluid-structure interaction . Progress in Aerospace Sciences, 2014, 65: 1-21 [本文引用: 1]
[5]
PeskinCS. Flow patterns around heart valves: A numerical method . Journal of Computational Physics, 1972, 10(2): 252-271 [本文引用: 1]
[6]
GoldsteinD, HandlerR, SirovichL. Modeling a no-slip flow boundary with an external force field . Journal of Computational Physics, 1993, 105(2): 354-366 [本文引用: 2]
[7]
Saiki EM, BiringenS. Numerical simulation of a cylinder in uniform flow: Application of a virtual boundary method . Journal of Computational Physics, 1996, 123(2): 450-465 [本文引用: 1]
[8]
LeeC. Stability characteristics of the virtual boundary method in three-dimensional applications . Journal of Computational Physics, 2003, 184(2): 559-591 [本文引用: 1]
[9]
KolomenskiyD, SchneiderK. A Fourier spectral method for the Navier-Stokes equations with volume penalization for moving solid obstacles . Journal of Computational Physics, 2009, 228(16): 5687-5709 [本文引用: 1]
[10]
KadochB, KolomenskiyD, AngotP, et al.A volume penalization method for incompressible flows and scalar advection-diffusion with moving obstacles . Journal of Computational Physics, 2012, 231(12): 4365-4383 [本文引用: 1]
Mohd-YosufJ. Combined immersed boundary/B-spline methods for simulations of flow in complex geometries, Annual Research Briefs, Center for Turbulence Research, NASA Ames Research Center/Stanford University1997, 317-328 [本文引用: 2]
[13]
FadlunEA, VerziccoR, OrlandiP, et al.Combined immersed-boundary finite-difference methods for three-dimensional complex flow simulations . Journal of Computational Physics, 2000, 161(1): 35-60 [本文引用: 2]
[14]
MittalR, DongH, BozkurttasM, et al.A versatile sharp interface immersed boundary method for incompressible flows with complex boundaries . Journal of Computational Physics, 2008, 227(10): 4825-4852 [本文引用: 1]
[15]
MajumdarS, IaccarinoG, DurbinP, RANS solvers with adaptive structured boundary non-conforming grids, Annual Research Briefs, NASA Ames Research Center/Stanford University. Stanford, CA, 2001, 353-366 [本文引用: 2]
[16]
TsengYH, FerzigerJH. A ghost-cell immersed boundary method for flow in complex geometry . Journal of Computational Physics, 2003, 192(2): 593-623 [本文引用: 1]
[17]
QuY, ShiR, BatraRC. An immersed boundary formulation for simulating high-speed compressible viscous flows with moving solids . Journal of Computational Physics, 2018, 354: 672-691 [本文引用: 2]
[18]
GhiasR, MittalR, DongH. A sharp interface immersed boundary method for compressible viscous flows . Journal of Computational Physics, 2007, 225(1): 528-553 [本文引用: 1]
[19]
ChaudhuriA, HadjadjA, ChinnayyaA. On the use of immersed boundary methods for shock/obstacle interactions . Journal of Computational Physics, 2011, 230(5): 1731-1748 [本文引用: 3]
[20]
BouchonF, DuboisT, JamesN. A second-order cut-cell method for the numerical simulation of 2D flows past obstacles . Computers & Fluids, 2012, 65: 80-91 [本文引用: 1]
[21]
IkenoT, KajishimaT. Finite-difference immersed boundary method consistent with wall conditions for incompressible turbulent flow simulations . Journal of Computational Physics, 2007, 226(2): 1485-1508 [本文引用: 1]
[22]
SeoJH, MittalR. A high-order immersed boundary method for acoustic wave scattering and low-Mach number flow-induced sound in complex geometries . Journal of Computational Physics, 2011, 230(4): 1000-1019 [本文引用: 1]
[23]
HaugenNEL, KragsetS. Particle impaction on a cylinder in a crossflow as function of Stokes and Reynolds numbers . Journal of Fluid Mechanics, 2010, 661: 239-261 [本文引用: 1]
[24]
LuoK, ZhuangZ, FanJ, et al.A ghost-cell immersed boundary method for simulations of heat transfer in compressible flows under different boundary conditions . International Journal of Heat and Mass Transfer, 2016, 92: 708-717 [本文引用: 2]
[25]
LuoK, MaoC, ZhuangZ, et al.A ghost-cell immersed boundary method for the simulations of heat transfer in compressible flows under different boundary conditions Part-II: Complex geometries . International Journal of Heat and Mass Transfer, 2017, 104: 98-111 [本文引用: 1]
[26]
LiXL, FuDX, MaYW, et al.Direct numerical simulation of shock/turbulent boundary layer interaction in a supersonic compression ramp . Science China Physics, Mechanics & Astronomy, 2010, 53(9): 1651-1658 [本文引用: 1]
[27]
LiX, FuD, MaY. Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack . Physics of Fluids, 2010, 22(2): 025105 [本文引用: 1]
(TongFulin, LiXinliang, TangZhigong. Numerical analysis of unsteady motion in shock wave/transitional boundary layer interaction . Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(1): 93-104 (in Chinese)) [本文引用: 2]
[29]
StegerJL, WarmingRF. Flux vector splitting of the inviscid gasdynamic equations with application to finite-difference methods . Journal of Computational Physics, 1981, 40(2): 263-293 [本文引用: 1]
[30]
JiangGS, ShuCW. Efficient implementation of weighted ENO schemes . Journal of Computational Physics, 1996, 126(1): 202-228 [本文引用: 2]
[31]
FrankeR. Scattered data interpolation: Tests of some methods . Mathematics of Computation, 1982, 38(157): 181-200 [本文引用: 1]
[32]
YangJ, LiuY, LomaxH. Computation of shock wave reflection by circular cylinders . AIAA Journal, 1987, 25(5): 683-689 [本文引用: 1]
[33]
GiepmanRHM, SchrijerFFJ, Van OudheusdenBW. Flow control of an oblique shock wave reflection with micro-ramp vortex generators: Effects of location and size . Physics of Fluids, 2014, 26(6): 066101 [本文引用: 1]
[34]
IgraD, TakayamaK, IgraO. Shock wave reflection over roughened wedges// 30th International Symposium on Shock Waves, Cham: Springer, 2017: 621-625 [本文引用: 1]
[35]
FangJ, YaoY, ZheltovodovAA, et al.Investigation of three dimensional shock wave turbulent boundary layer interaction initiated by a single fin . AIAA Journal, 2016, 55(2): 509-523 [本文引用: 1]
[36]
BrysonAE, GrossRWF. Diffraction of strong shocks by cones, cylinders, and spheres . Journal of Fluid Mechanics, 1961, 10(1): 1-16 [本文引用: 4]
[37]
KacaJ. An interferometric investigation of the diffraction of a planar shock wave over a semicircular cylinder . NASA STI/Recon Technical Report N, 1988, 89 [本文引用: 1]
[38]
ConstantinescuG, SquiresK. Numerical investigations of flow over a sphere in the subcritical and supercritical regimes . Physics of Fluids, 2004, 16(5): 1449-1466 [本文引用: 1]
[39]
RodriguezI, BorellR, LehmkuhlO, et al.Direct numerical simulation of the flow over a sphere at Re = 3700. Journal of Fluid Mechanics, 2011, 679: 263-287 [本文引用: 1]
(HuZhengmin, XiaoZhijian, TaoZhu, et al.LES investigations of flow over a sphere in the subcritical and supercritical regimes . Chinese Journal of Hydrodynamics, 2016, 31(4): 496-502 (in Chinese)) [本文引用: 1]
[41]
Von TerziDA, SandbergRD, FaselHF. Identification of large coherent structures in supersonic axisymmetric wakes . Computers & Fluids, 2009, 38(8): 1638-1650 [本文引用: 1]
[42]
YunG, KimD, ChoiH. Vortical structures behind a sphere at subcritical Reynolds numbers . Physics of Fluids, 2006, 18(1): 015102 [本文引用: 1]
(LinMengda, CuiGuixiang, ZhangZhaoshun, et al.Large eddy simulation on the evolution and the fast-time prediction of aircraft wake vortices . Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(6): 1185-1200 (in Chinese)) [本文引用: 1]
(JiChunning, HuaYang, XuDong, et al.Numerical simulation of vortex-induced vibration of a flexible cylinder exposed to shear flow at different shear rates . Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(1): 21-31 (in Chinese)) [本文引用: 1]
(RenAnlu, LiGuangwang, ZouJianfeng. Numerical study of uniform flow over sphere at intermediate Reynolds numbers . Journal of Zhejiang University (Engineering Science ), 2004, 38(5): 644-648 (in Chinese)) [本文引用: 1]
(ZouJianfeng, RenAnlu, DengJian. Numerical investigations of wake and force for flow past a sphere . Acta Aerodynamica Sinica, 2004, 22(3): 303-308 (in Chinese)) [本文引用: 1]
[47]
AchenbachE. Experiments on the flow past spheres at very high Reynolds numbers . Journal of Fluid Mechanics, 1972, 54(3): 565-575 [本文引用: 1]