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

doi:10.1007/s10483-022-2895-6

Abstract

Abstract: 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 Padé (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.

Key words: double diffusive convection (DDC), high resolution, heat and mass transfer, upwind compact difference, streamfunction formulation

2010 MSC Number:

Yaping YAN, Shuang WU, Kaiyuan TIAN, Zhenfu TIAN. Numerical simulation for 2D double-diffusive convection (DDC) in rectangular enclosures based on a high resolution upwind compact streamfunction model I: numerical method and code validation. Applied Mathematics and Mechanics (English Edition), 2022, 43(9): 1431-1448.

TrendMD

References

[1] TABRIZINEJADAS, S., FAHS, M., ATAIE-ASHTIANI, B., SIMMONS, C. T., ROUPERT, R. D. C., and YOUNES, A. A Fourier series solution for transient three-dimensional thermohaline convection in porous enclosures. Water Resources Research, 56, e2020WR028111(2020)
[2] TURNER, J. S. Double-diffusive phenomena. Annual Review of Fluid Mechanics, 6, 37-56(1974)
[3] HUPPTER, H. E. and TURNER, J. S. Double-diffusive convection. Journal of Fluid Mechanics, 106, 299-329(1981)
[4] SCHMITT, R. W. Double diffusion in oceanography. Annual Review of Fluid Mechanics, 26, 255-285(1994)
[5] CHAMKHA, A. J. Double-diffusive convection in a porous enclosure with cooperating temperature and concentration gradients and heat generation or absorption effects. Numerical Heat Transfer; Part A:Applications, 41, 65-87(2002)
[6] ABIDI, A., KOLSI, L., BORJINI, M. N., and BEN-AISSIA, H. Effect of heat and mass transfer through diffusive walls on three-dimensional double-diffusive natural convection. Numerical Heat Transfer; Part A:Applications, 53, 1357-1376(2008)
[7] CHEN, C. F. and CHAN, C. L. Stability of buoyancy and surface tension driven convection in a horizontal double-diffusive fluid layer. International Journal of Heat and Mass Transfer, 53, 1563-1569(2010)
[8] YU, P. X., XIAO, Z. C., WU, S., TIAN, Z. F., and CHENG, X. Z. High accuracy numerical investigation of double-diffusive convection in a rectangular cavity under a uniform horizontal magnetic field and heat source. International Journal of Heat and Mass Transfer, 110, 613-628(2017)
[9] RADKO, T. Double-Diffusive Convection, Cambridge University Press, Cambridge (2013)
[10] NISHIMURA, T., WAKAMATSU, M., and MOREGA, A. M. Oscillatory double-diffusive convection in a rectangular enclosure with combined horizontal temperature and concentration gradients. International Journal of Heat and Mass Transfer, 41, 1601-1611(1998)
[11] REDDY, N. and MURUGESAN, K. Magnetic field influence on double-diffusive natural convection in a square cavity-a numerical study. Numerical Heat Transfer; Part A:Applications, 71, 448-475(2017)
[12] SENTHIL KUMAR, D., MURUGESAN, K., and THOMAS, H. R. Numerical simulation of double diffusive mixed convection in a lid-driven square cavity using velocity-vorticity formulation. Numerical Heat Transfer; Part A:Applications, 54, 837-865(2008)
[13] WANG, J., YANG, M., and ZHANG, Y. Onset of double-diffusive convection in horizontal cavity with Soret and Dufour effects. International Journal of Heat and Mass Transfer, 78, 1003-1031(2014)
[14] DENG, Q. H., ZHOU, J., MEI, C., and SHEN, Y. M. Fluid, heat and contaminant transport structures of laminar double-diffusive mixed convection in a two-dimensional ventilated enclosure. International Journal of Heat and Mass Transfer, 47, 5257-5269(2004)
[15] GHORAYEB, K., KHALLOUF, H., and MOJTABI, A. Onset of oscillatory flows in doublediffusive convection. International Journal of Heat and Mass Transfer, 42, 629-643(1999)
[16] QIN, Q., XIA, Z. A., and TIAN, Z. F. High accuracy numerical investigation of double diffusive convection in a rectangular enclosure with horizontal temperature and concentration gradients. International Journal of Heat and Mass Transfer, 71,405-423(2014)
[17] LIANG, X., LIANG, X. L., FU, D. X., and MA, Y. W. Complex transition of double-diffusive convection in a rectangular enclosure with height-to-length ratio equal to 4:part I. Communications in Computational Physics, 6, 247-268(2009)
[18] LIANG, X., PENG, B., and TIAN, Z. F. Complex transition of double-diffusive convection in a rectangular enclosure with height-to-length ratio equal to 4:part II. International Journal of Heat and Mass Transfer, 135, 247-261(2019)
[19] MAHAPATRA, T. R., PAL, D., and MONDAL, S. Effects of buoyancy ratio on double-diffusive natural convection in a lid-driven cavity. International Journal of Heat and Mass Transfer, 57, 771-785(2013)
[20] CORCIONE, M., GRIGNAFFINI, S., and QUINTINO, A. Correlations for the double-diffusive natural convection in square enclosures induced by opposite temperature and concentration gradients. International Journal of Heat and Mass Transfer, 81, 811-819(2015)
[21] KRAMER, J., JECL, R., and SKERGET, L. Boundary domain integral method for the study of double diffusive natural convection in porous media. Engineering Analysis with Boundary Elements, 31, 897-905(2007)
[22] VERHAEGHE, F., BLANPAIN, B., and WOLLANTS, P. Lattice Boltzmann method for double diffusive natural convection. Physical Review E, 75, 046705(2007)
[23] WANG, L., SHI, B. C., CHAI, Z. H., and YANG, X. G. Regularized lattice Boltzmann model for double-diffusive convection in vertical enclosures with heating and salting from below. Applied Thermal Engineering, 103, 365-376(2016)
[24] LIU, Q., FENG, X. B., XU, X. T., and HE, Y. L. Multiple-relaxation-time lattice Boltzmann model for double-diffusive convection with Dufour and Soret effects. International Journal of Heat and Mass Transfer, 139, 713-719(2019)
[25] SHAO, Q., FAHS, M., YOUNES, A., and MAKRADI, A. A high-accurate solution for DarcyBrinkman double-diffusive convection in saturated porous media. Numerical Heat Transfer; Part B:Fundamentals, 69, 26-47(2016)
[26] BEN-ARTZI, M., CROISILLE, J. P., FISHELOV, D., and TRACHTENBERG, S. A pure-compact scheme for the streamfunction formulation of Navier-Stokes equations. Journal of Computational Physics, 205, 640-664(2005)
[27] YOUNES, A., FAHS, M., ZIDANE, A., HUGGENBERGER, P., and ZECHNER, E. A new benchmark with high accurate solution for hot-cold fluids mixing. Heat and Mass Transfer, 51, 1321-1336(2015)
[28] TIAN, Z. F. and YU, P. X. An effcient compact difference scheme for solving the streamfunction formulation of the incompressible Navier-Stokes equations. Journal of Computational Physics, 230, 6404-6419(2011)
[29] PANDIT, S. K. and KARMAKAR, H. An effcient implicit compact streamfunction velocity formulation of two dimensional flows. Journal of Scientific Computing, 68, 653-688(2016)
[30] YU, P. X. and TIAN, Z. F. Compact computations based on a stream-function-velocity formulation of two dimensional steady laminar natural convection in a square cavity. Physical Review E, 85, 036703(2012)
[31] YU, P. X., QIU, J. X., and TIAN, Z. F. Numerical investigation of natural convection in a rectangular cavity under different directions of uniform magnetic field. International Journal of Heat and Mass Transfer, 67, 1131-1144(2013)
[32] YU, P. X. and TIAN, Z. F. An upwind compact difference scheme for solving the streamfunctionvelocity formulation of the unsteady incompressible Navier-Stokes equations. Computers and Mathematics with Applications, 75, 3224-3243(2018)
[33] SHAH, A., GUO, H., and YUAN, L. A third-order upwind compact scheme on curvilinear meshes for the incompressible Navier-Stokes equations. Communication in Computational Physcis, 5, 712-729(2009)
[34] TIAN, Z. F., LIANG, X., and YU, P. X. A higher order compact finite difference algorithm for solving the incompressible Navier-Stokes equations. International Journal for Numerical Methods in Engineering, 88, 511-532(2011)
[35] CARPENTER, M. H., GOTTLIEB, D., and ABARBANEL, S. The stability of numerical boundary treatments for compact high-order schemes finite difference schemes. Journal of Computational Physics, 108, 272-295(1993)
[36] MOREGA, A. and NISHIMURA, T. Double-diffusive convection by a Chebyshev collocation method. Technology Reports of the Yamaguchi University, 5, 259-276(1996)