BZPEERReading, trust and collaboration for research
Menu

Discover

SearchFind papers, people, topics and opportunities from one place.GraphOpen maps around papers, topics and saved views.

Workspaces

HomeOpen the ranked research feed.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Source author record

Suncica Canic

Suncica Canic appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.APphysics.flu-dyn

Catalog footprint

What is connected

9works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2016arXiv

Comparison of reduced models for blood flow using Runge-Kutta discontinuous Galerkin methods

One-dimensional blood flow models take the general form of nonlinear hyperbolic systems but differ greatly in their formulation. One class of models considers the physically conserved quantities of mass and momentum, while another class describes mass and velocity. Further, the averaging process employed in the model derivation requires the specification of the axial velocity profile; this choice differentiates models within each class. Discrepancies among differing models have yet to be investigated. In this paper, we systematically compare several reduced models of blood flow for physiologically relevant vessel parameters, network topology, and boundary data. The models are discretized by a class of Runge-Kutta discontinuous Galerkin methods.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Existence of a weak solution to a fluid-elastic structure interaction problem with the Navier slip boundary condition

We study a nonlinear, moving boundary fluid-structure interaction problem between an incompressible, viscous Newtonian fluid, modeled by the 2D Navier-Stokes equations, and an elastic structure modeled by the shell or plate equations. The fluid and structure are coupled via the {\em Navier slip boundary condition} and balance of contact forces at the fluid-structure interface. The slip boundary condition is more realistic than the classical no-slip boundary condition in situations, e.g., when the structure is "rough", and in modeling dynamics near, or at a contact. Cardiovascular tissue and cell-seeded tissue constructs, which consist of grooves in tissue scaffolds that are lined with cells, are examples of "rough" elastic interfaces interacting with and incompressible, viscous fluid. The problem of heart valve closure is an example of a fluid-structure interaction problem with a contact. We prove the existence of a weak solution to this class of problems by designing a constructive proof based on the time discretization via operator splitting. This is the first existence result for fluid-structure interaction problems involving elastic structures satisfying the Navier slip boundary condition

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A partitioned scheme for fluid-composite structure interaction problems

We present a loosely-coupled partitioned scheme for a benchmark problem in fluid-composite structure interaction. The benchmark problem proposed here consists of an incompressible, viscous fluid interacting with a composite structure that consists of two layers: a thin elastic layer which is in contact with the fluid, and a thick elastic layer which sits on top of the thin layer. The motivation comes from fluid-structure interaction (FSI) in hemodyam- ics. The equations of linear elasticity are used to model the thick structural layer, while the Koiter member/shell equations are used to model the thin structural layer which serves as fluid-structure interface with mass. An effi- cient, modular, operator-splitting scheme is proposed to simulate solutions to the coupled, nonlinear FSI problem. The operator splitting scheme sepa- rates the elastodynamics structure problem, from a fluid problem in which the thin structure inertia is included as a Robin-type boundary condition to achieve unconditional stability, without requiring any sub-iterations within time-steps. An energy estimate associated with unconditional stability is derived for the fully nonlinear FSI problem defined on moving domains. Two instructive numerical examples are presented to test the performance of the scheme, where it is shown numerically, that the scheme is at least first-order accurate in time. The second example reveals a new phenomenon in FSI problems: the presence of a thin fluid-structure interface with mass regularizes solutions to the full FSI problem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Runge-Kutta Discontinuous Galerkin Method for Traffic Flow Model on Networks

We propose a bound-preserving Runge-Kutta (RK) discontinuous Galerkin (DG) method as an efficient, effective and compact numerical approach for numerical simulation of traffic flow problems on networks, with arbitrary high order accuracy. Road networks are modeled by graphs, composed of a finite number of roads that meet at junctions. On each road, a scalar conservation law describes the dynamics, while coupling conditions are specified at junctions to define flow separation or convergence at the points where roads meet. We incorporate such coupling conditions in the RK DG framework, and apply an arbitrary high order bound preserving limiter to the RK DG method to preserve the physical bounds on the network solutions (car density). We showcase the proposed algorithm on several benchmark test cases from the literature, as well as several new challenging examples with rich solution structures. Modeling and simulation of Cauchy problems for traffic flows on networks is notorious for lack of uniqueness or (Lipschitz) continuous dependence. The discontinuous Galerkin method proposed here deals elegantly with these problems, and is perhaps the only realistic and efficient high-order method for network problems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stability of the kinematically coupled β-scheme for fluid-structure interaction problems in hemodynamics

It is well-known that classical Dirichlet-Neumann loosely coupled partitioned schemes for fluid-structure interaction (FSI) problems are unconditionally unstable for certain combinations of physical and geometric parameters that are relevant in hemodynamics. It was shown in \cite{causin2005added} on a simple test problem, that these instabilities are associated with the so called ``added-mass effect''. By considering the same test problem as in \cite{causin2005added}, the present work shows that a novel, partitioned, loosely coupled scheme, recently introduced in \cite{MarSun}, called the kinematically coupled $β$-scheme, does not suffer from the added mass effect for any $β\in [0,1]$, and is unconditionally stable for all the parameters in the problem. Numerical results showing unconditional stability are presented for a full, nonlinearly coupled benchmark FSI problem, first considered in \cite{formaggia2001coupling}.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Modular, Operator Splitting Scheme for Fluid-Structure Interaction Problems with Thick Structures

We present an operator-splitting scheme for fluid-structure interaction (FSI) problems in hemodynamics, where the thickness of the structural wall is comparable to the radius of the cylindrical fluid domain. The equations of linear elasticity are used to model the structure, while the Navier-Stokes equations for an incompressible viscous fluid are used to model the fluid. The operator splitting scheme, based on Lie splitting, separates the elastodynamics structure problem, from a fluid problem in which structure inertia is included to achieve unconditional stability. We prove energy estimates associated with unconditional stability of this modular scheme for the full nonlinear FSI problem defined on a moving domain, without requiring any sub-iterations within time steps. Two numerical examples are presented, showing excellent agreement with the results of monolithic schemes. First-order convergence in time is shown numerically. Modularity, unconditional stability without temporal sub-iterations, and simple implementation are the features that make this operator-splitting scheme particularly appealing for multi-physics problems involving fluid-structure interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Existence of a solution to a fluid-multi-layered-structure interaction problem

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Existence of a weak solution to a nonlinear fluid-structure interaction problem modeling the flow of an incompressible, viscous fluid in a cylinder with deformable walls

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 \cite{GioSun} to 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fluid-structure interaction in blood flow capturing non-zero longitudinal structure displacement

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 advantages of partitioned schemes, such as modularity, simple implementation, and low computational costs.

0 citations0 reviewsNo code linkedOpen access