41 research outputs found
An unfitted discontinuous Galerkin scheme for conservation laws on evolving surfaces
Motivated by considering partial differential equations arising from conservation laws posed on evolving surfaces, a new numerical method for an advection problem is developed and simple numerical tests are performed. The method is based on an unfitted discontinuous Galerkin approach where the surface is not explicitly tracked by the mesh which means the method is extremely flexible with respect to geometry. Furthermore, the discontinuous Galerkin approach is well-suited to capture the advection driven by the evolution of the surface without the need for a space-time formulation, back-tracking trajectories or streamline diffusion. The method is illustrated by a one-dimensional example and numerical results are presented that show good convergence properties for a simple test problem
The interface for functions in the dune-functions module
The dune-functions dune module introduces a new programmer interface for discrete and non-discrete functions. Unlike the previous interfaces considered in the existing dune modules, it is based on overloading operator(), and returning values by-value. This makes user code much more readable, and allows the incorporation of newer C++ features such as lambda expressions. Run-time polymorphism is implemented not by inheritance, but by type erasure, generalizing the ideas of the std::function class from the C++11 standard library. We describe the new interface, show its possibilities, and measure the performance impact of type erasure and return-by-value
Function space bases in the dune-functions module
The dune-functions Dune module provides interfaces for functions and function space bases. It forms one abstraction level above grids, shape functions, and linear algebra, and provides infrastructure for full discretization frameworks like dune-pdelab and dune-fem. This document describes the function space bases provided by dune-functions. These are based on an abstract description of bases for product spaces as trees of simpler bases. From this description, many different numberings of degrees of freedom by multi-indices can be derived in a natural way. We describe the abstract concepts, document the programmer interface, and give a complete example program that solves the stationary Stokes equation using Taylor-Hood elements
Co-assembly of chitosan and phospholipids into hybrid hydrogels
Novel hybrid hydrogels were formed by adding chitosan (Ch) to phospholipids (P) self-assembled particles in lactic acid. The effect of the phospholipid concentration on the hydrogel properties was investigated and was observed to affect the rate of hydrogel formation and viscoelastic properties. A lower concentration of phospholipids (0.5% wt/v) in the mixture, facilitates faster network formation as observed by Dynamic Light Scattering, with lower elastic modulus than the hydrogels formed with higher phospholipid content. The nano-porous structure of Ch/P hydrogels, with a diameter of 260±20 nm, as observed by cryo-scanning electron microscopy, facilitated the penetration of water and swelling. Cell studies revealed suitable biocompatibility of the Ch/P hydrogels that can be used within life sciences applications
The Dune framework: Basic concepts and recent developments
This paper presents the basic concepts and the module structure of the Distributed and Unified Numerics Environment and reflects on recent developments and general changes that happened since the release of the first Dune version in 2007 and the main papers describing that state Bastian etal. (2008a, 2008b). This discussion is accompanied with a description of various advanced features, such as coupling of domains and cut cells, grid modifications such as adaptation and moving domains, high order discretizations and node level performance, non-smooth multigrid methods, and multiscale methods. A brief discussion on current and future development directions of the framework concludes the paper
Cell-scale degradation of peritumoural extracellular matrix fibre network and its role within tissue-scale cancer invasion
Local cancer invasion of tissue is a complex, multiscale process which plays
an essential role in tumour progression. Occurring over many different temporal
and spatial scales, the first stage of invasion is the secretion of matrix
degrading enzymes (MDEs) by the cancer cells that consequently degrade the
surrounding extracellular matrix (ECM). This process is vital for creating
space in which the cancer cells can progress and it is driven by the activities
of specific matrix metalloproteinases (MMPs). In this paper, we consider the
key role of two MMPs by developing further the novel two-part multiscale model
introduced in [33] to better relate at micro-scale the two micro-scale
activities that were considered there, namely, the micro-dynamics concerning
the continuous rearrangement of the naturally oriented ECM fibres within the
bulk of the tumour and MDEs proteolytic micro-dynamics that take place in an
appropriate cell-scale neighbourhood of the tumour boundary. Focussing
primarily on the activities of the membrane-tethered MT1-MMP and the soluble
MMP-2 with the fibrous ECM phase, in this work we investigate the MT1-MMP/MMP-2
cascade and its overall effect on tumour progression. To that end, we will
propose a new multiscale modelling framework by considering the degradation of
the ECM fibres not only to take place at macro-scale in the bulk of the tumour
but also explicitly in the micro-scale neighbourhood of the tumour interface as
a consequence of the interactions with molecular fluxes of MDEs that exercise
their spatial dynamics at the invasive edge of the tumour
Toward transient finite element simulation of thermal deformation of machine tools in real-time
Finite element models without simplifying assumptions can accurately describe the spatial and temporal distribution of heat in machine tools as well as the resulting deformation. In principle, this allows to correct for displacements of the Tool Centre Point and enables high precision manufacturing. However, the computational cost of FE models and restriction to generic algorithms in commercial tools like ANSYS prevents their operational use since simulations have to run faster than real-time. For the case where heat diffusion is slow compared to machine movement, we introduce a tailored implicit–explicit multi-rate time stepping method of higher order based on spectral deferred corrections. Using the open-source FEM library DUNE, we show that fully coupled simulations of the temperature field are possible in real-time for a machine consisting of a stock sliding up and down on rails attached to a stand
Structured models of cell migration incorporating molecular binding processes
The dynamic interplay between collective cell movement and the various
molecules involved in the accompanying cell signalling mechanisms plays a
crucial role in many biological processes including normal tissue development
and pathological scenarios such as wound healing and cancer. Information about
the various structures embedded within these processes allows a detailed
exploration of the binding of molecular species to cell-surface receptors
within the evolving cell population. In this paper we establish a general
spatio-temporal-structural framework that enables the description of molecular
binding to cell membranes coupled with the cell population dynamics. We first
provide a general theoretical description for this approach and then illustrate
it with two examples arising from cancer invasion