Hostname: page-component-76fb5796d-wq484 Total loading time: 0 Render date: 2024-04-25T08:09:34.663Z Has data issue: false hasContentIssue false

Viscoelastic Immersed Boundary Methods for Zero Reynolds Number Flow

Published online by Cambridge University Press:  20 August 2015

Wanda Strychalski*
Affiliation:
Department of Mathematics, University of California, Davis CA 95616, USA
Robert D. Guy*
Affiliation:
Department of Mathematics, University of California, Davis CA 95616, USA
*
Corresponding author.Email:wanda@math.ucdavis.edu
Get access

Abstract

The immersed boundary method has been extensively used to simulate the motion of elastic structures immersed in a viscous fluid. For some applications, such as modeling biological materials, capturing internal boundary viscosity is important. We present numerical methods for simulating Kelvin-Voigt and standard linear viscoelastic structures immersed in zero Reynolds number flow. We find that the explicit time immersed boundary update is unconditionally unstable above a critical boundary to fluid viscosity ratio for a Kelvin-Voigt material. We also show there is a severe time step restriction when simulating a standard linear boundary with a small relaxation time scale using the same explicit update. A stable implicit method is presented to overcome these computation challenges.

Type
Research Article
Copyright
Copyright © Global Science Press Limited 2012

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Alberts, B., Johnson, A., Lewis, J., Raff, M., Roberts, K., and Walter, P.. Molecular Biology of the Cell. Garland Science, New York, 4th edition, 2002.Google Scholar
[2]Batchelor, G.K.. An Introduction to Fluid Dynamics. Cambridge University Press, Cambridge, 1967.Google Scholar
[3]Beyer, Richard P. Jr., A computational model of the cochlea using the immersed boundary method. J. Comput. Phys., 98:145–162, 1992.Google Scholar
[4]Bottino, D. C.. Modeling viscoelastic networks and cell deformation in the context of the immersed boundary method. J. Comput. Phys., 147(1):86–113, 1998.CrossRefGoogle Scholar
[5]Bray, D.. Cell Movements: From Molecules to Motility. Garland Science, New York, 2nd edition, 2001.Google Scholar
[6]Ceniceros, H. D., Fisher, J. E., and Roma, A. M.. Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method. J. Comput. Phys., 228:7137–7158, 2009.Google Scholar
[7]Chrispell, J.C., Cortez, R., Khismatullin, D.B., and Fauci, L.J.. Shape oscillations of a droplet in an oldroyd-b fluid. Physica D: Nonlinear Phenomena, In Press, 2011.Google Scholar
[8]Cortez, R.. The method of regularized stokeslets. SIAM J. Sci. Comput., 23(4):1204–1225, 2001.Google Scholar
[9]Dong, C., Skalak, R., Sung, K. L., Schmid-Schönbein, G. W., and Chien, S.. Passive deformation analysis of human leukocytes. J. Biomech. Eng., 110(1):27–36, 1988.Google Scholar
[10]Fogelson, A. L.. A mathematical model and numerical method for studying platelet adhesion and aggregation during blood clotting. J. Comput. Phys., 56(1):111–134, 1984.Google Scholar
[11]Fung., Y. CBiomechanics: mechanical properties of living tissues. Springer-Verlag, New York, 2nd edition, 1993.CrossRefGoogle Scholar
[12]Goldstein, D., Handler, R., and Sirovich, L.. Modeling a no-slip flow boundary with an external force field. J. Comput. Phys., 105(2):354–366, 1993.Google Scholar
[13]Griffith, B.E. and Peskin, C.S.. On the order of accuracy of the immersed boundary method: Higher order convergence rates for sufficiently smooth problems. J. Comput. Phys., 208(1):75–105, 2005.Google Scholar
[14]Guy, R.D. and Philip, B.. A multigrid method for a model of the implicit immersed boundary equations. Commun. Comput. Phys. in press, 2012.Google Scholar
[15]Hou, T. Y. and Shi, Z.. An efficient semi-implicit immersed boundary method for the navier-stokes equations. J. Comput. Phys., 227:8968–8991, 2008.CrossRefGoogle Scholar
[16]Huang, W. X. and Sung, H. J.. An immersed boundary method for fluid-flexible structure interaction. Comput. Method. Appl. M., 198(33-36):2650–2661, 2009.Google Scholar
[17]LeVeque, R. J. and Li, Z.. Immersed interface methods for stokes flow with elastic boundaries or surface tension. SIAM J. Sci. Comput., 18(3):709–735, 1997.Google Scholar
[18]Mayer, M., Depken, M., Bois, J. S., Jülicher, F., and Grill, S. W.. Anisotropies in cortical tension reveal the physical basis of polarizing cortical flows. Nature, 467(7315):617–21, 2010.Google Scholar
[19]Mayo, A. A. and Peskin, C. S.. An implicit numerical method for fluid dynamics problems with immersed elastic boundaries. Fluid Dynamics in Biology: Proc. AMS–IMS–SIAM Joint Summer Research Conf. on Biofluiddynamics, 141 of Contemporary Mathematics, AMS:261–277, 1993.Google Scholar
[20]Miller, L. A. and Peskin, C. S.. When vortices stick: An aerodynamic transition in tiny insect flight. J. Exp. Biol., 207(17):3073–3088, 2004.Google Scholar
[21]Mori, Y. and Peskin., C. S.Implicit second-order immersed boundary methods with boundary mass. Comp. Method. Appl. M., 197(25-28):2049–2067, 2008. Immersed Boundary Method and Its Extensions.Google Scholar
[22]Newren, E. P., Fogelson, A. L., Guy, R. D., and Kirby, R. M.. Unconditionally stable discretizations of the immersed boundary equations. J. Comput. Phys., 222(2):702–719, 2007.Google Scholar
[23]Peskin, C. S.. Numerical analysis of blood flow in the heart. J. Comput. Phys., 25(3):220–252, 1977.Google Scholar
[24]Pozrikidis, C.. Boundary Integral and Singularity Methods for Linearized Viscous Flow, volume 7. Cambridge University Press, Cambridge, 1992.Google Scholar
[25]Schmid-Schönbein, G. W., Sung, K. L., Tözeren, H., Skalak, R., and Chien, S.. Passive mechanical properties of human leukocytes. Biophys. J., 36(1):243–56, 1981.Google Scholar
[26]Stamenović, D.. Rheological behavior of mammalian cells. Cell. Mol. Life. Sci., 65(22):3592–3605, 2008.Google Scholar
[27]Stockie, J. M. and Wetton, B. R.. Analysis of stiffness in the immersed boundary method and implications for time-stepping schemes. J. Comput. Phys, 154:41–64, 1998.Google Scholar
[28]Strychalski, W. and Guy, R. D.. A computational model of bleb formation. Math. Med. Biol. in press, 2012.Google Scholar
[29]Teran, J., Fauci, L., and Shelley, M.. Viscoelastic fluid response can increase the speed and efficiency of a free swimmer. Phys. Rev. Lett., 104(3):038101, 2010.Google Scholar
[30]Theret, D. P., Levesque, M. J., Sato, M., Nerem, R. M., and Wheeler, L. T.. The application of a homogeneous half-space model in the analysis of endothelial cell micropipette measurements. J. Biomech. Eng., 110(3):190–199, 1988.Google Scholar
[31]Tu, C. and Peskin, C. S.. Stability and instability in the computation of flows with moving immersed boundaries: A comparison of three methods. SIAM J. Sci. Stat. Comput., 13:1361–1376, 1992.Google Scholar
[32]Wottawah, F., Schinkinger, S., Lincoln, B., Ananthakrishnan, R., Romeyke, M., Guck, J., and Kas, J.. Optical rheology of biological cells. Phys. Rev. Lett., 94(9), 2005.Google Scholar
[33]Yeung, A. and Evans, E.. Cortical shell-liquid core model for passive flow of liquid-like spherical cells into micropipets. Biophys. J., 56(1):139–149, 1989.Google Scholar