Introduction to immersed boundary method
|
|
- Phillip Walker
- 5 years ago
- Views:
Transcription
1 Introduction to immersed boundary method Ming-Chih Lai Department of Applied Mathematics Center of Mathematical Modeling and Scientific Computing National Chiao Tung University 1001, Ta Hsueh Road, Hsinchu Taiwan Workshop on Fluid Motion Driven by Immersed Structure Fields Institute, August 9-13, 2010
2 Review of the IB method A general computational framework for simulations of problems involving the interaction of fluids and immersed elastic structures Original example: blood interacts with valve leaflet (Charles S. Peskin, 1972, flow patterns around heart valves) Applications Computer-assisted design of prothetic valve (Peskin & McQueen) Platelet aggregation during blood clotting (Fogelson) Flow of particle suspensions (Fogelson & Peskin; Sulsky & Brackbill) Wave propagation in the cochlea (Beyer)
3 Swimming organism (Fauci, Dillon, Cortez) Arteriolar flow (Arthurs, et. al.) Cell and tissue deformation under shear flow (Bottino, Stockie & Green; Eggleton & Popel) Flow around a circular cylinder (Lai & Peskin; Su, Lai & Lin), Interfacial flow with insoluble surfactant (Lai, Tseng & Huang) Valveless pumping (Jung & Peskin) Flapping filament in a flowing soap film (Zhu & Peskin), 2D dry foam (Kim, Lai & Peskin) Falling papers, sails, parachutes (Kim & Peskin), insect flight, (Jane Wang, Miller)... Many others (Lim,..)
4 IB features Mathematical formulation: Treat the elastic material as a part of fluid The material acts force into the fluid The material moves along with the fluid Numerical method: Finite difference discretization Eulerian grid points for the fluid variables Lagrangian markers for the immersed boundary The fluid-boundary are linked by a smooth version of Dirac delta function
5 Mathematical formulation Consider a massless elastic membrane Γ immersed in viscous incompressible fluid domain Ω. Γ : X(s, t), 0 s L b L b : unstressed length
6 Equations of motion ( ) u ρ + (u )u + p = µ u + f t X(s, t) t F (s, t) = (T τ ), s T u = 0 f(x, t) = = u(x(s, t), t) = Γ Ω = σ( X s F (s, t)δ(x X(s, t))ds u(x, t)δ(x X(s, t))dx ; s, t), τ = X/ s X/ s FLUID u(x, t) : velocity p(x, t) : pressure ρ : density µ : viscosity IMMERSED BOUNDARY X(x, t) : boundary configuration F (x, t) : boundary force T : tension τ : unit tangent
7 Why F (s, t) = s (T τ )? total force = force of fluid on boundary seg. + (T τ ) s2 s 1 boundary is massless total force = 0 force of boundary seg. on fluid = (T τ ) s2 s 1 = F = T s τ + T τ s = T s τ + T X s κn s2 s 1 (T τ )ds s where σ 0 : the stiffness constant τ / s κ = X/ s τ / s n = τ / s T = σ 0 ( X s boundary curvature unit normal ) 1
8 The force density is singular f(x, t) = Γ F (s, t)δ(x X(s, t))ds. f behaves like a one-dimensional delta function. f, w = f(s, t)w(x, t)dx Ω = F (s, t)δ(x X(s, t))ds w(x, t)dx Ω Γ = F (s, t) w(x, t)δ(x X(s, t))dxds Γ Ω = F (s, t)w(x(s, t), t)ds Γ If w(x, t) = u(x, t), the the total work done by the boundary is equal to the total work done on the fluid. Thus, the solution is NOT smooth. In fact the pressure and velocity derivatives are discontinuous across the boundary.
9 Theorem (1) The pressure and the velocity normal derivatives across the boundary satisfy µ [ u n [p] = ] F n X/ s, = F τ X/ s τ. Theorem (2) The normal derivatives of the pressure across the boundary satisfies ( ) [ ] F τ p s X/ s =. n X/ s Physical meaning of the jump conditions. Ref: Peskin & Printz 1993, LeVeque & Li 1997 (2D), Lai & Li 2001 (3D)
10 Test Example u t + (u )u + p = u + f + g u = 0 f = X(s, t) t Γ F (s, t)δ(x X(s, t))ds = u(x(s, t), t) ( ) ( ) X X(s), Y (s) = 0.5 cos s, 0.5 sin s, s = 0.5 n = ( cos s, sin s ), τ = ( sin s, cos s )
11 Test Example (cont.) u(x, y, t) = v(x, y, t) = p(x, y, t) = g 1 (x, y, t) = g 2 (x, y, t) = { e ( t 2y y ) r r r < 0.5, { e ( ) t 2x + x r r 0.5 [ ] u = 2e t sin s n [ ] v = 2e t cos s n 0 r < 0.5, { 0 r 0.5, [p] = 1 1 r < 0.5 { e ( t y r 2y y ) ( r + e 2t 4x 3 r 4x ) x r r r < 0.5 { e ( t x r 2x ) ( x r + e 2t 4y 3 r 4y y ) r r r < 0.5 Fn = 0.5 F 1 = 0.5 cos s e t sin s [g n] = 0 Fτ = e t F 2 = 0.5 sin s + e t cos s
12 Numerical Method FLUID-Eulerian N N lattice points x = (ih, jh) u n u(x, n t), f n f(x, n t), p n p(x, n t), IMMERSED BOUNDARY-Lagrangian M moving points X k U n k u(x k, n t) F n k F (X k, n t)
13 Finite difference operators on lattice points and immersed boundary points The finite difference approximations to the derivative are D + x φ i,j = φ i+1,j φ i,j h, Dx φ i,j = φ i,j φ i 1,j, h Dxφ 0 i,j = φ i+1,j φ i 1,j. 2h The upwind difference approximation to the advection term is { ψ i,j D x ± ψ i,j Dx φ i,j if ψ i,j > 0, φ i,j = ψ i,j D x + φ i,j if ψ i,j < 0. The centered difference approximation to the derivative on the Lagrangian boundary can be defined by D s χ k = χ k+1/2 χ k 1/2, s where χ k+1/2 is some boundary quantity defined on s = (k ) s.
14 Numerical algorithm How to march (X n, u n ) to (X n+1, u n+1 )? Step1 Compute the boundary force T n k+1/2 = σ( D sx n k ), τ n k+1/2 = D sx n k D s X n k, F n k = D s (T τ ) n k+1/2, where Tk+1/2 n and τ n k+1/2 are both defined on s = (k ) s, and F n k is defined on s = k s. Step2 Apply the boundary force to the fluid f n (x) = k F n kδ h (x X n k) s. Step3 Solve the Navier-Stokes equations with the force to update the velocity ( ) u n+1 u n 2 2 ρ + u n i D ± i t un = D 0 p n+1 + µ D + i D i un+1 + f n, i=1 D 0 u n+1 = 0, where T n k+1/2 and τ n k+1/2 are both defined on s = (k ) s, and F n k is defined on s = k s. i=1
15 Numerical algorithm, cont. Step4 Interpolate the new velocity on the lattice into the boundary points and move the boundary points to new positions. U n+1 k = x u n+1 δ h (x X n k)h 2 X n+1 k = X n k + tu n+1 k
16 Discrete delta function δ h δ h (x) = δ h (x)δ h (y) 1. δ h is positive and continuous function. 2. δ h (x) = 0, for x 2h. 3. j δ h(x j α)h = 1 for all α. ( δ h(x j α)h = δ h (x j α)h = 1 2 ) j even j odd 4. j (x j α)δ h (x j α)h = 0 for all α. 5. j (δ h(x j α)h) 2 = C for all α. ( j δ h(x j α)δ h (x j β)h 2 C) Uniquely determined: C = h (3 2 x /h + ) x /h 4( x /h) 2 x h, δ h (x) = 1 8h (5 2 x /h ) x /h 4( x /h) 2 h x 2h, 0 otherwise.
17 Numerical issues of IB method Simple and easy to implement Embedding complex structure into Cartesian domain, no complicated grid generation Fast elliptic solver (FFT) on Cartesian grid can be applied Numerical smearing near the immersed boundary First-order accurate, accuracy of 1D IB model (Beyer & LeVeque 1992, Lai 1998), formally second-order scheme (Lai & Peskin 2000, Griffith & Peskin 2005) Adaptive IB method, (Roma, Peskin & Berger 1999, Griffith et. al. 2007) Immersed Interface Method (IIM, LeVeque & Li 1994), 3D jump conditions (Lai & Li 2001) High-order discrete delta function in 2D, 3D Numerical stability tests, different semi-implicit methods (Tu & Peskin 1992, Mayo & Peskin 1993, Newren, Fogelson, Guy & Kirby 2007, 2008) Stability analysis (Stockie & Wetton 1999, Hou & Shi 2008) Convergence analysis (Y. Mori 2008)
18 Review articles C. S. Peskin, The immersed boundary method, Acta Numerica, pp 1-39, (2002) R. Mittal & G. Iaccarino, Immersed boundary methods, Annu. Rev. Fluid Mech. 37: , (2005), Flow around solid bodies Accuracy, stability and convergence R. P. Beyer & R. J. LeVeque, Analysis of a one-dimensional model for the immersed boundary method, SIAM J. Numer. Anal., 29: , (1992) E. P. Newren, A. L. Fogelson, R. D. Guy & R. M. Kirby, Unconditionally stable discretizations of the immersed boundary equations, JCP, 222: , (2007) Y. Mori, Convergence proof of the velocity field for a Stokes flow immersed boundary method, CPAM, vol LXI , (2008)
19 Application: Flow around a solid body Figure: The fluid feels the force along the surface body to stop it. u t + (u )u + p = 1 u = 0 f(x, t) = 0 = u b (X(s), t) = Re u + f Lb 0 Ω F (X(s), t)δ(x X(s))ds u(x, t)δ(x X(s))dx u(x, t) u as x
20 u n+1 u n + (u n )u n = p n t Re un+1 + f n+1, u n+1 = 0, M f n+1 (x) = F n+1 (X j )δ h (x X j ) s, for all x, j=1 U n+1 b (X k ) = x u n+1 (x)δ h (x X k )h 2, k = 1,..., M. One can solve the above linear system directly or approximately (A projection approach, Taira & Colonius, 2007).
21 Feedback forcing approach Goldstein, et. al. 1993; Saiki & Biringen 1996; Lai & Peskin 2000 u t + (u )u + p = 1 Re u + F (s, t)δ(x X(s, t))ds F (s, t) = K t 0 u = 0 (U e (s, t ) U(s, t ))dt + R(U e (s, t) U(s, t)) or F (s, t) = K(X e (s) X(s, t)), K 1 X(s, t) = U(s, t) = u(x, t)δ(x X(s, t))dx t Ω u(x, t) u as x Γ Treat the body surface as a nearly rigid boundary (stiffness is large) embedded in the fluid Allow the boundary to move a little but the force will bring it back to the desired location Be able to handle complicated geometry with known boundary velocity
22 Interpolating forcing approach Su, Lai & Lin 2007 Idea: To compute the boundary force F (X k ) accurately so the prescribed boundary velocity can be achieved in the intermediate step of projection method. 1. Distribute the force to the grid by the discrete delta function M f (x) = F (X j )δ h (x X j ) s 2. j=1 u u = f, thus, u (X k ) = u b (X k ). t u b (X k ) u (X k ) t = f (x)δ h (x X k )h 2 x = M F (X j )δ h (x X j ) s δ h (x X k )h 2 x j=1 ( ) M = δ h (x X j )δ h (x X k )h 2 s F (X j ) j=1 x
23 Re Su et al. Lai & Peskin Silva et al. Ye et al. Willimson(exp) Table: The comparison of Strouhal number for different Reynolds numbers.
An Immersed Boundary Method with Formal Second-Order Accuracy and Reduced Numerical Viscosity 1
Journal of Computational Physics 160, 705 719 (2000) doi:10.1006/jcph.2000.6483, available online at http://www.idealibrary.com on An Immersed Boundary Method with Formal Second-Order Accuracy and Reduced
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 informationAn Immersed Boundary Method for Restricted Diffusion with Permeable Interfaces
An Immersed Boundary Method for Restricted Diffusion with Permeable Interfaces Huaxiong Huang Kazuyasu Sugiyama Shu Takagi April 6, 009 Keywords: Restricted diffusion; Permeable interface; Immersed boundary
More informationImmersed boundary method:
Immersed boundary method: Recent developments in analysis, algorithms and applications John Stockie Department of Mathematics Simon Fraser University Burnaby, British Columbia, Canada http://www.math.sfu.ca/
More informationAn Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains
An Immersed Interace Method or the Incompressible Navier-Stokes Equations in Irregular Domains Duc-Vinh Le, Boo Cheong Khoo, Jaime Peraire Singapore-MIT Alliance Department o Mechanical Engineering, National
More informationAn immersed boundary method for interfacial flows with insoluble surfactant
An immersed boundary method for interfacial flows with insoluble surfactant Ming-Chih Lai Yu-Hau Tseng Huaxiong Huang Abstract In this paper, an immersed boundary method is proposed for the simulation
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 informationSimulation of a Flapping Flexible Filament in a Flowing Soap Film by the Immersed Boundary Method
Journal of Computational Physics 179, 452 468 (2002) doi:10.1006/jcph.2002.7066 Simulation of a Flapping Flexible Filament in a Flowing Soap Film by the Immersed Boundary Method Luoding Zhu 1 and Charles
More informationModelling and simulation of porous immersed boundaries
Modelling and simulation of porous immersed boundaries John M. Stockie Department of Mathematics, Simon Fraser University, 8888 University Drive, Burnaby, BC, V5A 1S6, Canada Abstract The immersed boundary
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 informationTETHER FORCE CONSTRAINTS IN STOKES FLOW BY THE IMMERSED BOUNDARY METHOD ON A PERIODIC DOMAIN
SIAM J. SCI. COMPUT. Vol. 31, No. 5, pp. 3404 3416 c 2009 Society for Industrial and Applied Mathematics TETHER FORCE CONSTRAINTS IN STOKES FLOW BY THE IMMERSED BOUNDARY METHOD ON A PERIODIC DOMAIN JOSEPH
More informationTWO DIMENSIONAL IMMERSED BOUNDARY SIMULATIONS OF SWIMMING JELLYFISH
TWO DIMENSIONAL IMMERSED BOUNDARY SIMULATIONS OF SWIMMING JELLYFISH by Haowen Fang B.Eng., Nanjing University of Science and Technology, 2006 a Thesis submitted in partial fulfillment of the requirements
More informationNumerical method for osmotic water flow and solute diffusion in moving cells
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
More informationNew Finite Difference Methods Based on IIM for Inextensible Interfaces in Incompressible Flows
East Asian Journal of Applied Mathematics Vol. x, No. x, pp. -7 doi:.48/eajam xxx x New Finite Difference Methods Based on IIM for Inextensible Interfaces in Incompressible Flows Zhilin Li, and Ming-Chih
More informationA new 3D immersed boundary method with application
A new 3D immersed boundary method with application Luoding Zhu Department of Mathematical Sciences Indiana Univ - Purdue Univ Indianapolis Runan Hua, Xiyun Lu Department of Modern Mechanics University
More information2 Z. Li and C. Wang where u = (u; v) is the velocity, p is the pressure, ν is the viscosity. We assume that the boundary of the is piecewise
A fast finite difference method for solving Navier-Stokes Equations on irregular domains Zhilin Li Λ Cheng Wang y April 17, 2002 keywords: Navier-Stokes equations, irregular domains, vorticity stream-function
More informationEfficient solutions to robust, semi-implicit discretizations of the Immersed Boundary Method
Efficient solutions to robust, semi-implicit discretizations of the Immersed Boundary Method Hector D. Ceniceros Department of Mathematics, University of California Santa Barbara, CA 93106 Jordan E. Fisher
More informationAN AUGMENTED IIM-LEVEL SET METHOD FOR STOKES EQUATIONS WITH DISCONTINUOUS VISCOSITY
Sixth Mississippi State Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conference 15 (2007), pp. 193 210. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More information2-D Parachute Simulation by the Immersed Boundary Method
-D Parachute Simulation by the Immersed Boundary Method Yongsam Kim and Charles S. Peskin Abstract Parachute aerodynamics involves an interaction between the flexible, elastic, porous parachute canopy
More informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for this
More informationEducation: Ph.D. Mathematics, Courant Institute of Mathematical Sciences, New York University, 1998/09.
Ming-Chih Lai Department of Applied Mathematics, National Chiao Tung University, Hsinchu 30050, Taiwan Tel : +886-3-5131361 E-mail: mclai@math.nctu.edu.tw http://www.math.nctu.edu.tw/ mclai/ Education:
More informationc 2003 International Press
COMMUNICATIONS IN MATH SCI Vol. 1, No. 1, pp. 181 197 c 2003 International Press A FAST FINITE DIFFERENCE METHOD FOR SOLVING NAVIER-STOKES EQUATIONS ON IRREGULAR DOMAINS ZHILIN LI AND CHENG WANG Abstract.
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 informationSpontaneous Symmetry Breaking of Hinged Flapping Filament Generates Lift
Spontaneous Symmetry Breaking of Hinged Flapping Filament Generates Lift Shervin Bagheri (KTH, Stockholm) Andrea Mazzino & Alessandro Bottaro (Genova University, Italy) Fluid & Elasticity 2012 La Jolla,
More informationA Computational Model of the Mechanics of Growth of the Villous Trophoblast Bilayer
Bulletin of Mathematical Biology (24) 66, 99 232 doi:.6/j.bulm.23.6. A Computational Model of the Mechanics of Growth of the Villous Trophoblast Bilayer KATARZYNA A. REJNIAK Mathematical Biosciences Institute,
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 informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
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 informationAdvances in the mathematical theory of the finite element immersed boundary method
Advances in the mathematical theory of the finite element immersed boundary method Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 12, 2014
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 informationThe immersed boundary method
Acta Numerica (2002), pp. 479 517 c Cambridge University Press, 2002 DOI: 10.1017/S0962492902000077 Printed in the United Kingdom The immersed boundary method Charles S. Peskin Courant Institute of Mathematical
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 informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationA STRONG COUPLING SCHEME FOR FLUID-STRUCTURE INTERACTION PROBLEMS IN VISCOUS INCOMPRESSIBLE FLOWS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2005 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2005 A STRONG COUPLING SCHEME
More informationVesicle dynamics in Newtonian and non-newtonian fluids
Vesicle dynamics in Newtonian and non-newtonian fluids Ming-Chih Lai mclai@math.nctu.edu.tw Department of Applied Mathematics National Chiao Tung University Taiwan Modelling and Simulation of Interface
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 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 informationABSTRACT NUMERICAL SIMULATION OF LIQUID - VAPOR INTERFACES IN THE SHARP INTERFACE LIMIT
ABSTRACT Title of Thesis: NUMERICAL SIMULATION OF LIQUID - VAPOR INTERFACES IN THE SHARP INTERFACE LIMIT Vamsee K. Yerramilli, Master of Science, 2005 Thesis directed by: Dr. Elias Balaras Department of
More informationInteraction of Incompressible Fluid and Moving Bodies
WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics
More informationImmersed Boundary Method Based on the Implementation of Conservation Equations along the Boundary using Control-Volume Finite-Element Scheme
doi:.58/jatm.v9i3.67 Immersed Boundary Method Based on the Implementation of Conservation Equations along the Boundary using Control-Volume Finite-Element Scheme Seyedeh Nasrin Hosseini, Seyed Mohammad
More informationABSTRACT. coefficients. The main motivation is to get not only a second order accurate solution but
ABSTRACT CHEN, GUANYU. Accurate Gradient Computation for Elliptic Interface Problems with Discontinuous and Variable Coefficients. (Under the direction of Dr. Zhilin Li.) A new numerical method is proposed
More informationANALYSIS AND NUMERICAL METHODS FOR SOME CRACK PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 2, Number 2-3, Pages 155 166 c 2011 Institute for Scientific Computing and Information ANALYSIS AND NUMERICAL METHODS FOR SOME
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationA Numerical Method for Two Phase Flows
A Numerical Method for Two Phase Flows with Insoluble and Soluble Surfactants by Shilpa Madangopal Khatri A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of
More informationOn the Solution of the Elliptic Interface Problems by Difference Potentials Method
On the Solution of the Elliptic Interface Problems by Difference Potentials Method Yekaterina Epshteyn and Michael Medvinsky Abstract Designing numerical methods with high-order accuracy for problems in
More informationThere are a host of physical problems that involve interconnected
The Voronoi Implicit Interface Method for computing multiphase physics Robert I. Saye a,b and James A. Sethian a,b,1 a Department of Mathematics; and b Lawrence Berkeley National Laboratory, University
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 informationFast Direct Solver for Poisson Equation in a 2D Elliptical Domain
Fast Direct Solver for Poisson Equation in a 2D Elliptical Domain Ming-Chih Lai Department of Applied Mathematics National Chiao Tung University 1001, Ta Hsueh Road, Hsinchu 30050 Taiwan Received 14 October
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationA Fluctuating Immersed Boundary Method for Brownian Suspensions of Rigid Particles
A Fluctuating Immersed Boundary Method for Brownian Suspensions of Rigid Particles Aleksandar Donev Courant Institute, New York University APS DFD Meeting San Francisco, CA Nov 23rd 2014 A. Donev (CIMS)
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 informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
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 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 informationVORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS
THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1395-1399 1395 VORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS by Li-Zhong HUANG a* and De-Ming NIE b a State Key Laboratory
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
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 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 informationMany of the smallest flying insects clap their wings together at the end of each upstroke
DRAFT Miller, L. A. and Peskin, C. S. Title: Flexible clap and fling in tiny insect flight. Abstract Many of the smallest flying insects clap their wings together at the end of each upstroke and fling
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationWhen vortices stick: an aerodynamic transition in tiny insect flight
The Journal of Experimental Biology 7, 7-88 Published by The Company of Biologists 4 doi:.4/jeb.8 7 When vortices stick: an aerodynamic transition in tiny insect flight Laura A. Miller* and Charles S.
More informationThe Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations
The Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations George Biros Lexing Ying and Denis Zorin Courant Institute of Mathematical Sciences, New York University, NY 10012,
More informationModelling of interfaces and free boundaries
University of Regensburg Regensburg, March 2009 Outline 1 Introduction 2 Obstacle problems 3 Stefan problem 4 Shape optimization Introduction What is a free boundary problem? Solve a partial differential
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
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 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 informationA Cartesian Grid Method for Solving the Two-Dimensional Streamfunction-Vorticity Equations in Irregular Regions
Journal of Computational Physics 176, 231 275 (2002) doi:10.1006/jcph.2001.6970, available online at http://www.idealibrary.com on A Cartesian Grid Method for Solving the Two-Dimensional Streamfunction-Vorticity
More informationAn Eulerian-Lagrangian approach for fluid-structure coupled systems using X-FEM Antoine Legay / CNAM-Paris
An Eulerian-Lagrangian approach for fluid-structure coupled systems using X-FEM Antoine Legay / CNAM-Paris Collaborations: A. Tralli, P. Gaudenzi, Università di Roma La Sapienza, Italia 8th USNCCM, 25-27th
More information14. FLUID DYNAMICS OF THE HEART AND ITS VALVES
14. FLUID DYNAMICS OF THE HEART AND ITS VALVES Charles S. Peskin and David M. McQueen 14.1 INTRODUCTION Blood is a viscous incompressible fluid which is propelled through the arteries, capillaries, and
More informationA fast, robust, and non-stiff Immersed Boundary Method
A fast, robust, and non-stiff Immersed Boundary Method Hector D. Ceniceros Department of Mathematics, University of California Santa Barbara, CA 93106 Jordan E. Fisher Department of Mathematics, University
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
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 informationAugmenting the Immersed Boundary Method with Radial Basis Functions (RBFs) for the Modeling of Platelets in Hemodynamic Flows
Boise State University ScholarWorks Mathematics Faculty Publications and Presentations Department of Mathematics 12-10-2015 Augmenting the Immersed Boundary Method with Radial Basis Functions (RBFs) for
More informationFluid Flow in Open Helical Channels
Hayden Tronnolone Supervisor: Yvonne Stokes School of Mathematical Sciences September 24, 2010 Much research has been conducted into flows in closed, helical pipes. Such flows have application to modelling
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationMarch Algebra 2 Question 1. March Algebra 2 Question 1
March Algebra 2 Question 1 If the statement is always true for the domain, assign that part a 3. If it is sometimes true, assign it a 2. If it is never true, assign it a 1. Your answer for this question
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 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 informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationComputing Fluid-Structure Interaction by the Partitioned Approach with Direct Forcing
Commun. Comput. Phys. doi: 10.4208/cicp.080815.090516a Vol. xx, No. x, pp. 1-29 xxx 201x Computing Fluid-Structure Interaction by the Partitioned Approach with Direct Forcing Asim Timalsina 1, Gene Hou
More informationExplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code
Eplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code V. Rutka A. Wiegmann November 25, 2005 Abstract The Eplicit Jump Immersed Interface method is a powerful tool to solve elliptic
More informationNew York Journal of Mathematics. Fluid Flow in Collapsible Elastic Tubes: A Three-Dimensional Numerical Model
New York Journal of Mathematics New York J. Math. 7 (2001) 281 302. Fluid Flow in Collapsible Elastic Tubes: A Three-Dimensional Numerical Model M. E. Rosar and Charles S. Peskin Abstract. A three-dimensional
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationA COMPUTATIONAL FLUID DYNAMICS STUDY OF CLAP AND FLING IN THE SMALLEST INSECTS. Laura A. Miller* and Charles S. Peskin**
A COMPUTATIONAL FLUID DYNAMICS STUDY OF CLAP AND FLING IN THE SMALLEST INSECTS Laura A. Miller* and Charles S. Peskin** *Department of Mathematics, University of Utah, 155 South 1400 East, Salt Lake City,
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 informationA FLUID-STRUCTURE INTERACTION MODEL BASED ON PERIDYNAMICS AND NAVIER-STOKES EQUATIONS FOR HYDRAULIC FRACTURE PROBLEMS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK A FLUID-STRUCTURE INTERACTION MODEL BASED ON PERIDYNAMICS
More informationFrom Immersed Boundary Method to Immersed Continuum Method
From Immersed Boundary Method to Immersed Continuum Method X. Sheldon Wang Department of Mathematical Sciences New Jersey Institute of Technology Newark, NJ 07102 CAMS Report 0506-7, Spring 2006 Center
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 informationSimulation of free surface fluids in incompressible dynamique
Simulation of free surface fluids in incompressible dynamique Dena Kazerani INRIA Paris Supervised by Pascal Frey at Laboratoire Jacques-Louis Lions-UPMC Workshop on Numerical Modeling of Liquid-Vapor
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 informationHydro-elastic Wagner impact using variational inequalities
Hydro-elastic Wagner impact using variational inequalities Thomas GAZZOLA, Alexander KOROBKIN, Šime MALENICA Introduction A simple model of water impact has been introduced by Wagner [6]. This model is
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 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 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 informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationA MOTION OF FREELY OSCILLATING DROPLET OF A YIELD STRESS FLUID: ANALYSIS AND NUMERICAL STUDIES
A MOTION OF FREELY OSCILLATING DROPLET OF A YIELD STRESS FLUID: ANALYSIS AND NUMERICAL STUDIES A Dissertation Presented to the Faculty of the Department of Mathematics University of Houston In Partial
More informationAccepted Manuscript. A study of self-propelled elastic cylindrical micro-swimmers using modeling and computation
Accepted Manuscript A study of self-propelled elastic cylindrical micro-swimmers using modeling and computation Lingling Shi, Sunčica Čanić, Annalisa Quaini, Tsorng-Whay Pan PII: S1-9991(16)143-1 DOI:
More informationIFE for Stokes interface problem
IFE for Stokes interface problem Nabil Chaabane Slimane Adjerid, Tao Lin Virginia Tech SIAM chapter February 4, 24 Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 / 2 Problem statement
More informationarxiv: v2 [math.na] 23 Nov 2015
An Immersed Boundary Method for Rigid Bodies Bakytzhan Kallemov, 1 Amneet Pal Singh Bhalla, 2 Boyce E. Griffith, 3 and Aleksandar Donev 1, 1 Courant Institute of Mathematical Sciences, New York University,
More information