Application of the log-conformation tensor to three-dimensional time-dependent free surface flows

The numerical simulation of flows of highly elastic fluids has been the subject of intense research over the past decades with important industrial applications. Therefore, many efforts have been made to improve the convergence capabilities of the numerical methods employed to simulate viscoelastic fluid flows. An important contribution for the solution of the High-Weissenberg Number Problem has been presented by Fattal and Kupferman [J. Non-Newton. Fluid. Mech. 123 (2004) 281-285] who developed the matrix-logarithm of the conformation tensor technique, henceforth called log-conformation tensor. Its advantage is a better approximation of the large growth of the stress tensor that occur in some regions of the flow and it is doubly beneficial in that it ensures physically correct stress fields, allowing converged computations at high Weissenberg number flows. In this work we investigate the application of the log-conformation tensor to three-dimensional unsteady free surface flows. The log-conformation tensor formulation was applied to solve the Upper-Convected Maxwell (UCM) constitutive equation while the momentum equation was solved using a finite difference Marker-and-Cell type method. The resulting developed code is validated by comparing the log-conformation results with the analytic solution for fully developed pipe flows. To illustrate the stability of the log-conformation tensor approach in solving three-dimensional free surface flows, results from the simulation of the extrudate swell and jet buckling phenomena of UCM fluids at high Weissenberg numbers are presented. (C) 2012 Elsevier B.V. All rights reserved.

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.

Capillary break-up rheometry of low-viscosity elastic fluids

We investigate the dynamics of the capillary thinning and break-up process for low viscosity elastic fluids such as dilute polymer solutions. Standard measurements of the evolution of the midpoint diameter of the necking fluid filament are augmented by high speed digital video images of the break up dynamics. We show that the successful operation of a capillary thinning device is governed by three important time scales (which characterize the relative importance of inertial, viscous and elastic processes), and also by two important length scales (which specify the initial sample size and the total stretch imposed on the sample). By optimizing the ranges of these geometric parameters, we are able to measure characteristic time scales for tensile stress growth as small as 1 millisecond for a number of model dilute and semi-dilute solutions of polyethylene oxide (PEO) in water and glycerol. If the final aspect ratio of the sample is too small, or the total axial stretch is too great, measurements are limited, respectively, by inertial oscillations of the liquid bridge or by the development of the well-known beads-on-a-string morphology which disrupt the formation of a uniform necking filament. By considering the magnitudes of the natural time scales associated with viscous flow, elastic stress growth and inertial oscillations it is possible to construct an operability diagram characterizing successful operation of a capillary break-up extensional rheometer. For Newtonian fluids, viscosities greater than approximately 70 mPas are required; however for dilute solutions of high molecular weight polymer, the minimum Viscosity is substantially lower due to the additional elastic stresses arising from molecular extension. For PEO of molecular weight 2.10(6) g/mol, it is possible to measure relaxation times of order 1 ms in dilute polymer solutions with zero-shear-rate viscosities on the order of 2-10 mPas.

The atomisation of detergent slurries

A study was made to determine the conditions under which the optimum droplet size distribution (ie., narrowest size range with a minimum of fines and over-sized agglomerates), is generated in sprays from centrifugal pressure nozzles. A range of non-Newtonian detergent slurries were tested but the results are of wider application and parallel work was undertaken with water, ionic solutions and chalk slurries. Six centrifugal pressure nozzles were used and the drop-size distributions correlated as a function of fluid properties, pressure, fiowrate, feed temperature, and nozzle characteristics. Measurements were made using a Malvern Particle Size Anayser slung across a specially-designed transparent tower section of approximately 1.2m diameter in order to reduce obscuration caused by the spray and improve existing droplet sizing techniques. The results obtained were based upon the Rosin-Rammler distribution model and the Size Analyser provided a direct print-out of size distribution and the parameters characterising it. A Spraying System nozzle, AAASSTC8-8, produced the optimum spray distribution with the detergent slurry at a temperature of 60°C whilst operating at 1200 psi. With other fluids the Delevan 2.2SJ nozzle produced the optimum spray distribution operating at 1200 psi but with the Spraying Systems nozzles there was no clear-cut optimum set of conditions, ie. the nozzle and pressure varied depending upon the fluid being sprayed. The mechanisms of liquid sheet break-up and droplet dispersion were investigated in specially-constructed, scaled-up, transparent nozzles. A mathematical model of centrifugal pressure nozzle atomisation was developed based upon fundamental operating parameters rather than resorting to empirical correlations. This enabled theoretical predictions to be made over a wide range of operating conditions and nozzle types. The model predictions for volumetric fiowrate, liquid sheet length and air core diameter showed good agreement with the experimentally determined results. However, the model predicted smaller droplet sizes than were produced experimentally due to inaccuracies identified in the initial assumptions.

Drag coefficient in a crystal suspension

A single air bubble rising in xanthan gum crystal suspension has been studied experimentally. The suspension was made by different concentrations of xanthan gum solutions with 0.23 mm polystyrene crystal particles. Drag co-efficient data and a new correlation of drag coefficient is presented for spherical and nonspherical bubbles in non-Newtonian crystal suspension. The correlation is developed in terms of the Reynolds number, Re and the bubble shape factor, � (the ratio between the surface equivalent sphere diameter to the volume equivalent sphere diameter). The experimental drag coefficient was found to be consistent with this new predicted correlation and published data over the ranges, 0.1

Suspended sediment properties and suspended sediment flux estimates in an inundated urban environment during a major flood event

During a major flood event, the inundation of urban environments leads to some complicated flow motion most often associated with significant sediment fluxes. In the present study, a series of field measurements were conducted in an inundated section of the City of Brisbane (Australia) about the peak of a major flood in January 2011. Some experiments were performed to use ADV backscatter amplitude as a surrogate estimate of the suspended sediment concentration (SSC) during the flood event. The flood water deposit samples were predominantly silty material with a median particle size about 25 μm and they exhibited a non-Newtonian behavior under rheological testing. In the inundated urban environment during the flood, estimates of suspended sediment concentration presented a general trend with increasing SSC for decreasing water depth. The suspended sediment flux data showed some substantial sediment flux amplitudes consistent with the murky appearance of floodwaters. Altogether the results highlighted the large suspended sediment loads and fluctuations in the inundated urban setting associated possibly with a non-Newtonian behavior. During the receding flood, some unusual long-period oscillations were observed (periods about 18 min), although the cause of these oscillations remains unknown. The field deployment was conducted in challenging conditions highlighting a number of practical issues during a natural disaster.

Analytical Solution of a Walters’ Liquid B Flow Over a Linear Stretching Sheet in a Porous Medium

This chapter represents the analytical solution of two-dimensional linear stretching sheet problem involving a non-Newtonian liquid and suction by (a) invoking the boundary layer approximation and (b) using this result to solve the stretching sheet problem without using boundary layer approximation. The basic boundary layer equations for momentum, which are non-linear partial differential equations, are converted into non-linear ordinary differential equations by means of similarity transformation. The results reveal a new analytical procedure for solving the boundary layer equations arising in a linear stretching sheet problem involving a non-Newtonian liquid (Walters’ liquid B). The present study throws light on the analytical solution of a class of boundary layer equations arising in the stretching sheet problem.

Analytical solution of unsteady three-dimensional MHD boundary layer flow and heat transfer due to impulsively stretched plane surface

The unsteady magnetohydrodynamic viscous flow and heat transfer of Newtonian fluids induced by an impulsively stretched plane surface in two lateral directions are studied by using an analytic technique, namely, the homotopy method. The analytic series solution presented here is highly accurate and uniformly valid for all time in the entire region. The effects of the stretching ratio and the magnetic field on the surface shear stresses and heat transfer are studied. The surface shear stresses in x- and y-directions and the surface heat transfer are enchanced by increasing stretching ratio for a fixed value of the magnetic parameter. For a fixed stretching ratio, the surface shear stresses increase with the magnetic parameter, but the heat transfer decreases. The Nusselt number takes longer time to reach the steady state than the skin friction coefficients. There is a smooth transition from the initial unsteady state to the steady state.

Flow of micropolar fluid through a tube with stenosis

mathematical model for the steady flow of non-Newtonian fluid through a stenotic region is presented. The results indicate that the general shape and size of the stenosis together with rheological properties of blood are important in understanding the flow characteristics and the presence of flow separation.

Free convection heat transfer to mercury in vertical annuli

Data on free convection heat transfer to water and mercury are collected using a test rig in vertical annuli of three radii ratios, the walls of which are maintained at uniform temperatures. A theoretical analysis of the boundary layer equations has been attempted using local similarity transformation and double boundary layer approach. Correlations derived from the present theoretical analysis are compared with the analysis and the experimental data available in literature for non-metallic fluids and also with the present experimental data on water and mercury. Generalised correlations are set up for expressing the ratio of heat transferred by convection to the heat transferred by pure conduction and Nusselt's number, in terms of Grashof, Rayleigh and Prandtl numbers, based on the theoretical analysis and the present data on mercury and water. The present generalised correlations agree with the reported and present data for non-metallic fluids and liquid metals with an average deviation of 9% and maximum deviation of ± 13.7%.

Generalized aerodynamic noise equation

An exact aerodynamic noise equation is formulated for Newtonian fluids. The cause−effect problem is discussed. Finally, the importance of external additions of mass, momentum, and energy is examined. Physics of Fluids is copyrighted by The American Institute of Physics.

Unsteady mixed convection flow of a thermomicropolar fluid on a long thin vertical cylinder

The unsteady laminar mixed convection boundary layer flow of a thermomicropolar fluid over a long thin vertical cylinder has been studied when the free stream velocity varies with time. The coupled nonlinear partial differential equations with three independent variables governing the flow have been solved numerically using an implicit finite difference scheme in combination with the quasilinearization technique. The results show that the buoyancy, curvature and suction parameters, in general, enhance the skin friction, heat transfer and gradient of microrotation, but the effect of injection is just opposite. The skin friction and heat transfer for the micropolar fluid are considerably less than those for the Newtonian fluids. The effect of microrotation parameter is appreciable only on the microrotation gradient. The effect of the Prandtl number is appreciable on the skin friction, heat transfer and gradient of microtation.

Steady flow of couple stress fluid through tubes of slowly varying cross-sections--application to blood flows

Steady laminar flow of a non-Newtonian fluid based on couple stress fluid theory, through narrow tubes of varying cross-sections has been studied theoretically. Asymptotic solutions are obtained for the basic equations and the expressions for the velocity field and the wall shear stress are derived for a general cross-section. Computation and discussions are carried out for the geometries which occur in the context of physiological flows or in particular blood flows. The tapered tubes and constricted tubes are of special importance. It is observed that increase in certain parameters results in erratic flow behaviour proximal to the constricted areas which is further enhanced by the increase in the geometric parameters. This elucidates the implications of the flow in the development of vascular lesions.

Sodium-alginate-based proton-exchange membranes as electrolytes for DMFCs

Novel mixed-matrix membranes prepared by blending sodium alginate (NaAlg) with polyvinyl alcohol (PVA) and certain heteropolyacids (HPAs), such as phosphomolybdic acid (PMoA), phosphotungstic acid (PWA) and silicotungstic acid (SWA), followed by ex-situ cross-linking with glutaraldehyde (GA) to achieve the desired mechanical and chemical stability, are reported for use as electrolytes in direct methanol fuel cells (DMFCs). NaAlg-PVA-HPA mixed matrices possess a polymeric network with micro-domains that restrict methanol cross-over. The mixed-matrix membranes are characterised for their mechanical and thermal properties. Methanol cross-over rates across NaAlg-PVA and NaAlg-PVA-HPA mixed-matrix membranes are studied by measuring the mass balance of methanol using a density meter. The DMFC using NaAlg-PVA-SWA exhibits a peak power-density of 68 mW cm(-2) at a load current-density of 225 mA cm(-2), while operating at 343 K. The rheological properties of NaAlg and NaAlg-PVA-SWA viscous solutions are studied and their behaviour validated by a non-Newtonian power-law.

Spray characterization of straight vegetable oils at high injection pressures

We present results of high pressure spray characterization of Straight Vegetable Oils (SVOs) which are potential diesel fuel substitutes. SVO sprays are visualized at high injection pressures (up to 1600 bar) to study their atomization characteristics. Spray structure studies are reported for the first time for Jatropha and Pongamia vegetable oils, under atmospheric conditions. Jatropha and Pongamia SVO sprays are found to be poorly atomized and intact liquid cores are observed even at an injection pressure of 1600 bar. Non-Newtonian behavior of Jatropha and Pongamia oil is shown to be the reason for observed spray structure. (C) 2012 Elsevier Ltd. All rights reserved.