Equibiaxial extensional flow of polymer melts via lubricated squeezing flow. II. Flow modeling

A model for lubricated squeezing flow of a viscoelastic fluid is developed in order to study the viability of this flow as a rheological technique for generating equibiaxial extensional deformations in polymer melts. In this simple flow model, the melt, described by an upper-convected Maxwell fluid, is squeezed between thin films of a Newtonian fluid. Comparisons of the model predictions for constant strain rate and constant stress flows are made with experimental results presented in the first paper. Predictions from the model are able to describe the effects of lubricant viscosity and experimental configuration and indicate the technique fails for these flows at Hencky strains of approximately one. The cause for this failure is lubricant thinning, which leads to significant errors in both the measured stress difference and the strain. Received: 31 January 2000 Accepted: 31 May 2000     

A Statistical Model for Predicting the Heat Transfer of Solid Particles in Turbulent Flows

The objective of this paper is twofold: (1) to present a statistical model of particle transport and heat transfer in turbulent flows and (2) to examine the performance of this model in various turbulent flows going from a simple flow to a more complicated one. This model is based on a kinetic equation for the probability density function of the particle velocity and temperature distributions in anisotropic turbulent flow. The model predictions compare reasonable well with numerical simulations and properly reproduce the crucial trends of computations performed in various turbulent flows.     

On unsteady unidirectional flows of a second grade fluid

Some properties of unsteady unidirectional flows of a fluid of second grade are considered for flows produced by the sudden application of a constant pressure gradient or by the impulsive motion of one or two boundaries. Exact analytical solutions for these flows are obtained and the results are compared with those of a Newtonian fluid. It is found that the stress at the initial time on the stationary boundary for flows generated by the impulsive motion of a boundary is infinite for a Newtonian fluid and is finite for a second grade fluid. Furthermore, it is shown that initially the stress on the stationary boundary, for flows started from rest by sudden application of a constant pressure gradient is zero for a Newtonian fluid and is not zero for a fluid of second grade. The required time to attain the asymptotic value of a second grade fluid is longer than that for a Newtonian fluid. It should be mentioned that the expressions for the flow properties, such as velocity, obtained by the Laplace transform method are exactly the same as the ones obtained for the Couette and Poiseuille flows and those which are constructed by the Fourier method. The solution of the governing equation for flows such as the flow over a plane wall and the Couette flow is in a series form which is slowly convergent for small values of time. To overcome the difficulty in the calculation of the value of the velocity for small values of time, a practical method is given. The other property of unsteady flows of a second grade fluid is that the no-slip boundary condition is sufficient for unsteady flows, but it is not sufficient for steady flows so that an additional condition is needed. In order to discuss the properties of unsteady unidirectional flows of a second grade fluid, some illustrative examples are given.     

Numerical investigation of fully developed turbulent fluid flow and heat transfer in a square duct

Three-dimensional fully developed turbulent fluid flow and heat transfer in a square duct are numerically investigated with the author's anisotropic low-Reynolds-number k-ε turbulence model. Special attenton has been given to the regions close to the wall and the corner, which are known to influence the characteristics of secondary flow a great deal. Hence, instead of the common wall function approach, the no-slip boundary condition at the wall is directly used. Velocity and temperature profiles are predicted for fully developed turbulent flows with constant wall temperature. The predicted variations of both local wall shear stress and local wall heat flux are shown to be in close agreement with available experimental data. The present paper also presents the budget of turbulent kinetic energy equation and the systematic evaluation for existing wall function forms. The commonly adopted wall function forms that are valid for two-dimensional flows are found to be inadequate for three-dimensional turbulent flows in a square duct.     

Exact solutions for the unsteady flow of a Burger’s fluid in a duct induced by time-dependent prescribed volume flow rate

In this investigation, some unsteady flows in a circular duct have been studied. The fluid obeys viscoelastic non-Newtonian model with the Burgers’ constitutive equation and all fluid properties are constant. The flows in a duct are due to the prescribed arbitrary time dependent inlet volume flow rates. Four types of flow situations are considered. The governing equations are first developed and then solved using Laplace transform technique. Results indicate the strong effect of Burgers’ fluid parameter on the velocity fields and pressure gradients.     

Laws of flow with a limiting gradient in anisotropic porous media

The equations of viscoplastic fluid flow through a porous medium are written for all types of anisotropy. It is shown that in anisotropic media the flows with a limiting gradient are characterized by two material tensors: the tensor of permeability (flow resistance) coefficients and the tensor of limiting gradients. A complex of laboratory measurements for determining the tensors of permeability coefficients and limiting gradients is considered for all types of anisotropic media. It is shown that the tensors of permeability coefficients and limiting gradients are coaxial. Conditions of flow onset and fluid flow laws are formulated for media with monoclinic and triclinic symmetries of flow characteristics.     

Leak-off During Placement of Cr(III)-Partially Hydrolyzed Polyacrylamide Gelling Solution in Fractured Porous Media

During waterflooding of a fractured formation, water may channel through the fracture or interconnected network of fractures, leaving a large portion of oil bearing rock unswept. One remedial practice is injection of a gelling solution into the fracture. Such placement of a gelling mixture (referred as gelant) is associated with leak-off from the face of the fracture into the adjoining matrix. As the gelant gets more crosslinked, the gelant encounters more resistance in flowing into the porous matrix. This article addresses the build-up of flow resistance as the Cr(III)-partially hydrolyzed polyacrylamide gelant, at various stages of crosslinking flows into the matrix. Flow experiments were conducted at constant injection pressure in unfractured Berea rocks that represent a matrix adjoining a fracture. Before entering the core, gelants underwent post-mixing delays, shorter than their gel time. On continued displacement, flow resistance developed that reduced the flow rate further. More delay, after mixing of gelant hastened, the build-up of resistance to flow and the resistance was contained nearer to the inlet face. Effect of flow over fracture face on the build-up of flow resistance in the matrix was also evaluated by conducting displacement of gelant in two fractured slabs. In one case, a part of the injected fluid came out of the fracture outlet with the rest leaking off into matrix. In the other case, all the fluid that entered into the fracture leaked off into the matrix. Build-up of flow resistances in the matrix for the two cases was compared. A simple conceptual model is presented that could explain the flow of gelant and build-up of resistance in porous rock at constant injection pressure.     

A film flow model for analysing gravity-driven,thin wavy fluid films

To analyse the physics underlying gravity-driven runoff of thin wavy films, a film flow model is developed, and is solved with computational fluid dynamics. This model is based on the lubrication theory, and takes into account the gravitational, wall shear and surface tension forces. A key characteristic of the model is that it assumes only one computational cell over the film height, which enables studying film flow on larger computational domains. A main aim of this study is to perform a detailed validation of the numerical model. The film flow model is validated against several experiments of gravity-driven, thin fluid films on smooth surfaces. The time-averaged film thickness and the fluid speed profiles predicted by the model show very good agreement with experimental results. Similarly, the film flow model is able to predict the wave speeds with sufficient accuracy. The energy spectra of the waves, where higher frequency waves are present in film flows at higher Reynolds numbers, show an exponentially decaying trend at these high frequencies. The model performs better than the Nusselt equation for film flows, which under-predicts the time-averaged film thickness and over-predicts the time-averaged fluid speeds, even for flows at low Reynolds numbers. The film flow model is compared qualitatively for fingering behaviour. This model also allows to investigate film flows on large surfaces, which can be rough, curved and of complex geometrical shape.     

Exact solutions for generalized Burgers’ fluid in an annular pipe

This paper deals with some unsteady unidirectional transient flows of generalized Burgers’ fluid in an annular pipe. Exact solutions of some unsteady flows of generalized Burgers’ fluid in an annular pipe are obtained by using Hankel transform and Laplace transform. The following two problems have been studied: (1) Poiseuille flow due to a constant pressure gradient; (2) axial Couette flow in a annulus. The well known solutions for Navier-Stokes fluid, as well as those corresponding to a Maxwell fluid, a second grade fluid and an Oldroyd-B fluid appear as limiting cases of our solutions.     

An unsymmetric constitutive equation for anisotropic viscoelastic fluid

A continuum constitutive theory of corotational derivative type is developed for the anisotropic viscoelastic fluid–liquid crystalline (LC) polymers. A concept of anisotropic viscoelastic simple fluid is introduced. The stress tensor instead of the velocity gradient tensor D in the classic Leslie–Ericksen theory is described by the first Rivlin–Ericksen tensor A and a spin tensor W measured with respect to a co-rotational coordinate system. A model LCP-H on this theory is proposed and the characteristic unsymmetric behaviour of the shear stress is predicted for LC polymer liquids. Two shear stresses thereby in shear flow of LC polymer liquids lead to internal vortex flow and rotational flow. The conclusion could be of theoretical meaning for the modern liquid crystalline display technology. By using the equation, extrusion–extensional flows of the fluid are studied for fiber spinning of LC polymer melts, the elongational viscosity vs. extension rate with variation of shear rate is given in figures. A considerable increase of elongational viscosity and bifurcation behaviour are observed when the orientational motion of the director vector is considered. The contraction of extrudate of LC polymer melts is caused by the high elongational viscosity. For anisotropic viscoelastic fluids, an important advance has been made in the investigation on the constitutive equation on the basis of which a series of new anisotropic non-Newtonian fluid problems can be addressed. The project supported by the National Natural Science Foundation of China (10372100, 19832050) (Key project). The English text was polished by Yunming Chen.     

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.     

An Analysis of the Movement of Wetting and Nonwetting Fluids in Homogeneous Porous Media

The movement of wetting and nonwetting fluid flow in columns packed with glass beads is used to understand the more complicated flows in homogeneous porous media. The motion of two immiscible liquids (oil and water) is observed with different surfactants. Through dimensional analyses, fluid velocity is well correlated with interfacial tension and less dependent on particle size. In water–oil (W/O) experiments, finger pattern flows are observed if water is the displacing fluid that flows in an oil-filled porous media, whereas oil ganglia tend to form if oil is the displacing fluid in the water-wetted porous media. The results are well described by a simple model based on an earlier theory of flow in a tube.     

Effect of anisotropic resistance characteristic on boundary-layer transitional flow

It is observed that the feather surface exhibits anisotropic resistances for the streamwise and spanwise flows. To obtain a qualitative understanding about the effect of this anisotropic resistance feature of surface on the boundary-layer transitional flow over a flat plate, a simple phenomenological model for the anisotropic resistance is established in this paper. By means of the large eddy simulation (LES) with high-order accurate finite difference method, the numerical investigations are conducted. The numerical results show that with the spanwise resistance hindering the formation of vortexes, the transition from laminar flow to turbulent flow can be delayed, and turbulence is weakened when the flow becomes fully turbulent, which leads to significant drag reduction for the plate. On the contrary, the streamwise resistance renders the flow less stable, which leads to the earlier transition and enhances turbulence in the turbulent region, causing a drag increase for the plate. Thus, it is indicated that a surface with large resistance for spanwise flow and small resistance for streamwise flow can achieve significant drag reduction. The present results highlight the anisotropic resistance characteristic near the feather surface for drag reduction, and shed a light on the study of bird’s efficient flight.     

A streaming birefringence study of the flow at the junction of the aorta and the renal arteries

Streaming birefringence with an organic dye (Milling Yellow) was used to investigate the flow near the junction of the renal arteries and the descending aorta in a model of human vessels. The dye concentration was adjusted to give fluid rheological properties, typical of blood. Steady and pulsatile flow were investigated at branch-to-trunk flow ratios of 0.050–0.350. The flow ratio range over which flow separation and simple secondary flows were identified during systole near the renal ostia are reported. Streaming birefringence has the advantage of allowing visualization of the entire flow field. Also, the fluid rather than suspended particles are observed. An important disadvantage, however, is that three-dimensional flows make interpretation difficult.     

The relevance of entry flow measurements for the estimation of extensional viscosity of polymer melts

The extensional viscosity of some flexible chain polymers and a thermotropic liquid crystalline polymer was measured in uniaxial extensional flow at constant extension rate. Power law functions were found for the dependence of the extensional viscosity at constant accumulated strain on strain rate. The stress growth curves were compared with measurements in axisymmetric entry flow, where both elongation and shear occur. The comparison showed that the values of the extensional viscosity calculated from the measurements in the entry flow correspond to the ones calculated from the viscosity growth measured in uniaxial elongation and averaged over extensional strain equal to what is accumulated on the fluid as it flows from the barrel into the capillary.     

Boundary element analysis of three-dimensional mixing flow of Newtonian and viscoelastic fluids

The boundary element method (BEM) is implemented for the simulation of three-dimensional transient flows of typical relevance to mixing. Creeping Newtonian and viscoelastic fluids of the Maxwell type are examined. A boundary-only formulation in the time domain is proposed for linear viscoelastic flows. Special emphasis is placed on cavity flows involving simple- and multiple-connected moving domains. The BEM becomes particularly suited in multiple-connected flows, where part of the boundary (stirrer or rotor) is moving, and the remaining outer part (cavity or barrel) is at rest. In this case, conventional methods, such as the finite element method (FEM), generally require remeshing or mesh refinement of the three-dimensional fluid volume as the flow evolves and the domain of computation changes with time. The BEM is shown to be much easier to implement since the kinematics of the elements bounding the fluid is known (imposed). It is found that, for simple cavity flow induced by a rotating vane at constant angular velocity, the tractions at the vane tip and cavity face exhibit non-linear periodic dynamical behavior with time for fluids obeying linear constitutive equations. © 1998 John Wiley & Sons, Ltd.     

On selection of reference temperature of heat transfer coefficient for complicated flows

Convective heat transfer coefficient is closely related with flow and thermal conditions. To define heat transfer coefficient, a reference temperature needs to be properly selected, which can be the fluid bulk mean temperature for internal flows or the temperature at the far field for external flows. For complicated flows, the adiabatic wall temperature is commonly adopted as the reference temperature, while other options can also be applied. This paper analyzed some of the potential selections of the reference temperature for different flow settings, including film cooling, jet impingement with cross flows, and a mixing flow in a straight duct with or without internal heat source. It is observed that heat transfer coefficient changes dramatically with selection of reference temperatures. In case of constant wall temperature, using adiabatic wall temperature as reference temperature can result in negative heat transfer coefficient, which means the heat flux has a different direction with the defined driving temperature difference. To avoid the inconsistency due to the reference temperature, an innovative method is proposed to calculate the heat transfer coefficient of complicated flows.     

Simple,coupled algorithms for solving creeping flows and their application to electro‐osmotic flows

In this study we developed simple, coupled algorithms for solving low‐Reynolds‐number flows applicable to micro‐scale flows such as electro‐osmotic flows. The most popular scheme, i.e. the projection method, is not suitable for such flows because of its undesirable slip effect on boundaries at low‐Reynolds‐numbers. In our method, the velocity and pressure are strongly coupled, and the momentum and pressure equations are solved iteratively by using the successive over relaxation (SOR) method while exchanging the unknown variables as soon as they have been updated. The developed methods are applied to a model flow for evaluating their performance. It was found that the coupled schemes are indeed superior to a projection method, i.e. the fractional‐step method, in both numerical accuracy and CPU time. The code is then applied to a dc electro‐osmotic flow within a cavity driven by electrical force acting on the ions spread in the fluid. In this application, the system of equations for the fluid flow and that for the ion transport are solved in a decoupled way, but each system is solved by using fully implicit schemes. From the simulations and by introducing the concept of vorticity source, we can identify two roles of the body force, one contributing to build‐up of the osmotic pressure and the other to the fluid flow. The interesting reverse flow occurring after the external potentials applied on the electrodes have been shut off is also investigated in terms of the vorticity source. Copyright © 2010 John Wiley & Sons, Ltd.     

Vortex Formations in Plane Flows of Viscous and Ideal Fluids

An analytic representation of individual vortex formations with different cross-sections is obtained for plane-parallel flows of viscous and ideal fluids. A flow in which one of the streamlines has cusp is also found. In these flows the vorticity is constant; therefore, the solutions obtained for the equations describing the kinematic variables in ideal fluid flows must simultaneously satisfy the Navier-Stokes equations.     

Correcting Anisotropic Gradient Transport of k

The gradient transport model for k is extended to classes of turbulent flows for which the gradient transport hypothesis is relevant but the anisotropy of the Reynolds stress, to which the eddy diffusivity is proportional, is large and variable. In highly anisotropic turbulence the standard isotropic model used in engineering practice is fundamentally wrong and the conventional anisotropic approximation inadequate. The work is motivated by the important observations that the eddy diffusivity coefficient for a standard gradient transport model for various transported quantities is a factor of 3–10 times larger in highly anisotropic turbulence than that used in standard engineering models. While the conventional anisotropic eddy diffusivity approximation appears adequate for material conserved scalars it is inadequate for k. The problem is solved by addressing the anisotropy of the turbulent transport of k at the level of the underlying third order tensor. It is shown that, unlike the traditional transport models for k, the orientation of the anisotropy with respect to the direction of the gradient plays a crucial role not accounted for in conventional models used in engineering calculations. The new anisotropic eddy diffusivity tensor is quadratic in the anisotropy (the traditional model is linear in the anisotropy). It is shown that the new more rigorous anisotropic eddy diffusivity varies 300% more than the standard model comparing the isotropic limit to the value for the two-dimensional limit. The two-dimensional limit is important in strongly stably stratified flows, in pressure gradient or shock driven flows and in rotating flows. Using the simple shear and the homogeneous non-equilibrium Rayleigh Taylor turbulence the new anisotropic diffusivity tensor is validated in inhomogeneous Rayleigh Taylor turbulence at early and late times.