@article{WeberFischerDamerauetal., author = {Weber, Marie-Christin and Fischer, Lisa and Damerau, Alexandra and Ponomarev, Igor and Pfeiffenberger, Moritz and Gaber, Timo and G{\"o}tschel, Sebastian and Lang, Jens and R{\"o}blitz, Susanna and Buttgereit, Frank and Ehrig, Rainald and Lang, Annemarie}, title = {In vitro and in silico modeling of cellular and matrix-related changes during the early phase of osteoarthritis}, series = {BioRxiv}, journal = {BioRxiv}, doi = {10.1101/725317}, abstract = {Understanding the pathophysiological processes of osteoarthritis (OA) require adequate model systems. Although different in vitro or in vivo models have been described, further comprehensive approaches are needed to study specific parts of the disease. This study aimed to combine in vitro and in silico modeling to describe cellular and matrix-related changes during the early phase of OA. We developed an in vitro OA model based on scaffold-free cartilage-like constructs (SFCCs), which was mathematically modeled using a partial differential equation (PDE) system to resemble the processes during the onset of OA. SFCCs were produced from mesenchymal stromal cells and analyzed weekly by histology and qPCR to characterize the cellular and matrix-related composition. To simulate the early phase of OA, SFCCs were treated with interleukin-1β (IL-1β), tumor necrosis factor α (TNFα) and examined after 3 weeks or cultivated another 3 weeks without inflammatory cytokines to validate the regeneration potential. Mathematical modeling was performed in parallel to the in vitro experiments. SFCCs expressed cartilage-specific markers, and after stimulation an increased expression of inflammatory markers, matrix degrading enzymes, a loss of collagen II (Col-2) and a reduced cell density was observed which could be partially reversed by retraction of stimulation. Based on the PDEs, the distribution processes within the SFCCs, including those of IL-1β, Col-2 degradation and cell number reduction was simulated. By combining in vitro and in silico methods, we aimed to develop a valid, efficient alternative approach to examine and predict disease progression and new therapeutic strategies.}, language = {en} } @article{WeberFischerDamerauetal., author = {Weber, Marie-Christin and Fischer, Lisa and Damerau, Alexandra and Ponomarev, Igor and Pfeiffenberger, Moritz and Gaber, Timo and G{\"o}tschel, Sebastian and Lang, Jens and R{\"o}blitz, Susanna and Buttgereit, Frank and Ehrig, Rainald and Lang, Annemarie}, title = {Macroscale mesenchymal condensation to study cytokine-driven cellular and matrix-related changes during cartilage degradation}, series = {Biofabrication}, volume = {12}, journal = {Biofabrication}, number = {4}, doi = {10.1088/1758-5090/aba08f}, abstract = {Understanding the pathophysiological processes of cartilage degradation requires adequate model systems to develop therapeutic strategies towards osteoarthritis (OA). Although different in vitro or in vivo models have been described, further comprehensive approaches are needed to study specific disease aspects. This study aimed to combine in vitro and in silico modeling based on a tissue-engineering approach using mesenchymal condensation to mimic cytokine-induced cellular and matrix-related changes during cartilage degradation. Thus, scaffold-free cartilage-like constructs (SFCCs) were produced based on self-organization of mesenchymal stromal cells (mesenchymal condensation) and i) characterized regarding their cellular and matrix composition or secondly ii) treated with interleukin-1β (IL-1β) and tumor necrosis factor α (TNFα) for 3 weeks to simulate OA-related matrix degradation. In addition, an existing mathematical model based on partial differential equations was optimized and transferred to the underlying settings to simulate distribution of IL-1β, type II collagen degradation and cell number reduction. By combining in vitro and in silico methods, we aim to develop a valid, efficient alternative approach to examine and predict disease progression and effects of new therapeutics.}, language = {en} } @misc{LangMerz, author = {Lang, Jens and Merz, Wilhelm}, title = {Two-Dimensional Adaptive Simulation of Dopant Diffusion in Silicon}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-5716}, number = {00-03}, abstract = {One important step in the fabrication of silicon-based integrated circuits is the creation of semiconducting areas by diffusion of dopant impurities into silicon. Complex models have been developed to investigate the redistribution of dopants and point defects. In general, numerical analysis of the resulting PDEs is the central tool to assess the modelling process. We present an adaptive approach which is able to judge the quality of the numerical approximation and which provides an automatic mesh improvement. Using linearly implicit methods in time and multilevel finite elements in space, we are able to integrate efficiently the arising reaction-drift-diffusion equations with high accuracy. Two different diffusion processes of practical interest are simulated.}, language = {en} } @misc{ErdmannLangRoitzsch, author = {Erdmann, Bodo and Lang, Jens and Roitzsch, Rainer}, title = {KARDOS - User"s Guide}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7096}, number = {02-42}, abstract = {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}}}, language = {en} } @misc{ErdmannKoberLangetal., author = {Erdmann, Bodo and Kober, Cornelia and Lang, Jens and Sader, Robert and Zeilhofer, Hans-Florian and Deuflhard, Peter}, title = {Efficient and Reliable Finite Element Methods for Simulation of the Human Mandible}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6403}, number = {01-14}, abstract = {By computed tomography data (CT), the individual geometry of the mandible is quite well reproduced, also the separation between cortical and trabecular bone. Using anatomical knowledge about the architecture and the functional potential of the masticatory muscles, realistic situations were approximated. The solution of the underlying partial differential equations describing linear elastic material behaviour is provided by an adaptive finite element method. Estimations of the discretization error, local grid refinement, and multilevel techniques guarantee the reliability and efficiency of the method.}, language = {en} } @misc{KoberErdmannLangetal., author = {Kober, Cornelia and Erdmann, Bodo and Lang, Jens and Sader, Robert and Zeilhofer, Hans-Florian}, title = {Adaptive Finite Element Simulation of the Human Mandible Using a New Physiological Model of the Masticatory Muscles}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7917}, number = {04-16}, abstract = {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.}, language = {en} } @misc{Lang, author = {Lang, Jens}, title = {Adaptive FEM for Reaction-Diffusion Equations}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-2393}, number = {SC-96-28}, abstract = {An integrated time--space adaptive finite element method for solving mixed systems of nonlinear parabolic, elliptic, and differential algebraic equations is presented. The approach is independent of the spatial dimension. For the discretization in time we use singly diagonally linearly implicit Runge--Kutta methods of Rosenbrock type. Local time errors for the step size control are defined by an embedded strategy. A multilevel finite element Galerkin method is subsequently applied for the discretization in space. A posteriori estimates of local spatial discretization errors are obtained solving local problems with higher order approximation. Superconvergence arguments allow to simplify the required computations. Two different strategies to obtain the start grid of the multilevel process are compared. The devised method is applied to a solid--solid combustion problem.}, language = {en} } @misc{FroehlichLang, author = {Fr{\"o}hlich, Jochen and Lang, Jens}, title = {Twodimensional Cascadic Finite Element Computations of Combustion Problems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-2167}, number = {SC-96-05}, abstract = {We present an integrated time--space adaptive finite element method for solving systems of twodimensional nonlinear parabolic systems in complex geometry. The partial differential system is first discretized in time using a singly linearly implicit Runge--Kutta method of order three. Local time errors for the step size control are defined by an embedding strategy. These errors are used to propose a new time step by a PI controller algorithm. A multilevel finite element method with piecewise linear functions on unstructured triangular meshes is subsequently applied for the discretization in space. The local error estimate of the finite element solution steering the adaptive mesh refinement is obtained solving local problems with quadratic trial functions located essentially at the edges of the triangulation. This two--fold adaptivity successfully ensures an a priori prescribed tolerance of the solution. The devised method is applied to laminar gaseous combustion and to solid--solid alloying reactions. We demonstrate that for such demanding applications the employed error estimation and adaption strategies generate an efficient and versatile algorithm.}, language = {en} } @misc{Lang, author = {Lang, Jens}, title = {KARDOS - KAskade Reaction Diffusion One-dimensional System.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-5019}, number = {TR-93-09}, abstract = {A software package for the adaptive solution of time--dependent reaction--diffusion systems and linear elliptic systems in one space dimension is presented. The used algorithm is based on fundamental arguments in J.~Lang, A.~Walter: {\it A Finite Element Method Adaptive in Space and Time for Nonlinear Reaction--Diffusion Systems.} IMPACT of Computing in Science and Engineering, 4, p.~269--314 (1992). Here, only brief outlines of the algorithm are given. This software package is based on the KASKADE toolbox B.~Erdmann, J.~Lang, R.~Roitzsch: {\it KASKADE -- Manual.} To appear as Technical Report TR 93--5, Konrad--Zuse--Zentrum (ZIB) (1993).}, language = {en} } @misc{LangErdmannRoitzsch, author = {Lang, Jens and Erdmann, Bodo and Roitzsch, Rainer}, title = {Three-Dimensional Fully Adaptive Solution of Thermo-Diffusive Flame Propagation Problems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-2879}, number = {SC-97-18}, abstract = {In this paper we present a self--adaptive finite element method to solve flame propagation problems in 3D. An implicit time integrator of Rosenbrock type is coupled with a multilevel approach in space. The proposed method is applied to an unsteady thermo--diffusive combustion model to demonstrate its potential for the solution of complicated problems.}, language = {en} } @misc{ErdmannLangRoitzsch, author = {Erdmann, Bodo and Lang, Jens and Roitzsch, Rainer}, title = {Kaskade Manual - Version 2.0.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-4970}, number = {TR-93-05}, abstract = {The KASKADE toolbox defines an interface to a set of C subroutines which can be used to implement adaptive multilevel Finite Element Methods solving systems of elliptic equations in two and three space dimensions. The manual contains the description of the data structures and subroutines. The main modules of the toolbox are a runtime environment, triangulation and node handling, assembling, direct and iterative solvers for the linear systems, error estimators, refinement strategies, and graphic utilities. Additionally, we included appendices on the basic command language interface, on file formats, and on the definition of the partial differential equations which can be solved. The software is available on the ZIB ftp--server {\tt elib} in the directory {\tt pub/kaskade}. TR 93--5 supersedes TR 89--4 and TR 89--05.}, language = {en} } @misc{Lang, author = {Lang, Jens}, title = {An Adaptive Finite Element Method for Convection-Diffusion Problems by Interpolation Techniques.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-4795}, number = {TR-91-04}, abstract = {For adaptive solution of convection- difussion problems with the streamline-diffusion finite element method, an error estimator based on interpolation techniques is developed. It can be shown that for correctness of this error estimator a restriction of the maximum angle is to be sufficient. Compared to usual methods, the adaptive process leads to more accurate solutions at much less computational cost. Numerical tests are enclosed. {\bf Keywords: } Adaptive finite elements, convection- diffusion equation, internal and boundary layers, streamline-diffusion. {\bf Subject Classifications:} AMS(MOS): 65N15, 65N30}, language = {en} } @misc{SeewaldPolleiKraumeetal., author = {Seewald, Frank and Pollei, Anja and Kraume, Matthias and Mittelbach, Walter and Lang., Jens}, title = {Numerical Calculation of the Heat Transfer in an Adsorption Energy Storage with KARDOS}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3927}, number = {SC-99-04}, abstract = {A new seasonal energy storage for thermal solar systems has been developed on the basis of an adsorption-desorption process. Design and optimization of this storage will be supported by numerical simulations of heat and mass transfer with KARDOS. This paper focuses on the unsteady heat transfer during the major operating step of energetic discharge of the storage, which is characterized by conductive heat transfer in the fixed bed and a strong heat source caused by the adsorption enthalpy. Results are interpreted concerning the influence of variations in the parameter set. The method of implementation of the differential equation will be shown as well as the post-processing and gridwriting programs.}, language = {en} } @misc{ErdmannLangSeebass, author = {Erdmann, Bodo and Lang, Jens and Seebass, Martin}, title = {Adaptive Solutions of Nonlinear Parabolic Equations with Application to Hyperthermia Treatments}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3130}, number = {SC-97-44}, abstract = {We present a self-adaptive finite element method to solve nonlinear evolution problems in 3D. An implicit time integrator of Rosenbrock type is coupled with a multilevel approach in space. The proposed method is applied to hyperthermia treatments to demonstrate its potential for the solving of complicated problems.}, language = {en} } @misc{LangMerz, author = {Lang, Jens and Merz, Wilhelm}, title = {Numerical Simulation of Single Species Dopant Diffusion in Silicon under Extrinsic Conditions}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3163}, number = {SC-97-47}, abstract = {In this article we consider a general model for phosphorus diffusion in silicon under extrinsic doping conditions. At such high concentrations we have to include the charged species and the internal electric field of the crystal, both of which can have profound effects on diffusion. In principle, this leads to a very large number of drift--diffusion--reaction equations: one for each charge state of every species, plus one Poisson equation to describe the internal electric field (in terms of the electron/hole concentration). The number of equations can be reduced substantially by making additional assumptions on the distribution of charge states and local equilibrium assumptions concerning the reaction terms. The resulting model turns out to be very interesting for numerical investigation. We solve the problem numerically in two space dimensions with the adaptive finite element program KARDOS and describe the numerical method used here to treat the resulting drift--diffusion--reaction problem.}, language = {en} } @misc{DeuflhardLangNowak, author = {Deuflhard, Peter and Lang, Jens and Nowak, Ulrich}, title = {Adaptive Algorithms in Dynamical Process Simulation}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-1823}, number = {SC-95-16}, abstract = {Dynamical simulation of industrially relevant processes strongly advises the use of algorithms, which are {\em adaptive} both in time and in space discretization. The paper presents two alternatives: (a) a fully adaptive method of lines approach, which is based on finite difference methods and essentially applicable to 1D problems; (b) a fully adaptive Rothe method, which is based on a fast multilevel finite element method and applicable to 1D up to 3D.}, language = {en} } @misc{LourencoRosaCastroetal., author = {Lourenco, Maria Jos{\´e} and Rosa, Samuel Costa S. and Castro, Carlos Alberto Nieto de and Albuquerque, C. and Erdmann, Bodo and Lang, Jens and Roitzsch, Rainer}, title = {Simulation of the Transient Heating in an Unsymmetrical Coated Hot--Strip Sensor with a Self--Adaptive Finite Element Method}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3656}, number = {SC-98-22}, abstract = {The transient heating in an unsymmetrical coated hot--strip sensor was simulated with a self--adaptive finite element method. The first tests of this model show that it can determine with a small error the thermal conductivity of liquids, from the transient temperature rise in the hot--strip, deposited in a substrate and coated by an alumina spray.}, language = {en} } @misc{LangErdmannRoitzsch, author = {Lang, Jens and Erdmann, Bodo and Roitzsch, Rainer}, title = {Adaptive Time-Space Discretization for Combustion Problems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3094}, number = {SC-97-40}, abstract = {We present a self--adaptive finite element method to solve combustion problems in 1D, 2D, and 3D. An implicit time integrator of Rosenbrock type is coupled with a multilevel approach in space. A posteriori error estimates are obtained by constructing locally higher order solutions involving all variables of the problem. Adaptive strategies such as step size control, spatial refinement and coarsening allow us to get economically an accurate solution. Various examples are presented to demonstrate practical applications of the proposed method.}, language = {en} } @misc{LangWalter, author = {Lang, Jens and Walter, Artur}, title = {A Finite Element Method Adaptive in Space and Time for Nonlinear Reaction-Diffusion- Systems.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-755}, number = {SC-92-05}, abstract = {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.}, language = {en} } @misc{LangWalter, author = {Lang, Jens and Walter, Artur}, title = {An Adaptive Discontinuous Finite Element Method for the Transport Equation.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-579}, number = {SC-91-07}, abstract = {In this paper we introduce a discontinuous finite element method. In our approach, it is possible to combine the advantages of finite element and finite difference methods. The main ingredients are numerical flux approximation and local orthogonal basis functions. The scheme is defined on arbitrary triangulations and can be easily extended to nonlinear problems. Two different error indicators are derived. Especially the second one is closely connected to our approach and able to handle arbitrary variing flow directions. Numerical results are given for boundary value problems in two dimensions. They demonstrate the performance of the scheme, combined with the two error indicators. {\bf Key words:} neutron transport equation, discontinuous finite element, adaptive grid refinement. {\bf Subject classifications:} AMS(MOS) 65N30, 65M15.}, language = {en} }