Research Papers: Flows in Complex Systems

Numerical and Experimental Analysis of the Fluid-Structure Interaction in Presence of a Hyperelastic Body

[+] Author and Article Information
H. Esmailzadeh

Mechanical Engineering Department,
Ferdowsi University of Mashhad,
Mashhad 91775-1111, Iran
e-mail: esmailzadeh_hamed@yahoo.com

M. Passandideh-Fard

Associate Professor
Mechanical Engineering Department,
Ferdowsi University of Mashhad,
Mashhad 91775-1111, Iran
e-mail: mpfard@um.ac.ir

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received July 28, 2013; final manuscript received June 18, 2014; published online September 4, 2014. Assoc. Editor: Zhongquan Charlie Zheng.

J. Fluids Eng 136(11), 111107 (Sep 04, 2014) (12 pages) Paper No: FE-13-1457; doi: 10.1115/1.4027893 History: Received July 28, 2013; Revised June 18, 2014

In this study, a numerical algorithm is developed for simulating the interaction between a fluid and a 2D/axisymmetric hyperelastic body based on a full Eulerian fluid-structure interaction (FSI) method. In this method, the solid volume fraction is used for describing the multicomponent material and the deformation tensor for describing the deformation of the hyperelastic body. The core elements of the simulation method are the constitutive law in the Cauchy stress form and an equation for the transport of the deformation tensor field. A semi-implicit formulation is used for the elastic stress to avoid instability especially for solid with high stiffness. The strain rate has a discontinuity across the fluid/solid interface. For improving the accuracy in capturing the interface, solid is treated as a highly viscous fluid. The viscosity term has the effect of smoothing the velocity and keeping the simulation stable. An experimental setup is used to validate the numerical results. The movement of a sphere made of silicone in air and its impact on a rigid substrate are investigated. The images are captured using a high speed CCD camera and the image processing technique is employed to obtain the required data from the images. For all cases considered, the results are in good agreement with those of the experiment performed in this study and other numerical results reported in the literature.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Belytschko, T., and Kennedy, J., 1975, “Finite Element Approach to Pressure Wave Attenuation by Reactor Fuel Subassemblies,” ASME J. Pressure Wave Technol., 97(3), pp. 172–177. [CrossRef]
Donea, J., Fasoli-Stella, P., and Giuliani, S., 1976, “Finite Element Solution of Transient Fluid-Structure Problems in Lagrangian Coordinates,” International Meeting on Fast Reactor Safety and Related Physics, Chicago, IL, Oct. 8, Vol. 3, pp. 1427–1435.
Hirt, C., Amsden, A., and Cook, J., 1974, “An Arbitrary Lagrangian–Eulerian Computing Method for all Flow Speeds,” J. Comput. Phys., 14(3), pp. 227–253. [CrossRef]
Nitikitpaiboon, C., and Bathe, K., 1993, “An Arbitrary Lagrangian–Eulerian Velocity Potential Formulation for Fluid-Structure Interaction,” Comput. Struct., 47(5), pp. 871–891. [CrossRef]
Figueroa, C., Vignon, I., Jansen, K., Hughes, T., and Taylor, C., 2006, “A Coupled Momentum Method for Modeling Blood Flow in Three-Dimensional Deformable Arteries,” J. Comput. Methods Appl. Mech. Eng., 195(41–43), pp. 1–22. [CrossRef]
Hu, H., 1996, “Direct Simulation of Flows of Solid-Liquid Mixtures,” Int. J. Multiphase Flow, 22(2), pp. 335–352. [CrossRef]
Gao, T., and Hu, H., 2009, “Deformation of Elastic Particles in Viscous Shear Flow,” J. Comput. Phys., 228, pp. 2132–2151. [CrossRef]
Haug, E. J., 1992, Intermediate Dynamics, Prentice Hall, Englewood Cliffs, NJ.
Wang, X., Wang, C., and Zhang, L. T., 2012, “Semi Implicit Formulation of the Immersed Finite Element Method,” J. Comput. Mech., 49, pp. 421–430. [CrossRef]
Peskin, C., 1972, “Flow Patterns Around Heart Valves: A Numerical Method,” J. Comput. Phys., 10, pp. 252–271. [CrossRef]
Wang, X. S., and Liu, W. K., 2004, “Extended Immersed Boundary Method Using FEM and RKPM,” Comput. Methods Appl. Mech. Eng., 193, pp. 1305–1321. [CrossRef]
Gilmanov, A., and Acharya, S., 2008, “A Hybrid Immersed Boundary and Material Point Method for Simulating 3D Fluid-Structure Interaction Problems,” Int. J. Numer. Methods Fluids, 56, pp. 151–177. [CrossRef]
Le, D. V., Khoo, B. C., and Peraire, J., 2006, “An Immersed Interface Method for Viscous Incompressible Flows Involving Rigid and Flexible Boundaries,” J. Comput. Phys., 220, pp. 109–138. [CrossRef]
Glowinski, R., Pan, T., Hesla, T., and Joseph, D., 1999, “A Distributed Lagrange Multiplier/Fictitious Domain Method for Particulate Flows,” Int. J. Multiphase Flow, 25, pp. 755–794. [CrossRef]
Hoogstraten, P., Slaats, P., and Baijens, F., 1991, “An Eulerian Approach to the Finite Element Modeling of Neo-Hookean Rubber Material,” Appl. Sci. Res., 48, pp. 193–210. [CrossRef]
Liu, C., and Walkington, N., 2001, “An Eulerian Description of Fluids Containing Visco-Elastic Particles,” Arch. Ratio. Mech. Anal., 159, pp. 229–252. [CrossRef]
Sugiyama, K., Ii, S., Takeuchi, S., Takagi, S., and Matsumoto, Y., 2010, “Full Eulerian Simulation of Biconcave Neo-Hookean Particles in a Poiseuille Flow,” J. Comput. Mech., 46, pp. 147–157. [CrossRef]
Nagano, N., Sugiyama, K., Takeuchi, S., Ii, S., Takagi, S., and Matsumoto, Y., 2010, “Full Eulerian Finite Difference Simulation of Fluid Flow in Hyperelastic Wavy Channel,” J. Fluid Sci. Technol., 5(3), pp. 475–490. [CrossRef]
Sugiyama, K., Ii, S., Takeuchi, S., Takagi, S., and Matsumoto, Y., 2011, “A Full Eulerian Finite Difference Approach for Solving Fluid-Structure Coupling Problems,” J. Comput. Phys., 230, pp. 596–627. [CrossRef]
Ii, S., Sugiyama, K., Takeuchi, S., Takagi, S., and Matsumoto, Y., 2011, “An Implicit Full Eulerian Method for the Fluid-Structure Interaction Problem,” Int. J. Numer. Methods Fluids, 65, pp. 150–165. [CrossRef]
Sugiyama, K., Nagano, N., Takeuchi, S., Ii, S., Takagi, S., and Matsumoto, Y., 2011, “Particle-in-Cell Method for Fluid-Structure Interaction Simulations of Neo-Hookean Tube Flows,” Theor. Appl. Mech., 59, pp. 245–256.
Ii, S., Gong, X., Sugiyama, K., Wu, J., Huang, H., and Takagi, S., 2011, “A Full Eulerian Fluid-Membrane Coupling Method With a Smoothed Volume of Fluid Approach,” J. Comput. Phys., 69, pp. 1–33. [CrossRef]
Takagi, S., Sugiyama, K., Ii, S., and Matsumoto, Y., 2012, “A Review of Full Eulerian Methods for Fluid-Structure Interaction Problems,” ASME J. Appl. Mech., 79, pp. 1–18. [CrossRef]
Ii, S., Sugiyama, K., Takagi, S., and Matsumoto, Y., 2012, “A Computational Blood Flow Analysis in a Capillary Vessel Including Multiple Red Blood Cells and Platelets,” J. Biomech. Sci. Eng., 7(1), pp. 72–83. [CrossRef]
Mirzaii, I., and Passandideh-Fard, M., 2012, “Modeling Free Surface Flows in Presence of an Arbitrary Moving Object,” Int. J. Multiphase Flow, 39, pp. 216–226. [CrossRef]
Carlson, M., Mucha, P., Brooks, R., and Turk, G., 2002, “Melting and Flowing,” ACM SIGGRAPH, 2002 Symposium on Computer Animation, San Antonio, TX, July 21–22, pp. 167–174.
Youngs, D. L., 1984, “An Interface Tracking Method for a 3D Eulerian Hydrodynamics Code,” AWRE, Technical Report No. 44/2092/35.
Bonet, J., and Wood, R., 1997, Nonlinear Continuum Mechanics for Finite Element Analysis, 2nd ed., Cambridge University, Cambridge, New York.
Chang, Y., Bao, K., Zhu, J., and Wu, E., 2011, “High Viscosity Fluid Simulation Using Particle-Based Method,” IEEE International Symposium on Virtual Reality Innovation, Singapore, March 19–20, pp. 199–205.
Kershaw, D. S., 1978, “The Incomplete Cholesky-Conjugate Gradient Method for the Iterative Solution of Systems of Linear Equations,” J. Comput. Phys., 26(1), pp. 43–65. [CrossRef]
Sardarabadi, M., Passandideh-Fard, M., and Zeinali Heris, S., 2014, “Experimental Investigation of the Effects of Silica/Water Nanofluid on PV/T (Photovoltaic Thermal Units),” J. Energy, 66, pp. 264–272. [CrossRef]
Darkner, S., and Erleben, K., 2011, “A Hyperelasticity Method for Interactive Virtual Design of Hearing Aids,” J. Visual Comput., 27, pp. 645–653. [CrossRef]
Marckmann, G., and Verron, E., 2006, “Comparison of Hyperelastic Models for Rubber-Like Materials,” Rubber Chem. Technol., 79, pp. 835–858. [CrossRef]
Wang, X., and Zhang, L., 2010, “Interpolation Functions in the Immersed Boundary and Finite Element Methods,” J. Comput. Mech., 45, pp. 321–334. [CrossRef]
Zhao, H., Freund, J. B., and Moser, R. D., 2008, “A Fixed-Mesh Method for Incompressible Flow-Structure Systems With Finite Solid Deformations,” J. Comput. Phys., 227, pp. 3114–3140. [CrossRef]
Robinson, A., Schroeder, C., and Fedkiw, R., 2011, “A Symmetric Positive Definite Formulation for Monolithic Fluid-Structure Interaction,” J. Comput. Phys., 230, pp. 1547–1566. [CrossRef]


Grahic Jump Location
Fig. 1

Flowchart of the sequence of the computational cycle for velocity, deformation tensor, and solid volume fraction

Grahic Jump Location
Fig. 2

Schematic of the experimental setup for the motion and impact of a hyperelastic sphere onto a rigid substrate

Grahic Jump Location
Fig. 3

Iamge processing technique using the Matlab software to obtain the position of the center of sphere volume: (a) image of sphere and a certain piece with known dimension for calibration, (b) the image after being processed using a threshold to reveal the sphere and the calibration piece by reducing the color values into black and white, (c) the top and bottom pixels of the sphere to obtain the position of its center while moving in air, and (d) the detection of the fluid/solid and solid/substrate interfaces to obtain the center during the impact

Grahic Jump Location
Fig. 4

Schematic of a sphere during its motion in air and the initial and boundary conditions for the simulation

Grahic Jump Location
Fig. 5

The result of mesh refinement study for (a) the position of the center of the sphere volume versus time and (b) the axial velocity versus time. The mesh size is characterized based on the number of CPR of the solid.

Grahic Jump Location
Fig. 6

The evolution of the motion and impact of a hyperelastic sphere with a radius 19.2 mm and a density 1106 kg/m3 in air onto a rigid substrate from the present model for (a) cross-sectional images, (b) 3D views, and (c) experimental results performed in this study

Grahic Jump Location
Fig. 7

Comparison of the numerical results with those of the experiments for (a) the position of the center of sphere volume and (b) the axial velocity

Grahic Jump Location
Fig. 8

Comparison of the numerical results with those of the experiments for various initial positions and velocities of the sphere: (a) the position of center of the sphere volume and (b) the axial velocity

Grahic Jump Location
Fig. 9

The simulation results for a soft wall in a lid-driven cavity flow from (a) present study and (b) those of Wang et al. [34]. (Reprinted with permission from Springer Science and Business Media.)

Grahic Jump Location
Fig. 10

The simulation results for a deformable solid motion in a lid-driven cavity flow from (a) present study and (b) those of Zhao et al. [35]. (Reprinted with permission from Elsevier.)

Grahic Jump Location
Fig. 11

The simulation results for an oscillating disk in fluid from (a) present study and (b) those of Robinson et al. [36]. (Reprinted with permission from Elsevier.)

Grahic Jump Location
Fig. 12

The hyperelastic rectangle in a time-varying shear flow: (a) present study and (b) those of Sugiyama et al. [19]. (Reprinted with permission from Elsevier.)




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In