Flow of Particulate-Fluid Suspension in a Channel with Porous Walls

The problem of the dispersed particulate-fluid two-phase flow in a channel with permeable walls under the effect of the Beavers and Joseph slip boundary condition is concerned in this paper. The analytical solution has been derived for the longitude pressure difference, stream functions, and the velocity distribution with the perturbation method based on a small width to length ratio of the channel. The graphical results for pressure, velocity, and stream function are presented and the effects of geometrical coefficients, the slip parameter and the volume fraction density on the pressure variation, the streamline structure and the velocity distribution are evaluated numerically and discussed. It is shown that the sinusoidal channel, accompanied by a higher friction factor, has higher pressure drop than that of the parallel-plate channel under fully developed flow conditions due to the wall-induced curvature effect. The increment of the channel’s width to the length ratio will remarkably increase the flow rate because of the enlargement of the flow area in the channel. At low Reynolds number ranging from 0 to 65, the fluids move forward smoothly following the shape of the channel. Moreover, the slip boundary condition will notably increase the fluid velocity and the decrease of the slip parameter leads to the increment of the velocity magnitude across the channel. The fluid-phase axial velocity decreases with the increment of the volume fraction density.     

MHD peristaltic flow of a third order fluid in an asymmetric channel

This study is concerned with peristaltic flow of a magnetohydrodynamic (MHD) fluid in an asymmetric channel. Asymmetry in the flow is induced by waves on the channel walls having different amplitudes and phase. A systematic approach based on an expansion of Deborah number is used for the solution series. Analytic expressions have been developed for the stream function, axial velocity and axial pressure gradient. The pressure rise over a wavelength has been addressed through numerical integration. Particular attention has been given to the effects of Hartman number and Deborah number on the pressure rise over a wavelength and the trapping phenomenon. Several limiting solutions of interest are obtained as the special cases of the presented analysis by taking the appropriate parameter(s) to be zero. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical modeling of turbulent compound channel flow using the lattice Boltzmann method

The flow of water in a straight compound channel with prismatic cross section is investigated with a relatively new tool, the lattice Boltzmann method. The large eddy simulation model is added in the lattice Boltzmann model for nonlinear shallow water equations (LABSWE^TM) so that the turbulence, caused by lateral exchange of momentum in the shear layer between the main channel and floodplain, can be taken into account and modeled efficiently. To validate the numerical model, a symmetrical compound channel with trapezoidal main channel and flat floodplain is tested. Similar to most natural watercourses, the floodplain has higher roughness values than the main channel. Different relative depths, D_r (the ratio of the depth of flow on the floodplain to that in the main channel), are considered. The Reynolds number is set at 30 000 in the main channel. The lateral distributions of the longitudinal velocity, the boundary shear stress, the Reynolds stress and the apparent shear stress across the channel are obtained after the large eddy simulation is performed. The results of numerical simulations are compared with the available experiment data, which show that the LABSWE^TM is capable of modeling the features of flow turbulence in compound channels and is sufficiently accurate for practical applications in engineering. Copyright © 2008 John Wiley & Sons, Ltd.     

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

Dispersion of particles in wall-bounded particle-laden turbulent flows with high wall permeability

A series of numerical simulations were performed to investigate the distribution and deposition properties of particles in turbulent flows bounded by permeable walls using the Large Eddy Simulation (LES) with a Lagrangian trajectory approach. The wall permeation speeds were taken from 10⁻⁴ to 10⁻² of the bulk velocity. The directions of the permeation speed were the same at both walls, and they were inward on one wall but outward on the other wall to reserve the fluid mass. Particles with Stokes number (respecting viscous time scale) around 0.1, 1 and 10 were released in the fully developed turbulent channel flow. The particle–particle interaction and the retroaction from particles to the fluid were neglected. The fluid-phase turbulence statistical properties and particle's transport characteristics by vortexes were then analyzed in details. If the wall permeation exists, the turbulence intensities will be depressed close to the outward permeable wall but increased near the inward permeable wall. Not influenced by the wall permeation, the suspended particles with St⁺ ∼O(1) tend to accumulate in the less vortical zones away from the wall, while those particles in the flow regions near the outward permeable wall will distribute disregarding of the vorticity. The turbulence structures near the outward permeable wall are found to exert promotional effects on the particle deposition rate, but such effects are different for particles with various Stokes number. A distribution tendency of streamwise streaks for the deposited particles is also found on the wall imposed by the high outward permeation speed and the clustering deposition pattern is more obvious with increasing particle Stokes number.     

强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.      

Investigation of the development of laminar flow at the entrance of a two-dimensional channel,using a laser-doppler velocimeter

The results of an experimental investigation of the velocity structure of Poiseuille flow at the entrance section of a two-dimensional horizontal channel with a discontinuous velocity profile at the entrance, at Re=300 and 500, are given. The experimental facility and the technique for measuring the velocity field in this channel using a laser-Doppler velocimeter (LDV) are discussed. Velocity profiles are obtained which are concave on the channel axis, and this agrees with the results of the numerical computations of [1].     

基于格子Boltzmann两相模型的粘弹流体挤出胀大分析

挤出胀大的数值模拟是非牛顿流体研究中具有挑战性的问题．本文运用格子Boltzmann方法(LBM)分析Oldroyd-B和多阶松弛谱PTT粘弹流体的挤出胀大现象，采用颜色模型模拟出口处粘弹流体和空气的两相流动，通过重新标色获得两种流体的界面，并最终获得胀大的形状．Navier-Stokes方程和本构方程的求解采用双分布函数模型．将胀大的结果与解析解、实验解和单相自由面LBM结果进行了比较，发现格子Boltzmann两相模型结果与解析解和实验结果相吻合，相比于单相模型，收敛速度更快，解的稳定性更高．研究了流道尺寸对胀大率的影响，并对挤出胀大的内在机理进行了分析．     

A practical implementation of high‐order RKDG models for the 1D open‐channel flow equations

This paper comprises an implementation of a fourth‐order Runge–Kutta discontinuous Galerkin (RKDG4) scheme for computing the open‐channel flow equations. The main features of the proposed methodology are simplicity and easiness in the implementation, which may be of possible interest to water resources numerical modellers. A version of the RKDG4 is blended with the Roe Riemann solver, an adaptive high‐order slope limiting procedure, and high‐order source terms approximations. A comparison of the performance of the proposed method with lower‐order RKDG models is performed showing a benefit of considering the RKDG4 model. The scheme is applied to computerize the Saint Venant system by considering benchmark tests that have exact solutions. Finally, numerical results are illustrated discussing the performance of the proposed high‐order model. Copyright © 2008 John Wiley & Sons, Ltd.     

A similarity solution for the flow and heat transfer over a moving permeable flat plate in a parallel free stream

This paper presents both a numerical and analytical study in connection with the steady boundary layer flow and heat transfer induced by a moving permeable semi-infinite flat plate in a parallel free stream. Both the velocities of the flat plate and the free stream are proportional to x ^1/3. The surface temperature is assumed to be constant. The governing partial differential equations are converted into ordinary differential equations by a new similarity transformation. Numerical results for the flow and heat transfer characteristics are obtained for various values of the moving parameter, transpiration parameter and the Prandtl number. Approximate analytical solutions are also obtained when the suction or injection parameter is very large. It is found that dual solutions exist for the case when the fluid and the plate move in the opposite directions.     

Numerical simulation of vortical ideal fluid flow through curved channel

A numerical algorithm to study the boundary‐value problem in which the governing equations are the steady Euler equations and the vorticity is given on the inflow parts of the domain boundary is developed. The Euler equations are implemented in terms of the stream function and vorticity. An irregular physical domain is transformed into a rectangle in the computational domain and the Euler equations are rewritten with respect to a curvilinear co‐ordinate system. The convergence of the finite‐difference equations to the exact solution is shown experimentally for the test problems by comparing the computational results with the exact solutions on the sequence of grids. To find the pressure from the known vorticity and stream function, the Euler equations are utilized in the Gromeka–Lamb form. The numerical algorithm is illustrated with several examples of steady flow through a two‐dimensional channel with curved walls. The analysis of calculations shows strong dependence of the pressure field on the vorticity given at the inflow parts of the boundary. Plots of the flow structure and isobars, for different geometries of channel and for different values of vorticity on entrance, are also presented. Copyright © 2003 John Wiley & Sons, Ltd.     

环形激波绕射、反射和聚焦流场的间断有限元模拟

采用间断有限元方法对环形激波在圆柱形激波管内绕射、反射和聚焦流场进行了数值模拟。将二维守恒方程的间断有限元方法发展到轴对称Euler方程,并对环形激波绕后台阶流动进行了数值计算。计算结果表明,采用间断有限元方法能够有效地捕捉运动激波在圆柱形激波管内传播的复杂流场结构;在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力,在聚焦点附近不会产生振荡或抹平间断现象。     

One‐dimensional simulation of supercritical flow at a confluence by means of a nonlinear junction model applied with the RKDG2 method

We investigate the one‐dimensional computation of supercritical open‐channel flows at a combining junction. In such situations, the network system is composed of channel segments arranged in a branching configuration, with individual channel segments connected at a junction. Therefore, two important issues have to be addressed: (a) the numerical solution in branches, and (b) the internal boundary conditions treatment at the junction. Going from the advantageous literature supports of RKDG methods to a particular investigation for a supercritical benchmark, the second‐order Runge–Kutta discontinuous Galerkin (RKDG2) scheme is selected to compute the water flow in branches. For the internal boundary handling, we propose a new approach by incorporating the nonlinear model derived from the conservation of the momentum through the junction. The nonlinear junction model was evaluated against available experiments and then applied to compute the junction internal boundary treatment for steady and unsteady flow applications. Finally, a combining flow problem is defined and simulated by the proposed framework and results are illustrated for many choices of junction angles. Copyright © 2007 John Wiley & Sons, Ltd.     

Multiple-relaxation-time lattice Boltzmann method for study of two-lid-driven cavity flow solution multiplicity

As a fundamental subject in fluid mechanics, sophisticated cavity flow patterns due to the movement of multi-lids have been routinely analyzed by the computational fluid dynamics community. Unlike those reported computational studies that were conducted using more conventional numerical methods, this paper features employing the multiple-relaxation-time (MRT) lattice Boltzmann method (LBM) to numerically investigate the two-dimensional cavity flows generated by the movements of two adjacent lids. The obtained MRT-LBM results reveal a number of important bifurcation flow features, such as the symmetry and steadiness of cavity flows at low Reynolds numbers, the multiplicity of stable cavity flow patterns when the Reynolds number exceeds its first critical value, as well as the periodicity of the cavity flow after the second critical Reynolds number is reached. Detailed flow characteristics are reported that include the critical Reynolds numbers, the locations of the vortex centers, and the values of stream function at the vortex centers. Through systematic comparison against the simulation results obtained elsewhere by using the lattice Bhatnagar–Gross–Krook model and other numerical schemes, not only does the MRT-LBM approach exhibit fairly satisfactory accuracy, but also demonstrates its remarkable flexibility that renders the adjustment of its multiple relaxation factors fully manageable and, thus, particularly accommodates the need of effectively investigating the multiplicity of flow patterns with complex behaviors.     

Peristaltic motion of third grade fluid in curved channel

＂ Analysis is performed to study the slip effects on the peristaltic flow of non-Newtonian fluid in a curved channel with wall properties. The resulting nonlinear partial differential equations are transformed to a single ordinary differential equation in a stream function by using the assumptions of long wavelength and low Reynolds number. This differential equation is solved numerically by employing the built-in routine for solving nonlinear boundary value problems （BVPs） through the software Mathematica. In addition, the analytic solutions for small Deborah number are computed with a regular perturbation technique. It is noticed that the symmetry of bolus is destroyed in a curved channel. An intensification in the slip effect results in a larger magnitude of axial velocity. Further, the size and circulation of the trapped boluses increase with an increase in the slip parameter. Different from the case of planar channel, the axial velocity profiles are tilted towards the lower part of the channel. A comparative study between analytic and numerical solutions shows excellent agreement.     

Flow past a plate lattice with developed cavitation

Problems of streamline cavitation flow past a lattice were studied in [1–8] using the Kirchhoff scheme. In this scheme the magnitude of the velocity at the free surface is equal to the stream velocity behind the lattice, and the cavitation number is zero (for a lattice the relative velocity and the cavitation number are defined from the stream velocity behind the lattice). In [4, 7] a solution is given of the problem of flow past a lattice using a scheme with an Efros-Gilbargreturn streamline, which permits considering arbitrary cavitation numbers; however, a unique solution is not given. Some other streamline schemes are mentioned in [8].In the following we consider the cavitational flow of an ideal incompressible inviscid and weightless fluid past an infinite lattice of flat plates, using the streamline wake model previously utilized by Wu [9] in studying cavitational flow past an isolated obstacle. In accordance with this model, the streamlines which separate from the body and bound the cavity behind it pass into two curvilinear infinitely long walls, along which the pressure increases and approaches the pressure in the undisturbed stream.It is further assumed that in the hodograph plane there corresponds to the curvilinear walls a cut along some line and that the complex potential takes the same values at points lying on opposite sides of the cut. In particular, at the points of contact of the streamlines with the curvilinear walls the complex potential is the same. In the Wu scheme the latter condition leads to vanishing of the velocity circulation along the contour CABC₁ (Fig. 1).In conclusion the author wishes to thank N. V. Yurtaeva for the accurately performed numerical work.     

Unsteady turbulent flow of a gas suspension in a channel under conditions of injection and forced pressure oscillations

A turbulent flow of a suspension of solid particles in a gas is considered. The suspension is located in a channel with permeable walls (the pressure at the left end face of the channel follows a sinusoidal law). The flow considered here reflects the principal features of the flow in the combustion chamber of a solid-propellant rocket motor. The unsteady flow of the gas suspension is described by using the Eulerian-Lagrangian approach. A stochastic variant of the discrete-trajectory approach is used for modeling the particle motion. The influence of the condensed phase on the turbulence characteristics and acoustic oscillations of the parameters of the working medium in the channel in the case of injection is discussed. The calculated results are compared with data obtained in a physical experiment.     

Finite element analysis of air flow around permeable sand fences

A numerical study of the turbulent air flow in a trench trap and the turbulent flow around a permeable sand fence is reported in this paper. The two-dimensional modified k–ε turbulence model proposed by Kato and Launder is used to predict the turbulent characteristics of the air flow. The discretization method for the governing equations is the three-step Taylor/Galerkin finite element method proposed by the authors. For the flow in a trench trap the numerical results are compared with experimental data obtained under realistic conditions using a large wind tunnel. For the air flow around a permeable sand fence a pressure loss model is used to represent the effect of the porosity of the fence on the flow field. © 1997 John Wiley & Sons, Ltd.