We study a nonlinear, unsteady, moving boundary, fluid-structure interaction (FSI) problem arising in modeling blood flow through elastic and viscoelastic arteries. The fluid flow, which is driven by the time-dependent pressure data, is governed by 2D incompressible Navier-Stokes equations, while the elastodynamics of the cylindrical wall is modeled by the 1D cylindrical Koiter shell model. Two cases are considered: the linearly viscoelastic and the linearly elastic Koiter shell. The fluid and structure are fully coupled (2-way coupling) via the kinematic and dynamic lateral boundary conditions describing continuity of velocity (the no-slip condition), and balance of contact forces at the fluid-structure interface. We prove existence of weak solutions to the two FSI problems (the viscoelastic and the elastic case) as long as the cylinder radius is greater than zero.The proof is based on a novel semi-discrete, operator splitting numerical scheme, known as the kinematically coupled scheme, introduced in [32] to numerically solve the underlying FSI problems.The backbone of the kinematically coupled scheme is the well-known Marchuk-Yanenko scheme, also known as the Lie splitting scheme. We effectively prove convergence of that numerical scheme to a solution of the corresponding FSI problem.

No abstract

We present a new model and a novel loosely coupled partitioned numerical scheme modeling fluid-structure interaction (FSI) in blood flow allowing non-zero longitudinal displacement. Arterial walls are modeled by a linearly viscoelastic, cylindrical Koiter shell model capturing both radial and longitudinal displacement. Fluid flow is modeled by the Navier-Stokes equations for an incompressible, viscous fluid. The two are fully coupled via kinematic and dynamic coupling conditions. Our numerical scheme is based on a new modified Lie operator splitting that decouples the fluid and structure sub-problems in a way that leads to a loosely coupled scheme which is unconditionally stable. This was achieved by a clever use of the kinematic coupling condition at the fluid and structure sub-problems, leading to an implicit coupling between the fluid and structure velocities. The proposed scheme is a modification of the recently introduced "kinematically coupled scheme" for which the newly proposed modified Lie splitting significantly increases the accuracy. The performance and accuracy of the scheme were studied on a couple of instructive examples including a comparison with a monolithic scheme. It was shown that the accuracy of our scheme was comparable to that of the monolithic scheme, while our scheme retains all the main ad-$ vantages of partitioned schemes, such as modularity, simple implementation, and low computational costs.

The broad research thematic of flows on networks was addressed in recent years by many researchers, in the area of applied mathematics, with new models based on partial differential equations. The latter brought a significant innovation in a field previously dominated by more classical techniques from discrete mathematics or methods based on ordinary differential equations. In particular, a number of results, mainly dealing with vehicular traffic, supply chains and data networks, were collected in two monographs: Traffic flow on networks, AIMSciences, Springfield, 2006, and Modeling, simulation, and optimization of supply chains, SIAM, Philadelphia, 2010. The field continues to flourish and a considerable number of papers devoted to the subject is published every year, also because of the wide and increasing range of applications: from blood flow to air traffic management. The aim of the present survey paper is to provide a view on a large number of themes, results and applications related to this broad research direction. The authors cover different expertise (modeling, analysis, numeric, optimization and other) so to provide an overview as extensive as possible. The focus is mainly on developments which appeared subsequently to the publication of the aforementioned books. 1 The author acknowledges partial support of 2013 GNAMPA project "Leggi di Conservazione: Teoria e Applicazioni". 2 The author acknowledges support by BMBF KinOpt, DFG Cluster of Excellence EXC128 and DAAD 54365630, 55866082. 3 The author acknowledges partial support of NSF Research Network in the Mathematical Sciences KI-Net "Kinetic description of emerging challenges in multiscale problems of natural sciences" Grant # : 1107444.

We study a nonlinear, unsteady, moving boundary, fluid-structure (FSI) problem in which the structure is composed of two layers: a thin layer which is in contact with the fluid, and a thick layer which sits on top of the thin structural layer. The fluid flow, which is driven by the time-dependent dynamic pressure data, is governed by the 2D Navier-Stokes equations for an incompressible, viscous fluid, defined on a 2D cylinder. The elastodynamics of the cylinder wall is governed by the 1D linear wave equation modeling the thin structural layer, and by the 2D equations of linear elasticity modeling the thick structural layer. The fluid and the structure, as well as the two structural layers, are fully coupled via the kinematic and dynamic coupling conditions describing continuity of velocity and balance of contact forces. The thin structural layer acts as a fluid-structure interface with mass. The resulting FSI problem is a nonlinear moving boundary problem of parabolic-hyperbolic type. This problem is motivated by the flow of blood in elastic arteries whose walls are composed of several layers, each with different mechanical characteristics and thickness. We prove existence of a weak solution to this nonlinear FSI problem as long as the cylinder radius is greater than zero. The proof is based on a novel semi-discrete, operator splitting numerical scheme, known as the kinematically coupled scheme. We effectively prove convergence of that numerical scheme to a solution of the nonlinear fluid-multi-layered-structure interaction problem. The spaces of weak solutions presented in this manuscript reveal a striking new feature: the presence of a thin fluid-structure interface with mass regularizes solutions of the coupled problem. 1 arXiv:1305.5310v1 [math.AP] 23 May 20133)