Steady-state solutions for some simple flows of generalized Burgers fluids

The steady-state solutions for three types of unsteady oscillating flows of generalized Burgers fluids are determined by means of the Fourier sine transforms. These solutions are also presented in equivalent forms in terms of elementary functions exp, sine, cosine, hyperbolic sine and hyperbolic cosine. The similar solutions for Burgers, Oldroyd-B, Maxwell, Second grade and Navier-Stokes fluids can be also obtained as limiting cases of our solutions.     

On the unsteady unidirectional flows generated by impulsive motion of a boundary or sudden application of a pressure gradient

Some exact solutions of the time-dependent Navier-Stokes equations are discussed for flows impulsively started from rest by the motion of a boundary or two boundaries or by sudden application of a pressure gradient. It is shown that the expressions of the quantities such as velocity, flux and skin friction obtained in a series form for large times can also be used to find their values for small times or vice versa. Furthermore, conditions for which the values of the quantities at large times are nearly the same with those at small times are derived.     

On Solutions of Some Non-Linear Differential Equations Arising in Newtonian and Non-Newtonian Fluids

Some steady as well as unsteady solutions of the equations of motion for an incompressible Newtonian and non-Newtonian (second-grade) fluids are obtained by applying different methods including the Lie symmetry group method. The flows considered are axially symmetric with the swirling motion, and the governing equations for second-grade fluid flow have been modeled. Expressions for streamlines, velocity and vorticity components are constructed explicitly in each case. Exact analytical solutions in second-grade fluid are obtained and compared with the corresponding viscous solutions.     

An algorithm for rotating axisymmetric flows: model,validation and industrial applications

The study of axisymmetric flows is of interest not only from an academic point of view, due to the existence of exact solutions of Navier–Stokes equations, but also from an industrial point of view, since these kind of flows are frequently found in several applications. In the present work the development and implementation of a finite element algorithm to solve Navier–Stokes equations with axisymmetric geometry and boundary conditions is presented. Such algorithm allows the simulation of flows with tangential velocity, including free surface flows, for both laminar and turbulent conditions. Pseudo‐concentration technique is used to model the free surface (or the interface between two fluids) and the k–ε model is employed to take into account turbulent effects. The finite element model is validated by comparisons with analytical solutions of Navier–Stokes equations and experimental measurements. Two different industrial applications are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

A case of strong nonlinearity: Intermittency in highly turbulent flows

It has long been suspected that flows of incompressible fluids at large or infinite Reynolds number (namely at small or zero viscosity) may present finite time singularities. We review briefly the theoretical situation on this point. We discuss the effect of a small viscosity on the self-similar solution to the Euler equations for inviscid fluids. Then we show that single-point records of velocity fluctuations in the Modane wind tunnel display correlations between large velocities and large accelerations in full agreement with scaling laws derived from Leray's equations (1934) for self-similar singular solutions to the fluid equations. Conversely, those experimental velocity–acceleration correlations are contradictory to the Kolmogorov scaling laws.     

Unsteady flows of Oldroyd-B fluids in a channel of rectangular cross-section

The aim of this note is to present the exact solutions corresponding to two types of unsteady flows of an Oldroyd-B fluid in a channel of rectangular cross-section. The solutions that have been obtained satisfy both the associate partial differential equations and all imposed initial and boundary conditions. For λ_r or λ→0 they tend toward similar solutions for a Maxwell or second-grade fluid. If both λ_r and λ→0, the solutions for Navier-Stokes fluids are recovered.     

Theoretical analysis of non-uniform extensional flows in relation to processing rheology and predictions of some constitutive equations

In this paper, a characteristic equation involving the stream function, already given by one of the authors in a previous work for classifying axisymmetric incompressible flows, is re-considered. Non-uniform nearly extensional flows are derived as particular solutions from this equation. Using experimental data in the literature for polymer solutions and melts, it is proved that particular solutions of the characteristic equation lead to kinematics very close to those encountered in the fiber-spinning process. The kinematic equations satisfactorily correlating the fiber-spinning data are used in order to determine the ability of constitutive equations to predict realistic stresses in the flow domain. The rheological parameters of the fluids, obtained from experiments, are used for computation of differential and integral constitutive equations in the spinning conditions. Comparisons with the stress response of adequate constitutive equations are given and discussed.Also affiliated to: Université Joseph Fourier Grenoble I and Institut National Polytechnique de Grenoble, Associé au CNRS (URA 1510)     

On variational principles for coherent vortex structures

Different approaches are discussed of variational principles characterizing coherent vortex structures in two-dimensional flows. Turbulent flows seem to form ordered structures in the large scales of the motion and the self-organization principle predicts asymptotic states realizing an extremal value of the energy or a minimum of enstrophy. On the other hand the small scales take care of the increase of entropy, and asymptotic results can be obtained by applying the theory of equilibrium statistical mechanics.     

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.     

Comparing analytical and numerical solution of nonlinear two and three‐dimensional hydrostatic flows

New test cases for frictionless, three‐dimensional hydrostatic flows have been derived from some known analytical solutions of the two‐dimensional shallow water equations. The flow domain is a paraboloid of revolution and the flow is determined by the initial conditions, the nonlinear advective terms, the Coriolis acceleration and by the hydrostatic pressure. Wetting and drying is also included. Some specific properties of the exact solutions are discussed under different hypothesis and relative importance of the forcing terms. These solutions are proposed for testing the stability, the accuracy and the efficiency of numerical models to be used for simulating environmental hydrostatic flows. The computed solutions obtained with a semi‐implicit finite difference—finite volume algorithm on unstructured grid are compared with the corresponding analytical solutions in both two and three space dimension. Excellent agreement are obtained for the velocity and for the resulting water surface elevation. Comparison of the computed inundation area also shows a good agreement with the analytical solution with degrading accuracy observed when the inundation area becomes relatively large and for long simulation time. Copyright © 2006 John Wiley & Sons, Ltd.     

Global Analysis of the Flows of Fluids with Pressure-Dependent Viscosities

 To describe the flows of fluids over a wide range of pressures, it is necessary to take into account the fact that the viscosity of the fluid depends on the pressure. That the viscosity depends on the pressure has been verified by numerous careful experiments. While the existence of solutions local-in-time to the equations governing the flows of such fluids are available for small, special data and rather unrealistic dependence of the viscosity on the pressure, no global existence results are in place. Our interest here is to establish the existence of weak solutions for spatially periodic three-dimensional flows that are global in time, for a large class of physically meaningful viscosity-pressure relationships. (Accepted May 1, 2002) Published online November 15, 2002 Communicated by S. S. ANTMAN     

Numerical Analysis of the Influence of Thermal Stresses on the Nonlinear Thermomolecular Pressure Difference Effect

On the basis of a numerical analysis of the non-Navier-Stokes gas-dynamic equations for slow non-isothermal gas flows, the nonlinear thermomolecular pressure difference effect due to a large temperature gradient along the lateral surface of a capillary is investigated. It is shown that the magnitude of the effect is substantially different from the values calculated using the Navier-Stokes equations. For two models of molecular interaction (Maxwell molecules and hard spheres), the possibility of a quasi-one-dimensional interpretation of the effect for experimental estimation purposes is demonstrated. The solutions of the relaxation kinetic equation for flow in a plane capillary at small Knudsen numbers and the gas-dynamic equations for slow non-isothermal flows are compared and the range of their applicability is estimated.     

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.     

Creation of very-low-Reynolds-number chaotic fluid motions in microchannels using viscoelastic surfactant solution

Solutions of flexible high-molecular-weight polymers or some kinds of surfactant are viscoelastic fluids. The elastic stress is induced in such viscoelastic fluid flows and grows nonlinearly with the flow-rate resulting in many particular flow phenomena, including purely elastic instability. The purely elastic instability can even result in a kind of chaotic fluid motion, the so-called elastic turbulence, which is a recently discovered flow phenomenon and arises at arbitrarily small Reynolds number. By using viscoelastic surfactant solution, we attempted to create the peculiar chaotic fluid motions in several specially designed microchannels in which flows with curvilinear streamlines can be generated. The viscoelastic working fluids were aqueous solutions of surfactant, CTAC/NaSal (cetyltrimethyl ammonium chloride/sodium salicylate). CTAC solutions with weight concentration of 200 ppm (part per million) and 1000 ppm, respectively, at room temperature were tested. For comparison, water flows in the same microchannels were also visualized. The Reynolds numbers for all the microchannel flows were quite small (for solution flows, the Reynolds numbers were the order of or smaller than one) and the flow should be definitely laminar for Newtonian fluid. It was found that the regular laminar flow patterns for low-Reynolds-number Newtonian fluid flow in different microchannels were strongly deformed in solution flows: either asymmetrical flow structures or time-dependent vortical fluid motions appeared. These chaotic flow phenomena were considered to be induced by the viscoelasticity of the CTAC solutions. Discussions about the potential applications using such kind of chaotic fluid motions were also made.     

Properties and Averaged Equations for Flows of Bubbly Liquids

Averaged properties of bubbly liquids in the limit of large Reynolds and small Weber numbers are determined as functions of the volume fraction, mean relative velocity, and velocity variance of the bubbles using numerical simulations and a pair interaction theory. The results of simulations are combined with those obtained recently for sheared bubbly liquids [19] and the mixture momentum and continuity equations to propose a complete set of averaged equations and closure relations for the flows of bubbly liquids at large Reynolds and small Weber numbers.     

Unsteady flows of a generalized second grade fluid with the fractional derivative model between two parallel plates

The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced. Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus. The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates. In addition, the solutions of the shear stresses at the plates are also determined. The project supported by the National Natural Science Foundation of China (10372007, 10002003) and CNPC Innovation Fund     

Transient solutions for three-dimensional lid-driven cavity flows by a least-squares finite element method

A time-accurate least-squares finite element method is used to simulate three-dimensional flows in a cubic cavity with a uniform moving top. The time- accurate solutions are obtained by the Crank-Nicolson method for time integration and Newton linearization for the convective terms with extensive linearization steps. A matrix-free algorithm of the Jacobi conjugate gradient method is used to solve the symmetric, positive definite linear system of equations. To show that the least-squares finite element method with the Jacobi conjugate gradient technique has promising potential to provide implicit, fully coupled and time-accurate solutions to large-scale three-dimensional fluid flows, we present results for three-dimensional lid-driven flows in a cubic cavity for Reynolds numbers up to 3200.     

Computation of Viscoelastic flow in a cylindrical tank with a rotating lid

An orthogonal collocation method is used to compute steady flows of viscoelastic fluids in a cylindrical tank covered by a rotating disk. The stability of these flows with respect to small disturbances is also analyzed. The Criminale-Ericksen-Filbey constitutive equation is used, since the primary flow is viscometric and the secondary flow is small. The solutions satisfy the equation of continuity and the boundary conditions exactly, including the velocity discontinuity at the edge of the disk. The computed flows for aqueous solutions of Separan 30 exhibit single or double vortices, according to the concentration and the rotation speed. Reasonable agreement is found with the data of Hill [7,8].