Studies on the dynamics and rheology of dilute suspensions of Brownian particles in simple shear flow under the action of an external force

We present a novel approach to computing the orientation moments and rheological properties of a dilute suspension of spheroids in a simple shear flow at arbitrary Peclct number based on a generalised Langevin equation method. This method differs from the diffusion equation method which is commonly used to model similar systems in that the actual equations of motion for the orientations of the individual particles are used in the computations, instead of a solution of the diffusion equation of the system. It also differs from the method of 'Brownian dynamics simulations' in that the equations used for the simulations are deterministic differential equations even in the presence of noise, and not stochastic differential equations as in Brownian dynamics simulations. One advantage of the present approach over the Fokker-Planck equation formalism is that it employs a common strategy that can be applied across a wide range of shear and diffusion parameters. Also, since deterministic differential equations are easier to simulate than stochastic differential equations, the Langevin equation method presented in this work is more efficient and less computationally intensive than Brownian dynamics simulations.We derive the Langevin equations governing the orientations of the particles in the suspension and evolve a procedure for obtaining the equation of motion for any orientation moment. A computational technique is described for simulating the orientation moments dynamically from a set of time-averaged Langevin equations, which can be used to obtain the moments when the governing equations are harder to solve analytically. The results obtained using this method are in good agreement with those available in the literature.The above computational method is also used to investigate the effect of rotational Brownian motion on the rheology of the suspension under the action of an external force field. The force field is assumed to be either constant or periodic. In the case of con- I stant external fields earlier results in the literature are reproduced, while for the case of periodic forcing certain parametric regimes corresponding to weak Brownian diffusion are identified where the rheological parameters evolve chaotically and settle onto a low dimensional attractor. The response of the system to variations in the magnitude and orientation of the force field and strength of diffusion is also analyzed through numerical experiments. It is also demonstrated that the aperiodic behaviour exhibited by the system could not have been picked up by the diffusion equation approach as presently used in the literature.The main contributions of this work include the preparation of the basic framework for applying the Langevin method to standard flow problems, quantification of rotary Brownian effects by using the new method, the paired-moment scheme for computing the moments and its use in solving an otherwise intractable problem especially in the limit of small Brownian motion where the problem becomes singular, and a demonstration of how systems governed by a Fokker-Planck equation can be explored for possible chaotic behaviour.

The Inertio-Elastic Planar Entry Flow of Low-Viscosity Elastic Fluids in Micro-fabricated Geometries

The non-Newtonian flow of dilute aqueous polyethylene oxide (PEO) solutions through microfabricated planar abrupt contraction-expansions is investigated. The contraction geometries are fabricated from a high-resolution chrome mask and cross-linked PDMS gels using the tools of soft-lithography. The small length scales and high deformation rates in the contraction throat lead to significant extensional flow effects even with dilute polymer solutions having time constants on the order of milliseconds. The dimensionless extra pressure drop across the contraction increases by more than 200% and is accompanied by significant upstream vortex growth. Streak photography and videomicroscopy using epifluorescent particles shows that the flow ultimately becomes unstable and three-dimensional. The moderate Reynolds numbers (0.03 ≤ Re ≤ 44) associated with these high Deborah number (0 ≤ De ≤ 600) microfluidic flows results in the exploration of new regions of the Re-De parameter space in which the effects of both elasticity and inertia can be observed. Understanding such interactions will be increasingly important in microfluidic applications involving complex fluids and can best be interpreted in terms of the elasticity number, El = De/Re, which is independent of the flow kinematics and depends only on the fluid rheology and the characteristic size of the device.

Estudi i simulació aerodinàmica amb programari de dinàmica de fluids computacional (CFD) en un model bàsic de vehicle de competició de la categoria CM (curses de muntanya)

Aquest projecte té com a objectiu la simulació numérica de la carrosseria d’ un vehicle de curses de muntanya de categoria CM

QUAGMIRE v1.3: a quasi-geostrophic model for investigating rotating fluids experiments

QUAGMIRE is a quasi-geostrophic numerical model for performing fast, high-resolution simulations of multi-layer rotating annulus laboratory experiments on a desktop personal computer. The model uses a hybrid finite-difference/spectral approach to numerically integrate the coupled nonlinear partial differential equations of motion in cylindrical geometry in each layer. Version 1.3 implements the special case of two fluid layers of equal resting depths. The flow is forced either by a differentially rotating lid, or by relaxation to specified streamfunction or potential vorticity fields, or both. Dissipation is achieved through Ekman layer pumping and suction at the horizontal boundaries, including the internal interface. The effects of weak interfacial tension are included, as well as the linear topographic beta-effect and the quadratic centripetal beta-effect. Stochastic forcing may optionally be activated, to represent approximately the effects of random unresolved features. A leapfrog time stepping scheme is used, with a Robert filter. Flows simulated by the model agree well with those observed in the corresponding laboratory experiments.

Linear and nonlinear microrheology of dense colloidal suspensions

The length and time scales accessible to optical tweezers make them an ideal tool for the examination of colloidal systems. Embedded high-refractive-index tracer particles in an index-matched hard sphere suspension provide 'handles' within the system to investigate the mechanical behaviour. Passive observations of the motion of a single probe particle give information about the linear response behaviour of the system, which can be linked to the macroscopic frequency-dependent viscous and elastic moduli of the suspension. Separate 'dragging' experiments allow observation of a sample's nonlinear response to an applied stress on a particle-by particle basis. Optical force measurements have given new data about the dynamics of phase transitions and particle interactions; an example in this study is the transition from liquid-like to solid-like behaviour, and the emergence of a yield stress and other effects attributable to nearest-neighbour caging effects. The forces needed to break such cages and the frequency of these cage breaking events are investigated in detail for systems close to the glass transition.

Concentration of aluminium in breast cyst fluids collected from women affected by gross cystic breast disease

Gross cystic breast disease (GCBD) is the most common benign breast disorder, but the molecular basis of cyst formation remains to be identified. If the use of aluminium-based antiperspirant salts is involved in the etiology of gross breast cyst formation, it might be expected that aluminium would be at elevated levels in human breast cyst fluid (BCF). Aluminium was measured by ICP-MS in 48 samples of BCF, 30 samples of human blood serum and 45 samples of human breast milk at different stages of lactation (colostrum, intermediate, mature). The median level of aluminium in apocrine type I BCF (n:= 27, 150 mu g I-1) was significantly higher than in transuclative type II BCF (n = 21, 32 mu g I-1; P < 0.0001). By comparison, aluminium measurements gave a median concentration of 6 mu g I-1 in human serum and 25 mu g I-1 in human breast milk, with no difference between colostrum, intermediate and mature milk. Levels of aluminium were significantly higher in both types of BCF than in human serum (P < 0.0001). However when compared with human breast milk, aluminium levels were only significantly higher in apocrine type I BCF (P < 0.0001) and not in transudative type II BCF (P = 0.152). It remains to be identified why such high levels of aluminium were found in the apocrine type I BCF and from where the aluminium originated. However, if aluminium-based antiperspirants are found to be the source and to play any causal role in development of breast cysts, then it might become possible to prevent this common breast disorder. Copyright (C) 2008 John Wiley & Sons, Ltd.

Analysis of aluminium content and iron homeostasis in nipple aspirate fluids from healthy women and breast cancer-affected patients

Aluminium is not a physiological component of the breast but has been measured recently in human breast tissues and breast cyst fluids at levels above those found in blood serum or milk. Since the presence of aluminium can lead to iron dyshomeostasis, levels of aluminium and iron-binding proteins (ferritin, transferrin) were measured in nipple aspirate fluid (NAF), a fluid present in the breast duct tree and mirroring the breast microenvironment. NAFs were collected noninvasively from healthy women (NoCancer; n = 16) and breast cancer-affected women (Cancer; n = 19), and compared with levels in serum (n = 15) and milk (n = 45) from healthy subjects. The mean level of aluminium, measured by ICP-mass spectrometry, was significantly higher in Cancer NAF (268.4 ± 28.1 μg l−1; n = 19) than in NoCancer NAF (131.3 ± 9.6 μg l−1; n = 16; P < 0.0001). The mean level of ferritin, measured through immunoassay, was also found to be higher in Cancer NAF (280.0 ± 32.3 μg l−1) than in NoCancer NAF (55.5 ± 7.2 μg l−1), and furthermore, a positive correlation was found between levels of aluminium and ferritin in the Cancer NAF (correlation coefficient R = 0.94, P < 0.001). These results may suggest a role for raised levels of aluminium and modulation of proteins that regulate iron homeostasis as biomarkers for identification of women at higher risk of developing breast cancer. The reasons for the high levels of aluminium in NAF remain unknown but possibilities include either exposure to aluminium-based antiperspirant salts in the adjacent underarm area and/or preferential accumulation of aluminium by breast tissues.

Thermogravimetric analysis of water release from wheat flour and wheat bran suspensions

Bran is hygroscopic and competes actively for water with other key components in baked cereal products like starch and gluten. Thermogravimetric analysis (TGA) of flour–water mixtures enriched with bran at different incorporation levels was performed to characterise the release of compartmentalised water. TGA investigations showed that the presence of bran increased compartmentalised water, with the measurement of an increase of total water loss from 58.30 ± 1.93% for flour only systems to 71.80 ± 0.37% in formulations comprising 25% w/w bran. Deconvolution of TGA profiles showed an alteration of the distribution of free and bound water, and its interaction with starch and gluten, within the formulations. TGA profiles showed that water release from bran-enriched flour is a prolonged event with respect to the release from non-enriched flour, which suggests the possibility that bran may interrupt the normal characteristic processes of texture formation that occur in non-enriched products.

Available potential energy and exergy in stratified fluids

Lorenz’s theory of available p otential energy (APE) remains the main framework for studying the atmospheric and oceanic energy cycles. Because the APE generation rate is the volume integral of a thermodynamic efficiency times the local diabatic heating/cooling rate, APE theory is often regarded as an extension of the theory of heat engines. Available energetics in classical thermodynamics, however, usually relies on the concept of exergy, and is usually measured relative to a reference state maximising entropy at constant energy, whereas APE’s reference state minimises p otential energy at constant entropy. This review seeks to shed light on the two concepts; it covers local formulations of available energetics, alternative views of the dynamics/thermodynamics coupling, APE theory and the second law, APE production/dissipation, extensions to binary fluids, mean/eddy decomp ositions, APE in incompressible fluids, APE and irreversible turbulent mixing, and the role of mechanical forcing on APE production.

A general method for finding extremal states of Hamiltonian dynamical systems, with applications to perfect fluids

In addition to the Hamiltonian functional itself, non-canonical Hamiltonian dynamical systems generally possess integral invariants known as ‘Casimir functionals’. In the case of the Euler equations for a perfect fluid, the Casimir functionals correspond to the vortex topology, whose invariance derives from the particle-relabelling symmetry of the underlying Lagrangian equations of motion. In a recent paper, Vallis, Carnevale & Young (1989) have presented algorithms for finding steady states of the Euler equations that represent extrema of energy subject to given vortex topology, and are therefore stable. The purpose of this note is to point out a very general method for modifying any Hamiltonian dynamical system into an algorithm that is analogous to those of Vallis etal. in that it will systematically increase or decrease the energy of the system while preserving all of the Casimir invariants. By incorporating momentum into the extremization procedure, the algorithm is able to find steadily-translating as well as steady stable states. The method is applied to a variety of perfect-fluid systems, including Euler flow as well as compressible and incompressible stratified flow.