500 Naturwissenschaften und Mathematik
Refine
Document Type
- Doctoral Thesis (41)
- Article (23)
- Bachelor Thesis (5)
- Master's Thesis (5)
- Conference Proceeding (2)
- Habilitation (1)
- Preprint (1)
- Study Thesis (1)
Institute
- Institut für Strukturmechanik (ISM) (33)
- F. A. Finger-Institut für Baustoffkunde (FIB) (6)
- Professur Angewandte Mathematik (5)
- Professur Werkstoffe des Bauens (4)
- Bauhaus-Institut für zukunftsweisende Infrastruktursysteme (b.is) (3)
- Professur Bauchemie und Polymere Werkstoffe (3)
- Professur Modellierung und Simulation - Konstruktion (3)
- Professur Modellierung und Simulation - Mechanik (3)
- Professur Stochastik und Optimierung (3)
- Institut für Konstruktiven Ingenieurbau (IKI) (2)
Keywords
- Finite-Elemente-Methode (6)
- Isogeometric Analysis (6)
- Isogeometrische Analyse (5)
- Maschinelles Lernen (5)
- Gestaltoptimierung (4)
- Machine learning (4)
- OA-Publikationsfonds2020 (4)
- Peridynamik (4)
- Beton (3)
- Modellierung (3)
Ein aktuelles Thema in der Forschung der Betonindustrie ist die gezielte Steuerung des Erstarrens und der Entwicklung der (Früh)Festigkeit von Betonen und Mörteln. Aus ökonomischer Sicht sind außerdem die Reduktion der CO2-Emission und die Schonung von Ressourcen und Energie wichtige Forschungsschwerpunkte. Eine Möglichkeit zum Erreichen dieser Ziele ist es, die Reaktivität/Hydratation der silikatischen Klinkerphasen gezielt anzuregen. Neben den bereits bekannten Möglichkeiten der Hydratationsbeschleunigung (u.a. Wärmebehandlung, Zugabe von Salzen) bietet die Anwendung von Power-Ultraschall (PUS) eine weitere Alternative zur Beschleunigung der Zementhydratation. Da bis zum jetzigen Zeitpunkt noch keine Erfahrungen zum Einsatz von PUS in der Zementchemie vorliegen, sollen mit der vorliegenden Arbeit grundlegende Kenntnisse zum Einfluss von PUS auf das Fließ- und Erstarrungsverhalten von Zementsuspensionen erarbeitet werden.
Dazu wurde die Arbeit in fünf Hauptuntersuchungsabschnitte aufgeteilt.
Im ersten Teil wurden optimale PUS-Parameter wie Amplitude und Energieeintrag ermittelt, die eine effiziente Beschleunigung der Portlandzement(CEM I)hydratation bei kurzen Beschallzeiten und begrenzter Zementleimtemperaturerhöhung erlauben. Mit Hilfe unabhängiger Untersuchungsmethoden (Bestimmung des Erstarrungsbeginns, der Festigkeitsentwicklung, zerstörungsfreier Ultraschallprüfung, isothermer Wärmeflusskalorimetrie, hochauflösender Rasterelektronmikroskopie (REM) wurde die Wirkung von PUS auf den Hydratationsverlauf von CEM I-Suspensionen charakterisiert. Die Ergebnisse zeigen, dass die Behandlung von CEM I-Suspensionen mit PUS grundsätzlich ein beschleunigtes Erstarren und eine beschleunigte (Früh)Festigkeitsentwicklung hervorruft.
Anhand von REM-Untersuchungen konnte eindeutig nachgewiesen werden, dass die Beschleunigung der CEM I-Hydratation mit einer beschleunigten Hydratation der Hauptklinkerphase Alit korreliert. Auf Grundlage dieser Erkenntnisse wurden die Ursachen der Aktivierung der Alithydratation untersucht. Dazu wurden Untersuchungen an Einzelsystemen des CEM I (silikatische Klinkerphase) durchgeführt.
Es ist bekannt, das die Hydratation der Hauptklinkerphase Alit (in der reinen Form Tricalciumsilikat 3CaO*SiO2; C3S) durch Lösungs-/Fällungsreaktionen (Bildung von Calcium-Silikat-Hydrat Phasen, C-S-H Phasen) bestimmt wird. Mit Hilfe von Untersuchungen zur Auflösung (C3S) und Kristallbildung (C-S-H Phasen) in Lösungen und Suspensionen (Aufzeichnung der elektrischen Leitfähigkeit sowie Bestimmung der Ionenkonzentrationen der wässrigen Phase, REM-Charakterisierung der Präzipitate) wurde die Beeinflussung dieser durch eine PUS-Behandlung charakterisiert. Die Ergebnisse zeigen, dass in partikelfreien Lösungen (primäre Keimbildung) eine PUS-Behandlung keinen Einfluss auf die Kinetik der Kristallisation von C-S-H Phasen hervorruft. Das heißt, auch die durch PUS eingetragene Energie reicht offensichtlich nicht aus, um in Abwesenheit von Oberflächen die C-S-H Phasen Bildung zu beschleunigen. Das weist darauf hin, dass die Bildung von C-S-H Phasen nicht durch eine Beschleunigung von Ionen in der Lösung (erhöhte Diffusion durch Anwendung von PUS) hervorgerufen wird. Eine Beschleunigung des Kristallisationsprozesses (Keimbildung und Wachstum von C-S-H Phasen) durch PUS wird nur in Anwesenheit von Partikeln in der Lösung (Suspension) erzielt. Das belegen Ergebnisse, bei denen die Bildung erster C-S-H Phasen bei geringer Übersättigung (heterogene Keimbildung, in Anwesenheit von Oberflächen) erfolgt. Unter diesen Bedingungen konnte gezeigt werden, dass PUS innerhalb der ersten 30 Minuten der Hydratation eine erhöhte Fällung von ersten C-S-H Phasen bewirkt. Diese fungieren dann vermutlich während der Haupthydratation als Keim bzw. geeignete Oberfläche zum beschleunigten Aufwachsen von weiteren C-S-H Phasen. Weiterhin ist vorstellbar, dass (in Analogie zu anderen Bereichen der Sonochemie) PUS durch Kavitation Schockwellen hervorruft, welche Partikel und wässriges Medium beschleunigen und damit erhöhte Partikelbewegungen und -kollisionen induziert. Dies wiederum bewirkt, dass die anfänglich auf der C3S-Oberfläche gebildeten C-S-H Phasen teilweise wieder entfernt werden. Damit ist das Inlösunggehen von Ca- und Si-Ionen aus dem C3S weiterhin möglich. Um den genauen Mechanismus weiter zu charakterisieren sollten mit geeigneten Methoden weitere Untersuchungen durchgeführt werden.
Im zweiten Teil der Arbeit wurde der Einfluss von PUS auf das Fließverhalten von CEM I-Suspensionen untersucht. Aus der Anwendung von PUS in anderen technischen Bereichen sind unter anderem Effekte wie das Entlüften, das Homogenisieren und das Dispergieren von Suspensionen und Emulsionen mittels PUS bekannt. Mit Hilfe der Bestimmung des Luftporengehaltes, Sedimentationsversuchen und cryo-SEM Untersuchungen wurde der Einfluss von PUS auf CEM I-Suspensionen charakterisiert. Die Ergebnisse belegen, dass durch PUS eine verbesserte Homogenität und Dispergierung der CEM I-Suspension erzielt wird. Damit wird für CEM I-Suspensionen unterschiedlichster w/z-Werte eine verbesserte Fließfähigkeit festgestellt. Ergebnisse der Bestimmung von Ausbreitmaßen und Trichterauslaufzeiten zeigen, dass PUS einen direkten Einfluss vor allem auf die Viskosität der CEM I-Suspensionen besitzt. Werden Fließmitteln (FM) der CEM I-Suspension zugegeben, wird nicht in jedem Fall eine verbesserte Fließfähigkeit festgestellt. Hier scheint unter bestimmten Voraussetzungen (w/z-Wert, FM-Gehalt, PUS) die Reaktion zwischen Aluminat- und Sulfatphase des Klinkers gestört. Zur eindeutigen Klärung dieses Sachverhaltes bedarf es jedoch weiterer quantitativer Untersuchungen zum Reaktionsumsatz.
Im dritten Teil der Arbeit wurden die am CEM I gewonnenen Erkenntnisse zum Einfluss von PUS auf die Hydratation an Portland-Hüttensand(HÜS)-Zement-Systemen verifiziert. Dafür
wurden auch in diesem Teil der Arbeit zunächst die optimalen PUS-Parameter festgelegt und der Einfluss auf das Erstarrung- und Erhärtungsverhalten dokumentiert. Untersuchungsmethoden sind unter anderem die Bestimmung des Erstarrungsbeginns und der (Früh)Festigkeitsentwicklung, Temperaturaufzeichnungen und isothermale Wärmeflusskalorimetrie sowie REM. Die Ergebnisse zeigen, dass auch die Reaktion von HÜS-Zementen durch PUS beschleunigt wird. Weiterführende Untersuchungen belegen, dass die erzielte Beschleunigung vorwiegend auf der Beschleunigung der Alitkomponente des CEM I beruht.
Im Fokus der Teile vier und fünf dieser Arbeit stand die Anwendbarkeit der PUS-Technik unter praktischen Bedingungen. Zum einen wurde die Anwendbarkeit von PUS in fertig gemischten Mörteln beurteilt. Anhand des Vergleichs wichtiger Frisch- und Festmörteleigenschaften unterschiedlich hergestellter Mörtel (beschallt im Anschluss an konventionelle Mischtechnik, beschallt im Anschluss an Suspensionsmischtechnik mit anschließender Zumischung der Gesteinskörnung und nicht beschallt) wird gezeigt, dass im Fall von Mörteln mit hohem Leimanteil eine durch PUS induzierte beschleunigte Festigkeitsentwicklung auch mit herkömmlichen Mischabläufen (ohne aufwendige Umstellung des Mischprozesses) möglich ist.
Abschließend wird untersucht, ob der Herstellungsprozess von Wandbauteilen im Fertigteilwerk durch den Einsatz von PUS optimiert werden kann und ob eine Einbindung der PUS-Technik in den Fertigungsprozess ohne größeren Aufwand möglich ist. Dazu wurden in einem ersten Schritt die Frisch- und Festbetoneigenschaften eines aktuell angewendeten selbstverdichtenden Betons im Labormaßstab (Mörtel) in Abhängigkeit einer PUS-Behandlung dokumentiert und mit der seiner unbeschallten Referenz verglichen. Aufgrund der durch PUS verursachten verbesserten Fließ- und Festigkeitseigenschaften kann die beschallte Mörtelrezeptur hinsichtlich Fließmittelgehalt und Dauer der Wärmebehandlung optimiert werden. Somit werden ca. 30 % der Fließmittelzugabe und 40 % der Dauer der Wärmebehandlung eigespart. Eine Einbindung der PUS-Technik in das betrachtete Fertigteilwerk ist nach Überprüfung der konstruktiven Gegebenheiten der Fertigungsstrukturen ohne größeren Aufwand möglich.
Bei einem marktüblichen Calciumsulfat-Fließestrich wurden in der Praxis schädigende Volu-menexpansionen festgestellt. Diese sind ein Resultat aus dem Zusammenwirken des einge-setzten Bindemittel-Compounds und einer kritischen Gesteinskörnung.
Das Ziel dieser Arbeit ist es, ein Calciumsulfat-Bindemittelsystem zu konfektionieren, welches in der Lage ist, die im Mörtel festgestellten Volumenexpansionen zu unterbinden. Es sollen verschiedene Bindemittel- und Additivzusammensetzungen untersucht werden, welche in Verbindung mit der kritischen Gesteinskörnung die Herstellung eines volumenstabilen Fließestrichs ermöglichen. Dazu soll folgende Fragestellung beantwortet werden: Welche Ursachen hat die Volumenzunahme und wie ist diese zu minimieren bzw. unterbinden?
Dabei werden unterschiedliche Bindemittelrezepturen aus α-Halbhydrat, Thermoanhydrit und Naturanhydrit, sowie verschiedene Additivzusammensetzungen hergestellt und untersucht.
Durch Längenänderungsmessungen in der Schwindrinne werden die Einflüsse der Binde-mittel, der Additivzusammensetzungen und der Wasser/Bindemittel-Werte auf das Län-genänderungsverhalten untersucht. Mittels Variation der einzelnen Compound-Bestandteile kann festgestellt werden, dass der Stabilisierer die Längenänderung negativ beeinflusst. Dieser bindet freies Wasser, welches für eine Reaktion zwischen Bindemittel und Gesteins-körnung im plastischen Zustand nicht mehr zur Verfügung steht. Diese Reaktion kann folglich erst im erhärteten Zustand ablaufen und verursacht die schädigende Volumenexpansion.
Abschließend wurde ein Bindemittel-Compound konfektioniert, welcher ohne Zusatz von Stabilisierern in Zusammenhang mit der kritischen Gesteinskörnung volumenstabil ist und keine Schäden auslöst.
We present a stochastic deep collocation method (DCM) based on neural architecture search (NAS) and transfer learning for heterogeneous porous media. We first carry out a sensitivity analysis to determine the key hyper-parameters of the network to reduce the search space and subsequently employ hyper-parameter optimization to finally obtain the parameter values. The presented NAS based DCM also saves the weights and biases of the most favorable architectures, which is then used in the fine-tuning process. We also employ transfer learning techniques to drastically reduce the computational cost. The presented DCM is then applied to the stochastic analysis of heterogeneous porous material. Therefore, a three dimensional stochastic flow model is built providing a benchmark to the simulation of groundwater flow in highly heterogeneous aquifers. The performance of the presented NAS based DCM is verified in different dimensions using the method of manufactured solutions. We show that it significantly outperforms finite difference methods in both accuracy and computational cost.
The vibration control of the tall building during earthquake excitations is a challenging task due to their complex seismic behavior. This paper investigates the optimum placement and properties of the Tuned Mass Dampers (TMDs) in tall buildings, which are employed to control the vibrations during earthquakes. An algorithm was developed to spend a limited mass either in a single TMD or in multiple TMDs and distribute them optimally over the height of the building. The Non-dominated Sorting Genetic Algorithm (NSGA – II) method was improved by adding multi-variant genetic operators and utilized to simultaneously study the optimum design parameters of the TMDs and the optimum placement. The results showed that under earthquake excitations with noticeable amplitude in higher modes, distributing TMDs over the height of the building is more effective in mitigating the vibrations compared to the use of a single TMD system. From the optimization, it was observed that the locations of the TMDs were related to the stories corresponding to the maximum modal displacements in the lower modes and the stories corresponding to the maximum modal displacements in the modes which were highly activated by the earthquake excitations. It was also noted that the frequency content of the earthquake has significant influence on the optimum location of the TMDs.
In the last two decades, Peridynamics (PD) attracts much attention in the field of fracture mechanics. One key feature of PD is the nonlocality, which is quite different from the ideas in conventional methods such as FEM and meshless method. However, conventional PD suffers from problems such as constant horizon, explicit algorithm, hourglass mode. In this thesis, by examining the nonlocality with scrutiny, we proposed several new concepts such as dual-horizon (DH) in PD, dual-support (DS) in smoothed particle hydrodynamics (SPH), nonlocal operators and operator energy functional. The conventional PD (SPH) is incorporated in the DH-PD (DS-SPH), which can adopt an inhomogeneous discretization and inhomogeneous support domains. The DH-PD (DS-SPH) can be viewed as some fundamental improvement on the conventional PD (SPH). Dual formulation of PD and SPH allows h-adaptivity while satisfying the conservations of linear momentum, angular momentum and energy. By developing the concept of nonlocality further, we introduced the nonlocal operator method as a generalization of DH-PD. Combined with energy functional of various physical models, the nonlocal forms based on dual-support concept are derived. In addition, the variation of the energy functional allows implicit formulation of the nonlocal theory. At last, we developed the higher order nonlocal operator method which is capable of solving higher order partial differential equations on arbitrary domain in higher dimensional space. Since the concepts are developed gradually, we described our findings chronologically.
In chapter 2, we developed a DH-PD formulation that includes varying horizon sizes and solves the "ghost force" issue. The concept of dual-horizon considers the unbalanced interactions between the particles with different horizon sizes. The present formulation fulfills both the balances of linear momentum and angular momentum exactly with arbitrary particle discretization. All three peridynamic formulations, namely bond based, ordinary state based and non-ordinary state based peridynamics can be implemented within the DH-PD framework. A simple adaptive refinement procedure (h-adaptivity) is proposed reducing the computational cost. Both two- and three- dimensional examples including the Kalthoff-Winkler experiment and plate with branching cracks are tested to demonstrate the capability of the method.
In chapter 3, a nonlocal operator method (NOM) based on the variational principle is proposed for the solution of waveguide problem in computational electromagnetic field. Common differential operators as well as the variational forms are defined within the context of nonlocal operators. The present nonlocal formulation allows the assembling of the tangent stiffness matrix with ease, which is necessary for the eigenvalue analysis of the waveguide problem. The present formulation is applied to solve 1D Schrodinger equation, 2D electrostatic problem and the differential electromagnetic vector wave equations based on electric fields.
In chapter 4, a general nonlocal operator method is proposed which is applicable for solving partial differential equations (PDEs) of mechanical problems. The nonlocal operator can be regarded as the integral form, ``equivalent'' to the differential form in the sense of a nonlocal interaction model. The variation of a nonlocal operator plays an equivalent role as the derivatives of the shape functions in the meshless methods or those of the finite element method. Based on the variational principle, the residual and the tangent stiffness matrix can be obtained with ease. The nonlocal operator method is enhanced here also with an operator energy functional to satisfy the linear consistency of the field. A highlight of the present method is the functional derived based on the nonlocal operator can convert the construction of residual and stiffness matrix into a series of matrix multiplications using the predefined nonlocal operators. The nonlocal strong forms of different functionals can be obtained easily via the concept of support and dual-support. Several numerical examples of different types of PDEs are presented.
In chapter 5, we extended the NOM to higher order scheme by using a higher order Taylor series expansion of the unknown field. Such a higher order scheme improves the original NOM in chapter 3 and chapter 4, which can only achieve one-order convergence. The higher order NOM obtains all partial derivatives with specified maximal order simultaneously without resorting to shape functions. The functional based on the nonlocal operators converts the construction of residual and stiffness matrix into a series of matrix multiplication on the nonlocal operator matrix. Several numerical examples solved by strong form or weak form are presented to show the capabilities of this method.
In chapter 6, the NOM proposed as a particle-based method in chapter 3,4,5, has difficulty in imposing accurately the boundary conditions of various orders. In this paper, we converted the particle-based NOM into a scheme with interpolation property. The new scheme describes partial derivatives of various orders at a point by the nodes in the support and takes advantage of the background mesh for numerical integration. The boundary conditions are enforced via the modified variational principle. The particle-based NOM can be viewed a special case of NOM with interpolation property when nodal integration is used. The scheme based on numerical integration greatly improves the stability of the method, as a consequence, the operator energy functional in particle-based NOM is not required. We demonstrated the capabilities of current method by solving the gradient solid problems and comparing the numerical results with the available exact solutions.
In chapter 7, we derived the DS-SPH in solid within the framework of variational principle. The tangent stiffness matrix of SPH can be obtained with ease, and can be served as the basis for the present implicit SPH. We proposed an hourglass energy functional, which allows the direct derivation of hourglass force and hourglass tangent stiffness matrix. The dual-support is {involved} in all derivations based on variational principles and is automatically satisfied in the assembling of stiffness matrix. The implementation of stiffness matrix comprises with two steps, the nodal assembly based on deformation gradient and global assembly on all nodes. Several numerical examples are presented to validate the method.
Im Rahmen der Dissertation wurde ein Toolboxmodell für transdisziplinäres Wasserressourcenmanagement entwickelt. Das Modell liefert den methodischen Rahmen Wasserressourcen nachhaltig und transdisziplinär zu bewirtschaften.
Der Begriff der Nachhaltigkeit und eine Konkretisierung der nachhaltigen Bewirtschaftung globaler Wasserressourcen scheinen unüberschaubar und suggerieren die Forderung nach einer neuen Weltformel. Die globale Bedeutung der Wasserressourcen, die für Regionen spezifischen Besonderheiten des natürlichen Wasserhaushalts und der anthropogenen Nutzung, die Zeitskala und die Kontextualisierung in alle betroffenen und benachbarten Disziplinen deuten auf die Komplexität der Thematik hin. Es wird eine Systematisierung des Planungsprozesses von Wasserressourcen notwendig, anhand derer eine holistische Herangehensweise mit einer Strategieentwicklung für Regionen spezifischer Schwerpunktprobleme erfolgt. Ziel der Arbeit ist die Erarbeitung einer Strategie zur Systematisierung nach diesen Forderungen und die Bereitstellung eines Toolboxmodelles als Planungswerkzeug für das transdisziplinäre Wasserressourcenmanagement.
Das Toolboxmodell stellt den konzeptionellen Rahmen für die Bewirtschaftung von Wasserressourcen mit der Anwendung transdisziplinärer Forschungsmethoden bereit. Wesentliche Herausforderung bei der Anwendung der transdisziplinären Methode sind die Implementierung verschiedener Skalenbereiche, der Umgang mit der Komplexität von Daten, das Bewahren von Transparenz und Objektivität sowie die Ermöglichung eines auf andere Regionen übertragbaren Planungsprozesses.
Die theoretischen Grundlagen naturwissenschaftlicher Forschung zur Nachhaltigkeit haben ihren Ursprung in den biologischen und geographischen Disziplinen. Das Ineinandergreifen naturräumlicher Zusammenhänge und der Einfluss anthropogener Nutzung und technischer Innovationen auf den Naturhaushalt sind Kern der Kausalität übergreifenden Denkens und Verstehens. Mit dem Ansatz des integrierten Wasserressourcenmanagements (IWRM) erfolgt die Berücksichtigung wirtschaftlicher und sozioökonomischer Ziele in den Planungsprozess für ökologisch nachhaltige Wasserwirtschaft. Das Instrument der Wasserrahmenrichtlinie (EU-WRRL) ist auf eine Gewässerökologie ausgerichtete Richtlinie, welche die Integration verschiedener Interessenvertreter in den Planungsprozess vorsieht. Das Konzept der neuartigen Sanitärsysteme basiert auf Stoffflüssen zwischen konkurrierenden Handlungsbereichen, wie Abfall-, Ressourcen- und Landwirtschaft.
Den integrierten Ansätzen fehlt eine übergeordnete gemeinsame Zielstrategie – eine sogenannte Phase Null. Diese Phase Null – das Lernen aller 7 Zusammenfassung 157 relevanten, konkurrierenden und harmonisierenden Handlungsfelder eines Planungshorizontes wird durch eine transdisziplinäre Perspektive ermöglicht. Während bei der integralen Perspektive eine disziplinorientierte Kooperation im Vordergrund steht, verlangt die transdisziplinäre Perspektive nach einer problemorientierten Kooperation zwischen den Interessenvertretern (Werlen 2015). Die bestehenden Konzepte und Richtlinien für das nachhaltige Management von Wasserressourcen sind etabliert und evaluiert. Der Literatur zur Folge ist eine Weiterentwicklung nach der Perspektive der Transdisziplinarität erforderlich. Das Toolboxmodell für integrales Wasserressourcenmanagement entspricht einem Planungstool bestehend aus Werkzeugen für die Anwendung wissenschaftlicher Methoden. Die Zusammenstellung der Methoden/Werkzeuge erfüllt im Rahmen die Methode transdisziplinärer Forschung. Das Werkzeug zum Aufstellen der relevanten Handlungsfelder umfasst die Charakterisierung eines Untersuchungsgebietes und Planungsrahmens, die kausale Verknüpfung des Bewirtschaftungskonzeptes und konkurrierender sowie sich unterstützender Stakeholder. Mit dem Werkzeug der Kontextualisierung und Indikatorenaufstellung wird eine Methode der stufenweisen und von einer Skala unabhängigen Bewertung des Umweltzustandes für die Zielpriorisierung vorgenommen. Damit wird das Toolboxmodell dem Problem der Komplexität und Datenverfügbarkeit gerecht. Anhand der eingesetzten ABC Methode, werden die Bewertungsgrößen differenziert strukturiert auf verschiedene Skalen und Datenressourcen (A=Ersterkennung,B=Zeigerwerte, C=Modell/Index). Die ABC-Methode ermöglicht die Planung bereits mit unsicherer und lückenhafter Datengrundlage, ist jederzeit erweiterbar und bietet somit eine operative Wissensgenerierung während des Gestaltungsprozesses.
Für das Werkzeug zur Bewertung und Priorisierung wird der Algorithmus der Composite Programmierung angewandt. Diese Methode der Mehrfachzielplanung erfüllt den Anspruch der permanenten Erweiterbarkeit und der transparenten und objektiven Entscheidungsfindung. Die Komplexität des transdisziplinären Wasserressourcenmanagements kann durch die Methode der Composite Programmierung systematisiert werden. Das wesentliche Ergebnis der Arbeit stellt die erfolgreiche Erarbeitung und Anwendung des Tool-boxmodells für das transdisziplinäre Wasserressourcenmanagement im Untersuchungsgebiet Stadt Darkhan in der Mongolei dar. Auf Grund seiner besonderen hydrologischen und strukturellen Situa-tion wird die Relevanz eines nachhaltigen Bewirtschaftungskonzeptes deutlich. Im Rahmen des Querschnittsmoduls des MoMo-Projektes wurde eine für das Toolboxmodell geeignete Datengrundlage erarbeitet. Planungsrelevante Handlungsfelder wurden im Rahmen eines Workshops mit verschiedenen Interessenvertretern erarbeitet. Im Ergebnis dessen wurde die Systematik eines Zielbaumes mit Hauptzielen und untergeordneten Teilzielen als Grundlage der Priorisierung nach den holistischen Anspruch der transdisziplinären Forschung aufgestellt. Für die Messbarkeit, in-wieweit Teilziele erreicht sind oder Handlungsbedarf besteht, wurden Indikatoren erarbeitet. Die Indikatoren-Aufstellung erfolgte exemplarisch für das Handlungsfeld Siedlungswasserwirtschaft in allen Skalen des ABC-Systems. Die im BMBF-MoMo Projekt generierte umfassende Datengrundlage ermöglichte die Anwendung und Evaluierung des Toolboxmodells mit unterschiedlichem quantitativem und qualitativem Dateninput. Verschiedene Kombination von A (Ersterkennung), B (Zeigerwerte) und C (Modell/Index) als Grundlage der Priorisierung mit der Compostite Programmierung ermöglichten die Durchführung und Bewertung des transdisziplinären Planungstools. Die er-mittelten Rangfolgen von Teilzielen mit unterschiedlichen Bewertungsvarianten ergaben ähnliche
Tendenzen. Das ist ein Hinweis dafür, dass für die zukünftige Anwendung des Toolboxmodells die operative Wissensgenerierung, d.h. das schrittweise Hinzufügen neu ermittelter, gesicherterer Daten, funktioniert. Eine schwierige Datenverfügbarkeit oder eine noch im Prozess befindliche wissenschaftliche Analyse sollen keine Hindernisse für eine schrittweise und erweiterbare Zielpriorisierung und Maßnahmenplanung sein. Trotz der Komplexität des transdisziplinären Ansatzes wird durch die Anwendung des Toolboxmodells eine effiziente und zielorientierte Handlungspriorisierung ermöglicht. Die Effizienz wird erreicht durch ressourcenschonende und flexible, Ziel fokussierte Datenermittlung. Zeit und Kosten im Planungsprozess können eingespart werden. Die erzielte Priorisierung von letztlich Handlungsempfehlungen erfolgt individuell auf die Eigenart des Untersuchungsgebietes angepasst, was hinsichtlich seiner Wirkung als erfolgsversprechend gilt.
Das Kernthema dieser Arbeit ist die Beschäftigung mit den Folgen des Uranbergbaus in dem Gebiet um die ehemalige Abbauregion der Wismut SAG/SDAG in Ronneburg (Ostthüringen). Dieses Thema wird unter historischen, sozialen, kulturanthropologischen und künstlerischen Aspekten betrachtet und in den Zusammenhang mit den weltweiten Voraussetzungen der Nuklearindustrie und Auswirkungen des Uranbergbaus und seiner Folgen gestellt. Die Arbeit legt dar, wie eine Uranbergbaufolgelandschaft entsteht und welches Wissen ist für ein angemessenes Verständnis des Phänomens wichtig ist. Es wird untersucht, ob Kunst bezüglich der Uranbergbaufolgelandschaft einen relevanten Beitrag leisten kann bzw. in welcher Form dies versucht wurde, bzw. stellte Arbeiten vor, die verwandete Themen bearbeitet haben. In Kombination dieser beiden Hauptaspekte geht die Arbeit der Frage nach, welche Faktoren die Uranbergbaufolgelandschaft prägen und ob es sinnvolle Beteiligungsfelder für künstlerisches Forschen oder Handeln gibt sowie welche Bedingungen hierfür erfüllt werdenmüssten. Die Kernthese der Arbeit ist, dass künstlerische Arbeiten im Themenfeld des Uranbergbaus unter bestimmten Bedingungen relevante Beiträge leisten können.
The K-nearest neighbors (KNN) machine learning algorithm is a well-known non-parametric classification method. However, like other traditional data mining methods, applying it on big data comes with computational challenges. Indeed, KNN determines the class of a new sample based on the class of its nearest neighbors; however, identifying the neighbors in a large amount of data imposes a large computational cost so that it is no longer applicable by a single computing machine. One of the proposed techniques to make classification methods applicable on large datasets is pruning. LC-KNN is an improved KNN method which first clusters the data into some smaller partitions using the K-means clustering method; and then applies the KNN for each new sample on the partition which its center is the nearest one. However, because the clusters have different shapes and densities, selection of the appropriate cluster is a challenge. In this paper, an approach has been proposed to improve the pruning phase of the LC-KNN method by taking into account these factors. The proposed approach helps to choose a more appropriate cluster of data for looking for the neighbors, thus, increasing the classification accuracy. The performance of the proposed approach is evaluated on different real datasets. The experimental results show the effectiveness of the proposed approach and its higher classification accuracy and lower time cost in comparison to other recent relevant methods.
In this paper, an artificial neural network is implemented for the sake of predicting the thermal conductivity ratio of TiO2-Al2O3/water nanofluid. TiO2-Al2O3/water in the role of an innovative type of nanofluid was synthesized by the sol–gel method. The results indicated that 1.5 vol.% of nanofluids enhanced the thermal conductivity by up to 25%. It was shown that the heat transfer coefficient was linearly augmented with increasing nanoparticle concentration, but its variation with temperature was nonlinear. It should be noted that the increase in concentration may cause the particles to agglomerate, and then the thermal conductivity is reduced. The increase in temperature also increases the thermal conductivity, due to an increase in the Brownian motion and collision of particles. In this research, for the sake of predicting the thermal conductivity of TiO2-Al2O3/water nanofluid based on volumetric concentration and temperature functions, an artificial neural network is implemented. In this way, for predicting thermal conductivity, SOM (self-organizing map) and BP-LM (Back Propagation-Levenberq-Marquardt) algorithms were used. Based on the results obtained, these algorithms can be considered as an exceptional tool for predicting thermal conductivity. Additionally, the correlation coefficient values were equal to 0.938 and 0.98 when implementing the SOM and BP-LM algorithms, respectively, which is highly acceptable. View Full-Text
Material properties play a critical role in durable products manufacturing. Estimation of the precise characteristics in different scales requires complex and expensive experimental measurements. Potentially, computational methods can provide a platform to determine the fundamental properties before the final experiment. Multi-scale computational modeling leads to the modeling of the various time, and length scales include nano, micro, meso, and macro scales. These scales can be modeled separately or in correlation with coarser scales. Depend on the interested scales modeling, the right selection of multi-scale methods leads to reliable results and affordable computational cost. The present dissertation deals with the problems in various length and time scales using computational methods include density functional theory (DFT), molecular mechanics (MM), molecular dynamics (MD), and finite element (FE) methods.
Physical and chemical interactions in lower scales determine the coarser scale properties. Particles interaction modeling and exploring fundamental properties are significant challenges of computational science. Downscale modelings need more computational effort due to a large number of interacted atoms/particles. To deal with this problem and bring up a fine-scale (nano) as a coarse-scale (macro) problem, we extended an atomic-continuum framework. The discrete atomic models solve as a continuum problem using the computationally efficient FE method. MM or force field method based on a set of assumptions approximates a solution on the atomic scale. In this method, atoms and bonds model as a harmonic oscillator with a system of mass and springs. The negative gradient of the potential energy equal to the forces on each atom. In this way, each bond's total potential energy includes bonded, and non-bonded energies are simulated as equivalent structural strain energies. Finally, the chemical nature of the atomic bond is modeled as a piezoelectric beam element that solves by the FE method.
Exploring novel materials with unique properties is a demand for various industrial applications. During the last decade, many two-dimensional (2D) materials have been synthesized and shown outstanding properties. Investigation of the probable defects during the formation/fabrication process and studying their strength under severe service life are the critical tasks to explore performance prospects. We studied various defects include nano crack, notch, and point vacancy (Stone-Wales defect) defects employing MD analysis. Classical MD has been used to simulate a considerable amount of molecules at micro-, and meso- scales. Pristine and defective nanosheet structures considered under the uniaxial tensile loading at various temperatures using open-source LAMMPS codes. The results were visualized with the open-source software of OVITO and VMD.
Quantum based first principle calculations have been conducting at electronic scales and known as the most accurate Ab initio methods. However, they are computationally expensive to apply for large systems. We used density functional theory (DFT) to estimate the mechanical and electrochemical response of the 2D materials. Many-body Schrödinger's equation describes the motion and interactions of the solid-state particles. Solid describes as a system of positive nuclei and negative electrons, all electromagnetically interacting with each other, where the wave function theory describes the quantum state of the set of particles. However, dealing with the 3N coordinates of the electrons, nuclei, and N coordinates of the electrons spin components makes the governing equation unsolvable for just a few interacted atoms. Some assumptions and theories like Born Oppenheimer and Hartree-Fock mean-field and Hohenberg-Kohn theories are needed to treat with this equation. First, Born Oppenheimer approximation reduces it to the only electronic coordinates. Then Kohn and Sham, based on Hartree-Fock and Hohenberg-Kohn theories, assumed an equivalent fictitious non-interacting electrons system as an electron density functional such that their ground state energies are equal to a set of interacting electrons. Exchange-correlation energy functionals are responsible for satisfying the equivalency between both systems. The exact form of the exchange-correlation functional is not known. However, there are widely used methods to derive functionals like local density approximation (LDA), Generalized gradient approximation (GGA), and hybrid functionals (e.g., B3LYP). In our study, DFT performed using VASP codes within the GGA/PBE approximation, and visualization/post-processing of the results realized via open-source software of VESTA.
The extensive DFT calculations are conducted 2D nanomaterials prospects as anode/cathode electrode materials for batteries. Metal-ion batteries' performance strongly depends on the design of novel electrode material. Two-dimensional (2D) materials have developed a remarkable interest in using as an electrode in battery cells due to their excellent properties. Desirable battery energy storage systems (BESS) must satisfy the high energy density, safe operation, and efficient production costs. Batteries have been using in electronic devices and provide a solution to the environmental issues and store the discontinuous energies generated from renewable wind or solar power plants. Therefore, exploring optimal electrode materials can improve storage capacity and charging/discharging rates, leading to the design of advanced batteries.
Our results in multiple scales highlight not only the proposed and employed methods' efficiencies but also promising prospect of recently synthesized nanomaterials and their applications as an anode material. In this way, first, a novel approach developed for the modeling of the 1D nanotube as a continuum piezoelectric beam element. The results converged and matched closely with those from experiments and other more complex models. Then mechanical properties of nanosheets estimated and the failure mechanisms results provide a useful guide for further use in prospect applications. Our results indicated a comprehensive and useful vision concerning the mechanical properties of nanosheets with/without defects. Finally, mechanical and electrochemical properties of the several 2D nanomaterials are explored for the first time—their application performance as an anode material illustrates high potentials in manufacturing super-stretchable and ultrahigh-capacity battery energy storage systems (BESS). Our results exhibited better performance in comparison to the available commercial anode materials.
Modern digital material approaches for the visualization and simulation of heterogeneous materials allow to investigate the behavior of complex multiphase materials with their physical nonlinear material response at various scales. However, these computational techniques require extensive hardware resources with respect to computing power and main memory to solve numerically large-scale discretized models in 3D. Due to a very high number of degrees of freedom, which may rapidly be increased to the two-digit million range, the limited hardware ressources are to be utilized in a most efficient way to enable an execution of the numerical algorithms in minimal computation time. Hence, in the field of computational mechanics, various methods and algorithms can lead to an optimized runtime behavior of nonlinear simulation models, where several approaches are proposed and investigated in this thesis.
Today, the numerical simulation of damage effects in heterogeneous materials is performed by the adaption of multiscale methods. A consistent modeling in the three-dimensional space with an appropriate discretization resolution on each scale (based on a hierarchical or concurrent multiscale model), however, still contains computational challenges in respect to the convergence behavior, the scale transition or the solver performance of the weak coupled problems. The computational efficiency and the distribution among available hardware resources (often based on a parallel hardware architecture) can significantly be improved. In the past years, high-performance computing (HPC) and graphics processing unit (GPU) based computation techniques were established for the investigationof scientific objectives. Their application results in the modification of existing and the development of new computational methods for the numerical implementation, which enables to take advantage of massively clustered computer hardware resources. In the field of numerical simulation in material science, e.g. within the investigation of damage effects in multiphase composites, the suitability of such models is often restricted by the number of degrees of freedom (d.o.f.s) in the three-dimensional spatial discretization. This proves to be difficult for the type of implementation method used for the nonlinear simulation procedure and, simultaneously has a great influence on memory demand and computational time.
In this thesis, a hybrid discretization technique has been developed for the three-dimensional discretization of a three-phase material, which is respecting the numerical efficiency of nonlinear (damage) simulations of these materials. The increase of the computational efficiency is enabled by the improved scalability of the numerical algorithms. Consequently, substructuring methods for partitioning the hybrid mesh were implemented, tested and adapted to the HPC computing framework using several hundred CPU (central processing units) nodes for building the finite element assembly. A memory-efficient iterative and parallelized equation solver combined with a special preconditioning technique for solving the underlying equation system was modified and adapted to enable combined CPU and GPU based computations.
Hence, it is recommended by the author to apply the substructuring method for hybrid meshes, which respects different material phases and their mechanical behavior and which enables to split the structure in elastic and inelastic parts. However, the consideration of the nonlinear material behavior, specified for the corresponding phase, is limited to the inelastic domains only, and by that causes a decreased computing time for the nonlinear procedure. Due to the high numerical effort for such simulations, an alternative approach for the nonlinear finite element analysis, based on the sequential linear analysis, was implemented in respect to scalable HPC. The incremental-iterative procedure in finite element analysis (FEA) during the nonlinear step was then replaced by a sequence of linear FE analysis when damage in critical regions occured, known in literature as saw-tooth approach. As a result, qualitative (smeared) crack initiation in 3D multiphase specimens has efficiently been simulated.
Maturation and Structure Formation Processes in Binders with Aqueous Alkali-Silicate Solutions
(2017)
Maturation and structure formation processes can lead to crack formation in silicate and aluminosilicate binders (e.g. for coating materials...) through restricted deformation, loss of strength and thus to loss of durability. These processes are evaluated with silicate materials with an outlook on aluminosilicate binders.
Problem definition and research objectives
1. The production of Portland cement clinker causes approx. 5% to 8% of the annual man-made CO2 emissions. This is due to the usage of mainly fossil fuel (approx. 40 % of the total CO2) and because of the decarbonation of limestone as a main component of the raw meal (approx. 60 % of the total CO2).
2. Various strategies are applied in order to reduce the green-house gas-emissions, such as optimizing the process of clinker production, the use of alternative fuel and the partial substitution of the clinker in blended cement by so-called SCM (supplementary cementitious materials). Hereby blast-furnace slag, fly ash and limestone are the most used materials.
3. Quaternary systems containing three SCM simultaneously besides Portland cement contribute to the reduction of CO2 emissions due to the decrease of the clinker content. In addition, such systems allow to use blast-furnace slag and fly ash in the most economical way and provide the possibility to account for shortages of SCM on the market.
4. Blast-furnace slag and fly ash show similarities in their principal chemical compositions such that similar hydrates are formed during their reaction in presence of Portland cement. Compared to ternary systems based on blast-furnace slag or fly ash besides limestone, quaternary systems that contain both, blast-furnace slag and fly ash, simultaneously besides limestone, are expected to perform similar in terms of phase assemblage and strength development.
5. The use of SCM as cement replacing materials is limited due to their generally slower reaction compared to neat cement which also leads to lower strength development, especially in the early stage of the hydration up to 28 d. To account for this it is necessary to study the reactivity of SCM such as blast-furnace slag and fly ash in detail in order to develop strategies to enhance the reactivity and thereby the strength development of SCM-containing systems.
6. The early hydration of clinker phases is studied in detail, mainly in diluted systems. It is unclear if processes that were found to control the reaction of such model systems are also prevailing in concentrated cement pastes under realistic water-to-solid ratios. Deeper insight to this aspect is needed to better understand interactions of neat Portland cement and SCM in the first hours of hydration.
State-of-the-art
7. Increasing Ca-concentrations lead to decreasing dissolution rates of C3S and C2S in diluted systems.
8. The hydration kinetics of C3S is controlled by the interplay of undersaturation with respect to C3S and oversaturation with respect to C-S-H.
9. Increasing Al-concentrations lead to a retardation of the hydration of C3S. It is unclear if the uptake of aluminum in C-S-H to form C-(A)-S-H which has a significantly lower growth rate than pure C-S-H or a retarding effect of Al on the dissolution of C3S causes this phenomenon.
10. The surface of limestone provides excellent conditions for the nucleation and growth of C-S-H such that significantly more C-S-H nuclei are formed in presence of limestone compared to other SCM.
11. The reactivity of blast-furnace slag and fly ash depends on the particle size as well as on the intrinsic reactivity of especially the amorphous phases.
12. An increase in network modifying oxides (e.g. CaO) in the chemical composition of amorphous (calcium)aluminosilicates leads to an increasingly depolymerized network which in turn causes increasing reactivity. The role of amphoteric oxides (Al2O3, Fe2O3) that can be present as network modifying oxides as well as network forming oxides is not completely solved.
13. CO2-containing AFm-phases are thermodynamically more stable than monosulfoaluminate. This indirectly stabilizes the voluminous ettringite which causes a higher total volume of hydrates and lower porosity whereby higher compressive strength is reached.
14. Only a few percent of limestone in blended cement reacts chemically dependent on the Al2O3 available for reaction. Al2O3 that is provided by the reaction of Portland cement but also by the dissolution of SCM, especially by fly ash, reacts to form hemicarboaluminate which is transformed to monocarboaluminate as the hydration proceeds.
Methodology
15. The influence of SCM on the early hydration of Portland cement in binary (including blast-furnace slag or fly ash or limestone or quartz) and ternary (including fly ash and limestone) systems was investigated applying isothermal calorimetry and analysis of the pore solution chemistry. Calculated saturation indices and solubility products of relevant phases were correlated with heat development. Based on the gained data it was reviewed if mechanisms that control the hydration of pure phases in diluted systems are also prevailing in cement pastes under realistic conditions.
16. The influence of the chemical composition of synthetic glasses on their dissolution at high pH was investigated in highly diluted systems using ion chromatography. Pozzolanity tests were conducted on pastes using simplified model systems and glass-blended Portland cements. The process of the glass dissolution was investigated by isothermal calorimetry and by thermogravimetry. Correlation of experimentally determined total bound water with bound water determined by mass balance calculations as a function of amount of glass reacted allowed to estimate the degree of glass reaction in the pastes. Further on selective dissolution experiments were carried out to crosscheck the results of the bound water/mass balance approach.
17. The reaction kinetics of quaternary pastes containing blast-furnace slag and fly ash simultaneously in the presence of limestone were investigated up to 28 d using isothermal calorimetry and chemical shrinkage measurements. In addition strength tests on mortar bars were carried out.
18. Pastes of quaternary blends were also investigated in terms of hydrate assemblage at ages of up to 182 d. Thermodynamic calculations regarding total volume of hydrates as a function of limestone and fly ash/blast-furnace slag content were conducted. The calculations were supported by thermogravimetric determination of bound water and portlandite content as well as qualitative X-ray diffraction. The results were correlated with strength tests on mortar bars.
19. The pore solutions of hydrated quaternary blends were extracted and investigated by means of ion chromatography at ages of up to 728 d. Based on the ion concentrations in the solutions saturation indices were calculated for relevant phases. In order to gain better insight to the blast-furnace slag reaction sulphate speciation was carried out at two blast furnace slag levels (20 and 30 wt.%) for selected samples up to 91 d of hydration and at 91 d for the whole matrix under investigation.
Main results
20. Investigations on the early hydration kinetics of binary systems showed a higher heat flow in presence of SCM compared to neat Portland cement. This is caused by the higher surface area that is available for the nucleation of hydrates and by the lower (over)saturation with respect to C-S-H. An increase in the Ca-concentration in the pore solution did not cause lower dissolution rates of C3S as was reported for pure phases in diluted systems. The highest dissolution was observed in the presence of limestone, i.e. at the highest Ca-concentration. The general trend of the reaction rate is inversely related to the degree of undersaturation with respect to C3S. The more undersaturated the faster the observed reaction. The presence of increasing Al-concentrations caused a retardation of the reaction which is in line with investigations on pure phases in diluted systems. Higher sulphate concentrations could be detected for the fly ash containing blend which possibly hindered ettringite precipitation and results in higher Al-concentrations. Correspondingly the low sulphate concentrations lead to lower Al-concentrations in the presence of quartz, blast-furnace slag and limestone compared to fly ash.
21. The early hydration kinetics of quaternary systems is significantly accelerated in the presence of limestone while fly ash leads to retardation. Compared to reference systems containing inert quartz, investigations by means of isothermal calorimetry and chemical shrinkage revealed an acceleration caused by blast-furnace slag. Additions of fly ash, limestone or mixtures thereof introduced another acceleration but differences are too small to be significant and clear distinguishing between the various SCM is not possible.
22. Investigations on the reactivity of synthetic glasses showed that increasing amounts of network modifying oxides caused an increase in reactivity and dissolution rates. The results reveal that Al2O3 acts mainly as network modifying oxide in all investigated glasses. Experimentally determined bound water (thermogravimetric experiments) in model systems and blended cements can be compared with bound water determined by mass balance calculations carried out as a function of the amount of glass reacted. This enables to estimate the degree of glass reaction.
23. The actual content of blast-furnace slag, fly ash or limestone does not exert significant influence on the development of compressive strength up to 7 d. At later ages thermodynamic calculations predict a degree of CaCO3 reaction of 2 to 5 wt.%. This leads to the formation of hemicarbonate and monocarbonate whereby ettringite is indirectly stabilized. As a result the total amount of solids is increased and compressive strength shows a slight maximum. Hereby the degree of CaCO3 reaction depends on the Al2O3 available for reaction which is not only provided by the dissolution of Portland cement but especially by the dissolution of the fly ash.
24. In general the presence of blast-furnace slag and fly ash in the presence of limestone exerts little influence on the hydrate assemblage. The substitution of some of the blast-furnace slag by fly ash leads to a slight decrease of portlandite and C-S-H and gives rise to the formation of more monocarbonate and hemicarbonate. Portlandite is consumed in a pozzolanic reaction with the fly ash whereby C-S-H is formed. However, the low reactivity of the fly ash causes a decrease in the amount of C-S-H formed. Thereby a lower total volume of hydrates is generated which is in line with slightly lower compressive strength in case of increasing fly ash content. The overall influence is small and all systems investigated reach strength class 42.5 N according to EN 196-1.
25. Corresponding to the investigations of the hydrate assemblage the pore solution chemistry of quaternary systems showed only small variations. Depending on the fly ash content the highest variations are observed for aluminium, i.e. increasing fly ash content leads to higher Al-concentrations. Another effect of increasing fly ash contents is an increasing undersaturation with respect to portlandite and a decreasing undersaturation with respect to strätlingite indicating the dissolution of portlandite.
26. The total concentration of sulfur in the pore solution is controlled by sulphate (SO42–) while the concentrations of sulphite (SO32–) and thiosulphate (S2O32–) were very low. Up to 2 d of hydration about 90 % of the total sulphur is present as SO42–. After 91 d this value is reduced to about 36 % while about 28 % are present as S2O32–. In general higher blast-furnace slag content leads to higher concentrations of sulphite and thiosulphate after 7 d.
Bentonite sind quellfähige Tone, die häufig in der Umwelttechnik (in Abdichtungsbauwerken oder in der Bodensanierung) eingesetzt werden. Ziel der Arbeit war die Klärung, wie eine unterschiedliche Kationenbelegung mit Cu2+ und NH4+ die Eigenschaften der Bentonite bei Raumtemperatur und nach moderater Wärmebehandlung (300 – 450°C) beeinflusst. Im Blickpunkt stand insbesondere die gleichzeitige Präsenz von Kupfer- und Ammoniumionen, die als Vertreter für häufig auftretende Inhaltsstoffe von Wässern in der Umgebung von technischen Bentoniten ausgewählt wurden.
Die Untersuchungen zur Cu2+-Sorption bei Raumtemperatur und nach moderater Wärmebehandlung (300 – 450°C) erfolgten an Pulverproben von zwei technischen Bentoniten, die sich in der ursprünglichen Kationenbelegung, Art und Anteil an Nebengemengteilen, sowie der Schichtladungsverteilung in den Montmorilloniten unterscheiden. Vor der Wärmebehandlung wurden die Bentonite durch Kontakt mit verschieden konzentrierten Kupfer- und Ammoniumlösungen mit unterschiedlichen Gehalten der Kationen Cu2+, NH4+, Na+, Ca2+, Mg2+ belegt.
Der Eintrag von Kupferionen in die Bentonite durch Kationenaustausch bei Raumtemperatur wurde erwartungsgemäß durch präsente Nebengemengteile (wie Carbonat) beeinflusst, so dass die Kupferionen zusätzlich spezifisch adsorbiert und in festen Phasen angereichert wurden.
Die Cu2+-Fixierung infolge der Wärmebehandlung wurde vom Cu2+-Totalgehalt in den Bentoniten, der Präsenz von Nebengemengteilen und die Schichtladungsverteilung in den Montmorilloniten beeinflusst. Es waren generell Behandlungstemperaturen von > 400°C erforderlich, um Cu2+-Fixierungsraten von > 95% zu erzielen.
Waren während der Wärmebehandlung neben Cu2+-Ionen gleichzeitig NH4+-Ionen in den Bentoniten präsent, konnte die Cu2+-Fixierungstemperatur herabgesetzt werden. Die Deammonisierung (NH4+ --> NH3 + H+) der NH4+-belegten Bentonite fand größtenteils unterhalb der Dehydroxylierungstemperatur der Bentonite statt.
Durch Untersuchungen (XRD, FTIR, NMR, ESR) zum Mechanismus der Cu2+-Einbindung in die Bentonite an speziell aufbereiteten Proben (carbonatfrei, < 2 µm) konnte nachgewiesen werden, dass in den Cu2+-belegten Montmorilloniten die Cu2+-Ionen infolge der Wärmebehandlung nicht bis in die Oktaederschicht der Tonminerale vordringen, sondern nur bis in die Tetraederschicht wandern. In den NH4+-belegten Montmorilloniten treten im Zusammenhang mit der Deammonisierung keine zusätzlichen Strukturänderungen (wie Auflösung der Oktaederschicht) infolge der Wärmebehandlung auf.
In this thesis, a new approach is developed for applications of shape optimization on the time harmonic wave propagation (Helmholtz equation) for acoustic problems. This approach is introduced for different dimensional problems: 2D, 3D axi-symmetric and fully 3D problems. The boundary element method (BEM) is coupled with the isogeometric analysis (IGA) forming the so-called (IGABEM) which speeds up meshing and gives higher accuracy in comparison with standard BEM. BEM is superior for handling unbounded domains by modeling only the inner boundaries and avoiding the truncation error, present in the finite element method (FEM) since BEM solutions satisfy the Sommerfeld radiation condition automatically. Moreover, BEM reduces the space dimension by one from a volumetric three-dimensional problem to a surface two-dimensional problem, or from a surface two-dimensional problem to a perimeter one-dimensional problem. Non-uniform rational B-splines basis functions (NURBS) are used in an isogeometric setting to describe both the CAD geometries and the physical fields.
IGABEM is coupled with one of the gradient-free optimization methods, the Particle Swarm Optimization (PSO) for structural shape optimization problems. PSO is a straightforward method since it does not require any sensitivity analysis but it has some trade-offs with regard to the computational cost. Coupling IGA with optimization problems enables the NURBS basis functions to represent the three models: shape design, analysis and optimization models, by a definition of a set of control points to be the control variables and the optimization parameters as well which enables an easy transition between the three models.
Acoustic shape optimization for various frequencies in different mediums is performed with PSO and the results are compared with the benchmark solutions from the literature for different dimensional problems proving the efficiency of the proposed approach with the following remarks:
- In 2D problems, two BEM methods are used: the conventional isogeometric boundary element method (IGABEM) and the eXtended IGABEM (XIBEM) enriched with the partition-of-unity expansion using a set of plane waves, where the results are generally in good agreement with the linterature with some computation advantage to XIBEM which allows coarser meshes.
-In 3D axi-symmetric problems, the three-dimensional problem is simplified in BEM from a surface integral to a combination of two 1D integrals. The first is the line integral similar to a two-dimensional BEM problem. The second integral is performed over the angle of revolution. The discretization is applied only to the former integration. This leads to significant computational savings and, consequently, better treatment for higher frequencies over the full three-dimensional models.
- In fully 3D problems, a detailed comparison between two BEM methods: the conventional boundary integral equation (CBIE) and Burton-Miller (BM) is provided including the computational cost. The proposed models are enhanced with a modified collocation scheme with offsets to Greville abscissae to avoid placing collocation points at the corners. Placing collocation points on smooth surface enables accurate evaluation of normals for BM formulation in addition to straightforward prediction of jump-terms and avoids singularities in $\mathcal{O} (1/r)$ integrals eliminating the need for polar integration. Furthermore, no additional special treatment is required for the hyper-singular integral while collocating on highly distorted elements, such as those containing sphere poles. The obtained results indicate that, CBIE with PSO is a feasible alternative (except for a small number of fictitious frequencies) which is easier to implement. Furthermore, BM presents an outstanding treatment of the complicated geometry of mufflers with internal extended inlet/outlet tube as an interior 3D Helmholtz acoustic problem instead of using mixed or dual BEM.
A novel combination of the ant colony optimization algorithm (ACO)and computational fluid dynamics (CFD) data is proposed for modeling the multiphase chemical reactors. The proposed intelligent model presents a probabilistic computational strategy for predicting various levels of three-dimensional bubble column reactor (BCR) flow. The results prove an enhanced communication between ant colony prediction and CFD data in different sections of the BCR.
Die vorliegende Arbeit richtet sich an Ingenieur*innen und Wissenschaftler*innen der technischen Gebäudeausrüstung. Sie greift einen sich abzeichnenden Änderungsbedarf in der Umwelt- und Nachhaltigkeitsbewertung von Gebäuden und wärmetechnischen Anlagen auf. Der aktuell genutzte nicht erneuerbare Primärenergiebedarf wird insbesondere hinsichtlich künftiger politischer Klima- und Umweltschutzziele als alleinige Bewertungsgröße nicht ausreichend sein. Die mit dieser Arbeit vorgestellte Ökoeffizienzbewertungsmethode kann als geeignetes Instrument zur Lösung der Probleme beitragen. Sie ermöglicht systematische, ganzheitliche Bewertungen und reproduzierbare Vergleiche wärmetechnischer Anlagen bezüglich ihrer ökologischen und ökonomischen Nachhaltigkeit. Die wesentlichsten Neuentwicklungen sind die spezifische Umweltleistung, in Erweiterung zum genutzten Primärenergiefaktor, und der Ökoeffizienzindikator UWI.
The development of a hydro-mechanically coupled Coupled-Eulerian–Lagrangian (CEL) method and its application to the back-analysisof vibratory pile driving model tests in water-saturated sand is presented. The predicted pile penetration using this approachis in good agreement with the results of the model tests as well as with fully Lagrangian simulations. In terms of pore water pressure, however, the results of the CEL simulation show a slightly worse accordance with the model tests compared to the Lagrangian simulation. Some shortcomings of the hydro-mechanically coupled CEL method in case of frictional contact problems and pore fluids with high bulk modulus are discussed. Lastly, the CEL method is applied to the simulation of vibratory driving of open-profile piles under partially drained conditions to study installation-induced changes in the soil state. It is concluded that the proposed method is capable of realistically reproducing the most important mechanisms in the soil during the driving process despite its addressed shortcomings.
The world society faces a huge challenge to implement the human right of “access to sanitation”. More and more it is accepted that the conventional approach towards providing sanitation services is not suitable to solve this problem. This dissertation examines the possibility to enhance “access to sanitation” for people who are living in areas with underdeveloped water and wastewater infrastructure systems. The idea hereby is to follow an integrated approach for sanitation, which allows for a mutual completion of existing infrastructure with resource-based sanitation systems.
The notion “integrated sanitation system (iSaS)” is defined in this work and guiding principles for iSaS are formulated. Further on the implementation of iSaS is assessed at the example of a case study in the city of Darkhan in Mongolia. More than half of Mongolia’s population live in settlements where yurts (tents of Nomadic people) are predominant. In these settlements (or “ger areas”) sanitation systems are not existent and the hygienic situation is precarious.
An iSaS has been developed for the ger areas in Darkhan and tested over more than two years. Further on a software-based model has been developed with the goal to describe and assess different variations of the iSaS. The results of the assessment of material-flows, monetary-flows and communication-flows within the iSaS are presented in this dissertation. The iSaS model is adaptable and transferable to the socio-economic conditions in other regions and climate zones.
Die Gase Sauerstoff und Stickstoff werden für eine Vielzahl an technischen, industriellen, biologischen und medizinischen Einsatzzwecken benötigt. So liegen Anwendungsgebiete dieser Gase neben der klassischen metallverarbeitenden und der chemischen Industrie bei Sauerstoff vor allem in der Medizin, Verbrennungs- und Kläranlagenoptimierung sowie der Fischzucht und bei Stickstoff als Schutz- beziehungsweise Inertgas in der Kunststoffindustrie, der Luft- und Raumfahrt sowie dem Brandschutz.
Die Bereitstellung der Gase Sauerstoff und Stickstoff wird nahezu ausschließlich durch die Abtrennung aus der Umgebungsluft realisiert, welche aus ca. 78 Vol.-% Stickstoff, 21 Vol.-% Sauerstoff und 1 Vol.-% Spurengasen (Ar, CO2, Ne, He, ...) besteht. Am Markt etablierte Verfahren der Luftzerlegung sind das Linde-, das PSA- (pressure swing adsorption/Druckwechseladsorption) oder verschiedene Membran-Verfahren. Hierdurch werden die benötigten Gase entweder direkt vor Ort beim Verbraucher erzeugt (PSA- und Polymer-Membranverfahren: geringe Reinheiten) oder zentral in großen Anlagen hergestellt (Linde-Verfahren: hohe Reinheiten) und anschließend zum Verbraucher in Form von Flaschen- oder Tankgasen geliefert (Tansportkosten).
Für kleinere Verbraucher mit hohen Ansprüchen an die Reinheit des benötigten Sauerstoffs beziehungsweise Stickstoffs ergibt sich nur die Möglichkeit, die Gase als kostenintensive Transportgase zentraler Gaseversorger zu beziehen und sich somit in eine Abhängigkeit (Lieferverträge, Flaschen-/Tankmieten, ...) zu diesen zu begeben sowie eine eigene Lagerhaltung für die benötigten Gase (Mehraufwand, Lagerkosten, Platzbedarf) zu betreiben.
Ziel dieser Arbeit ist es, keramische Material-Systeme auf Basis chemischer Hochtemperatur-Reaktionen als Reaktive Oxidkeramiken zu entwickeln und diese hinsichtlich eines möglichen Einsatzes für die Sauerstoffseparation in neuartigen Luftzerlegungsanlagen zu untersuchen.
Derartige Anlagen sollen in ihrem Prinzip an die regenerative Sauerstoffseparation angelehnt sein und in ihren Reaktoren die Reaktiven Oxidkeramiken als Festbett-Material abwechselnd mit Luft be- und Vakuum oder O2-armen Atmosphären entladen.
Die Verwendung Reaktiver Oxidkeramiken, welche im Vergleich zu den bisherigen Materialien höhere Sauerstoffaustauschmengen und -raten bei gleichzeitig hoher Lebensdauer und Korrosionsbeständigkeit sowie relativ einfacher Handhabe aufweisen würden, soll ein Schritt in Richtung einer effizienten alternativen Luftzerlegungstechnologie sein.
Mit den Reaktiven Oxidkeramiken in einer Luftzerlegungsanlage sollte es im besten Fall möglich sein, in kleinen Anlagen sehr reinen Sauerstoff und zugleich sauerstofffreies Inertgas zu erzeugen sowie eine Sauerstoffan- oder -abreicherung von Luft, Prozess- oder Abgasen zu generieren.
Somit besäße eine solche, auf Reaktiven Oxidkeramiken basierende Technologie sehr weit gefächerte Einsatzgebiete und demzufolge ein enormes wirtschaftliches Potential.
Im Rahmen dieser Arbeit wird das Charakterisieren struktureller Veränderungen zementgebundener Baustoffe durch zwei auf dem Ultraschall-Transmissionsverfahren beruhenden Methoden der zerstörungsfreien Prüfung (ZfP) mit mechanischen Wellen vorgenommen.
Zur kontinuierlichen Charakterisierung der Erstarrung und Erhärtung frischer zementgebundener Systeme wird ein auf Ultraschallsensoren für Longitudinal- und Scherwellen basierendes Messsystem in Kombination mit zugehörigen Verfahrensweisen zur Datenauswertung konzipiert, charakterisiert und angewandt. Gegenüber der bislang üblichen alleinigen Bewertung der Verfestigung anhand indirekter Ultraschallparameter wie Ausbreitungsgeschwindigkeit, Signalenergie oder Frequenzgehalt der Longitudinalwelle lässt sich damit eine direkte, sensible Erfassung der sich während der Strukturbildung entwickelnden dynamischen elastischen Eigenschaften auf der Basis primärer physikalischer Werkstoffparameter erreichen. Insbesondere Scherwellen und der dynamische Schubmodul sind geeignet, den graduellen Übergang zum Festkörper mit Überschreiten der Perkolationsschwelle sensibel und unabhängig vom Luftgehalt zu erfassen. Die zeitliche Entwicklung der dynamischen elastischen Eigenschaften, die Strukturbildungsraten sowie die daraus extrahierten diskreten Ergebnisparameter ermöglichen eine vergleichende quantitative Charakterisierung der Strukturbildung zementgebundener Baustoffe aus mechanischer Sicht. Dabei lassen sich typische, oft unvermeidbare Unterschiede in der Zusammensetzung der Versuchsmischungen berücksichtigen.
Der Einsatz laserbasierter Methoden zur Anregung und Erfassung von mechanischen Wellen und deren Kombination zu Laser-Ultraschall zielt darauf ab, die mit der Anwendung des konventionellen Ultraschall-Transmissionsverfahrens verbundenen Nachteile zu eliminieren. Diese resultieren aus der Sensorgeometrie, der mechanischen Ankopplung und bei einer Vielzahl von Oberflächenpunkten aus einem hohen prüftechnischen Aufwand. Die laserbasierte, interferometrische Erfassung mechanischer Wellen ist gegenüber Ultraschallsensoren rauschbehaftet und vergleichsweise unsensibel. Als wesentliche Voraussetzung der scannenden Anwendung von Laser-Ultraschall auf zementgebundene Baustoffe erfolgen systematische experimentelle Untersuchungen zur laserinduzierten ablativen Anregung. Diese sollen zum Verständnis des Anregungsmechanismus unmittelbar auf den Oberflächen von zementgebundenen Baustoffen, Gesteinskörnungen und metallischen Werkstoffen beitragen, relevante Einflussfaktoren aus den charakteristischen Materialeigenschaften identifizieren, geeignete Prozessparameter gewinnen und die Verfahrensgrenzen aufzeigen. Unter Einsatz von Longitudinalwellen erfolgt die Anwendung von Laser-Ultraschall zur zeit- und ortsaufgelösten Charakterisierung der Strukturbildung und Homogenität frischer sowie erhärteter Proben zementgebundener Baustoffe. Während der Strukturbildung wird erstmals eine simultane berührungslose Erfassung von Longitudinal- und Scherwellen vorgenommen. Unter Anwendung von tomographischen Methoden (2D-Laufzeit¬tomo¬graphie) werden überlagerungsfreie Informationen zur räumlichen Verteilung struktureller Gefügeveränderungen anhand der longitudinalen Ausbreitungsgeschwindigkeit bzw. des relativen dynamischen Elastizitätsmoduls innerhalb von virtuellen Schnittebenen geschädigter Probekörper gewonnen. Als beton-schädigende Mechanismen werden exemplarisch der kombinierte Frost-Tausalz-Angriff sowie die Alkali-Kieselsäure-Reaktion (AKR) herangezogen.
Die im Rahmen dieser Arbeit entwickelten Verfahren der zerstörungsfreien Prüfung bieten erweiterte Möglichkeiten zur Charakterisierung zementgebundener Baustoffe und deren strukturellen Veränderungen und lassen sich zielgerichtet in der Werkstoffentwicklung, bei der Qualitätssicherung sowie zur Analyse von Schadensprozessen und -ursachen einsetzen.
Isogeometric analysis (IGA) is a numerical method for solving partial differential equations (PDEs), which was introduced with the aim of integrating finite element analysis with computer-aided design systems. The main idea of the method is to use the same spline basis functions which describe the geometry in CAD systems for the approximation of solution fields in the finite element method (FEM). Originally, NURBS which is a standard technology employed in CAD systems was adopted as basis functions in IGA but there were several variants of IGA using other technologies such as T-splines, PHT splines, and subdivision surfaces as basis functions. In general, IGA offers two key advantages over classical FEM: (i) by describing the CAD geometry exactly using smooth, high-order spline functions, the mesh generation process is simplified and the interoperability between CAD and FEM is improved, (ii) IGA can be viewed as a high-order finite element method which offers basis functions with high inter-element continuity and therefore can provide a primal variational formulation of high-order PDEs in a straightforward fashion. The main goal of this thesis is to further advance isogeometric analysis by exploiting these major advantages, namely precise geometric modeling and the use of smooth high-order splines as basis functions, and develop robust computational methods for problems with complex geometry and/or complex multi-physics.
As the first contribution of this thesis, we leverage the precise geometric modeling of isogeometric analysis and propose a new method for its coupling with meshfree discretizations. We exploit the strengths of both methods by using IGA to provide a smooth, geometrically-exact surface discretization of the problem domain boundary, while the Reproducing Kernel Particle Method (RKPM) discretization is used to provide the volumetric discretization of the domain interior. The coupling strategy is based upon the higher-order consistency or reproducing conditions that are directly imposed in the physical domain. The resulting coupled method enjoys several favorable features: (i) it preserves the geometric exactness of IGA, (ii) it circumvents the need for global volumetric parameterization of the problem domain, (iii) it achieves arbitrary-order approximation accuracy while preserving higher-order smoothness of the discretization. Several numerical examples are solved to show the optimal convergence properties of the coupled IGA–RKPM formulation, and to demonstrate its effectiveness in constructing volumetric discretizations for complex-geometry objects.
As for the next contribution, we exploit the use of smooth, high-order spline basis functions in IGA to solve high-order surface PDEs governing the morphological evolution of vesicles. These governing equations are often consisted of geometric PDEs, high-order PDEs on stationary or evolving surfaces, or a combination of them. We propose an isogeometric formulation for solving these PDEs. In the context of geometric PDEs, we consider phase-field approximations of mean curvature flow and Willmore flow problems and numerically study the convergence behavior of isogeometric analysis for these problems. As a model problem for high-order PDEs on stationary surfaces, we consider the Cahn–Hilliard equation on a sphere, where the surface is modeled using a phase-field approach. As for the high-order PDEs on evolving surfaces, a phase-field model of a deforming multi-component vesicle, which consists of two fourth-order nonlinear PDEs, is solved using the isogeometric analysis in a primal variational framework. Through several numerical examples in 2D, 3D and axisymmetric 3D settings, we show the robustness of IGA for solving the considered phase-field models.
Finally, we present a monolithic, implicit formulation based on isogeometric analysis and generalized-alpha time integration for simulating hydrodynamics of vesicles according to a phase-field model. Compared to earlier works, the number of equations of the phase-field model which need to be solved is reduced by leveraging high continuity of NURBS functions, and the algorithm is extended to 3D settings. We use residual-based variational multi-scale method (RBVMS) for solving Navier–Stokes equations, while the rest of PDEs in the phase-field model are treated using a standard Galerkin-based IGA. We introduce the resistive immersed surface (RIS) method into the formulation which can be employed for an implicit description of complex geometries using a diffuse-interface approach. The implementation highlights the robustness of the RBVMS method for Navier–Stokes equations of incompressible flows with non-trivial localized forcing terms including bending and tension forces of the vesicle. The potential of the phase-field model and isogeometric analysis for accurate simulation of a variety of fluid-vesicle interaction problems in 2D and 3D is demonstrated.
Tropical coral reefs, one of the world’s oldest ecosystems which support some of the highest levels of biodiversity on the planet, are currently facing an unprecedented ecological crisis during this massive human-activity-induced period of extinction. Hence, tropical reefs symbolically stand for the destructive effects of human activities on nature [4], [5]. Artificial reefs are excellent examples of how architectural design can be combined with ecosystem regeneration [6], [7], [8]. However, to work at the interface between the artificial and the complex and temporal nature of natural systems presents a challenge, i.a. in respect to the B-rep modelling legacy of computational modelling.
The presented doctorate investigates strategies on how to apply digital practice to realise what is an essential bulwark to retain reefs in impossibly challenging times. Beyond the main question of integrating computational modelling and high precision monitoring strategies in artificial coral reef design, this doctorate explores techniques, methods, and linking frameworks to support future research and practice in ecology led design contexts.
Considering the many existing approaches for artificial coral reefs design, one finds they often fall short in precisely understanding the relationships between architectural and ecological aspects (e.g. how a surface design and material composition can foster coral larvae settlement, or structural three-dimensionality enhance biodiversity) and lack an integrated underwater (UW) monitoring process. Such a process is necessary in order to gather knowledge about the ecosystem and make it available for design, and to learn whether artificial structures contribute to reef regeneration or rather harm the coral reef ecosystem.
For the research, empirical experimental methods were applied: Algorithmic coral reef design, high precision UW monitoring, computational modelling and simulation, and validated through parallel real-world physical experimentation – two Artificial Reef Prototypes (ARPs) in Gili Trawangan, Indonesia (2012–today). Multiple discrete methods and sub techniques were developed in seventeen computational experiments and applied in a way in which many are cross valid and integrated in an overall framework that is offered as a significant contribution to the field. Other main contributions include the Ecosystem-aware design approach, Key Performance Indicators (KPIs) for coral reef design, algorithmic design and fabrication of Biorock cathodes, new high precision UW monitoring strategies, long-term real-world constructed experiments, new digital analysis methods and two new front-end web-based tools for reef design and monitoring reefs. The methodological framework is a finding of the research that has many technical components that were tested and combined in this way for the very first time.
In summary, the thesis responds to the urgency and relevance in preserving marine species in tropical reefs during this massive extinction period by offering a differentiated approach towards artificial coral reefs – demonstrating the feasibility of digitally designing such ‘living architecture’ according to multiple context and performance parameters. It also provides an in-depth critical discussion of computational design and architecture in the context of ecosystem regeneration and Planetary Thinking. In that respect, the thesis functions as both theoretical and practical background for computational design, ecology and marine conservation – not only to foster the design of artificial coral reefs technically but also to provide essential criteria and techniques for conceiving them.
Keywords: Artificial coral reefs, computational modelling, high precision underwater monitoring, ecology in design.
Scalarization methods are a category of multiobjective optimization (MOO) methods. These methods allow the usage of conventional single objective optimization algorithms, as scalarization methods reformulate the MOO problem into a single objective optimization problem. The scalarization methods analysed within this thesis are the Weighted Sum (WS), the Epsilon-Constraint (EC), and the MinMax (MM) method. After explaining the approach of each method, the WS, EC and MM are applied, a-posteriori, to three different examples: to the Kursawe function; to the ten bar truss, a common benchmark problem in structural optimization; and to the metamodel of an aero engine exit module.
The aim is to evaluate and compare the performance of each scalarization method that is examined within this thesis. The evaluation is conducted using performance metrics, such as the hypervolume and the generational distance, as well as using visual comparison.
The application to the three examples gives insight into the advantages and disadvantages of each method, and provides further understanding of an adequate application of the methods concerning high dimensional optimization problems.
This paper presents the development of an assessment scheme for a visual qualitative evaluation of nailed connections in existing structures, such as board trusses. In terms of further use and preservation, a quick visual inspection will help to evaluate the quality of a structure regarding its load-bearing capacity and deformation behaviour. Tests of old and new nailed joints in combination with a rating scheme point out the correlation between the load-bearing capacity and condition of a joint. Old joints of comparatively good condition tend to exhibit better results than those of poor condition. Moreover, aged joints are generally more load-bearing than newly assembled ones.
Finite Element Simulations of dynamically excited structures are mainly influenced by the mass, stiffness, and damping properties of the system, as well as external loads. The prediction quality of dynamic simulations of vibration-sensitive components depends significantly on the use of appropriate damping models. Damping phenomena have a decisive influence on the vibration amplitude and the frequencies of the vibrating structure. However, developing realistic damping models is challenging due to the multiple sources that cause energy dissipation, such as material damping, different types of friction, or various interactions with the environment.
This thesis focuses on thermoelastic damping, which is the main cause of material damping in homogeneous materials. The effect is caused by temperature changes due to mechanical strains. In vibrating structures, temperature gradients arise in adjacent tension and compression areas. Depending on the vibration frequency, they result in heat flows, leading to increased entropy and the irreversible transformation of mechanical energy into thermal energy.
The central objective of this thesis is the development of efficient simulation methods to incorporate thermoelastic damping in finite element analyses based on modal superposition. The thermoelastic loss factor is derived from the structure's mechanical mode shapes and eigenfrequencies. In subsequent analyses that are performed in the time and frequency domain, it is applied as modal damping.
Two approaches are developed to determine the thermoelastic loss in thin-walled plate structures, as well as three-dimensional solid structures. The realistic representation of the dissipation effects is verified by comparing the simulation results with experimentally determined data. Therefore, an experimental setup is developed to measure material damping, excluding other sources of energy dissipation.
The three-dimensional solid approach is based on the determination of the generated entropy and therefore the generated heat per vibration cycle, which is a measure for thermoelastic loss in relation to the total strain energy. For thin plate structures, the amount of bending energy in a modal deformation is calculated and summarized in the so-called Modal Bending Factor (MBF). The highest amount of thermoelastic loss occurs in the state of pure bending. Therefore, the MBF enables a quantitative classification of the mode shapes concerning the thermoelastic damping potential.
The results of the developed simulations are in good agreement with the experimental results and are appropriate to predict thermoelastic loss factors. Both approaches are based on modal superposition with the advantage of a high computational efficiency. Overall, the modeling of thermoelastic damping represents an important component in a comprehensive damping model, which is necessary to perform realistic simulations of vibration processes.
In this study, we propose a nonlocal operator method (NOM) for the dynamic analysis of (thin) Kirchhoff plates. The nonlocal Hessian operator is derived based on a second-order Taylor series expansion. The NOM does not require any shape functions and associated derivatives as ’classical’ approaches such as FEM, drastically facilitating the implementation. Furthermore, NOM is higher order continuous, which is exploited for thin plate analysis that requires C1 continuity. The nonlocal dynamic governing formulation and operator energy functional for Kirchhoff plates are derived from a variational principle. The Verlet-velocity algorithm is used for the time discretization. After confirming the accuracy of the nonlocal Hessian operator, several numerical examples are simulated by the nonlocal dynamic Kirchhoff plate formulation.
Material failure can be tackled by so-called nonlocal models, which introduce an intrinsic length scale into the formulation and, in the case of material failure, restore the well-posedness of the underlying boundary value problem or initial boundary value problem. Among nonlocal models, peridynamics (PD) has attracted a lot of attention as it allows the natural transition from continuum to discontinue and thus allows modeling of discrete cracks without the need to describe and track the crack topology, which has been a major obstacle in traditional discrete crack approaches. This is achieved by replacing the divergence of the Cauchy stress tensor through an integral over so-called bond forces, which account for the interaction of particles. A quasi-continuum approach is then used to calibrate the material parameters of the bond forces, i.e., equating the PD energy with the energy of a continuum. One major issue for the application of PD to general complex problems is that they are limited to fairly simple material behavior and pure mechanical problems based on explicit time integration. PD has been extended to other applications but losing simultaneously its simplicity and ease in modeling material failure. Furthermore, conventional PD suffers from instability and hourglass modes that require stabilization. It also requires the use of constant horizon sizes, which drastically reduces its computational efficiency. The latter issue was resolved by the so-called dual-horizon peridynamics (DH-PD) formulation and the introduction of the duality of horizons.
Within the nonlocal operator method (NOM), the concept of nonlocality is further extended and can be considered a generalization of DH-PD. Combined with the energy functionals of various physical models, the nonlocal forms based on the dual-support concept can be derived. In addition, the variation of the energy functional allows implicit formulations of the nonlocal theory. While traditional integral equations are formulated in an integral domain, the dual-support approaches are based on dual integral domains. One prominent feature of NOM is its compatibility with variational and weighted residual methods. The NOM yields a direct numerical implementation based on the weighted residual method for many physical problems without the need for shape functions. Only the definition of the energy or boundary value problem is needed to drastically facilitate the implementation. The nonlocal operator plays an equivalent role to the derivatives of the shape functions in meshless methods and finite element methods (FEM). Based on the variational principle, the residual and the tangent stiffness matrix can be obtained with ease by a series of matrix multiplications. In addition, NOM can be used to derive many nonlocal models in strong form.
The principal contributions of this dissertation are the implementation and application of NOM, and also the development of approaches for dealing with fractures within the NOM, mostly for dynamic fractures. The primary coverage and results of the dissertation are as follows:
-The first/higher-order implicit NOM and explicit NOM, including a detailed description of the implementation, are presented. The NOM is based on so-called support, dual-support, nonlocal operators, and an operate energy functional ensuring stability. The nonlocal operator is a generalization of the conventional differential operators. Combining with the method of weighted residuals and variational principles, NOM establishes the residual and tangent stiffness matrix of operate energy functional through some simple matrix without the need of shape functions as in other classical computational methods such as FEM. NOM only requires the definition of the energy drastically simplifying its implementation. For the sake of conciseness, the implementation in this chapter is focused on linear elastic solids only, though the NOM can handle more complex nonlinear problems. An explicit nonlocal operator method for the dynamic analysis of elasticity solid problems is also presented. The explicit NOM avoids the calculation of the tangent stiffness matrix as in the implicit NOM model. The explicit scheme comprises the Verlet-velocity algorithm. The NOM can be very flexible and efficient for solving partial differential equations (PDEs). It's also quite easy for readers to use the NOM and extend it to solve other complicated physical phenomena described by one or a set of PDEs. Several numerical examples are presented to show the capabilities of this method.
-A nonlocal operator method for the dynamic analysis of (thin) Kirchhoff plates is proposed. The nonlocal Hessian operator is derived from a second-order Taylor series expansion. NOM is higher-order continuous, which is exploited for thin plate analysis that requires $C^1$ continuity. The nonlocal dynamic governing formulation and operator energy functional for Kirchhoff plates are derived from a variational principle. The Verlet-velocity algorithm is used for time discretization. After confirming the accuracy of the nonlocal Hessian operator, several numerical examples are simulated by the nonlocal dynamic Kirchhoff plate formulation.
-A nonlocal fracture modeling is developed and applied to the simulation of quasi-static and dynamic fractures using the NOM. The phase field's nonlocal weak and associated strong forms are derived from a variational principle. The NOM requires only the definition of energy. We present both a nonlocal implicit phase field model and a nonlocal explicit phase field model for fracture; the first approach is better suited for quasi-static fracture problems, while the key application of the latter one is dynamic fracture. To demonstrate the performance of the underlying approach, several benchmark examples for quasi-static and dynamic fracture are solved.
The Variability of the Void Ratio of Sand and its Effect on Settlement and Infinite Slope Stability
(2018)
The uncertainty of a soil property can significantly affect the physical behavior of soil, so as to influence geotechnical practice. The uncertainty can be expressed by its stochastic parameters, including the mean, the standard deviation, and the spatial correlation length. These stochastic parameters are regarded as constant value in most of the former studies. The main aim of this thesis is to prove whether they are depth-dependent, and to evaluate the effect of this depth-dependent character on both the settlement and the infinite slope stability during rainwater infiltration.
A stochastic one-dimensional settlement simulation is carried out using random finite element method with the von Wolffersdorff hypoplastic model, so as to evaluate the effect of stress level on the stochastic parameters of void ratio related parameters of sand. It is found that these stochastic parameters are both stress-dependent and depth-dependent.
The non-stationary random field, considering the depth-dependent character of these stochastic parameters, can be generated through the distortion of the stationary random field.
The one-dimensional settlement analysis is carried out to evaluation the effect of the depth-dependent character of the stochastic parameters of void ratio on the strain. It is found that the depth-dependent character has low effect on the strain.
The deterministic analysis of infinite slope stability during rainwater infiltration is simulated.
The transient seepage is carried out using finite difference method, while the steady state seepage is simulated using the analytical solution. The saturated hydraulic conductivity (ks) is taken as the only variable. The results show that the depth-dependent ks has a significant influence on the stability of the slope when the negative flux is high. Without considering the depth-dependent character, can overestimate the factor of safety of the slope. A slope can fail if the depth-dependent character is considered, while it is stable if the depth-dependent character is neglected. The failure time of the slope with a greater depth-dependent ks is earlier during transient infiltration.
Meanwhile, the stochastic infinite slope stability analysis during infiltration, is also carried out to highlight the effect of the depth-dependent character of the stochastic parameters of ks. The results show that: the probability of failure is significantly increased if the depth-dependent character of mean is considered, while, it is moderately reduced if the depth-dependent character of the standard deviation is accounted. If the depth-dependent character of both the mean and standard deviation of ks is considered, the depth-dependent mean value plays a dominant influence on the results. Furthermore, the depth-dependent character of the spatial correlation length can slightly reduce the probability of failure.