Modeling and Simulation of Particulate Processes

Particulate processes can be modeled by means of populations balances. This is an important class of nonlinear partial differential equations with many applications in chemical and biochemical engineering. Major challenges are multidimensional problems, coupling with nonideal flow fields and feedbac...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Электронное моделирование
Datum:2016
Hauptverfasser: Kienle, A., Palis, S., Mangold, M., Durr, R.
Format: Artikel
Sprache:English
Veröffentlicht: Інститут проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України 2016
Schlagworte:
Online Zugang:https://nasplib.isofts.kiev.ua/handle/123456789/115839
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Modeling and Simulation of Particulate Processes / A. Kienle, S. Palis, M. Mangold, R. Durr // Электронное моделирование. — 2016. — Т. 38, № 5. — С. 23-33. — Бібліогр.: 19 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id nasplib_isofts_kiev_ua-123456789-115839
record_format dspace
spelling Kienle, A.
Palis, S.
Mangold, M.
Durr, R.
2017-04-14T09:18:32Z
2017-04-14T09:18:32Z
2016
Modeling and Simulation of Particulate Processes / A. Kienle, S. Palis, M. Mangold, R. Durr // Электронное моделирование. — 2016. — Т. 38, № 5. — С. 23-33. — Бібліогр.: 19 назв. — англ.
0204-3572
https://nasplib.isofts.kiev.ua/handle/123456789/115839
Particulate processes can be modeled by means of populations balances. This is an important class of nonlinear partial differential equations with many applications in chemical and biochemical engineering. Major challenges are multidimensional problems, coupling with nonideal flow fields and feedback control. Possible solution approaches to these problems are presented and illustrated with different types of process applications including fluidized bed spray granulation, crystallization and influenza vaccine production processes.
Процессы в макрочастицах можно моделировать, используя популяционный баланс. Он представляет собой важный класс нелинейных дифференциальных уравнений в частных производных и широко применяются в химической и биохимической инженерии. Основными проблемами при этом являются многомерные задачи, взаимосвязь с неидеальными полями течения и управление с обратными связями. В работе представлены возможные подходы к решению этих задач на примере различных процессов, таких как грануляция в кипящем слое, кристаллизация и процессы производства вакцин от гриппа.
Процеси обробки часток можна моделювати, використовуючи популяційний баланс. Він є важливим класом нелінійних диференціальних рівнянь з частинними похідними та широко застосовується в хімічній та біохімічної інженерії. Основними проблемами тут є багатовимірні задачі, взаємозв’язок з неідеальними полями течії та керування зі зворотними зв’язками. У роботі представлені можливі підходи до вирішення цих задач на приклад і різних процесів, таких як грануляція в киплячому шарі, кристалізація і процеси виробництва вакцин від грипу.
The financial support of the German Science Foundation within the priority program SPP 1679, the Federal Ministry for Education and Research within the projects SimParTurS (grant 03KIPAA2) and CellSys (grant 0316189A) as well as the German Competence Network of Chemical Engineering Pro3 is greatly acknowledged. Special thanks goes to Prof. Svjatnyj and his students for many years of excellent cooperation.
en
Інститут проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України
Электронное моделирование
Математическое моделирование и вычислительные методы
Modeling and Simulation of Particulate Processes
Article
published earlier
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
title Modeling and Simulation of Particulate Processes
spellingShingle Modeling and Simulation of Particulate Processes
Kienle, A.
Palis, S.
Mangold, M.
Durr, R.
Математическое моделирование и вычислительные методы
title_short Modeling and Simulation of Particulate Processes
title_full Modeling and Simulation of Particulate Processes
title_fullStr Modeling and Simulation of Particulate Processes
title_full_unstemmed Modeling and Simulation of Particulate Processes
title_sort modeling and simulation of particulate processes
author Kienle, A.
Palis, S.
Mangold, M.
Durr, R.
author_facet Kienle, A.
Palis, S.
Mangold, M.
Durr, R.
topic Математическое моделирование и вычислительные методы
topic_facet Математическое моделирование и вычислительные методы
publishDate 2016
language English
container_title Электронное моделирование
publisher Інститут проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України
format Article
description Particulate processes can be modeled by means of populations balances. This is an important class of nonlinear partial differential equations with many applications in chemical and biochemical engineering. Major challenges are multidimensional problems, coupling with nonideal flow fields and feedback control. Possible solution approaches to these problems are presented and illustrated with different types of process applications including fluidized bed spray granulation, crystallization and influenza vaccine production processes. Процессы в макрочастицах можно моделировать, используя популяционный баланс. Он представляет собой важный класс нелинейных дифференциальных уравнений в частных производных и широко применяются в химической и биохимической инженерии. Основными проблемами при этом являются многомерные задачи, взаимосвязь с неидеальными полями течения и управление с обратными связями. В работе представлены возможные подходы к решению этих задач на примере различных процессов, таких как грануляция в кипящем слое, кристаллизация и процессы производства вакцин от гриппа. Процеси обробки часток можна моделювати, використовуючи популяційний баланс. Він є важливим класом нелінійних диференціальних рівнянь з частинними похідними та широко застосовується в хімічній та біохімічної інженерії. Основними проблемами тут є багатовимірні задачі, взаємозв’язок з неідеальними полями течії та керування зі зворотними зв’язками. У роботі представлені можливі підходи до вирішення цих задач на приклад і різних процесів, таких як грануляція в киплячому шарі, кристалізація і процеси виробництва вакцин від грипу.
issn 0204-3572
url https://nasplib.isofts.kiev.ua/handle/123456789/115839
citation_txt Modeling and Simulation of Particulate Processes / A. Kienle, S. Palis, M. Mangold, R. Durr // Электронное моделирование. — 2016. — Т. 38, № 5. — С. 23-33. — Бібліогр.: 19 назв. — англ.
work_keys_str_mv AT kienlea modelingandsimulationofparticulateprocesses
AT paliss modelingandsimulationofparticulateprocesses
AT mangoldm modelingandsimulationofparticulateprocesses
AT durrr modelingandsimulationofparticulateprocesses
first_indexed 2025-11-24T23:38:57Z
last_indexed 2025-11-24T23:38:57Z
_version_ 1850500655954788352
fulltext A. Kienle 1,2 , Prof. Dr.-Ing., S. Palis 2 , Jun. Prof. Dr.-Ing., M. Mangold 1 , Prof. Dr.-Ing., R. D��urr 2 , Dipl.-Ing. 1 Max Planck Institute for Dynamics of Complex Technical Systems (Sandtorstr. 1, 39106 Magdeburg, Germany), 2 Otto von Guericke Universt��at (Universit��atsplatz 2, 39106 Magdeburg, Germany) (Tel. +49 391 67 58523, e-mail: achim.kienle@ovgu.de) Modeling and Simulation of Particulate Processes Particulate processes can be modeled by means of populations balances. This is an important class of nonlinear partial differential equations with many applications in chemical and biochemi- cal engineering. Major challenges are multidimensional problems, coupling with nonideal flow fields and feedback control. Possible solution approaches to these problems are presented and il- lustrated with different types of process applications including fluidized bed spray granulation, crystallization and influenza vaccine production processes. Ïðîöåññû â ìàêðî÷àñòèöàõ ìîæíî ìîäåëèðîâàòü, èñïîëüçóÿ ïîïóëÿöèîííûé áàëàíñ.Îí ïðåäñòàâëÿåò ñîáîé âàæíûé êëàññ íåëèíåéíûõ äèôôåðåíöèàëüíûõ óðàâíåíèé â ÷àñòíûõ ïðîèçâîäíûõ è øèðîêî ïðèìåíÿåüñÿ â õèìè÷åñêîé è áèîõèìè÷åñêîé èíæåíåðèè. Îñíîâ- íûìè ïðîáëåìàìè ïðè ýòîì ÿâëÿþòñÿ ìíîãîìåðíûå çàäà÷è, âçàèìîñâÿçü ñ íåèäåàëüíûìè ïîëÿìè òå÷åíèÿ è óïðàâëåíèå ñ îáðàòíûìè ñâÿçÿìè.  ðàáîòå ïðåäñòàâëåííû âîçìîæíûå ïîäõîäû ê ðåøåíèþ ýòèõ çàäà÷ íà ïðèìåðå ðàçëè÷íûõ ïðîöåññîâ, òàêèõ êàê ãðàíóëÿöèÿ â êèïÿùåì ñëîå, êðèñòàëèçàöèÿ è ïðîöåññû ïðîèçâîäñòâà âàêöèí îò ãðèïïà. K e y w o r d s: partial differential equations, population balances, control, model reduction, proper orthogonal decomposition, direct quadrature method of moments. Introduction. Particulate products like crystals, granules, and powders play a major role in process industries. Typical examples are pharmaceuticals, deter- gents, pigments, polymers etc. They represent about 60 % of the produced value in the chemical industry. Typical production processes comprise crystallization, granulation, and polymerization among others. Function and effectiveness of particulate products often depend on particle properties - such as size, porosity, morphology or composition. Main objective of our research in this field of appli- cation is to devise new methods and tools for modeling, simulation and control of particulate processes aiming at the directed adjustment of desired product properties. This is a challenging issue due to nonuniformity of particle systems, where particles differ with respect to individual properties, and product properties are ISSN 0204–3572. Ýëåêòðîí. ìîäåëèðîâàíèå. 2016. Ò. 38. ¹ 5 23 � A. Kienle, S. Palis, M. Mangold, R. D��urr, 2016 represented by the collective behavior of the particle population. From the theoretical point of view particulate processes belong to a special class of distributed parameter systems, so called population balance sys- tems. They are described by nonlinear partial differential equations (PDE) often coupled to integro differential equations describing the surrounding medium (see e.g. [1]). The paper will address main challenges for modeling, simulation and control of par- ticulate processes and present recent solution approaches developed in our group at the Max Planck Institute and the Otto von Guericke University in Magdeburg. Theoretical concepts will be illustrated with practical application examples in- cluding granulation, crystallization and vaccine production in cell cultures. Modeling of particulate processes. The generic structure of particulate processes is illustrated in Fig. 1 for a crystallization process. It consists of a con- tinuous phase — also sometimes called the medium — and a disperse particle phase, which interacts with the continuous phase through an exchange of mass, energy and momentum. Particles are characterized by their specific properties like particle size, morphology or composition, for example. These properties are summarized in a vector of internal coordinates x. Typically, particles are nonuniform (see also Fig. 1) giving rise to a property distribution like a particle size distribution, for example. In case of spatial inhomogeneity of the continuous phase an additional vector of external space coordinates r is required for the mathematical desc- ription of the system. Kinetic processes in the particle phase may comprise nu- cleation, i.e. formation of new particles, particle growth, aggregation and brea- kage processes. Microscopic models based on discrete element methods are computatio- nally very expensive and therefore currently not well suited for model based pro- cess design and process control purposes [2]. Therefore, focus in this paper is on a more macroscopic approach using population balances [1]. In general form the population balance equation reads � �t n t n n h t dV dV i ( , , ) ( � ) ( � ) ( , , )x r r x x x rr x x x �� �� � � �� �, �� �, ,n ni i iin out , where n t( , , )x r is a number density, which depends on external and internal variables and evolves over time. The equation consists of an accumulation term, A. Kienle, S. Palis, M. Mangold, R. D��urr 24 ISSN 0204–3572. Electronic Modeling. 2016. V. 38. ¹ 5 Fig. 1. Crystallization process convection in the direction of the external (space) variables r and convection in the direction of the internal state variables x. The latter corresponds to the kinetic processes contributing to the continuous change of particle properties like growth processes, for example. The convolution type of integral on the right hand side represents aggregation and breakage processes. Furthermore, material transport across the system boundaries is taken into account. The population ba- lance is coupled with corresponding material, energy and momentum balances of the continuous phase, including spatial gradients in r and integral interaction with the particle phase. Specific examples will be given later in this paper. Major challenges to be discussed subsequently comprise: • problems with multiple internal coordinates; • coupling with nonideal flow fields leading to problems with internal and external coordinates; • process control. Nonlinear dynamics and control of fluidized bed spray granulation processes. First, focus is on process control. As an application example continu- ous fluidized bed spray granulation as illustrated in Fig. 2 is considered. The pro- cess consists of a granulation chamber, where the particles are fluidized with a stream of hot air. From the top a liquid suspension is sprayed onto the particles giving rise to particle growth. Particles are removed continuously from the granu- lation chamber and classified into a product fraction, fractions of undersized and oversized particles. Oversized particles are ground with a mill and returned to- Modeling and Simulation of Particulate Processes ISSN 0204–3572. Ýëåêòðîí. ìîäåëèðîâàíèå. 2016. Ò. 38. ¹ 5 25 Fig. 2. Fluidized bed spray granulation process with external product classification gether with the undersized particles to the granulation chamber giving rise to new nuclei for particle growth. On the one hand, this mode of operation is very economic since it avoids the supply of external nuclei. On the other hand, it may lead to dynamic instability in the from of self sustained oscillations [3]. This im- plies oscillating product properties like particle size L, which is usually not ac- ceptable from the practical point of view. Therefore, a new method was develo- ped for stabilization by means of feedback control. The model equations are given by: � � � � � n t G n L n n nP O M M� � � � �� � � ( ) , � ( )( ( )) �n T L T L nP � �2 11 , � ( ) �n T L nO � 1 , (1) with growth rate G m m L ndL e e( � ) � � � � 2 0 2�� , A. Kienle, S. Palis, M. Mangold, R. D��urr 26 ISSN 0204–3572. Electronic Modeling. 2016. V. 38. ¹ 5 Fig. 3. Open loop simulation results: a — stable (coarse grinding with� M = 0.9 mm); b — un- stable (fine grinding with � M = 0.7 mm) which assumes that the injected liquid suspension �me is equally distributed over the available total particle surface. The withdrawal �n follows from a total mate- rial balance assuming constant bed mass. T L1( ) and T L2( ) are the sigmoid classi- fying functions of the sieves. For the nomenclature of the flow rates � , � , �n n nP O we refer to Fig. 2. Although the model is rather simple it contains all necessary ingredients to reproduce the nonlinear oscillations. It further shows the important influence of the mill grade � M on process stability as illustrated in Fig. 3. For coarse milling the process is settling down to a stable steady state. In contrast to this, the fine milling triggers instability [4]. This further indicates that the mill grade is a suit- able handle to stabilize the process by means of feedback control. For the present class of processes a new control method was developed. It is a Lyapunov type of approach based directly on the nonlinear PDE (1) and uses the mill grade as a lumped handle. Stability is proven in the integral sense, for the norm of the difference of the overall particle surface wrto to its given reference. It should be noted that this quantity is not a norm of the number density n L t( , ) but a weaker error measure called a discrepancy. It was shown that this integral stability implies pointwise convergence of n L t( , ), if the corresponding zero dy- namics are stable as shown in Fig. 4 for the moments � i iL n L t dL� � � 0 ( , ) , i = 0. 2 and the number density n L t( , ). The figure also shows the control error e and the Modeling and Simulation of Particulate Processes ISSN 0204–3572. Ýëåêòðîí. ìîäåëèðîâàíèå. 2016. Ò. 38. ¹ 5 27 Fig. 4. Closed loop simulation results: without control; ···· with control manipulated variable � M . For unstable zero dynamics an additional robust pa- rallel compensator can be used for stabilization. The details are given in [5]. The method is rather simple and intuitive and applies to a large class of problems. Application was also demonstrated for fluidized bed spray granula- tion with internal product classification [6], spray coating in a Wurster coater [7] and continuous crystallization [8]. Experimental validation is in progress. Nonlinear model reduction for particulate processes with nonideal flow fields. Models are becoming more involved, when spatial inhomogeneities play a major role leading to problems with internal and external coordinates. A typi- cal example is the crystallization process in a tube described in [9]. Model equa- tions are given by the population balance coupled to the Navier Stokes equation, i.e. the momentum balance, as well as energy and material balances of the con- tinuous phase. Internal coordinate is again the particle size L which is comple- mented with two external coordinates summarized in the vector r to describe the spatial position in the crystallizer. These kinds of models are usually too complex for online applications in the process control or process design. Therefore, nonlinear model reduction is inevitable. Proper orthogonal decomposition was used for the first time for this class of processes in [9, 10]. The main idea is to approximate the solution of the coupled system of equations by means of a series expansion n L t t L i n i i( , , ) ( ) ( , )r r� � 1 � � with basis functions � i which depend on external and internal coordinates only and time dependent weights� i . The basis functions were generated from offline simulations with the full reference model using a finite element approach. The reduced model equations for the time dependent weights were obtained by Galerkin projection. Best point interpolation was used for the approximation of nonlinearities to avoid expensive online evaluation of multidimensional inte- grals during runtime. With this approach the system order could be reduced by four orders of magnitude compared to the discretized PDE reference model. The computation time for online application could be reduced by three orders of magnitude, once expensive offline calculations are completed. Recent work is concerned with automatic generation of reduced order mo- dels by means of proper orthogonal decomposition [11]. As an application a novel type of fluidized bed crystallizer as described in [12, 13] is considered among others. Multidimensional modeling of influenza vaccine production processes. The third application example is concerned with biotechnological processes for A. Kienle, S. Palis, M. Mangold, R. D��urr 28 ISSN 0204–3572. Electronic Modeling. 2016. V. 38. ¹ 5 the production of influenza vaccines. Virus influenza is a serious disease killing many people every year. Vac- cination is currently the best option to protect against infection. Since influ- enza viruses are changing continu- ously, flexible production processes are required. Therefore production in mammalian cell cultures is becoming more and more important [14]. In this type of process, first a culture of mam- malian cells is grown. At the begin- ning of the production process the cell culture is infected with virus seeds. Typi- cally, in technical production processes, the amount of virus is low compared to the number of cells in the culture. Therefore in the beginning of the process only a part of the cells is infected. The virus replicates inside the infected cells and new virus is released to the medium and thereby infects step by step also the other cells. With time, active virus will also degrade and exhausted cells will die. The production process is stopped when the virus yield goes through a maxi- mum. Then the bioreactor is harvested and the vaccine can be produced from the virus yield. As indicated by measurements with flow cytometry the cell to cell variabi- lity is a major feature of this process which is addressed here with distributed population balance models [15]. For this purpose two different approaches are available: 1. Top down approach using global kinetics for virus infection, virus repli- cation and virus release, which are fitted to experimental data using an inverse problem approach [15, 16]. 2. Bottom up approach starting from detailed single cell kinetics. The latter accounts in more detail for available biological knowledge and therefore provides more insight into the underlying biological processes and their interaction. Dynamics on the single cell level are described by a vector of physiologically relevant components such as a viral genome or viral proteins, which change over time. Since the amount of these quantities differs from cell to cell, these state variables on the single cell level translate into internal coordi- nates on the population level leading to multidimensional population balance equations. External coordinates are neglected here due to perfect mixing of the medium in the bioreactor. Numerical solution of multidimensional population balance models of cel- lular systems can be very challenging, if N the number of internal coordinates is Modeling and Simulation of Particulate Processes ISSN 0204–3572. Ýëåêòðîí. ìîäåëèðîâàíèå. 2016. Ò. 38. ¹ 5 29 Infection Virus release gen v1 v2 tem str Fig. 5. Scheme for viral replication on the single cell level large. Standard discretization techniques using finite differences, finite volumes or finite elements are typically limited to N in the order of 5. Monte Carlo simu- lation which is another popular solution method in this field is typically limited to N in the order of 10. For higher dimensional problems a new quadrature method of moments was developed. It combines monomial cubatures with the method of characteristics [17, 18]. With this approach the computational costs increase only polynomially with N, in the simplest case even only linearly. For large N, this is a drastic reduction in computational costs compared to standard Gaussian cubatures, where the computational costs increase exponentially with N. In the remainder, application is demonstrated for a relatively simple five di- mensional problem, which was adopted from [19] and modified slightly. The replication kinetics are illustrated in Fig. 5. Relevant species are the viral ge- nome gen, the template of the viral genome tem, which is used for synthesizing new genome and structural proteins str. Generation of template and structural proteins is catalyzed by viral enzymes v1 and v2. Viral genome and structural proteins are merged to build new viruses which are released to the medium. The model equations are. • Infected cells � � i t t i t k U t V t I kc c c ic ( , ) { � ( , )} ( ) ( ) ( ) x x x xx� �� � �inf cd, ( ) ( , )x xi tc , (2) where �x in the first term on the r.h.s. corresponds to the reaction rates on the sin- gle cell level according to � [ ][ ] [ ] [ ] [ ][ ] [ ][ x � � � � k v gen k tem k tem k v gen k gen 1 1 6 3 1 1 5 str k gen k tem v k str k gen str f f v v ] [ ] [ ][ ] [ ] [ ][ ] � � � 7 2 2 4 5 1 2 � � � � � � � � � � � � � � � � . The second term on the r.h.s of equation (2) corrpesponds to infection and is proportional to the number of uninfected cells and the number of free viruses. The third term corresponds to cell death. • Uninfected cells dU t t k U t V t k U t k U tc c U c U cc c ( ) ( ) ( ) ( ) ( ) � � � � �inf gro, cd, . In contrast to the infected cells, uninfected cells Uc and free virus V are de- scribed by ordinary differential equations. The number of uninfected cells de- creases due to infection and cell death. A. Kienle, S. Palis, M. Mangold, R. D��urr 30 ISSN 0204–3572. Electronic Modeling. 2016. V. 38. ¹ 5 • Free virus in the medium dV t t r i t d k U t V t k V tc c V ( ) ( ) ( , ) ( ) ( ) ( � � � �� x x x xrel inf deg, ). The number of free viruses increases due to release of newly formed virus from infected cells and decreases due to infection and virus degradation. Results for the first two moments of viral genome and template are shown in Fig. 6 for different types of cubatures. The reference solution is obtained by Monte Carlo sampling with 104 samples. Compared to a monomial and a Gaus- sian cubature with 66 and 243 abscissas, respectively. Both show similar perfor- mance with some deviation from the reference solution. By far the best match was obtained with a monomial cubature with mixed Gaussian model densities with 121 abscissas. Currently the approach is applied to a detailed model of the single cell kine- tics with about 30 state variables translating into 30 internal coordinates. The model aims at identifying targets or genetic modification of the cell line to im- Modeling and Simulation of Particulate Processes ISSN 0204–3572. Ýëåêòðîí. ìîäåëèðîâàíèå. 2016. Ò. 38. ¹ 5 31 Fig. 6. Simulation results for the viral replication processes in Fig. 5: m t x i t dl k k l c, ( ) ( ) ( , )� � x x x x; — Monomial (N = 66); Gauss (N = 243); Monomial+GMD (N = 121); � Monte-Carlo (N = 104); prove influenza vaccine production within a joint research project with partners from academia and industry. Conclusions. It was shown that population balance systems are an impor- tant class of distributed parameter systems with many applications in chemical and biochemical engineering. Modeling, simulation and control of population balance systems turns out to be challenging, in particular for multidimensional problems with multiple internal and/or external coordinates. Several solution ap- proaches were presented which were developed in our research group at the Max Planck Institute and the Otto von Guericke University in Magdeburg during the last years. Nevertheless, there are still many opportunities for productive re- search in this field. Acknowledgement. The financial support of the German Science Founda- tion within the priority program SPP 1679, the Federal Ministry for Education and Research within the projects SimParTurS (grant 03KIPAA2) and CellSys (grant 0316189A) as well as the German Competence Network of Chemical En- gineering Pro3 is greatly acknowledged. Special thanks goes to Prof. Svjatnyj and his students for many years of excellent cooperation. Ïðîöåñè îáðîáêè ÷àñòîê ìîæíà ìîäåëþâàòè, âèêîðèñòîâóþ÷è ïîïóëÿö³éíèé áàëàíñ. Â³í º âàæëèâèì êëàñîì íåë³í³éíèõ äèôåðåíö³àëüíèõ ð³âíÿíü ç ÷àñòèííèìè ïîõ³äíèìè òà øèðîêî çàñòîñîâóºòüñÿ â õ³ì³÷í³é òà á³îõ³ì³÷íî¿ ³íæåíåð³¿. Îñíîâíèìè ïðîáëåìàìè òóò º áàãàòîâèì³ðí³ çàäà÷³, âçàºìîçâ’ÿçîê ç íå³äåàëüíèìè ïîëÿìè òå÷³¿ òà êåðóâàííÿ ç³ çâîðîò- íèìè çâ’ÿçêàìè. Ó ðîáîò³ ïðåäñòàâëåí³ ìîæëèâ³ ï³äõîäè äî âèð³øåííÿ öèõ çàäà÷ íà ïðèê- ëàä³ ð³çíèõ ïðîöåñ³â, òàêèõ ÿê ãðàíóëÿö³ÿ â êèïëÿ÷îìó øàð³, êðèñòàë³çàö³ÿ ³ ïðîöåñè âèðîáíèöòâà âàêöèí â³ä ãðèïó. REFERENCES 1. Ramkrishna, D. (2000), Population balances: Theory and applications to particulate systems in engineering, Academic Press, New York, USA. 2. Deen, N.G., Van Sint Annaland, M., Van der Hoef, M.A. and Kuipers, J.A.M. (2007), Re- view of discrete particle modeling of fluidized beds, Chem. Engng. Sci., Vol. 62, pp. 28-44. 3. Heinrich, S., Peglow, M., Ihlow, M., Henneberg, M. and L. M�rl, L. (2002), Analysis of start-up process in continuous fluidized bed spray granulation by population balance model- ing, Chem. Engng. Sci., Vol. 57, pp. 4369-4390. 4. Radichkov, R., M�ller, T., Kienle, A., Heinrich, S., Peglow, M. and M�rl, L. (2006), A nu- merical bifurcation analysis of fluidized bed spray granulation with external classification, Chem. Engng. Proc., Vol. 45, pp. 826-837. 5. Palis, S. and Kienle, A. (2014), Discrepancy based control of particulate processes, J. Proc. Contr., Vol. 24, pp. 33-46. 6. Palis, S. and Kienle, A. Discrepancy based control of continuous fluidized bed spray granu- lation with internal product classification, Proc. 8th IFAC International Symposium on Ad- vanced Control of Chemical Processes, Singapore, July 10-13, 2012, pp. 756-761. 7. Palis, S., B�ck, A. and Kienle, A. (2013), Discrepancy based control of systems of popula- tion balances, Proc. 1st IFAC Workshop on Control of Systems Modeled by Partial Differen- tial Equations, Paris, September 25-27, 2013, pp. 172-176. A. Kienle, S. Palis, M. Mangold, R. D��urr 32 ISSN 0204–3572. Electronic Modeling. 2016. V. 38. ¹ 5 8. Palis, S. and Kienle, A. (2012), Diskrepanz basierte Regelung der kontinuierlichen Fl�s- sigkristallisation, AT-Automatisierungstechnik, Vol. 60, pp. 145-154. 9. Krasnyk, M., Mangold, M., Ganesan, S. and L. Tobiska, L. (2012), Numerical reduction of a crystallizer model with internal and external coordinates by proper orthogonal decomposi- tion, Chem. Engng. Sci., Vol. 70, pp. 77-86. 10. Krasnyk, M., Mangold, M. and Kienle, M. (2010), Extensions of the POD model reduction to multi-parameter domains, Chem. Engng. Sci., Vol. 65, pp. 6238-6246. 11. Khlopov, D. and Mangold, M. (2015), Automatic model reduction of linear population ba- lance models by proper orthogonal decomposition, Proc. Vienna Conference on Mathemati- cal Modeling, Vienna, February 18-20, DOI 10.1016/j.ifacol.2015.05.019, 2015. 12. Mangold, M., Feng, L.H., Khlopov, D., Palis, S., Benner, P., Binev, D. and Seidel-Mor- genstern, A. (2015), Nonlinear model reduction of a continuous fluidized bed crystallizer, J. Comp. Appl. Math., Vol. 89, pp. 253-266. 13. Mangold, M., Khlopov, D., Danker, G., Palis, S., Sviatnyi, V. and Kienle, A. (2014), Deve- lopment and nonlinear analysis of dynamic plant models in ProMoT/DIANA, Chemie- Ing.-Techn., Vol. 86, pp. 1-12. 14. Genzel, Y. and Reichl, U. (2009), Continuous cell lines as a production system for influenza vaccines, Expert Rev. Vaccines, Vol. 8, pp. 1681-1692. 15. M�ller, T., D�rr, R., Isken, B., Schulze-Horsel, J., Reichl, U. and Kienle, A. (2013), Distri- buted modeling of human influenza a virus-host cell interactions during vaccine production, Biotechnol. Bioengng., Vol. 110, pp. 2252-2266. 16. D�rr, R., M�ller, T., Isken, B., Schulze-Horsel, J., Reichl, U. and Kienle, A. (2012), Distri- buted modeling and parameter estimation of influenza virus replication during vaccine pro- duction, Proc. Vienna Conference on Mathematical Modeling, Vienna, February 15-17, 2012. 17. D�rr, R. and Kienle, A. (2014), An efficient method for calculating the moments of multi- di- mensional growth processes in population balance systems, Can. J. Chem. Eng., Vol. 92, pp. 2088-2097. 18. D�rr, R., M�ller, T. and Kienle, A. (2015), Efficient DQMOM for multivariate population balance equations and application to virus replication in cell cultures, Proc. Vienna Con- ference on Mathematical Modeling, Vienna, February 18-20, DOI 10.1016/j.ifacol.2015. 05.045, 2015. 19. Haseltine, E.L., Yin, J. and Rawlings, J.B. (2005), Dynamics of viral infections: incorporating both the intracellular and extracellular levels, Comput. Chem. Engng., Vol. 29, pp. 675-686. Received 13.06.16 Modeling and Simulation of Particulate Processes ISSN 0204–3572. Ýëåêòðîí. ìîäåëèðîâàíèå. 2016. Ò. 38. ¹ 5 33