Refine
Document Type
- Doctoral thesis (17)
Has Fulltext
- yes (17)
Is part of the Bibliography
- no (17)
Year of publication
Keywords
- Braunkohle (5)
- Lignite (4)
- Trocknung (4)
- Wirbelschicht (4)
- Wirbelschichtverfahren (4)
- Fluidized bed (3)
- Heat transfer (3)
- Kraftwerk (3)
- Modellierung (3)
- Numerische Strömungssimulation (3)
Institute
PEM water electrolysis is a clean technology for hydrogen production. In spite of its many advantages, the costs of the conventional PEM electrolysis cell makes it commercially less competitive vis-à-vis its peers. An alternative cell design has been proposed which has up to a 25 % costs advantage over the conventional cell. In this alternative cell design, the flow channel plate which bears the most costs in the conventional cell design has been replaced with a 3-D Porous Transport Layer (PTL) structure. It has however, been observed that the conventional cell by far out performs the low cost cell at high current density operations, due to increased mass transport limitation in the later. Industrial and commercial hydrogen production efforts are focused towards high current density operation (> 3 A/cm²), so the alternative cell design must be optimized for mass transport limitation.
PEM water electrolysis is a clean technology for hydrogen production. In spite of its many advantages, the costs of the conventional PEM electrolysis cell makes it commercially less competitive vis-à-vis its peers. An alternative cell design has been proposed which has up to a 25 % costs advantage over the conventional cell. In this alternative cell design, the flow channel plate which bears the most costs in the conventional cell design has been replaced with a 3-D Porous Transport Layer (PTL) structure. It has however, been observed that the conventional cell by far out performs the low cost cell at high current density operations, due to increased mass transport limitation in the later. Industrial and commercial hydrogen production efforts are focused towards high current density operation (> 3 A/cm²), so the alternative cell design must be optimized for mass transport limitation.
This work seeks to understand the source of, and to eliminate the mass transport losses in the alternative cell design to get it performing at least as good as the conventional cell at current densities up to 5 A/cm². A 2-D non-isothermal semi-empirical fully-coupled models of both cell designs have been developed and experimentally validated. The developed validated models were then used as tools to simulate and predict the best operating conditions, design parameters and micro-structural properties of the PTL at which the mass transport issues in the alternate cell will be at its minimum, at high current densities. The models are based on a multi-physics approach in which thermodynamic, electrochemical, thermal and mass transport sub-models are coupled and solved numerically, to predict the cell polarization and individual overpotentials, as well as address heat and water management issues. The most unique aspect of this work however, is the development of own semi-empirical equations for predicting the mass transport overpotential imposed by the gas phase (bubbles) at high current densities. For the very first time, calculated polarization curves up to 5 A/cm² have been validated by own experimental data. The results show that, the temperature and pressure, water flowrate and thickness of the PTL are the critical parameters for mitigating mass transport limitation. It was found that, for the size of the cells studied (25 cm² active area each), when both cells are operating at the same temperature of 60 °C, alternative design will have a comparable performance to the conventional designed cell even at 5 A/cm² current density when; the operating pressure is ≥ 5 bar, the feed water flowrate is ≥ 0.024l/min∙cm², PTL porosity is 50 %, PTL pore size is ≥ 11 µm and PTL thickness is 0.5 mm. At these operating, design and micro-structural conditions, the predicted difference between the polarizations of both cells will be only ~10 mV at 5 A/cm² operating current density.
With the rapid growth of renewable energy sources (RES) in the power generation mix in accordance with the German energy transition policy (‘Energiewende’), fewer baseload coal power plants will be required. Future power generation will be supplied through decentralized power utilities such as off-shore wind parks and also through high operational flexibility of existing conventional coal power units. High operational flexibility means conventional power plants have to increase cyclic operations to cope with feed-ins from variable-RES such as wind and solar.
Unlike medium and peak load power plants that can react quickly to load changes and power ramps, baseload power plants are not suited for such operations. Important technical requirements for flexible operation include among others; frequent start-ups and shut-downs, a minimum downtime, shorter startup time and short operational periods. Baseload coal power plants however do not meet these requirements.
This increased cyclic mode of operation can have severe impacts on vital power plant components such as superheater and reheater tubes resulting in high temperature cyclic oxidation/corrosion especially because these plants were not designed for frequent cyclic operations. To optimize plant operations, minimize material damage and reduce operational and maintenance cost, it is therefore important to understand the oxidation and corrosion risk to plants materials associated with this flexible mode of operation.
In this context, thermochemical modeling in FactSage 6.4ᵀᴹ as well as experimental investigations were carried out. For the experimental investigations, five commercial coal boiler superheater and reheater materials, namely T91, VM12-SHC, TP347-HFG, DMV304 HCu and DMV310 N were exposed for 1000 hours under discontinuous isothermal oxidation conditions and 1000 hours thermo–cyclic oxidation conditions at a metal surface temperature of 650 °C. The synthetic corrosive flue gas consisted of a mixture of CO₂, O₂, SO₂, N₂ and H₂O. The test material samples were partly covered in fly ash to investigate the effect ash deposits on the corrosion and oxidation behavior of the test materials. After exposure metallographic analysis by means of light microscopy and scanning electron microscopy (LOM and SEM–EDS) were carried out to study the oxide morphology and micro–structural properties of the materials. The oxidation kinetics (weight change) results showed significant oxide growth rates (weight gain) under cyclic oxidation conditions especially in the martensitic alloys – T91, VM12-SHC.
Furthermore, metallographic analysis revealed severe oxide spallation in the ash covered sections of these alloys. The austenitic materials (TP374-HFG, DMV310 N) with the exception of DMV304 HCu showed good oxidation behavior with minimal oxide growth both under isothermal and thermal cyclic conditions. However, severe grain boundary attack and internal sulphidation were found in these alloys. DMV310 N showed the best corrosion and oxidation performance. The thermochemical modeling calculations supported the experimental results.
Im Rahmen der vorliegenden Arbeit wurde eine modellgestützte Prozessoptimierung an einer 400 kWth Verbrennungsanlage entwickelt, mit dem Ziel, die Betriebsoptimierung in Echtzeit zu realisieren. Die Grundlage bildet ein physikalisch-netzwerkgekoppeltes Modell der Verbrennungsanlage, das mit einem neuronalen Netzwerk zur Analyse der Betriebsdaten des Brenners hinsichtlich der Verbrennungsstabilität kombiniert wird. Der Optimierungsalgorithmus basiert auf einer Zielfunktion, die die Flammenstabilität und den Luftanteil berücksichtigt. Im Überblick zum aktuellen Stand der Technik wurde zunächst die Problemstellung der Flammenstabilität auf das untersuchte Feststoff-Brennersystem übertragen und der erforderliche Zugriff auf physikalisch-chemische Parameter der Flammenzone zur Zustandsbewertung des Systems erläutert. Das Modell wurde auf der Grundlage experimenteller Untersuchungen am Gas-/Kohle-Kombibrenner des Verbrennungsversuchsstands VVS400 an der BTU Cottbus-Senftenberg entwickelt. Dabei wurden verschiedene Messmethoden zur Erfassung instationärer Verbrennungseffekte überprüft. Insbesondere erwies sich die UV-Flammenstrahlung und die Chemolumineszenz des OH*-Radikals als geeignet für die Bewertung der primären Verbrennungszone. Diese Methodik ermöglichte es, das Verbrennungsverhalten zu analysieren und ein Modell zur Bewertung der Flammenstabilität zu entwickeln. Eine wesentliche Erkenntnis war, dass das Flammenwächtersignal, unterstützt durch eine Frequenzanalyse, eine zuverlässige Datenquelle für die Modellierung und Optimierung des Brennersystems darstellt. Das entwickelte Bewertungssystem ermöglicht eine präzise Beurteilung des Verbrennungszustands sowie die Berechnung von Optimierungsmaßnahmen, die als Handlungsempfehlungen für den Brennerbetrieb dienen. Im weiteren Verlauf wurde ein quasistationäres Modell entwickelt, das physikalische Modellierungsansätze mit neuronalen Netzwerken kombiniert, um komplexe Betriebsparameter abzubilden. Diese Modellierung zeigte eine hohe Übereinstimmung mit den experimentellen Messdaten und bildete somit eine effektive Grundlage für die Prozessoptimierung. Das Optimierungsverfahren basiert auf dem Pattern-Search-Verfahren, wobei die Flammenstabilität und der Luftanteil als Optimierungsziele definiert wurden. Die Ergebnisse der Untersuchungen belegen, dass das integrierte Verfahren zur Flammenbewertung und -optimierung erfolgreich zur Echtzeit-Analyse und Steuerung von Brennersystemen eingesetzt werden kann. Anhand von Validierungsmessungen wird das Potenzial zur Modellierung und Optimierung des Verbrennungszustands auf Basis der dynamischen Flammenstabilität am Beispiel eines Kohlebrenners diskutiert. Das vorgestellte Verfahren soll als Hilfsmittel für Brennersysteme dienen, um den steigenden Anforderungen an Emissionsgrenzen und Lastflexibilität gerecht zu werden. Die Übertragbarkeit vom experimentellen Rahmen auf die Industrieanwendung wird diskutiert, und weiterer Entwicklungsbedarf abgeleitet.
Im Rahmen der vorliegenden Arbeit wird der Transport von CO₂ aus energietechnischen Anlagen untersucht. Zunächst werden die physikalischen und technischen Randbedingungen erläutert, welche ein Transportsystem beeinflussen und beschränken, insbesondere werden die Einflüsse der thermo-dynamischen Fluideigenschaften von CO₂ und CO₂-Gemischen auf die Transportprozesse analysiert. Des Weiteren erörtert die Arbeit die Besonderheiten, die bei der Auslegung und dem Betrieb von CO₂-Pipelines (im Gegensatz beispielsweise zu Erdgaspipelines) auftreten.
Der Hauptteil der Arbeit befasst sich mit der Modellierung und Simulation zweier potentieller Pipeline-transportsysteme, welche CO₂ vom Kraftwerksstandort Jänschwalde in der Lausitz a) zu einem CO₂-Hub im Hafen von Rotterdam und b) zu einem Offshore-Aquifer unter der Nordsee transportieren. Anhand dieser konkreten Fallbeispiele werden die Auslegung des Transportsystems dargestellt und die statischen und die dynamischen Betriebskennlinien der Pipelines und des Injektionsbrunnens aufgenommen und hinsichtlich der Auswirkungen auf den Betrieb unter flexiblen Lastbedingungen untersucht.
Die Modellierung der Pipelines wird mithilfe der Modellierungs- und Simulationssoftware OLGA Dynamic Multiphase Flow Simulator von Schlumberger umgesetzt. Die Validierung der genutzten Submodelle der Software erfolgt über die Simulation historischer Lastfälle der Canyon-Reef-Carrier-Pipeline in Texas und der CO₂-Injektionspilotanlage Ketzin westlich von Potsdam.
Die Hauptmodelle umfassen a) die Pipeline vom Kraftwerkszaun bis zum Übergabepunkt am CO₂-Hub in Rotterdam bzw. b) die Pipeline vom Kraftwerkszaun bis zum Meeresgrund sowie den Injektionsbrunnen und ein Modell des Aquiferspeichers. Es ist mit den Modellen möglich, die thermo-fluiddynamischen Prozesse des CO₂ in der Rohrleitung sowie im Injektionsbrunnen bis zum Speicher nachzubilden. Es wurden insgesamt 66 Szenarien simuliert und ausgewertet.
Die simulierten Szenarien umfassen statische Kennlinien bei konstanter Last, definierte Lastzyklen mit Sprüngen konstanter Höhe, reale Lastgänge auf Basis historischer Kraftwerksmessdaten sowie typische Betriebsabläufe wie An- und Abfahren der Pipeline, Absperrung einzelner Abschnitte der Pipeline bis hin zu den Auswirkungen eines Risses in der Pipelinewand.
Zu den wesentlichen Ergebnissen der Simulationen zählen die auftretenden Druckverläufe, Wärme-übergangs- und Wärmespeichervorgänge sowie das dynamische Verhalten bei wechselnden Lastbedingungen.
Ausgehend von den Simulationsergebnissen werden die Pipelinesysteme charakterisiert und die dynamischen Prozesse hinsichtlich ihrer Auswirkungen auf Betriebsabläufe und Auslegungskriterien untersucht.
Es kann gezeigt werden, dass Pipelines für die gestellte Transportaufgabe geeignet sind und mit einem gewissen Maß an Prozessüberwachung und -regelung auch dauerhaft sicher betrieben werden können. Problematische Prozessverläufe werden eingehend diskutiert und Regelstrategien aufgezeigt, welche einen Betrieb der Pipelines sicherstellen können.
In der vorliegenden Arbeit werden die Modellierung und die numerische Simulation von Strömung, Wärme- und Stofftransport zur Abbildung eines Trocknungsprozesses von Braunkohlepartikeln in einer blasenbildenden Wirbelschicht behandelt. Hierbei wird das Euler-Euler Two Fluid Model zugrunde gelegt und die Beschreibung des Fließverhaltens der Feststoffphase mittels der Kinetic Theory of Granular Flow realisiert. Ausgehend von einem Überblick zum aktuellen Stand der Forschung wird der Bedarf nach einer geeigneten Konfiguration von Modellparametern und Submodellen hinsichtlich der Strömungsmodellierung sowie nach der Implementierung eines Trocknungsmodells abgeleitet und als Zielsetzung formuliert.
Experimentelle Voruntersuchungen an einer Wirbelschichtanlage im Labormaßstab schaffen die Datenbasis für die spätere Validierung des Modells. Dabei erfolgt zunächst eine Betrachtung einzelner Zustandspunkte von diskontinuierlichen Trocknungsprozessen, um die darin auftretenden, strömungstechnischen Veränderungen zu quantifizieren. Zur Bewertung werden die makroskopischen Eigenschaften sowie Charakteristiken meso- und mikroskaliger Strömungstrukturen herangezogen, die sich aus der Analyse von Druckfluktuationen ergeben. Des Weiteren werden die Trocknungsverläufe mehrerer Chargen für verschiedene Betriebsparameter aufgezeichnet.
In einem ersten Untersuchungsschwerpunkt werden ausschließlich strömungsmechanische Aspekte fokussiert. Auf Basis umfangreicher Sensitivitätsanalysen zu den Einflüssen der rheologisch relevanten Parameter, der Impulsaustauschfunktion einschließlich der Partikelsphärizität und der Randbedingungen für Behälterwand und Gaseinlass, sowie durch Vergleiche mit den experimentellen Daten wird eine Parameterkonfiguration des Strömungsmodells vorgeschlagen. Die damit erzielten Ergebnisse sind insgesamt zufriedenstellend und geben die beobachteten Veränderungen im Trocknungsprozess korrekt wieder. Ursachen für bestehende Abweichungen zum Experiment werden diskutiert.
Im zweiten Untersuchungsschwerpunkt wird der gesamte Trocknungsprozess betrachtet. Hierzu werden Wärme- und Stofftransportmechanismen innerhalb des Modells berücksichtigt, wobei die eigens implementierte Trocknungskinetik auf Ebene der Partikel ansetzt. Die Einflüsse der Sphärizität und zwei verschiedener Formulierungen zur Berechnung der Wärme- und Stoffübergangskoeffizienten werden untersucht. Für den favorisierten Parametersatz werden die Simulationen mit einem expliziten Vorwärtsverfahren gekoppelt, wodurch die Trocknungsverläufe auf makroskopischer Zeitskale approximiert und somit mit den experimentellen Daten verglichen werden können. Es wird gezeigt, dass das vorgeschlagene Gesamtmodell die Trocknungsverläufe bei niedrigen bis mäßigen Leerrohrgeschwindigkeiten zuverlässig abbilden kann. Auftretende Diskrepanzen werden diskutiert und weiterer Entwicklungsbedarf abgeleitet.
Im Rahmen dieser Arbeit wurden Dampferzeuger-Werkstoffe mit Chromgehalten von 0,1 bis 12 Ma.-% hinsichtlich ihres Korrosionsverhaltens unter konventionellen und Oxyfuel-Betriebsbedingungen untersucht. Dabei stand neben der Gasphasenkorrosion auch die Korrosion unter Beteiligung von Aschedepositionen im Fokus der Untersuchungen. Die Auslagerung der Werkstoffe erfolgte durch eine kombinierte Versuchsanordnung aus Technikumanlage und einem eigens errichteten Labor-Korrosionsversuchsstand, wobei eine Gesamtauslagerungsdauer von 1.110 Stunden realisiert wurde. Dabei konnte das in der Technikumanlage begonnene Schichtwachstum im Labor-Korrosionsversuchsstand erfolgreich weitergeführt werden. Mit Hilfe der kombinierten Versuchsanordnung, insbesondere des Labor-Korrosionsversuchsstandes, konnten reale Betriebsbedingungen erfolgreich simuliert werden. Nach erfolgter Auslagerung in der Feuerungsanlage zeigten die Proben beginnende Korrosionserscheinungen an den Werkstoffoberflächen. Die Untersuchungsergebnisse der einzelnen Werkstoffproben hinsichtlich des Korrosionsverhaltens ließen jedoch kaum Unterschiede zwischen höher und niedrig legierten Werkstoffen erkennen. Wie bei anderen Autoren konnte nur am Werkstoff mit 9-12 Ma.-% Chromgehalt (Werkstoff VM12SHC) Aufkohlung festgestellt werden. Dabei hat die Aufkohlungstiefe gegenüber den Proben nach 110 Stunden Auslagerungszeit unter allen Betriebsbedingungen zugenommen. Die werkstofftechnischen Analysen nach der Auslagerung im Labor-Korrosionsversuchsstand zeigte eine deutliche Abnahme großer Poren in den Oxidschichten. Die Anwesenheit von Aschedepositionen führt zur Ausbildung dickerer Oxidschichten. Ursache sind Reaktionen zwischen den gebildeten Oxidschichten und den angelagerten Aschepartikeln unter Ausbildung einer Reaktionsschicht. Die bei steigendem Chromgehalt in den Werkstoffen zu erwartende Abnahme der Korrosionsaktivität durch die Bildung einer schützenden Eisen-Chrom-Spinellschicht ist aufgrund der geringen Auslagerungstemperatur von 490°C nicht eingetreten. Mit Hilfe der in dieser Arbeit durchgeführten Untersuchungen konnte gezeigt werden, dass bei niedrigen Temperaturen (T = 490°C) auch für Werkstoffe mit Chromgehalten von 0,1-6 Ma.-% Chrom konstante parabolische Wachstumsraten auftreten. Es konnten keine Vorteile im Hinblick auf geringere Korrosionsraten analysiert werden. Inwiefern die Werkstoffe bei der Temperatur von 490°C auch ein anormales Oxidationsverhalten aufweisen, könnte im Rahmen einer separaten Arbeit erforscht werden. Sowohl die berechneten als auch die gemessenen Abzehrungsraten an den untersuchten Werkstoffen liegen über den zu erwartenden, normalbedingten Abzehrungsraten, sind jedoch als nicht kritisch einzustufen. Für die realisierten Auslagerungsbedingungen (T = 490°C, p = 1 bar) ist unter Oxyfuel-Bedingungen, sowohl mit als auch ohne Aschedepositionen, keine erhöhte Korrosion feststellbar. Die ermittelten parabolischen Wachstumskonstanten weisen die gleiche Größenordnung wie die in der Literatur verfügbaren parabolischen Wachstumsraten auf. Dementsprechend befanden sich die untersuchten Korrosionsprozesse sich noch nicht in einem Langzeitgleichgewicht. Die untersuchten Werkstoffe wurden im Hinblick auf ihre Korrosionsbeständigkeit anhand der Untersuchungsergebnisse eingestuft. Darüber hinaus erfolgte unter Berücksichtigung der Beschaffungskosten für die untersuchten Werkstoffe eine Erweiterung in Form eines kommerziellen Rankings. Die Werkstoffe 13CrMo4-5 und 16Mo3 weisen das beste Preis- / Leistungsverhältnis auf. Diese Werkstoffe können unter den realisierten Bedingungen ohne Einschränkungen hinsichtlich der Lebensdauer und der Kosten eingesetzt werden. Der Werkstoff 10CrMo9-10 weist teilweise höhere Korrosionsraten als die Werkstoffe 7CrMoVTiB10-10 und VM12SHC auf, ist jedoch deutlich kostengünstiger zu beschaffen. Deshalb erfolgt im kommerziellen Ranking die Einordnung des Werkstoffes 10CrMo9-10 vor den Werkstoffen 7CrMoVTiB10-10 und VM12SHC.
In der vorliegenden Arbeit wird die Kondensation von Wasserdampf im horizontalen Rohr bei Drücken bis 10 bar im Hinblick auf das Wärmeübertragungsverhalten und die entstehenden Druckdifferenzen untersucht. Um hierfür eine entsprechende Ausgangsbasis zu schaffen, werden im ersten Teil der Arbeit die bereits veröffentlichten experimentellen Untersuchungen zum Thema zusammengefasst und diskutiert. Anschließend werden die vorhandenen theoretischen Grundlagen der Hydrodynamik von Zweiphasenströmungen sowie zum Wärmeübertragungsverhalten dargelegt. Dabei wird jeweils auf die bestehenden Schwierigkeiten bei der physikalischen Beschreibung der Phänomene und damit einhergehenden Unsicherheiten bei der Berechnung von Zielgrößen aufmerksam gemacht. Zur Berechnung des Reibungsdruckverlustes werden die Modelle von Friedel und Lockhart/Martinelli und zur Berechnung des Wärmeübergangskoeffizienten die Modelle von Thome, Shah, Cavallini, Numrich/Müller, Dobson/Chato sowie Huhn gegenübergestellt.
Der eigentliche Schwerpunkt der Arbeit liegt daran anschließend bei der experimentellen Bestimmung der beiden Zielgrößen - spezifischer Druckverlust und Wärmeübergangskoeffizient - bei absoluten Sattdampfdrücken zwischen 4…10 bar und Eintrittsgeschwindigkeiten bis 21,5 m/s. Die Entwicklung und die Errichtung der hierzu erforderlichen Versuchsapparatur erfolgten im Rahmen dieser Arbeit. Mit acht geometrisch identischen Doppelrohrwärmeübertragern (Innendurchmesser Kondensationsrohr: 20,5 mm) lassen sich die benötigten lokalen Temperaturinformationen in Abhängigkeit des Kondensationsfortschrittes generieren, indem unter Konstanthalten der Ein- bzw. Austrittsbedingungen (Sattdampfdruck und -massenstrom, Kühlwassertemperatur und -durchsatz) die Anzahl durchströmter Wärmeübertrager variiert wird. Mit in Summe 63 Versuchspunkten bei drei Druckstufen und jeweils fünf Eintrittsgeschwindigkeiten werden 15 Wärmeübergangs- und Druckverlustverläufe als Funktion des Kondensatanteils bzw. der Rohrlängskoordinate ermittelt. Der Vergleich der experimentellen Befunde mit den vorgestellten und bekannten Theorien zeigt zum Teil starke Abweichungen und nicht wiedergegebene Abhängigkeiten. Die mittleren relativen Abweichungen zwischen experimentellen und von den bestehenden Modellvorstellungen berechneten Wärmeübergangskoeffizienten liegen zwischen 16…42 %.
Ausgehend von den experimentellen Befunden dieser Arbeit werden Berechnungsgleichungen für den lokalen Wärmeübergangskoeffizienten sowie den spezifischen Druckverlust entwickelt. Der darüber hinaus durchgeführte Abgleich mit experimentellen Befunden aus der Literatur führt dazu, dass ein breiterer Parameterbereich erfasst wird, als der in dieser Arbeit experimentell untersuchte. Die entwickelten Berechnungsmodelle ermöglichen eine deutlich bessere Vorhersage von Wärmeübergang und Druckverlust als die bestehenden Modellvorstellungen. Aufgrund der komplexen Phänomenologie des Kondensationsprozesses im horizontalen Rohr bleibt eine Extrapolation dennoch mit Risiken verbunden. Die empfohlene Vorgehensweise zur Weiterentwicklung des Wissensstandes auf diesem Gebiet wird im Ausblick der Arbeit gegeben.
In bubbling fluidized beds, bubble characteristics such as size, shape and velocity have a vital influence on the hydrodynamics of the bed and hence on its performance as a chemical reactor and/or a heat exchange unit. In many industrial applications such as lignite dryer heat transfer tubes are usually inserted to enhance the rate of heat and mass transfer and chemical conversion. However, their presence strongly influences the bubbling behavior of the beds. Therefore, reliable design and scale-up of these systems come only after fundamental understanding of the bubbling behavior is achieved. Therefore, in this research work the influences of tube bank geometries and particle size on bubble characteristics were thoroughly investigated. In this research work, both numerical and experimental studies were employed. For the experimental measurements, a new nonintrusive digital image analysis technique was developed. The technique allowed for the simultaneous measurements of bed expansion and various bubble properties. An in-house software was developed to fully automate the image acquisition and data processing procedure. For the numerical studies, the Eulerian-Eulerian two-fluid model based on the kinetic theory of granular flow was used. Though, this CFD model has been considered as a fundamental tool for modeling gas-solid fluidized beds, its quantitative validation remains insufficient for a wide range of reactor geometries and operating conditions. Therefore, in this work validation of the model using experimental measurements of bed expansion and bubble properties obtained from a pseudo-two-dimensional fluidized bed was performed. The influence of two-dimensional simulations and different modeling parameters such as the friction packing limit, drag model and solid-wall boundary conditions were investigated. The two-fluid model generally showed reasonable agreement with the experimental measurements of pressure drop, bed expansion and bubble properties in bubbling regime. However, as the gas superficial velocity is increased and the bed moved towards slugging and turbulent regimes a big deviation arose and the two-fluid model failed to predict reasonably the fluidized bed hydrodynamics for the freely bubbling bed. The mean bubble properties predicted by two-dimensional simulations were in reasonable agreement with experiments at lower superficial velocities. They deviated at higher bed height and this was more pronounced at higher gas superficial velocities. The results from three-dimensional simulations were in better agreement with the experimental measurements; however, the computational effort need was very high making them impractical for parametric studies and sensitivity analyses. It was also showed that the choice of friction packing limits, drag laws and specularity coefficients have little influence on the bubble properties. For a bubbling bed, both experimental measurements and numerical simulations showed that inserting horizontal tube banks had either no or marginal influence on the static bed pressure drop and bed expansion. On the other hand, bubble hydrodynamics were strongly influenced and controlled by the geometry of the immersed tubes. In freely bubbling fluidized beds bubble size as well as rise velocity increased with bed height and superficial velocity. In beds with immersed tubes, such general trends were completely disturbed. Tubes appeared to restrict rapidly growing bubbles. As a result the mean bubble diameter and rise velocity were lower in the vicinity of the tube banks than in the freely bubbling bed. Results from different particle sizes showed that in a freely bubbling bed increasing the mean particle size increased the bubble diameter and rise velocity as well as bed expansion. In fluidized beds with dense horizontal tubes on the other hand, the mean bubble properties were almost independent of the particle sizes.
In der vorliegenden Arbeit wird die numerische Modellierung der Hydrodynamik und Wärmeübertragung in blasenbildenden Wirbelschichten bestehend aus Braunkohlepartikeln untersucht. Zur Validierung der aus den Simulationen extrahierten Daten stehen zwei verschiedene Versuchsanlagen zur Verfügung. In einem quasi-2D Wirbelschichtkaltmodell können hydrodynamische Größen gemessen werden. Mithilfe transparenter Vorder- und Rückwände der Wirbelschicht kann eine automatische Bildanalyse per Digital Image Analysis Technique (DIAT) durchgeführt werden, um die Eigenschaften der aufsteigenden Blasen ermitteln zu können. Eine weitere zylindrische Wirbelschichtanlage mit einem einzelnen, horizontal eingetauchten Rohr dient zur Bestimmung des Wärmeübergangs von dem Rohr in die Wirbelschicht. Der Wärmeübergangskoeffizient wird dabei über eine Sonde gemessen, welche die Ermittlung der lokalen Verteilung der Werte über dem Gesamtumfang ermöglicht. Als Wirbelgut wurde in beiden Versuchsanlagen abgesiebte Trockenbraunkohle verwendet. Dem eigentlichen Schwerpunkt der Arbeit wird eine Untersuchung zu Simulationsparametern und -gittern vorangestellt, um deren Einfluss auf die Ergebnisse abschätzen und passende Einstellungen ableiten zu können. Dies beinhaltet die Diskretisierung der konvektiven Terme in den Erhaltungsgleichungen, die gewählten Konvergenzkriterien und Zeitschrittweiten, die mindestens notwendigen Simulationsdauern und Mittelungszeiträume sowie die Auflösung und die Dimension des Rechengitters. Zur Modellierung der Hydrodynamik der Braunkohlewirbelschichten wird das Euler-Euler Two Fluid Model (TFM) unter Einsatz der Kinetic Theory of Granular Flow (KTGF) zur Modellierung der Rheologie der Feststoffphase angewendet und die wichtigsten partikelrelevanten Modellparameter und Submodelle untersucht. Die in Braunkohlewirbelschichten sowohl ohne als auch mit dicht gepackten Tauchheizflächen durchgeführten Vergleiche werden anhand der experimentell ermittelten Versuchsergebnisse validiert. Die individuelle Wirkung der unterschiedlichen Parameterwerte und Submodelle ist im Allgemeinen in freien Wirbelschichten größer und deutlicher voneinander abgrenzbar als in den Wirbelschichten mit den Strömungswiderständen. Deren Einfluss ist aufgrund ihrer dichten Anordnung dominierend. Von den untersuchten Größen stellt die Impulsaustauschfunktion den bedeutendsten Faktor zur adäquaten Beschreibung der Hydrodynamik der Wirbelschicht dar. Für die Untersuchung der Wärmeübertragung in Braunkohlewirbelschichten werden verschiedene Ansätze zur Modellierung der Wärmeleitfähigkeiten von Gas- und Feststoffphase betrachtet. Neben einem Standardmodell für die Wirbelschicht gehören dazu ebenfalls zwei unterschiedliche Methoden zur Einbeziehung von Effekten, die in unmittelbarer Nähe eingetauchter Oberflächen zu beachten sind und die lokale Wärmeleitfähigkeit der Gas-Partikel-Suspension beeinflussen. Die Simulationen mit den genannten Modellen zur Beschreibung der Wärmeleitfähigkeiten werden mittels globaler und lokaler Mittelwerte der sich ergebenden Wärmeübergangskoeffizienten an einem eingetauchten, horizontalen Einzelrohr miteinander und den experimentell gewonnenen Daten verglichen. Zusätzlich erfolgt eine Variation der Leerrohrgeschwindigkeit sowie eine getrennte Diskussion der gas- und partikelkonvektiven Anteile. Die Resultate zeigen, dass ein Wärmeleitungsansatz in der Nähe einer Übertragerfläche zu deutlich niedrigeren Wärmeübergangskoeffizienten als die anderen beiden führt, da er dort für die Feststoffphase geringere Wärmeleitfähigkeiten annimmt. Weiterhin kann in den Simulationen unabhängig von Partikeldurchmesser und Geschwindigkeit ein starker Abfall der erzielten Wärmeübertragung im Bereich der Rohrflanken beobachtet werden, welcher die größte qualitative Abweichung zwischen den Daten aus den Simulationen und Experimenten darstellt. In einer Diskussion werden die Gründe für diese Diskrepanz beleuchtet.
Die vorliegende Arbeit behandelt den Wärmeübergang in Wirbelschichten mit horizontalen Wärmetauscherflächen in Form von Einzelrohren oder Rohrbündeln. Die zugrunde liegende Anwendung ist die (Druckaufgeladene) Dampfwirbelschicht-Trocknung von Braunkohle, welche auf eine Korngröße von 0-2 mm aufgemahlen ist ("Feinkorn", Sauterdurchmesser < 85 Mikrometer), was den Anforderungen von Staubbrennern im Kraftwerk genügt. In diesen Wirbelschichten dominiert deutlich der partikelkonvektive Wärmeübergangsmechanismus. Die verschiedenen (D)DWT- und andere Trocknungsverfahren werden verglichen und anhand von Energieeffizienz und großtechnischer Anwendbarkeit bewertet. Derzeit existieren in der Literatur große Unsicherheiten und Widersprüche bezüglich der Vorhersage des Wärmeübergangs bei unterschiedlichen Einzelrohr- und Rohrbündelgeometrien. Untersuchungsbereiche in der Arbeit sind die Wärmetauschergeometrie, der Einfluss von Korngrößenverteilung und beheizten Nachbarrohren im Rohrbündel sowie die Effekte der Anwesenheit und gleichzeitigen Verdampfung von Wasser. Die Vorhersagegüte der anerkannten Theorien von Martin und Molerus/Dietz wird anhand des (D)DWT-Verfahrens bezüglich des Wärmeübergangs vom Heizrohr in die Wirbelschicht überprüft. Dabei stellt sich die Verwendung des Martindurchmessers zur Beschreibung des stark polydispersen Wirbelgutes als geeignet heraus. Das zusätzlich eingeführte 40-my-Kriterium wurde als geeigneter Standard zum Vergleich von feinkörnig aufgemahlenen Braunkohlen unterschiedlicher Herkunft identifiziert. Die Auswirkungen des Rohrdurchmessers auf den Wärmeübergang wurden über den weiten Bereich von 12 - 200 mm untersucht und in die Theorie implementiert. Den Schwerpunkt der Arbeit bildet die Untersuchung unterschiedlicher Rohrbündelgeometrien, insbesondere mit versetzter Rohranordnung, zur Herleitung eines auf geometrischen dimensionslosen Kennzahlen basierenden Rohrbündelfaktors, welcher die Reduktion des Wärmeübergangskoeffizienten am Rohrbündel im Vergleich zum Einzelrohr vorhersagt. Die verwendeten Rohrbündel umfassten Rohrdurchmesser von 12 - 35 mm und Rohrteilungen von 1,5 bis 4,0. Die gefundenen Haupteinflussgrößen auf den Rohrbündelfaktor sind die horizontale und diagonale Teilung sowie der kleinste freie Rohrabstand zwischen den Rohren. Wichtige Randeffekte wurden ebenfalls betrachtet und diskutiert. Ergänzend wurden die Rohrbündelgeometrien anhand des spezifischen, volumenbezogenen Wärmeübergangskoeffizienten verglichen. In der Arbeit konnten relativ präzise und detaillierte Auslegungsgrundlagen für die Konstruktion von Rohrbündelwärmetauschern in Wirbelschichten erarbeitet werden. Diese umfassen den Anwendungsfall der (D)DWT, sind aber auch auf andere Feinkornwirbelschichten übertragbar.
