Computing Viscous Flow in an Elastic Tube  

作　　者：Yi Li Ioannis Sgouralis Anita T.Layton 

机构地区：[1]Department of Mathematics,Duke University,Durham,NC,USA

出　　处：《Numerical Mathematics(Theory,Methods and Applications)》2014年第4期555-574,共20页高等学校计算数学学报（英文版）

基　　金：supported in part by the National Institutes of Health,National Institute of Diabetes and Digestive and Kidney Diseases,through grant DK089066(to Layton).

摘　　要：We have developed a numerical method for simulating viscous flow through a compliant closed tube,driven by a pair of fluid source and sink.As is natural for tubular flow simulations,the problem is formulated in axisymmetric cylindrical coordinates,with fluid flow described by the Navier-Stokes equations.Because the tubular walls are assumed to be elastic,when stretched or compressed they exert forces on the fluid.Since these forces are singularly supported along the boundaries,the fluid velocity and pressure fields become unsmooth.To accurately compute the solution,we use the velocity decomposition approach,according to which pressure and velocity are decomposed into a singular part and a remainder part.The singular part satisfies the Stokes equations with singular boundary forces.Because the Stokes solution is unsmooth,it is computed to second-order accuracy using the immersed interface method,which incorporates known jump discontinuities in the solution and derivatives into the finite difference stencils.The remainder part,which satisfies the Navier-Stokes equations with a continuous body force,is regular.The equations describing the remainder part are discretized in time using the semi-Lagrangian approach,and then solved using a pressure-free projection method.Numerical results indicate that the computed overall solution is secondorder accurate in space,and the velocity is second-order accurate in time.

关 键 词：Navier-Stokes flow AXISYMMETRY immersed boundary fluid-structure interaction 

分 类 号：O17[理学—数学]