Populationsbilanzgleichungen (PBE) modellieren die Zeitentwicklung eines Systems von Entitäten, die unabhängig von ihrer eigentlichen Natur ”Partikel” genannt werden. Der Zustand des Systems ändert sich durch Entstehung und Vernichtung von Partikeln sowie durch Modifikationen bestehender Partikel. In klassischen mit PBE modellierten Systemen handelt es sich bei den Entitäten um physikalische Partikel wie Ruß, Kristalle oder Wassertröpfchen. In solchen Zusammenhängen beschreibt die PBE den Nettoeffekt (Bilanz) von Ereignissen wie Nukleation, Auflösung, Koagulation und Entfernung von Partikeln durch Aussieben. Populationsbilanzgleichungen sind jedoch nicht auf solche physikalischen Partikel beschränkt. Sie können zum Beispiel auch für Fragestellungen aus der Biologie genutzt werden. Dann sind die ”Partikel” Gene oder Zellen, oder sogar einzelne Individuen einer Herde von Tieren.

Wenn die Interaktion der Partikel untereinander, beispielsweise durch Koagulation, in einer PBE berücksichtigt wird, so stellt sie sich als Integro-Differentialgleichung dar, was spezifische Schwierigkeiten der Analysis und Numerik mit sich bringt. Noch herausforderdernder sind Systeme, in die auch die Interaktion von Partikeln mit ihrer Umgebung einfließt. Im Fall physikalischer Partikel beinhaltet dies Phänomene wie Partikelwachstum, Transport oder Oberflächenreaktionen. Die Umgebung der Partikel kann selbst wieder gewissen Differentialgleichungen unterworfen sein. Werden jene mit der Populationsbilanzgleichung für die Partikel kombiniert, so erhält man ein Populationsbilanzsystem. Pobpulationsbilanzsysteme sind ein wesentliches Forschungsfeld am WIAS.

Insbesondere für klassische Fluid-Partikel-Systeme wie sie im Chemieingenieurwesen auftreten, werden am WIAS Lösungsalgorithmen entwickelt, implementiert und analysiert. Außerdem wird die Qualität solcher Algorithmen in Bezug auf Stabilität, Rechenzeit und Anwendbarkeit bewertet. Durch Anwendung von PBS-Methoden auf die Modellierung und Simulation wissenschaftlicher Experimente überbrückt das Institut die Lücke zwischen Theorie und Anwendung von PBS. Zu diesem Zweck wurden bereits einige, mit nationalen und internationalen Partnern entwickelte, Forschungs-Codes erfolgreich angewandt: ParMooN, sein Vorgänger MooNMD, und der Partikel-Löser Brush.



Abbildung 1: Drei Momentaufnahmen der Dichteverteilung von Wassertröpfchen am Auslass einer turbulenten Kanalströmung. Numerische Simulation mit PBS. Tröpfchenwachstum und -transport in Wolken ist ein praktisches Anwendungsfeld für Numerik und Analysis von Populationsbilanzsystemen.


Aktuelle Aspekte der Forschung zu PBS am WIAS sind:

  • stochastische Simulationsalgorithmen für PBE,
  • höher-dimensionale Partikel,
  • Kopplung stochastischer Simulation und Finite Element Methode,
  • Partikel in turbulenten Strömungen,
  • Vergleich von QMoM und Finite Elemente Ansätzen bei PBE.



Besondere Beiträge des WIAS

Effiziente Numerik für Populationsbilanzsysteme

Populationsbilanzgleichungen numerisch zu lösen kann eine große Herausforderung darstellen. Abhängig von der Dimension des räumlichen Gebiets (äußere Koordinaten) und der Anzahl betrachteter Partikeleigenschaften (innere Koordinaten), können sehr hoch-dimensionale Gleichungen auftreten. Um diesem Effekt entgegenzutreten können unterschiedliche Methoden der Ordnungsreduktion eingesetzt werden. So kann zum Beispiel ein achsensymmetrisches Gebiet in 3D, wie es Strömungsrohre häufig sind, auf ein rechteckiges 2D-Gebiet transformiert werden. Auch können Operator-Splitting Schemata zum Einsatz kommen. Diese und ähnliche Methoden werden bezüglich ihrer Genauigkeit und Effizienz bewertet und zur Simulation realer Experimente genutzt.


Abbildung 2: Vergleich der Genauigkeit verschiedener numerischer Methoden bei der Approximation der Verteilung von Harnstoff-Kristallen am Auslass eines experimentellen Rohrkristallers.


Stochastik und Numerik kombiniert

Ein neuer Ansatz zur Kombination stochastischer Partikelsimulation und Numerik partieller Differentialgleichungen wird auf die Modellierung und Simulation eines ASS (Acetylsalicylsäure) Rohrkristallers angewandt. Eine gesättigte Lösung von ASS in Ethanol wird bei niedriger Geschwindigkeit durch einen dünnen Schlauch gepumpt. Die Strömung transportiert ungelöste ASS Kristalle, Kühlung von außen führt zu Übersättigung, welche wiederum Wachstum der Kristalle anregt. Zusätzlich führt die Verklumpung (Koagulation) von Partikeln zur Bildung größerer Aggregate. Rohrkristaller werden in der Industrie für die Produktion von Kristallen mit möglichst regelmäßiger Form und Größe eingesetzt, was unter anderem in der Erzeugung von Pharmazeutika relevant ist. In dem benutzten Simulationsalgorithmus wird der stochastische Partikel-Löser Brush an den Finite Elemente-Löser ParMooN gekoppelt. Ergebnisse in 2D (äußere Koordinaten) und einer inneren Koordinate (Partikelmasse) sind vielversprechend. Es besteht die Absicht, Probleme mit mehreren inneren Koordinaten in 3D auf die gleiche Weise anzugehen.




Abbildung 3: Momentaufnahme einer gekoppelten Simulation eines ASS Strömungskristallers. Links: Temperatur in K. Mitte: Molare Konzentration von ASS. Rechts: Massenkonzentration ungelöster ASS Kristalle, mit dem stochastischen Löser Brush berechnet.

Publikationen

  Artikel in Referierten Journalen

  • V. Wiedmeyer, F. Anker, C. Bartsch, A. Voigt, V. John, K. Sundmacher, Continuous crystallization in a helically-coiled flow tube: Analysis of flow field, residence time behavior and crystal growth, Industrial and Engineering Chemistry Research, 56 (2017), pp. 3699--3712, DOI 10.1021/acs.iecr.6b04279 .

  • R.I.A. Patterson, Properties of the solutions of delocalised coagulation and inception problems with outflow boundaries, Journal of Evolution Equations, 16 (2016), pp. 261--291.
    Abstract
    Well posedness is established for a family of equations modelling particle populations undergoing delocalised coagulation, advection, inflow and outflow in a externally specified velocity field. Very general particle types are allowed while the spatial domain is a bounded region of $d$-dimensional space for which every point lies on exactly one streamline associated with the velocity field. The problem is formulated as a semi-linear ODE in the Banach space of bounded measures on particle position and type space. A local Lipschitz property is established in total variation norm for the propagators (generalised semi-groups) associated with the problem and used to construct a Picard iteration that establishes local existence and global uniqueness for any initial condition. The unique weak solution is shown further to be a differentiable or at least bounded variation strong solution under smoothness assumptions on the parameters of the coagulation interaction. In the case of one spatial dimension strong differentiability is established even for coagulation parameters with a particular bounded variation structure in space. This one dimensional extension establishes the convergence of the simulation processes studied in [Patterson, Stoch. Anal. Appl. 31, 2013] to a unique and differentiable limit.

  • F. Anker, S. Ganesan, V. John, E. Schmeyer, A comparative study of a direct discretization and an operator-splitting solver for population balance systems, Comput. Chem. Engng., 75 (2015), pp. 95--104.
    Abstract
    A direct discretization approach and an operator-splitting scheme are applied for the numerical simulation of a population balance system which models the synthesis of urea with a uni-variate population. The problem is formulated in axisymmetric form and the setup is chosen such that a steady state is reached. Both solvers are assessed with respect to the accuracy of the results, where experimental data are used for comparison, and the efficiency of the simulations. Depending on the goal of simulations, to track the evolution of the process accurately or to reach the steady state fast, recommendations for the choice of the solver are given.

  • E. Schmeyer, R. Bordás, D. Thévenin, V. John, Numerical simulations and measurements of a droplet size distribution in a turbulent vortex street, Meteorologische Zeitschrift, 23 (2014), pp. 387--396.
    Abstract
    A turbulent vortex street in an air flow interacting with a disperse droplet population is investigated in a wind tunnel. Non-intrusive measurement techniques are used to obtain data for the air velocity and the droplet velocity. The process is modeled with a population balance system consisting of the incompressible Navier--Stokes equations and a population balance equation for the droplet size distribution. Numerical simulations are performed that rely on a variational multiscale method for turbulent flows, a direct discretization of the differential operator of the population balance equation, and a modern technique for the evaluation of the coalescence integrals. After having calibrated two unknown model parameters, a very good agreement of the experimental and numerical results can be observed.

    Eine turbulente Wirbelstraße in einer Luftströmung mit einer dispergierten Tröpfchenpopulation wird in einem Windkanal untersucht. Nichtintrusive Messtechniken werden verwendet, um Daten bezüglich der Luft-- und Tröpfchengeschwindigkeiten zu gewinnen. Der zu Grunde liegende Prozess wird mit einem Populationsbilanzsystem modelliert, welches aus den inkompressiblen Navier--Stokes--Gleichungen und einer Populationsbilanzgleichung für die Tröpfchenverteilungsdichte besteht. Numerische Simulationen werden durchgeführt, welche ein variationelle Mehrskalenmethode für turbulente Strömungen, eine direkte Diskretisierung des Differentialoperators der Populationsbilanzgleichung und ein modernes Verfahren zur Berechnung der Koaleszensintegrale verwenden. Nachdem zwei unbekannte Modellparameter kalibriert worden sind, kann eine sehr gute Übereinstimmung der experimentellen und numerischen Ergebnisse beobachtet werden.

  • V. John, C. Suciu, Direct discretizations of bi-variate population balance systems with finite difference schemes of different order, Chemical Engineering Sciences, 106 (2014), pp. 39--52.
    Abstract
    The accurate and efficient simulation of bi-variate population balance systems is nowadays a great challenge since the domain spanned by the external and internal coordinates is five-dimensional. This report considers direct discretizations of this equation in tensor-product domains. In this situation, finite difference methods can be applied. The studied model includes the transport of dissolved potassium dihydrogen phosphate (KDP) and of energy (temperature) in a laminar flow field as well as the nucleation and growth of KDP particles. Two discretizations of the coupled model will be considered which differ only in the discretization of the population balance equation: a first order monotone upwind scheme and a third order essentially non-oscillatory (ENO) scheme. The Dirac term on the right-hand side of this equation is discretized with a finite volume method. The numerical results show that much different results are obtained even in the class of direct discretizations.

  • R. Bordás, V. John, E. Schmeyer, D. Thévenin, Numerical methods for the simulation of a coalescence-driven droplet size distribution, Theoretical and Computational Fluid Dynamics. Springer-Verlag, Berlin., 27 (2013), pp. 253--271.
    Abstract
    A droplet size distribution in a turbulent flow field is considered and modeled by means of a population balance system. This paper studies different numerical methods for the 4D population balance equation and their impact on an output of interest, the time-space-averaged droplet size distribution at the outlet which is known from experiments. These methods include different interpolations of the experimental data at the inlet, various discretizations in time and space, and different schemes for computing the aggregation integrals. It will be shown that notable changes in the output of interest might occur. In addition, the efficiency of the studied methods is discussed.

  • L.G.M. DE Souza, G. Janiga, V. John, D. Thévenin, Reconstruction of a distribution from a finite number of moments with an adaptive spline-based algorithm, Chemical Engineering Sciences, 65 (2010), pp. 2741--2750.

  • V. John, M. Roland, On the impact of the scheme for solving the higher-dimensional equation in coupled population balance systems, International Journal for Numerical Methods in Engineering, 82 (2010), pp. 1450--1474.

  • V. John, T. Mitkova, M. Roland, K. Sundmacher, L. Tobiska, A. Voigt, Simulations of population balance systems with one internal coordinate using finite element methods, Chemical Engineering Sciences, 64 (2009), pp. 733--741.

  Beiträge zu Sammelwerken

  • V. John, M. Roland, Simulations of 3D/4D precipitation processes in a turbulent flow field, in: Numerical Mathematics and Advanced Applications 2009, G. Kreiss, P. Lötstedt, A. Målqvist, M. Neytcheva, eds., Springer, Heidelberg et al., 2010, pp. 479-- 487.
    Abstract
    Precipitation processes are modeled by population balance systems. A very expensive part of the simulation of population balance systems is the solution of the equation for the particle size distribution (PSD) since this equation is defined in a higher dimensional domain than the other equations in the system. This paper studies different approaches for the solution of this equation: two finite difference upwind schemes and a linear finite element flux--corrected transport method. It is shown that the different schemes lead to qualitatively different solutions for an output of interest.

  Preprints, Reports, Technical Reports

  • A. Jha, V. John, On basic iteration schemes for nonlinear AFC discretizations, Preprint no. 2533, WIAS, Berlin, 2018, DOI 10.20347/WIAS.PREPRINT.2533 .
    Abstract, PDF (814 kByte)
    Algebraic flux correction (AFC) finite element discretizations of steady-state convection-diffusion-reaction equations lead to a nonlinear problem. This paper presents first steps of a systematic study of solvers for these problems. Two basic fixed point iterations and a formal Newton method are considered. It turns out that the fixed point iterations behave often quite differently. Using a sparse direct solver for the linear problems, one of them exploits the fact that only one matrix factorization is needed to become very efficient in the case of convergence. For the behavior of the formal Newton method, a clear picture is not yet obtained.

  • C. Bartsch, V. John, R.I.A. Patterson, Simulations of an ASA flow crystallizer with a coupled stochastic-deterministic approach, Preprint no. 2483, WIAS, Berlin, 2018, DOI 10.20347/WIAS.PREPRINT.2483 .
    Abstract, PDF (378 kByte)
    A coupled solver for population balance systems is presented, where the flow, temperature, and concentration equations are solved with finite element methods, and the particle size distribution is simulated with a stochastic simulation algorithm, a so-called kinetic Monte-Carlo method. This novel approach is applied for the simulation of an axisymmetric model of a tubular flow crystallizer. The numerical results are compared with experimental data.

  Vorträge, Poster

  • R. Ahrens, F. Anker, C. Bartsch, A. Voigt, V. Wiedmeyer, K. Sundmacher, V. John, S. Le Borne, Advanced numerical methods for the simulation of population balance systems, 6th International Conference on Population Balance Modelling (PBM2018), Belgium, May 7 - 9, 2018.

  • C. Bartsch, V. John, R.I.A. Patterson, A new mixed stochastic-deterministic simulation approach to particle populations in fluid flows, 6th International Conference on Population Balance Modelling (PBM2018), Belgium, May 7 - 9, 2018.

  • V. John, A new mixed stochastic-deterministic simulation approach for particle populations in fluid flows, 6th European Seminar on Computing (ESCO), June 3 - 8, 2018, University of West Bohemia, Pilsen, Czech Republic, June 6, 2018.

  • R.I.A. Patterson, Confidence intervals for coagulation--advection simulations, Clausthal-Göttingen International Workshop on Simulation Science, April 27 - 28, 2017, Georg-August-Universität Göttingen, Institut für Informatik, April 28, 2017.

  • R.I.A. Patterson, Coagulation --- Transport Simulations with Stochastic Particles, CIM-WIAS Workshop ``Topics in Applied Analysis and Optimisation'', December 6 - 8, 2017, University of Lisbon, International Center for Mathematics, Lisboa, Portugal, December 7, 2017.

  • R.I.A. Patterson, Simulation of particle coagulation and advection, Numerical Methods and Applications of Population Balance Equations, October 13, 2017, GRK 1932, Technische Universität Kaiserslautern, Fachbereich Mathematik, October 13, 2017.

  • R.I.A. Patterson, Population balance simulation, University of Cambridge, Department for Chemical Engineering and Biotechnology, UK, May 5, 2016.

  • V. Wiedmeyer, F. Anker, A. Voigt, V. John, K. Sundmacher, Crystal shape evolution in a continuous helically coiled flow tube crystallizer (HCFT), 10th European Congress of Chemical Engineering (ECCE10), Nice, France, September 28 - 29, 2015.

  • F. Anker, A comparative study of the accuracy and efficiency of numerical techniques for the solution of population balance systems, 10th European Congress of Chemical Engineering (ECCE10), September 27 - October 1, 2015, Nice, France, September 28, 2015.

  • V. John, S. Le Borne, K. Sundmacher, Numerische Lösungsverfahren für gekoppelte Populationsbilanzsysteme zur dynamischen Simulation multivariater Feststoffprozesse am Beispiel der formselektiven Kristallisation, Evaluation Colloquium of the Priority Program (SPP) 1679 ``Dyn-Sim-FP -- Dynamische Simulation vernetzter Feststoffprozesse'', Frankfurt/Main, January 13, 2015.

  • V. John, On the numerical simulation of population balance systems, Karlsruher Institut für Technologie, Fakultät für Mathematik, December 9, 2009.