Numerische Mathematik und Wissenschaftliches Rechnen
Seminar Numerische Mathematik / Numerical mathematics seminars
aktuelles Programm / current program
Donnerstag, 03. 05. 2018, 14:00 Uhr (ESH)
Dr. H. Stephan (WIAS Berlin)
One million perrin pseudo primes including a few giants
Pseudoprimes are integers that are no primes but behave like them in some sense. Suppose we have a theorem like the following: If n is a prime, then statement A(n) holds. In general, the opposite is not true: It may be that A(n) holds, but n is a composite number, a so-called pseudo prime with respect to statement A. Pseudoprimes are interesting if they are very rare, as for instance Perrin's pseudoprimes, the smallest of which is 271441. The talk introduces pseudoprimes which are based on recurrent sequences. In addition, some new numerical results on Perrin's pseudoprimes and a fast algorithm for their calculation are presented.
Mittwoch, 18. 04. 2018, 13:30 Uhr (ESH)
Prof. P. Deuring (Universitè du Littoral ''Côte d'Opale'', Frankreich)
Stabilitäts- und Fehlerschranken einer FE-FV-Diskretisierung von Konvektions-Diffusionsgleichung: exponentielle Abhängigkeit von Parametern
Wir betrachten eine FE-FV-Methode zur Bestimmung von Näherungslösungen von Konvektions-Diffusionsgleichungen. Bei dieser Methode wird der Diffusionsterm mittels Crouzeix-Raviart-Elementen diskretisiert, während man den Konvektionsterm durch baryzentrische finite Volumen approximiert. Wir gehen der Frage nach, ob sich Stabilitäts- und Fehlerschranken für diese Methode finden lassen, die von keiner Größe exponentiell abhängen, insbesondere nicht vom Kehrwert des Diffusionskoeffizienten.
Donnerstag, 12. 04. 2018, 14:00 Uhr (ESH)
Prof. M. J. Neilan (University of Pittsburgh, USA)
Inf-sup stable Stokes pairs on barycentric refinements producing divergence-free approximations
We construct several stable finite element pairs for the Stokes problem on barycentric refinements in arbitrary dimensions and for any polynomial degree. A key feature of the spaces is that the divergence maps the discrete velocity space onto the discrete pressure space; thus, when applied to models of incompressible flows, the pairs yield divergence-free velocity approximations. The key ingredients to prove these results are local inf-sup stability estimates and a modification of Bernardi-Raugel bubble functions. This is joint work with Johnny Guzman.
Donnerstag, 01. 03. 2018, 14:00 Uhr (ESH)
Dr. F. Dassi (Politecnico di Milano, Italien)
Recent advancements of the Virtual Element Method in 3D
The Virtual Element Method is a novel way to discretize a partial differential equation. It avoids the explicit integration of shape functions and introduces an innovative construction of the stiffness matrix so that it acquires very interesting properties and advantages. One among them is the possibility to apply the VEM to general polygonal/polyhedral domain decomposition, also characterized by non-conforming and non-convex elements.
In this talk we focus on the definition/construction of the Virtual Element functional spaces in three dimensions and how apply this new strategy to set a standard Laplacian problem in 3D.
Finally we test the method to show its robustness with respect to element distortion and the polynomial approximation degree $k$. Then, we move to more involved cases: convection-diffusion-reaction problems with variable coefficients and magnetostatic Maxwell equations.
Donnerstag, 01. 02. 2018, 14:00 Uhr (ESH)
Dr. E. Sinibaldi (Italian Institute of Technology)
Selected modeling approaches for biomedical applications and biorobotics tools
The effective deployment of theranostic agents (including drugs), smart?materials?based systems and interfaces to target regions (prospectively) in the human body where to perform the sought actions also requires theoretical models and physical tools. Modeling the delivery and the (remote) actuation/stimulation of the aforementioned agents/systems/interfaces, in particular, permits to compensate for experimental conditions hard to characterize and to better interpret the experimental results, thus paving the way for quantitative therapy design and control. This is complemented by model-based design of related experimental devices, flexible tools able to safely navigate anatomical pathways, and miniaturized effectors. In this talk I will firstly present the solution to an inverse problem, namely to determine the velocity profile in a vessel cross-section starting from the flow-rate, which is relevant to pulsatile biological flows (blood and cerebro-spinal fluid flows), with application, e.g., to magnetic particle targeting. Then I will address numerical models for intrathecal drug-delivery (drug infusion in the cerebro-spinal fluid, also accounting for transport to the spinal cord) and simplified analytical models for intra-tissue drug-delivery (particle transport into porous/poroelastic media, with application to intratumoral thermotherapy). Finally, I will overview model-based approaches for: piezoelectric nanoparticle-mediated cell stimulation; real-scale physical models of the blood-brain barrier; flexible biorobotics probes; miniature (bioinspired) actuators and effectors.und 15:00 Uhr
L. Blank (WIAS Berlin)
Towards a simple and robust finite element method for the numerical simulation of porous media flow
The topic of this talk is the numerical simulation of flow through porous media based on the Brinkman model as a unified framework that allows the transit between Darcy and Stokes problems. Therefore an unconditionally stable low order finite element approach, which is robust with respect to the physical parameters, is proposed. This approach is based on the combination of stabilized equal order finite elements with a non-symmetric (penalty-free) Nitsche method for the weak imposition of essential boundary conditions. Focusing on the two-dimensional case, optimal a priori error estimates in a mesh-dependent norm, which allows to extend the results also to the Stokes and Darcy limits, are obtained.
Donnerstag, 18. 01. 2018, 14:00 Uhr (ESH)
Prof. W. Dreyer (WIAS Berlin)
Modeling of non-Newtonian fluids without material frame indifference
In the fifties the modeling of non-Newtonian fluids initiated the search for invariant time derivatives with respect to certain space time transformations. Then Euclidean Transformations were selected to establish the Principle of Material Frame Indifference and moreover, the introduction of Nth grade fluids. However, the two concepts lead to serious inconsistencies in both thermodynamic modeling and experimental observations. In 1986 the subject has been resolved in a remarkable paper by I. Müller and K. Wilmanski. In this lecture we show how an appropriate model of non-Newtonian fluids may be incorporated in Continuum Thermodynamics as it was laid down by D. Bothe and W. Dreyer.
Donnerstag, 11. 01. 2018, 14:00 Uhr (ESH)
Prof. G. Lube (Georg-August-Universität Göttingen)
Why exactly divergence-free H(div)-conforming FEM for transient incompressible flows?
Recent results show that one can reconstruct classical inf-sup stable H1-conforming FEM for incompressible flow problems to be pressure-robust. Exactly divergence-free FEM are pressure-robust per construction. In particular, exactly divergence-free H(div)-conforming FEM combine excellent stability and conservation properties with minimal stabilization. We show that the numerical analysis of such methods allows to separate the (static) linear Stokes regime and the nonlinear dynamic regime. Semi-robust error estimates w.r.t. the Reynolds number will be given. In a hybridized form, H(div)-conforming FEM are suitable for large scale computation. Applications of the approach to two- and three-dimensional problems of vortex dynamics will be presented for high Reynolds numbers.und 15:00 Uhr
Dr. A. Rasheed (Lahore University of Management Sciences, Pakistan)
Influence of magnetic field on dendrites during solidification of binary mixtures
A phase field model has been proposed which incorporates convection and magnetic field in an isothermal environment. A numerical scheme is proposed and numerical analysis of model in two-dimensional geometry is performed. The numerical stability and error analysis of this approximation scheme which is based on mixed finite-element method are performed. An application of a nickel-copper binary alloy is considered. Influence of various magnetic fields on the dendrites during the solidification process has been discussed.