@misc{SanderKlapprothYouettetal.2012, author = {Sander, Oliver and Klapproth, Corinna and Youett, Jonathan and Kornhuber, Ralf and Deuflhard, Peter}, title = {Towards an Efficient Numerical Simulation of Complex 3D Knee Joint Motion}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-14512}, number = {12-06}, year = {2012}, abstract = {We present a time-dependent finite element model of the human knee joint of full 3D geometric complexity. Its efficient numerical simulation requires advanced numerical algorithms that have been developed just recently. Up to now, the model comprises bones, cartilage, and the major ligaments (patella and menisci are still missing). Bones (femur, tibia, and fibula) are modelled by linear elastic materials, cartilage by viscoelastic materials, ligaments by one-dimensional so-called Cosserat rods. In order to capture the dynamical contact problems correctly, we solve the full PDEs of elasticity in the presence of strict contact inequalities. For the total spatio-temporal discretization we apply a method of layers approach (first time, then space discretization). For the time discretization of the elastic and viscoelastic parts, we apply a new contact-stabilized Newmark method, while for the Cosserat rods we choose an energy-momentum method. For the space discretization, we use linear finite elements for the elastic and viscoelastic parts and novel geodesic finite elements for the Cosserat rods. The coupled system is solved by a Dirichlet-Neumann method, and the arising large algebraic systems are solved by a recent fast multigrid solver, the truncated non-smooth Newton multigrid method.}, language = {en} } @misc{NaegelDeuflhardWittum2018, author = {N{\"a}gel, Arne and Deuflhard, Peter and Wittum, Gabriel}, title = {Efficient Stiff Integration of Density Driven Flow Problems}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-70891}, year = {2018}, abstract = {The paper investigates the efficient use of a linearly implicit stiff integrator for the numerical solution of density driven flow problems. Upon choosing a one-step method of extrapolation type (code LIMEX), the use of full Jacobians and reduced approximations are discussed. Numerical experiments include nonlinear density flow problems such as diffusion from a salt dome (2D), a (modified) Elder problem (3D), the saltpool benchmark (3D) and a real life salt dome problem (2D). The arising linear equations are solved using either a multigrid preconditioner from the software package UG4 or the sparse matrix solver SuperLU.}, language = {en} }