Numerical simulation for 2D double-diffusive convection(DDC) in rectangular enclosures based on a high resolution upwind compact streamfunction model Ⅰ: numerical method and code validation  

作　　者：Yaping YAN Shuang WU Kaiyuan TIAN Zhenfu TIAN 

机构地区：[1]Institute of Mechanics and Computational Engineering,Department of Aeronautics and Astronautics,Fudan University,Shanghai 200433,China [2]College of Science,Liaoning Technical University,Fuxin 123000,Liaoning Province,China [3]Research Institute of Intelligent Complex Systems,Fudan University,Shanghai 200433,China [4]School of Life Sciences,Fudan University,Shanghai 200433,China

出　　处：《Applied Mathematics and Mechanics(English Edition)》2022年第9期1431-1448,共18页应用数学和力学（英文版）

基　　金：supported by the National Natural Science Foundation of China(Nos.11872151,11372075,and 91330112)。

摘　　要：A high resolution upwind compact streamfunction numerical algorithm for two-dimensional(2D)double-diffusive convection(DDC)is developed.The unsteady Navier-Stokes(N-S)equations in the streamfunction-velocity form and the scalar temperature and concentration equations are used.An optimized third-order upwind compact(UCD3 opt)scheme with a low dispersion error for the first derivatives is utilized to approximate the third derivatives of the streamfunction in the advection terms of the N-S equations and the first derivatives in the advection terms of the scalar temperature and concentration equations.The remaining first derivatives of the streamfunction(velocity),temperature,and concentration variables used in the governing equations are discretized by the fourth-order compact Pade(SCD4)schemes.With the temperature and concentration variables and their approximate values of the first derivatives obtained by the SCD4 schemes,the explicit fourth-order compact schemes are suggested to approximate the second derivatives of temperature and concentration in the diffusion terms of the energy and concentration equations.The discretization of the temporal term is executed with the second-order Crank-Nicolson(C-N)scheme.To assess the spatial behavior capability of the established numerical algorithm and verify the developed computer code,the DDC flow is numerically solved.The obtained results agree well with the benchmark solutions and some accurate results available in the literature,verifying the accuracy,effectiveness,and robustness of the provided algorithm.Finally,a preliminary application of the proposed method to the DDC is carried out.

关 键 词：double diffusive convection(DDC) high resolution heat and mass transfer upwind compact difference streamfunction formulation 

分 类 号：O241.8[理学—计算数学] O357.1[理学—数学]