Multilevel Domain Decomposition Algorithms for Monolithic Fluid-Structure Interaction Problems with Application to Haemodynamics

Finite element techniques for solving the problem of fluid-structure interaction of an elastic solid material in a laminar incompressible viscous flow are described. The mathematical problem consists of the Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian formulation coupled with a non-linear structure model, considering the problem as one continuum. The coupling between the structure and the fluid is enforced inside a monolithic framework which computes simultaneously for the fluid and the structure unknowns within a unique solver. We used the well-known Crouzeix-Raviart finite element pair for discretization in space and the method of lines for discretization in time. A stability result using the Backward-Euler time-stepping scheme for both fluid and solid part and the finite element method for the space discretization has been proved. The resulting linear system has been solved by multilevel domain decomposition techniques. Our strategy is to solve several local subproblems over subdomain patches using the Schur-complement or GMRES smoother within a multigrid iterative solver. For validation and evaluation of the accuracy of the proposed methodology, we present corresponding results for a set of two FSI benchmark configurations which describe the self-induced elastic deformation of a beam attached to a cylinder in a laminar channel flow, allowing stationary as well as periodically oscillating deformations, and for a benchmark proposed by COMSOL multiphysics where a narrow vertical structure attached to the bottom wall of a channel bends under the force due to both viscous drag and pressure. Then, as an example of fluid-structure interaction in biomedical problems, we considered the academic numerical test which consists in simulating the pressure wave propagation through a straight compliant vessel. All the tests show the applicability and the numerical efficiency of our approach to both two-dimensional and three-dimensional problems.

Decomposition Methods and Network Design Problems

Decomposition based approaches are recalled from primal and dual point of view. The possibility of building partially disaggregated reduced master problems is investigated. This extends the idea of aggregated-versus-disaggregated formulation to a gradual choice of alternative level of aggregation. Partial aggregation is applied to the linear multicommodity minimum cost flow problem. The possibility of having only partially aggregated bundles opens a wide range of alternatives with different trade-offs between the number of iterations and the required computation for solving it. This trade-off is explored for several sets of instances and the results are compared with the ones obtained by directly solving the natural node-arc formulation. An iterative solution process to the route assignment problem is proposed, based on the well-known Frank Wolfe algorithm. In order to provide a first feasible solution to the Frank Wolfe algorithm, a linear multicommodity min-cost flow problem is solved to optimality by using the decomposition techniques mentioned above. Solutions of this problem are useful for network orientation and design, especially in relation with public transportation systems as the Personal Rapid Transit. A single-commodity robust network design problem is addressed. In this, an undirected graph with edge costs is given together with a discrete set of balance matrices, representing different supply/demand scenarios. The goal is to determine the minimum cost installation of capacities on the edges such that the flow exchange is feasible for every scenario. A set of new instances that are computationally hard for the natural flow formulation are solved by means of a new heuristic algorithm. Finally, an efficient decomposition-based heuristic approach for a large scale stochastic unit commitment problem is presented. The addressed real-world stochastic problem employs at its core a deterministic unit commitment planning model developed by the California Independent System Operator (ISO).

Extracting evolutionary information from the spectral decomposition of early-type galaxies.

Holding the major share of stellar mass in galaxies and being also old and passively evolving, early-type galaxies (ETGs) are the primary probes in investigating these various evolution scenarios, as well as being useful means to provide insights on cosmological parameters. In this thesis work I focused specifically on ETGs and on their capability in constraining galaxy formation and evolution; in particular, the principal aims were to derive some of the ETGs evolutionary parameters, such as age, metallicity and star formation history (SFH) and to study their age-redshift and mass-age relations. In order to infer galaxy physical parameters, I used the public code STARLIGHT: this program provides a best fit to the observed spectrum from a combination of many theoretical models defined in user-made libraries. the comparison between the output and input light-weighted ages shows a good agreement starting from SNRs of ∼ 10, with a bias of ∼ 2.2% and a dispersion 3%. Furthermore, also metallicities and SFHs are well reproduced. In the second part of the thesis I performed an analysis on real data, starting from Sloan Digital Sky Survey (SDSS) spectra. I found that galaxies get older with cosmic time and with increasing mass (for a fixed redshift bin); absolute light-weighted ages, instead, result independent from the fitting parameters or the synthetic models used. Metallicities, instead, are very similar from each other and clearly consistent with the ones derived from the Lick indices. The predicted SFH indicates the presence of a double burst of star formation. Velocity dispersions and extinctiona are also well constrained, following the expected behaviours. As a further step, I also fitted single SDSS spectra (with SNR∼ 20), to verify that stacked spectra gave the same results without introducing any bias: this is an important check, if one wants to apply the method at higher z, where stacked spectra are necessary to increase the SNR. Our upcoming aim is to adopt this approach also on galaxy spectra obtained from higher redshift Surveys, such as BOSS (z ∼ 0.5), zCOSMOS (z 1), K20 (z ∼ 1), GMASS (z ∼ 1.5) and, eventually, Euclid (z 2). Indeed, I am currently carrying on a preliminary study to estabilish the applicability of the method to lower resolution, as well as higher redshift (z 2) spectra, just like the Euclid ones.

Bio-inspired polymer-supported nitrido molybdenum(VI) complexes for ammonia synthesis − reactivity towards protons and electrons : EPR investigations on paramagnetic molybdenum(V) and copper(II) complexes

Die biologische Stickstofffixierung durch Molybdän-haltige Nitrogenasen sowie die Erforschung des zugrundeliegenden komplexen Mechanismus (N2-Aktivierung an Metall-Zentren, 6-fache Protonierung und Reduktion, N–N Bindungsspaltung unter Bildung von Ammoniak) ist von erheblichem Interesse. Insbesondere Molybdän-Komplexe wurden bereits erfolgreich als Modellverbindungen für die Untersuchung elementarer Einzelschritte der N2-Aktivierung eingesetzt. Durch die Verwendung von Triamidoamin-Liganden ist es Schrock et al. sogar gelungen mehrere Katalysezyklen zu durchlaufen und einen Mechanismus zu formulieren. Trotz der sterisch anspruchsvollen Substituenten in den Schrock-Komplexen ist die Umsatzrate dieses homogenen Katalysators, aufgrund Komplex-Deaktivierung infolge intermolekularer Reaktionen wie Dimerisierung und Disproportionierung, limitiert. In der vorliegenden Arbeit wurden einige dieser Herausforderungen angegangen und die aktiven Spezies auf einer Festphase immobilisiert, um intermolekulare Reaktionen durch räumliche Isolierung der Komplexe zu unterdrücken.rnEin Polymer-verankertes Analogon des Schrock Nitrido-Molybdän(VI)-Komplexes wurde auf einem neuen Reaktionsweg synthetisiert. Dieser beinhaltet nur einen einzigen Reaktionsschritt, um die funktionelle Gruppe „MoN“ einzuführen. Protonierung des immobilisierten Nitrido-Molybdän(VI)-Komplexes LMoVIN (L = Polymer-verankerter Triamidoamin-Ligand) mit 2,6-Lutidinium liefert den entsprechenden Imido-Molybdän(VI)-Komplex. Durch anschließende Ein-Elektronen-Reduktion mit Cobaltocen wird der Polymer-angebundene Imido-Molybdän(V)-Komplex erhalten, bewiesen durch EPR-Spektroskopie (g1,2,3 = 1.989, 1.929, 1.902). Durch die Immobilisierung und die effektive räumliche Separation der Reaktionszentren auf der Festphase werden bimolekulare Nebenreaktionen, die oft in homogenen Systemen auftreten, unterdrückt. Dies ermöglicht zum ersten Mal die Darstellung des Imido-Molybdän(V)-Intermediates des Schrock-Zyklus.rnEPR-Spektren des als Spin-Label eingeführten immobilisierten Nitrato-Kupfer(II)-Komplexes wurden unter verschiedenen Bedingungen (Lösungsmittel, Temperatur) aufgenommen, wobei sich eine starke Abhängigkeit zwischen der Zugänglichkeit und Reaktivität der immobilisierten Reaktionszentren und der Art des Lösungsmittels zeigte. Somit wurde die Reaktivität von LMoVIN gegenüber Protonen und Elektronen, welches zur Bildung von NH3 führt, unter Verwendung verschiedener Lösungsmittel untersucht und optimiert. Innerhalb des kugelförmigen Polymers verläuft die Protonierung und Reduktion von LMoVIN stufenweise. Aktive Zentren, die sich an der „äußeren Schale“ des Polymers befinden, sind gut zugänglich und reagieren schnell nach H+/e− Zugabe. Aktive Zentren im „Inneren des Polymers“ hingegen sind schlechter zugänglich und zeigen langsame diffusions-kontrollierte Reaktionen, wobei drei H+/e− Schritte gefolgt von einer Ligandenaustausch-Reaktion erforderlich sind, um NH3 freizusetzen: LMoVIN  LMoVNH  LMoIVNH2  LMoIIINH3 und anschließender Ligandenaustausch führt zur Freisetzung von NH3.rnIn einem weiteren Projekt wurde der Bis(ddpd)-Kupfer(II)-Komplex EPR-spektroskopisch in Hinblick auf Jahn−Teller-Verzerrung und -Dynamik untersucht. Dabei wurden die EPR-Spektren bei variabler Temperatur (70−293 K) aufgenommen. Im Festkörperspektrum bei T < 100 K erscheint der Kupfer(II)-Komplex als gestreckter Oktaeder, wohingegen das EPR-Spektrum bei höheren Temperaturen g-Werte aufzeigt, die einer pseudo-gestauchten oktaedrischen Kupfer(II)-Spezies zuzuordnen sind. Diese Tatsache wird einem intramolekularen dynamischen Jahn−Teller Phänomen zugeschrieben, welcher bei 100 K eingefroren wird.

Comparison of Model Chemistry and Density Functional Theory Thermochemical Predictions with Experiment for Formation of Ionic Clusters of the Ammonium Cation Complexed with Water and Ammonia; Atmospheric Implications

The G2, G3, CBS-QB3, and CBS-APNO model chemistry methods and the B3LYP, B3P86, mPW1PW, and PBE1PBE density functional theory (DFT) methods have been used to calculate ΔH° and ΔG° values for ionic clusters of the ammonium ion complexed with water and ammonia. Results for the clusters NH4+(NH3)n and NH4+(H2O)n, where n = 1−4, are reported in this paper and compared against experimental values. Agreement with the experimental values for ΔH° and ΔG° for formation of NH4+(NH3)n clusters is excellent. Comparison between experiment and theory for formation of the NH4+(H2O)n clusters is quite good considering the uncertainty in the experimental values. The four DFT methods yield excellent agreement with experiment and the model chemistry methods when the aug-cc-pVTZ basis set is used for energetic calculations and the 6-31G* basis set is used for geometries and frequencies. On the basis of these results, we predict that all ions in the lower troposphere will be saturated with at least one complete first hydration shell of water molecules.