Paralleler Rechencluster
Zusammenfassung der Projektergebnisse
Effiziente Kopplungsalgorithmen für Fluid-Struktur Wechselwirkung (FSI): Im Bereich der Methodenentwicklung von Simulationsmodellen für FSI-Phänomene wurden vor allem folgende neuartige Ansätze entwickelt. a) Eine monolithische Formulierung des gekoppelten nichtlinearen 3-Feld-Problems. Ein spezielles Augenmerk lag u.a. auf der Konsistenz bei Verwendung unterschiedlicher Zeitintegratoren in den Einzelfeldern sowie bei der Zeitadaptivität. b) Effiziente iterative und parallele Gleichungslösungsverfahren für monolithische FSI-Probleme auf der Basis von algebraischen Mehrgittermethoden. Es wurden Verfahren entwickelt, die die Interaktionsgleichungen auch auf den algebraischen Grobleveln enthalten. c) Bei FSI Problemen mit Topologieänderungen sind ALE Ansätze nur bedingt brauchbar. Es wurden neuartige sog. Fixed-Grid Ansätze entwickelt, bei denen die Interface Gleichungen über einen XFEM Ansatz bei nichtkonformer Diskretisierung eingebracht werden. Zusätzlich ist es möglich über eine gemischte ALE/Fixed-Grid Diskretisierung bewegte Randschichten im Fluid gut aufzulösen und diese mit der XFEM Methodik im Rahmen einer nichtkonformen Fluid-Fluid Kopplung in ein Fixed-Grid einzubetten. Ein mehrskaliges virtuelles Lungenmodell: Es wurde ein virtuelles Lungenmodell entwickelt, das die relevanten mechanischen Effekte (Strömung, Transport, Gewebedehnung, etc.) auf den unterschiedlichen Skalen und biologischen Prozesse berücksichtigt und koppelt und sich zur Prädiktion von Therapien und Krankheitsverläufen eignet. Das Lungenmodell ist mit einer geschachtelt parallelen Software realisiert, die eine effiziente Betrachtung dieses hochkomplexen Systems überhaupt erst ermöglicht. Bei dem Modell handelt es sich um das weltweit derzeit umfassendste Computermodell der Lunge. Kontaktmechanik bei großen Deformationen: Neuartige Ansätze für reibungsbehaftete Kontaktmechanik auf der Basis von dualen Mortaransätzen sowie deren konsistente Linearisierung wurden entwickelt. Dabei wurden die auftretenden Ungleichheitsnebenbedingungen über nichtlineare Komplementärfunktionen (NCP) eingebracht sowie die Bestimmung des sog. Aktiven Set über eine semi-smooth Newton Iteration eingebettet. Parallele Aspekte der Lastverteilung und der Kontaktsuche wurden studiert und hierzu effiziente Algorithmen entwickelt. Rupturrisiko und Mechanobiologie des Abdominellen Aortenaneurysma: In Kooperation mit Chirurgen, Humanbiologen und med. Bildgebungsexperten wurden Simulationsmodelle zur Risikoabschätzung und Vorhersage des Krankheitsverlaufs beim AAA entwickelt. Dabei wurden Aspekte der Mechanobiologie über die Zusammenführung von simulationsbasierten mechanischen Größen und klinisch messbaren sowie experimentellen Informationen erfasst. Zur Quantifizierung von Unsicherheiten (UQ) der Modellbildung wurden parallele stochastische Finite Elemente Methoden entwickelt, die auf der Basis einer Zweiskalen Diskretisierung und eines Bayesischen Regressionsmodells quasi erstmals eine vollwertige UQ auch bei diesen großen, nichtlinearen Modellen erlaubt. Brownsche Dynamik biopolymerer Netzwerke: Biopolymere Netzwerke, wie bspw. das Zytoskelett in biologische Zellen, haben eine große Bedeutung in der Biophysik. Um die Brown’sche Bewegung solcher Netzwerke zu simulieren, wurde ein neuartiges Modell und eine neuartige Methode auf Basis eindimensionaler Kontinua entwickelt und für Höchstleistungsrechner implementier. Dieser Ansatz ist konkurrierenden Verfahren derzeit um einen Faktor von drei bis vier Zehnerpotenzen überlegen und erlaubt so die Beantwortung bislang unbeantworteter Fragestellungen. Galvanische Prozesse und Elektrochemie: Die Galvanotechnik ist ein weitverbreitetes industrielles Verfahren zur funktionalen und dekorativen Ummantelung elektrisch leitfähiger Bauteile mit metallischen Schichten. Um galvanische Prozesse numerisch simulieren zu können, wurde ein numerisches Mehrskalenmodell entwickelt, das auf stabilisierten finiten Elementen beruht. Sämtliche strömungsmechanischen, elektrischen und chemischen Feldgrößen werden in vollständig gekoppelter Form berechnet und es wurden sowohl natürliche Konvektion infolge von Dichtevariationen in Elektrolyten als auch nichtlineare kinetische Gesetze für die elektrochemischen Oberflächenreaktionen zwischen dem Elektrolyten und den Elektroden berücksichtigt. Turbulente Ein- und Mehrphasenströmungen: Turbulenzmodelle für Large-Eddy Simulationen für ein- und mehrphasige Strömungen einschließlich Skalartransport wurden entwickelt. Für den neuartigen Subskalenansatz basierend auf den multifraktalen Eigenschaften dieser Phänomene, erfolgte die erforderliche Skalenseparation unter Verwendung eines auf rein algebraischen Mehrgitteroperatoren basierendem Ansatz (AVM4). Bei den mehrphasigen Strömungen wurde eine XFEM Diskretisierung zur Beschreibung der nichtkonformen Zweiphasenkopplung entwickelt und eingesetzt.
Projektbezogene Publikationen (Auswahl)
-
An algebraic variational multiscale-multigrid method for large-eddy simulation: generalized-alpha time integration, Fourier analysis and Application to turbulent flow past a square-section cylinder.
Computational Mechanics, Vol. 47. 2011, Issue 2, pp 217–233.
Gravemeier V., Kronbichler M., Gee M.W., Wall W.A.
-
Local strain distribution in real three-dimensional alveolar
geometries. Annals of Biomedical Engineering, Vol. 39. 2011, pp. 2835-2843.
Rausch S.M.K., Haberthür D., Stampanoni M., Schittny J.C., Wall W.A.
-
Dual quadratic mortar finite element methods for 3D finite deformation contact. SIAM Journal on Scientific Computing, Vol. 34. 2012, Issue 4, B421-B446.
Popp A., Wohlmuth B.I., Gee M.W., Wall W.A.
-
Information-flux method: a meshfree maximum-entropy Petrov-Galerkin method including stabilised finite element methods. Computer Methods in Applied Mechanics and Engineering, Vol. 241–244. 2012, pp. 225-237.
Nissen K., Cyron C.J., Gravemeier V., Wall W.A.
-
A monolithic computational approach to thermo-structure interaction.
International Journal for Numerical Methods in Engineering, Vol. 95. 2013, Issue 13, pp. 1053–1078.
Danowski C., Gravemeier V., Yoshihara L., Wall W.A.
-
Measuring and modeling patient-specific distributions of material properties in abdominal aortic aneurysm wall. Biomechanics and Modeling in Mechanobiology, Vol. 12. 2013, Issue 4, pp 717–733.
Reeps C., Maier A., Pelisek J., Härtl F., Grabher-Meier V., Wall W.A., Essler M., Eckstein H.-H., Gee M.W.
-
A combined fluid-structure interaction and multi-field scalar
transport model for simulating mass transport in biomechanics.
International Journal for Numerical Methods in Engineering,
Vol. 100. 2014, Issue 4, pp. 277–299.
Yoshihara L., Coroneo M., Comerford A., Bauer G., Klöppel T., Wall W.A.
-
An Objective 3D Large Deformation Finite Element Formulation for Geometrically Exact Curved Kirchhoff Rods. Computer Methods in Applied Mechanics and Engineering, Vol. 278. 2014, pp. 445-478.
Meier C., Popp A., Wall W.A.
-
Numerical identification method for the non-linear viscoelastic compressible behaviour of soft tissue using uniaxial tensile tests and Image registration - Application to rat lung parenchyma. Journal of the Mechanical Behavior of Biomedical Materials,
Vol. 29. 2014, pp. 360-374.
Bel-Brunon A., Kehl S., Martin C., Uhlig S., Wall W.A.
-
Towards efficient uncertainty quantification in complex and large scale biomechanical problems based on a Bayesian multi fidelity scheme. Biomechanics and Modeling in Mechanobiology, accepted 2014
Biomechanics and Modeling in Mechanobiology, Vol. 14. 2015, Issue 3, pp 489–513.
Biehler J., Gee M.W., Wall W.A.