Numerical analysis of the coupling of free fluid with a poroelastic material

This paper presents a numerical solution of the coupled system of the time-dependent Stokes and fully dynamic Biot equations. The numerical scheme is based on standard inf-sup stable finite elements in space and the Backward Euler scheme in time. We establish stability of the scheme and derive error estimates for the fully discrete coupled scheme. To handle realistic parameters which may cause nonphysical oscillations in the pore fluid pressure, a heuristic stabilization technique is considered. Numerical errors and convergence rates for smooth problems as well as tests on realistic material parameters are presented.     

A fluid‐fluid interaction method using decoupled subproblems and differing time steps

A time stepping procedure is proposed for a coupled fluid model motivated by the dynamic core of the atmosphere‐ocean system. The method exploits properties of the atmosphere‐ocean system to obtain efficiency. The momentum equations for the two fluids may be solved in parallel with different time step sizes. Stability is maintained with large time steps via a balanced two‐way passing of momentum flux. Numerical tests are provided that demonstrate the efficiency of the method. Published 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

Strong solutions to a nonlinear fluid structure interaction system

In this paper, we prove the existence of local-in-time smooth solutions to the nonlinear fluid structure interaction model first introduced in [J.-L. Lions, Quelques méthodes de résolution des problèmes aux limites non linéaires, Dunod, 1969] and considered in [V. Barbu, Z. Gruji?, I. Lasiecka, A. Tuffaha, Existence of the energy-level weak solutions for a nonlinear fluid-structure interaction model, in: Fluids and Waves, in: Contemp. Math., vol. 440, Amer. Math. Soc., Providence, RI, 2007, pp. 55-82; V. Barbu, Z. Gruji?, I. Lasiecka, A. Tuffaha, Smoothness of weak solutions to a nonlinear fluid-structure interaction model, Indiana Univ. Math. J. 57 (3) (2008) 1173-1207]. In particular, the strong solutions here are obtained given initial datum for the Navier-Stokes equation in the space H¹, and initial data for the wave equation w₀ and w₁ in the spaces H²(Ωe) and H¹(Ωe) respectively.     

Elasto-plastic earthquake response of arch dams including fluid–structure interaction by the Lagrangian approach

Elasto-plastic earthquake response of arch dams including fluid–structure interaction by the Lagrangian approach is mainly investigated in this study. To this aim, three-dimensional eight-noded version of Lagrangian fluid finite element including the effects of compressible wave propagation and surface sloshing motion, and three-dimensional version of Drucker–Prager model based on associated flow rule assumption were programmed in FORTRAN language by authors and incorporated into the program NONSAP. Two new components added into the NONSAP were tested on a simple fluid tank and a simple fluid–structure system and obtained very reasonable results.     

Local projection stabilized and characteristic decoupled scheme for the fluid–fluid interaction problems

In this article, we propose and analyse a local projection stabilized and characteristic decoupled scheme for the fluid–fluid interaction problems. We use the method of characteristics type to avert the difficulties caused by the nonlinear term, and use the local projection stabilized method to control spurious oscillations in the velocities due to dominant convection, and use a geometric averaging idea to decouple the monolithic problems. The stability analysis is derived and numerical tests are performed to demonstrate the robustness of this new method. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 704–723, 2017     

A FE‐BE coupling for a fluid‐structure interaction problem: Hierarchical a posteriori error estimates

In this article, we establish a hierarchical a posteriori error estimate for a coupling of finite elements and boundary elements for a fluid‐structure interaction problem posed in two and three dimensions. These methods combine boundary elements for the exterior fluid and finite elements for the elastic structure. We consider two weak formulations, a nonsymmetric one and a symmetric one, which are both uniquely solvable. We present the reliability and efficiency of the error estimates. For the two dimensional case, we compute local error indicators which allow us to develop an adaptive mesh refinement strategy on triangles. For the three dimensional case, we use hexahedrons as elements. Numerical experiments underline our theoretical results. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

Analysis of the operator splitting scheme for the Cahn‐Hilliard equation with a viscosity term

In this paper, we consider a second‐order fast explicit operator splitting method for the viscous Cahn‐Hilliard equation, which includes a viscosity term αΔu_t (α ∈ (0, 1)) described the influences of internal micro‐forces. The choice α = 0 corresponds to the classical Cahn‐Hilliard equation whilst the choice α = 1 recovers the nonlocal Allen‐Cahn equation. The fundamental idea of our method is to split the original problem into linear and nonlinear parts. The linear subproblem is numerically solved using a pseudo‐spectral method, and thus an ordinary differential equation is obtained. The nonlinear one is solved via TVD‐RK method. The stability and convergence are discussed in L²‐norm. Numerical experiments are performed to validate the accuracy and efficiency of the proposed method. Besides, a detailed comparison is made for the dynamics and the coarsening process of the metastable pattern for various values of α. Moreover, energy degradation and mass conservation are also verified.     

A micropolar fluid model of blood flow through a tapered artery with a stenosis

A nonlinear two‐dimensional micropolar fluid model for blood flow in a tapered artery with a single stenosis is considered. This model takes into account blood rheology in which blood consists of microelements suspended in plasma. The classical Navier–Stokes theory is inadequate to describe the microrotations or particles' spin of such suspension in a viscous medium. The governing equations involving unsteady nonlinear partial differential equations are solved using a finite difference scheme. A quantitative analysis performed through numerical computation shows that the axial velocity profile and the flow rate decrease and the wall shear stress increases once the artery is narrower in the presence of the polar effect. Furthermore, the taper angle certainly bears the potential to influence the velocity and the flow characteristics to considerable extent. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical simulation of three-dimensional incompressible fluid in a box flow passage considering fluid–structure interaction by differential quadrature method

A numerical simulation scheme of 3D incompressible viscous fluid in a box flow passage is developed to solve Navier–Stokes (N–S) equations by firstly taking fluid–structure interaction (FSI) into account. This numerical scheme with FSI is based on the polynomial differential quadrature (PDQ) approximation technique, in which motions of both the fluid and the solid boundary structures are well described. The flow passage investigated consists of four rectangular plates, of which two are rigid, while another two are elastic. In the simulation the elastic plates are allowed to vibrate subjected to excitation of the time-dependent dynamical pressure induced by the unsteady flow in the passage. Meanwhile, the vibrating plates change the flow pattern by producing many transient sources and sinks on the plates. The effects of FSI on the flow are evaluated by running numerical examples with the incoming flow’s Reynolds numbers of 3000, 7000 and 10,000, respectively. Numerical computations show that FSI has significant influence on both the velocity and pressure fields, and the DQ method developed here is effective for modelling 3D incompressible viscous fluid with FSI.     

An efficient ghost fluid method for compressible multifluids in Lagrangian coordinate

In this paper, we present an extended ghost fluid method (GFM) for computations of inviscid compressible multifluids in Lagrangian coordinate. That is, we capture the appropriate interface conditions by defining a fluid that has the velocity and the pressure of the real fluid at each point, but the entropy or the internal energy of some other fluid. Meanwhile, a single-fluid solver, CWENO-type central-upwind scheme, is developed in Lagrangian coordinate. The high resolution and the non-oscillatory quality of the scheme can be verified by solving several numerical experiments.     

Strong stability for a fluid–structure interaction model

In this paper we consider the question of stabilization of a fluid–structure model that describes the interaction between a 3‐D incompressible fluid and a 2‐D plate, the interface, which coincides with a flat flexible part of the surface of the vessel containing the fluid. The mathematical model comprises the Stokes equations and the equations for the longitudinal deflections of the plate with inclusion of the shear stress, which the fluid exerts on the plate. We show that the energy associated with the model decays strongly when the interface is equipped with a locally supported dissipative mechanism. Our main tool is an abstract resolvent criterion due to Tomilov. Copyright © 2008 John Wiley & Sons, Ltd.     

Existence of local strong solutions to fluid–beam and fluid–rod interaction systems

We study an unsteady nonlinear fluid–structure interaction problem. We consider a Newtonian incompressible two-dimensional flow described by the Navier–Stokes equations set in an unknown domain depending on the displacement of a structure, which itself satisfies a linear wave equation or a linear beam equation. The fluid and the structure systems are coupled via interface conditions prescribing the continuity of the velocities at the fluid–structure interface and the action-reaction principle. Considering three different structure models, we prove existence of a unique local-in-time strong solution, for which there is no gap between the regularity of the initial data and the regularity of the solution enabling to obtain a blow up alternative. In the case of a damped beam this is an alternative proof (and a generalization to non zero initial displacement) of the result that can be found in [20]. In the case of the wave equation or a beam equation with inertia of rotation, this is, to our knowledge the first result of existence of strong solutions for which no viscosity is added. The key points consist in studying the coupled system without decoupling the fluid from the structure and to use the fluid dissipation to control, in appropriate function spaces, the structure velocity.     

A regularity result for a solid–fluid system associated to the compressible Navier–Stokes equations

In this paper we deal with a fluid-structure interaction problem for a compressible fluid and a rigid structure immersed in a regular bounded domain in dimension 3. The fluid is modelled by the compressible Navier–Stokes system in the barotropic regime with no-slip boundary conditions and the motion of the structure is described by the usual law of balance of linear and angular moment.     

On the Euler implicit/explicit iterative scheme for the stationary Oldroyd fluid

In this article, we consider the stationary Oldroyd fluid equations from the large time behavior research of the nonstationary equations. Thus, to obtain its numerical solution, we first solve the nonstationary Oldroyd fluid equations via the Euler implicit/explicit finite element method with the integral term discretized by the right‐hand rectangle rule, then increase the total time (i.e., number of time steps) to approximate the solution of the original stationary equations. Under a new uniqueness condition (A2), we prove the exponential stability of the solution pair for the stationary equations and the almost unconditional stability of the numerical method. Furthermore, we also obtain the uniform optimal and error estimates in time integral . Finally, several numerical experiments are provided to verify our theoretical results.     

Existence of weak solutions for an interaction problem between an elastic structure and a compressible viscous fluid

We prove an existence result of weak solutions for an interaction problem between an elastic structure and a compressible fluid in three space dimensions. Solutions are defined as long as there is no collision and as long as conditions of non-interpenetration and of preservation of orientation are satisfied by the displacement field of the structure. To cite this article: M. Boulakia, C. R. Acad. Sci. Paris, Ser. I 340 (2005).     

New lax pair for restricted multiple three wave interaction system, quasiperiodic solutions and bi-Hamiltonian structure

We study restricted multiple three wave interaction system by the inverse scattering method. We develop the algebraic approach in terms of classical r-matrix and give an interpretation of the Poisson brackets as linear r-matrix algebra. The solutions are expressed in terms of polynomials of theta functions. In particular case for n = 1 in terms of Weierstrass functions.      

Existence of weak solutions for an interaction problem between an elastic structure and a compressible viscous fluid

In this article, we are interested by the three-dimensional motion of an elastic structure immersed in a viscous compressible fluid. The fluid and the structure are contained in a fixed bounded set. To describe the structure motion, we choose an Eulerian point of view and we strongly regularize the equation of the solid motion in order to get additional estimates on the elastic deformations. Our maim result is an existence result of weak solutions defined as long as no collisions occur and as long as conditions of non-interpenetration and of preservation of orientation are satisfied.     

An operator approach for an oblique derivative boundary‐transmission problem

We consider a boundary‐transmission problem for the Helmholtz equation, in a Bessel potential space setting, which arises within the context of wave diffraction theory. The boundary under consideration consists of a strip, and certain conditions are assumed on it in the form of oblique derivatives. Operator theoretical methods are used to deal with the problem and, as a consequence, several convolution type operators are constructed and associated to the problem. At the end, the well‐posedness of the problem is shown for a range of non‐critical regularity orders of the Bessel potential spaces, which include the finite energy norm space. In addition, an operator normalization method is applied to the critical orders case. Copyright © 2004 John Wiley & Sons, Ltd.