Flow of viscoelastic fluids through a porous channel—I

The flow of viscoelastic fluids through a porous channel with one impermeable wall is computed. The flow is characterized by a boundary value problem in which the order of the differential equation exceeds the number of boundary conditions. Three solutions are developed: (i) an exact numerical solution, (ii) a perturbation solution for small R, the cross-flow Reynold's number and (iii) an asymptotic solution for large R. The results from exact numerical integration reveal that the solutions for a non-Newtonian fluid are possible only up to a critical value of the viscoelastic fluid parameter, which decreases with an increase in R. It is further demonstrated that the perturbation solution gives acceptable results only if the viscoelastic fluid parameter is also small. Two more related problems are considered: fluid dynamics of a long porous slider, and injection of fluid through one side of a long vertical porous channel. For both the problems, exact numerical and other solutions are derived and appropriate conclusions drawn.     

A hybrid method for computing the flow of viscoelastic fluids

A hybrid method for computing the flow of viscoelastic and second-order fluids is presented. It combines the features of the finite difference technique and the shooting method. The method is accurate because it uses central differences. Its convergence is at least superlinear. The method is applied to obtain the solutions to three problems of flow of Walters' B' fluid: (a) flow near a stagnation point, (b) flow over a stretching sheet and (c) flow near a rotating disk. Numerical results reveal some new characteristics of flows which are not easy to demonstrate using the perturbation technique.     

Laminar flow through a channel with uniformly porous walls of different permeability

Summary The steady laminar flow of a viscous incompressible fluid through a two-dimensional channel, having fluid sucked or injected with different velocities through its uniformly porous parallel walls is considered. A solution for small suction Reynolds number has been given by the authors in a previous paper. The purpose of this paper is to present a solution valid for large Reynolds numbers for the cases of (i) suction at both walls, and (ii) suction at one wall and injection at the other. A technique of matching outer and inner expansions is used to obtain an asymptotic solution for both of these cases. Further a perturbation solution for the case of suction at one wall and injection at the other is obtained by choosing the difference between two wall velocities as the perturbation parameter. Both asymptotic and perturbation solutions are confirmed by exact numerical solutions. As expected, the resulting solutions show the presence of the usual suction boundary layers in both types of flow considered in this paper.     

Application of differential constraint method to exact solution of second-grade fluid

A differential constraint method is used to obtain analytical solutions of a second-grade fluid flow. By using the first-order differential constraint condition, exact solutions of Poiseuille flows, jet flows and Couette flows subjected to suction or blowing forces, and planar elongational flows are derived. In addition, two new classes of exact solutions for a second-grade fluid flow are found. The obtained exact solutions show that the non-Newtonian second-grade flow behavior depends not only on the material viscosity but also on the material elasticity. Finally, some boundary value problems are discussed.     

四角点支撑纳米板稳态受迫振动解析解

本文基于非局部弹性理论及辛叠加方法，得到放置在黏弹性介质上四角点支撑矩形纳米板稳态受迫振动问题的解析解．将纳米板受迫振动问题导入哈密顿体系，得到哈密顿控制方程，在无需任何预设函数的情况下可直接对哈密顿控制方程进行求解，得到简支纳米板稳态受迫振动问题在辛空间展开形式的解析解．进而通过边界叠加，可求出四角点支撑纳米板稳态受迫振动的解析解．数值算例中验证了本文应用辛叠加方法得到解析解的准确性，并以石墨烯纳米板为例，分析了非局部参数和黏弹性介质参数对四角点支撑石墨烯纳米板稳态受迫振动的影响．结果表明，非局部参数和黏弹性介质参数的变化会影响石墨烯纳米板的共振频率及共振幅值．     

非牛顿流体在渐变管中压力和剪切应力的二次摄动解

本文利用双摄动方法求解缓慢变化管道中Johnson-Segalman(J-S)流体流动的渐近解．将管道的扩张（或收缩）角度和粘弹性参数分别作为双摄动的参数，由流函数和涡量函数的形式，推导出压力和壁面剪切应力的渐近解．在此基础上，分析了管道角度，粘弹性参数和雷诺数等参数对压力以及剪切力影响．主要结论如下：(1) 管道扩张角度增加时，流向同一位置处径向压力以及壁面剪切应力随扩张角度减小；(2) 在同一扩张管道中，径向压力随着流向位移减小，收缩管与之相反；(3) 扩张角度与雷诺数对流场起主导作用，粘弹性系数起次要作用．     

非等温黏弹性复杂流动的改进SPH方法模拟

非等温黏弹性流体广泛存在于自然界和工业生产中,准确预测黏弹性流体的非等温流动机理和复杂流变特性有着重要的应用价值.文章提出一种改进的光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法对非等温黏弹性复杂流动进行了数值模拟,其中流体的黏弹特性通过eXtended Pom-Pom本构模型来表征.为了提高模拟结果的精度,采用了一种核函数梯度的修正算法;为了灵活地施加边界条件,发展了边界粒子和虚拟粒子相联合的边界处理方法;为了消除流动过程中的拉伸不稳定性,施加了粒子迁移技术.运用改进SPH方法数值模拟了液滴撞击固壁和F型腔注塑成型问题,通过与Basilisk软件得到的结果进行比较验证了改进SPH方法求解非等温黏弹性流体的有效性.通过利用不同粒子初始间距进行计算,评价了改进SPH方法的数值收敛性.研究了非等温流动相较于等温流动的不同流动特征,深入分析了不同热流变参数对流动过程的影响.数值结果表明,文章提出的改进SPH方法可稳定、准确地描述非等温黏弹性复杂流动的传热机理、复杂流变特性和自由面变化特性.     

复合材料层合板的三维研究

本文提出了固支复合材料各向异性层合圆板受均布横向载荷作用下的满足三维弹性力学基本微分方程和边界条件的解析解答。文中采用一种发展的摄动方法进行求解，板中的每个应力和位移都展开为无量纲厚度参数ε的摄动级数，并采用二维板理论解答作为其相应三维摄动解答的一个基本解的形式，通过摄动方法逐级求解而获得完整的三维解答。文中以解析形式和数值形式给出了高精确度的三维应力和位移结果，结果表明，本文求解三维问题的解析方法是合理有效的。     

Lighthill transpiration velocity revisited: An exact formulation

The Lighthill transpiration-velocity correction is commonly used in the analysis of viscous flows of aeronautical interest, in order to take into account the perturbation to the potential flow caused by the presence of the vorticity in the boundary layer. This correction consists of considering the boundary as a permeable surface from which the fluid flows through the boundary surface, with a velocity (named the transpiration velocity) determined on the basis of the local boundary-layer characteristics. Here, we use a new potential-vorticity decomposition in order to derive an exact representation of the effects of the vorticity on the external flow. The relationship between approximate transpiration-velocity representation and the exact one presented here is analyzed: it is shown that, under typical boundary-layer assumptions, the new representation reduces to that by Lighthill, except for a corrective field term. Finally, in order to quantify, in a simple case, the contribution of the corrective terms which arise in the new formulation, we examine, as a numerical test case, the problem of an attached boundary-layer flow over a flat plate: the numerical results indicate that the corrective term is negligible for Reynolds numbers above 10⁴.A preliminary version of this paper was presented at VI Convegno Italiano di Meccanica Computazionale, Brescia, Italy, October 1991.     

Steady Flows of Viscoelastic Fluids in Domains with Outlets to Infinity

The equations governing the motion of incompressible viscoelastic fluids of Rivlin—Ericksen and Oldroyd type are investigated in domains with cylindrical and paraboloidal outlets to infinity. For sufficiently small fluxes, prescribed in each outlet, existence and uniqueness of solutions are proven in weighted Hölder spaces. In domains with paraboloidal outlets the solution is obtained as a perturbation of the corresponding Navier—Stokes solution and in domains with cylindrical outlets as a perturbation of a flux carrier, constructed by joining together the exact solutions found in each outlet. These exact solutions are shown to be either rectilinear flows of Poiseuille type or flows composed of a rectilinear and of a transverse secondary component.     

A low cost laser-speckle photographic technique for velocity measurement in slow flows

A low cost, low power laser-speckle photographic technique has been developed and is duscussed for the measurement of point velocities in slow laminar flows. The technique is particularly suitable for axisymmetric flows where the two velocity components can be easily measured. The accuracy of the technique is established by measurement of the velocity distribution for Poiseuille flow and from data obtained for acceleration of an inelastic Newtonian fluid through a four-to-one circular contraction. Preliminary results are also presented in the contracting flow field for a non shear-thinning highly elastic fluid. These data are particularly significant for verification of finite element numerical solutions currently being developed for viscoelastic fluids in circular entry flows.     

Rotating electroosmotic flow of an Eyring fluid

A perturbation analysis is presented in this paper for the electroosmotic(EO) flow of an Eyring fluid through a wide rectangular microchannel that rotates about an axis perpendicular to its own. Mildly shear-thinning rheology is assumed such that at the leading order the problem reduces to that of Newtonian EO flow in a rotating channel, while the shear thinning effect shows up in a higher-order problem.Using the relaxation time as the small ordering parameter,analytical solutions are deduced for the leading-as well as first-order problems in terms of the dimensionless Debye and rotation parameters. The velocity profiles of the Ekman–electric double layer(EDL) layer, which is the boundary layer that arises when the Ekman layer and the EDL are comparably thin, are also deduced for an Eyring fluid. It is shown that the present perturbation model can yield results that are close to the exact solutions even when the ordering parameter is as large as order unity. By this order of the relaxation time parameter, the enhancing effect on the rotating EO flow due to shear-thinning Eyring rheology can be significant.     

Application of a vortex method to free surface flows

Vortex methods have found wide applications in various practical problems. The use of vortex methods in free surface flow problems, however, is still very limited. This paper demonstrates a vortex method for practical computation of non-linear free surface flows produced by moving bodies. The method is a potential flow formulation which uses the exact non-linear free surface boundary condition at the exact location of the instantaneous free surface. The position of the free surface, on which vortices are distributed, is updated using a Lagrangian scheme following the fluid particles on the free surface. The vortex densities are updated by the non-linear dynamic boundary condition, derived from the Euler equations, with an iterative Lagrangian numerical scheme. The formulation is tested numerically for a submerged circular cylinder in unsteady translation. The iteration is shown to converge for all cases. The results of the unsteady simulations agree well with classical linearized solutions. The stability of the method is also discussed.     

Asymptotic solutions of the flow of a Johnson-Segalman fluid through a slowly varying pipe using double perturbation strategy

A double perturbation strategy is presented to solve the asymptotic solutions of a Johnson-Segalman (J-S) fluid through a slowly varying pipe. First, a small parameter of the slowly varying angle is taken as the small perturbation parameter, and then the second-order asymptotic solution of the flow of a Newtonian fluid through a slowly varying pipe is obtained in the first perturbation strategy. Second, the viscoelastic parameter is selected as the small perturbation parameter in the second perturbation strategy to solve the asymptotic solution of the flow of a J-S fluid through a slowly varying pipe. Finally, the parameter effects, including the axial distance, the slowly varying angle, and the Reynolds number, on the velocity distributions are analyzed. The results show that the increases in both the axial distance and the slowly varying angle make the axial velocity slow down. However, the radial velocity increases with the slowly varying angle, and decreases with the axial distance. There are two special positions in the distribution curves of the axial velocity and the radial velocity with different Reynolds numbers, and there are different trends on both sides of the special positions. The double perturbation strategy is applicable to such problems with the flow of a non-Newtonian fluid through a slowly varying pipe.     

Extended homotopy perturbation method and the axisymmetric flow past a porous stretching sheet

The extended homotopy perturbation method, which is an extension of the celebrated homotopy perturbation method (HPM), is applied to obtain a solution to the problem of the steady, laminar, axisymmetric flow of a viscous, incompressible fluid past a porous stretching sheet. The solution so obtained is totally analytical and is expressible in terms of the cross‐flow velocity of the fluid past the stretching sheet. Its hallmark is that it does not depend upon computation of any auxiliary parameter for enlarging the convergence region of the solution. Rather, it calculates the solution automatically adjusting the scaling factor of the independent similarity variable normal to the sheet. The results obtained by the extended HPM are in excellent agreement with the exact numerical solution. Also, an asymptotic solution valid for large suction parameter is developed, which matches well with the exact solution even for moderate values of the suction parameter. Copyright © 2011 John Wiley & Sons, Ltd.     

HLLC scheme for the preconditioned pseudo-compressibility Navier–Stokes equations for incompressible viscous flows

A new HLLC (Harten-Lax-van leer contact) approximate Riemann solver with the preconditioning technique based on the pseudo-compressibility formulation for numerical simulation of the incompressible viscous flows has been proposed, which follows the HLLC Riemann solver (Harten, Lax and van Leer solver with contact resolution modified by Toro) for the compressible flow system. In the authors' previous work, the preconditioned Roe's Riemann solver is applied to the finite difference discretisation of the inviscid flux for incompressible flows. Although the Roe's Riemann solver is found to be an accurate and robust scheme in various numerical computations, the HLLC Riemann solver is more suitable for the pseudo-compressible Navier--Stokes equations, in which the inviscid flux vector is a non-homogeneous function of degree one of the flow field vector, and however the Roe's solver is restricted to the homogeneous systems. Numerical investigations have been performed in order to demonstrate the efficiency and accuracy of the present procedure in both two- and three-dimensional cases. The present results are found to be in good agreement with the exact solutions, existing numerical results and experimental data.     

Prediction of periodic boundary layers

A relatively simple, yet efficient and accurate finite difference method is developed for the solution of the unsteady boundary layer equations for both laminar and turbulent flows. The numerical procedure is subjected to rigorous validation tests in the laminar case, comparing its predictions with exact analytical solutions, asymptotic solutions, and/or experimental results. Calculations of periodic laminar boundary layers are performed from low to very high oscillation frequencies, for small and large amplitudes, for zero as well as adverse time-mean pressure gradients, and even in the presence of significant flow reversal. The numerical method is then applied to predict a relatively simple experimental periodic turbulent boundary layer, using two well-known quasi-steady closure models. The predictions are shown to be in good agreement with the measurements, thereby demonstrating the suitability of the present numerical scheme for handling periodic turbulent boundary layers. The method is thus a useful tool for the further development of turbulence models for more complex unsteady flows.     

Hydrodynamic Investigation of Certain Groundwater Flows in Coastal Pressurized Aquifers

Mathematical models of certain fresh groundwater flows through a semi-infinite pressurized aquifer toward a salt-water sea (basin, reservoir, etc.) are examined within the framework of two-dimensional steady-state flow theory. In order to study the models, mixed boundary value problems of the theory of analytic functions are formulated and solved using the Polubarinova-Kochina method. Algorithms for calculating the flows are developed on the basis of these models for situations in which groundwater streams arrive at the sea from one side or from below. Both the structure and the characteristic features of the simulated processes and the effect of all the physical characteristics of the models on the nature of the flow are analyzed in detail using the exact analytic dependences obtained and numerical calculations. The calculation results for the two inflow schemes are compared and features of the flows that depend on the initial position of the fluid contact are discussed.     

PTT黏弹性流体的光滑粒子动力学方法模拟

运用光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)方法对基于PHan-Thien-Tanner (PTT)模型的黏弹性流动进行了数值模拟. 首先, 利用SPH方法模拟了基于PTT模型的平板 Poiseuille流, 通过与文献结果的比较, 验证了SPH方法模拟黏弹性流动的准确性和有效性; 随后, 基于PTT模型对黏弹性自由表面流-液滴碰撞问题进行了SPH模拟, 研 究了PTT模型中拉伸参数对碰撞过程的影响. 为了解决张力不稳定问题, 采用简化的 人工应力公式. 数值结果表明, SPH方法可有效而灵活地模拟黏弹自由表面流问题.