Verfahren zur Erstellung eines Hydrauliknetzwerkes für einen optimierten Wärmeübertragungs- und Stofftransport
Technisches Anwendungsgebiet
Die vorliegende Erfindung betrifft ein Verfahren zur Erstellung eines Hydrauliknetzwerkes zur Wärme- und/oder StoffÜbertragung, bei dem eine vorgegebene Übertragungsfläche oder ein vorgegebenes Übertragungs- volumen mit einem mehrfach verzweigten Hydrauliknetzwerk versehen wird.
Der Wärmetransport mit Fluiden stellt eine wichtige technische Aufgabe sowohl bei konventionellen Wärmetauschern als auch bei der Nutzung erneuerbarer Energien dar, so beispielsweise bei der Übertragung der absorbierten Solarstrahlung auf das Wärmeträgermedium in einem Sonnenkollektor als auch bei der Abgabe dieser Energie an einen Speicher. Bei der energetischen
Betrachtung solcher Systeme muss sowohl die thermische Effizienz bei der Wärmeübertragung als auch die für den Fluidtransport benötigte hydraulische Leistung berücksichtigt werden. Um Wärme von einer Wärmeübertragungs- fläche auf eine Kanalstruktur zu übertragen, müssen die Kanäle möglichst gleichmäßig auf diese Fläche verteilt werden.
Aus dem Stand der Technik sind bereits zahlreiche Hydrauliknetzwerke bekannt, bei denen die Kanäle relativ gleichmäßig auf die Wärmeübertragungsfläche verteilt sind. So sind beispielsweise Wärmeübertragungssysteme bekannt, bei denen die Fluidkanäle
mäanderförmig über die Fläche verlaufen. Eine derartige Ausgestaltung wird auch als Serpentinenabsorber bezeichnet. Weiterhin sind beispielsweise auf dem Gebiet der Solartechnik so genannte Harfenabsorber bekannt, bei denen die Kanäle parallel zueinander angeordnet sind und auf der Ein- und Austrittsseite in Sammelkanäle münden. Eine weitere bekannte Ausgestaltung ist, der so genannte Kissenabsorber, bei dem das Fluid eine die Wärmeübertragungsfläche begrenzende flache Kammer durchströmt, so dass es vollflächig mit der Wärmeübertragungsfläche in Kontakt ist .
Diese bisher bekannten Konstruktionen sind jedoch häufig mit Nachteilen verbunden. So können bei mäander- förmig angeordneten Fluidkanälen hohe Druckverluste entstehen, die eine hohe Pumpleistung erfordern. Eine Parallelanordnung von Kanälen kann aufgrund der Druckverluste beim Übergang des Fluids von den Sammelkanälen in die parallelen Kanäle bzw. umgekehrt zu einer ungleichmäßigen Volumenstromverteilung und damit zu entsprechenden Temperaturunterschieden auf der Wärmeübertragungsfläche führen. Bei vollflächigen Konstruktionen besteht schließlich die Gefahr, dass einzelne Bereiche der Kammer schlecht oder gar nicht durchströmt werden. Bei den letztgenannten Systemen ist zudem die Wärmekapazität aufgrund des größeren Fluid- inhalts in der Regel höher als bei kanalartigen Systemen, was sich negativ auf das dynamische Verhalten auswirken kann.
Aus der EP 0 857 554 A2 sind Formen und Werkzeuge mit einem feingliedrig verzweigten Temperierkanal- netzwerk bekannt, das in der Art eines arteriellen
Gefäßsystems aufgebaut ist. In einer in der Druckschrift dargestellten Ausgestaltung wird hierfür ein binär verzweigtes Hydrauliknetzwerk eingesetzt. Ein grundsätzliches Problem besteht jedoch in der Erstellung eines derartigen verzweigten Hydrauliknetzwerkes, so dass es die Wärmeübertragungsfläche bzw. das Wärmeübertragungsvolumen möglichst gleichmäßig durchsetzt. Durch die in der Druckschrift angegebene Vorgehensweise, in der Natur vorkommende Verzweigungs- arten als Vorbild zu nehmen, wird die Erstellung eines derartigen Netzwerkes deutlich erschwert.
Ausgehend von diesem Stand der Technik besteht die Aufgabe der vorliegenden Erfindung darin, ein Verfahren zur Erstellung eines Hydrauliknetzwerkes zur Wärme- und/oder Stoffübertragung anzugeben, das sich bei beliebigen Geometrien der Übertragungsfläche bzw. des Ubertragungsvolumens einfach durchführen lässt und einen optimierten Wärmeübertragungs- und Stofftransport ermöglicht.
Darstellung der Erfindung
Die Aufgabe wird mit dem Verfahren gemäß Patentanspruch 1 gelöst. Vorteilhafte Ausgestaltungen des Verfahrens sind Gegenstand der Unteransprüche oder lassen sich aus der nachfolgenden Beschreibung sowie dem Ausführungsbeispiel entnehmen.
Bei dem vorliegenden Verfahren zur Erstellung eines Hydrauliknetzwerkes zur Wärme- und/oder Stoff-
Übertragung, bei dem eine vorgegebene Übertragungs- fläche oder ein vorgegebenes Ubertragungsvolumen mit einem mehrfach, vorzugsweise binär, verzweigten
Hydrauliknetzwerk versehen wird, wird der Verlauf des Hydrauliknetzwerkes durch folgende Schritte festgelegt:
- Festlegen eines Fluideintrittspunktes und eines Fluidaustrittspunktes; - Legen eines ersten Pfades zwischen Fluideintritts- und Fluidaustrittspunkt mit annähernd gleichem Abstand zu seitlichen Begrenzungen der Übertragungsfläche oder des Ubertragungsvolumens;
- Festlegen von ersten Verzweigungsstartpunkten und - endpunkten an Stellen des ersten Pfades, an denen ein vorgebbarer Höchstabstand zu den seitlichen Begrenzungen überschritten bzw. ein vorgebbarer Mindest- abstand unterschritten wird, in einem ersten Iterationsschritt ; - Verbinden der ersten Verzweigungsstartpunkte und - endpunkte durch Legen jeweils zumindest zweier zweiter Pfade mit annähernd gleichem Abstand zu nachstliegenden seitlichen Begrenzungen der Übertragungsfläche oder des Ubertragungsvolumens und dem ersten Pfad; - Festlegen von weiteren Verzweigungsstartpunkten und - endpunkten an Stellen der jeweils unmittelbar vorangehend gelegten Pfade, an denen ein vorgebbarer Höchst- abstand zu nachstliegenden seitlichen Begrenzungen der Übertragungsfläche oder des Ubertragungsvolumens überschritten bzw. ein vorgebbarer Mindestabstand unterschritten wird, in einem weiteren Iterationsschritt;
- Verbinden der weiteren Verzweigungsstartpunkte und - endpunkte durch Legen jeweils zumindest zweier weiterer Pfade mit annähernd gleichem Abstand zu nachstliegenden seitlichen Begrenzungen der Übertragungsfläche oder des Ubertragungsvolumens und einem nachstliegenden der bisher gelegten Pfade; und
- gegebenenfalls Wiederholung des weiteren Iterationsschrittes, bis eine vorgebbare Anzahl von Iterations- schritten oder ein vorgebbarer Mindestabstand der als Letztes gelegten Pfade zu den nachstliegenden seit- liehen Begrenzungen der Übertragungsfläche oder des Ubertragungsvolumens erreicht ist . Die nach Beendigung der Iterationsschritte im letzten Iterationsschritt erhaltenen Pfade sowie deren Verbindung mit dem Fluideintritts- und Fluidaustrittspunkt über Teile der bisher gelegten Pfade legen die Struktur des Hydrauliknetzwerkes fest, das anschließend entsprechend gefertigt wird.
Vorzugsweise wird als Hydrauliknetzwerk bei einer vorgegebenen Übertragungsfläche ein binäres Netzwerk und bei einem vorgegebenen Ubertragungsvolumen ein ternäres Netzwerk gebildet, so dass an jeder Verzweigung ein Kanal in entweder zwei oder drei Pfade bzw. Kanäle verzweigt. Mit einer idealen Verzweigungs- Struktur teilt sich der Volumenstrom des Fluids an jedem Verzweigungspunkt jeweils zur Hälfte auf den rechten und linken Kanal auf. Dies geschieht in allen Iterationsstufen, also bei jeder weiteren Verzweigung, so dass eine gleichmäßige Volumenstromverteilung entsteht. Gleichzeitig sind die Verzweigungswinkel nur verhältnismäßig klein, so dass die Strömung nur sanft umgelenkt wird, die einzelnen Kanäle annähernd parallel zueinander verlaufen und der Druckverlust gering ist. Bei manchen konkreten Anwendungen gibt es aufgrund einer asymmetrischen Übertragungsfläche zwangsläufig leichte Abweichungen von diesem vorzugsweisen Idealfall.
Vorzugsweise werden die Durchmesser bzw. die Querschnitte der mit jeder Verzweigung neu gebildeten Kanäle der Hydraulikstruktur jeweils so gewählt, dass sich der Volumenstrom an jeder Verzweigung zu gleichen Teilen auf die neu gebildeten Kanäle aufteilt. Die Durchmesser der einzelnen Kanäle können dabei unterschiedlich ausfallen, falls unterschiedliche Längen der neu gebildeten Kanäle kompensiert werden müssen.
Selbstverständlich lässt sich jedoch auch eine andere Verteilung der Kanalquerschnitte an den Verzweigungen wählen, falls dies für eine Optimierung der Wärmeübertragung sowie eine Minimierung Druckverluste bei einem vorgegebenen Ubertragungsvolumen oder einer vorgegebenen Übertragungsfläche erforderlich ist.
Vorzugsweise werden die Vernetzungsparameter des Hydrauliknetzwerkes bei der Durchführung des vor- liegenden Verfahrens mit einem Optimierungsverfahren variiert, um die für die jeweilige Anwendung optimalen Vernetzungsparameter zu erhalten. Als Vernetzungs- parameter kommen hierbei beispielsweise die Mindest- und Höchstabstände, die Positionen von Fluideintritts- und Fluidaustrittspunkt , die Zahl der Iterations- schritte und/oder das Verhältnis der Durchmesser der Kanäle der Hydraulikstruktur vor und nach jeder Verzweigung in Frage. Hierbei können sowohl einzelne als auch beliebige Kombinationen dieser Vernetzungs- parameter dem Optimierungsverfahren unterzogen werden. Bei dem Optimierungsverfahren werden die Vernetzungsparameter ausgehend von einem vorgegebenen Start- parametersatz variiert und damit verschiedene
Geometrien der Hydraulikstruktur erzeugt . Diese unterschiedlichen Geometrien werden durch thermische und hydraulische Simulationen hinsichtlich ihrer Gesamteffizienz berechnet. Die Gesamteffizienz ergibt sich dabei aus einer Berechnung der thermischen Effizienz und der hydraulischen Leistung. Beide Berechnungen können mit einem so genannten CFD-Programm (Computational Fluid Dynamics) oder mit anderen vereinfachten Simulationswerkzeugen durchgeführt werden. Die Parameter der besten Variante oder eine
Kombination von Parametern aus mehreren der Varianten wird als neuer Startparametersatz verwendet. Von diesen werden wiederum neue Varianten erzeugt usw.. Durch dieses Optimierungsverfahren, beispielsweise das Optimierungsverfahren der „Evolutionsstrategie", können insbesondere vielparametrische Probleme, für die es keine analytische Gesamtbeschreibung gibt, gelöst werden.
Das Verfahren zur Erstellung der Hydraulikstrukturen lässt sich beispielsweise für Anwendungen in Sonnenkollektoren, Heiz- und Kühlelementen jeglicher Art und Größe, Latentwärmespeichern oder Brennstoffzellen, um nur einige zu nennen, einsetzen. Insbesondere dreidimensionale Hydrauliknetzwerke können beispielsweise in Sorptionswärmespeichern eingesetzt werden, bei denen sowohl der Wärme- als auch der Stofftransport möglichst gleichmäßig im Raum erfolgen soll. Dabei soll jedoch das Speichervolumen nur minimal eingeschränkt werden, so dass Strukturen benötigt werden, die bei geringem Eigenvolumen möglichst den gesamten Raum bedienen können. Ebenso kann das vorliegende Verfahren zum Erzeugen von Strömungskanälen
für den Stofftransport (Reaktionsgase und entstehendes Wasser) in Brennstoffzellen eingesetzt werden.
Ein besonderer Vorteil des vorliegenden Verfahrens besteht darin, dass es sich nach Vorgabe des Fluid- eintritts- und Fluidaustrittspunktes sowie der entsprechenden Höchst- bzw. Mindestabstände und/oder Iterationen ebenso wie beim zusätzlichen Einsatz des Optimierungsverfahrens vollkommen automatisiert durch- führen lässt.
Zur technischen Umsetzung sind verschiedene Herstellverfahren nutzbar, so beispielsweise Prägen, Gießen, Rollbond, Fräsen oder Laminieren. Auch bekannte Rapid-Prototyping-Verfahren eignen sich selbstverständlich zum Herstellen des Hydrauliknetzwerkes. Ebenso können diverse Materialien zum Einsatz kommen, wie beispielsweise Kunststoff, Metall, Epoxydharz, Gießkeramik und andere .
Kurze Beschreibung der Zeichnungen
Das vorliegende Verfahren wird nachfolgend anhand eines Ausführungsbeispiels in Verbindung mit den Zeichnungen nochmals kurz erläutert. Hierbei zeigen:
Fig. 1 eine Gegenüberstellung von zwei bekannten Hydraulikstrukturen des Standes der Technik mit einem nach dem vorliegenden Verfahren erstellten Hydrauliknetzwerk;
Fig. 2 verschiedene Schritte bei der Durchführung des vorliegenden Verfahrens;
Fig. 3 ein Beispiel zur Darstellung des in einer Ausgestaltung des Verfahrens eingesetzten Glättungsalgorithmus; und
Fig. 4 unterschiedliche Hydraulikstrukturen, die mit unterschiedlichen Vernetzungsparametern des vorliegenden Verfahrens erstellt wurden.
Wege zur Ausführung der Erfindung Fig. 1 zeigt beispielhaft eine Gegenüberstellung zwischen zwei bekannten Hydraulikstrukturen und einer mit dem vorliegenden Verfahren erstellten Hydraulikstruktur. In Teilabbildung a) ist die bekannte mäanderförmige Anordnung eines Fluidkanals 2 an einer
Wärmeübertragungsfläche 1 zu erkennen. In einer derartigen Struktur können allerdings hohe Druckverluste entstehen. Teilabbildung b) zeigt eine parallele Anordnung von Stromungskanalen 2, die jeweils in Sammelkanäle 3 münden. Bei dieser Ausgestaltung können aufgrund der Druckverluste beim Übergang von den Sammelkanälen 3 bzw. in die Sammelkanäle 3 ungleichmäßige Volumenstromverteilungen auftreten.
Teilabbildung c) zeigt schließlich ein mehrfach verzweigtes Hydrauliknetzwerk, das mit dem vorliegenden Verfahren erstellt wurde. Hierbei sind sehr gut die
einzelnen Verzweigungen 4 zu erkennen, an denen der Fluidstrom jeweils in zwei Teilkanäle 2 aufgeteilt wird. Aufgrund der kleinen Winkel dieser Verzweigungen sowie der guten Ausfüllung der Wärmeübertragungsfläche 1 durch die einzelnen Teilkanäle 2 kommt es einerseits zu sehr geringen Druckverlusten und andererseits zu einer hohen thermischen Effizienz bei der Wärmeübertragung von der oder auf die Wärmeübertragungs- flache 1.
Fig. 2 zeigt beispielhaft eine Möglichkeit zum Erstellen eines derartigen mehrfach verzweigten Hydrauliknetzwerkes gemäß dem vorliegenden Verfahren.
Zunächst werden der Fluideintritts- 5 und
Fluidaustrittspunkt 6 an der Übertragungsfläche 1 festgelegt (Teilabbildung 1) . Diese Punkte werden miteinander verbunden. Anschließend werden von beiden Punkten 5, 6 aus jeweils im Abstand s0 (vorgegebene Schrittweite) die Senkrechten auf der entstandenen
Strecke bestimmt (Teilabbildung 2) . Diese Senkrechten schneiden den Rand 11 der Übertragungsfläche 1 in jeweils zwei Punkten. Die Strecke zwischen diesen Schnittpunkten wird halbiert. Die Verbindung vom Fluideintritts- 5 und Fluidaustrittspunkt 6 zu diesem Mittelpunkt bilden jeweils das erste Teilstück des ersten Pfades 7 (von zwei Seiten aus betrachtet) . Nun werden diese Verfahrensschritte - ausgehend von den neu entstandenen Punkten, also den Mittelpunkten - wiederholt: Bildung der Verbindungsstrecke, der
Senkrechten, Mittelpunktbestimmung, neues Pfadstück des ersten Pfades 7 (Teilabbildung 3)).
Diese Vorgehensweise wird wiederholt, bis der Abstand der Punkte gering genug geworden ist (Teilabbildungen 4 bis 6) . Schließlich werden die letzten Punkte miteinander verbunden. Das Ziel bei der Ermittlung dieses ersten Pfades 7 besteht darin, einen möglichst gleichen Abstand zum rechten und linken Rand 11 der Übertragungsfläche 1 zu erreichen.
Im nächsten Schritt (Teilabbildung 7) wird der Start- und Endpunkt für die nächste Iteration, d.h. die erste Verzweigung, bestimmt. Hierzu verfolgt man in Strömungsrichtung Teilstück für Teilstück des erzeugten ersten Pfades 7 (0. Iterationsstufe) und ermittelt jeweils die Abstände zum rechten und linken Rand 11 (dre bzw. dχι) . Sobald einer dieser Abstände über einem vorgegebenen Maximalwert dmaX(0 liegt, wird der Startpunkt 8 für die erste Verzweigung gesetzt. Nun läuft man auf dem Pfad weiter, bis einer der beiden Abstände unter einem vorgegebenen Minimalwert dmin,0 (z.B. dmin,o = 0,8 -dmax^o) liegt. Dort wird der Endpunkt 9 der Verzweigung gesetzt. Je nach Form der Fläche können auch mehrere Verzweigungen auftreten, beispielsweise vor und nach einer Einschnürung. Nach Bestimmung der Verzweigungspunkte 8, 9 wird der erste Pfad 7 geglättet (Teilabbildung 8) .
Nun wird das Verfahren analog zu den beschriebenen Schritten in der vorangehenden 0. Iterationsstufe - beginnend bei den zuvor berechneten Start- und Endpunkten 8, 9 - durchgeführt (Teilabbildung 9) .
Zunächst werden dazu wiederum die Punkte 8, 9 verbunden und die Senkrechten im Abstand Si bestimmt (der Index bezeichnet die Iterationsstufe; die Schrittweite s kann
sich in jeder Iteration ändern) . Die jeweilige Senkrechte schneidet den ersten Pfad 7 der 0. Iterationsstufe. Mit diesem entstehendem Schnittpunkt und den Schnittpunkten mit dem Rand 11 wird die Senkrechte in eine rechte und eine linke Strecke geteilt. Die Verbindung vom Start- 8 bzw. Endpunkt 9 zu den Mittelpunkten dieser Strecken bilden die ersten Teilstücke auf den Pfaden 10 der ersten Iterationsstufe (rechter und linker Pfad 10) . Das Verfahren wird entsprechend fortgesetzt, wobei nun die rechte bzw. linke Teilfläche als Begrenzung berücksichtigt werden. Sobald die Pfade 10 der ersten Iterationsstufe erzeugt und die neuen Start- und Endpunkte ermittelt worden sind, wird wiederum - unter Berücksichtigung des Übergangs zum ersten Pfad 7 der 0. Iterationsstufe - der Glättungsalgorithmus angewandt (Teilabbildung 10) .
Die beschriebenen Schritte werden fortgesetzt, bis eine vorgegebene Iterationstiefe oder ein Abbruch- kriterium (z.B. ein minimaler Pfadabstand) erreicht wurde .
Die Änderung des Kanalquerschnitts von einer Iterationsstufe zur nächsten stellt einen weiteren Parameter dar, der das Strömungsverhalten beeinflusst. Gegebenenfalls kann er sich auch innerhalb einer Iterationsstufe - d.h. bei parallelen Kanälen - unterscheiden, um unterschiedliche Kanallängen auszugleichen (hydraulischer Abgleich) .
Die Funktionsweise des in diesem Beispiel eingesetzten Glättungsalgorithmus soll exemplarisch anhand von Fig. 3 erklärt werden. Ausgehend von dem in
Teilabbildung 1 dargestellten Pfad, der aus Teilstücken der Längen aι-a4 besteht, werden zunächst für jedes Teilstück neue Punkte bestimmt, die jeweils den Abstand f • ai vom Anfangs- bzw. Endpunkt jedes Teilstücks haben, wobei f ein voreingestellter Faktor ist (günstig ist z.B. f = 0,3) . Auf diese Weise wird jedes Teilstück in drei Teile geteilt. Im ersten Glättungsschritt werden nun diese neuen Punkte miteinander verbunden, so dass nur die Mittelteile der ursprünglichen Teilstücke übernommen und die bisherigen Ecken ignoriert werden (Teilabbildung 2) . Dieses Verfahren wird mit dem jeweils neu entstandenen Pfad wiederholt, bis die gewünschte Glättung erreicht ist (Teilabbildüngen 3 und 4) .
Fig. 4 zeigt exemplarisch den Einfluss der Vernetzungsparameter auf das Aussehen der erzeugten Hydraulikstrukturen. So wurden in den Varianten A) bis D) die Schrittweite s0, der Faktor fs (Si+i= fs • Sj.) , der Maximalabstand dmaX/0, der Faktor f (dmax,i+ι =f • dmax,i) sowie die Anzahl der Iterationen n variiert. Der oben beschriebene Faktor f für den Glättungsalgorithmus war in allen Fällen gleich (f=0,3). Es ist erkennbar, dass je nach Wahl der Parameter entweder die Flächen- ausnutzung (Variante A) oder die günstige Hydraulik durch kleine Winkel bei Strömungsumlenkungen (Variante D) im Mittelpunkt stehen. Für die jeweilige Anwendung gilt es, die besten Parameter im Hinblick auf die Gesamteffizienz zu finden.
Auch wenn in dem vorangehenden Ausführungsbeispiel lediglich eine zweidimensionale Hydraulikstrutur erläutert wurde, so sind selbstverständlich mit dem
vorliegenden Verfahren auch dreidimensionale, baumartige Strukturen realisierbar. Hierbei kann beispielsweise analog zum Mittelpunkt einer Strecke im zweidimensionalen Fall der Schwerpunkt einer Schnittfläche im dreidimensionalen Fall bestimmt werden. Statt der rechten und linken Teilstrecke im zweidimensionalen Fall, wie dies anhand der Fig. 2 erläutert wurde, werden hier beispielsweise drei flächengleiche Teilflächen berücksichtigt. Entsprechend verzweigt sich in diesem Beispiel dann ein Kanal an einem Verzweigungspunkt in drei Teilkanäle. Die damit erzeugbare dreidimensionale Hydraulikstruktur erzeugt die gleichen Vorteile an einem Ubertragungsvolumen wie dies mit der zweidimensionalen Hydraulikstruktur an einer Wärmeübertragungsfläche erreicht wird.
BEZUGSZEICHENLISTE