Accuracy Analysis of the Immersed Boundary Method Using the Method of Manufactured Solutions

GONG Zhao-xin; LU Chuan-jing; HUANG Hua-xiong

doi:10.3879/j.issn.1000-0887.2010.10.001

Volume 31 Issue 10

Turn off MathJax

Article Contents

Article Navigation > Applied Mathematics and Mechanics > 2010 > 31(10): 1141-1151

GONG Zhao-xin, LU Chuan-jing, HUANG Hua-xiong. Accuracy Analysis of the Immersed Boundary Method Using the Method of Manufactured Solutions[J]. Applied Mathematics and Mechanics, 2010, 31(10): 1141-1151. doi: 10.3879/j.issn.1000-0887.2010.10.001

Citation:

PDF( 3231 KB)

Accuracy Analysis of the Immersed Boundary Method Using the Method of Manufactured Solutions

doi: 10.3879/j.issn.1000-0887.2010.10.001

1.
Department of Mechanics, Shanghai Jiaotong University, Shanghai 200240, P. R. China;

Received Date: 1900-01-01
Rev Recd Date: 2010-06-17
Publish Date: 2010-10-15

Abstract

Abstract

The immersed boundary method is an effective technique for modeling and simulating fluid-structure interactions especially in the area of biomechanics.The principle work was to analyze the accuracy of the immersed boundary method.The whole procedure contained two parts:code verification and accuracy analysis.The code verification part provided the confidence that the code used here was mistake free and the accuracy analysis part gave the order of accuracy of the immersed boundary method.The method of manufactured solutions was taken as the research means for both parts.In the first part, the numerical code employed a second order discretization scheme, i.e., the theoretical order of accuracy was second order.It was matched by the calculated order of accuracy obtained by numerical calculation for all variables.This meant that the code contained no mistake, which was the premise of the following work.The second part introduced a jump in the manufactured solution for pressure and added corresponding singular forcing terms in the momentum equations.By analyzing the discretization errors, the accuracy of the immersed boundary method was proved to be first order, even though the discretization scheme was second order.In the meantime, it was found that the coarser mesh might be not sensitive enough to capture the influence of the immersed boundary and that refinement on the Lagrangian markers barely had any effects on the numerical calculation.
- manufactured solutions,
- immersed boundary method,
- order of accuracy,
- codeveri-fication,
- discretization error

FullText(HTML)

References(29)

References

[1]	Peskin C S. Numerical analysis of blood flow in the heart[J]. J Comput Phys，1977, 25(3): 220-252. doi: 10.1016/0021-9991(77)90100-0
[2]	Peskin C S， McQueen D M. A three-dimensional computational method for blood flow in the heart 1—immersed elastic fibers in a viscous incompressible fluid[J]. J Comput Phys，1989, 81(2): 372-405. doi: 10.1016/0021-9991(89)90213-1
[3]	Dillion R, Fauci L J, Gaver D. A microscale model of bacteria swimming, chemotaxis, and substrate transport[J]. J Theor Biol，1995, 177(4):325-340. doi: 10.1006/jtbi.1995.0251
[4]	Fauci L J，McDonald A. Sperm motility in the presence of boundaries[J]. B Math Biol， 1995, 57(5):679-699.
[5]	Fauci L J，Peskin C S. A computational model of aquatic animal locomotion[J]. J Comput Phys，1988, 77(1): 85-108. doi: 10.1016/0021-9991(88)90158-1
[6]	Bottino D C. Modeling viscoelastic networks and cell deformation in the context of the immersed boundary method[J]. J Comput Phys，1998, 147(1):86-113. doi: 10.1006/jcph.1998.6074
[7]	Fogelson A L. Continumm models of platelet aggregation: formulation and mechanical properties[J]. SIAM J Appl Math， 1992, 52(4):1089-1110. doi: 10.1137/0152064
[8]	Fauci L J，Fogelson A L. Truncated Newton’s methods and the modeling of complex immersed elastic structures[J]. Comm Pur Appl Math，1993, 46(6):787-818. doi: 10.1002/cpa.3160460602
[9]	Eggleton C D，Popel A S. Large deformation of red blood cell ghosts in a simple shear flow[J]. Phys Fluids， 1998, 10(8):1834-1845. doi: 10.1063/1.869703
[10]	Tu C，Peskin C S. Stability and instability in the computation of flows with moving immersed boundaries: a comparison of three methods[J]. SIAM J Sci Stat Comput， 1992, 13(6):1361-1376. doi: 10.1137/0913077
[11]	Stockie J M，Wetton B R. Analysis of stiffness in the immersed boundary method and implications for time-stepping schemes[J]. J Comput Phys，1999, 154(1): 41-64. doi: 10.1006/jcph.1999.6297
[12]	Stockie J M，Wetton B R. Stability analysis for the immersed fiber problem[J]. SIAM J Appl Math，1995, 55(6): 1577-1591. doi: 10.1137/S0036139994267018
[13]	Gong Z X, Huang H X，Lu C J. Stability analysis of the immersed boundary method for a two-dimensional membrane with bending rigidity[J]. Communications in Computational Physics, 2008, 3(3): 704-723.
[14]	Beyer R P，Leveque R J. Analysis of a one-dimensional model for the immersed boundary method[J]. SIAM J Num Anal，1992, 29(2):332-364. doi: 10.1137/0729022
[15]	Lai M C，Peskin C S. An immersed boundary method with formal second-order accuracy and reduced numerical viscosity[J]. J Comput Phys, 2000, 160(2): 705-719. doi: 10.1006/jcph.2000.6483
[16]	Griffith B E, Peskin C S. On the order of accuracy of the immersed boundary method: higher order convergence rates for sufficiently smooth problems[J], J Comput Phys, 2005, 208(1): 75-105.
[17]	Lai M C. Simulations of the flow past an array of circular cylinders as a test of the immersed boundary method[D]. Ph D thesis.New York: Courant Institute of Mathematical Sciences, New York University, 1998.
[18]	Leveque R J，Li Z. The immersed interface method for elliptic equations with discontinuous coefficients and singular sources[J]. SIAM J Numer Anal, 1994, 31(4): 1019-1044. doi: 10.1137/0731054
[19]	Steinberg S，Roache P J. Symbolic manipulation and computation fluid dynamics[J]. J Comput Phys, 1985, 57(2):251-284. doi: 10.1016/0021-9991(85)90045-2
[20]	Roache P J. Verification and Validation in Computational Science and Engineering[M]. Albuquerque, NM: Hermosa Publishers, 1998.
[21]	Roache P J. Code verification by the method of manufactured solutions [J]. J Fluid Eng，2002, 124(1):4-10. doi: 10.1115/1.1436090
[22]	Oberkampf W L, Trucano T G. Validation methodology in computational fluid dynamics[R]. AIAA Paper 2000-2549, Denver, CO. 19-22, June, 2000.
[23]	Oberkampf W L, Trucano T G. Verification and validation in computational fluid dynamics[J]. AIAA Progress in Aerospace Sciences, 2002, 38(3): 209-272. doi: 10.1016/S0376-0421(02)00005-2
[24]	Roy C J, Nelson C C, Smith T M，Ober C C. Verification of Euler/Navier-Stokes codes using the method of manufactured solutions[J]. Int J Numer Mech Fluids, 2004, 44(6):599-620. doi: 10.1002/fld.660
[25]	Bond R B, Ober C C，Knupp P M. A manufactured solution for verifying CFD boundary conditions—part Ⅲ[C]36th AIAA Fluid Dynamics Conference and Exhibit. San Francisco, CA, 2006, 1966-1982.
[26]	Brunner T A. Development of a grey nonlinear thermal radiation diffusion verification problem[J]. Transactions of the American Nuclear Society， 2006, 95:876-878.
[27]	Eca L, Hoekstra M, Hay A，Pelletier D. On the construction of manufactured solutions for one- and two-equation eddy-viscosity models[J]. Int J Num Mech Fluids， 2007, 54(2):119-154. doi: 10.1002/fld.1387
[28]	Tremblay D, Etienne S, Pelletier D. Code verification and the method of manufactured solutions for fluid-structure interaction problems[C]36th AIAA Fluid Dynamics Conference, San Francisco, CA. 2. 2006, 882-892.
[29]	Peskin C S. The immersed boundary method [J]. Acta Numerica，2002, 11:435-477. doi: 10.1017/CBO9780511550140.006