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.     

Homotopy analysis solution for micropolar fluid flow through porous channel with expanding or contracting walls of different permeabilities

The flow of a micropolar fluid through a porous channel with expanding or contracting walls of different permeabilities is investigated. Two cases are considered, in which opposing walls undergo either uniform or non-uniform motion. In the first case,the homotopy analysis method (HAM) is used. to obtain the expressions for the velocity and micro-rotation fields. Graphs are sketched for some parameters. The results show that the expansion ratio and the different permeabilities have important effects on the dynamic characteristics of the fluid. Following Xu's model, in the second case which is more general, the wall expansion ratio varies with time. Under this assumption, the governing equations are transformed into nonlinear partial differential equations that can also be solved analytically by the HAM. In the process, both algebraic and exponential models are considered to describe the evolution of α(t) from the initial state α0 to the final state α1. As a result, the time-dependent solutions are found to approach the steady state very rapidly. The results show that the time-dependent variation of the wall expansion ratio can be ignored because of its limited effects.     

On the comparison of two different solutions in the form of series of the governing equation of an unsteady flow of a second grade fluid

In this paper, two different solutions in the form of series of the governing equation of unsteady flow of a second grade fluid are considered. These are series expansions with respect to inverse power of time and a perturbation expansion. Two illustrative examples are given. One of them is the unsteady flow of a second grade fluid over a plane wall suddenly set in motion and the other is the diffusion of a line vortex in a fluid of second grade. It is a remarkable fact that the expression of the series expansion with respect to inverse power of time is exactly in the same form as that of the perturbation expansion. Thus, it is possible to replace a series expansion with respect to inverse power of time with a perturbation expansion.     

Analytic solution for flow of Sisko fluid through a porous medium

This paper presents the analytic solution for flow of a magnetohydrodynamic (MHD) Sisko fluid through a porous medium. The non-linear flow problem in a porous medium is formulated by introducing the modified Darcy’s law for Sisko fluid to discuss the flow in a porous medium. The analytic solutions are obtained using homotopy analysis method (HAM). The obtained analytic solutions are explicitly expressed by the recurrence relations and can give results for all the appropriate values of material parameters of the examined fluid. Moreover, the well-known solutions for a Newtonian fluid in non-porous and porous medium are the limiting cases of our solutions.     

Unsteady mixed convective heat transfer of two immiscible fluids confined between long vertical wavy wall and parallel flat wall

The combined effects of thermal and mass convection of viscous incompressible and immiscible fluids through a vertical wavy wall and a smooth flat wall are analyzed. The dimensionless governing equations are perturbed into a mean part (the zeroth-order) and a perturbed part (the first-order). The first-order quantities are obtained by the perturbation series expansion for short wavelength, in which the terms of the exponential order arise. The analytical expressions for the zeroth-order, the first-order, and the total solutions are obtained. The numerical computations are presented graphically to show the salient features of the fluid flow and the heat transfer characteristics. Separate solutions are matched at the interface by using suitable matching conditions. The shear stress and the Nusselt number are also analyzed for variations of the governing parameters. It is observed that the Grashof number, the viscosity ratio, the width ratio, and the conductivity ratio promote the velocity parallel to the flow direction. A reversal effect is observed for the velocity perpendicular to the flow direction.     

Stokes＇ first problem for a viscoelastic fluid with the generalized Oldroyd-B model

The flow near a wall suddenly set in motion for a viscoelastic fluid with the generalized Oldroyd-B model is studied. The fractional calculus approach is used in the constitutive relationship of fluid model. Exact analytical solutions of velocity and stress are obtained by using the discrete Laplace transform of the sequential fractional derivative and the Fox H-function. The obtained results indicate that some well known solutions for the Newtonian fluid, the generalized second grade fluid as well as the ordinary Oldroyd-B fluid, as limiting cases, are included in our solutions. The project supported by the National Natural Science Foundation of China (10272067), the Doctoral Program Foundation of the Education Ministry of China (20030422046), the Natural Science Foundation of Shandong Province, China (Y2006A14) and the Research Foundation of Shandong University at Weihai. The English text was polished by Keren Wang.     

Laws of Incompressible-Fluid Turbulent Near-Wall Flows Following from the Limiting Asymptotic Properties

The averaged layered turbulent wall flows of an incompressible fluid are considered for an arbitrary wall roughness and a friction coefficient tending to zero. It is assumed that in the asymptotic limit the flow may be divided into two regions: an outer, only slightly disturbed region, and a thin wall region. A similarity law for the velocity disturbances in the outer region is established for cylindrical tubes with an arbitrary cross section. Using the asymptotically matched solutions in the inner and outer regions, the existence of a logarithmic velocity profile in the matching zone is proved. On the basis of a unique asymptotic approach, the existing semi-empirical solutions for tube flows are brought into consistency. A number of semi-empirical formulas for the velocity profiles, which refine and expand the range of applicability of existing relations, are proposed.     

Axisymmetric flow and heat transfer of the Carreau fluid due to a radially stretching sheet: Numerical study

The prime objective of this article is to study the axisymmetric flow and heat transfer of the Carreau fluid over a radially stretching sheet. The Carreau constitutive model is used to discuss the characteristics of both shear-thinning and shear-thickening fluids. The momentum equations for the two-dimensional flow field are first modeled for the Carreau fluid with the aid of the boundary layer approximations. The essential equations of the problem are reduced to a set of nonlinear ordinary differential equations by using local similarity transformations. Numerical solutions of the governing differential equations are obtained for the velocity and temperature fields by using the fifth-order Runge–Kutta method along with the shooting technique. These solutions are obtained for various values of physical parameters. The results indicate substantial reduction of the flow velocity as well as the thermal boundary layer thickness for the shear-thinning fluid with an increase in the Weissenberg number, and the opposite behavior is noted for the shear-thickening fluid. Numerical results are validated by comparisons with already published results.     

Hydromagnetic flow through uniform channel bounded by porous media

The combined effects of the magnetic field, permeable walls, Darcy velocity, and slip parameter on the steady flow of a fluid in a channel of uniform width are studied. The fluid flowing in the channel is assumed to be homogeneous, incompressible,and Newtonian. Analytical solutions are constructed for the governing equations using Beavers-Joseph slip boundary conditions. Effects of the magnetic field, permeability,Darcy velocity, and slip parameter on the axial velocity, slip velocity, and shear stress are discussed in detail. It is shown that the Hartmann number, Darcy velocity, porous parameter, and slip parameter play a vital role in altering the flow and in turn the shear stress.     

A high-accuracy temporal–spatial pseudospectral method for time-periodic unsteady fluid flow and heat transfer problems

A temporal–spatial pseudospectral (TSP) method is proposed for the high-accuracy solutions of time-periodic unsteady fluid flow and heat transfer problems. In this method, both the spatial and temporal derivative terms in the governing equations are computed by pseudospectral method. The spatial derivatives are computed through Chebyshev and Lagrange polynomials while the time derivatives are computed by Fourier series. The TSP method is capable of directly finding out the periodic state solutions without the necessity to resolve the initial transient state solutions, hence holds high computational efficiency and high numerical accuracy properties for the time-periodic problems. This method is validated by three 2D benchmark problems: the time-periodic incompressible flow with exact solutions; the natural convection in enclosure with time-periodic temperature on one sidewall, and on both sidewalls. The TSP results fit well the exact solutions or the benchmark solutions and the TSP accuracy is much higher than the time marching spatial pseudospectral accuracy. Some time-dependent fluid flow and heat transfer characteristic parameters are analysed. The proposed TSP method could be further extended to more complex time-periodic unsteady fluid flow and heat transfer problems where high-accuracy results are required.     

Gravity-capillary waves in the presence of constant vorticity

Periodic and solitary gravity-capillary waves propagating at a constant velocity at the surface of a fluid of finite depth are considered. The vorticity in the fluid is assumed to be constant. Analytical solutions are presented for waves of small amplitude. For waves of large amplitude, numerical solutions are computed by boundary integral equation methods. The results unify previous findings for irrotational gravity capillary waves and gravity waves with constant vorticity. In particular solitary waves with oscillatory tails and branches of solutions which exist only for waves of large amplitude are found.     

爆炸冲击波绕流的三维数值模拟研究

对爆炸冲击波绕流问题采用可压缩流体模型建立了质量、动量和能量的偏微分守恒方程组,基于多物质流体Euler型算法, 采用算子分裂格式, 运用体积份额法处理多介质界面, 用自行开发的三维数值模拟程序MMIC3D模拟计算了三维空中爆炸的爆源附近和爆炸场中冲击波流场的发展规律, 研究了挡墙拐角处绕流的形成和变化情况, 分析比较了挡墙的位置、形状对挡墙后爆炸冲击波的影响, 并与经验公式进行了比较研究. 数值计算基本符合物理规律,说明文中采用的模型和算法是合理的, 为工程设计减压设施提供了数值计算依据．      

Computation of compressible flows with high density ratio and pressure ratio

The WENO method, RKDG method, RKDG method with original ghost fluid method, and RKDG method with modified ghost fluid method are applied to singlemedium and two-medium air-air, air-liquid compressible flows with high density and pressure ratios： We also provide a numerical comparison and analysis for the above methods. Numerical results show that, compared with the other methods, the RKDG method with modified ghost fluid method can obtain high resolution results and the correct position of the shock, and the computed solutions are converged to the physical solutions as themesh is refined.     

Magnetohydrodynamic transient flows of a non-Newtonian fluid

Some exact solutions of the time-dependent partial differential equations are discussed for flows of an Oldroyd-B fluid. The fluid is electrically conducting and incompressible. The flows are generated by the impulsive motion of a boundary or by application of a constant pressure gradient. The method of Laplace transform is applied to obtain exact solutions. It is observed from the analysis that the governing differential equation for steady flow in an Oldroyd-B fluid is identical to that of the viscous fluid. Several results of interest are obtained as special cases of the presented analysis.     

Dynamic coupling between horizontal vessel motion and two-layer shallow-water sloshing

Numerical and analytical results are presented for fluid sloshing, of a two-layer inviscid, incompressible and immiscible fluid with thin layers and a rigid lid, coupled to a vessel which is free to undergo horizontal motion governed by a nonlinear spring. Exact analytical results are obtained for the linear problem, giving the natural frequencies and the resonance structure, particularly between the fluid and vessel. A numerical method for the linear and nonlinear equations is developed based on the high-resolution f-wave-propagation finite volume methods due to Bale et al. (2002) [SIAM Journal on Scientific Computing 24, 955–978], adapted to include the pressure gradient at the rigid-lid, and coupled to a Runge–Kutta solver for the vessel motion. The numerical simulations in the linear limit are compared with the exact analytical solutions. The coupled nonlinear numerical solutions with simulations near the internal 1:1 resonance are presented. Of particular interest is the partition of energy between the vessel and fluid motion.     

Study on the 3D Green's functions of the fluid and piezoelectric bimaterials

Because most piezoelectric devices have interfaces with fluid in engineering,it is valuable to study the coupled field between fluid and piezoelectric media.As the fundamental problem,the 3 D Green's functions for point forces and point charge loaded in the fluid and piezoelectric bimaterials are studied in this paper.Based on the 3 D general solutions expressed by harmonic functions,we constructed the suitable harmonic functions with undetermined constants at first.Then,the couple field in the fluid and piezoelectric bimaterials can be derived by substitution of harmonic functions into general solutions.These constants can be obtained by virtue of the compatibility,boundary,and equilibrium conditions.At last,the characteristics of the electromechanical coupled fields are shown by numerical results.     

New analytical solutions to water wave radiation by vertical truncated cylinders through multi-term Galerkin method

New analytical solutions to water wave radiation by vertical truncated circular cylinders are developed based on linear potential flow theory. Two typical cylinder configurations of a surface-piercing cylinder and a submerged floating cylinder are considered. The multi-term Galerkin method is employed in the solution procedure, in which the fluid velocity on the interface between different regions is expanded into a set of basis function involving the Gegenbauer polynomials, and the cube-root singularity of fluid velocity at the side edges of the truncated cylinders is correctly modeled. The present solutions have the merits of very rapid convergence. The results with six-figure accuracy for added mass and radiation damping can be obtained using a few truncated numbers in the basis function for three motions (surge, heave and roll). The calculated results of the present solutions agree well with that by a higher-order boundary element method solution. Calculation examples are presented to investigate the influence of the motion frequency on the added mass and the radiation damping of the truncated cylinders with different geometric parameters. The present solutions can be used as a reliable benchmark for numerical solutions to water wave radiation by complicated structures.

     

Propagation of a penny-shaped fluid-driven fracture in an impermeable rock: asymptotic solutions

This paper presents an analysis of the propagation of a penny-shaped hydraulic fracture in an impermeable elastic rock. The fracture is driven by an incompressible Newtonian fluid injected from a source at the center of the fracture. The fluid flow is modeled according to lubrication theory, while the elastic response is governed by a singular integral equation relating the crack opening and the fluid pressure. It is shown that the scaled equations contain only one parameter, a dimensionless toughness, which controls the regimes of fracture propagation. Asymptotic solutions for zero and large dimensionless toughness are constructed. Finally, the regimes of fracture propagation are analyzed by matching the asymptotic solutions with results of a numerical algorithm applicable to arbitrary toughness.     

A variational approach to non-steady non-newtonian flow in a circular pipe

The transient response of a non-Newtonian power-law fluid to several assumed forms of pressure pulse in a circular tube is analysed by the semi-direct variational method of Kanntovorich. Velocity profiles are shown for several power-law indices, and by comparing the results for the Newtonian case with the exact solution given by Szymanski, it is observed that the results are good to 5%. More accurate solutions have been found for the case involving Newtonian fluid flow. New results are reported concerning the effect of a triangular pressure pulse on the development and transient response of the flow field of a non-Newtonian fluid.     

Nonsteady flow of groundwater in aquifer system with viscoelastic properties

The transient flow behavior of groundwater in aquifer-aquitard system with viscoelastic properties is studied. On the basis of previous works (Hantush, Neuman, Brutsaert, Corapcioglu), the new partial differential-integral equations are derived. The well-known equations (Hantush, Brutsaert) are the special cases of the new equations. The new equations describe the flow of a slightly compressible groundwater in layers with viscoelastic properties.Analytical solutions of the partial differential-integral equations are obtained by using the method of Laplace transform. The viscoelastic properties enhance the heterogeneities of elastic aquifer system which have delay and feed qualities. The agreements between the numericl inversion results of Laplace transform and the analytical solutions are good. The formulae predict the transient flow behavior of groundwater in this heterogeneous layers.