Fluid–solid interaction problems with thermal convection using the immersed element‐free Galerkin method

In this work, the immersed element‐free Galerkin method (IEFGM) is proposed for the solution of fluid–structure interaction (FSI) problems. In this technique, the FSI is represented as a volumetric force in the momentum equations. In IEFGM, a Lagrangian solid domain moves on top of an Eulerian fluid domain that spans over the entire computational region. The fluid domain is modeled using the finite element method and the solid domain is modeled using the element‐free Galerkin method. The continuity between the solid and fluid domains is satisfied by means of a local approximation, in the vicinity of the solid domain, of the velocity field and the FSI force. Such an approximation is achieved using the moving least‐squares technique. The method was applied to simulate the motion of a deformable disk moving in a viscous fluid due to the action of the gravitational force and the thermal convection of the fluid. An analysis of the main factors affecting the shape and trajectory of the solid body is presented. The method shows a distinct advantage for simulating FSI problems with highly deformable solids. Copyright © 2009 John Wiley & Sons, Ltd.     

A time splitting fictitious domain algorithm for fluid–structure interaction problems (A fictitious domain algorithm for FSI)

A Finite Element Method in mixed Eulerian and Lagrangian formulation is developed to allow direct numerical simulations of dynamical interaction between an incompressible fluid and a hyper-elastic incompressible solid. A Fictitious Domain Method is applied so that the fluid is extended inside the deformable solid volume and the velocity field in the entire computational domain is resolved in an Eulerian framework. Solid motion, which is tracked in a Lagrangian framework, is imposed through the body force acting on the fluid within the solid boundaries. Solid stress smoothing on the Lagrangian mesh is performed with the Zienkiewicz–Zhu patch recovery method. High-order Gaussian integration quadratures over cut elements are used in order to avoid sub-meshing within elements in the Eulerian mesh that are intersected by the Lagrangian grid. The algorithm is implemented and verified in two spatial dimensions by comparing with the well validated simulations of solid deformation in a lid driven cavity and periodic elastic wall deformation driven by a time-dependent flow. It shows good agreement with the numerical results reported in the literature. In 3-D the method is validated against previously reported numerical simulations of 3-D rhythmically contracting alveolated ducts.     

Numerical analysis of gas bubbles in close proximity to a movable or deformable body

The growth and collapse of gaseous bubbles near a movable or deformable body are investigated numerically using the boundary element method and fluid–solid coupling technique. The fluid is treated as inviscid, incompressible and the flow irrotational. The unsteady Bernoulli equation is applied on the bubble surface as one of the boundary conditions of the Laplace’s equation for the potential. Good agreements between the numerical and experimental results demonstrate the robustness and accuracy of the present method. The translation and rotation of the rigid body due to the bubble evolution are captured by solving the six-degrees-of-freedom equations of motion for the rigid body. The fluid–solid coupling is achieved by matching the normal component of the velocity and the pressure at the fluid–solid interface. Compared to a fixed rigid body, the expansion of the bubble is not affected too much but much faster collapsing velocities during the collapsing phase of bubble can be observed when considering the motion of the rigid body. The rigid body is pushed away as the bubble grows and moved toward the bubble as the bubble collapses. The motion of two bubbles near a movable cylinder is also simulated. The large rotation of the cylinder and obvious deformation and distortion for the bubble in close proximity to a curved wall are observed in our codes. Finally, the growth and collapse of bubble near a deformable ellipsoid shell are also simulated using the combination of boundary element method (BEM) and finite element method (FEM) techniques. The oscillations of the ellipsoid shell can be observed during the growth and collapse of bubble, which much differs from the results obtained by only considering effects of a rigidly movable body on the bubble evolution.     

Simulating 3D deformable particle suspensions using lattice Boltzmann method with discrete external boundary force

A method for direct numerical analysis of three‐dimensional deformable particles suspended in fluid is presented. The flow is computed on a fixed regular ‘lattice’ using the lattice Boltzmann method (LBM), where each solid particle is mapped onto a Lagrangian frame moving continuously through the domain. Instead of the bounce‐back method, an external boundary force (EBF) is used to impose the no‐slip boundary condition at the fluid–solid interface for stationary or moving boundaries. The EBF is added directly to the lattice Boltzmann equation. The motion and orientation of the particles are obtained from Newtonian dynamics equations. The advantage of this approach is outlined in comparison with the standard and higher‐order interpolated bounce‐back methods as well as the LBM immersed‐boundary and the volume‐of‐fluid methods. Although the EBF method is general, in this application, it is used in conjunction with the lattice–spring model for deformable particles. The methodology is validated by comparing with experimental and theoretical results. Copyright © 2009 John Wiley & Sons, Ltd.     

Interaction of a deformable solid with two-phase flows: An Eulerian-based numerical model for fluid-structure interaction using the level contour reconstruction method

We describe the formulation of a method for fluid-structure interaction involving the coupling of moving and/or flexible solid structures with multiphase flows in the framework of the Level Contour Reconstruction Method. We present an Eulerian-based numerical procedure for tracking the motion and interaction of a liquid-gas interface with a fluid-solid interface in the Lagrangian frame together with the evaluation of the fluid transport equations coupled to those for the solid transport, namely the left Cauchy-Green strain tensor field, in the Eulerian frame. To prevent excessive dissipation due to the convective nature of the solid transport equation, a simple incompressibility constraint for the strain field is enforced. A single grid structure is used for both the fluid and solid phases which allows for a simple and natural coupling of the fluid and solid dynamics. Several benchmark tests are performed to show the accuracy of the numerical method and which demonstrate accurate results compared to several of those in the existing literature. In particular we show that surface tension effects including contact line dynamics on the deforming solid phase can be properly simulated. The three-phase interaction of a droplet impacting on a flexible cantilever is investigated in detail. The simulations follow the detailed motion of the droplet impact (and subsequent deformation, breakup, and fall trajectory) along with the motion of the deformable solid cantilever due to its own weight as well as due to the force of the droplet impact.     

Modelling Pressurized Water Reactor cores in terms of porous media

The aim of this study is to develop a tractable model of a nuclear reactor core taking the complexity of the structure (including its nonlinear behaviour) and fluid flow coupling into account. The mechanical behaviour modelling includes the dynamics of both the fuel assemblies and the fluid. Each rod bundle is modelled in the form of a deformable porous medium; then, the velocity field of the fluid and the displacement field of the structure are defined over the whole domain. The fluid and the structure are first modelled separately, before being linked together. The equations of motion for the structure are obtained using a Lagrangian approach and, to be able to link up the fluid and the structure, the equations of motion for the fluid are obtained using an arbitrary Lagragian Eulerian approach. The finite element method is applied to spatially discretize the equations. Simulations are performed to analyse the effects of the characteristics of the fluid and of the structure. Finally, the model is validated with a test involving two fuel assemblies, showing good agreement with the experimental data.     

PIV measurements combined with the motion tracking technique to analyze flow around a moving porous structure

To gain a better understanding of the fluid–structure interaction and especially when dealing with a flow around an arbitrarily moving body, it is essential to develop measurement tools enabling the instantaneous detection of moving deformable interface during the flow measurements. A particularly useful application is the determination of unsteady turbulent flow velocity field around a moving porous fishing net structure which is of great interest for selectivity and also for the numerical code validation which needs a realistic database. To do this, a representative piece of fishing net structure is used to investigate both the Turbulent Boundary Layer (TBL) developing over the horizontal porous moving fishing net structure and the turbulent flow passing through the moving porous structure. For such an investigation, Time Resolved PIV measurements are carried out and combined with a motion tracking technique allowing the measurement of the instantaneous motion of the deformable fishing net during PIV measurements. Once the two-dimensional motion of the porous structure is accessed, PIV velocity measurements are analyzed in connection with the detected motion. Finally, the TBL is characterized and the effect of the structure motion on the volumetric flow rate passing though the moving porous structure is clearly demonstrated.     

STUDY ON THE DYNAMIC BEHAVIOR OF AXIALLY MOVING RECTANGULAR PLATES PARTIALLY SUBMERSED IN FLUID

The natural frequencies, complex modes and critical speeds of an axially moving rectangular plate, which is partially immersed in a fluid and subjected to a pretension, are investigated. The effects of free surface waves, compressibility and viscidity of the fluid are neglected in the analysis. The subsection functions are used to describe the discontinuous characteristics of the system due to partial immersion. The classical thin plate theory is adopted to formulate the equations of motion of a vibrating plate. The velocity potential and Bernoulli's equation are used to describe the fluid pressure acting on the moving plate. The effect of fluid on the vibrations of the plate may be equivalent to the added mass on the plate. The effects of distance ratio, moving speed, immersed-depth ratio, boundary conditions, stiffness ratio and aspect ratio of the plate as well as the fluid-plate density ratios on the free vibrations of the moving plate-fluid system are investigated.     

斑马鱼C型起动中动力学特性的活体实验研究

提出了一种从鱼类自主游动的运动学实验数据出发推算其动力学特性的实验研究方法。该方法基于变形体动力学方程,不仅可以计算出变形鱼体的整体转动角速度以完善其运动学数据,还可以计算出作用在自主游动的鱼体上的流体合力和流体合力矩,进而分析鱼体的力能学特征。本文运用此方法研究了斑马鱼的两种典型C型起动,对比分析了它们的运动学及力能学特征。结果表明,在相仿的C型弯曲变形下,逃逸型C型起动与非逃逸型C型起动相比,前者推力较大,导致其在前进方向上的质心速率较大;前者侧向力较大,导致其转弯半径较小;前者转矩较大,但因在起动中后期的S型摆动产生了反向转矩,最终导致转动角度小于后者。通过对两种典型C型起动的实验研究发现,斑马鱼会因不同的运动需求而表现出不同的机动性能。     

On the stress state of shells penetrating into a deformable solid

The motion of an axisymmetric shell in a deformable solid medium is considered. It is assumed that the medium resistance is described by a two-term expression containing a constant term (the rigidity characteristic) and an inertial term quadratic with respect to the penetration velocity. A model of the impactor penetration with the normal interactions with the resisting medium taken into account is proposed. The membrane forces and the arising stresses are determined for decelerated motions of the impactor.     

Vortex-induced vibrations of pipes conveying fluid in the subcritical and supercritical regimes

In this paper, the vortex-induced vibrations of a hinged–hinged pipe conveying fluid are examined, by considering the internal fluid velocities ranging from the subcritical to the supercritical regions. The nonlinear coupled equations of motion are discretized by employing a four-mode Galerkin method. Based on numerical simulations, diagrams of the displacement amplitude versus the external fluid reduced velocity are constructed for pipes transporting subcritical and supercritical fluid flows. It is shown that when the internal fluid velocity is in the subcritical region, the pipe is always vibrating periodically around the pre-buckling configuration and that with increasing external fluid reduced velocity the peak amplitude of the pipe increases first and then decreases, with jumping phenomenon between the upper and lower response branches. When the internal fluid velocity is in the supercritical region, however, the pipe displays various dynamical behaviors around the post-buckling configuration such as inverse period-doubling bifurcations, periodic and chaotic motions. Moreover, the bifurcation diagrams for vibration amplitude of the pipe with varying internal fluid velocities are constructed for each of the lowest four modes of the pipe in the lock-in conditions. The results show that there is a significant difference between the vibrations of the pipe around the pre-buckling configuration and those around the post-buckling configuration.     

Nonlinear Forced Vibration of Cantilevered Pipes Conveying Fluid

The nonlinear forced vibrations of a cantilevered pipe conveying fluid under base excitations are explored by means of the full nonlinear equation of motion, and the fourthorder Runge-Kutta integration algorithm is used as a numerical tool to solve the discretized equations. The self-excited vibration is briefly discussed first, focusing on the effect of flow velocity on the stability and post-flutter dynamical behavior of the pipe system with parameters close to those in previous experiments. Then, the nonlinear forced vibrations are examined using several concrete examples by means of frequency response diagrams and phase-plane plots. It shows that, at low flow velocity, the resonant amplitude near the first-mode natural frequency is larger than its counterpart near the second-mode natural frequency. The second-mode frequency response curve clearly displays a softening-type behavior with hysteresis phenomenon, while the first-mode frequency response curve almost maintains its neutrality. At moderate flow velocity,interestingly, the first-mode resonance response diminishes and the hysteresis phenomenon of the second-mode response disappears. At high flow velocity beyond the flutter threshold, the frequency response curve would exhibit a quenching-like behavior. When the excitation frequency is increased through the quenching point, the response of the pipe may shift from quasiperiodic to periodic. The results obtained in the present, work highlight the dramatic influence of internal fluid flow on the nonlinear forced vibrations of slender pipes.     

雨刮器非线性摩擦振动的稳定性与分岔特性

针对雨刮器建立2自由度非线性摩擦振动动力学模型,基于复模态理论计算复特征值并进行稳定性及其对刮刷速度的依赖性分析;通过数值计算分析摩擦振动对刮刷速度的分岔特性,并利用相轨迹、庞加莱映射、频谱特性分析不同刮刷速度下的非线性振动现象.研究发现:摩擦-速度特性的负斜率是导致系统不稳定的根本原因,增大刮刷速度有利于提高系统的稳定性;在高、低刮速区,随着刮刷速度的下降,系统振动形态遵循周期→准周期→混沌的演化规律,并会伴随显著的粘滑振动;仅高速区的周期振动和非振动条件下,刮刷时无附加的粘滑振动.     

Parametric vibrations of flexible hoses excited by a pulsating fluid flow,Part I: Modelling,solution method and simulation

The article presents an analysis of a model describing lateral vibrations of a pipe induced by fluid flow velocity pulsation. The motion has been described with a set of two non-linear partial differential equations with periodically variable coefficients. In the analysis Galerkin method has been applied using orthogonal polynomials as shape function. To determine instability regions Floquet theory has been employed. The effect of selected parameters on parametric resonance ranges and regions of increased vibration level has been investigated. The character and form of vibrations have been investigated indicating the possibility of excitation of sub-harmonic and quasi-periodic vibrations in the combination resonance ranges.     

紊动流场中悬浮颗粒分布的随机理论

通过分析固体颗粒在紊动流场中的随机运动,建立了二维流场中垂直于时均流动的方向上颗粒随机位移的概率密度分布函数所满足的方程。由该方程解出的分布函数在一定条件下即相当于颗粒浓度分布函数。运用这一方法研究了[1]、[2]中报道的壁面附近颗粒浓度降低的现象。     

Wave propagation in 3-D poroelastic media including gradient effects

In the framework of the theory of mixtures, the governing equations of motion of a fluid-saturated poroelastic medium including microstructural (for both the solid and the fluid) and micro-inertia (for the solid) effects are derived. This is accomplished by appropriately combining the conservation of mass and linear momentum equations with the constitutive equations for both the solid and the fluid constituents. The solid is assumed to be gradient elastic, that is, its stress tensor depends on the strain and the second gradient of strain tensor. The fluid is assumed to have an analogous behavior, that is, its stress tensor depends on the pressure and the second gradient of pressure. A micro-inertia term in the form of the second gradient of the acceleration of the solid is also included in the equations of motion. The equations of motion in three dimensions are seven equations with seven unknowns, the six displacement components for the solid and the fluid and the pore-fluid pressure. Because of the microstructural effects, the order of these equations is two degrees higher than in the classical case. Application of the divergence and the rot operations on these equations enable one to study the propagation of plane harmonic waves in the infinitely extended medium separately in the form of dilatational and rotational dispersive waves. The effects of the microstructure and the micro-inertia on the dispersion curves are determined and discussed.     

Basic properties of natural vibrations of an extended segment of a pipeline

Natural transverse vibrations of an extended segment of a pipeline containing a uniformly moving fluid are considered. The mechanical model under study takes into account the inertial forces of the pipe and environment and the moment of Coriolis and centrifugal forces arising because of the medium motion. It is proved that all natural frequencies of the pipeline rigidly clamped at both ends are real (and hence no flutter can arise in this model). For the first three modes, the dependence of the eigenvalues on the fluid flow velocity (varying from zero to the buckling velocity) are constructed, and their properties depending on the inertia parameter are studied. Families of vibration mode shapes of the pipeline are constructed and investigated.     

Simulation of forced deformable bodies interacting with two-dimensional incompressible flows: Application to fish-like swimming

We present an efficient algorithm for simulation of deformable bodies interacting with two-dimensional incompressible fluid flows. The temporal and spatial discretizations of the Navier–Stokes equations in vorticity stream-function formulation are based on classical fourth-order Runge–Kutta scheme and compact finite differences, respectively. Using a uniform Cartesian grid we benefit from the advantage of a new fourth-order direct solver for the Poisson equation to ensure the incompressibility constraint down to machine zero over an optimal grid. For introducing a deformable body in fluid flow, the volume penalization method is used. A Lagrangian structured grid with prescribed motion covers the deformable body which is interacting with the surrounding fluid due to the hydrodynamic forces and the torque calculated on the Eulerian reference grid. An efficient law for controlling the curvature of an anguilliform fish, swimming toward a prescribed goal, is proposed which is based on the geometrically exact theory of nonlinear beams and quaternions. Validation of the developed method shows the efficiency and expected accuracy of the algorithm for fish-like swimming and also for a variety of fluid/solid interaction problems.     

Plane Couette flow of viscoplastic materials along a slippery vibrating wall

This study looks at the influence of slip at the wall on plane Couette flows of viscous and yield stress fluids with ultrasonic wall motion. These fluids are used in coating processes. A constant speed V at one wall creates the flow, and vibrations and slip take place at the other wall. Isothermal conditions and arbitrary (longitudinal or transverse) vibrations are considered, with negligible vibrational inertia.For the Bingham model, due to its nonlinearity, whatever the vibration direction and the wall slipperiness, significant decreases occur in the average stress as soon as moderate values of the dimensionless vibration velocity amplitude are involved. Such effects are associated with adherent or slippery walls, even with linear friction laws. They do not occur with linear viscous (Newtonian) models.Average stress reductions can reach nearly 100% for very high Oldroyd numbers, i.e. for stress values without vibration close to the yield limit. Slip velocity also decreases. The cost in terms of the power dissipated remains relatively less than in the Newtonian case, and may contribute to a change in the temperature field. Even when the flow without vibration is a pure slip one, large enough amplitude vibrations, either longitudinal or transverse, applied at the wall can reduce the average shear stress and slip velocity, giving rise to an average axial shear flow.Hence vibrations of moderate or high-velocity amplitude applied to adherent or slippery walls enhance plane Couette flow rates for viscoplastic materials. With moderate values of this amplitude, longitudinal vibrations may be 1.5–2 times more efficient than transverse vibrations with an equivalent cost. However, if for technological reasons transverse vibrations have to be preferred, they can also produce significant results. In any case, coating flows should benefit from an adequate application of ultrasound at the wall.