Tim Ricken

Multigenerational interstitial growth of biological tissues

This study formulates a theory for multigenerational interstitial growth of biological tissues whereby each generation has a distinct reference configuration determined at the time of its deposition. In this model, the solid matrix of a growing tissue consists of a multiplicity of intermingled porous permeable bodies, each of which represents a generation, all of which are constrained to move together in the current configuration. Each generation’s reference configuration has a one-to-one mapping with the master reference configuration, which is typically that of the first generation. This mapping is postulated based on a constitutive assumption with regard to that generations’ state of stress at the time of its deposition. For example, the newly deposited generation may be assumed to be in a stress-free state, even though the underlying tissue is in a loaded configuration. The mass content of each generation may vary over time as a result of growth or degradation, thereby altering the material properties of the tissue. A finite element implementation of this framework is used to provide several illustrative examples, including interstitial growth by cell division followed by matrix turnover.

A biphasic model for sinusoidal liver perfusion remodeling after outflow obstruction

Liver resection can lead to focal outflow obstruction due to transection of hepatic veins. Outflow obstruction may cause additional damage to the small remnant liver. Drainage of the obstructed territories is reestablished via dilatation of sinusoids. Subsequently, sinusoidal canals are formed draining the blood from the obstructed territory to the neighboring unobstructed territories. We raised the phenomenological hypothesis that the blood pressure gradient is the main driving force for the formation of sinusoidal vascular canals. We generated a biphasic mechanical model to describe this vascular remodeling process in relation to the variable pressure gradient. Therefore, we introduced a transverse isotropic permeability relation as well as an evolutional optimization rule to describe the relationship between pressure gradient and the direction of the sinusoidal blood flow in the fluid phase. As a next step, we developed a framework for the calculation concept including the representation of the governing weak formulations. Then, we examined a representative numerical example with simulation of the blood flow under both conditions, the physiological situation as well as after outflow obstruction. Doing so, we were able to reproduce numerically the experimentally observed process of reestablishing hepatic venous drainage via redirection of blood flow and formation of new vascular structures in respect to the fluid flow. The calculated results support the hypothesis that the reorientation of blood flow mainly depends on the pressure gradient. Further investigations are needed to determine the micromechanical influences on the reorientation of the sinusoids.

A hyperelastic biphasic fibre-reinforced model of articular cartilage considering distributed collagen fibre orientations: continuum basis, computational aspects and applications

Cartilage is a multi-phase material composed of fluid and electrolytes (68–85% by wet weight), proteoglycans (5–10% by wet weight), chondrocytes, collagen fibres and other glycoproteins. The solid phase constitutes an isotropic proteoglycan gel and a fibre network of predominantly type II collagen, which provides tensile strength and mechanical stiffness. The same two components control diffusion of the fluid phase, e.g. as visualised by diffusion tensor MRI: (i) the proteoglycan gel (giving a baseline isotropic diffusivity) and (ii) the highly anisotropic collagenous fibre network. We propose a new constitutive model and finite element implementation that focus on the essential load-bearing morphology: an incompressible, poroelastic solid matrix reinforced by an inhomogeneous, dispersed fibre fabric, which is saturated with an incompressible fluid residing in strain-dependent pores of the collagen–proteoglycan solid matrix. The inhomogeneous, dispersed fibre fabric of the solid further influences the fluid permeability, as well as an intrafibrillar portion that cannot be ‘squeezed out’ from the tissue. Using representative numerical examples on the mechanical response of cartilage, we reproduce several features that have been demonstrated experimentally in the cartilage mechanics literature.

Structural Analysis of Articular Cartilage Using Multiphoton Microscopy: Input for Biomechanical Modeling

The 3-D morphology of chicken articular cartilage was quantified using multiphoton microscopy (MPM) for use in continuum-mechanical modeling. To motivate this morphological study we propose aspects of a new, 3-D finite strain constitutive model for articular cartilage focusing on the essential load-bearing morphology: an inhomogeneous, poro-(visco)elastic solid matrix reinforced by an anisotropic, (visco)elastic dispersed fiber fabric which is saturated by an incompressible fluid residing in strain-dependent pores. Samples of fresh chicken cartilage were sectioned in three orthogonal planes and imaged using MPM, specifically imaging the collagen fibers using second harmonic generation. Employing image analysis techniques based on Fourier analysis, we derived the principal directionality and dispersion of the collagen fiber fabric in the superficial layer. In the middle layer, objective thresholding techniques were used to extract the volume fraction occupied by extracellular collagen matrix. In conjunction with information available in the literature, or additional experimental testing, we show how this data can be used to derive a 3-D map of the initial solid volume fraction and Darcy permeability.

Modeling function-perfusion behavior in liver lobules including tissue, blood, glucose, lactate and glycogen by use of a coupled two-scale PDE-ODE approach.

This study focuses on a two-scale, continuum multicomponent model for the description of blood perfusion and cell metabolism in the liver. The model accounts for a spatial and time depending hydro-diffusion–advection–reaction description. We consider a solid-phase (tissue) containing glycogen and a fluid-phase (blood) containing glucose as well as lactate. The five-component model is enhanced by a two-scale approach including a macroscale (sinusoidal level) and a microscale (cell level). The perfusion on the macroscale within the lobules is described by a homogenized multiphasic approach based on the theory of porous media (mixture theory combined with the concept of volume fraction). On macro level, we recall the basic mixture model, the governing equations as well as the constitutive framework including the solid (tissue) stress, blood pressure and solutes chemical potential. In view of the transport phenomena, we discuss the blood flow including transverse isotropic permeability, as well as the transport of solute concentrations including diffusion and advection. The continuum multicomponent model on the macroscale finally leads to a coupled system of partial differential equations (PDE). In contrast, the hepatic metabolism on the microscale (cell level) was modeled via a coupled system of ordinary differential equations (ODE). Again, we recall the constitutive relations for cell metabolism level. A finite element implementation of this framework is used to provide an illustrative example, describing the spatial and time-depending perfusion–metabolism processes in liver lobules that integrates perfusion and metabolism of the liver.

A microstructurally based continuum model of cartilage viscoelasticity and permeability incorporating measured statistical fiber orientations.

The remarkable mechanical properties of cartilage derive from an interplay of isotropically distributed, densely packed and negatively charged proteoglycans; a highly anisotropic and inhomogeneously oriented fiber network of collagens; and an interstitial electrolytic fluid. We propose a new 3D finite strain constitutive model capable of simultaneously addressing both solid (reinforcement) and fluid (permeability) dependence of the tissue’s mechanical response on the patient-specific collagen fiber network. To represent fiber reinforcement, we integrate the strain energies of single collagen fibers—weighted by an orientation distribution function (ODF) defined over a unit sphere—over the distributed fiber orientations in 3D. We define the anisotropic intrinsic permeability of the tissue with a structure tensor based again on the integration of the local ODF over all spatial fiber orientations. By design, our modeling formulation accepts structural data on patient-specific collagen fiber networks as determined via diffusion tensor MRI. We implement our new model in 3D large strain finite elements and study the distributions of interstitial fluid pressure, fluid pressure load support and shear stress within a cartilage sample under indentation. Results show that the fiber network dramatically increases interstitial fluid pressure and focuses it near the surface. Inhomogeneity in the tissue’s composition also increases fluid pressure and reduces shear stress in the solid. Finally, a biphasic neo-Hookean material model, as is available in commercial finite element codes, does not capture important features of the intra-tissue response, e.g., distributions of interstitial fluid pressure and principal shear stress.

A finite element simulation of biological conversion processes in landfills

Landfills are the most common way of waste disposal worldwide. Biological processes convert the organic material into an environmentally harmful landfill gas, which has an impact on the greenhouse effect. After the depositing of waste has been stopped, current conversion processes continue and emissions last for several decades and even up to 100years and longer. A good prediction of these processes is of high importance for landfill operators as well as for authorities, but suitable models for a realistic description of landfill processes are rather poor. In order to take the strong coupled conversion processes into account, a constitutive three-dimensional model based on the multiphase Theory of Porous Media (TPM) has been developed at the University of Duisburg-Essen. The theoretical formulations are implemented in the finite element code FEAP. With the presented calculation concept we are able to simulate the coupled processes that occur in an actual landfill. The models theoretical background and the results of the simulations as well as the meantime successfully performed simulation of a real landfill body will be shown in the following.

Ice Formation in Porous Media

Ice formation in porous media results from coupled heat and mass transport and is accompanied by ice expansion. The volume increase in space and time corresponds to the moving freezing front inside the porous solid. In this contribution, a macroscopic model based on the Theory of Porous Media (TPM) is presented toward the description of freezing and thawing processes in saturated porous media. Therefore, a quadruple model consisting of the constituents solid, ice, liquid and gas is used. Attention is paid to the description of capillary suction, liquid- and gas pressure on the surrounding surfaces, volume deformations due to ice formation, temperature distribution as well as influence of heat of fusion under thermal loading. For detection of energetic effects regarding the control of phase transition of water and ice, a physically motivated evolution equation for the mass exchange based on the local divergence of the heat flux is used. Numerical examples are presented to the applications of the model.

Characterization of methane oxidation in a simulated landfill cover system by comparing molecular and stable isotope mass balances

Biological methane oxidation may be regarded as a method of aftercare treatment for landfills to reduce climate relevant methane emissions. It is of social and economic interest to estimate the behavior of bacterial methane oxidation in aged landfill covers due to an adequate long-term treatment of the gas emissions. Different approaches assessing methane oxidation in laboratory column studies have been investigated by other authors recently. However, this work represents the first study in which three independent approaches, ((i) mass balance, (ii) stable isotope analysis, and (iii) stoichiometric balance of product (CO2) and reactant (CH4) by CO2/CH4-ratio) have been compared for the estimation of the biodegradation by a robust statistical validation on a rectangular, wide soil column. Additionally, an evaluation by thermal imaging as a potential technique for the localization of the active zone of bacterial methane oxidation has been addressed in connection with stable isotope analysis and CO2/CH4-ratios. Although landfills can be considered as open systems the results for stable isotope analysis based on a closed system correlated better with the mass balance than calculations based on an open system. CO2/CH4-ratios were also in good agreement with mass balance. In general, highest values for biodegradation were determined from mass balance, followed by CO2/CH4-ratio, and stable isotope analysis. The investigated topsoil proved to be very suitable as a potential cover layer by removing up to 99% of methane for CH4 loads of 35-65gm-2d-1 that are typical in the aftercare phase of landfills. Finally, data from stable isotope analysis and the CO2/CH4-ratios were used to trace microbial activity within the reactor system. It was shown that methane consumption and temperature increase, as a cause of high microbial activity, correlated very well.

Computational Modeling in Liver Surgery

The need for extended liver resection is increasing due to the growing incidence of liver tumors in aging societies. Individualized surgical planning is the key for identifying the optimal resection strategy and to minimize the risk of postoperative liver failure and tumor recurrence. Current computational tools provide virtual planning of liver resection by taking into account the spatial relationship between the tumor and the hepatic vascular trees, as well as the size of the future liver remnant. However, size and function of the liver are not necessarily equivalent. Hence, determining the future liver volume might misestimate the future liver function, especially in cases of hepatic comorbidities such as hepatic steatosis. A systems medicine approach could be applied, including biological, medical, and surgical aspects, by integrating all available anatomical and functional information of the individual patient. Such an approach holds promise for better prediction of postoperative liver function and hence improved risk assessment. This review provides an overview of mathematical models related to the liver and its function and explores their potential relevance for computational liver surgery. We first summarize key facts of hepatic anatomy, physiology, and pathology relevant for hepatic surgery, followed by a description of the computational tools currently used in liver surgical planning. Then we present selected state-of-the-art computational liver models potentially useful to support liver surgery. Finally, we discuss the main challenges that will need to be addressed when developing advanced computational planning tools in the context of liver surgery.