Coupled Problems of Dynamics in Materially Incompressible Saturated Porous Media

The mechanical behavior of saturated porous materials is largely governed by the interaction between the solid skeleton and the pore fluid. This interaction is particularly strong in dynamic problems and leads to numerical challenges especially in the case of incompressible constituents. In fact, the permeability plays a significant role in this coupling and influences the choice of a proper time integration scheme. Proceeding from the macroscopic Theory of Porous Media (TPM) within the isothermal and geometrical linear regime, the governing balance equations of the dynamic binary solid–fluid model are the solid and fluid momentum balances, and the overall volume balance of the biphasic mixture. This set of coupled partial differential equations (PDEs) is solved within the framework of the mixed Finite Element Method (FEM) applying two different time solution methods, viz., a monolithic implicit and a splitted implicit–explicit scheme. The time stepping algorithms are implemented into the FE program PANDAS and a Scilab FE routine and compared on a one–dimensional wave propagation example. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

How an added mass matrix estimation may dramatically improve FSI calculations for moving foils

This paper presents a corrected partitioned scheme for investigating fluid–structure interaction (FSI) that may be encountered by lifting devices immersed in heavy fluid such as liquids. The purpose of this model is to counteract the penalizing impact of the added mass effect on the classical partitioned FSI coupling scheme. This work is based on an added mass corrected version of the classical strongly coupled partitioned scheme presented in Song et al. (2013). Results show that this corrected version systematically allows convergence to the coupled solution. The fluid flow model considered here uses a non-stationary potential approach, commonly termed the Panel Method. The advantage of this kind of approach is twofold: first, in restricting itself to a boundary method and, second, in allowing an added mass matrix to be estimated as a post-processing phase. Whereas the classical scheme encounters an acceptable (no numerical oscillation) convergence limit for fluid densities higher than 8 kg/m³ for the considered case, our corrected scheme is not dependent on fluid density and converges with only 6 iterations. This makes it possible to investigate the dynamic behavior of a 2D foil immersed in heavy fluids such as water. For example, it recognizes that frequency shifting may occur as the consequence of a strong added mass effect.     

A coupled NMM-SPH method for fluid-structure interaction problems

The main challenges in the numerical simulation of fluid–structure interaction (FSI) problems include the solid fracture, the free surface fluid flow, and the interactions between the solid and the fluid. Aiming to improve the treatment of these issues, a new coupled scheme is developed in this paper. For the solid structure, the Numerical Manifold Method (NMM) is adopted, in which the solid is allowed to change from continuum to discontinuum. The Smoothed Particle Hydrodynamics (SPH) method, which is suitable for free interface flow problem, is used to model the motion of fluids. A contact algorithm is then developed to handle the interaction between NMM elements and SPH particles. Three numerical examples are tested to validate the coupled NMM-SPH method, including the hydrostatic pressure test, dam-break simulation and crack propagation of a gravity dam under hydraulic pressure. Numerical modeling results indicate that the coupled NMM-SPH method can not only simulate the interaction of the solid structure and the fluid as in conventional methods, but also can predict the failure of the solid structure.     

Vibration of submerged floating tunnels due to moving loads

The concept of submerged floating tunnel (SFT) has become an increasingly attractive idea to cross the straits. The structural solution in this scheme includes buoyancy force on tunnel body plus tension in mooring tethers. This paper investigates the effect of submergence on the dynamic response of SFT due to moving load. The inertial effect of the fluid is accounted for by evaluating the added mass of tunnel using two and three dimensional models. It is found that fluid–structure interaction increases dynamic amplification of the tunnel deflection (in some cases very significantly). The results show that although the 3D model predicts lesser inertial contribution for surrounding fluid, it is not always possible to associate the larger response with the 2D or 3D models. The discrepancy between the results of the two models decreases as the tether stiffness increases. This indicates that the adoption of Morison’s equation for evaluating the fluid loading on the tunnel is a reasonable assumption when the tether stiffness is high. It is also found that by increasing the tether stiffness, it is possible to introduce a major reduction in the dynamic amplification of the response and by this way control the dynamic response of the SFT.     

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.     

Coupling of mixed finite element and stabilized boundary element methods for a fluid–solid interaction problem in 3D

We introduce and analyze the coupling of a mixed finite element and a boundary element for a three‐dimensional time‐harmonic fluid–solid interaction problem. We consider a formulation in which the Cauchy stress tensor and the rotation are the main variables in the elastic structure and use the usual pressure formulation in the acoustic fluid. The mixed variational formulation in the solid is completed with boundary integral equations relating the Cauchy data of the acoustic problem on the coupling interface. A crucial point in our formulation is the stabilization technique introduced by Hiptmair and coworkers to avoid the well‐known instability issue appearing in the boundary element method treatment of the exterior Helmholtz problem. The main novelty of this formulation, with respect to a previous approach, consists in reducing the computational domain to the solid media and providing a more accurate treatment of the far field effect. We show that the continuous problem is well‐posed and propose a conforming Galerkin method based on the lowest‐order Arnold–Falk–Winther mixed finite element. Finally, we prove that the numerical scheme is convergent with optimal order.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1211–1233, 2014     

A Lagrangian-based SPH-DEM model for fluid–solid interaction with free surface flow in two dimensions

A Lagrangian-based SPH-DEM coupling model is proposed to study fluid–solid interaction (FSI) problems with free-surface flow. In this model, SPH uses an incompressible divergence-free scheme for simulating complex flow problems. Based on the Mohr–Coulomb criterion with tension cut, the DEM describes the characteristics of solid deformation and failure by means of contact models between particles. The coupling mechanism between SPH and DEM is realised by the decoupling of the force field during the process of fluid–solid interaction. That is, the motions of fluid and solid particles are reflected by the Navier–Stokes equations and interactions among solid particles are determined by Newton's second law in the DEM. To demonstrate the applicability of the SPH-DEM model, three case studies are used to verify the different fluid interaction situations with rigid bodies, deformable objects, and granular assemblies, respectively. The results of the proposed model shows good agreement with experimental data and indicates that it is capable of capturing the features of solid movement, deformation and failure under complex flow conditions with convincing accuracy and high efficiency.     

计入气穴的水下爆炸作用下结构流固耦合三维数值模拟

水下爆炸在结构物面附近产生的气穴现象,严重影响水下爆炸作用下的流固耦合动响应,是舰船水下爆炸领域的难点,传统的边界元方法、有限元方法(FEM)难以解决水下爆炸气穴现象这类强非线性问题．针对此问题,计及流体中的气穴现象,考虑流体的可压缩型,忽略流体粘性,建立水下爆炸瞬态强非线性流固耦合三维数值模型,采用流体谱单元方法(SEM)和结构有限元方法求解该模型．计算结果表明:相对有限元法,谱单元法具有更高的计算精度,且谱单元解与解析解、试验值吻合良好．在此基础上,结合ABAQUS软件,分别探讨三维球壳、船体板架在水下爆炸作用下的瞬态流固耦合机理,给出气穴区域及其对水中结构物动响应的影响特征,旨在为舰船水下爆炸瞬态流固耦合问题的相关研究提供参考．     

Finite element method to fluid–solid interaction problems with unbounded periodic interfaces

Consider a time‐harmonic acoustic plane wave incident onto a doubly periodic (biperiodic) surface from above. The medium above the surface is supposed to be filled with a homogeneous compressible inviscid fluid of constant mass density, whereas the region below is occupied by an isotropic and linearly elastic solid body characterized by its Lamé constants. This article is concerned with a variational approach to the fluid–solid interaction problems with unbounded biperiodic Lipschitz interfaces between the domains of the acoustic and elastic waves. The existence of quasiperiodic solutions in Sobolev spaces is established at arbitrary frequency of incidence, while uniqueness is proved only for small frequencies or for all frequencies excluding a discrete set. A finite element scheme coupled with Dirichlet‐to‐Neumann mappings is proposed and the convergence analysis is performed. The Dirichlet‐to‐Neumann mappings are approximated by truncated Rayleigh series expansions. Finally, numerical tests in 2D are presented to confirm the convergence of solutions and the energy balance formula. In particular, the frequency spectrum of normally reflected signals is plotted for water–brass and water–brass–water interfaces. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 5–35, 2016     

2D simulation of fluid-structure interaction using finite element method

This paper deals with pressure-based finite element analysis of fluid–structure systems considering the coupled fluid and structural dynamics. The present method uses two-dimensional fluid elements and structural line elements for the numerical simulation of the problem. The equations of motion of the fluid, considered inviscid and compressible, are expressed in terms of the pressure variable alone. The solution of the coupled system is accomplished by solving the two systems separately with the interaction effects at the fluid–solid interface enforced by an iterative scheme. Non-divergent pressure and displacement are obtained simultaneously through iterations. The Galerkin weighted residual method-based FE formulation and the iterative solution procedure are explained in detail followed by some numerical examples. Numerical results are compared with the existing solutions to validate the code for sloshing with fluid–structure coupling.     

Numerical investigation of the behavior of thin flexible membranes in interaction with a fluid domain

With the increasing concern over massive damage ensuing from flood, the application of flood protection systems has aroused growing interest. As part of an interdisciplinary project regarding the development of a temporary installable textile protective system against flood, we study, how the thin flexible membrane behaves in interaction with a fluid domain numerically and experimentally. The numerical research was carried out by means of the FV solver available in Star CCM+ based on both morphing motion model and on the basis of overset mesh and morphing motion model, with the help of a FEM solver available in Abaqus and finally through FSI Star CCM+ and Abaqus Co-Simulation (interconnection between FV solver for fluid domain analysis and FEM solver for solid domain analysis). A key challenge facing us was the implementation of proper damping models to address the problem with oscillation of the membrane. We have examined, how the change in water level affects the behavior of the thin flexible membrane. A comparable physical validation with respect to experimental model has been achieved through numerical analysis. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A time discretization scheme of a characteristics method for a fluid–rigid system with discontinuous density

We propose a new characteristics method for the time discretization of a fluid–rigid system in the case when the densities of the fluid and the solid are different. This method is based on a global weak formulation involving only terms defined on the whole fluid–rigid domain. The main idea is to construct a characteristic function which preserves the rigidity of the solid at the discrete time levels. A convergence result for this semi-discrete scheme is then given.     

A projection algorithm for fluid–structure interaction problems with strong added-mass effect

This Note aims at introducing a semi-implicit coupling scheme for fluid–structure interaction problems with a strong added-mass effect. Our main idea relies on the splitting of added-mass, viscous effects and geometrical/convective non-linearities, through a Chorin–Temam projection scheme within the fluid. We state some theoretical stability results, in the linear case, and provide some numerical experiments. The main interest of the proposed scheme is its efficiency compared to the implicit approach. To cite this article: M.A. Fernández et al., C. R. Acad. Sci. Paris, Ser. I 342 (2006).     

A new approach to the stabilization of a fluid–structure interaction model

This article is a continuation of our work on a linear fluid–structure interaction model [Grobbelaar-Van Dalsen, On a fluid–structure model in which the dynamics of the structure involves the shear stress due to the fluid, J. Math. Fluid Mech. 10(3) (2008), pp. 388–401; Grobbelaar-Van Dalsen, Strong stability for a fluid––structure model, Math. Methods Appl. Sci., 32(2009) pp. 1452–1466]. The model 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 the inclusion of the shear stress that the fluid exerts on the plate. A dissipative damping mechanism of Kelvin–Voigt type is applied to the interior of the plate. While our earlier work shows that weak solutions in a space of finite energy are strongly asymptotically stable under no-slip transmission conditions at the interface with uniform exponential stability only attainable under an additional domination condition, the present research is directed at achieving uniform exponential stability of weak solutions without imposing the domination condition. Using energy methods we establish uniform exponential decay under a modified transmission condition at the interface. This condition entails that the fluid velocity at the interface is coupled to a linear combination of the plate velocity and displacement.     

Convergence analysis of an approximation to miscible fluid flows in porous media by combining mixed finite element and finite volume methods

This article deals with development and analysis of a numerical method for a coupled system describing miscible displacement of one incompressible fluid by another through heterogeneous porous media. A mixed finite element (MFE) method is employed to discretize the Darcy flow equation combined with a conservative finite volume (FV) method on unstructured grids for the concentration equation. It is shown that the FV scheme satisfies a discrete maximum principle. We derive L^∞ and BV estimates under an appropriate CFL condition. Then we prove convergence of the approximate solutions to a weak solution of the coupled system. Numerical results are presented to see the performance of the method in two space dimensions. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

An added-mass free semi-implicit coupling scheme for fluid–structure interaction

In this Note we propose a semi-implicit coupling scheme for the numerical simulation of fluid–structure interaction systems involving a viscous incompressible fluid. The scheme is stable irrespectively of the so-called added-mass effect and allows for conservative time-stepping within the structure. The efficiency of the scheme is based on the explicit splitting of the viscous effects and geometrical/convective non-linearities, through the use of the Chorin–Temam projection scheme within the fluid. Stability relies on the implicit treatment of the pressure stresses and on the Nitsche's treatment of the viscous coupling. To cite this article: M. Astorino et al., C. R. Acad. Sci. Paris, Ser. I 347 (2009).     

A kinematically coupled time-splitting scheme for fluid–structure interaction in blood flow

We present a new time-splitting scheme for the numerical simulation of fluid–structure interaction between blood flow and vascular walls. This scheme deals in a successful way with the problem of the added mass effect. The scheme is modular and it embodies the stability properties of implicit schemes at the low computational cost of loosely coupled ones.     

Dynamic fractal unifying interaction confirmed with magnetospheric behavior and orbital data

The solar wind makes the magnetosphere to expand and contract as indicated by the expansions and contractions of the auroral oval due to balancing of the dynamic pressure of the ambient space plasma at inner and outer magnetic lines. This self‐similar magnetospheric behavior elucidates the controversial magnetic storm‐substorm relationship and reveals the 3D‐spiral structure of magnetic interaction. The found self‐similarly evolving structure of one seen as fundamental interaction suggests dynamic fractal unifying interaction that builds a firework universe having 3D‐spiral code. The unifying interaction is described with equation drawn in new fundamental dynamic fractal framework. The equation of unifying interaction converges to the inverse square laws and the principle of uncertainty at laboratory scales. The dynamic fractal fundamental framework is made of one 3D‐spirally‐faster‐inward contracting and expanding, oscillating, basic matter. It simply accounts for observed constant speed of light and for the creation of bright and dark bands on a screen behind a tiny slit. The dynamic fractal framework is quantitatively confirmed with the orbital data for the Milky Way Galaxy, the Sun, the Earth, and the triple asteroid system 87 Sylvia. Many testable predictions are also made. The presented new fundamental dynamic fractal framework allows qualitative and quantitative modeling and simplification. © 2007 Wiley Periodicals, Inc. Complexity 12: 61–76, 2007     

Divergence-free finite elements for the numerical solution of a hydroelastic vibration problem

In this paper, we analyze a divergence-free finite element method to solve a fluid–structure interaction spectral problem in the three-dimensional case. The unknowns of the resulting formulation are the fluid and solid displacements and the fluid pressure on the interface separating both media. The resulting mixed eigenvalue problem is approximated by using appropriate basis of the divergence-free lowest order Raviart–Thomas elements for the fluid, piecewise linear elements for the solid and piecewise constant elements for the interface pressure. It is proved that eigenvalues and eigenfunctions are correctly approximated and some numerical results are reported in order to assess the performance of the method.