Numerical method for osmotic water flow and solute diffusion in moving cells
|
|
- Isaac Jennings
- 5 years ago
- Views:
Transcription
1 Numerical method for osmotic water flow and solute diffusion in moving cells Lingxing Yao and Yoichiro Mori Case Western Reserve University and University of Minnesota IMA Workshop Electrohydrodynamics and Electrodiffusion in Material Sciences and Biology March 15, 018 Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
2 Cell structure and movement Torsten Wittmann, Scripps Research Institute. Actin filiments (purple), microtubules (yellow), and nuclei (green) 3 A. Schwab, A. Fabian, P. J. Hanley, and C. Stock, Physiol. Rev., 9, , (01) Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, 018 / 5
3 Osmotic engine experiment Cancer cell can move through confined extracellular matrix, including longitudinal tracks of varying length, with actin polymerization and depolymerization inhibited Osmotic shock, or polarized distribution of ion channels and aquaporins can lead to sustained cell migration in confined collagen matrix For cancer cells, actin polymerization cycle and water permeation due to osmosis may be the two mechanisms to drive the migration 1 1 K. M. Stroka, H. Jian, S. Chen, Z. Tong, D. Wirtz, S. Sun, and K. Konstantopoulos, Cell, 157, , (014) Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
4 Our goal and model assumptions We want to verify the role of chemical osmosis in driving cell migration by studying the diffusion of chemicals, osmotic water flow, and actin network, in a model cell system, which includes: Cell membrane is an elastic structure dividing the intra- and extra-cellular space and allows water and chemicals to pass through Chemicals diffuse inside and outside the cell, and cross cell membrane via passive or active channels Difference in solute concentrations across the membrane, together with mechanical forces from the membrane, lead to movement of the membrane and transmembrane water flow Actin network is present over the entire intra-cellular space as another phase and is able to polymerize/depolymerize (currently ongoing with some results) Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
5 Model equation part 1: fluid flow 1 The variables in this module include fluid velocity u, pressure p. We have incompressible Stokes (or Navier-Stokes) equation: u t + u u = Σ m(u, p), u = 0, where ν is fluid viscosity, and viscous stress Σ m (u, p) ν( u + ( u) T ) pi. Fluid velocity is continuous across membrane: [u] = 0 The membrane may not be moving with fluid velocity: u X t = f w n, with water flux f w = k w [ψ], and water potential ψ = c (Σ m (u, p)n) n. Here c is the chemical concentration in fluid. 1 L. Yao and Y. Mori, A numerical method for osmotic water flow and solute diffusion with deformable membrane boundaries in two spatial dimension, Journal of Computational Physics, 350, , 017 Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
6 Model equation part : membrane elasticity and force The membrane is denoted by X(s) (Lagrangian description of the membrane), which has elastic forces: (( ) ) X F elas = k elas s s l τ, τ = X 1 X s s, F bend = k bend 4 X s 4, Along cell membrane, viscous stress difference is balanced by the mechanical force of the membrane (n is outward normal along the membrane): [Σ m (u, p)n] = F mem X s 1, F mem = F elas + F bend Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
7 Model equation part 3: chemical concentration The chemical is defined over entire computational domain, with concentration c: c + (uc) = (D c), t D is diffusion coefficient. Along cell membrane, chemical boundary condition satisfies: (cu D c) n = c X t n + f c + f p, on Γ I, Γ E, where f c = k c [c] is transmembrane passive chemical flux, and f p is active flux H(s, c i, c e ) = f p = k p H(s, c i, c e ), k p = k p (s) { c i if k p (s) 0, c e if k p (s) < 0. X s 1, Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
8 The free energy identity 1 The model proposed so far is thermodynamically consistent, as it has the following energy identity (for Stokes flow): d dt (E bulk + E mem ) = I J, E bulk = ωdx, ω = RT (c ln c c), Ω i Ω e ( ( ) X E mem = k elas Γ ref s l X ) + k bend s ds, I = (ν S u + DRT ) c µ dx, µ = RT ln c, Ω i Ω e J = ([ψ]f w + [µ](f c + f p )) dm Γ, dm Γ = X s ds. Γ The entropic and membrane elastic free energy are dissipated through bulk flux (viscous and solute dissipation) and membrane flux (water and chemical flux). 1 Based on Y. Mori,, C. Liu, and R. S. Eisenberg, Physica D, 40, , (011) Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
9 The design of numerical method Fluid equations are rewritten into immersed boundary (IB) form: u t where δ is Dirac delta function. + u u = ν u p + f, u = 0, f(x, t) = F mem (X)δ(x X(s, t))ds, Γ Membrane location X is updated semi-implicitly using water permeable velocity boundary condition in the IB framework Chemical concentrations c ij are defined at all cell centers over a fixed computational grid, and updated with a Cartesian grid method. Chemical concentrations c b i and cb e are also defined along Γ i and Γ e to enforce chemical boundary condition, and increase numerical stability. Fluid structure interaction and chemicals update are decoupled on the cell membrane using fractional steps in time. Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
10 Diagram of computational grids and variables Ω 1 Γ 1 Ω 3 Γ 3 Ω Γ Ω 0 u = (u, v), p are in MAC arrangement, X i, i = 1,..., n ring are IB points to discretize cell membrane Γ, chemical on computational grid c ij are at all computational cell centers, and c b i,e are defined at crossings of the fixed Eulerian grid and Γ. Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
11 Numerical algorithm 1 Given u n, v n, p n and X n, use explicit IB method to compute u n+1, v n+1 and p n+1. Once this is found, use X n, c b,n i,e as well as the newly found u n+1, v n+1 to get the new IB locations X n+1 : f n+1 w X n+1 = X n + (u n+1 fw n+1 n) t = k w [ψ w ] = k w ( [c] n F mem (X n+1 ) Xn s 1 n) Given X n+1 i and u n+1, update chemical c n to c n+1 at all cell centers by using Cartesian grid method (need modification near membrane) c n+1 c n + (u n+1 c n+1 ) = D c n+1 t with chemical boundary condition, which is enforced with help from c b,n+1 i,e fw n+1 c D c n = k c [c] + f p Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
12 Numerical details for IB step: fluid equations u n+1 u n i,j+ 1 i,j+ 1 +[(u t i,j+ 1 C x + E x E y v i,j+ 1 C y)u i,j+ 1 ]n+ 1 = ν L(un+1 + u n i,j+ 1 i,j+ ) D 1 x p n+1 + f n i+ 1,j+ 1 x,i,j+, 1 v n+1 v n i,j+ 1 i,j+ 1 +[(E x E y u t i+ 1,j C x + v i+ 1,j C y )v i+ 1,j ] n+ 1 = ν L(vn+1 + i+ 1,j vn i+ 1,j) D y p n+1 + f n i+ 1,j+ 1 y,i+ 1,j, 0 = Dx u i+1,j+ 1 + D y v i+ 1,j+1. D ± x w α,β = ± w α±1,β w α,β h, D y ± w α,β = ± w α,β±1 w α,β, h Lw α,β = w α+1,β + w α,β+1 + w α 1,β + w α,β 1 4w α,β h, E x w α,β = w α+ 1,β + w α 1,β, C x w α,β = w α+1,β w α 1,β h Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
13 Continue IB step: Elastic force and Delta function f n x,i,j+ 1 = n ring k=1 n ring f n y,i+ 1,j = k=1 where delta function is given by F n x,kδ h (x i X n k )δ h (y j+ 1 Y n k ) s F n y,kδ h (x i+ 1 Xn k )δ h (y j Y n k ) s δ h (r) = 1 ( r ) h φ, h 1 8 (3 r r 4r ) r 1, 1 φ(r) = 8 (5 r r 4r ) 1 < r, 0 < r. Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
14 Numerical details for chemicals: auxiliary variables Extrapolations are used to ensure chemicals are updated using chemicals from proper side: c p1 = (1 θ) + θ c 3(1 θ) p θ c 6 p 5 + (1 + θ)( + θ) cb i p p1 va (xp, yp) n (x b 1, y1) b C p3 B A θ y vo (x b, y) b p4 (X n, Y n ): IB point (Marker); (x b 1, y b 1)& (x b, y b ): grid crossings of interface Γ and lines connecting cell centers between p 1 to A and p 3 to A; A, B, C: irregular cell centers Chemical boundary condition is p7 p5 p6 (Xn+1, Yn+1) (Xn, Yn) (Xn 1, Yn 1) f n+1 w c D c n = k c [c] + f p, where the gradient will be enforced using c b i (v o / v o ) 3 v o cb i + v o c B 1 v o [(1 θ)c C+θc p7 ], (v o is off cell center direction that n is decomposed into) Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
15 Numerical details of freshly cleared chemicals At freshly cleared point F (on extracellular side) is replaced by c F t c n+1 F cn+1 F c n e,p F t=n t t c n F, t ũ F D 0 xc n+1 F ṽ F D 0 yc n+1 F, t n t n+1 p1 py pf F px p where with D 0 x y w = 1 (ũ F, ṽ F ) = (D +x y w + D x y w ), ( xf x F t, y F y F t ). Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
16 Test of Cartesian grid method for updating chemicals Computational domain [0, 1] [0, 1] with initial membrane in circle: X(s) = cos s, sin s, and chemical field: c(x, y, 0) = 0.5(1.5 + sin(π(x 0.5)))(1 + sin(π(y 0.5))) The fluid velocity and cell motion are prescribed by u = u(x, y), v(x, y) = 1 4 (y 1 ), 0 ; dx dt = ( 1 4 (y 1 ) ) 1 cos t, 0 Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
17 Snapshots of chemicals and cell locations t=0.00 t=0.50 t=1.00 t=1.50 t=.00 t= Chemical Field Chemical Field Chemical Field x t =0.05 t =0.5 Test chemical module t =.5 Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
18 Convergence results of in chemical module t e I L e I L e E L e E L ɛ { x, t}/ /ɛ { x, t}/ ɛ { x, t}/ /ɛ { x, t}/ ɛ { x, t}/ /ɛ { x, t}/ ɛ { x, t}/ /ɛ { x, t}/ ɛ { x, t}/ /ɛ { x, t}/ At different times, convergence rates of chemical c on cell centers for the analytical IB locations simulation. x = 1/64, and t = Superscript I and E indicates the intracellular and extracellular cell centers, and subscript L and L indicate errors norms used. Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
19 π Test of full coupling of fluid flow and chemicals Computational domain [0, 1] [0, 1] with initial membrane configuration X(s) = cos s, sin s, and chemical field is constant 1 everywhere. Parameters are listed here k c k w D k elas k bend l r ν Case 1a kp(s) 0 3π π π s with active pump strength: twid = 0.5hwid 1 s k p (s) = k h (e (h wid ) + e (s π) (h wid ) ) + k t e (s π) where k h = k t =, h wid = t wid = 0.π, and s [0, π] (t wid ), Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
20 Velocity and chemical snapshots of full coupling Velocity Chemical Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
21 Convergence results for all variables Rate on variables t = 0.5 t = 1.5 t =.5 t = 3.75 ɛ(c E ) { x, t}/ / ɛ(c E ) { x, t}/ ɛ(c I ) { x, t}/ / ɛ(c I ) { x, t}/ ɛ(c L E ) { x, t}/ / ɛ(c L E ) { x, t}/ ɛ(c L I ) { x, t}/ / ɛ(c L I ) { x, t}/ ɛ(x) { x, t}/ / ɛ(x) { x, t}/ ɛ(u) { x, t}/ / ɛ(u) { x, t}/ At different times, convergence rates of chemicals at cell centers on intracellular side c I and extracellular side c E of Γ, chemicals at IB points on intracelular side c L I side and extracellular side cl E of Γ, and locations of IB points X are recorded. Here x = 1/64, t = 0.005, and n ring = 160 for this pair of x and t. Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
22 Full model predictions with varying membrane elasticities t=0.0 t=3. t=6.4 t= t=0.0 t=3. t=6.4 t= t=0.0 t=3. t=6.4 t= t=0.0 t=3. t=6.4 t= Comparison of cell development with different membrane elasticity while all other parameters and pump are the same. Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, 018 / 5
23 Full model predictions in permeability and chemical flux speed kc = 0.1, D = 0.1 kc = 0., D = 0. kc = 0.4, D = 0.4 speed kw = 1.11, D = 0.1 kw = 1.11, D = 0. kw = 1.11, D = 0.4 kw =., D = 0.1 kw =., D = 0. kw =., D = kw kc 1 Speed is measure by using the distance traveled by the center of mass for the cell membrane. 1 Left panel: the speed (unit µm/s) calculated to t = 40s, vs water flux coefficient k w (unit µm s/kg). Right panel: the speed calculated to t = 40s vs water passive chemical flux constant k c (unit µm/s), at different diffusion coefficient D (unit 10 3 µm /s) and k w. Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
24 Conclusions, current and future work Conclusions: We construct a thermodynamically consistent model for osmotic water flow and cell movement, and design a numerical method that is capable of handling fluid structure interactions and chemical osmosis and their coupling. Fluid structure interaction is computed using IB method with semi-implicit treatment on advancing IB points Chemical update is obtained by using a Cartesian grid method for advection diffusions in moving domains, with special cares given to enforce boundary condition and deal with freshly cleared point test cases Current and future work Include new phase of actin network inside the cell and adapt our numerical scheme to handle the full coupling of chemicals, network, and fluid flow (done implementation, testing with applications) Include charges of chemicals and electro diffusion into our model and numerical framework Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
25 Acknowledgement Yizeng Li, Sean Sun Johns Hopkins University Alex Mogilner New York University Aaron Fogelson University of Utah Funding support Yoichiro Mori: NSF DMS Lingxing Yao: NSF DMS Thank you! Any questions? Yao, Mori (CWRU, UMN) Simulating osmotic flow in cells IMA, March 15, / 5
Introduction to immersed boundary method
Introduction to immersed boundary method Ming-Chih Lai mclai@math.nctu.edu.tw Department of Applied Mathematics Center of Mathematical Modeling and Scientific Computing National Chiao Tung University 1001,
More informationElectroMechanics of Living Cells
ElectroMechanics of Living Cells Sean Sun Department of Mechanical Engineering, and Department of Biomedical Engineering, Johns Hopkins University. IMA, Minneapolis, 218 Acknowledgements Yizeng Li Flori
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
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 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 informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationInvestigating platelet motion towards vessel walls in the presence of red blood cells
Investigating platelet motion towards vessel walls in the presence of red blood cells (Complex Fluids in Biological Systems) Lindsay Crowl and Aaron Fogelson Department of Mathematics University of Utah
More informationFREE BOUNDARY PROBLEMS IN FLUID MECHANICS
FREE BOUNDARY PROBLEMS IN FLUID MECHANICS ANA MARIA SOANE AND ROUBEN ROSTAMIAN We consider a class of free boundary problems governed by the incompressible Navier-Stokes equations. Our objective is to
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationNonlinear elasticity and gels
Nonlinear elasticity and gels M. Carme Calderer School of Mathematics University of Minnesota New Mexico Analysis Seminar New Mexico State University April 4-6, 2008 1 / 23 Outline Balance laws for gels
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 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 informationNernst Equilibrium Potential. p. 1
Nernst Equilibrium Potential p. 1 Diffusion The conservation law for a compound with concentration c: rate change of c = local production + accumulation due to transport. Model: d c dv = p dv J n da dt
More information2. Conservation of Mass
2 Conservation of Mass The equation of mass conservation expresses a budget for the addition and removal of mass from a defined region of fluid Consider a fixed, non-deforming volume of fluid, V, called
More informationA Momentum Exchange-based Immersed Boundary-Lattice. Boltzmann Method for Fluid Structure Interaction
APCOM & ISCM -4 th December, 03, Singapore A Momentum Exchange-based Immersed Boundary-Lattice Boltzmann Method for Fluid Structure Interaction Jianfei Yang,,3, Zhengdao Wang,,3, and *Yuehong Qian,,3,4
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationChapter 13. Eddy Diffusivity
Chapter 13 Eddy Diffusivity Glenn introduced the mean field approximation of turbulence in two-layer quasigesotrophic turbulence. In that approximation one must solve the zonally averaged equations for
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationASTR 320: Solutions to Problem Set 2
ASTR 320: Solutions to Problem Set 2 Problem 1: Streamlines A streamline is defined as a curve that is instantaneously tangent to the velocity vector of a flow. Streamlines show the direction a massless
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 informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationA level set projection model of lipid vesicles in general flows
, A level set projection model of lipid vesicles in general flows D. Salac a, M. Miksis a a Northwestern University, Engineering Sciences and Applied Mathematics, Evanston, IL, 60208 Abstract A new numerical
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 informationA monolithic FEM solver for fluid structure
A monolithic FEM solver for fluid structure interaction p. 1/1 A monolithic FEM solver for fluid structure interaction Stefan Turek, Jaroslav Hron jaroslav.hron@mathematik.uni-dortmund.de Department of
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationSplash singularity for a free-boundary incompressible viscoelastic fluid model
Splash singularity for a free-boundary incompressible viscoelastic fluid model Pierangelo Marcati (joint work with E.Di Iorio, S.Spirito at GSSI) Workshop 2016 Modeling Computation of Shocks Interfaces
More informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
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 informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More informationDrift-Diffusion Simulation of the Ephaptic Effect in the Triad Synapse of the Retina
Drift-Diffusion Simulation of the Ephaptic Effect in the Triad Synapse of the Retina Jeremiah Jones PhD Thesis Defense, Applied Mathematics SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES April 5, 2013
More informationAn arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem
An arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem Yifan Wang University of Houston, Department of Mathematics
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationFlow simulation on moving boundary-fitted grids and application to fluid-structure interaction problems
Flow simulation on moving boundary-fitted grids and application to fluid-structure interaction problems Martin Engel and Michael Griebel Institute of Numerical Simulation, University of Bonn, Wegelerstr.
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
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 informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.306 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Monday March 12, 2018. Turn it in (by 3PM) at the Math.
More informationSupplementary Figure 1
Supplementary Figure 1 Activation of P2X2 receptor channels in symmetric Na + solutions only modestly alters the intracellular ion concentration. a,b) ATP (30 µm) activated P2X2 receptor channel currents
More informationMechanical Simulations of cell motility
Mechanical Simulations of cell motility What are the overarching questions? How is the shape and motility of the cell regulated? How do cells polarize, change shape, and initiate motility? How do they
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
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 informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationIMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS
Approved for public release; distribution is unlimited IMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS Sean L. Gibbons Captain, United
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More information3 Generation and diffusion of vorticity
Version date: March 22, 21 1 3 Generation and diffusion of vorticity 3.1 The vorticity equation We start from Navier Stokes: u t + u u = 1 ρ p + ν 2 u 1) where we have not included a term describing a
More information2.6 The Membrane Potential
2.6: The Membrane Potential 51 tracellular potassium, so that the energy stored in the electrochemical gradients can be extracted. Indeed, when this is the case experimentally, ATP is synthesized from
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationThe lattice Boltzmann method for contact line dynamics
The lattice Boltzmann method for contact line dynamics Sudhir Srivastava, J.H.M. ten Thije Boonkkamp, Federico Toschi April 13, 2011 Overview 1 Problem description 2 Huh and Scriven model 3 Lattice Boltzmann
More informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationLagrangian acceleration in confined 2d turbulent flow
Lagrangian acceleration in confined 2d turbulent flow Kai Schneider 1 1 Benjamin Kadoch, Wouter Bos & Marie Farge 3 1 CMI, Université Aix-Marseille, France 2 LMFA, Ecole Centrale, Lyon, France 3 LMD, Ecole
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
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 informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationSpatial discretization scheme for incompressible viscous flows
Spatial discretization scheme for incompressible viscous flows N. Kumar Supervisors: J.H.M. ten Thije Boonkkamp and B. Koren CASA-day 2015 1/29 Challenges in CFD Accuracy a primary concern with all CFD
More informationComputational Modelling of Mechanics and Transport in Growing Tissue
Computational Modelling of Mechanics and Transport in Growing Tissue H. Narayanan, K. Garikipati, E. M. Arruda & K. Grosh University of Michigan Eighth U.S. National Congress on Computational Mechanics
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationA Multidomain Model for Ionic Electrodiffusion and Osmosis with an Application to Cortical Spreading Depression
A Multidomain Model for Ionic Electrodiffusion and Osmosis with an Application to Cortical Spreading Depression Yoichiro Mori School of Mathematics, University of Minnesota, MN 55455, U.S.A. Abstract Ionic
More informationAnswers to Problem Set # 01, MIT (Winter-Spring 2018)
Answers to Problem Set # 01, 18.306 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Contents 1 Nonlinear solvable ODEs 2 1.1 Statement:
More informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
More informationTodd Arbogast. Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences (ICES)
CONSERVATIVE CHARACTERISTIC METHODS FOR LINEAR TRANSPORT PROBLEMS Todd Arbogast Department of Mathematics and, (ICES) The University of Texas at Austin Chieh-Sen (Jason) Huang Department of Applied Mathematics
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationA Comparison of Implicit Solvers for the Immersed Boundary Equations
A Comparison of Implicit Solvers for the Immersed Boundary Equations Elijah P. Newren a, Aaron L. Fogelson,a,b, Robert D. Guy c, Robert M. Kirby d,b a Department of Mathematics, University of Utah, Salt
More informationRegularity Theory a Fourth Order PDE with Delta Right Hand Side
Regularity Theory a Fourth Order PDE with Delta Right Hand Side Graham Hobbs Applied PDEs Seminar, 29th October 2013 Contents Problem and Weak Formulation Example - The Biharmonic Problem Regularity Theory
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationA PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM Hyun Geun LEE 1, Jeong-Whan CHOI 1 and Junseok KIM 1 1) Department of Mathematics, Korea University, Seoul
More informationA Poroelastic Immersed Boundary Method with Applications to Cell Biology
A Poroelastic Immersed Boundary Method with Applications to Cell Biology Wanda Strychalski a,, Calina A. Copos b, Owen L. Lewis b, Robert D. Guy b a Department of Mathematics, Applied Mathematics, and
More informationEstimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations.
Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations. May 6, 2009 Motivation Constitutive Equations EnKF algorithm Some results Method Navier Stokes equations
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
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 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 informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationA Second Order Energy Stable Scheme for the Cahn-Hilliard-Hele-Shaw Equations
A Second Order Energy Stable Scheme for the Cahn-Hilliard-Hele-Shaw Equations Wenbin Chen Wenqiang Feng Yuan Liu Cheng Wang Steven M. Wise August 9, 07 Abstract We present a second-order-in-time finite
More informationENGI 9420 Lecture Notes 1 - ODEs Page 1.01
ENGI 940 Lecture Notes - ODEs Page.0. Ordinary Differential Equations An equation involving a function of one independent variable and the derivative(s) of that function is an ordinary differential equation
More informationAn introduction to the Immersed boundary method and its finite element approximation
FE An introduction to te and its finite element approximation Collaborators: Daniele offi, Luca Heltai and Nicola Cavallini, Pavia Pavia, 27 maggio 2010 FE History Introduced by Peskin Flow patterns around
More information2 Law of conservation of energy
1 Newtonian viscous Fluid 1 Newtonian fluid For a Newtonian we already have shown that σ ij = pδ ij + λd k,k δ ij + 2µD ij where λ and µ are called viscosity coefficient. For a fluid under rigid body motion
More informationGauge finite element method for incompressible flows
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 34: 701 710 Gauge finite element method for incompressible flows Weinan E a, *,1 and Jian-Guo Liu b,2 a Courant Institute
More information.u= 0 ρ( u t +(u. )u)= ρ g p+.[µ( u+ t u)]
THETIS is a numerical simulation tool developed by University of Bordeaux. It is a versatile code to solve different problems: fluid flows, heat transfers, scalar transports or porous mediums. The potential
More informationA dynamic model of polyelectrolyte gels. A dissertation submitted to the faculty of the graduate school of the university of minnesota by.
A dynamic model of polyelectrolyte gels A dissertation submitted to the faculty of the graduate school of the university of minnesota by Haoran Chen In partial fulfillment of the requirements for the degree
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF
More informationSupplement to Molecular Gas TitleBirkhäuser, Boston, 007 Dynamic Version Authors Sone, Yoshio Citation Yoshio Sone. 008 Issue Date 008-09-0 URL http://hdl.handle.net/433/66098 Right Type Book Textversion
More informationDaniel J. Jacob, Models of Atmospheric Transport and Chemistry, 2007.
1 0. CHEMICAL TRACER MODELS: AN INTRODUCTION Concentrations of chemicals in the atmosphere are affected by four general types of processes: transport, chemistry, emissions, and deposition. 3-D numerical
More informationINCOMPRESSIBLE FLUIDS IN THIN DOMAINS WITH NAVIER FRICTION BOUNDARY CONDITIONS (II) Luan Thach Hoang. IMA Preprint Series #2406.
INCOMPRESSIBLE FLUIDS IN THIN DOMAINS WITH NAVIER FRICTION BOUNDARY CONDITIONS II By Luan Thach Hoang IMA Preprint Series #2406 August 2012 INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY OF
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 informationFree energy concept Free energy approach LBM implementation Parameters
BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH
More informationIntroduction to Physiology II: Control of Cell Volume and Membrane Potential
Introduction to Physiology II: Control of Cell Volume and Membrane Potential J. P. Keener Mathematics Department Math Physiology p.1/23 Basic Problem The cell is full of stuff: Proteins, ions, fats, etc.
More informationQuantifying Intermittent Transport in Cell Cytoplasm
Quantifying Intermittent Transport in Cell Cytoplasm Ecole Normale Supérieure, Mathematics and Biology Department. Paris, France. May 19 th 2009 Cellular Transport Introduction Cellular Transport Intermittent
More informationSmoluchowski Navier-Stokes Systems
Smoluchowski Navier-Stokes Systems Peter Constantin Mathematics, U. of Chicago CSCAMM, April 18, 2007 Outline: 1. Navier-Stokes 2. Onsager and Smoluchowski 3. Coupled System Fluid: Navier Stokes Equation
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 informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
More informationComputational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. # 02 Conservation of Mass and Momentum: Continuity and
More informationSurface phase separation and flow in a simple model of drops and vesicles
Surface phase separation and flow in a simple model of drops and vesicles Tutorial Lecture 4 John Lowengrub Department of Mathematics University of California at Irvine Joint with J.-J. Xu (UCI), S. Li
More information