Modelling oedemous limbs and venous ulcers using partial differential equations

Background Oedema, commonly known as tissue swelling, occurs mainly on the leg and the arm. The condition may be associated with a range of causes such as venous diseases, trauma, infection, joint disease and orthopaedic surgery. Oedema is caused by both lymphatic and chronic venous insufficiency, which leads to pooling of blood and fluid in the extremities. This results in swelling, mild redness and scaling of the skin, all of which can culminate in ulceration. Methods We present a method to model a wide variety of geometries of limbs affected by oedema and venous ulcers. The shape modelling is based on the PDE method where a set of boundary curves are extracted from 3D scan data and are utilised as boundary conditions to solve a PDE, which provides the geometry of an affected limb. For this work we utilise a mixture of fourth order and sixth order PDEs, the solutions of which enable us to obtain a good representative shape of the limb and associated ulcers in question. Results A series of examples are discussed demonstrating the capability of the method to produce good representative shapes of limbs by utilising a series of curves extracted from the scan data. In particular we show how the method could be used to model the shape of an arm and a leg with an associated ulcer. Conclusion We show how PDE based shape modelling techniques can be utilised to generate a variety of limb shapes and associated ulcers by means of a series of curves extracted from scan data. We also discuss how the method could be used to manipulate a generic shape of a limb and an associated wound so that the model could be fine-tuned for a particular patient

3D mesh processing using GAMer 2 to enable reaction-diffusion simulations in realistic cellular geometries

Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer simulations of such processes using numerical techniques such as the Finite Element Method. In this paper, we describe the use of our recently rewritten mesh processing software, GAMer 2, to bridge the gap between poorly conditioned meshes generated from segmented micrographs and boundary marked tetrahedral meshes which are compatible with simulation. We demonstrate the application of a workflow using GAMer 2 to a series of electron micrographs of neuronal dendrite morphology explored at three different length scales and show that the resulting meshes are suitable for finite element simulations. This work is an important step towards making physical simulations of biological processes in realistic geometries routine. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery at the interface of geometry and cellular processes. We posit that a new frontier at the intersection of computational technologies and single cell biology is now open.Comment: 39 pages, 14 figures. High resolution figures and supplemental movies available upon reques

Advances in reduced order methods for parametric industrial problems in computational fluid dynamics

Reduced order modeling has gained considerable attention in recent decades owing to the advantages offered in reduced computational times and multiple solutions for parametric problems. The focus of this manuscript is the application of model order reduction techniques in various engineering and scientific applications including but not limited to mechanical, naval and aeronautical engineering. The focus here is kept limited to computational fluid mechanics and related applications. The advances in the reduced order modeling with proper orthogonal decomposition and reduced basis method are presented as well as a brief discussion of dynamic mode decomposition and also some present advances in the parameter space reduction. Here, an overview of the challenges faced and possible solutions are presented with examples from various problems

Multilevel Solvers for Unstructured Surface Meshes

Parameterization of unstructured surface meshes is of fundamental importance in many applications of digital geometry processing. Such parameterization approaches give rise to large and exceedingly ill-conditioned systems which are difficult or impossible to solve without the use of sophisticated multilevel preconditioning strategies. Since the underlying meshes are very fine to begin with, such multilevel preconditioners require mesh coarsening to build an appropriate hierarchy. In this paper we consider several strategies for the construction of hierarchies using ideas from mesh simplification algorithms used in the computer graphics literature. We introduce two novel hierarchy construction schemes and demonstrate their superior performance when used in conjunction with a multigrid preconditioner