Erweiterte Diskrete-Elemente-Methode
Die Erweiterte Diskrete-Elemente-Methode, (englisch Extended Discrete Element Method) oder XDEM ist in der numerischen Mathematik ein Ansatz, der die Bewegung von granularen Materialien oder Partikel, wie sie durch die klassische Diskrete-Elemente-Methode (DEM) (Cundall[1] and Allen[2]) beschrieben wird, um weitere Eigenschaften wie den thermodynamischen Zustand, mechanische Spannung/Dehnung oder elektromagnetische Felder für jedes Partikel erweitert. Im Gegensatz zur Kontinuumsmechanik löst die XDEM die Partikelphase mit den verschiedenen zu einem Partikel gehörenden Prozessen auf. Während die Diskrete Element Methode die Position und Orientierung in Raum und Zeit der Partikel beschreibt, berechnet die Erweiterte Diskrete-Elemente-Methode zusätzlich Eigenschaften wie Temperatur und/oder Speziesverteilungen oder mechanische Wechselwirkung mit Strukturen.
Geschichtlicher Rückblick
Die Molekulardynamik, wie sie in den späten 50er Jahren von Alder et al.[3] und den frühen 60er von Rahman[4] entwickelt wurde, kann als ersten Schritt in Richtung der Extended Discrete Element Methode gesehen werden, obwohl die Kräfte aufgrund von Kontakten zwischen Partikel nur durch Potentiale wie beispielsweise Lennard-Jones Potentiale von Molekülen oder Atomen zur Beschreibung der Wechselwirkung ersetzt wurden. In ähnlicher Weise wurde die fluiddynamische Wechselwirkung zwischen Teilchen, die sich in einem Fluid bewegen, untersucht. Die strömungsmechanischen Widerstandskräfte aufgrund der Relativgeschwindigkeit zwischen Partikel und umgebenden Fluid wurden als zusätzliche Kräfte betrachtet, die auf die Partikel wirken. Mit diesem mehrphasigen Strömungsansatz bestehend aus einer festen und gasförmigen Phase wird die Partikelphase mit diskreten Methoden erfasst, während die gasförmige oder flüssige Phase mit kontinuumsmechanischen Methoden beschrieben wird. Dieser Ansatz wird deshalb als combined continuum and discrete model (CCDM) bezeichnet wie er von Kawaguchi et al.,[5] Hoomans,[6] Xu 1997[7] and Xu 1998.[8] verwendet wurde. Da die feste Phase bestehend aus Partikeln diskret beschrieben wird, entfallen konstitutive Gesetze, was zu einem besseren Verständnis der fundamentalen Physik führt. Dies wird von Zhu 2007 et al.[9] und Zhu 2008 et al.[10] in ihrer Literaturübersicht zu partikulären Strömungen modelliert mit CCDM bestätigt. CCDM wurde in den beiden letzten Jahrzehnten entwickelt und beschreibt die granulare Phase mit der Diskrete-Elemente-Methode (DEM) während die gasförmige oder flüssige Phase mit den Navier-Stokes Gleichungen behandelt wird. Damit hat sie sich als effektive Methode etabliert, um Wechselwirkungen zwischen einer Partikel- und einer Fluidphase zu untersuchen wie in der Literaturübersicht von Yu und Xu,[11] Feng und Yu[12] und Deen et al.[13] dargestellt wird. Basierend auf der CCDM Methode wurden die charakteristischen Eigenschaften von spouted oder fluidisierten Betten von Gryczka et al.[14] erfasst.
Die theoretischen Grundlagen zu XDEM wurden von Peters[15] gelegt, der die Verbrennung von Holzhackschnitzel auf einem Vorschubrost berechnete.[16] Dieses Konzept wurde später von Simsek et al.[17] übernommen, um eine Rostfeuerung zu beschreiben. Anwendungen für den Hochofen wurden durch Shungo et al.[18] durchgeführt. Heute ist die numerische Simulation von Einspritzvorgängen von Flüssigkeiten wie Diesel in ein gasförmiges Medium ein Standardwerkzeug für eine große Zahl von CFD-Kodes wie Star-CD von CD-adapco, Ansys und AVL-Fire. Tropfen, die bei der Zerstäubung entstehen, werden mit einem null-dimensionalen Ansatz beschrieben, um Stoff- und Wärmeübergang in die Gasphase zu bestimmen.
Methodologie
Zahlreiche Anforderungen in den Ingenieurswissenschaften existieren und bilden sich heraus, die gleichzeitig eine kontinuierliche und diskrete Phase enthalten und deshalb nicht mit entweder kontinuierlichen oder diskreten Methoden alleine in hinreichender Genauigkeit gelöst werden können. Aus diesem Grunde wurde XDEM entwickelt, dass eine ausgezeichnete Plattform darstellt, mit der diskrete und kontinuierliche Phasen für eine sehr große Anzahl von Anwendungen gekoppelt werden. Obwohl Forschung und Entwicklung der numerischen Methoden auf den Gebieten der diskreten und kontinuierlichen Lösungsmethoden anhält, haben diese Werkzeuge mittlerweile einen hohen Entwicklungsstand erreicht. Um diskrete und kontinuierliche Löser zu koppeln existieren zwei Konzepte:
- Monolithisches Konzept: Die die mehrphasige Anwendung beschreibenden Gleichungen werden simultan behandelt und ergeben eine umfassende Lösung des Problems.
- Sequentielles Konzept: Die Gleichungen, die die diskrete und kontinuierliche Phase beschreiben, werden separat und sequentiell durch speziell entwickelte Methoden gelöst, wobei die Ergebnisse eines Gleichungssatzes für die Verbleibenden zur Verfügung gestellt wird.
Das erste Konzept erfordert einen Löser, der die gesamten physikalischen Eigenschaften des Systems erfasst und deshalb große Anforderungen an die Implementierung stellt. Allerdings existieren Szenarien, für die die Matrixkoeffizienten von gekoppelten Differentialgleichungen nicht oder nur unter sehr hohem Aufwand bestimmt werden können. Aus diesem Grunde zeigt das sequentielle Konzept als eine Kopplung zwischen verschiedenen Lösern, die die diskreten und kontinuierlichen Phasen repräsentieren, große Vorteile gegenüber dem monolithischen Ansatz. Er bietet einen großen Grad an Freiheit für die Kopplung von einer Vielzahl von Lösern. Außerdem bleibt eine mehr modulare Software erhalten, die individuell adaptierte Lösungsmethoden zulässt. Allerdings erfordert das sequentielle Konzept stabile und ausreichend genaue Kopplungen zwischen den einzelnen Modulen. Innerhalb des sequentiellen Konzepts für XDEM werden die Ergebnisse für kontinuierliche Probleme durch Lösen der relevanten Erhaltungsgleichungen berechnet. Eigenschaften wie Temperatur von individuellen Partikeln werden ebenfalls durch Lösen der dazugehörigen Gleichungen für das Partikel gefunden, das die interne Verteilung der relevanten Variablen in Raum und Zeit repräsentiert. Wesentliche Erhaltungsprinzipien mit den assoziierten Variablen, die für die Partikel gelten, sind in der folgenden Tabelle zusammengestellt:
Erhaltungsgesetz | Gleichung | Variable |
---|---|---|
Masse (kompressibles Medium) | Kontinuität | Druck/Dichte |
Linearer Impuls | Navier-Stokes | Geschwindigkeit |
Energie | Energie | Temperatur |
Spezies | Spezies Transport | Massenanteile |
Ladung, Strom | Maxwell | elektrisches/magnetisches Feld, elektrisches Verschiebungsfeld |
Die Lösung dieser Gleichungen definiert das dreidimensionale und instationäre Feld für die relevanten Variablen wie Temperatur oder Spezies. Allerdings wird die Auflösung im Speziellen in räumliche Dimensionen dadurch begrenzt, dass die Gleichungen für eine große Zahl von Partikeln gelöst werden muss. Folglich reduziert sich die räumlich Auflösung in der Regel meist auf eine repräsentative Dimension, was ausreichende Genauigkeit ergibt und durch Man and Byeong,[19] belegt wird, während der instationäre Charakter von Lee et al.[20] hervorgehoben wird.
Anwendungen
Anwendungen, die eine kontinuierliche als auch diskrete Phase beinhalten, treten in der pharmazeutischen Industrie wie beispielsweise Tablettenherstellung, Nahrungsmittelindustrie, Bergbau, Bau- und Landmaschinen, Metallgewinnung, Energiegewinnung und Systembiologie auf. Einige hervorstechende Beispiele sind Kaffeebohnen, Corn Flakes, Nüsse, Kohle, Sand, regenerative Kraftstoffe wie Biomasse und Düngemittel. Ursprünglich waren solche Untersuchungen laut Hoomans[21] auf einfache Strömungskonfigurationen begrenzt. Allerdings konnten Chu and Yu[22] demonstrieren, dass die Methode auf komplexe Strömungskonfigurationen wie eine Wirbelschicht, Förderband oder Zyklon angewendet werden kann. In ähnlicher Weise hat Zhou et al.[23] den CCDM Ansatz auf die komplexe Geometrie einer Kohlepulververbrennung angewendet und Chu et al.[24] modellierte die Strömung von Luft, Kohlestaub und Magnetitpartikeln verschiedener Größe.
Der CCDM Ansatz wurde ebenfalls für Wirbelschichten eingesetzt wie von Rowe und Nienow[25] und Feng and Yu[26] recherchiert und von Feng and Yu[27] zur Untersuchung der chaotischen Bewegung von Partikeln in Wirbelschichten untersucht wurde. Kafuia et al.[28] beschreiben ebenfalls Partikel-Kontinuum Fluid Modellierung für Wirbelschichten. Weitere Anwendungen von XDEM beinhalten thermische Konvertierung von Biomasse auf Vor- und Rückschubrosten. Wärmeübergang in einem Festbettreaktor wurde ebenfalls für durch ein Festbett strömende heiße Luft untersucht, wobei die Partikel durch Verteilung und Position unterschiedlichen Wärmeübergang erfahren. Die Deformation eines Förderbandes aufgrund von auftreffenden granularen Materialien stellt eine Anwendung aus der Spannungs- und Dehnungsanalyse dar.
- Verteilung der Oberflächentemperatur von Partikeln auf einem Rückschubrost
- Reaktionsfortschritt der Pyrolyse von Stroh auf einem Vorschubrost
- Verteilung der Porosität und Partikeltemperatur in einem Festbettreaktor
Literaturübersicht
- P. A. Cundall, O. D. L. Strack: A discrete numerical model for granular assemblies. In: Geotechnique. 29, 1979, S. 47–65.
- M. P. Allen, D. J. Tildesley: Computer Simulation of Liquids. Clarendon Press Oxford, 1990.
- B. J. Alder, T. E. Wainwright: Studies in Molecular Dynamics. I. General Method. In: J. Chem. Phys.. 31, 1959, S. 459.
- A. Rahman: Correlations in the Motion of Atoms in Liquid Argon. In: Phys. Rev.. 136, 1964.
- T. Kawaguchi, Y. Tsuji, T. Tanaka: Discrete particle simulation of two-dimensional fluidized bed. In: Powder Technol.. 77, 1993.
- B. P. B. Hoomans, J. A. M. Kuipers, W. J. Briels, W. P. M. Van Swaaij: Discrete particle simulation of bubble and slug formation in a two-dimensional gas-fluidized bed: A hard-sphere approach. In: Chem. Eng. Sci.. 51, 1996.
- B. H. Xu, A. B. Yu: Numerical simulation of the gas-solid flow in a fluidized bed by combining discrete particle method with computational fluid dynamics. In: Chemical Engineering Science. 52, 1997, S. 2785.
- B. H. Xu, A. B. Yu: Comments on the paper numerical simulation of the gas-solid flow in a fluidized bed by combining discrete particle method with computational fluid dynamics. In: Chemical Engineering Science. 53, 1998, S. 2646–2647.
- H. P. Zhu, Z. Y. Zhou, R. Y. Yang, A. B. Yu: Discrete particle simulation of particulate systems: Theoretical developments. In: Chemical Engineering Science. 62, 2007, S. 3378–3396.
- H. P. Zhu, Z. Y. Zhou, R. Y. Yang, A. B. Yu: Discrete particle simulation of particulate systems: A review of major applications and findings. In: Chemical Engineering Science. 63, 2008, S. 5728–5770.
- B. H. Xu, A. B. Yu: Particle-scale modelling of gas-solid flow in fluidisation. In: Journal of Chemical Technology and Biotechnology. 78, 2003, S. 111–121.
- Y. Q. Feng, A. B. Yu, A. B. Yu, A. Vince: Assessment of model formulations in the discrete particle simulation of gas-solid flow. In: Industrial & Engineering Chemistry Research. 43, 2004, S. 8378–8390.
- N. G. Deen, M. V. S. Annaland, M. A. Van Der Hoef, J. A. M. Kuipers: Review of discrete particle modeling of fluidized beds. In: Chemical Engineering Science. 62, 2007, S. 28–44.
- O. Gryczka, S. Heinrich, N. S. Deen, M. van Sint Annaland, J. A. M. Kuipers, M. Mörl: CFD modeling of a prismatic spouted bed with two adjustable gas inlets. In: Canadian Journal of Chemical Engineering. 87, 2009, S. 318–328.
- B. Peters: Classification of combustion regimes in a packed bed based on the relevant time and length scales. In: Combustion and Flame. 116, 1999, S. 297–301.
- B. Peters: Measurements and application of a discrete particle model (DPM) to simulate combustion of a packed bed of individual fuel particles. In: Combustion and Flame. 131, 2002, S. 132–146.
- E. Simsek, B. Brosch, S. Wirtz, V. Scherer, F. Krüll: Numerical simulation of grate firing systems using a coupled CFD/Discrete Element Method (DEM). In: Powder Technology. 193, 2009, S. 266–273.
- Shungo Natsui, Shigeru Ueda, Zhengyun Fan, Nils Andersson, Junya Kano, Ryo Inoue, Tatsuro Ariyama: Characteristics of solid flow and stress distribution including asymmetric phenomena in blast furnace analyzed by discrete element method. In: ISIJ International. 50, 2010, S. 207–214.
- Y. H. Man, R. C. Byeong: A numerical study on the combustion of a single carbon particle entrained in a steady flow. In: Combustion and Flame. 97, 1994, S. 1–16.
- J. C. Lee, R. A. Yetter, F. L. Dryer: Numerical simulation of laser ignition of an isolated carbon particle in quiescent environment. In: Combustion and Flame. 105, 1996, S. 591–599.
- B. P. B. Hoomans, J. A. M. Kuipers, W. J. Briels, W. P. M. Van Swaaij: Discrete particle simulation of bubble and slug formation in a two-dimensional gas-fluidized bed: A hard-sphere approach. In: Chem. Eng. Sci.. 51, 1996.
- K. W. Chu, A. B. Yu: Numerical simulation of complex particle-fluid flows. In: Powder Technology. 179, 2008, S. 104–114.
- H. Zhou, G. Mo, J. Zhao, K. Cen: DEM-CFD simulation of the particle dispersion in a gas-solid two-phase flow for a fuel-rich/lean burner. In: Fuel. 90, 2011, S. 1584–1590.
- K. W. Chu, B. Wang, A. B. Yu, A. Vince, G. D. Barnett, P. J. Barnett: CFD-DEM study of the effect of particle density distribution on the multiphase flow and performance of dense medium cyclone. In: Minerals Engineering. 22, 2009, S. 893–909.
- P. N. Rowe, A. W. Nienow: Particle mixing and segregation in gas fluidized beds: A review. In: Powder Technology. 15, 1976, S. 141–147.
- Y. Q. Feng, A. B. Yu, A. B. Yu, A. Vince: Assessment of model formulations in the discrete particle simulation of gas-solid flow. In: Industrial & Engineering Chemistry Research. 43, 2004, S. 8378–8390.
- Y. Q. Feng, A. B. Yu: An analysis of the chaotic motion of particles of different sizes in a gas fluidized bed. In: Particuology. 6, 2008, S. 549–556.
- K. D. Kafuia, C. Thornton, M. J. Adams: Discrete particle-continuum fluid modelling of gas-solid fluidised beds. In: Chemical Engineering Science. 57, 2002, S. 2395–2410.