A new 3D immersed boundary method with application

Transcription

1 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 of Science & Technology of China

2 Talk Schedule 1) Background/motivation 2)Mathematical formulation 3)Numerical methods and algorithm 4)Verification/validation 5)Application 6)Summary 7)Future work

3 Immersed Boundary Methods C.S. Peskin, J. Comput. Phys. 25, pp.220 (1977) C.S. Peskin, Acta Numerica 11, 479 (2002). R. Mittal and G. Iaccarino, Annu. Rev. Fluid Mech. 37, pp , (2005). Immersed Interface Methods R.J. LeVeque and Z.L. Li, SIAM J. Numer. Anal. 31, p (1994). R.J. LeVeque and Z.L. Li, SIAM J. Sci. Comput. 18, p (1997). Z.L. Li, SIAM press, Philadelphia, (2006). Immersed Finite Element Methods Fictitious Domain Methods Arbitrary Lagrangian Eulerian Methods Phase Field Methods..

4 The immersed boundary method (Charles Peskin 1972) is a practical and effective method for problems of this kind. The boundary can be active (beating heart, swimming sperm) or passive (flag-in-wind) The boundary can be neutrally buoyant (swimming fish) or can have higher or lower density than surrounding fluids (aggregated RBCs in flowing blood) May be a body (eel, RBC) or a surface (flag, paper) May be open (flag) or closed (balloon) The boundary may be modeled by a collection of discrete elastic springs/fibers or by continuum/solid mechanics

5 Different versions of the IB method (incomplete list) Original version (Peskin 1972,1977,Peskin & McQueen 1993,1995,1996) Vortex-method version (McCracken & Peskin 1980) Volume-conserved version (Peskin & Printz 1993, Rosar & Peskin 2001) Adaptive mesh version (Roma, Peskin & Berger 1999, Griffth & Peskin 2007) Second-order version (Lai & Peskin 2000, Griffth & Peskin 2005) Multigrid version (Fogelson & Zhu, Zhu & Peskin 2002) Penalty version (Kim & Peskin 2006, Huang, Shin, Sung 2007, Huang Chang, Sung 2011) Stochastic version (Atzberger, Kramer & Peskin 2006) Kirchhoff rod version (Lim, Ferent, Wang & Peskin 2008) Viscoelastic fluid version (Chrispell, Cortez, Khismatullin, Fauci 2011, Chrispell, Fauci, Shelley 2013) Implicit version (Tu & Peskin 1992,Mayo & Peskin 1993, Fauci & Fogelson 1993,Taira & Colonius 2006,Mori & Peskin 2008, Hao & Zhu 2010,2011) Finite element version (B Griffith & XY Luo 2014) Lattice Boltzmann version (Zhu, He, et al. 2010)

6 Motivation FSI involving thin-walledstructures

7 Challenge of modeling thin-walled-structure by existing IB methods Elastic spring/fiber or 3D finite elements Lagrangian grid size for solid in the thickness direction is very small Accordingly the fluid grid has to be as small as the finest solid grid size to capture the thin yet finite thickness; may be too costly for computing in 3D A quick fix --- : a collection of thin flat plates model the thin-walled structure

8 Mathematical formulation of the IB method for fluid-thin-walled-structure interaction FLUID STRUCTURE ρ s Dd Dt = S FT + ρ s f s S -- 2 nd Piola-Kirchhoff stress tensor, F -- deformation gradient tensor

9 FSI by the penalty approach (Goldstein, Handler, Sirovich 1993, Kim & Peskin 2007, Huang, Shin, Sung 2007, Huang, Chang, Sung 2011) FSI inner B.C. Velocity: Force: Fluid virtual IB X IB real IB X compute interacting force by

10 existence and uniqueness of solutions of the mathematical formulation

11 Hybrid numerical methods (FEM + LBM) Lattice Boltzmann (D3Q19) for the Navier Stokes Finite element for the solid equations (Co-rotational scheme) Finite difference for other eqns

12 finite element method for thin plate motion large deformation large displacement large rotation nonlinearities in material & geometry small deformation (linear elasticity) co-rotational scheme for geometrical nonlinearity global local

13 FEM for in-plane & out-of-plane motions stiff matrix/load vector are computed in local coordinates in-plane motion by MRT (u,v,ωz ) out-of-plane motion by DKT (w,ωx,ωy) stiff matrix & load vector are transformed to and assembled in global coordinates by change of coordinates via transformation matrix is formed by element triads =

14 Newmark for time discretization Assembly stiff matrix & load vector in global coordinate system semi-discretization by FEM Newmark method for ODE system Newton s method for the nonlinear algebraic system LDL for linear algebraic system of each Newton step A=LDL

15 N-S eqns by lattice Boltzmann method A) Why LB method for Navier-Stokes? 1) Lattice Boltzmann method solves a series of linear scalar differential equations 2) Relationship between pressure and density 3) Clear physical interpretation of the scheme 4) Natural for parallelization 5)easier to model extra physics in a flow problem B) A brief introduction to the LB method C) Lattice-Boltzmann based IB method

16 Some References for the LBM Books on the Lattice Boltzmann Method: D. A. Wolf-Gladrow,``Lattice-gas cellular automata and lattice Boltzmann Models -- an introduction'', Springer, Berlin, (2000). S. Succi, ``The lattice Boltzmann equation'', Oxford Univ Press, Oxford (2001). M.C. Sukop and D.T. Thorne, Jr., ``Lattice Boltzmann Modeling: an introduction for geoscientists and engineers, Springer, Berlin, (2006). Zhaoli Guo and Cuguang Zheng, ``Theory and applications of Lattice Boltzmann method, Chinese Science Publisher, Beijing (2008). Widely cited paper on the LBM: L.-S. Luo. ``Unified Theory of the lattice Boltzmann models for nonideal gases, Phys. Rev. Lett. 81: 1618 (1998). S.Y. Chen, G.D. Doolen, ``Lattice Boltzmann Method for fluid flows'', Annu Rev. Fluid Mech., 30, p329, (1998). LBM is a fast growing area

17

18 Two widely used lattice Boltzmann models (left: D2Q9 right:d3q19)

19

20

21 Physical interpretation of lattice Boltzmann method 1) Collision 2) Streaming 3)Bounce-back BC

22 Literature on using LB method in the IB method Z.G. Feng and E.E. Michaelides, ``The immersed boundary-lattice Boltzmann method for solving fluid-particles interaction problems'', J. Comput.Phys. 195, (2004). Z.G. Feng and E.E. Michaelides, ``Proteus: a direct forcing method in the simulations of particulate flows'', J. Comput. Phys. 202, (2005). X.D. Niu, C. Shu, Y.T. Chew and Y. Peng,``A momentum exchange-based immersed boundary-lattice Boltzmann method for simulating incompressible viscous flows'', Physics Letters A, 354, p (2006). Y.Peng, C. Shu, Y.T. Chew, X.D. Niu, and X.Y. Lu,``Application of multi-block approach in the immersed boundary-lattice Boltzmann method for viscous fluid flows, J. Comp. Phys., 218, p (2006). Y. Sui, Y.T. Chew, P. Roy and H.T. Low,``A hybrid immersed boundary and multi-block lattice Boltzmann method for simulating fluid and moving-boundaries interactions, Int. J. Numer. Meth. Fluids, 53, p (2007). Y. Peng and L.-S. Luo, A comparative study of immersed-boundary and interpolated bounce-back methods in LBE, Prog. in Comput. Fluid. Dyn. 8 (1-4), pp (2008). Zhu, He, et al. 2010, Hao & Zhu 2010,2011 (both explicit and implicit) X Wang, C. Shu, J. Wu & LM Yang, Computer & Fluids, 100, pp (2014)

23 Some features of our new LB-IB method 1) Our new method is formulated for a thin-walled structure, and the motion of the structure is described by thin-plate eqns. 2) The fluid equations are solved by the D3Q19 LB model and structure equations are solved by an especial finite element method. Other equations in the formulation are solved by the finite difference methods. 3) The fluid-structure interaction forces are computed by the penalty approach.

24 Lattice Boltzmann (D3Q19) for N-S eqns

25 Summary: algorithm (from t to t+1) X(α, t+1) & Vs(α, t+1) next time step t+2

26 Verification/Validation in lieu of numerical analysis Comparision with Turek & Hron (2006)

27 Verification/Validation (II) Cd Rear vortex length flow past rigid plate Flow past a flag Re=100,M=1.0, K=0.0001,S=1000 rear edge bottom point rear edge middle point

28 Application Schouveiler, L. and C. Eloy (2013). "Flow-Induced Draping." Phys. Rev. Lett. 111:

29 Modeling and Simulation R Hua, L Zhu, X Lu, J Fluid Mech (2014), vol 759, pp Spatial: h=0.025 D Temporal: Re=100 & 200 M=0.1 & 2.0, d= K= ,S=

30 Deformation modes 6F 5F 4F 3F 2F 0F Folding number increases with the decreases of Kb Re=100,M=0.1,d=0.1,S=100

31 Deformation mode on S-K plane Re=100 Re=200 S increases, K decreases ====> n increases M=0.1,d=0.1

32 2 typical plate shapes and streamlines

33 2 typical plate shapes and vorticity contours

34 Vortex structure in the wake Re=100 Re=200 K=10-4 4F 4F K=10-4 K=10-2 2F 2F K=10-2 Re=100, steady for all modes Re=200,steady for 3F ~ 6F, unsteady for 2F and 0F M=0.1,d=0.1,S=100

35 Cd of the deformable plate Re=100 Re=200 M=0.1,d=0.1,S=100 Cd decreases with the decrease of K or increase of flexibility

36 Strain energy of the deformable plate Re=100 Re=200 M=0.1,d=0.1,S=100 Kb dominates Kp

37 Influence of fixed region size on deformation modes d=0.1 d=0.2 d=0.3 d=0.5 Re=100,M=0.1,S=100,K=10-4 n increases with d Computed n is consistent with theoretical n Cerda & Mahadevan 2003, PRL 90,

38 Influence of M on deformation mode M=0.1 (water tunnel) M=2.0 (wind tunnel) Re=100,d=0.1,S=100,K=10-4 M only prolongs the evolution process to the final stage, no effects on the final mode

39 Developing process of draping pattern Re=100, M=2.0, S=100, K=0.0001

40 Effect of Reynolds number on draping pattern M=2.0, S=1000, K= Re=

41 Summary 1) A new 3D IB method for fluid-thin-walled-structure-interaction is introduced where immersed structure is modeled by thin plate equations and solved by the co-rotational finite element scheme and the fluid equations are solved by the lattice Boltzmann method. The new method is verified and validated via comparisons with existing works. 2)The new method is applied to a 3D viscous flow past a deformable thin circular plate. Our main findings are: a) There exist six different plate deformation patterns; b)bending energy Kb dominates the plate strain energy; c) mass ratio M does not alter the deformation pattern of the plate; d)drag coefficient Cd decreases with the increase of structure flexibility; e)folding number increases with the size of the fixed region.

42 Ongoing and Future Work 1) Spectral/hp elements to model the thinwalled structures for extension to large deformation, displacement, and rotation, i.e. consider both material and geometry nonlinearities. (with S Dong and F Song) 2) Applications to Blood flows

43 Acknowledgement Thanks to supports from NSF of China ( ) and NSF of USA (DMS ). Thanks to my collaborators: Dr XY LU (USTC), Dr RN HUA, Thanks to JF Doyle for open source Fortran codes for MRT and DKT elements in his book: Doyle, James F. Nonlinear analysis of thin-walled structures: statics, dynamics, and stability. Springer Science & Business Media, Thanks to the workshop organizers and the invitation of Dr. Zhilin Li, the support from the IMS and Dr Sheng Xu via NSF conference grant.

44