TY - GEN A1 - Erdmann, Bodo A1 - Lang, Jens A1 - Roitzsch, Rainer T1 - KARDOS - User"s Guide N2 - The adaptive finite element code {\sc Kardos} solves nonlinear parabolic systems of partial differential equations. It is applied to a wide range of problems from physics, chemistry, and engineering in one, two, or three space dimensions. The implementation is based on the programming language C. Adaptive finite element techniques are employed to provide solvers of optimal complexity. This implies a posteriori error estimation, local mesh refinement, and preconditioning of linear systems. Linearely implicit time integrators of {\em Rosenbrock} type allow for controlling the time steps adaptively and for solving nonlinear problems without using {\em Newton's} iterations. The program has proved to be robust and reliable. The user's guide explains all details a user of {\sc Kardos} has to consider: the description of the partial differential equations with their boundary and initial conditions, the triangulation of the domain, and the setting of parameters controlling the numerical algorithm. A couple of examples makes familiar to problems which were treated with {\sc Kardos}. We are extending this guide continuously. The latest version is available by network: {\begin{rawhtml} Downloads. \end{rawhtml}} T3 - ZIB-Report - 02-42 KW - nonlinear parabolic systems KW - Rosenbrock methods KW - adaptivity in time and space Y1 - 2002 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7096 ER - TY - GEN A1 - Kober, Cornelia A1 - Erdmann, Bodo A1 - Sader, Robert A1 - Zeilhofer, Hans-Florian T1 - Simulation of the Human Mandible: Comparison of Bone Mineral Density and Stress/Strain Profiles due to Masticatory Muscles' Traction N2 - The correlation of the inner architecture of bone and its functional loading was already stated by Wolff in 1892. Our objective is to demonstrate this interdependence in the case of the human mandible. For this purpose, stress/strain profiles occuring at a human lateral bite were simulated. Additionally, by a combination of computer graphics modules, a three--dimensional volumetric visualization of bone mineral density could be given. Qualitative correspondences between the density profile of the jaw and the simulated stress/strain profiles could be pointed out. In the long run, this might enable the use of the simulation for diagnosis and prognosis. The solution of the underlying partial differential equations describing linear elastic material behaviour was provided by an adaptive finite element method. Estimates of the discretization errors, local grid refinement, and multilevel techniques guaranteed the reliability and efficiency of the method. T3 - ZIB-Report - 03-23 KW - linear elasticity KW - finite element method KW - adaptive grid refinement KW - stress/strain analysis KW - human mandible Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7458 ER - TY - GEN A1 - Kober, Cornelia A1 - Erdmann, Bodo A1 - Lang, Jens A1 - Sader, Robert A1 - Zeilhofer, Hans-Florian T1 - Adaptive Finite Element Simulation of the Human Mandible Using a New Physiological Model of the Masticatory Muscles N2 - Structural mechanics simulation of bony organs is of general medical and biomechanical interest, because of the interdependence of the inner architecture of bone and its functional loading already stated by Wolff in 1892. This work is part of a detailed research project concerning the human mandible. By adaptive finite element techniques, stress/strain profiles occurring in the bony structure under biting were simulated. Estimates of the discretization errors, local grid refinement, and multilevel techniques guarantee the reliability and efficiency of the method. In general, our simulation requires a representation of the organ's geometry, an appropriate material description, and the load case due to teeth, muscle, or joint forces. In this paper, we want to focus on the influence of the masticatory system. Our goal is to capture the physiological situation as far as possible. By means of visualization techniques developed by the group, we are able to extract individual muscle fibres from computed tomography data. By a special algorithm, the fibres are expanded to fanlike (esp. for the musc. temporalis) coherent vector fields similar to the anatomical reality. The activity of the fibres can be adapted according to compartmentalisation of the muscles as measured by electromyological experiments. A refined sensitivity analysis proved remarkable impact of the presented approach on the simulation results. T3 - ZIB-Report - 04-16 KW - adaptive finite elements KW - error estimation KW - local grid refinement KW - isotropic linear elasticity KW - human mandible KW - masticatory muscles Y1 - 2004 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7917 ER - TY - GEN A1 - Kober, Cornelia A1 - Erdmann, Bodo A1 - Hellmich, Christian A1 - Sader, Robert A1 - Zeilhofer, Hans-Florian T1 - Anisotropic Simulation of the Human Mandible N2 - We focus on the role of anisotropic elasticity in the simulation of the load distribution in a human mandible due to a lateral bite on the leftmost premolar. Based on experimental evidence, we adopt ``local''" orthotropy of the elastic properties of the bone tissue. Since the trajectories of anisotropic elasticity are not accessible from Computer Tomographic (CT) data, they will be reconstructed from (i) the organ's geometry and (ii) from coherent structures which can be recognized from the spatial distribution of the CT values. A sensitivity analysis comprising various 3D FE simulations reveals the relevance of elastic anisotropy for the load carrying behavior of a human mandible: Comparison of the load distributions in isotropic and anisotropic simulations indicates that anisotropy seems to ``spare''" the mandible from loading. Moreover, a maximum degree of anisotropy leads to kind of an load minimization of the mandible, expressed by a minimum of different norms of local strain, evaluated throughout the organ. Thus, we may suggest that anisotropy is not only relevant, but also in some sense ``optimal''. T3 - ZIB-Report - 04-12 KW - anisotropy KW - elasticity KW - bony organ simulation KW - human mandible KW - visualization of tissue quality KW - finite element analysis Y1 - 2004 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7873 ER - TY - GEN A1 - Weiser, Martin A1 - Röllig, Mathias A1 - Arndt, Ralf A1 - Erdmann, Bodo T1 - Development and test of a numerical model for pulse thermography in civil engineering N2 - Pulse thermography of concrete structures is used in civil engineering for detecting voids, honeycombing and delamination. The physical situation is readily modeled by Fourier's law. Despite the simplicity of the PDE structure, quantitatively realistic numerical 3D simulation faces two major obstacles. First, the short heating pulse induces a thin boundary layer at the heated surface which encapsulates all information and therefore has to be resolved faithfully. Even with adaptive mesh refinement techniques, obtaining useful accuracies requires an unsatisfactorily fine discretization. Second, bulk material parameters and boundary conditions are barely known exactly. We address both issues by a semi-analytic reformulation of the heat transport problem and by parameter identification. Numerical results are compared with measurements of test specimens. T3 - ZIB-Report - 08-45 KW - finite element model KW - KARDOS KW - pulse thermography KW - civil engineering KW - inverse solution Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10980 SN - 1438-0064 ER - TY - GEN A1 - Bornemann, Folkmar A. T1 - Adaptive Solution of One-Dimensional Scalar Conservation Laws with Convex Flux. N2 - A new adaptive approach for one-dimensional scalar conservation laws with convex flux is proposed. The initial data are approximated on an adaptive grid by a problem dependent, monotone interpolation procedure in such a way, that the multivalued problem of characteristic transport can be easily and explicitly solved. The unique entropy solution is chosen by means of a selection criterion due to LAX. For arbitrary times, the solutions is represented by an adaptive monotone spline interpolation. The spatial approximation is controlled by local $L^1$-error estimated. As a distinctive feature of the approach, there is no discretization in time. The method is monotone on fixed grids. Numerical examples are included, to demonstrate the predicted behavior. {\bf Key words.} method of characteristics, adaptive grids, monotone interpolation, $L^1$-error estimates {\bf AMS(MOS) subject classification.} 65M15, 65M25, 65M50. T3 - ZIB-Report - SC-92-18 KW - method of characteristics KW - adaptive grids KW - monotone interpolation KW - L1-error estimates Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-880 ER - TY - GEN A1 - Lang, Jens A1 - Walter, Artur T1 - A Finite Element Method Adaptive in Space and Time for Nonlinear Reaction-Diffusion- Systems. N2 - Large scale combustion simulations show the need for adaptive methods. First, to save computation time and mainly to resolve local and instationary phenomena. In contrast to the widespread method of lines, we look at the reaction- diffusion equations as an abstract Cauchy problem in an appropriate Hilbert space. This means, we first discretize in time, assuming the space problems solved up to a prescribed tolerance. So, we are able to control the space and time error separately in an adaptive approach. The time discretization is done by several adaptive Runge-Kutta methods whereas for the space discretization a finite element method is used. The different behaviour of the proposed approaches are demonstrated on many fundamental examples from ecology, flame propagation, electrodynamics and combustion theory. {\bf Keywords:} initial boundary value problem, Rothe- method, adaptive Runge-Kutta method, finite elements, mesh refinement. {\bf AMS CLASSIFICATION:} 65J15, 65M30, 65M50. T3 - ZIB-Report - SC-92-05 KW - initial boundary value problem KW - Rothe method KW - adaptive Runge-Kutta method KW - finite elements KW - mesh refinement Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-755 ER - TY - GEN A1 - Bornemann, Folkmar A. T1 - An Adaptive Multilevel Approach to Parabolic Equations III. N2 - Part III of the paper is devoted to the construction of an adaptive FEM solver in two spatial dimensions, which is able to handle the singularly perturbed elliptic problems arising from discretization in time. The problems of error estimation and multilevel iterative solution of the linear systems - both uniformly well behaved with respect to the time step - can be solved simultaneously within the framework of preconditioning. A multilevel nodal basis preconditioner able to handle highly nonuniform meshes is derived. As a numerical example an application of the method to the bioheat-transfer equation is included. {\bf AMS CLASSIFICATION:} 65F10, 65F35, 65M50, 65M60, 65N30. T3 - ZIB-Report - SC-91-01 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-519 ER - TY - GEN A1 - Weiser, Martin T1 - Pointwise Nonlinear Scaling for Reaction-Diffusion-Equations N2 - Parabolic reaction-diffusion systems may develop sharp moving reaction fronts which pose a challenge even for adaptive finite element methods. We propose a method to transform the equation into an equivalent form that usually exhibits solutions which are easier to discretize, giving higher accuracy for a given number of degrees of freedom. The transformation is realized as an efficiently computable pointwise nonlinear scaling that is optimized for prototypical planar travelling wave solutions of the underlying reaction-diffusion equation. The gain in either performance or accuracy is demonstrated on different numerical examples. T3 - ZIB-Report - 07-45 KW - reaction-diffusion equations KW - travelling waves KW - nonlinear scaling KW - discretization error Y1 - 2007 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10493 SN - 1438-0064 ER - TY - GEN A1 - Weiser, Martin A1 - Scacchi, Simone T1 - Spectral Deferred Correction methods for adaptive electro-mechanical coupling in cardiac simulation N2 - We investigate spectral deferred correction (SDC) methods for time stepping and their interplay with spatio-temporal adaptivity, applied to the solution of the cardiac electro-mechanical coupling model. This model consists of the Monodomain equations, a reaction-diffusion system modeling the cardiac bioelectrical activity, coupled with a quasi-static mechanical model describing the contraction and relaxation of the cardiac muscle. The numerical approximation of the cardiac electro-mechanical coupling is a challenging multiphysics problem, because it exhibits very different spatial and temporal scales. Therefore, spatio-temporal adaptivity is a promising approach to reduce the computational complexity. SDC methods are simple iterative methods for solving collocation systems. We exploit their flexibility for combining them in various ways with spatio-temporal adaptivity. The accuracy and computational complexity of the resulting methods are studied on some numerical examples. T3 - ZIB-Report - 14-22 Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-50695 SN - 1438-0064 ER - TY - GEN A1 - Erdmann, Bodo A1 - Lang, Jens A1 - Matera, Sebastian A1 - Wilmanski, Krzysztof T1 - Adaptive Linearly Implicit Methods for Linear Poroelastic Equations N2 - Adaptive numerical methods in time and space are introduced and studied for linear poroelastic models in two and three space dimensions. We present equivalent models for linear poroelasticity and choose both the {\em displacement--pressure} and the {\em stress--pressure} formulation for our computations. Their discretizations are provided by means of linearly implicit schemes in time and linear finite elements in space. Our concept of adaptivity opens a way to a fast and reliable simulation of different loading cases defined by corresponding boundary conditions. We present some examples using our code {\sf Kardos} and show that the method works efficiently. In particular, it could be used in the simulation of some bone healing models. T3 - ZIB-Report - 06-37 KW - Poroelasticity KW - Biot's model KW - bone healing KW - adaptive finite elements KW - adaptive time integration KW - Rosenbrock methods Y1 - 2006 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-9312 ER - TY - GEN A1 - Fackeldey, Konstantin T1 - Coupling Meshbased and Meshfree Methods by a Transfer Operator Approach N2 - In contrast to the well known meshbased methods like the finite element method, meshfree methods do not rely on a mesh. However besides their great applicability, meshfree methods are rather time consuming. Thus, it seems favorable to combine both methods, by using meshfree methods only in a small part of the domain, where a mesh is disadvantageous, and a meshbased method for the rest of the domain. We motivate, that this coupling between the two simulation techniques can be considered as saddle point problem and show the stability of this coupling. Thereby a novel transfer operator is introduced, which interacts in the transition zone, where both methods coexist. T3 - ZIB-Report - 10-12 KW - Gitterlose Methoden KW - inf-sup-Bedingung KW - Kopplung KW - meshfree KW - meshbased KW - inf-sup-Condition Y1 - 2010 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11755 SN - 1438-0064 ER -