Rheological characterization of complex fluids in electro-magnetic fields

The paper is focused on the experimental investigations and rheological characterisations in magnetic and electric fields of liquids based on water in crude oils emulsions, added with ferrofluids (two types of crude oils are used in experiments: asphaltic and paraffinic, respectively). The final samples disclose weakly effects in the presence of magnetic field (saturated magnetization: M_n < 300 [G]) and behave almost as isolators in electric field (conductivity: σ < 10⁻⁵ [S/m]). The main goal of the study is to explore to what extent rheometry of complex fluids in electric and magnetic fields is able to offer value information about the internal structure of the samples. The experimental results prove that anomalous rheological behaviour (thixotropy, non-monotonic flow curve or viscosity function) of a complex fluid (in our case, emulsions based on paraffinic oil) generate also thixotropic properties and non-monotonic answers in the presence ferrofluids, under low magnetic and/or electric fields intensity. Our prospective study suggests that novel experimental procedures based on interaction: electro-magnetic field–complex fluids can be developed, in order to determine indirectly some relevant rheological properties of the complex fluids with internal network structure.     

Validation of a discrete element model using magnetic resonance measurements

The discrete element model （DEM） is a very promising modelling strategy for two-phase granular systems. However, owing to a lack of experimental measurements, validation of numerical simulations of two-phase granular systems is still an important issue. In this study, a small two-dimensional gas- fluidized bed was simulated using a discrete element model. The dimensions of the simulated bed were 44mm × 10mm × 120 mm and the fluidized particles had a diameter dp = 1.2 mm and density ρp = 1000 kg/m^3. The comparison between DEM simulations and experiments are performed on the basis of time-averaged voidage maps. The drag-law of Beetstra et al. [Beetstra, R., van der Hoef, M.A., ＆ Kuipers,J. A. M. （2007b）. Drag force of intermediate Reynolds number flow past mono- and bidispersed arrays of spheres. AIChE Journal, 53,489-501 ] seems to give the best results. The simulations are fairly insensitive to the coefficient of restitution and the coefficient of friction as long as some route of energy dissipation during particle-particle and particle-wall contact is provided. Changing the boundary condition of the gas phase at the side-walls from zero-slip to full-slip does not affect the simulation results. Care is to be taken that the cell sizes are chosen so that a reasonable number of particles can be found in a fluid cell.     

Numerical simulations of negatively buoyant jets in an immiscible fluid using the Particle Finite Element Method

Negatively buoyant jets consist in a dense fluid injected vertically upward into a lighter ambient fluid. The numerical simulation of this kind of buoyancy‐driven flows is challenging as it involves multiple fluids with different physical properties. In the case of immiscible fluids, it requires, in addition, to track the motion of the interface between fluids and accurately represent the discontinuities of the flow variables. In this paper, we investigate numerically the injection of a negatively buoyant jet into a homogenous immiscible ambient fluid using the Particle Finite Element Method and compare the two‐dimensional numerical results with experiments on the injection of a jet of dyed water through a nozzle in the base of a cylindrical tank containing rapeseed oil. In both simulations and experiments, the fountain inlet flow velocity and nozzle diameter have been varied to cover a wide range of Froude Fr and Reynolds Re numbers ( 0.1 < Fr < 30, 8 < Re < 1350), reproducing both weak and strong laminar fountains. The flow behaviors observed for the different numerical simulations fit in the regime map based on the Re and Fr values of the experiments, and the maximum fountain height is in good agreement with the experimental observations, suggesting that particle finite element method is a useful tool for the study of immiscible two‐fluid systems. Copyright © 2011 John Wiley & Sons, Ltd.     

Series solution of the temperature distribution in the Falkner–Skan wedge flow by the homotopy analysis method

A new analytical method, namely the homotopy analysis method (HAM), has been applied to investigate the temperature field associated with the Falkner–Skan boundary-layer problem, and a series solution is provided in this paper. The results of the present work show agreement with those of numerical solutions in a large range of Prandtl numbers (0 < Pr ≤ 100), which demonstrates the validity of the present analysis.     

Cohesive properties of nickel-alumina interfaces determined via simulation of ductile bridging experiments

A combined experimental and computational study is carried out to characterize a nickel-alumina interface in terms of the two parameter (σ̂, Γ₀) computational cohesive zone (CCZ) model of Tvergaard and Hutchinson. Experiments were performed using a sandwich specimen consisting of a thin nickel foil bonded between two pre-cracked alumina plates. The specimen was loaded in tension with the nickel foil bridging the cracks in the ceramic. Numerical simulations of the experiments were used to extract the parameters for the CCZ model.Effects of various parameters of the CCZ model are investigated and it is found that the most dominant parameter is the interface strength, σ̂. Effects of the residual thermal stresses are also investigated and it is shown that these stresses can enhance the specimen fracture toughness by almost 16%. The parameters for the nickel-alumina interface are found to be σ̂ = 148 MPa and Γ₀ = 11 J m⁻². It is observed that for the foil thicknesses tested, the work of rupture does not vary linearly with the thickness as predicted by many theoretical models. We found that interfaces which are neither too strong nor too weak contribute most to the overall fracture toughness of such a composite. Although the macroscopic loading at the nickel-alumina interface is shear, the failure is primarily tensile due to the thinning that occurs in the metal as it is stretched.     

Heat transfer to non-Newtonian flows over a cylinder in cross flow

Over a range of 102<Re^*<5800, 6.5<Pr^*<79, and 0.6<n<1, circumferential wall temperatures for water and aqueous polymer (purely viscous) solution flows over a smooth cylinder were measured experimentally. The cylinder was heated by passing direct electric current through it. Aqueous solutions of Carbopol 934 and EZ1 were used as power-law non-Newtonian fluids. The peripherally averaged heat transfer coefficient for purely viscous non-Newtonian fluids, at any fixed flow rate, decreases with increasing polymer concentration. A new correlation is proposed for predicting the peripherally averaged Nusselt number for power-law fluid flows over a heated cylinder in cross flow.     

A hybrid Lagrangian–Eulerian particle‐level set method for numerical simulations of two‐fluid turbulent flows

A coupled Lagrangian interface‐tracking and Eulerian level set (LS) method is developed and implemented for numerical simulations of two‐fluid flows. In this method, the interface is identified based on the locations of notional particles and the geometrical information concerning the interface and fluid properties, such as density and viscosity, are obtained from the LS function. The LS function maintains a signed distance function without an auxiliary equation via the particle‐based Lagrangian re‐initialization technique. To assess the new hybrid method, numerical simulations of several ‘standard interface‐moving’ problems and two‐fluid laminar and turbulent flows are conducted. The numerical results are evaluated by monitoring the mass conservation, the turbulence energy spectral density function and the consistency between Eulerian and Lagrangian components. The results of our analysis indicate that the hybrid particle‐level set method can handle interfaces with complex shape change, and can accurately predict the interface values without any significant (unphysical) mass loss or gain, even in a turbulent flow. The results obtained for isotropic turbulence by the new particle‐level set method are validated by comparison with those obtained by the ‘zero Mach number’, variable‐density method. For the cases with small thermal/mass diffusivity, both methods are found to generate similar results. Analysis of the vorticity and energy equations indicates that the destabilization effect of turbulence and the stability effect of surface tension on the interface motion are strongly dependent on the density and viscosity ratios of the fluids. Copyright © 2007 John Wiley & Sons, Ltd.     

Experimental investigations and numerical simulations for an open channel flow of a weak elastic polymer solution around a T-profile

The present paper is concerned with experimental and numerical investigations of planar complex flows of weak elastic polymer solutions (whose concentration are below the critical overlap concentration), characterised by small relaxation times (<0.1 s) and almost constant shear viscosities for small and medium shear rates. The main aim of the study is to detect to what extent a very small amount of elasticity present in a viscous fluid can influence its behaviour in complex flows, without introducing major modifications of classical rheological tests. The samples are polymer solutions of low PIB molecular weight dissolved in highly viscous Newtonian mineral oil. The analysed motion is steady, and takes place in an open channel around a T profile. Maximum values of the characteristic parameters for the experiments, the Reynolds and Weissenberg numbers, were 45 and 0.1, respectively. The experiments show a decrease of the wake length downstream the profile for weak elastic solutions in comparison to the Newtonian solvent. Actually, the same wake length as in the Newtonian case was obtained for tested polymer solutions, but at higher Re numbers. Numerical simulations using the Giesekus model predict the same behaviour and are consistent with experiments from both qualitative and quantitative point of views. The results of research conclude that, even in small amounts, the presence of elasticity in pure viscous liquids induces quantitative changes from Newtonian flow in complex dominant elongational flows, at elongational rates for which the sudden thickening of extensional viscosity is remarkable. The study is important, since it should enable better understanding and modelling of viscoelastic flows that involve dilute polymer solutions, or fluids with similar rheology; biofluid mechanics being one area of application of this research. Corroboration of experimental flow visualization with numerical simulation is currently a feasible method used to characterise weak elastic polymer solutions, since classical rheological techniques generally fail to obtain realistic values of relaxation time for these particular viscoelastic fluids. Corneliu Balan dedicates this paper to the anniversary of one hundred years from the birth of Academician Dumitru Dumitrescu (1904–1983), charismatic personality of the Romanian school of fluid mechanics.

---

On corner flows of Oldroyd-B fluids

Exact series solutions for planar creeping flows of Oldroyd-B fluids in the neighbourhood of sharp corners are presented and discussed. Both reentrant and non-reentrant sectors are considered. For reentrant sectors it is shown that more than one type of series solution can exist formally, one type exhibiting Newtonian-like asymptotic behaviour at the corner, away from walls, and another type exhibiting the same kind of asymptotics as an Upper Convected Maxwell (UCM) fluid. The solutions which are Newtonian-like away from walls are shown to develop non-integrable stress singularities at the walls when the no-slip velocity boundary condition is imposed. These mathematical solutions are therefore inadmissible from the physical viewpoint under no-slip conditions. An inadmissible solution, with stress singularities which are not everywhere integrable, is identified among the solutions of UCM-type. For a 270° reentrant sector the radial behaviour of the normal stress is everywhere r^−0.613. In the viscometric region near a wall, the radial normal stress σ^rr behaves like (rε)^−0.613, where ε is the angle made with the wall. In addition σ^rθ is infinite (not integrable) at the wall even when r is non-zero. Another UCM-type solution has a normal stress behaviour away from walls which is r^−0.985 for 270° sector. Again, this solution has a non-integrable stress singularity and is therefore inadmissible. Finally, for non-reentrant sectors it is shown that the flow is always Newtonian-like away from walls.     

Dynamics of nanofibres conveyed by low Reynolds number flow in a microchannel

In this paper we aim to create an experimental and numerical model of nano and micro filaments suspended in a confined Poiseuille flow. The experimental data obtained for short nanofibres will help to elucidate fundamental questions concerning mobility and deformation of biological macromolecules due to hydrodynamic stresses from the surrounding fluid motion. Nanofibres used in the experiments are obtained by electrospinning polymer solutions. Their typical dimensions are 100–1000 μm (length) and 0.1–1 μm (diameter). The nanofibre dynamics is followed experimentally under a fluorescence microscope. A precise multipole expansion method of solving the Stokes equations, and its numerical implementation are used to construct a bead-spring model of a filament moving in a Poiseuille flow between two infinite parallel walls. Simulations show typical behaviour of elongated macromolecules. Depending on the parameters, folding and unfolding sequences of a flexible filament are observed, or a rotational and translation motion of a shape-preserving filament. An important result of our experiments is that nanofibres do not significantly change their shape while interacting with a micro-flow. It appeared that their rotational motion is better reproduced by the shape-preserving Stokesian bead model with all pairs of beads connected by springs, omitting explicit bending forces.     

Numerical results for the conservation of energy of Maxwell media for the Rayleigh–Stokes problem

This publication continues our studies of analytical solutions of the Rayleigh–Stokes problem for Maxwell fluids [J. Zierep, C. Fetecau, Energetic balance for the Rayleigh–Stokes problem of a Maxwell fluid, Int. J. Eng. Sci. 45 (2007) 617–627]. We start from the Fourier sine transform. The numerical result is given and discussed for the velocity u, the power of the wall shear stresses L, the dissipation Φand the boundary layer thickness δ. These new results are important for nature and technology.     

Three-dimensional flow of Newtonian and Boger fluids in square–square contractions

The flow of a Newtonian fluid and a Boger fluid through sudden square–square contractions was investigated experimentally aiming to characterize the flow and provide quantitative data for benchmarking in a complex three-dimensional flow. Visualizations of the flow patterns were undertaken using streak-line photography, detailed velocity field measurements were conducted using particle image velocimetry (PIV) and pressure drop measurements were performed in various geometries with different contraction ratios. For the Newtonian fluid, the experimental results are compared with numerical simulations performed using a finite volume method, and excellent agreement is found for the range of Reynolds number tested (Re₂ ≤ 23). For the viscoelastic case, recirculations are still present upstream of the contraction but we also observe other complex flow patterns that are dependent on contraction ratio (CR) and Deborah number (De₂) for the range of conditions studied: CR = 2.4, 4, 8, 12 and De₂ ≤ 150. For low contraction ratios strong divergent flow is observed upstream of the contraction, whereas for high contraction ratios there is no upstream divergent flow, except in the vicinity of the re-entrant corner where a localized atypical divergent flow is observed. For all contraction ratios studied, at sufficiently high Deborah numbers, strong elastic vortex enhancement upstream of the contraction is observed, which leads to the onset of a periodic complex flow at higher flow rates. The vortices observed under steady flow are not closed, and fluid elasticity was found to modify the flow direction within the recirculations as compared to that found for Newtonian fluids. The entry pressure drop, quantified using a Couette correction, was found to increase with the Deborah number for the higher contraction ratios.     

Dynamic interaction of a spherical radiator in a fluid-filled cylindrical borehole within a poroelastic formation

Harmonic acoustic radiation from a modally oscillating spherical source positioned at the center of a fluid-filled cylindrical cavity embedded within a fluid-saturated porous elastic formation is studied in an exact manner. The formulation utilizes the Biot theory of dynamic poroelasticity along with the cylindrical to spherical wave-field transformations, and the pertinent boundary conditions to obtain a closed-form series solution. The analytical results are illustrated with a numerical example in which the spherical source, with its polar axis oriented along the main axis of a water-filled borehole and embedded within a water-saturated Ridgefield sandstone formation, is excited in vibrational modes of various orders. The magnitude of the reflected component of acoustic pressure along the axis of the borehole for a pulsating (n = 0), an oscillating (n = 1), and also a multipole (n = 0–3) spherical source as a function of the excitation frequency is calculated and discussed for representative values of the parameters characterizing the system. Special attention is paid to the effects of source excitation frequency, size, surface velocity profile, and internal impedance as well as borehole interface permeability condition on the reflected pressure magnitudes. Limiting cases are considered and fair agreements with well-known solutions are obtained.     

Three-dimensional numerical simulation of thermocapillary flow in a shallow cylindrical pool

In order to understand the characteristics of thermocapillary flow, we conducted a series of unsteady three-dimensional (3D) numerical simulations of thermocapillary flow for fluids with Prandtl number Pr = 0.011 and 1.0 in a shallow cylindrical pool with an azimuthal temperature non-uniform, an adiabatic solid bottom and free surface. The simulation results indicate that thermocapillary flow is steady 3D flow at the small temperature non-uniform. When temperature non-uniform exceeds some critical value, the flow will undergo a transition to 3D oscillatory flow, which is characterized by traveling oscillatory wave reversely along the flow direction on the free surface and diffusing to two sides at the same time. Further, the critical conditions for the flow transition are determined.     

Nonlinear resonances and self-excited oscillations of a rotor caused by radial clearance and collision

This paper investigates oscillations in a flexible rotor system with radial clearance between an outer ring of the bearing and a casing by experiments and numerical simulations. The mathematical model considers the collisions of the bearing with the casing. The following phenomena are found: (1) Nonlinear resonances of subharmonic, super-subharmonic and combination oscillation occur. (2) Self-excited oscillation of a forward whirling mode occurs in a wide range above the major critical speed. (3) Entrainment phenomena from self-excited oscillation to nonlinear forced oscillation occur at these nonlinear resonance ranges. Moreover, this study analyzes periodic solutions of the mathematical model by the Harmonic Balance Method (HBM). As the results, the nonlinear resonances of subharmonic oscillation and its entrainment phenomenon can be explained theoretically by investigating the stability of the periodic solutions. The influence of the static force and the bearing damping on these oscillation are also clarified.     

Numerical simulation on micromixer based on synthetic jet

This paper studied a concept of micromixer with a synthetic jet placed at the bottom of a rectangular channel. Due to periodic ejections from and suctions into the channel, the fluids are mixed effectively. To study the effects of the inlet velocity, the jet intensity and frequency, and the jet location on the mixing efficiency, 3-D numerical simulations of the micromixer have been carried out. It has been found that when the jet intensity and the frequency are fixed, the mixing efficiency increases when Re 〈 50, and decreases when Re 〉 50 with the best mixing efficiency achieved at Re = 50. When the ratio of the jet velocity magnitude to the inlet velocity is taken as 10 and the jet frequency is 100 Hz, the mixing index reaches the highest value. It has also been found that to get better mixing efficiency, the orifice of the synthetic jet should be asymmetrically located away from the channel＇s centerline.     

A numerical investigation of the effects of the spanwise length on the 3-D wake of a circular cylinder

The numerical prediction of vortex-induced vibrations has been the focus of numerous investigations to date using tools such as computational fluid dynamics. In particular, the flow around a circular cylinder has raised much attention as it is present in critical engineering problems such as marine cables or risers. Limitations due to the computational cost imposed by the solution of a large number of equations have resulted in the study of mostly 2-D flows with only a few exceptions. The discrepancies found between experimental data and 2-D numerical simulations suggested that 3-D instabilities occurred in the wake of the cylinder that affect substantially the characteristics of the flow. The few 3-D numerical solutions available in the literature confirmed such a hypothesis. In the present investigation the effect of the spanwise extension of the solution domain on the 3-D wake of a circular cylinder is investigated for various Reynolds numbers between 40 and 1000. By assessing the minimum spanwise extension required to predict accurately the flow around a circular cylinder, the infinitely long cylinder is reduced to a finite length cylinder, thus making numerical solution an effective way of investigating flows around circular cylinders. Results are presented for three different spanwise extensions, namely πD/2, πD and 2πD. The analysis of the force coefficients obtained for the various Reynolds numbers together with a visualization of the three-dimensionalities in the wake of the cylinder allowed for a comparison between the effects of the three spanwise extensions. Furthermore, by showing the different modes of vortex shedding present in the wake and by analysing the streamwise components of the vorticity, it was possible to estimate the spanwise wavelengths at the various Reynolds numbers and to demonstrate that a finite spanwise extension is sufficient to accurately predict the flow past an infinitely long circular cylinder.     

Stokes paradox for power-law flow around a cylinder

A simple analysis for power-law fluids shows that the Stokes paradox for creeping flow around a cylinder is removed for shear-thinning (n < 1) but not for shear-thickening (n 1) fluids. An approximate drag value is found for n < 1 and is compared with computed results.