^{1}

^{1}

^{*}

^{2}

^{1}

^{3}

On the basis of a three-dimensional non stationary model of a convective cloud with detailed description of dynamic, thermodynamic and microphysical processes, numerical experiments were conducted to study the formation of parameters of convective clouds under unstable stratification of the atmosphere. Numerical experiments have been carried out to study the formation of convective processes in the atmosphere. The thermo hydrodynamic parameters in the zone of a thunderstorm cloud are determined, and regions with a vortex motion of air are identified. The main flows feeding the convective cloud in the mature stage are determined. Due to the means of visualization, the areas of formation and growth of precipitation particles are identified. In a three-dimensional form, the interaction of dynamic and thermodynamic processes is analyzed. The interaction of fields is manifested in the form of deformation of fields of thermodynamic parameters under the influence of dynamic processes. Trajectories of air streams around a cloud and the trajectories of drops in a cloud are determined. The results of numerical experiments confirm that dynamic processes significantly influence the formation of fields of thermodynamic parameters in the cloud, which also determine the course of microphysical processes and the nature of the growth of precipitation particles.

The development of physics of convective clouds and cloud seeding at the present stage requires the solution of a number of large and interrelated tasks of fundamental and applied character, which are complex research. Cloud formation and development processes are main of them, as include cloud clusters formation and interaction of clouds with surrounding atmosphere [ 1 - 3 ]. This is due to the fact that a convective cloud is a highly complex thermodynamic and microphysics system, important features of which are the nonstationarity, three-dimensionality and nonlinearity [ 3 - 7 ]. Therefore, the solution of problems of physics of convective clouds requires a holistic approach and the use of more effective methods, the most important among which is numerical modeling [ 1 - 10 ].

The aim of this work is the study on the basis of three-dimensional numerical model of convective clouds the formation of dynamic, thermodynamic and microstructural characteristics, and interaction of physical processes in clouds.

The subject of study is a three-dimensional convective cloud parameters: potential temperature, specific humidity, air flow, turbulence and distribution of drops and ice particles.

In the work briefly described is a three-dimensional non stationary model of convective cloud with detailed account for the thermodynamic, microphysical and electric processes.

The model includes the equations of hydrothermodynamics, microphysics and electrical processes that were used in the works [ 1 - 3 ].

Thermodynamic block of model consists of the equations of motion that describe moist convection in the Boussinesq approximation, which takes into account advective and turbulent transport, buoyancy force, friction and pressure gradients:

∂ u ∂ t + ( V ⋅ ∇ ) u = − ∂ π ′ ∂ x + Δ ′ u + l v ,

∂ v ∂ t + ( V ⋅ ∇ ) v = − ∂ π ′ ∂ y + Δ ′ v − l u , (1)

∂ w ∂ t + ( V ⋅ ∇ ) w = − ∂ π ′ ∂ z + Δ ′ w + g ( θ ′ / θ 0 + 0.61 s ′ − Q S ) ,

the equation of continuity

∂ u ∂ x + ∂ ν ∂ y + ∂ w ∂ z = σ w , (2)

the equations of thermodynamics

∂ θ ∂ t + ( V ⋅ ∇ ) θ = L k C p θ T δ M k δ t + L C C p θ T δ M C δ t + L З C p θ T δ M З δ t + Δ ′ θ , (3)

∂ s ∂ t + ( V ⋅ ∇ ) s = − δ M k δ t − δ M C δ t + Δ ′ s ,

where ( V ⋅ ∇ ) ≡ u ∂ ∂ x + v ∂ ∂ y + w ∂ ∂ z , Δ ′ = ∂ ∂ x K ∂ ∂ x + ∂ ∂ y K ∂ ∂ y + ∂ ∂ z K ∂ ∂ z ,

V = { u , v , w } ―velocity vector, u ( r , t ) , v ( r , t ) , w ( r , t ) ―components of the velocity vector of air flow in the cloud; θ ( r , t ) ―potential temperature; π ( r , t ) = c p θ ¯ ( P ( r , t ) / 1000 ) R / C p the dimensionless pressure; θ ¯ the average potential temperature; R―gas constant; s ( r , t ) the specific humidity of air, Q S ( r , t ) the total mixing ratio of liquid and solid phases in the cloud; σ ( z ) ―a parameter that takes into account the change of air density with height; P ( r , t ) and T ( r , t ) the pressure and the temperature, respectively; C_{p} the heat capacity of air at constant pressure; L C , L S , L F the specific heat of condensation, sublimation and freezing, respectively; π ′ ( r , t ) , θ ′ ( r , t ) , s ′ ( r , t ) ―deflection dimensionless pressure, potential temperature and specific humidity from their background values in ambient atmosphere

π 0 ( r ) , θ 0 ( r ) , s 0 ( r ) ; δ М C δ t , δ М S δ t the changes in specific humidity due to diffusion of vapor on the drops and crystals; δ М F δ t the mass of frozen water droplets in unit time in unit volume of air; K ( r , t ) the turbulent diffusion coefficient; l the Coriolis parameter. The vector r corresponds to the coordinates ( x , y , z ) . For the boundaries of the spatial domain symbols 0, L_{x}, 0, L_{y} 0, L_{z} are used.

Equations for the distribution functions of droplets and crystals by mass:

∂ f 1 ∂ t + u ∂ f 1 ∂ x + v ∂ f 1 ∂ y + ( w − V 1 ) ∂ f 1 ∂ z = ( ∂ f 1 ∂ t ) C + ( ∂ f 1 ∂ t ) C G + ( ∂ f 1 ∂ t ) A C + ( ∂ f 1 ∂ t ) D R + ( ∂ f 1 ∂ t ) F R + Δ ′ f 1 + I 1 ,

∂ f 2 ∂ t + u ∂ f 2 ∂ x + v ∂ f 2 ∂ y + ( w − V 2 ) ∂ f 2 ∂ z = ( ∂ f 2 ∂ t ) S B + ( ∂ f 2 ∂ t ) A C + ( ∂ f 2 ∂ t ) F R + Δ ′ f 2 + I 2 , (4)

( ∂ f 1 ∂ t ) C , ( ∂ f 1 ∂ t ) C G , ( ∂ f 1 ∂ t ) A C , ( ∂ f 1 ∂ t ) D R , ( ∂ f 1 ∂ t ) F R ―changes of the distribution functions of droplets due to microphysical condensation processes, coagulation of droplets, accretion of droplets and crystals, crushing and freezing of drops, respectively; ( ∂ f 2 ∂ t ) S B , ( ∂ f 2 ∂ t ) A C , ( ∂ f 2 ∂ t ) F R ―changes of the distribution functions of crystals due to sublimation, accretion and freezing of droplets; I 1 and I 2 ―sources of droplets and crystals.

Initial and boundary conditions for Equations (1)-(4) are as follows:

u ( r , 0 ) = u 0 ( r ) , v ( r , 0 ) = v 0 ( r ) , w ( r , 0 ) = w 0 ( r ) , θ ( r , 0 ) = θ 0 ( r ) , s ( r , 0 ) = s 0 ( r ) , f 1 ( r , m , 0 ) = f 2 ( r , m , 0 ) = f 3 ( r , m , 0 ) = 0 ,

u ( r , t ) = u 0 ( r ) , v ( r , t ) = v 0 ( r ) , w ( r , t ) = w 0 ( r ) , θ ( r , t ) = θ 0 ( r ) , s ( r , t ) = s 0 ( r ) | x = 0 , L x ; y = 0 , L y ; z = L z (5)

u ( r , t ) = v ( r , t ) = w ( r , t ) = 0 , θ ( r , t ) = θ 0 ( r ) , s ( r , t ) = s 0 ( r ) | z = 0

f 1 ( r , m , t ) = f 2 ( r , m , t ) = 0 | x = 0 , L x ; y = 0 , L y ; z = L z

∂ f 1 ( r , m , t ) / ∂ z = ∂ f 2 ( r , m , t ) / ∂ z = 0 | z = 0

Microphysical block of the model describes the processes of nucleation, condensation, coagulation of the droplets with the droplets, sublimation, accretion, freezing of drops, sedimentation of cloud particles in a gravity field, the airflow transfer [ 2 ].

The system of equations of model was solved by the method of splitting into physical processes and component splitting. In the first stages of the procedure of splitting the equations of hydro- thermodynamics are solved. One common approach to the numerical solution of the dynamics equations is based on the solution of an elliptic equation for the pressure. Firstly the system of equations of the wind and turbulent change of fields of airflows was solved. Derived fields are not adapted. To obtain the final values the Poisson equation for the pressure field is solve, we used an iterative method [ 3 ].

Let consider some of the results of studies of dynamic and microstructural parameters of convective clouds under different stratifications of the atmosphere.

The size of the spatial region in the calculations was set 40 × 40 km horizontal and 16 km in the vertical. The grid spacing in the X, Y coordinates was 500 m, along the Z―250 m. The X axis is directed eastward, Y―northward. The cloud was initiated by the set pulse at the earth surface with the overheating ΔT = 1.0˚C. The shape and size of the pulse was varied.

For the analysis of calculation results the software for three-dimensional data visualization was developed, which allows one to see the modeling object, rotate it in the vertical and horizontal plane, to change the form of representation: isosurfaces, isolines, vector field, etc. One can also quickly change the displayed value of the parameter by moving the slider in the program menu.

The first phase of research was aimed at modeling of the development of axisymmetric clouds in a calm atmosphere to analyze the quality of the model and algorithms of calculations. The data showed high accuracy of numerical algorithms developed for the model.

Hereinafter let we consider the results of calculations of dynamic and microstructure parameters of clouds in different stages.

The simulation results show that the dynamic, thermodynamic and microphysical parameters in convective clouds exert a mutual influence on each other.

During convective clouds simulations, downward streams of air around the clouds were observed. On the windward side, they can be amplified by an external wind, and on the leeward side they may be weaker due to the hydrodynamic shade in this part of the cloud.

Numerical experiments based on the upper air balloon sounding showed that the characteristics of the result clouds are sensitive to the vertical profiles of temperature and humidity. Quite powerful clouds were observed upon significant instability of the atmosphere, the tops of which consisted of crystals. When model clouds were compared to actual ones, qualitative and quantitative agreement was noticed.

The lines of airflows that run into the cloud in the maturity stage (

developed visualization tools. The figure shows several preferential flows from several directions: from the right side in the direction of cloud movement, from the rear part of the cloud, a small beam in the direction of the average wind in the atmosphere.

The numerical simulation of cloud let us to make space and time analysis of vertical and horizontal velocity of airflows, the turbulent diffusion coefficient, water content, ice content and other parameters on different stages of the cloud.

Three-dimensional non-stationary numerical model of a convective cloud with detailed description of thermodynamic and microphysical processes was developed.

The special software is developed for 3-dimensional visualization of simulation data. The software includes convenient tools for detailed analysis of the thermodynamic parameters and microstructure of the cloud.

Numerical experiments on research of formation of convective processes in the atmosphere under different stratifications of the atmosphere and the background wind were conducted. Thermodynamic parameters in a zone of thunderstorm clouds were determined in the areas with the vortex motion of the air. The main streams that feed the deep cloud in the mature stage were defined.

By means of visualization tools, regions of the formation and growth of precipitation particles were identified.

The interaction of dynamic and thermodynamic processes was analyzed in 3D. The interaction of the different physical processes is manifested in the form deformation of fields of thermodynamic parameters under the influence of dynamic processes. Cloud simulation furthers to track air flows around the cloud and the trajectories of the droplets in the cloud.

The results of numerical experiments confirm that dynamic processes significantly affect the formation of fields of thermodynamic parameters in the cloud that in turn determine the microphysical processes and the nature of the growth of precipitation particles.