The article is devoted to the development of software for calculating the hydrodynamic conditions and kinetic characteristics of granulation process in vortex devices. In the basis software Vortex Granulator© and Classification in vortex flow© original mathematical model for calculating the flow rate of gas and granules classification and separation processes of granules in a vortex granulator, kinetics of granules heating and removing moisture from the granules was put. The theoretical model, structure of software and algorithm of their work is shown. In the article an algorithm for calculation of the granulation process in the vortex granulator using the developed software is shown. Software designed in JavaFx platform. Vortex Granulator© and Classification in vortex flow© allow to conduct an optimization calculation of vortex granulator according the criteria of minimum required residence time of granules in device workspace.
A large number of chemical industries use heterogeneous processes, that take place in
«gas-liquid-solid» system and have a special place among other processes, because
the speed of occurrence is determined by the laws of heat and mass transfer in phases
that interact. These processes include various ways of granulation [
Among the variety of ways of chemical technology processes intensification working
environments flows twisting is one of the simplest and most common ways. This is
caused by the fact, that the application of swirling flows leads to improving of
heatmass transfer efficiency and temperature irregularities equalization and flows
stabilization [
Currently the software for the optimization calculation granulation process in vortex devices is not designed. This is explained by the fact that, due to insufficient knowledge of hydrodynamic view of traffic flows and kinetics of granulation in fluidized bed difficulties arise in design of vortex granulators. Vortex granulators in the most cases are patented by industrial enterprises which can conduct experimental research. The theoretical basis for calculation of such devices is limited by basic knowledge of fluid dynamics and heat and mass transfer in a classical fluid bed granulation. Calculation of vortex granulators is based on known algorithms of classical fluidized bed devices. Calculation by these techniques leads to significant error in results. A new software product that will be created is based on the original mathematical model for calculating exactly vortex granulators.
The purpose of the work is a creation of special software for calculating granulators of vortex type. One of the tasks of creating of software for system calculation of hydraulic and thermodynamic conditions for creating granules with porous structure, as well as the calculation of classification process and separation of granules in vortex granulator is to find the optimal programming language. It should allow to quickly calculate the values from the given formulas, it has a set of tools for developing client interface and visualization of values obtained in form of plots. It is also important to search for optimal programming language is the necessity to ensure application performance on different operating systems. 2
We write the Navier-Stokes equation for the motion of the real gas flow in the
diffuser (fig. 1), complementing it with flow continuity equations [
HVr −Vϕ2 =Fr − 1 ∂p +E(∆υr −Vr2 − 22 ∂Vϕ r r
∂ϕ r
V V HVϕ + r ϕ
r HVz =F z − 1 ∂p ρ ∂z ρ ∂r +E ∆V ,
z =Fϕ +E(∆Vϕ −Vϕ +
2
r
r
2
2 ∂Vr ),
∂ϕ
∂Vr + 1 ∂Vϕ + ∂Vz + Vr = 0, (2)
∂r r ∂ϕ ∂z r
where E – a coefficient of turbulent viscosity according to the Boussinesq hypothesis
[
∂ r ∂r +Vϕ ∂ r ∂ϕ
+V + 1 ∂ + 1 ∂2 r ∂r r2 ∂ϕ2 +
∂ z ∂z , ∂2 ∂z2 .
(3) (4) nents of the gas rate, respectively; α – half opening angle of the diffuser. Decision of system (1) - (2) has allowed to determine gas flow velocity components at each point of height and radius of vortex granulator, depending on its workspace design and gas flow rate.
Theoretical calculations of components of the granule rate were based on
mathematical instrument describing the hydrodynamic characteristics of the granules in the
working space of the vortex granulator [
= − +ψ ⋅ (Vϕ − Wϕ ) , (5) dτ r 8 ⋅ m
π ⋅µ gs ⋅ dgr dWz = −g +ψ ⋅ (Vz − Wz ) , dτ 8 ⋅ m where m – mass of granules, τ – time; r – the current radius of the working space of the vortex granulator; – radial, circumferential and axi
Wr ,Wϕ ,Wz al/longitudinal/perpendicular components of the velocity of the granules, respectively; g – acceleration of gravity; ψ – linear coefficient of the granule’s resistance to the gas flow; µ g – viscosity of the gas stream; dgr – diameter of the granule. To determine the granules trajectory in equations (5) we conduct the replacement given that during timeτ granule passes a way in radial Sr , circular Sϕ and vertical S z directions
ddWτr = ddτ2S2r ; ddWτϕ = dd2τS2ϕ ; ddWτz = dd2τS2z . (6) Solving system of equations (5) for the variable S in each direction occurs for some granules movement time inside the granulator.
With the use of the present system of differential equations of the granule rate in cylindrical coordinates and the numerical solution by the RungeKutta method (fig. 2) the trajectories were obtained according to the structure of the vortex granulator and properties of the granules Fig. 2. Design scheme for nu- The solid phase flow, which is modeled with Namerical determination of the vier-Stokes equations system (1) and flow contitrajectory of the granules by the nuity equation (2), gives part of amount of motion Runge-Kutta method. momentum to the granules. In the case of granules in device working volume appearing it is pulled into rotary motion by the energy of gas flow [19]. This is dispersed phase input to solid will lead to tangible change in value angular velocity component of gas flow. The gas flow circumferential velocity after the interaction with dispersed phase Vϕ' = Vϕ − ρ gr ⋅ Q gr Wϕ , ρ g Qg (7) where ρ gr – density of granule; ρ g – density of gas flow;
Based on power analysis (fig.3) calculation formulas [20] for the determination of granules classification conditions in the working space of device with a variable cross-sectional area were obtained: − gas flow velocity that meets the granules equilibrium conditions with size rgr in gas flow where rgr - radius of granule. − current location height of granule with size rgr
Z = 1, 584 − current height of location of granule with size rgr in case of mass changing (due to humidification and subsequent drying, for example in porous ammonium nitrate production): where U – moisture content of granule. − radius of lower cross-section of granulator workspace in condition of location on it granules with maximum size rmax in polydispersed system (8) (9) (10) (11) (12) (13) − heights according to fig. 3
Distribution of temperature in granule is described by the differential equation [21]: ddτ (rT (r,τ )) = a ddr22 (rT (r,τ )) , (3) where r – current value of radius, а – thermal diffusivity coefficient.
Given that U is a function of r and drying time τ, granules mass is determined by dependence
Vop = 1, 63⋅ ρ gr ⋅ g ⋅ rgr ,
ψ ⋅ ρ g tgα 2 ⋅ tgα 2 ⋅
Qg ρ g ⋅ g ⋅ rgr ψ ⋅ ρ g
.
Qg m ⋅ (1+ U ) ⋅ g ψ ⋅ ρ g ⋅ rgr
,
.
Z1 = Z − Z0 ,
Z0 = R0 / tgα .
R nπRr e− nπRmτ ∫Ueq + (U0 −Ueq ) ∑∞ 2(sin(nπ ) − nπ cos(nπ ))R sin 0 n=1 (nπ − sin(nπ ) cos(nπ ))rnπ dr (8)
π R3 R where m – diffusion coefficient; U0 – granules initial moisture content; Ueq – granules moisture content in kinetic equilibrium state.
The law of granules humidity changing in time
R 2RsinπRr e−π R2m2τ Rsin 2πRr e−4πR22mτ 0∫ U p + (U0 −U p π r − π r )dr
R (9) 3
Vortex Granulator©.
In the program Vortex Granulator© nine classes and one stylesheet were applied. Program model is presented in form of UML diagrams (fig. 4).
Main.java class is responsible for displaying the main menu and connection to other program files.
Class ControllerInput.java provides input to calculate the gas flow velocity components and the total granules and the rates of these phases.
Class Theory.java presents theoretical information about Vortex granulators, its advantages and disadvantages.
SpeedGas.java and Granuly.java classes are displayed graphically based velocity components (vertical, radial and circumferential) of gas flow and granules. Class Traektory.java displays the total rate of velocity.
Class Functional.java retains the input values of program and includes methods of forming solutions and data sets to display graphic dependencies.
Classification in vortex flow©.
Program Classification in vortex flow© has seven classes and one stylesheet were applied.
Program model is presented in form of UML diagrams (fig. 5).
Main.java class is responsible for displaying the main menu, and connects to other program files.
Vortex Granulator©.
The main class of Vortex Granulator© program is Main.java which is responsible for displaying the main menu and connection to other classes. Default program after starting first connects class, which is responsible for data entry. After entering the input parameters (fig. 6), the program checks them for validity using Error.java class. If there is an input error, the program calls AlertBox.java class responsible for the output error information. Checked data Functional.java fall into a class that is responsible for the storage and calculation.
Program has three calculation areas: components of velocity of gas flow, components of granules velocity, total velocity of gas flow and granules. When activated menu items the program respectively calculates the radial, vertical, circumferential (circular) components of gas flow rates (fig. 7) and granules (fig. 8).
Class DataInput.java retains the input values, and the program includes methods of forming solutions and data sets to display graphic dependencies.
ControllerInputSpeedGas.java class contains the class for displaying the graphical interface.
Class Graph.java responsible for the output values obtained in the form of plots. Class Excelgenerator.java generates a document electronically on the obtained values in the table format. • public Wr - calculates the radial component of velocity of granules along the radius of vortex granulator; • public RungeKuttWfi - implementing the Runge-Kutta algorithm for finding the value of the circumferential velocity component of granules; • public Wfi - calculates the circumferential component of granules velocity along the radius of vortex granulator; • public RungeKuttWz - implementing the Runge-Kutta algorithm for finding the value of the vertical component of velocity of granules; • public Wz - calculates the vertical component of velocity of granules adjustment vortex granulator workspace; • public resultSummSpeed - calculates the total gas flow velocity along the radius of vortex granulator; • public resultSummGranul - calculates the total velocity of granules along the radius of vortex granulator. • calculation of the distribution of granules on fractions in working space of granulator - shows the distribution of granules of different diameters at the height of working space of vortex granulator (fig.12). 3. Calculation of mass and kinetics of granules heating: • calculation of granules mass (fig.13); • calculation of kinetics of granules heating (distribution of temperature according to granules radius); • calculation of kinetics of granules heating (granule temperature change in time).
Fig.11. Calculation of geometry of working space of vortex granulator at a variable angle of disclosure of diffuser.
For calculation of the appropriate values, it is necessary to select a specific item, enter data and press the "Calculate" button. The program will check the correctness of the data, and then calculate the corresponding values and visualizes the values obtained in form of plots. The program Classification in vortex flow© has the ability to save the results in a spreadsheet format.
Fig.12. Calculation of distribution of granules on fractions in working space of granulator.
The main methods of calculation are follows: • public double masSuchGranuly - calculates the weight of dry granules; • public double masVologGranuly - calculates the weight of wet granules; • public double speedSuchGranuly - calculates the speed of dry granules; • public double speedVologGranuly - calculates the speed of wet granules; • public double heightVitanSuchGranuly - calculates the height of the flight of dry granules; • public double heightVitanVologGranuly - calculates the altitude moist granules; • public double radiusBottomPererizSuchGranuly - calculates the radius of bottom section of granulator for dry granules; • public double radiusBottomPererizVologGranuly - calculates the radius of lower section of granulator for wet granules; • public double simpsonLanst - implements Simpson algorithm and is designed to calculate integrals; • public double sumTemperatur - calculates the kinetics of granules heating (granules temperature distribution along the radius); • public double sumTemperaturtau - calculates the kinetics of heating of granules (granules change over time temperature). Creating of programs Vortex Granulator© and Classification in vortex flow© based on the author's mathematical model is a new stage in the construction of design algorithm of vortex granulator "theoretical calculation - industrial design". The use of these programs allows for a rational choice configuration workspace vortex granulator. Automation of hydrodynamic and thermodynamic calculation of vortex granulator allows to conduct multivariate experiment without the use of experimental and industrial installations. Such an approach in modeling the vortex granulator is cost-effective since equipment design optimization calculation is carried out on the stage before the stage of pilot implementation in production. Thus, proposed software allows to quickly calculate the hydro- and thermodynamic conditions for creating granules of porous structure, as well as the calculation of classification and separation process of granules in a vortex granulator with results visualization.
Algorithm of vortex granulator engineering calculation consists of separate blocks that are logically connected to each other (fig. 14). The feature of algorithm is presence of optimization calculation block with selection of optimal hydrodynamic vortex granulator working conditions, that make possible to remove moisture from the granules for the minimum required time. In order to automate the calculations and facilitate optimization calculation each blocks of algorithm is implemented in copyright software products. Software products have the function of the output of research results as files *.xls for further transfer to the program of creation of design and technological documentation. 6
This work was carried out under the project «Improving the efficiency of granulators and dryers with active hydrodynamic regimes for obtaining, modification and encapsulation of fertilizers», state registration No. 0116U006812. While preparing the article it was used information about author`s software products (certificates of authorship No 62692, 65140, 67472). The authors thank researchers of Processes and
Equipment of Chemical and Refining In
Fig.14. Algorithm of calculation of dustries Department, Sumy State Universiprocess of receiving porous surface ty and Department of Chemical Technololayer on granule of ammonium nitrate. gy, National University of Río Cuarto, for their valuable comments during the article preparation. 19. Artyukhov, A.E., Sklabinskyi, V.I.: Experimental and industrial implementation of porous ammonium nitrate producing process in vortex granulators.
Naukovyi Visnyk Natsionalnoho Hirnychoho Universytetu, 6, 42-48 (2013). 20. Artyukhov, A.E., Fursa, A.S., Moskalenko, K.V.: Classification and separation of granules in vortex granulators. Chemical and Petroleum Engineering, 51, Nos 5-6, 311-318 (2015). 21. Artyukhov, A.E., Sklabinskyi, V.I.: 3D Nanostructured Porous Layer of Ammonium Nitrate: Influence of the Moisturizing Method on the Layer's Structure. Journal of Nano- and Electronic Physics, 8, No 4(1), 04051(2016).