Boundary conditions at the solid-liquid surface over the multiscale channel size from nanometer to micron

The boundary condition at the solid surface is one of the important problems for the microfluidics. In this paper we study the effects of the channel sizes on the boundary conditions (BC), using the hybrid computation scheme adjoining the molecular dynamics (MD) simulations and the continuum fluid mechanics. We could reproduce the three types of boundary conditions (slip, no-slip and locking) over the multiscale channel sizes. The slip lengths are found to be mainly dependent on the interfacial parameters with the fixed apparent shear rate. The channel size has little effects on the slip lengths if the size is above a critical value within a couple of tens of molecular diameters. We explore the liquid particle distributions nearest the solid walls and found that the slip boundary condition always corresponds to the uniform liquid particle distributions parallel to the solid walls, while the no-slip or locking boundary conditions correspond to the ordered liquid structures close to the solid walls. The slip, no-slip and locking interfacial parameters yield the positive, zero and negative slip lengths respectively. The three types of boundary conditions existing in "microscale" still occur in "macroscale". However, the slip lengths weakly dependent on the channel sizes yield the real shear rates and the slip velocity relative to the solid wall traveling speed approaching those with the no-slip boundary condition when the channel size is larger than thousands of liquid molecular diameters for all of the three types of interfacial parameters, leading to the quasi-no-slip boundary conditions.

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

The research of rheology and thermodynamics of organic-inorganic hybrid membrane during the membrane formation

To clarify the mechanism of organic-inorganic hybrid membrane formation by phase-inversion method, the thermodynamical and theological properties of PSF/TiO2 casting solution were investigated by the viscosity measurement and the triangle phase diagram, respectively. TiO2 introduction decreased the non-solvent tolerance of casting solution with non-solvent 20% ethanol aqueous solution, which caused thermodynamic enhancement of phase separation, and also resulted in the change of theological properties from Newtonian fluid to non-Newtonian fluid and the viscosity increase of casting solution, which induced rheological hindrance in demixing process

The mixing of bubbles in two-dimensional bidisperse foams under extensional shear

Cox, S.J. (2006) The mixing of bubbles in two-dimensional bidisperse foams under extensional shear. Journal of Non-Newtonian Fluid Mechanics . 137:39-45.

Contraction/expansion flows: The pressure drop and related issues

Binding, David; Phillips, P.M.; Philips, T.N., (2006) 'Contraction/expansion flows: The pressure drop and related issues', Journal of Non-Newtonian Fluid Mechanics 137 pp.31-38 RAE2008

Viscoelastic flow through a planar contraction using a semi-Lagrangian finite volume method

A semi-Lagrangian finite volume scheme for solving viscoelastic flow problems is presented. A staggered grid arrangement is used in which the dependent variables are located at different mesh points in the computational domain. The convection terms in the momentum and constitutive equations are treated using a semi-Lagrangian approach in which particles on a regular grid are traced backwards over a single time-step. The method is applied to the 4 : 1 planar contraction problem for an Oldroyd B fluid for both creeping and inertial flow conditions. The development of vortex behaviour with increasing values of We is analyzed.

Computational modelling of metal extrusion and forging processes

The computational modelling of extrusion and forging processes is now well established. There are two main approaches: Lagrangian and Eulerian. The first has considerable complexities associated with remeshing, especially when the code is parallelised. The second approach means that the mould has to be assumed to be entirely rigid and this may not be the case. In this paper, a novel approach is described which utilises finite volume methods on unstructured meshes. This approach involves the solution of free surface non-Newtonian fluid flow equations in an Eulerian context to track the behaviour of the workpiece and its extrusion/forging, and the solution of the solid mechanics equations in the Lagrangian context to predict the deformation/stress behaviour of the die. Test cases for modelling extrusion and forging problems using this approach will be presented.

A mixed Eulerian-Lagrangian approach for metal extrusion

In this paper a mixed Eulerian-Lagrangian approach for the modelling metal extrusion processes is presented. The approach involves the solution of non-Newtonian fluid flow equations in an Eulerian context, using a free-surface algorithm to track the behaviour of the workpiece and its extrusion. The solid mechanics equations associated with the tools are solved in Lagangrian context. Thermal interactions between the workpiece are modelled and a fluid-structure interaction technique is employed to model the effect of the fluid traction load imposed by the workpiece on the tools. Two extrusion test cases are investigated and the results obtained show the potential of the model with regard to representing the physics of the process and the simulation time.

A cell-centred finite volume method for modelling viscoelastic flow

An unstructured cell-centred finite volume method for modelling viscoelastic flow is presented. The method is applied to the flow through a planar channel and the 4:1 planar contraction for creeping flow of an Oldroyd-B fluid. The results are presented for a range of Weissenberg numbers. In the case of the planar channel results are compared with analytical solutions. For the 4:1 planar contraction benchmark problem the convection terms in the constitutive equations are approximated using both first and second order differencing schemes to compare the techniques and the effect of mesh refinement on the solution is investigated. This is the first time that a fully unstructured, cell-centredfinitevolume technique has been used to model the Oldroyd-B fluid for the test cases presented in this paper.

Rheological and heat transfer behaviour of the ionic liquid, [C(4)mim][NTf2]

Systematic experiments have been carried out on the thermal and rheological behaviour of the ionic liquid, 1-butyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl} imide, [C(4)mim][NTf2], and, for the first time, on the forced convective heat transfer of an ionic liquid under the laminar flow conditions. The results show that the thermal conductivity of the ionic liquid is similar to 0.13 W m(-1) K-1, which is almost independent of temperature between 25 and 40 degrees C. Rheological measurements show that the [C(4)mim][NTf2] liquid is a Newtonian fluid with its shear viscosity decreasing with increasing temperature according to the exponential law over a temperature range of 20-90 degrees C. The convective heat transfer experiments demonstrate that the thermal entrance length of the ionic liquid is very large due to its high viscosity and low thermal conductivity. The convective heat transfer coefficient is observed to be much lower than that of distilled water under the same conditions. The convective heat transfer data are also found to fit well to the convectional Shah's equation under the conditions of this work. (C) 2007 Elsevier Inc. All rights reserved.

Avaluació del procés de desintegració de papers recuperats

La desintegració és una etapa important en la recuperació de paper vell, ja que té importants conseqüències en consum d'energia i en el comportament de les etapes posteriors. Per això els objectius es centren en analitzar la desintegració des del punt de vista del temps de desintegració, els aspectes energètics, modelització de la màquina de desintegració utilitzada i anàlisi dels factors de cisallament calculats com a mesura global de les forces implicades en la desintegració. Els autors que hi han treballat donen diferents explicacions a aquestes forces. Fins avui només s'ha pogut avaluar qualitativament la influència que tenen cada un dels mecanismes en el temps necessari per a desintegrar i en el consum energètic. Les característiques reològiques de les suspensions papereres, i el seu comportament no newtonià tenen una clara influència en el consum energètic i les forces de desfibrat en el desintegrador. Els experiments de desintegració s'han realitzat en un púlper convencional, amb tres tipus de paper recuperat: paper estucat d'alta qualitat imprès offset (PQ), paper revista estucat imprès en color (PR), paper blanc imprès en impresora làsser (PF). Anàlisi del temps de desintegració Per cada un del papers estudiats (PQ, PR i PF), les fraccions màssiques des de 0.06 fins a la màxima que estat possible per cada paper (de 0.14 a 0.18), i a dues velocitats d'agitació diferents, s'ha determinat el temps de desintegració (tD) fins a aconseguir un índex de Sommerville de 0.01%. S'obté que en augmentar la fracció màssica disminueix potencialment el temps de desintegració. S'ha estudiat la velocitat de desintegració, la producció teòrica del púlper en cada cas, i la seva relació amb les forces d'impacte i de fregament que produeixen la desintegració. Aspectes energètics El consum específic d'energia (SEC), definit com l'energia consumida per a desintegrar 1 kg de paper recuperat, disminueix molt en augmentar Xm, ja que a més de disminuir l'energia consumida en cada desintegració, el contingut en paper és més elevat. Pel disseny de desintegradors, cal tenir en compte que en augmentar Xm i en augmentar la velocitat, sempre augmenta la potència consumida. Però així com els beneficis de treballar a Xm alt són de 10 vegades en termes de SEC i de producció, l'augment de potència és només de l'ordre de 2 vegades la necessària respecte de la Xm baixa. Viscositat aparent i energia de fluidització S'estudia la relació entre el temps de desintegració, les forces de fregament i els valors de viscositat aparent de la bibliografia. Per cada paper i velocitat s'ha observat que el consum específic d'energia disminueix en funció de la viscositat aparent. Reologia del púlper Utilitzant el mètode de Metzner i Otto (1957) per determinar la viscositat aparent mitjana de les suspensions papereres, modificat per Roustan, s'ha caracteritzat el pulper mitjançant el model: Np= K· Rex·Fry S'han utilitzat dissolucions de glicerina com a fluid newtonià per a calcular les constants d'ajust, i a partir d'aquí, aïllar la viscositat aparent en funció de la potència neta i els paràmetres d'agitació. La viscositat aparent, d'acord amb Fabry (1999) es substitueix pel concepte de factor de cisallament. Factor de cisallament Calculat el factor de cisallament per a cada paperot i condicions d'agitació, s'ha relacionat amb Xm, SEC, tD, consum de potència, potència instal·lada i fracció cel·lulòsica. El factor de cisallament és un paràmetre útil per a quantificar les forces globals implicades en la desintegració.

The physics of baking: rheological and polymer molecular structure-function relationships in breadmaking

Molecular size and structure of the gluten polymers that make up the major structural components of wheat are related to their rheological properties via modem polymer rheology concepts. Interactions between polymer chain entanglements and branching are seen to be the key mechanisms determining the rheology of HMW polymers. Recent work confirms the observation that dynamic shear plateau modulus is essentially independent of variations in MW amongst wheat varieties of varying baking performance and is not related to variations in baking performance, and that it is not the size of the soluble glutenin polymers, but the structural and rheological properties of the insoluble polymer fraction that are mainly responsible for variations in baking performance. The rheological properties of gas cell walls in bread doughs are considered to be important in relation to their stability and gas retention during proof and baking, in particular their extensional strain hardening properties. Large deformation rheological properties of gas cell walls were measured using biaxial extension for a number of doughs of varying breadmaking quality at constant strain rate and elevated temperatures in the range 25-60 degrees C. Strain hardening and failure strain of cell walls were both seen to decrease with temperature, with cell walls in good breadmaking doughs remaining stable and retaining their strain hardening properties to higher temperatures (60 degrees C), whilst the cell walls of poor breadmaking doughs became unstable at lower temperatures (45-50 degrees C) and had lower strain hardening. Strain hardening measured at 50 degrees C gave good correlations with baking volume, with the best correlations achieved between those rheological measurements and baking tests which used similar mixing conditions. As predicted by the Considere failure criterion, a strain hardening value of I defines a region below which gas cell walls become unstable, and discriminates well between the baking quality of a range of commercial flour blends of varying quality. This indicates that the stability of gas cell walls during baking is strongly related to their strain hardening properties, and that extensional rheological measurements can be used as predictors of baking quality. (C) 2004 Elsevier B.V. All rights reserved.

Whither tube theory: From believing to measuring

This short contribution examines the difficulties that have not yet been fully overcome in the many developments made from the simplest (and original) tube model for entangled polymers. It is concluded that many more length scales have to be considered sequentially when deriving a continuum rheological model from molecular considerations than have been considered in the past. In particular, most unresolved issues of the tube theory are related to the length scales of tube diameter, and molecular dynamics simulations is the perfect route to resolve them. The power of molecular simulations is illustrated by two examples: stress contributions from bonded and non-bonded interaction, and the inter-chain coupling, which is usually neglected in the tube theory.

A novel X-ray rheometer for in situ studies of polymer melts and solutions during shear flow

A novel X-ray rheometer based on a parallel plate geometry is described. This system allows time-resolved X-ray scattering intensity data to be obtained from polymeric samples subjected to shear flow. The range of quantitative structural parameters, such as molecular orientation and inter chain correlations, which can be obtained from the data is highlighted. Examples of the utility of X-ray scattering in examining optically opaque samples and the extraction of 〈P2〉 and 〈P4〉 orientation parameters are given using anisotropic hydroxypropylcellulose solutions as the sample.