Numerical simulation of transient planar flow of a viscoelastic material with two moving free surfaces

In this study, we examine the numerical simulation of transient viscoelastic flows with two moving free surfaces. A modified Galerkin finite element method is implemented to the two-dimensional non-steady motion of the fluid of the Oldroyd-B type. The fluid is initially placed between two parallel plates and bounded by two straight free boundaries. In this Lagrangian finite element method, the spatial mesh deforms in time along with the moving free boundaries. The unknown shape of the free surfaces is determined with the flow field u, v, τ, p by the deformable finite element method, combined with a predictor-corrector scheme in an uncoupled fashion. The moving free surfaces and fluid motion of both Newtonian and non-Newtonian flows are investigated. The results include the influence of surface tension, fluid inertia and elasticity.     

Periodic solution of turbulent oscillating channel flows

The time-dependent turbulent Navier–Stokes equations are solved numerically by a finite element method with an algebraic eddy viscosity model (Baldwin–Lomax formulation) for oscillating turbulent channel flows. The method of averaging is used to analyse the resulting periodic motion of the fluid. Numerical results are obtained for various Strouhal numbers and relative amplitudes. A comparison is made between the numerical and published experimental results. It appears that for low relative amplitudes in a certain range of frequencies the agreement is satisfactory.     

Numerical analysis of finite amplitude motion of waves and a moored floating body under severe storm conditions

The motion of a moored floating body under the action of wave forces, which is influenced by fluid forces, shape of the floating body and mooring forces, should be analysed as a complex coupled motion system. Especially under severe storm conditions or resonant motion of the floating body it is necessary to consider finite amplitude motions of the waves, the floating body and the mooring lines as well as non-linear interactions of these finite amplitude motions. The problem of a floating body has been studied on the basis of linear wave theory by many researchers. However, the finite amplitude motion under a correlated motion system has rarely been taken into account. This paper presents a numerical method for calculating the finite amplitude motion when a floating body is moored by non-linear mooring lines such as chains and cables under severe storm conditions.     

流固偶合运动的边界元法

本文给出了流固偶合运动(包括物体散射辐射及偶合运动)的边界元法理论和应用.对于散射问题,求出了物体引起的散射势及入射波作用于物体的载荷.对于辐射问题,求出了辐射势及物体在流体中运动的附加质量和附加阻尼.偶合问题包括求其中包含的散射势和辐射势以及作用于物体之上的散射力、物体的附加质量、附加阻尼、物体在入射波作用下的运动.在偶合运动问题中,本文采取了边界积分方程与物体在流体中的运动方程联立求解的方法,并将其运用到边界元法的数值过程中.所编制的程序有较高的精度.最后给出了数值计算结果与理论解的比较.     

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.     

Forces on oscillating bodies in viscous fluid

This paper describes a method for determining the fluid forces on oscillating bodies in viscous fluid when the corresponding flow problem has been solved using the finite element method. These forces are characterized by the concept of added mass, added damping and added force. Numerical results are obtained for several example body shapes. Comparison is made with exact analytical results and other finite element results for the limiting cases of Stoke's flow and inviscid flow, and good agreement is obtained. The results for finite values of the body amplitude parameter β show the appearance of added force from the steady streaming component of the flow for asymmetric bodies. Results are also obtained for the associated flow where the fluid remote from a fixed body is oscillating.     

Finite element analysis of non-Newtonian fluid flow in 2-d branching channel

This paper presents finite element analysis of non-Newtonian fluid flow in 2-d branching channel. The Galerkin method and mixed finite element method are used. Here the fluid is considered as incompressible, non-Newtonian fluid with Oldyord differential-type constitutive equation. The non-linear algebraic equation system which is formulated with finite element method is solved by means of continuous differential method. The results show that finite element method is suitable for the analysis of non-Newtonian fluid flow with complex geometry.     

圆板非线性振动有限元分析的一种迭代方法

同时考虑横向振动和板平面内的运动,用3节点有限元研究均匀圆板的轴对称大振幅非线性振动,构造了一个避免发散加速收敛的平均迭代法,并将计算结果与文献的已有结果做了比较。     

Internal-external resonance of beams on non-linear viscoelastic foundation traversed by moving load

Vibration of a finite Euler–Bernoulli beam, supported by non-linear viscoelastic foundation traversed by a moving load, is studied and the Galerkin method is used to discretize the non-linear partial differential equation of motion. Subsequently, the solution is obtained for different harmonics using the Multiple Scales Method (MSM) as one of the perturbation techniques. Free vibration of a beam on non-linear foundation is investigated and the effects of damping and non-linear stiffness of the foundation on the responses are examined. Internal-external resonance condition is then stated and the frequency responses of different harmonics are obtained by MSM. Different conditions of the external resonance are studied and a parametric study is carried out for each case. The effects of damping and non-linear stiffness of the foundation as well as the magnitude of the moving load on the frequency responses are investigated. Finally, a thorough local stability analysis is performed on the system.     

A coupled fluid-elasticity model for the wave forcing of an ice-shelf

A mathematical model for predicting the vibrations of ice-shelves based on linear elasticity for the ice-shelf motion and potential flow for the fluid motion is developed. No simplifying assumptions such as the thinness of the ice-shelf or the shallowness of the fluid are made. The ice-shelf is modelled as a two-dimensional elastic body of an arbitrary geometry under plane-strain conditions. The model is solved using a coupled finite element method incorporating an integral equation boundary condition to represent the radiation of energy in the infinite fluid. The solution is validated by comparison with thin-beam theory and by checking energy conservation. Using the analyticity of the resulting linear system, we show that the finite element solution can be extended to the complex plane using interpolation of the linear system. This analytic extension shows that the system response is governed by a series of singularities in the complex plane. The method is illustrated through time-domain simulations as well as results in the frequency domain.     

A one-dimensional moving grid solution for the coupled non-linear equations governing multiphase flow in porous media. 1: Model development

A straightforward moving grid finite element method is developed to solve the one-dimensional coupled system of non-linear partial differential equations (PDEs) governing two- and three-phase flow in porous media. The method combines features from a number of self-adaptive grid techniques. These techniques are the equidistribution, the moving grid finite element and the local grid refinement/coarsening methods. Two equidistribution criteria, based on solution gradient and curvature, are employed and nodal distributions are computed iterativcly. Using the developed approach, an intermingle-free nodal distribution is guaranteed. The method involves examination of a single representative gradient to facilitate the application of moving grid algorithms to solve a non-linear coupled set of PDEs and includes a feature to limit mass balance error during nodal redistribution. The finite element part of the developed algorithm is verified against an existing finite difference model. A numerical simulation example involving a single-front two-phase flow problem is presented to illustrate model performance. Additional simulation examples are given in Part 2 of this paper. These examples include single and double moving fronts in two- and three-phase flow systems incorporating source/sink terms. Simulation sensitivity to the moving grid parameters is also explored in Part 2.     

The analysis of unsteady incompressible flows by a three-step finite element method

This paper describes a three-step finite element method and its applications to unsteady incompressible fluid flows. Stability analysis of the one-dimensional pure convection equation shows that this method has third-order accuracy and an extended numerical stability domain in comparison with the Lax--Wendroff finite element method. The method is cost-effective for incompressible flows because it permits less frequent updates of the pressure field with good accuracy. In contrast with the Taylor-Galerkin method, the present method does not contain any new higher-order derivatives, which makes it suitable for solving non-linear multidimensional problems and flows with complicated boundary conditions. The three-step finite element method has been used to simulate unsteady incompressible flows. The numerical results obtained are in good agreement with those in the literature.     

Parallel finite element simulation of mooring forces on floating objects

The coupling between the equations governing the free‐surface flows, the six degrees of freedom non‐linear rigid body dynamics, the linear elasticity equations for mesh‐moving and the cables has resulted in a fluid‐structure interaction technology capable of simulating mooring forces on floating objects. The finite element solution strategy is based on a combination approach derived from fixed‐mesh and moving‐mesh techniques. Here, the free‐surface flow simulations are based on the Navier–Stokes equations written for two incompressible fluids where the impact of one fluid on the other one is extremely small. An interface function with two distinct values is used to locate the position of the free‐surface. The stabilized finite element formulations are written and integrated in an arbitrary Lagrangian–Eulerian domain. This allows us to handle the motion of the time dependent geometries. Forces and momentums exerted on the floating object by both water and hawsers are calculated and used to update the position of the floating object in time. In the mesh moving scheme, we assume that the computational domain is made of elastic materials. The linear elasticity equations are solved to obtain the displacements for each computational node. The non‐linear rigid body dynamics equations are coupled with the governing equations of fluid flow and are solved simultaneously to update the position of the floating object. The numerical examples includes a 3D simulation of water waves impacting on a moored floating box and a model boat and simulation of floating object under water constrained with a cable. Copyright © 2003 John Wiley & Sons, Ltd.     

Dynamic coupling of fluid and structural mechanics for simulating particle motion and interaction in high speed compressible gas particle laden flow

In this work, structural finite element analyses of particles moving and interacting within high speed compressible flow are directly coupled to computational fluid dynamics and heat transfer analyses to provide more detailed and improved simulations of particle laden flow under these operating conditions. For a given solid material model, stresses and displacements throughout the solid body are determined with the particle–particle contact following an element to element local spring force model and local fluid induced forces directly calculated from the finite volume flow solution. Plasticity and particle deformation common in such a flow regime can be incorporated in a more rigorous manner than typical discrete element models where structural conditions are not directly modeled. Using the developed techniques, simulations of normal collisions between two 1 mm radius particles with initial particle velocities of 50–150 m/s are conducted with different levels of pressure driven gas flow moving normal to the initial particle motion for elastic and elastic–plastic with strain hardening based solid material models. In this manner, the relationships between the collision velocity, the material behavior models, and the fluid flow and the particle motion and deformation can be investigated. The elastic–plastic material behavior results in post collision velocities 16–50% of their pre-collision values while the elastic-based particle collisions nearly regained their initial velocity upon rebound. The elastic–plastic material models produce contact forces less than half of those for elastic collisions, longer contact times, and greater particle deformation. Fluid flow forces affect the particle motion even at high collision speeds regardless of the solid material behavior model. With the elastic models, the collision force varied little with the strength of the gas flow driver. For the elastic–plastic models, the larger particle deformation and the resulting increasingly asymmetric loading lead to growing differences in the collision force magnitudes and directions as the gas flow strength increased. The coupled finite volume flow and finite element structural analyses provide a capability to capture the interdependencies between the interaction of the particles, the particle deformation, the fluid flow and the particle motion.     

Non-linear free surface flows and an application of the Orlanski boundary condition

The focus of this paper is the analysis of spatially two-dimensional non-linear free surface problems. The critical aspects of the problem concern the treatment of the non-linear free surface, the body boundary condition for large motions and the imposition of suitable radiation conditions. To address such complexities, time domain simulation was chosen as the method of analysis. With the use of a finite domain for simulation, a major concern is with the radiation condition to be applied at the open or truncation boundary. For the two-dimensional problem at hand, no theoretical radiation conditions are known to exist. An extension of the Orlanski open boundary condition, based on phase velocity determination at the free surface, is proposed. Three categories of problems were analysed using numerical simulation-namely, freely moving steep waves, waves over a submerged body and forced body motion. Simulation results have been compared with linear theory and experiments.     

Dynamical analysis of axially moving plate by finite difference method

The complex natural frequencies for linear free vibrations and bifurcation and chaos for forced nonlinear vibration of axially moving viscoelastic plate are investigated in this paper. The governing partial differential equation of out-of-plane motion of the plate is derived by Newton’s second law. The finite difference method in spatial field is applied to the differential equation to study the instability due to flutter and divergence. The finite difference method in both spatial and temporal field is used in the analysis of a nonlinear partial differential equation to detect bifurcations and chaos of a nonlinear forced vibration of the system. Numerical results show that, with the increasing axially moving speed, the increasing excitation amplitude, and the decreasing viscosity coefficient, the equilibrium loses its stability and bifurcates into periodic motion, and then the periodic motion becomes chaotic motion by period-doubling bifurcation.     

A fictitious domain finite element method for simulations of fluid–structure interactions: The Navier–Stokes equations coupled with a moving solid

The paper extends a stabilized fictitious domain finite element method initially developed for the Stokes problem to the incompressible Navier–Stokes equations coupled with a moving solid. This method presents the advantage to predict an optimal approximation of the normal stress tensor at the interface. The dynamics of the solid is governed by Newton׳s laws and the interface between the fluid and the structure is materialized by a level-set which cuts the elements of the mesh. An algorithm is proposed in order to treat the time evolution of the geometry and numerical results are presented on a classical benchmark of the motion of a disk falling in a channel.     

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.     

A simple method for simulating viscoelastic fluid flows in slit channel

A low-cost semi-analysis finite element technique, named the finite piece method (FPM) is presented in this article. It aims to solve three-dimensional (3D) viscoelastic slit flows. The viscoelastic stress of the fluid is modelled using an K-BKZ integral constitutive equation of the Wagner type. Picard iteration is used to solve non-linear equations. The FPM is tested on flow problems in both planar and contraction channels. The accuracy of the method is assessed by comparing flow distributions and pressure with results obtained by 3D finite element method (FEM). It shows that the solution accuracy is excellent and a substantial amount of computing time and memory requirement can be saved.     

Stopping a Viscous Fluid by a Feedback Dissipative Field: I. The Stationary Stokes Problem

We consider a planar stationary flow of an incompressible viscous fluid in a semiinfinite strip governed by the Stokes system with a body forces field. We show how this fluid can be stopped at a finite distance of the entrance of the semi-infinite strip by means of a feedback field depending in a sub-linear way on the velocity field. This localization effect is proved reducing the problem to a non-linear bi-harmonic type one for which the localization of solutions is obtained by means of the application of a suitable energy method. Since the presence of the non-linear terms defined through the body forces field is not standard in the fluid mechanics literature, we establish also some results about the existence and uniqueness of weak solutions for this problem.