Arch. Min. Sci., Vol. 60 (2015), No 2, p. 507516 Electronic version (in color) of this paper is available: http://mining.archives.pl DOI 10.1515/amsc-2015-0033 VLADIMIR PAVLOVICH ZHUKOV*, HENRYK OTWINOWSKI**, ANTON NIKOLAEVICH BELYAKOV*, TOMASZ WYLECIAL***, VADIM EVGENEVICH MIZONOV* RÓWNANIE BOLTZMANNA W MODELOWANIU PROCESÓW PRZERÓBCZYCH The paper presents an application of the Boltzmann kinetic equation to the simultaneous modeling of multi-dimensional processes. This equation defines the evolution of the distribution of the probability density in a given phase space. In the case of a grinding process, the considered phase space is defined by the Cartesian coordinates of particle position, the components of particle velocity and the particle size. The theory of Markov processes is used in the paper to solve the Boltzmann equation for the multi-dimensional space of system states. In order to verify the presented model, research into the simultaneous comminution and movement of material in a drum ball mill was performed. The methodology developed to solve the Boltzmann equation significantly reduces the computational time, which is particularly important in the solution of multi-dimensional problems. Keywords: Boltzmann equation, mineral processing, comminution, classification, particles transport, matrix model, ball mill Równanie Boltzmanna jest podstawowym równaniem kinetycznej teorii gazów opisujcym ewolucj czstek w rozrzedzonym gazie. W równaniu tym wystpuje funkcja gstoci prawdopodobiestwa zmiennej losowej w trójwymiarowej przestrzeni fazowej (funkcja rozkladu). W artykule przedstawiono sposób wykorzystania równania Boltzmanna do analizy procesów przeróbki mechanicznej surowców mineralnych. Wynikiem tej analizy jest matematyczny model wystpujcych równoczenie procesów mielenia, klasyfikacji i transportu materialów ziarnistych. W tym przypadku równanie Boltzmanna opisuje ewolucj gstoci rozkladu ziaren wzgldem skladowych prdkoci, wspólrzdnych kartezjaskich oraz rozmiaru ziarna. W mlynie funkcja rozkladu zmienia si w wyniku rozdrabniania i ruchu ziaren, a w klasyfikatorze tylko w wyniku ruchu ziaren. W ogólnym przypadku funkcja rozkladu zaley od: czasu, ruchu ziaren, DEPARTMENT OF APPLIED MATHEMATICS, IVANOVO STATE POWER UNIVERSITY, RABFAKOVSKAYA 34, 153003 IVANOVO, RUSSIA. E-MAIL: zhukov@ispu.ru ** INSTITUTE OF THERMAL MACHINERY, CZESTOCHOWA UNIVERSITY OF TECHNOLOGY, ARMII KRAJOWEJ 21, 42-201 CZSTOCHOWA, POLAND. E-MAIL: otwinowski@imc.pcz.czest.pl *** DEPARTMENT OF INDUSTRIAL FURNACES AND ENVIRONMENTAL PROTECTION, CZESTOCHOWA UNIVERSITY OF TECHNOLOGY, ARMII KRAJOWEJ 19, 42-201 CZSTOCHOWA, POLAND. E-MAIL: wylecial@wip.pcz.pl * prdkoci ziaren i rozmiaru ziaren, który zmienia si w wyniku rozdrabniania. Uwzgldnienie zjawisk losowych wymaga wprowadzenia skladowej dyfuzyjnej do równania Boltzmanna. W artykule rozpatrzono zastosowanie równania Boltzmanna do rozdrabniania periodycznego i ciglego. W otrzymanych postaciach równania mona uwzgldni rzeczywiste warunki technologiczne, co pozwala opisa stanu ukladu podczas oddzielnych lub jednoczesnych procesów przeróbczych. Przy zaloeniu jednowymiarowoci procesów rozpatrywane zagadnienie sprowadza si do znanych przypadków, analizowanych jako oddzielne procesy. Obliczenia numeryczne wykonano metod macierzow z wykorzystaniem teorii lacuchów Markowa. Przedstawiono wyniki oblicze dla przypadku jednoczesnego rozdrabniania i ruchu ziaren w mlynie bbnowym kulowym. Analiza wyników oblicze wykazala, e przebieg ewolucji stanu ukladu ziaren jest prawidlowy. W przyszlych badaniach mona uwzgldni w równaniu Boltzmanna ksztalt ziaren, co oznacza wprowadzenie dodatkowych trzech wspólrzdnych do przestrzeni fazowej. Wspólrzdne te zwizane s ze zmian dlugoci, szerokoci i wysokoci ziarna. Slowa kluczowe: równanie Boltzmanna, procesy przeróbcze, rozdrabnianie, klasyfikacja, transport ziaren, model macierzowy, mlyn kulowy 1. Introduction The Boltzmann equation is one of the most important equations in non-equilibrium statistical mechanics. This equation describes the evolution of rarefied gas. In rarefied gas there is enough space for particles to have different velocities in one spatial volume element. In a microscopic description of a rarefied neutral gas, the gas particles are moving with a constant velocity until undergoing binary collisions. In a kinetic picture, the properties of the gas are described by the density function in the phase space f(x,v,) called the distribution function. The quantity f(x,v,)dxdv represents the number of particles in the phase-space volume of the element dxdv at time . Both x and v are three-dimensional independent variables. The distribution function satisfies the Boltzmann equation, an integro-differential equation, which describes the effect of the free flow and binary collisions between the particles (Cercignani, 1988; Kubo et al., 1991, 1991a; Huang, 1987; Kremer, 2010). The Boltzmann equation is use to analyse the thermophoresis of adhesive spherical particles immersed in fluid flows. Most theories for thermophoretic force are derived on the basis of the solution to the Boltzmann equation, or some approximations of the Boltzmann equation (Li et al., 2011). Nonlinear kinetic equations with an inflow term and rather general interactions were studied by Eibeck and Wagner (2003). A particular nonlinear kinetic equation with a gradient term was considered by Rjasanow and Wagner (2005). The Boltzmann equation allows to calculate flows of simple gases and gas mixtures in micro- and nanosize channels (Anikin et al., 2012). The kinetic approach was also applied to analyze the long time evolution of the viscous compressible weakly rarefied gas flow (Aristov & Rovenskaya, 2011). There is lack of using the Boltzmann equation to analyse of mechanical processing of mineral resources. However, it can only find Broek and co-authors' papers, in which the Boltzmann's distribution is used to modeling of the flotation process (Broek & Turno, 2005; Broek, 2010). In the work of Broek & Turno (2005) in the enrichment process it was assumed that the distribution of potential energy along the vertical axis of the separator chamber is in agreement with Boltzmann's distribution. The distribution of air bubble size in flotation cell obtained in heuristic considerations (Broek, 2010), is based upon an assumption that air flowing into the flotation cell is subject to dispersion in turbulent vortexes of the liquid and the newly formed bubble surfaces have energies whose correspond the Boltzmann's distribution. Broek and Mlynarczykowska assumed that in the pneumatic mechanical flotation machine a large air bubble is divided into smaller bubbles according to Boltzmann's law of energy distribution. 2. Multi-dimensional modeling of mineral processing The paper presents a concept for modeling the simultaneous processes of grinding, classification and transferring of granular materials based on the Boltzmann kinetic equation. This equation defines the variation in time of the probability density function of a random variable in the defined phase space. The sought function is the density of the particle distribution with r respect to: the Cartesian coordinates, velocity components and size f (r ,v ,x,), where: (z1,z2, z3) unit vector (versor) of the Cartesian coordinates, v (v1,v2, v3) vector of particle velocity com ponents, x particle size, time. The product of the function f(r ,v ,x,) and the phase space dV = dz1dz2dz3dv1dv2dv3dx express the probability of the following event: particle of a size from the elementary range (x, x + dx), staying in the position defined by the elementary coordinates (z1, z1 + dz), (z2, z2 + dz2), (z3, z3 + dz3) at the time (, + d) is moving at the velocity of the elementary components (v1, v1 + dv1), (v2, v2 + dv2), (v3, v3 + dv3). Integrating over all phase space gives the total number of particles in the granular material. In a mill, the distribution function varies as a result of particle size reduction and particle motion, while in a classifier only as a result of particle motion. In a general case, the distribution function depends on: the time, particle motion divr (v f ), · particle velocity divv (a f ) and the size variation caused by comminution f c. The Boltzmann equation in a differential form is as follows: f v1 f v2 f v3 f a1 f a2 f a3 f fc z1 z2 z3 v1 v2 v3 or f vk f ak f f c , (k 1,2,3) zk vk (2) (1) where: f(r ,v ,x,) the sought density of the particle distribution with respect to the Cartesian coordinates (z1,z2, z3), velocity components (v1,v2, v3) and particle size x; a (a1,a2, a3) particle r v acceleration; repetition of the index k in the numerator and in the denominator means summation · with respect to this index. The quantity f c on the right-hand side of the equation describes the variation of density of the particle distribution as a result of comminution of particles of the size of (x, x + dx) (which means the transition of those particles to finer classes) and the appearance of particles of the same sizes (x, x + dx), representing the product of comminution of particles of the larger size. By using the comminution model based on the mass balance of particle population (Mizonov et al., 1997), equation (2) takes the following form: f vk f ak f fS (y ) xk vk xmax fS (y )b(x,y )dy (3) where: S(y) selection function, b(x,y) breakage function, x particle size after comminution, y particle size prior to comminution, xmax maximum particle size. Further on in this paper, two cases of application of equation (3) will be considered: a) periodical comminution (bowl and roller mill, ring and ball mill) At zero particle velocities and accelerations, in the case of the periodical comminution, equation (3) will take the form of the known equation of the population mass balance: f fS (y ) xmax fS (y )b(x,y )dy (4) b) comminution by attrition In order to describe the comminution of material by attrition, it is necessary to transform equation (4). Attrition takes place by detaching fine particles from large ones and thereby reducing the particle size of the feed, which means the movement of particle along the coordinate x of its size at a velocity defined as the attrition speed: x dx dz4 v4 d d (5) In this case the particle motion is described in the four-dimensional space (z1,z2, z3, z4) (k = 1,2,3, 4). The formation of fine particles during attrition corresponds to the boundary values · of the function f c of the following equation: f vk f ak f fc , xk vk (k = 1,2,3,4) (6) Further generalization of the considered problem may result from taking into account the shape of particles. In that case, three additional coordinates, x4, x5, x6, are introduced into the phase space. These coordinates represent the change in the particle length, width and height, while the index k in equation (2) takes the values k = 1,2,...,6. 2 Dk f Consideration of the random phenomena needs to introduce the diffusion component 2 zk to the equation (1) (Huang, 1987; Bozhenko et al., 2011): f vk f ak f 2 Dk f fc , 2 zk vk zk (k =1,2,3) (7) In the obtained forms of the Boltzmann equation, specific real technological conditions can be considered, what allows the system state's evolution during either separate or simultaneous comminution, classification and transport processes to be described. Assuming one-dimensional movement or comminution processes, the problem is brought to the known cases which traditionally are solved separately. 3. Boltzmann equation solution A matrix method, using the Markov theory of chains (Berthiaux et al., 2005; Mizonov et al., 2008; Mizonov et al., 2011), was proposed to solve equation (7) numerically. Let N denotes the number of phase coordinates. Along each coordinate of process, the number of the selected cells (or steps) is ni, where the index i = 1,2,...,N denotes the ordinal number of the coordinate. The general number of the considered cells n is equal to the product N of the numbers: n = n1 · n2 · ... · nN. The vector of system state U = {Ui} with the dimension 1· n is composed of the probabilities of presence in the i-th cell, which are equal to the product of the distribution function f and the elementary volume of the phase space: Ui = fV. From the normalization condition (being equivalent to the mass conservation law), the sum of these probabilities is equal to unity U i 1 . For the description of the system evolution in time according to equation (7), a random quantity Yi is introduced into each cell, which defines possible transitions: Yi = {yi, y2, ..., ymi} (8) where: y1, y2, ..., ymi numbers of those cells of the state vector U, to which transition from the i-th cell is possible with the probability q1, q2, ..., qmi; mi number of cells, to which transition from the i-th cell is possible. These probabilities form the vector Q: Qi = {qi, q2, ..., qmi} (9) According to the normalization condition, the sum of the probabilities of all possible transitions equals 1: qm mi (10) If transition from the i-th cell is possible only to neighbouring cells, then the vector of possible transitions will take the following form: Yi i,i 1, i ni ,i ni n2 ,...,i ni n2 ...nmi 1 (11) In equation (11), the index i denotes transition to the considered cell, the indexes i ±1 transitions to the right-hand and left-hand adjacent cells, the indexes i ±n1 transition to the cells situated in the rows above and below the row with the i-numbered cell, etc. The probabilities Qi of the transitions, defined by the random quantity Yi, are determined on the basis of the physical laws applicable to the courses of the considered processes. For the one-dimensional motion of a particle at the velocity v, the probability of transition to the adjacent cell in time is determined from the following expression: qv where: x size of the cell along the x axis. The change in the particle velocity is defined by acceleration, while the probability of transition to the adjacent cell along the velocity coordinate in time is: qa a where: v size of the cell along the v axis. v (13) v x (12) The change in the particle size is connected with the comminution process. To determine the probability of particle transition to finer classes as a result of comminution, a relationship, based on the population balance, was proposed: qd x0 x1 (14) x The difference (x 0 1) denotes the change in the average particle size in time . In turn, the change in the average particle size for linear forms of the selection and breakage functions is defined by the following relationship: 1 exp x0 x1 where parameter of the selection function (Mizonov et al. 1997). The matrix operation was introduced to solve numerically the multi-dimensional equation (7): ^ U k 1 p Y ,U k (15) (16) ^ where the superscript k denotes the number of time step. The matrix operator p in equation (16) can be determined through the following identity summation of the probabilities of particle transition from all cells to the considered i-th cell of the system: U ik 1 U kj pij j 1 (17) where pij probability of transition from the j-th cell to the i-th cell. The probability pij is determined from the following expression: q j for y j1 i pij 1 0 for y j1 i where the index j1 = 1,2,...,mi. Comparing to the traditional multiplication of matrices, equation (16) allows computational time to be reduced due to the elimination of operations with zero elements. (18) 4. Results and summary According to the above algorithm the solution to the Boltzmann equation will be presented on the example of the problem of simultaneous material comminution and movement in a ball tumbling mill. The mill's z axis, the material velocity v along the mill axis and the particle size x were assumed as the coordinates; hence the number of coordinates, along which the process evolution is considered, is equal to three (N = 3). The three cells n1 = 3, n2 = 3, n3 = 3 were chosen along each coordinate. The considered phase space is presented in Figure 1. Fig. 1. The phase space in the problem of the simultaneous grinding and movement of material particles in the ball mill (the numbers denote the creation method of the state vector) Numeration of the cells in Figure 1 represents the method of creating the vector of system state U, which contains n = n1 · n2 · n3 = 27 elements. Comminution is regarded as the reduction of particle sizes or their transition to smaller-size cells with a specified probability d . At the initial moment, a unit sample of the coarse particles is fed into the mill. Considering the assumed order of cell numerations, the initial state can be described by the vector of system state in the following form: 1 for i 3 U i0 0 for i 3 (19) Each cell, situated inside the separated phase space, may border on 6 adjacent cells. If the cell is present at the phase space boundary, then the transitions through the boundary are defined by the adequate boundary conditions. In the considered example, an isolated phase volume was selected, from which no particles come out. If the material may only pass to adjacent cells during evolution, then the random quantity Y is expressed applying the following table: Y Q i 1d av i1 d i+1 0 i n1 0 i + n1 a i n1n2 0 i + n1n2 v The evolution of system state in time, calculated using equation (16) for the initial conditions (19) in the totally isolated phase space (Fig. 1), is presented in Figure 2. In time, the elementary sample of coarse particles passes to the cell with the minimum particle size, the maximum velocity and the maximum coordinate along the mill axis. The change in the contents of three classes in time is presented in Figure 3. The mass fraction of coarse class decreases and the mass fraction of fine class increases. The change in the contents at the outlet of the mill increases with the increasing of the convective mass transfer rate (Fig. 4). Fig. 2. An evolution of the state vector during material comminution and particle motion at a constant acceleration in the tumbling ball mill: the material passing to the cell at the maximum velocity, the minimum particle size and the maximum coordinate along the mill's axis. Designations of the axes: i cell number consistent with the numeration shown in Figure 1, TS time step, Ui value of the vector of system state Fig. 3. Change in the contents of a coarse (1), medium (2) and fine (3) classes in time Fig. 4. Change in the contents at the outlet of the mill for different convective mass transfer rates v: 1 v = 0.05, 2 v = 0.1, 3 v = 0.3 The model presented in the paper enables multi-dimensional simultaneous comminution, classification and transport processes to be analysed. The developed methodology to solve the Boltzmann equation, using operations only on non-zero elements, significantly reduces the computational time, which is particularly important in the solution of multi-dimensional problems.
Archives of Mining Sciences – de Gruyter
Published: Jun 1, 2015
You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!
Read and print from thousands of top scholarly journals.
Already have an account? Log in
Bookmark this article. You can see your Bookmarks on your DeepDyve Library.
To save an article, log in first, or sign up for a DeepDyve account if you don’t already have one.
Copy and paste the desired citation format or use the link below to download a file formatted for EndNote
All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.