05 Fakultät Informatik, Elektrotechnik und Informationstechnik

Permanent URI for this collectionhttps://elib.uni-stuttgart.de/handle/11682/6

Browse

Filters

2013 - 20183

Settings

Institute: search.filters.UBSInstitut.Institut für Parallele und Verteilte SystemeStart date: 2010End date: 2019Author: search.filters.author.Scheufele, Klaudius

Search Results

Now showing 1 - 3 of 3
  • Thumbnail Image
    ItemOpen Access
    Robust Quasi-Newton methods for partitioned fluid-structure simulations
    (2015) Scheufele, Klaudius
    In recent years, quasi-Newton schemes have proven to be a robust and efficient way for the coupling of partitioned multi-physics simulations in particular for fluid-structure interaction. The focus of this work is put on the coupling of partitioned fluid-structure interaction, where minimal interface requirements are assumed for the respective field solvers, thus treated as black box solvers. The coupling is done through communication of boundary values between the solvers. In this thesis a new quasi-Newton variant (IQN-IMVJ) based on a multi-vector update is investigated in combination with serial and parallel coupling systems. Due to implicit incorporation of passed information within the Jacobian update it renders the problem dependent parameter of retained previous time steps unnecessary. Besides, a whole range of coupling schemes are categorized and compared comprehensively with respect to robustness, convergence behaviour and complexity. Those coupling algorithms differ in the structure of the coupling, i.\,e., serial or parallel execution of the field solvers and the used quasi-Newton methods. A more in-depth analysis for a choice of coupling schemes is conducted for a set of strongly coupled FSI benchmark problems, using the in-house coupling library preCICE. The superior convergence behaviour and robust nature of the IQN-IMVJ method compared to well known state of the art methods such as the IQN-ILS method, is demonstrated here. It is confirmed that the multi-vector method works optimal without the need of tuning problem dependent parameters in advance. Furthermore, it appears to be especially suitable in conjunction with the parallel coupling system, in that it yields fairly similar results for parallel and serial coupling. Although we focus on FSI simulation, the considered coupling schemes are supposed to be equally applicable to various kinds of different volume- or surface-coupled problems.
  • Thumbnail Image
    ItemOpen Access
    Robuste Multilevel-Lösung elliptischer partieller Differentialgleichungen mit springenden Koeffizienten
    (2013) Scheufele, Klaudius
    In dieser Arbeit wird eine robuste Multilevel-Lösung für elliptische partielle Differentialgleichungen mit springenden Koeffizientenfunktionen im Kontext der Partition of Unity Methode realisiert und analysiert. Bei dieser gitterfreien Methode müssen die Koeffizientensprünge nicht auf dem gröbsten Level geometrisch aufgelöst sein, vielmehr kann durch geeignete Anreicherungsfunktionen mittels algebraischer Verfeinerung die Approximationsqualität verbessert werden. Die Implementierung eines stabilen und robusten Multilevel-Lösers sowie die Realisierung verschiedener Anreicherungsfunktionen sind Kernbereich dieser Arbeit. Insbesondere werden verschiedene Anreicherungsfunktionen entwickelt und deren Auswirkung im Hinblick auf die Robustheit des Lösers untersucht. Mögliche Ursachen für nicht robustes Verhalten des Lösers werden in diesem Zusammenhang detailliert diskutiert und Ansätze für Verbesserungen gegeben. Der realisierte Multilevel-Löser zeigt im Vergleich zu früheren Arbeiten effizienteres und für viele Fälle durchaus robustes Verhalten. Einfache Distanzfunktionen führen dabei i. A. zu den besten Ergebnissen jedoch lässt sich durch alleinige Verbesserung der Approximationsqualität der lokalen Anreicherungsräume die Effizienz und Robustheit des Lösers aufgrund schlechterer Glättungseigenschaften nicht beliebig steigern.
  • Thumbnail Image
    ItemOpen Access
    Coupling schemes and inexact Newton for multi-physics and coupled optimization problems
    (2018) Scheufele, Klaudius; Mehl, Miriam (Prof. Dr.)
    This work targets mathematical solutions and software for complex numerical simulation and optimization problems. Characteristics are the combination of different models and software modules and the need for massively parallel execution on supercomputers. We consider two different types of multi-component problems in Part I and Part II of the thesis: (i) Surface coupled fluid- structure interactions and (ii) analysis of medical MR imaging data of brain tumor patients. In (i), we establish highly accurate simulations by combining different aspects such as fluid flow and arterial wall deformation in hemodynamics simulations or fluid flow, heat transfer and mechanical stresses in cooling systems. For (ii), we focus on (a) facilitating the transfer of information such as functional brain regions from a statistical healthy atlas brain to the individual patient brain (which is topologically different due to the tumor), and (b) to allow for patient specific tumor progression simulations based on the estimation of biophysical parameters via inverse tumor growth simulation (given a single snapshot in time, only). Applications and specific characteristics of both problems are very distinct, yet both are hallmarked by strong inter-component relations and result in formidable, very large, coupled systems of partial differential equations. Part I targets robust and efficient quasi-Newton methods for black-box surface-coupling of parti- tioned fluid-structure interaction simulations. The partitioned approach allows for great flexibility and exchangeable of sub-components. However, breaking up multi-physics into single components requires advanced coupling strategies to ensure correct inter-component relations and effectively tackle instabilities. Due to the black-box paradigm, solver internals are hidden and information exchange is reduced to input/output relations. We develop advanced quasi-Newton methods that effectively establish the equation coupling of two (or more) solvers based on solving a non-linear fixed-point equation at the interface. Established state of the art methods fall short by either requiring costly tuning of problem dependent parameters, or becoming infeasible for large scale problems. In developing parameter-free, linear-complexity alternatives, we lift the robustness and parallel scalability of quasi-Newton methods for partitioned surface-coupled multi-physics simulations to a new level. The developed methods are implemented in the parallel, general purpose coupling tool preCICE. Part II targets MR image analysis of glioblastoma multiforme pathologies and patient specific simulation of brain tumor progression. We apply a joint medical image registration and biophysical inversion strategy, targeting at facilitating diagnosis, aiding and supporting surgical planning, and improving the efficacy of brain tumor therapy. We propose two problem formulations and decompose the resulting large-scale, highly non-linear and non-convex PDE-constrained optimization problem into two tightly coupled problems: inverse tumor simulation and medical image registration. We deduce a novel, modular Picard iteration-type solution strategy. We are the first to successfully solve the inverse tumor-growth problem based on a single patient snapshot with a gradient-based approach. We present the joint inversion framework SIBIA, which scales to very high image resolutions and parallel execution on tens of thousands of cores. We apply our methodology to synthetic and actual clinical data sets and achieve excellent normal-to-abnormal registration quality and present a proof of concept for a very promising strategy to obtain clinically relevant biophysical information. Advanced inexact-Newton methods are an essential tool for both parts. We connect the two parts by pointing out commonalities and differences of variants used in the two communities in unified notation.