Modelling Drug Delivery and Tumour Growth
|
|
- Mercy Thompson
- 5 years ago
- Views:
Transcription
1 Modelling Drug Delivery and Tumour Growth Matthew Hubbard and Christian Groh, Pamela Jones, Brian Sleeman, Steven Smye, Christopher Twelves (Leeds) Paul Loadman, Nagaraja Periasamy, Roger Phillips (Bradford) Helen Byrne (Oxford) (Partly funded by Yorkshire Cancer Research) University of Greenwich, 11th December 214
2 PART ONE DRUG DELIVERY (TO TUMOURS)
3 Drug Delivery Tumours contain poorly organised and dysfunctional vasculature. Hence, they are poorly perfused with blood. This difference in microenvironment has a profound effect on tumour response to chemotherapy. A major obstacle to effective chemotherapy is inadequate delivery of drug to cells. Can we model the penetration of drug from blood vessels, through the surrounding layers of cells?
4 Drug Delivery Blood Vessels (red): Drug (blue): Hypoxia (green) Primeau, Rendon, Hedley, Lilge, Tannock (25)
5 Drug Delivery Realistic tumour vasculatures are difficult to obtain. They would also be expensive to simulate. Start by considering a single blood vessel and the surrounding layers of cells. This set-up resembles a tumour cord.
6 Tumour Cords BV = Blood Vessels, N = Necrosis, Glut1 Hypoxia
7 Cell Model k 1 k 2 C 1 C 2 C 3 with C 3 C V 1 V 2 extracellular space intracellular space k 2 Transport of drug across cell membrane (k 1 ) Binding and unbinding of drug in the cell (k 2,k 2 ) Potential saturation of binding sites (C ) Interstitial drug diffusion will link the cells (D)
8 Mathematical Model V 1 dc 1 dt V 2 dc 2 dt V 2 dc 3 dt = ak 1 (C 2 C 1 ) = ak 1 (C 1 C 2 ) V 2 k 2 C 2 (C C 3 ) + V 2 k 2 C 3 = V 2 k 2 C 2 (C C 3 ) V 2 k 2 C 3 Biochemical parameters: k 1 k 2 k 2 C Determined from in vitro experiments. Physical parameters: V 1 V 2 a Determined by observation of tissue.
9 Binding Experiments The in vitro experiments considered the action of the drug (Doxorubicin) on cells (DLD-1) in suspension. The initial extracellular drug concentration (C 1 ) was specified. Measurements were taken at later times of: extracellular drug concentration (C 1 ), intracellular drug concentration (C 2 + C 3 ). All the required parameters are identifiable from the data collected.
10 Parameter Sensitivity In fact, the model is extremely insensitive to the transmembrane transport rate (k 1 ) as long as it is fast enough.
11 Timescales The fit to the data suggests two distinct time-scales, which are confirmed by mathematical analysis and can be linked to the literature. 11 Free extracellular drug, C 1 1 Free intracellular drug, C 2 2 Bound intracellular drug, C Numerical simulation Full asymptotic approximation Fast-scale asymptotic approximation Concentration (µm) 7 Concentration (µm) 6 4 Concentration (µm) Numerical simulation Full asymptotic approximation Fast-scale asymptotic approximation 2 Numerical simulation Full asymptotic approximation Fast-scale asymptotic approximation Time (min) Time (min) Time (min)
12 Radially Symmetric Models A n A 2 A 1 A L L V n V 2 V 1 l V l d d Compartment Model Continuum Model
13 Radially Symmetric Models A compartment model. Each cell layer constitutes 3 compartments. Drug is exchanged between layers and compartments within layers. Effectively a finite volume approach with the mesh size dictated by the biological cell size. A continuum model. Reaction-diffusion partial differential equations. Discretised using a spectral approach. Adaptive time-stepping is used (ode15s in Matlab).
14 Compartment Model dc (i) 1 δ 1 V i dt dc (i) 2 δ 2 V i dt dc (i) 3 δ 2 V i dt = A i 1 k (C (i 1) 1 C (i) 1 ) + A i k (C (i+1) 1 C (i) 1 ) + a i k 1 (C (i) 2 C (i) 1 ) = a i k 1 (C (i) 1 C (i) 2 ) δ 2 V i k 2 C (i) 2 (C C (i) 3 ) + δ 2 V i k 2 C (i) 3 = δ 2 V i k 2 C (i) 2 (C C (i) 3 ) δ 2 V i k 2 C (i) 3 k = D/d governs transport between layers. δ 1 and δ 2 are volume fractions (δ 1 + δ 2 = 1).
15 Continuum Model δ 1 C 1 t δ 2 C 2 t δ 2 C 3 t = D ( 2 C 1 r r + αk 1 (C 2 C 1 ) C 1 r ) = αk 1 (C 1 C 2 ) δ 2 [k 2 C 2 (C C 3 ) k 2 C 3 ] = δ 2 [k 2 C 2 (C C 3 ) k 2 C 3 ] Both models have the boundary conditions D C 1 r = k v (C v C 1 r=l ) D C 1 r=l r = r=l
16 Cell-Centre Model This is a multidimensional version of the compartment model. Tumour cord Discrete representation
17 Cell-Centre Model dc (i) 1 δ 1 V i dt dc (i) 2 δ 2 V i dt dc (i) 3 δ 2 V i dt = j N i A ij d ij D ( + a (i) k 1 C (i) 2 C (i) 1 ( ) ( C (j) 1 C (i) 1 + j Vi A ij k v C v (j) ) ( ) = a (i) k 1 C (i) 1 C (i) 2 [ ( ) δ 2 V i k 2 C (i) 2 C C (i) 3 = δ 2 V i [ k 2 C (i) 2 ( ) C C (i) 3 ] k 2 C (i) 3 ] k 2 C (i) 3 ) C (i) 1
18 Numerical Experiments To investigate: the differences between the compartment, continuum and cell-centre models; the effect of the pharmacokinetic profile; the influence of binding affinities.
19 Pharmacokinetic Profiles For a bolus injection the concentration of drug in the vessel (c v (t)) is given by: D τ D τ [ A α ( 1 e αt ) + B ( β 1 e βt ) + C ( γ 1 e γt )] t < τ [ A α (eατ 1) e αt + B ( β e βτ 1 ) ] e βt + C γ (eγτ 1) e γt t τ Robert, Illiadis, Hoerni, Cano, Durand, Lagarde (1982) A simplified form is given by c v (t) = Ae αt. Infusion can be simulated by c v (t) = K. The parameters are chosen to give the same area under curve, i.e. c v dt.
20 Bolus Injection.25 Free extracellular drug, C1 (t = 1h).25 Free extracellular drug, C1 (t = 6h).25 Free extracellular drug, C1 (t = 24h).25 Free extracellular drug, C1 (t = 72h).2 Cell-centre Model Continuum Model Compartment Model.2 Cell-centre Model Continuum Model Compartment Model.2 Cell-centre Model Continuum Model Compartment Model.2 Cell-centre Model Continuum Model Compartment Model Concentration (µm).15.1 Concentration (µm).15.1 Concentration (µm).15.1 Concentration (µm) Distance from vessel centre (µm) Distance from vessel centre (µm) Distance from vessel centre (µm) Distance from vessel centre (µm) 1.5 Bound intracellular drug, C3 (t = 1h) 1.5 Bound intracellular drug, C3 (t = 6h) 1.5 Bound intracellular drug, C3 (t = 24h) 1.5 Bound intracellular drug, C3 (t = 72h) Cell-centre Model Cell-centre Model Cell-centre Model Cell-centre Model Continuum Model Continuum Model Continuum Model Continuum Model Compartment Model Compartment Model Compartment Model Compartment Model Concentration (µm).5 Concentration (µm).5 Concentration (µm).5 Concentration (µm) Distance from vessel centre (µm) Distance from vessel centre (µm) Distance from vessel centre (µm) Distance from vessel centre (µm) Drug distribution: extracellular (top), bound intracellular (bottom). Snapshots at times 1h, 6h, 24h, 72h (left to right).
21 Spatial Variation.2 Free extracellular drug, C 1.2 Free extracellular drug, C 1.2 Free extracellular drug, C t = 1 h t = 6 h t = 24 h t = 72 h t = 1 h t = 6 h t = 24 h t = 72 h t = 1 h t = 6 h t = 24 h t = 72 h Concentration (µm) Concentration (µm) Concentration (µm) Distance from vessel centre (µm) Distance from vessel centre (µm) Distance from vessel centre (µm) 2 Bound intracellular drug, C 3 2 Bound intracellular drug, C 3 2 Bound intracellular drug, C t = 1 h 1.8 t = 1 h 1.8 t = 1 h 1.6 t = 6 h 1.6 t = 6 h 1.6 t = 6 h 1.4 t = 24 h t = 72 h 1.4 t = 24 h t = 72 h 1.4 t = 24 h t = 72 h Concentration (µm) Concentration (µm) Concentration (µm) Distance from vessel centre (µm) Distance from vessel centre (µm) Distance from vessel centre (µm) Drug distribution: extracellular (top), bound intracellular (bottom); for bolus (left), exponential (middle), infusion (right).
22 Temporal Variation.12 Free extracellular drug, C1.12 Free extracellular drug, C1.12 Free extracellular drug, C1 r 26 µm r 26 µm r 26 µm.1 r 18 µm.1 r 18 µm.1 r 18 µm r 19 µm r 19 µm r 19 µm Concentration (µm).6.4 Concentration (µm).6.4 Concentration (µm) Time (h) Time (h) Time (h) 2 Bound intracellular drug, C3 2 Bound intracellular drug, C3 2 Bound intracellular drug, C Concentration (µm) Concentration (µm) Concentration (µm) r 26 µm r 18 µm.4 r 26 µm r 18 µm.4 r 26 µm r 18 µm.2 r 19 µm.2 r 19 µm.2 r 19 µm Time (h) Time (h) Time (h) Drug evolution: extracellular (top), bound intracellular (bottom); for bolus (left), exponential (middle), infusion (right).
23 Exposure 45 Bound intracellular drug, C 3 (t = 24h) 9 Bound intracellular drug, C 3 (t = 72h) 4 tri-exponential mono-exponential 8 35 uniform Exposure (µm h) 25 2 Exposure (µm h) tri-exponential mono-exponential uniform Distance from vessel centre (µm) Distance from vessel centre (µm) Exposure to bound intracellular drug, C 3 dt, after 24h (left) and 72h (right).
24 Varying Binding Affinity 25 Bound intracellular drug, C 3 (t = 72h) Exposure (µm h) k 2 =.4 M 1 s 1 k 2 = 4 M 1 s 1 k 2 = 4 M 1 s Distance from vessel centre (µm) Exposure to bound intracellular drug, C 3 dt, after 72h, for three different binding affinities.
25 Preliminary Observations The two discrete models and the discretised continuum model produce very similar results. The pharmacokinetic profile has less effect on the exposure than expected at long times. The exposure is sensitive to the binding rates. The model is extremely simple. The experimental measurements contain the largest source of variation.
26 The Next Steps Convective transport, validated by transwell experiments and possibly artificial tumour cords. More (parameterised?) processes and cell types. Cell cycle (mitosis, apoptosis, necrosis). Cell response (and a tumour growth model). Drug clearance mechanisms, irreversible binding. With cell movement, realistic vasculature and three dimensions the scientific computing contribution becomes more significant.
27 PART TWO TUMOUR GROWTH
28 Vascular Tumour Characteristics Rapid cell proliferation Imbalance between mitosis and apoptosis Angiogenesis and co-option of vessels Can be spherical or irregular shapes Necrotic regions in tumour core and in avascular areas
29 A Continuum Multiphase Model Based on Breward, Byrne and Lewis (22,23) Phases: cells, extracellular material, blood vessels Incompressible, viscous, fluids with drag between phases Cell-cell interactions included in stress tensors Includes mitosis, apoptosis/necrosis, angiogenesis, vessel occlusion A diffusible nutrient Tumour interacts with external tissue
30 Mass Balance The tissue consists of N p interacting phases, whose volume fractions θ i are governed by mass conservation θ i t + (θ i u i ) = q i i = 1,...,N p. u i are phase velocities q i allow conversion between phases Densities assumed to be constant and equal There are no voids, i.e. N p i=1 θ i = 1 and N p i=1 q i =.
31 A Four-Phase Model Four phases are chosen, to represent healthy cells (α 1 ) tumour cells (α 2 ) extracellular material (β) blood vessels (γ) More could be added, e.g. different tumour cell phenotypes, mature/immature vessels. The q i are chosen to ensure θ i [, 1].
32 Mass Sources For the cell phases: q αj ( ) c = k 1,j α j β c p + c }{{} cell birth ( ) cc1 + c k 2,j α j c c2 + c }{{} cell death c/(cp + c) c 2 Birth rate is half maximum at c = c p k are rate constants c are concentration thresholds (cc1 + c)/(cc2 + c) Necrosis rate is half maximum at c = c c2 Baseline apoptosis rate c
33 Mass Sources For the blood vessel phase: q γ = k 3 γ H(α 1 p α1 + α 2 p α2 p crit,ǫ 3 ) }{{} occlusion ( )( ) β c + k 4 (α 1 + α 2 )γ ε + β (c a + c) }{{ 2 } angiogenesis where H(θ,ǫ) is a smooth switch. For the remaining material: (1 + tanh(p/ǫ3))/2 c/(ca + c) p = p crit p = θ 1 p 1 + θ 2 p q β = q α1 q α2 q γ 1.5 Rate is maximum at c = c a c
34 Momentum Balance Assuming low Reynolds number flow, (θ i σ i ) + F i = i = 1,...,N p. σ i are the phase stress tensors. F i are momentum sources and sinks. Note that mass conservation implies N p i=1 (θ i u i ) =.
35 Constitutive Laws For viscous Newtonian fluids σ i = p i I + µ i ( u i + u T i ) + λ i ( u i )I i = 1,...,N p, though λ i = 2 3 µ i is assumed (local thermodynamic equilibrium). Momentum sources include pressure effects and interphase drag (with coefficients d ij ): F i = p i I θ i + j i d ij θ i θ j ( u j u i ) i = 1,...,N p.
36 Pressure Relations Pressure due to the blood vessel phase is assumed independent of the local flow, i.e. p γ = p γ. Cells act like a fluid with additional cell-cell interactions, i.e. where p α1 = p α2 = p β + Σ α Σ α = max ( ) Λ(α α ) (1 α), 2 with α = α 1 + α 2 and α a natural equilibrium density. max(,(θ θ )/(1 θ) 2 ) θ = θ θ
37 Diffusible Species There are N d additional species whose concentrations c j are governed by D j 2 c j + q j = j = 1,...,N d. They are assumed to evolve much more rapidly than the tissue. They do not contribute to the volume. They provide non-local influence.
38 Diffusion Sources A single species representing nutrient is used here, for which q c = k 5 γ (c v c) }{{} replenishment k 6,1 α 1 c k 6,2 α 2 c }{{} baseline consumption ( ) ( ) c c k 7,1 α 1 β k 7,2 α 2 β c p + c c p + c }{{} consumption due to cell birth It would be simple to include others, e.g. VEGF, ph, glucose.
39 The Whole System θ i t + (θ i u i ) = q i N p i=1 θ i = 1 (θ i [µ i ( u i + u T i ) + λ i ( u i )I]) + j i d ij θ i θ j ( u j u i ) θ i (pi I) = N p (θ i u i ) = i=1 D j 2 c j + q j =
40 Boundary Conditions Mass balance: specify phase volume fractions at inflow. Momentum balance: u β = somewhere on the boundary; zero stress elsewhere. Diffusible species: c n = (zero flux).
41 The Whole System The unknowns are θ i phase volume fractions u i phase velocities p i phase pressures c j diffusible species concentrations Using unstructured triangular meshes in two space dimensions, evolve the phase volume fractions in time then update the other variables.
42 Numerical Approximation Phase volume fractions Hyperbolic conservation laws with sources/sinks Important to remain within [, 1] Cell-centre finite volume scheme on refined mesh Upwind with MUSCL slope limiting Forward Euler time-stepping
43 Numerical Approximation Phase velocities and pressures A form of multicomponent Stokes flow Stable finite elements quadratic velocities, linear pressures System is linear but not (quite) symmetric MUMPS (without the MP) Care is needed when θ i =
44 Numerical Approximation Nutrient concentration Standard reaction-diffusion system Galerkin finite elements with linear elements on refined mesh Newton iteration combined with MUMPS Allows nonlinear reaction terms Jacobian approximated numerically
45 Test Case Circular domain, radius 16 Triangular mesh: 2349 nodes, 4539 elements Gives DoFs in FEM solve Free boundaries: zero normal stress except β Seed healthy tissue in centre with α 2 =.5 cos 2 (πr/2) for r 1 otherwise (animations)
46 Phase Volume Fractions Tumour cells: θ2 at t = 3 Healthy cells: θ1 at t = 3 Cell pressure: p1,p2 at t = Extracellular material: θ4 at t = 3 Blood vessels: θ3 at t = 3 2 ECM pressure: p4 at t = p. 46/5
47 Phase Fluxes Healthy cells: θ 1 u 1 at t = 3 Tumour cells: θ 2 u 2 at t = 3 Extracellular material: θ 4 u 4 at t = 3
48 Varying Diffusion Healthy cells: θ1 at t = 225 Healthy cells: θ1 at t = 525 Healthy cells: θ1 at t = Tumour cells: θ2 at t = Tumour cells: θ2 at t = 8 Tumour cells: θ2 at t = p. 48/5
49 Varying Diffusion Extracellular material: θ4 at t = 225 Extracellular material: θ4 at t = 525 Extracellular material: θ4 at t = Blood vessels: θ3 at t = 225 Blood vessels: θ3 at t = 525 Blood vessels: θ3 at t = p. 49/5
50 Parameter Dependence Rmax(t) Rmax(t) Rmax(t) 6 3 D c =.1 D c =.1 D c =1. D c = Λ=.1 Λ=.1 Λ=.1 Λ=1. Λ= µ=1. µ=1. µ=1. µ= Time, t Time, t Time, t Nutrient diffusion rate (left), cell tension (middle), phase viscosity (right).
51 Front Irregularity Tumour cells: θ 2 at t = Tumour cells: θ 2 at t = Variable phase viscosities: low (left), high (right).
52 Model Behaviour The proliferating rim/necrotic core structure is reproduced. Shapes can be almost circular (in 2D) or irregular. The regularity is most sensitive to the phase viscosities and the cell tension constant. Tumour growth is approximately linear. The growth rate is increased by increasing nutrient diffusion rates, cell birth rates or cell tension, or decreasing phase viscosities or cell death rates.
53 Issues Speed and memory Iterative solvers, preconditioning, multigrid Better efficiency needed for 3D, more phases, etc. Physical/biological realism Measurable model parameters Numerical artefacts, resolution issues Adaptive approaches?
54 Current Developments Multiscale asymptotic analysis is allowing us to scale the effects of the capillary-level drug delivery model up to the scale of a full tumour. The large-scale drug delivery model will be linked with a tumour response model. The response model will investigate the combined effects of chemotherapy and radiotherapy. The ultimate aim is to devise optimal, personalised, treatment strategies.
Simulating Solid Tumor Growth Using Multigrid Algorithms
Simulating Solid Tumor Growth Using Multigrid Algorithms Asia Wyatt Applied Mathematics, Statistics, and Scientific Computation Program Advisor: Doron Levy Department of Mathematics/CSCAMM Abstract In
More informationBiotransport: Principles
Robert J. Roselli Kenneth R. Diller Biotransport: Principles and Applications 4 i Springer Contents Part I Fundamentals of How People Learn (HPL) 1 Introduction to HPL Methodology 3 1.1 Introduction 3
More informationMODELING THE ABSORPTION, CIRCULATION, AND METABOLISM OF TIRAPAZAMINE
MODELING THE ABSORPTION, CIRCULATION, AND METABOLISM OF TIRAPAZAMINE Olivia Nguyen and Xiaoyu Shi Problem Statement: Tirapazamine (TPZ; 1,2,4-benzotriazin-3-amine 1,4-dioxide) is a hypoxia-activated prodrug
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationSolving a non-linear partial differential equation for the simulation of tumour oxygenation
Solving a non-linear partial differential equation for the simulation of tumour oxygenation Julian Köllermeier, Lisa Kusch, Thorsten Lajewski MathCCES, RWTH Aachen Lunch Talk, 26.05.2011 J. Köllermeier,
More informationSolving a non-linear partial differential equation for the simulation of tumour oxygenation
Solving a non-linear partial differential equation for the simulation of tumour oxygenation Julian Köllermeier, Lisa Kusch, Thorsten Lajewski MathCCES, RWTH Aachen Talk at Karolinska Institute, 19.09.2011
More informationCode Verification of Multiphase Flow with MFIX
Code Verification of Multiphase Flow with MFIX A N I R U D D H A C H O U D H A R Y ( A N I R U D D @ V T. E D U ), C H R I S T O P H E R J. ROY ( C J R O Y @ V T. E D U ), ( V I R G I N I A T E C H ) 23
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationTABLE OF CONTENTS CHAPTER TITLE PAGE
v TABLE OF CONTENTS CHAPTER TITLE PAGE TABLE OF CONTENTS LIST OF TABLES LIST OF FIGURES LIST OF SYMBOLS LIST OF APPENDICES v viii ix xii xiv CHAPTER 1 INTRODUCTION 1.1 Introduction 1 1.2 Literature Review
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationChemotherapy in heterogeneous cultures of cancer cells with interconversion
Chemotherapy in heterogeneous cultures of cancer cells with interconversion Rui DILãO Institut des Hautes Études Scientifiques 35, route de Chartres 91440 Bures-sur-Yvette (France) Juin 2013 IHES/M/13/20
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationSimulating Solid Tumor Growth using Multigrid Algorithms
Simulating Solid Tumor Growth using Multigrid Applied Mathematics, Statistics, and Program Advisor: Doron Levy, PhD Department of Mathematics/CSCAMM University of Maryland, College Park September 30, 2014
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationInverse problems in tumor growth modelling
MC2 Damiano Lombardi Angelo Iollo Thierry Colin Olivier Saut Marseille, 18-09-09 Outline 1) Introduction: phenomenological model for tumor growth models 2) Models based on mixture theory: a simplified
More informationLevel Set Tumor Growth Model
Level Set Tumor Growth Model Andrew Nordquist and Rakesh Ranjan, PhD University of Texas, San Antonio July 29, 2013 Andrew Nordquist and Rakesh Ranjan, PhD (University Level Set of Texas, TumorSan Growth
More informationIntroduction to some topics in Mathematical Oncology
Introduction to some topics in Mathematical Oncology Franco Flandoli, University of Pisa y, Finance and Physics, Berlin 2014 The field received considerable attentions in the past 10 years One of the plenary
More informationUnsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe
Unsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe T S L Radhika**, M B Srinivas, T Raja Rani*, A. Karthik BITS Pilani- Hyderabad campus, Hyderabad, Telangana, India. *MTC, Muscat,
More informationStoichiometry Of Tumor Dynamics: Models And Analysis
Stoichiometry Of Tumor Dynamics: Models And Analysis Yang Kuang Department of Mathematics and Statistics Arizona State University Supported by NSF grant DMS-0077790 Y. Kuang, J. Nagy and J. Elser: Disc.
More informationA. Kovacevic N. Stosic I. Smith. Screw Compressors. Three Dimensional Computational Fluid Dynamics and Solid Fluid Interaction.
Screw Compressors A. Kovacevic N. Stosic I. Smith Screw Compressors Three Dimensional Computational Fluid Dynamics and Solid Fluid Interaction With 83 Figures Ahmed Kovacevic Nikola Stosic Ian Smith School
More informationProblem Set No. 4 Due: Monday, 11/18/10 at the start of class
Department of Chemical Engineering ChE 170 University of California, Santa Barbara Fall 2010 Problem Set No. 4 Due: Monday, 11/18/10 at the start of class Objective: To understand the thermodynamic and
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationNumerical Simulation of the Hagemann Entrainment Experiments
CCC Annual Report UIUC, August 14, 2013 Numerical Simulation of the Hagemann Entrainment Experiments Kenneth Swartz (BSME Student) Lance C. Hibbeler (Ph.D. Student) Department of Mechanical Science & Engineering
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationA Multigrid LCR-FEM solver for viscoelastic fluids with application to problems with free surface
A Multigrid LCR-FEM solver for viscoelastic fluids with application to problems with free surface Damanik, H., Mierka, O., Ouazzi, A., Turek, S. (Lausanne, August 2013) Page 1 Motivation Polymer melts:
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationScientific Computing I
Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationTaxis Diffusion Reaction Systems
Chapter 2 Taxis Diffusion Reaction Systems In this chapter we define the class of taxis diffusion reaction (TDR) systems which are the subject of this thesis. To this end, we start in Sec. 2. with the
More informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationBiophysics I. DIFFUSION
Biophysics I. DIFFUSION Experiment add a droplet of ink to a glass of water Observation: the stain spreads and eventually colours the entire fluid add a droplet of ink to HOT and COLD water Observation:
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationCompartmental Modelling: Eigenvalues, Traps and Nonlinear Models
Compartmental Modelling: Eigenvalues, Traps and Nonlinear Models Neil D. Evans neil.evans@warwick.ac.uk School of Engineering, University of Warwick Systems Pharmacology - Pharmacokinetics Vacation School,
More informationNumerical Model of the Influence of Shear Stress on the Adaptation of a Blood Vessel BMT 03-35
Numerical Model of the Influence of Shear Stress on the Adaptation of a Blood Vessel BMT 03-35 Mirjam Yvonne van Leeuwen Supervisor: Dr. Ir. M.C.M. Rutten Ir. N.J.B. Driessen TUE Eindhoven, The Netherlands
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More information2. Modeling of shrinkage during first drying period
2. Modeling of shrinkage during first drying period In this chapter we propose and develop a mathematical model of to describe nonuniform shrinkage of porous medium during drying starting with several
More informationHow are calculus, alchemy and forging coins related?
BMOLE 452-689 Transport Chapter 8. Transport in Porous Media Text Book: Transport Phenomena in Biological Systems Authors: Truskey, Yuan, Katz Focus on what is presented in class and problems Dr. Corey
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
More informationNANO 243/CENG 207 Course Use Only
L8. Drug Dispersion and Diffusion in Biological Systems (2) April 26, 2018 2. Diffusion in Water (3) Diffusion of small molecule in water drug molecules/solvents are identical to solvent Assume molecules
More informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
More informationThe Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011
The Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011 Derivation of Navier-Stokes Equation 1 The total stress tensor
More informationThe Navier-Stokes Equations
s University of New Hampshire February 22, 202 and equations describe the non-relativistic time evolution of mass and momentum in fluid substances. mass density field: ρ = ρ(t, x, y, z) velocity field:
More informationThe Polymers Tug Back
Tugging at Polymers in Turbulent Flow The Polymers Tug Back Jean-Luc Thiffeault http://plasma.ap.columbia.edu/ jeanluc Department of Applied Physics and Applied Mathematics Columbia University Tugging
More informationSolving a non-linear partial differential equation for the simulation of tumour oxygenation
Solving a non-linear partial differential equation for the simulation of tumour oxygenation Julian Köllermeier, Lisa Kusch, Thorsten Lajewski Contents Contents 3 1 Introduction 5 2 Application 7 2.1 Tumour
More informationA Mathematical Model for Tumor Growth and Angiogenesis
A Mathematical Model for Tumor Growth and Angiogenesis Wietse Boon July 11, 2011 1 Contents 1 Introduction 3 2 Angiogenesis 4 2.1 Discretization......................................... 5 2.2 Concentration
More informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Mats Nigam Copyright 2016 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their respective
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationNumerical study of blood fluid rheology in the abdominal aorta
Design and Nature IV 169 Numerical study of blood fluid rheology in the abdominal aorta F. Carneiro 1, V. Gama Ribeiro 2, J. C. F. Teixeira 1 & S. F. C. F. Teixeira 3 1 Universidade do Minho, Departamento
More informationSimulating Microbubble Flows Using COMSOL Mutiphysics
Presented at the COMSOL Conference 2008 Boston Simulating Microbubble Flows Using COMSOL Mutiphysics Xiaodong(Sheldon) Chen 1, Samir N. Ghadiali *2 1 Department of Mechanical Engineering, The Ohio State
More informationApplication of the immersed boundary method to simulate flows inside and outside the nozzles
Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,
More informationNumerical Simulation of Powder Flow
Numerical Simulation of Powder Flow Stefan Turek, Abderrahim Ouazzi Institut für Angewandte Mathematik, Univ. Dortmund http://www.mathematik.uni-dortmund.de/ls3 http://www.featflow.de Models for granular
More informationfluid mechanics as a prominent discipline of application for numerical
1. fluid mechanics as a prominent discipline of application for numerical simulations: experimental fluid mechanics: wind tunnel studies, laser Doppler anemometry, hot wire techniques,... theoretical fluid
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationUniversidad Nacional de Córdoba CIEM-CONICET, Argentina
Selçuk J. Appl. Math. Vol. 10. No. 1. pp. 147-155, 2009 Selçuk Journal of Applied Mathematics A Simple Discrete Model for the Growth Tumor Andrés Barrea, Cristina Turner Universidad Nacional de Córdoba
More informationComputer Modeling in Bioengineering
Computer Modeling in Bioengineering Theoretical Background, Examples and Software Milos Kojic Harvard School of Public Health, USA University of Kragujevac, Serbia University of Texas Health Science Center
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationIntroduction to stochastic multiscale modelling of tumour growth
Introduction to stochastic multiscale modelling of tumour growth Tomás Alarcón ICREA & Centre de Recerca Matemàtica T. Alarcón (ICREA & CRM, Barcelona, Spain) Lecture 1 CIMPA, Santiago de Cuba, June 2016
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationAGITATION AND AERATION
AGITATION AND AERATION Although in many aerobic cultures, gas sparging provides the method for both mixing and aeration - it is important that these two aspects of fermenter design be considered separately.
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
More informationSome notes about PDEs. -Bill Green Nov. 2015
Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can
More informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
More informationBiology in Time and Space: A Partial Differential Equation Approach. James P. Keener University of Utah
Biology in Time and Space: A Partial Differential Equation Approach James P. Keener University of Utah 26 December, 2015 2 Contents I The Text 9 1 Introduction 11 2 Conservation 13 2.1 The Conservation
More informationChristel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011
Notes by: Andy Thaler Christel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011 Many complex fluids are shear-thinning. Such a fluid has a shear
More informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationMagnetic Drug Targeting in Cancer Therapy
Magnetic Drug Targeting in Cancer Therapy SOLVED WITH COMSOL MULTIPHYSICS 3.5a COPYRIGHT 2008. All right reserved. No part of this documentation may be photocopied or reproduced in any form without prior
More informationProject 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions
Project 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions R. S. Sellers MAE 5440, Computational Fluid Dynamics Utah State University, Department of Mechanical and Aerospace Engineering
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationTumour angiogenesis as a chemo-mechanical surface instability. Supplementary Information Chiara Giverso 1, Pasquale Ciarletta 2
Tumour angiogenesis as a chemo-mechanical surface instability 1 Theoretical Derivation Supplementary Information Chiara Giverso 1, Pasquale Ciarletta 2 In order to close the description of the process
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationComparison of Heat and Mass Transport at the Micro-Scale
Comparison of Heat and Mass Transport at the Micro-Scale E. Holzbecher, S. Oehlmann Georg-August Univ. Göttingen *Goldschmidtstr. 3, 37077 Göttingen, GERMANY, eholzbe@gwdg.de Abstract: Phenomena of heat
More informationReliability of LES in complex applications
Reliability of LES in complex applications Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) DESIDER Symposium Corfu, June 7-8, 27 Sample of complex flow
More informationEPSRC Centre for Doctoral Training in Industrially Focused Mathematical Modelling
EPSRC Centre for Doctoral Training in Industrially Focused Mathematical Modelling Penetration of a liquid agent into a polymer Valentin Sulzer Contents 1. Introduction... 2 Background... 2 Modelling approach...
More informationDigital Simulation for the Behavior of the Flow of Non-Newtonian Fluids in 90 Pipe Bend
International Journal of Engineering and Technical Research (IJETR) ISSN: 2321-0869 (O) 2454-4698 (P), Volume-3, Issue-8, August 2015 Digital Simulation for the Behavior of the Flow of Non-Newtonian Fluids
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationA Mixed Finite Element Formulation for Solving Phase Change Problems with Convection
A Mixed Finite Element Formulation for Solving Phase Change Problems with Convection Youssef Belhamadia 1, Abdoulaye S. Kane 2, and André Fortin 3 1 University of Alberta, Campus Saint-Jean and Department
More informationA Finite-Element based Navier-Stokes Solver for LES
A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,
More informationA Computational Model of Vortex Ring Formation in Hydrogels
A Computational Model of Vortex Ring Formation in Hydrogels BEE/MAE 4530 Computer Aided Engineering: Applications to Biomedical Processes Shoshana Das Alyssa Fiore Dillon Hickman Varsha Vasthare Table
More informationA surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant
A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant Ashley J. James Department of Aerospace Engineering and Mechanics, University of Minnesota John Lowengrub
More informationAbsorption of gas by a falling liquid film
Absorption of gas by a falling liquid film Christoph Albert Dieter Bothe Mathematical Modeling and Analysis Center of Smart Interfaces/ IRTG 1529 Darmstadt University of Technology 4th Japanese-German
More informationTurbulence - Theory and Modelling GROUP-STUDIES:
Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence
More informationFluid Dynamics and Balance Equations for Reacting Flows
Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and
More information