A kinetic scheme for unsteady pressurised flows in closed water pipes

The aim of this paper is to present a kinetic numerical scheme for the computations of transient pressurised flows in closed water pipes. Firstly, we detail the mathematical model written as a conservative hyperbolic partial differentiel system of equations, and the we recall how to obtain the corresponding kinetic formulation. Then we build the kinetic scheme ensuring an upwinding of the source term due to the topography performed in a close manner described by Perthame et al. using an energetic balance at microscopic level for the Shallow Water equations. The validation is lastly performed in the case of a water hammer in a uniform pipe: we compare the numerical results provided by an industrial code used at EDF-CIH (France), which solves the Allievi equation (the commonly used equation for pressurised flows in pipes) by the method of characteristics, with those of the kinetic scheme. It appears that they are in a very good agreement

Interior feedback stabilization of wave equations with dynamic boundary delay

In this paper we consider an interior stabilization problem for the wave equation with dynamic boundary delay.We prove some stability results under the choice of damping operator. The proof of the main result is based on a frequency domain method and combines a contradiction argument with the multiplier technique to carry out a special analysis for the resolvent

Global existence and exponential growth for a viscoelastic wave equation with dynamic boundary conditions

The goal of this work is to study a model of the wave equation with dynamic boundary conditions and a viscoelastic term. First, applying the Faedo-Galerkin method combined with the fixed point theorem, we show the existence and uniqueness of a local in time solution. Second, we show that under some restrictions on the initial data, the solution continues to exist globally in time. On the other hand, if the interior source dominates the boundary damping, then the solution is unbounded and grows as an exponential function. In addition, in the absence of the strong damping, then the solution ceases to exist and blows up in finite time.Comment: arXiv admin note: text overlap with arXiv:0810.101

Air entrainment in transient flows in closed water pipes: a two-layer approach

In this paper, we first construct a model for free surface flows that takes into account the air entrainment by a system of four partial differential equations. We derive it by taking averaged values of gas and fluid velocities on the cross surface flow in the Euler equations (incompressible for the fluid and compressible for the gas). The obtained system is conditionally hyperbolic. Then, we propose a mathematical kinetic interpretation of this system to finally construct a two-layer kinetic scheme in which a special treatment for the "missing" boundary condition is performed. Several numerical tests on closed water pipes are performed and the impact of the loss of hyperbolicity is discussed and illustrated. Finally, we make a numerical study of the order of the kinetic method in the case where the system is mainly non hyperbolic. This provides a useful stability result when the spatial mesh size goes to zero

A pseudo active kinematic constraint for a biological living soft tissue: an effect of the collagen network

Recent studies in mammalian hearts show that left ventricular wall thickening is an important mechanism for systolic ejection and that during contraction the cardiac muscle develops significant stresses in the muscular cross-fiber direction. We suggested that the collagen network surrounding the muscular fibers could account for these mechanical behaviors. To test this hypothesis we develop a model for large deformation response of active, incompressible, nonlinear elastic and transversely isotropic living soft tissue (such as cardiac or arteries tissues) in which we include a coupling effect between the connective tissue and the muscular fibers. Then, a three-dimensional finite element formulation including this internal pseudo-active kinematic constraint is derived. Analytical and finite element solutions are in a very good agreement. The numerical results show this wall thickening effect with an order of magnitude compatible with the experimental observations

A model for unsteady mixed flows in non uniform closed water pipes and a well-balanced finite volume scheme

We present the derivation of a new unidirectional model for We present the derivation of a new unidirectional model for unsteady mixed flows in non uniform closed water pipes. We introduce a local reference frame to take into account the local perturbation caused by the changes of section and slope. Then an asymptotic analysis is performed to obtain a model for free surface flows and another one for pressurized flows. By coupling these models through the transition points by the use of a common set of variables and a suitable pressure law, we obtain a simple formulation called PFS-model close to the shallow water equations with source terms. It takes into account the changes of section and the slope variation in a continuous way through transition points. Transition point between the two types of flows is treated as a free boundary associated to a discontinuity of the gradient of pressure. The numerical simulation is performed by making use of a Roe-like finite volume scheme that we adapted to take into account geometrical source terms in the convection matrix. Finally some numerical tests are presented

A kinetic scheme for transient mixed flows in non uniform closed pipes: a global manner to upwind all the source terms

We present a numerical kinetic scheme for an unsteady mixed pressurised and free surface model. This model has a source term depending on both the space variable and the unknown, U, of the system. The source term is composed by a topography, a section variation, a curvature (also called corrective) and a friction term. Using the Finite Volume and Kinetic (FVK) framework, we propose an approximation of the source terms following the principle of interfacial upwind with a kinetic interpretation: the source term is not treated as a volumic term, but included in the numerical fluxes. Then, several numerical tests are presented

Asymptotic error distribution for the Ninomiya-Victoir scheme in the commutative case

In a previous work, we proved strong convergence with order $1$ of the Ninomiya-Victoir scheme $X^{NV}$ with time step $T/N$ to the solution $X$ of the limiting SDE when the Brownian vector fields commute. In this paper, we prove that the normalized error process $N \left(X - X^{NV}\right)$ converges to an affine SDE with source terms involving the Lie brackets between the Brownian vector fields and the drift vector field. This result ensures that the strong convergence rate is actually $1$ when the Brownian vector fields commute, but at least one of them does not commute with the drift vector field. When all the vector fields commute the limit vanishes. Our result is consistent with the fact that the Ninomiya-Victoir scheme solves the SDE in this case.Comment: arXiv admin note: text overlap with arXiv:1601.0526

Ninomiya-Victoir scheme: strong convergence, antithetic version and application to multilevel estimators

In this paper, we are interested in the strong convergence properties of the Ninomiya-Victoir scheme which is known to exhibit weak convergence with order 2. We prove strong convergence with order $1/2$. This study is aimed at analysing the use of this scheme either at each level or only at the finest level of a multilevel Monte Carlo estimator: indeed, the variance of a multilevel Monte Carlo estimator is related to the strong error between the two schemes used on the coarse and fine grids at each level. Recently, Giles and Szpruch proposed a scheme permitting to construct a multilevel Monte Carlo estimator achieving the optimal complexity $O\left(\epsilon^{-2}\right)$ for the precision $\epsilon$. In the same spirit, we propose a modified Ninomiya-Victoir scheme, which may be strongly coupled with order $1$ to the Giles-Szpruch scheme at the finest level of a multilevel Monte Carlo estimator. Numerical experiments show that this choice improves the efficiency, since the order $2$ of weak convergence of the Ninomiya-Victoir scheme permits to reduce the number of discretization levels