Three-dimensional numerical modeling of natural convection in underground cavities connected to the surface by a well
摘要
We present two- and three-dimensional numerical simulations of the natural convection flow within an artificial underground cavity. We consider a realistic geometry reconstructed from measurements (3D point cloud obtained with a photogrammetric method) of the Barcq cavity existing in Normandy region, France. The numerical model solves the Navier-Stokes equations coupled with the unsteady heat equation using a finite-element method for the space discretization. Time integration is based on a semi-implicit method and the Newton method to solve the resulting non-linear equation. To implement the method, we use the free software FreeFem++ offering efficient mesh refinement/adaptivity tools. The parallelization of the simulations relies upon the PETSc software library, to which FreeFem++ offers an easy-to-use interface. We analyze the onset of the natural convection and characterize the convective cells developing within the Barcq cavity for four Rayleigh (Ra) numbers (from