Численное моделирование конвективного тепломассопереноса в сферических координатах

Abstract

Целью исследования является построение дискретного аналога обобщенного диф-ференциального уравнения, описывающего конвекцию в вязкой несжимаемой жидкости в сферических координатах. Математическая модель конвективного тепломассопереноса в вязкой несжимаемой жидкости задается системой дифференциальных уравнений, полученных на основе уравнений гидродинамики, тепло- и массообмена. Эти уравнения подчиняются обобщенному закону сохранения, который описывается дифференциальным уравнением для обобщенной переменной. Для дискретизации дифференциального уравнения используется метод контрольного объема. Расчетная область разбивается на множество непересекающихся контрольных объемов с узловой точкой в каждом из них. Дифференциальное уравнение интегрируется по контрольным объемам. В результате получается дискретный аналог, связывающий значение обобщенной переменной в узловой точке с ее значениями в соседних узлах. Метод гарантирует строгое выполнение законов сохранения как во всей расчетной области, так и в любой ее части. Чтобы применять лучшие аппроксимации профилей обобщенной переменной, находятся точные решения уравнений сохранения отдельно по каждой координате. Кратко поясняется физический смысл точных решений. В итоге строится дискретный аналог для обобщенного дифференциального уравнения с использованием полученных аналитических решений. The aim of the research is to construct a discrete analogue of the generalized differential equation describing convection in a viscous incompressible fluid in spherical coordinates. The mathematical model of convective heat and mass transfer in a viscous incompressible fluid is given by a system of differential equations derived from the equations of hydrodynamics, heat and mass transfer. These equations satisfy the generalized conservation law, which is described by a differential equation for the generalized variable. The control volume method is used to obtain a discrete analogue of the differential equation. The computational domain is divided into a multiplicity of control volumes with a node in each of them. As a result, a discrete analogue is obtained that relates the value of the generalized variable at the node point to its values at neighboring nodes. The method guarantees strict compliance of conservation laws both in the entire calculation area and in any part of it. To apply the best approximation of the profiles of the generalized variable, there are exact solutions of the conservation equation separately for each coordinate. The physical meaning of exact solutions is briefly explained. As a result, a discrete analogue is constructed for the generalized differential equation using the obtained analytical solutions.