A unified compatibility method for exact solutions of non-linear flow models of Newtonian and non-Newtonian fluids

Criteria are established for higher order ordinary differential equations to be compatible with lower order ordinary differential equations. Necessary and sufficient compatibility conditions are derived which can be used to construct exact solutions of higher order ordinary differential equations subject to lower order equations. We provide the connection to generalized groups through conditional symmetries. Using this approach of compatibility and generalized groups, new exact solutions of non-linear flow problems arising in the study of Newtonian and non-Newtonian fluids are derived. The ansatz approach for obtaining exact solutions for non-linear flow models of Newtonian and non-Newtonian fluids is unified with the application of the compatibility and generalized group criteria.     

On a new exact solution to Stokes’ first problem for Maxwell fluids

Stokes’ first problem for Maxwell fluids is re-examined using integral transform methods and several forms of the exact solution are presented. In the process, we call attention to, and correct, a recent case in the literature where this solution has been stated incorrectly. In addition, we note a number of analytical results and give several velocity profile plots. Finally, connections are drawn to other areas of research.     

Symmetry reductions and exact solutions to the ill-posed Boussinesq equation

In this paper, using the Lie symmetry analysis method, we study the ill-posed Boussinesq equation which arises in shallow water waves and non-linear lattices. The similarity reductions and exact solutions for the equation are obtained. Then the exact analytic solutions are considered by the power series method, and the physical significance of the solutions is considered from the transformation group point of view.     

Analytical solutions of singularities moving with an arbitrary path when two fluids are present

The derivations are carried out for the velocity potentials of singularities moving withan arbitrary path either in the upper fluid or in the lower fluid with or without a horizontalbottom when two fluids are present.In such a case,the pressure distribution is no longerequal to a constant or zero at the free interface.Taking the influence of an upper fluid uponthe lower fluid into consideration,a series of fundamental solutions in closed forms arepresented in this paper.     

Axial annular flow of plastic fluids: dead zones and plug-free flow

In axial annular flow, the shear stress decreases from its value τ(κR) at the inner cylinder to 0 at r = λR and increases from then on to τ(R) at the outer cylinder. For plastic fluids with a yield stress τ _c, λ will be such that flow commences when τ(κR) = τ(R) = τ _c. For fluids with position-dependent yield stresses (electro- and magnetorheological fluids are examples), the situation is more complex. While it is possible that yielding and flow occur everywhere, it is also possible that flow occurs only in parts of the fluid-filled space, and a dead zone (region in which the fluid is at rest) close to one of the walls exists. In that case, the fluid will flow no matter how small the applied pressure difference is. If P is large enough, the dead zone ceases to exist and flow without any plug is possible. The fluid flows as if no yield stress exists.

MHD Flow and heat transfer of non-Newtonian power-law fluid with diffusion and chemical reaction on a moving cylinder

The problem of flow and heat transfer of an electrically conducting non-Newtonian fluid over a continuously moving cylinder in the presence of a uniform magnetic field is analyzed for the case of power-law variation in the temperature and concentration at the cylinder surface. A diffusion equation with a chemical reaction source term is taken into account. The governing non-similar partial differential equation are solved numerically by employing shooting method. The effects of various parameters on the velocity, temperature and concentration profiles as well as the heat and mass transfer rate from the cylinder surface to the surrounding fluid are presented graphically and in tabulated form.     

Vortex shedding in cylinder flow of shear-thinning fluids: II. Flow characteristics

A detailed experimental study on the flow characteristics of various vortex shedding regimes was carried out for the flow of non-Newtonian fluids around a cylinder. The fluids were aqueous solutions of carboxymethyl cellulose (CMC) and tylose at weight concentrations ranging from 0.1 to 0.6%, which had varying degrees of shear-thinning and elasticity. Two cylinders of 10 and 20 mm diameter were used in the experiments, defining an aspect ratio of 12 and 6 and producing blockages of 5 and 10%, respectively. The Reynolds number (Re) ranged from 50 to 9×10³.Shear-thinning gave rise to a decrease of the cylinder boundary-layer thickness and to a reduction of the diffusion length (l_d), which raised the Strouhal number, St. In the laminar shedding regime, a modified Strouhal number was successful at overlapping the shedding frequency variation with the Reynolds number for the various solutions. In contrast, fluid elasticity was found to increase the formation length (l_f), and this contributed to a decrease of the Strouhal number. The overall effect of shear-thinning and elasticity was an increase in the Strouhal number.The increase in polymer concentration and the corresponding increase in fluid elasticity were responsible for the reduction of the critical Reynolds number marking the sudden decrease of the formation length, Re_lf. In the shear layer transition regime, the formation length and Strouhal number data collapsed onto single curves as function of a Reynolds number difference, which confirmed Coelho and Pinho (J. Non-Newtonian Fluid Mech. (2003), accepted for publication) finding that an important effect of fluid rheology was in changing the demarcations of the various flow regimes.     

Existence, uniqueness and stability of steady flows of second and third grade fluids in an unbounded “pipe-like” domain

We study steady motions of viscous incompressible third-grade fluids in unbounded channels with arbitrary shape. Such flows exist for small fluxes, due to a pressure drop. We prove that they are asymptotically stable in time, provided the viscosity is sufficiently large, and the initial condition on the perturbation sufficiently small.     

The Effect of Double Dispersion on Natural Convection Heat and Mass Transfer in a Non-Newtonian Fluid Saturated Non-Darcy Porous Medium

The effect of power law index parameter of the non-Newtonian fluid on free convection heat and mass transfer from a vertical wall is analyzed by considering double dispersion in a non-Darcy porous medium with constant wall temperature and concentration conditions. The Ostwald–de Waele power law model is used to characterize the non-Newtonian fluid behavior. In this case a similarity solution is possible. The variation of heat and mass transfer coefficients with the governing parameters such as power law index, thermal and solutal dispersion parameters, inertia parameter, buoyancy ratio, and the Lewis number is discussed for a wide range of values of these parameters.     

Displacement of a Newtonian fluid by a non-Newtonian fluid in a porous medium

This paper presents an analytical Buckley-Leverett-type solution for one-dimensibnal immiscible displacement of a Newtonian fluid by a non-Newtonian fluid in porous media. The non-Newtonian fluid viscosity is assumed to be a function of the flow potential gradient and the non-Newtonian phase saturation. To apply this method to field problems a practical procedure has been developed which is based on the analytical solution and is similar to the graphic technique of Welge. Our solution can be regarded as an extension of the Buckley-Leverett method to Non-Newtonian fluids. The analytical result reveals how the saturation profile and the displacement efficiency are controlled not only by the relative permeabilities, as in the Buckley-Leverett solution, but also by the inherent complexities of the non-Newtonian fluid. Two examples of the application of the solution are given. One application is the verification of a numerical model, which has been developed for simulation of flow of immiscible non-Newtonian and Newtonian fluids in porous media. Excellent agreement between the numerical and analytical results has been obtained using a power-law non-Newtonian fluid. Another application is to examine the effects of non-Newtonian behavior on immiscible displacement of a Newtonian fluid by a power-law non-Newtonian fluid.     

Exact and numerical solutions of a fully developed generalized second-grade incompressible fluid with power-law temperature-dependent viscosity

We compute exact and numerical solutions of a fully developed flow of a generalized second-grade fluid, with power-law temperature-dependent viscosity (μ=θ^-M), down an inclined plane. Analytical solutions are found for the case when M=m+1, m≠1, m being a constant that models shear thinning (m<0) or shear thickening (m>0). The exact solutions are given in terms of Bessel functions. The numerical solutions indicate that both the velocity and the temperature increase with decreasing Froude number and that there is a critical value of Fr below which temperature “overshoots” its free surface value of unity. This phenomena is not reported in the work of Massoudi and Phuoc [Fully developed flow of a modified second grade fluid with temperature dependent viscosity, Acta Mech. 150 (2001) 23-37.] for viscosity that depends exponentially on temperature.     

A robust algorithm for moving interface of multi-material fluids based on Riemann solutions

In the paper, the numerical simulation of interface problems for multiple material fluids is studied. The level set function is designed to capture the location of the material interface. For multi-dimensional and multi-material fluids, the modified ghost fluid method needs a Riemann solution to renew the variable states near the interface. Here we present a new convenient and effective algorithm for solving the Riemann problem in the normal direction. The extrapolated variables are populated by Taylor series expansions in the direction. The anti-diffusive high order WENO difference scheme with the limiter is adopted for the numerical simulation. Finally we implement a series of numerical experiments of multi-material flows. The obtained results are satisfying, compared to those by other methods.The English text was polished by Boyi Wang.     

The yield surface of viscoelastic and plastic fluids in a vane viscometer

Vane viscometers are often used to investigate the low shear rate properties of plastic fluids. The shear stress is determined by assuming that the material is held in the space between the vane blades so that it behaves like a rigid cylinder. Experimental evidence supports this assumption and the aim of the present study is to model numerically the yield process in a vane rheometer using viscoelastic and plastic fluids. The finite element method has been used to model the behavior of Herschel-Bulkley (Bingham), Casson and viscoelastic (Maxwell type) fluids. The penalty function approach for the pressure approximation and a rotating reference frame are used together with fine meshes containing more than 1300 elements. The results show that for Herschel-Bulkley (Bingham), and Casson fluids a rotating rigid cylinder of fluid is trapped inside the periphery of the vane, the shear stress is uniformly distributed over the surface of the cylinder. Finally a modified second order fluid is used to simulate the viscoelastic behaviour, anticipated to be an intermediate between the elastic deformation and the plastic flow, to provide a more realistic simulation of the yield process about a vane. In this case, contrast with the concentration of the elastic strain rate at the blade tips, a nearly uniform distribution of the plastic shear rate is still found. This implies that the plastic shear always distributes uniformly during the entire yielding process. Evidently the assumption of uniform shear on a rotating cylinder of material occluded in the blades of a vane is a valid and useful model for many types of fluid possessing a yield stress.     

Symmetry groups and similarity solutions for a free convective boundary-layer problem

The free convective boundary-layer problem due to the motion of an elastic surface into an electrically conducting fluid is studied with group-theoretical methods. The symmetry groups admitted by the corresponding boundary value problem are obtained. Particular attention is paid on the group of scaling which provides the similarity solution of the problem. Also, the admissible form of the data, in order to be conformed to the obtained symmetries, is provided. Finally, with the use of the entailed similarity solution the problem is transformed into a boundary value problem of ODEs and is solved numerically.     

Similarity solutions for creeping flow and heat transfer in second grade fluids

The two-dimensional equations of motions for the slowly flowing and heat transfer in second grade fluid are written in cartesian coordinates neglecting the inertial terms. When the inertia terms are simply omitted from the equations of motions the resulting solutions are valid approximately for Re?1. This fact can also be deduced from the dimensionless form of the momentum and energy equations. By employing Lie group analysis, the symmetries of the equations are calculated. The Lie algebra consist of four finite parameter and one infinite parameter Lie group transformations, one being the scaling symmetry and the others being translations. Two different types of solutions are found using the symmetries. Using translations in x and y coordinates, an exponential type of exact solution is presented. For the scaling symmetry, the outcoming ordinary differential equations are more involved and only a series type of approximate solution is presented. Finally, some boundary value problems are discussed.     

Bifurcations of traveling wave solutions and exact solutions to generalized Zakharov equation and Ginzburg-Landau equation

This paper studies the dynamic behaviors of some exact traveling wave solutions to the generalized Zakharov equation and the Ginzburg-Landau equation. The effects of the behaviors on the parameters of the systems are also studied by using a dynamical system method. Six exact explicit parametric representations of the traveling wave solutions to the two equations are given.     

An exact analytical solution of the Falkner-Skan equation with mass transfer and wall stretching

In this paper, an exact analytical solution of the famous Falkner-Skan equation is obtained. The solution involves the boundary layer flow over a moving wall with mass transfer in presence of a free stream with a power-law velocity distribution. Multiple solution branches are observed. The effects of mass transfer and wall stretching are analyzed. Interesting velocity profiles including velocity overshoot and reversal flows are observed in the presence of both mass transfer and wall stretching. These solutions greatly enrich the analytical solution for the celebrated Falkner-Skan equation and the understanding of this important and interesting equation.     

Backlund transformation and exact solutions for whitham-broer-kaup equations in shallow water

I.IntroductionTodescribethepropagationofshallowwaterwave,manywell-knowncompletelyintegrablemodelsareintroduced,suchasKdVequatioll,Boussinesqequation,K-Pequation,WBKequation,etc.UnderBoussinesqapproximation,Whitham,BroerandKaupl"2'3]obtainednon-lilleurWBKequationwhereu=u(x,l)isthefieldofhorizontalvelocity;v=v(x,t)istheheightthatdeviatefromequilibriumpositionofliquid;a,gareconstantsthatrepresentdift'erentdispersivepower.TheEqs.(1.I),(l.2)areverygoodmodelstodescl.ibedispersivewave.Ifa=0,P…     

Analytical and numerical solutions of electro-osmotically driven flow of a third grade fluid between micro-parallel plates

Electro-osmotic flow of a third grade fluid between micro-parallel plates is considered. The equations of motion are derived and made dimensionless. Approximate analytical solutions are obtained by perturbation techniques. Constant viscosity and temperature dependent viscosity (Reynolds model) cases are treated separately. Numerical solutions of the equations are also obtained. Influences of non-Newtonian parameter, Joule heating effect, viscosity index and electro-kinetic effect on the velocity and temperature profiles are shown. Approximate and numerical solutions are contrasted.     

The arbitrarily and the periodically laminated elastic hollow sphere: exact solutions and homogenization

Summary The aim of this paper is (1) to develop a rational method for the analysis of an arbitrarily laminated elastic, isotropic or transversely isotropic hollow sphere under internal and/or external pressure, (2) to solve the problem of a periodically layered sphere consisting of many equal groups of n different thin layers. The transfer matrix method is used, and exact closed-form solutions are worked out, supplemented by a numerical example. It turns out that by means of the proposed homogenization an originally (periodically) inhomogeneous isotropic sphere is replaced by a homogeneous anisotropic one belonging to the type of spherical symmetric anisotropy. Received 20 October 1997; accepted for publication 15 December 1997