New method for three-dimensional convection on large grids

In the present work the method of solution of three-dimensional hydrodynamic equations in Bussinesq approach is presented. The finite-difference method is used for spatial derivative and implicit method is used for time derivative. Global nonlinear system of equations is linearized and solved by the stabilization method at each time step. To solve the linearized equations of convectional heat-mass transfer we apply the method using exponential transformation and the conjugate gradient method with preconditioning by the incomplete factorisation method.