NAJan 25, 2013
Pressure jump interface law for the Stokes-Darcy coupling: Confirmation by direct numerical simulationsThomas Carraro, Christian Goll, Anna Marciniak-Czochra et al.
It is generally accepted that the effective velocity of a viscous flow over a porous bed satisfies the Beavers-Joseph slip law. To the contrary, interface law for the effective stress has been a subject of controversy. Recently, a pressure jump interface law has been rigorously derived by Marciniak-Czochra and Mikelić. In this paper, we provide a confirmation of the analytical result using direct numerical simulation of the flow at the microscopic level.
NAMar 9, 2016
Coupling vs decoupling approaches for PDE/ODE systems modeling intercellular signalingThomas Carraro, Elfriede Friedmann, Daniel Gerecht
We consider PDE/ODE systems for the simulation of intercellular signaling in multicellular environments. The intracellular processes for each cell described here by ODEs determine the long-time dynamics, but the PDE part dominates the solving effort. Thus, it is not clear if commonly used decoupling methods can outperform a coupling approach. Based on a sensitivity analysis, we present a systematic comparison between coupling and decoupling approaches for this class of problems and show numerical results. For biologically relevant configurations of the model, our quantitative study shows that a coupling approach performs much better than a decoupling one.
NAMay 19, 2018
A level-set approach for a multi-scale cancer invasion modelThomas Carraro, Sven E. Wetterauer, Ana Victoria Ponce Bobadilla et al.
Central to the quest for a deeper understanding of the cancer growth and spread process, the naturally multiscale character of cancer invasion demands appropriate multiscale modelling and analysis approach. The cross-talk between the tissue scale (macro-scale) cancer cell population dynamics and the cell-scale (micro-scale) proteolytic molecular processes along the tumour boundary plays a particularly important role within the invasion processes, leading to dramatic changes in tumour morphology and influencing the overall pattern of cancer spread. Building on the multiscale moving boundary framework proposed in Trucu et al. (Multiscale Model. Simul 11(1): 309-335), in this work we propose a new formulation of this process involving a novel derivation of the macro-scale boundary movement law based on micro-dynamics, involving a transport equation combined with the level-set method. This is explored numerically in a novel finite element macro-micro framework based on cut-cells.
NAJul 15, 2015
On the implementation of the eXtended Finite Element Method (XFEM) for interface problemsThomas Carraro, Sven Wetterauer
The eXtended Finite Element Method (XFEM) is used to solve interface problems with an unfitted mesh. We present an implementation of the XFEM in the FEM-library deal.II. The main parts of the implementation are (i) the appropriate quadrature rule; (ii) the shape functions for the extended part of the finite element formulation; (iii) the boundary and interface conditions. We show how to handle the XFEM formulation providing a code that demonstrates the solution of two exemplary interface problems for a strong and a weak discontinuity respectively. In the weak discontinuity case, the loss of conformity due to the blending effect and its remedy are discussed. Furthermore, the optimal convergence of the presented unfitted method is numerically verified.