Exact solutions for the flow of second grade fluid in annulus between torsionally oscillating cylinders

The velocity field and the associated shear stress corresponding to the torsional oscillatory flow of a second grade fluid, between two infinite coaxial circular cylinders, are determined by means of the Laplace and Hankel transforms. At time t = 0, the fluid and both the cylinders are at rest and at t = 0 + , cylinders suddenly begin to oscillate around their common axis in a simple harmonic way having angular frequencies ω 1 and ω 2 . The obtained solutions satisfy the governing differential equation and all imposed initial and boundary conditions. The solutions for the motion between the cylinders, when one of them is at rest, can be obtained from our general solutions. Furthermore, the corresponding solutions for Newtonian fluid are also obtained as limiting cases of our general solutions.     

Effect of induced magnetic field on natural convection in vertical concentric annuli

In the present paper,we have considered the steady fully developed laminar natural convective flow in open ended vertical concentric annuli in the presence of a radial magnetic field.The induced magnetic field produced by the motion of an electrically conducting fluid is taken into account.The transport equations concerned with the considered model are first recast in the non-dimensional form and then unified analytical solutions for the velocity,induced magnetic field and temperature field are obtained for the cases of isothermal and constant heat flux on the inner cylinder of concentric annuli.The effects of the various physical parameters appearing into the model are demonstrated through graphs and tables.It is found that the magnitude of maximum value of the fluid velocity as well as induced magnetic field is greater in the case of isothermal condition compared with the constant heat flux case when the gap between the cylinders is less or equal to 1.70 times the radius of inner cylinder,while reverse trend occurs when the gap between the cylinders is greater than 1.71 times the radius of inner cylinder.These fields are almost the same when the gap between the cylinders is equal to 1.71 times the radius of inner cylinder for both the cases.It is also found that as the Hartmann number increases,there is a flattening tendency for both the velocity and the induced magnetic field.The influence of the induced magnetic field is to increase the velocity profiles.     

Unsteady flow of viscoelastic fluid between two cylinders using fractional Maxwell model

The unsteady flow of an incompressible fractional Maxwell fluid between two infinite coaxial cylinders is studied by means of integral transforms.The motion of the fluid is due to the inner cylinder that applies a time dependent torsional shear to the fluid.The exact solutions for velocity and shear stress are presented in series form in terms of some generalized functions.They can easily be particularized to give similar solutions for Maxwell and Newtonian fluids.Finally,the influence of pertinent parameters on the fluid motion,as well as a comparison between models,is highlighted by graphical illustrations.     

Shear flow induced vibrations of long slender cylinders with a wake oscillator model

A time domain model is presented to study the vibrations of long slender cylinders placed in shear flow. Long slender cylinders such as risers and tension legs are widely used in the field of ocean engineering. They are subjected to vortex-induced vibrations(VIV) when placed within a transverse incident flow. A three dimensional model coupled with wake oscillators is formulated to describe the response of the slender cylinder in cross-flow and in-line directions. The wake oscillators are distributed along the cylinder and the vortex-shedding frequency is derived from the local current velocity. A non-linear fiuid force model is accounted for the coupled effect between cross-flow and in-line vibrations. The comparisons with the published experimental data show that the dynamic features of VIV of long slender cylinder placed in shear flow can be obtained by the proposed model,such as the spanwise average displacement,vibration frequency,dominant mode and the combination of standing and traveling waves. The simulation in a uniform flow is also conducted and the result is compared with the case of nonuniform flow. It is concluded that the flow shear characteristic has significantly changed the cylinder vibration behavior.     

Hybrid algorithm for modeling of fluid-structure interaction in incompressible, viscous flows

The objective of this paper is to present and to validate a new hybrid coupling (HC) algorithm for modeling of fluid-structure interaction (FSI) in incompressible, viscous flows. The HC algorithm is able to avoid numerical instability issues associated with artificial added mass effects, which are often encountered by standard loosely coupled (LC) and tightly coupled (TC) algorithms, when modeling the FSI response of flexible structures in incompressible flow. The artificial added mass effect is caused by the lag in exchange of interfacial displacements and forces between the fluid and solid solvers in partitioned algorithms. The artificial added mass effect is much more prominent for light/flexible structures moving in water, because the fluid forces are in the same order of magnitude as the solid forces, and because the speed at which numerical errors propagate in an incompressible fluid. The new HC algorithm avoids numerical instability issues associated with artificial added mass effects by embedding Theodorsen’s analytical approximation of the hydroelastic forces in the solution process to obtain better initial estimates of the displacements. Details of the new HC algorithm are presented. Numerical validation studies are shown for the forced pitching response of a steel and a plastic hydrofoil. The results show that the HC algorithm is able to converge faster, and is able to avoid numerical instability issues, compared to standard LC and TC algorithms, when modeling the transient FSI response of a plastic hydrofoil. Although the HC algorithm is only demonstrated for a NACA0009 hydrofoil subject to pure pitching motion, the method can be easily extended to model general 3-D FSI response and stability of complex, flexible structures in turbulent, incompressible, multiphase flows.     

The SPH approach to the process of container filling based on non-linear constitutive models

In this work,the transient free surface of container filling with non-linear constitutive equation’s fluids is numerically investigated by the smoothed particle hydrodynamics(SPH) method.Specifically,the filling process of a square container is considered for non-linear polymer fluids based on the Cross model.The validity of the presented SPH is first verified by solving the Newtonian fluid and OldroydB fluid jet.Various phenomena in the filling process are shown,including the jet buckling,jet thinning,splashing or spluttering,steady filling.Moreover,a new phenomenon of vortex whirling is more evidently observed for the Cross model fluid compared with the Newtonian fluid case.     

Natural convection flow of a couple stress fluid between two vertical parallel plates with Hall and ion-slip effects

The Hall and ion-slip effects on fully developed electrically conducting couple stress fluid flow between vertical parallel plates in the presence of a temperature dependent heat source are investigated. The governing non-linear partial differential equations are transformed into a system of ordinary differential equations using similarity transformations. The resulting equations are then solved using the homotopy analysis method (HAM). The effects of the magnetic parameter, Hall parameter, ion-slip parameter and couple stress fluid parameter on velocity and temperature are discussed and shown graphically.     

Scale-adaptive simulation of flow past wavy cylinders at a subcritical Reynolds number

The Xu & Yan scale-adaptive simulation (XYSAS) model is employed to simulate the flows past wavy cylinders at Reynolds number 8 × 10 3.This approach yields results in good agreement with experimental measurements.The mean flow field and near wake vortex structure are replicated and compared with that of a corresponding circular cylinder.The effects of wavelength ratios λ/D m from 3 to 7,together with the amplitude ratios a /D m of 0.091 and 0.25,are fully investigated.Owing to the wavy configuration,a maximum reduction of Strouhal number and root-meansquare (r.m.s) fluctuating lift coefficients are up to 50% and 92%,respectively,which means the vortex induced vibration (VIV) could be effectively alleviated at certain larger values of λ/D m and a /D m.Also,the drag coefficients can be reduced by 30%.It is found that the flow field presents contrary patterns with the increase of λ/D m.The free shear layer becomes much more stable and rolls up into mature vortex only further downstream when λ/D m falls in the range of 5-7.The amplitude ratio a /D m greatly changes the separation line,and subsequently influences the wake structures.     

正方形排列四圆柱体绕流特性及互扰效应研究

基于计算流体动力学理论,运用大涡模拟方法对雷诺数Re=3900三维正方形排列四圆柱体结构群的绕流问题进行数值计算,主要分析来流攻角与间距比两个参数对四圆柱体结构群流体参数及流场模态的影响。结果表明:来流攻角与间距比均对四圆柱体结构群绕流特性有较强的影响;来流攻角θ=0°、22.5°、45°下,临界间距比分别为3.5、4.0、3.0;间距比的变化会导致下游圆柱表面压力系数分布发生改变;另一方面,间距比较小时,四圆柱体结构之间的互扰作用均以临近效应为主;随间距比增大,上游圆柱尾流对下游圆柱有显著影响,其互扰作用会转变尾激效应。     

Formulation of thermo-hydro-mechanical coupling behavior of unsaturated soils based on hybrid mixture theory

Thermo-Hydro-Mechanical （THM） coupling pro- cesses in unsaturated soils are very important in both theoretical researches and engineering applications. A coupled formulation based on hybrid mixture theory is derived to model the THM coupling behavior of unsaturated soils. The free-energy and dissipative functions for different phases are derived from Taylor＇s series expansions. Constitutive relations for THM coupled behaviors of unsaturated soils, which include deformation, entropy change, fluid flow, heat conduction, and dynamic compatibility conditions on the interfaces, are then established. The number of field equations is shown to be equal to the number of unknown variables; thus, a closure of this coupling problem is established. In addition to modifications of the physical conservation equations with coupling effect terms, the constitutive equations, which consider the coupling between elastoplastic deformation of the soil skeleton, fluid flow, and heat transfer, are also derived.     

涡激诱导并列双圆柱碰撞数值模拟研究

圆柱类结构物的涡激振动是工程中较为常见的一种现象,如果圆柱结构物之间的距离较小, 就会产生涡激诱导碰撞现象,而涡激碰撞会比涡激振动对结构物疲劳破坏产生更严重的威胁.采用浸入边界法模拟流体中的动边界问题,避免了传统贴体网格方法在求解流体中存在固体间碰撞问题时出现数值求解不稳定问题,采用有限元方法对圆柱的运动和碰撞进行求解,通过数据回归方法建立了流体流动条件下的润滑模型,对不同间隙比下涡激诱导并列双圆柱振动及碰撞过程进行了数值模拟, 数值结果表明,如果两圆柱产生了碰撞将会有连续的碰撞发生, 碰撞时出现了多阶频率,振动主频率要比无碰撞时大, 两圆柱碰撞时的相对速度比自由来流速度小;当两圆柱相互接近时, 随着涡环分离角度的逐渐倾斜, 横向流体力先逐渐减小,当两圆柱间涡环开始相互影响发生挤压时, 横向流体力开始逐渐增大;当两圆柱开始反弹时, 两圆柱间形成了低压区, 改变了横向流体阻力的方向,使两圆柱又产生了接近运动,如此反复从而产生了碰撞后横向流体力和圆柱速度的振荡现象.      

Simultaneous CFD evaluation of VIV suppression using smaller control cylinders

Two-degree-of-freedom vortex-induced vibrations (VIV) of a circular cylinder with and without two smaller control cylinders are investigated numerically by computational fluid dynamics (CFD) models coupling with a fluid–structure interaction (FSI) computational method. The numerical model is validated against experimental data of VIV of an isolated cylinder in uniform current. The study is aimed to investigate the effect of smaller control cylinders on VIV suppression. The trajectories of cylinder motion, amplitude response, and temporal evolution of vortex shedding and streamlines are obtained by conducting a series of simulations. And the effect of Reynolds number, located angle and rotational rate of small control cylinders are discussed in detail. It is found that placing small cylinders at 45° to the downstream vector can achieve a good suppression effect, but the effect is different at different Re. Rotating control cylinders with a reasonable rotation velocity can further enhance the VIV suppression by injecting enough momentum into the boundary layer of the main cylinder. The best effect is found at U_c=10, which has a 64.56% reduction in the transverse vibration response.     

Numerical predictions for unsteady viscous flow past an array of cylinders

The unsteady two-dimensional flow around an array of circular cylinders submerged in a uniform onset flow is analysed. The fluid is taken to be viscous and incompressible. The array of cylinders consists of two horizontal rows extending to infinity in the upstream and downstream directions. The centre-to-centre distance between adjacent cylinders is fixed at three diameters, and the rows are staggered. Advantage is taken of spatially periodic boundary conditions in the flow direction. This reduces the computational domain to a rectangular region surrounding a single circular cylinder. Two cases, for Reynolds numbers of 1000 and 10,000, are presented.     

Three-dimensional computation for flow-induced vibrations of an upstream circular cylinder in two tandem circular cylinders

It is well known from a lot of experimental data that fluid forces acting on two tandem circular cylinders are quite different from those acting on a single circular cylinder. Therefore, we first present numerical results for fluid forces acting on two tandem circular cylinders, which are mounted at various spacings in a smooth flow, and second we present numerical results for flow-induced vibrations of the upstream circular cylinder in the tandem arrangement. The two circular cylinders are arranged at close spacing in a flow field. The upstream circular cylinder is elastically placed by damper-spring systems and moves in both the in-line and cross-flow directions. In such models, each circular cylinder is assumed as a rigid body. On the other hand, we do not introduce a turbulent model such as the Large Eddy Simulation (LES) or Reynolds Averaged Navier-Stokes (RANS) models into the numerical scheme to compute the fluid flow. Our numerical procedure to capture the flow-induced vibration phenomena of the upstream circular cylinder is treated as a fluid-structure interaction problem in which the ideas of weak coupling is taken into consideration.     

Effects of thermophoresis particle deposition and of the thermal conductivity in a porous plate with dissipative heat and mass transfer

Network simulation method(NSM) is used to solve the laminar heat and mass transfer of an electricallyconducting,heat generating/absorbing fluid past a perforated horizontal surface in the presence of viscous and Joule heating problem. The governing partial differential equations are non-dimensionalized and transformed into a system of nonlinear ordinary differential similarity equations,in a single independent variable,η. The resulting coupled,nonlinear equations are solved under appropriate transformed boundary conditions. Computations are performed for a wide range of the governing flow parameters,viz Prandtl number,thermophoretic coeffcient(a function of Knudsen number),thermal conductivity parameter,wall transpiration parameter and Schmidt number. The numerical details are discussed with relevant applications. The present problem finds applications in optical fiber fabrication,aerosol filter precipitators,particle deposition on hydronautical blades,semiconductor wafer design,thermo-electronics and problems including nuclear reactor safety.     

THE FLUID AND MECHANICAL COUPLING BETWEEN TWO CIRCULAR CYLINDERS IN TANDEM ARRANGEMENT

The effect of the fluid and mechanical coupling between two flexible circular cylinders on their vortex-induced oscillations is examined in a wind tunnel. The cylinders are in a tandem arrangement and a mechanical coupling allows vibrations in phase and out of phase. The separation between the cylinders ranges between 7 and 25 diameters. When the cylinders are coupled only by fluid, the upstream cylinder behaves as an isolated one, while the downstream cylinder response depends on the separation and presents a hysteresis and three discontinuities. The hysteresis and its discontinuities agree with those of the upstream cylinder and the third discontinuity coincides with the end of synchronization of the windward cylinder. The addition of the mechanical coupling complicates the dynamic response which is greatly influenced by the separation and the coupling mode involved. The results are interpreted according to those dealing with the fluid coupling only. The effect of small-scale grid turbulence is also examined.     

Peristaltic flow of Sisko fluid in a uniform inclined tube

We have analyzed an incompressible Sisko fluid through an axisymmetric uniform tube with a sinusoidal wave propagating down its walls. The present analysis of non- Newtonian fluid is investigated under the considerations of long wavelength and low Reynolds number approximation. The analytic solution is obtained using （i） the regular perturbation method （ii） the Homotopy analysis method （HAM）. The comparison of both the solutions is presented graphically. The results for the pressure rise, frictional force and pressure gradient have been calculated numerically and the results are studied for various values of the physical parameters of interest, such as α （angle of inclination）, b^＊ （Sisko fluid parameter）, Ф （amplitude ratio） and n （power law index）. Trapping phenomena is discussed at the end of the article.     

On the pressure and stress singularities induced by steady flows of a pair of nonmiscible, incompressible, viscous fluids in contact with a wall

Design for structural integrity requires an appreciation of where stress singularities can occur in structural configurations. While there is a rich literature devoted to the identification of such singular behavior in solid mechanics, only of late has there been much in the way of corresponding identifications of flow-induced stress singularities in fluid mechanics. These recent asymptotic identifications are for a single incompressible viscous fluid： Here the asymptotic approach is extended to apply to a configuration entailing two such fluids, For this configuration, various specifications leading to power or log singularities are determined. These results demonstrate that flow-induced stress singularities can occur in a structural container at a location where no singularities are identified within solid mechanics alone.     

A class of exact solutions for the incompressible viscous magnetohydrodynamic flow over a porous rotating disk

The present paper is concerned with a class of exact solutions to the steady Navier-Stokes equations for the incompressible Newtonian viscous electrically conducting fluid flow due to a porous disk rotating with a constant angular speed.The three-dimensional hydromagnetic equations of motion are treated analytically to obtained exact solutions with the inclusion of suction and injection.The well-known thinning/thickening flow field effect of the suction/injection is better understood from the constructed closed form velocity equations.Making use of this solution,analytical formulas for the angular velocity components as well as for the permeable wall shear stresses are derived.Interaction of the resolved flow field with the surrounding temperature is further analyzed via the energy equation.The temperature field is shown to accord with the dissipation and the Joule heating.As a result,exact formulas are obtained for the temperature field which take different forms corresponding to the condition of suction or injection imposed on the wall.     

Numerical modeling and simulation of PEM fuel cells: Progress and perspective

This paper provides a comprehensive review on the research and development in multi-scale numerical modeling and simulation of PEM fuel cells. An overview of recent progress in PEM fuel cell modeling has been provided. Fundamental transport phenomena in PEM fuel cells and the corresponding mathematical formulation of macroscale models are analyzed. Various important issues in PEM fuel cell modeling and simulation are examined in detail, including fluid flow and species transport, electron and proton transport, heat transfer and thermal management, liquid water transport and water management, transient response behaviors, and cold-start processes. Key areas for further improvements have also been discussed.