Rheology of a dilute suspension of spheres in a viscoelastic fluid under LAOS
|
|
- Ethel Norman
- 5 years ago
- Views:
Transcription
1 Rheology of a dilute suspension of spheres in a viscoelastic fluid under LAOS GAETANO D AVINO 1, MARTIEN A. HULSEN 2, FRANCESCO GRECO 3 AND PIER LUCA MAFFETTONE 1 1 Dipartimento di Ingegneria Chimica, Università degli Studi di Napoli Federico II, Piazzale Tecchio 8, 8125 Napoli, Italia 2 Department of Mechanical Engineering, Eindhoven University of Technology, PO Box 513, 56 MB Eindhoven, The Netherlands 3 Istituto di Ricerche sulla Combustione, IRC, CNR, Piazzale Tecchio 8, 8125 Napoli, Italia Abstract: Dilute suspensions of spheres in a viscoelastic liquid under oscillatory shear flow are studied via numerical simulations. In linear viscoelastic regime, the 3D finite element simulation results agree with the analytical results from Palierne (199). The analysis is extended to Large Amplitude Oscillations (LAOS). The rotation rate of the sphere as well as the bulk rheology are investigated. The moduli increase when particles are added into the fluid, in agreement with experimental findings. In the nonlinear range both moduli are found to be decreasing functions of the forcing amplitude. The non-linear response of the material is also illustrated by performing a Fourier analysis of the bulk shear stress signal which shows the appearance of higher harmonic contributions. Keywords: Suspensions, Rheology, LAOS, Oscillatory shear flow, Simulations 1. INTRODUCTION Suspensions of particles in liquids are of significant industrial interest. Polymers with fillers, ceramic pastes, biomedical materials, food, cosmetics, and detergents are examples of typical applications of such materials. New applications with nanosized fillers are also emerging, with the suspending fluid being in most cases a nonnewtonian liquid (e.g., nanocomposite corresponding author. pierluca.maffettone@unina.it 1
2 polymer melts). It is then quite obvious that understanding of the rheological and mechanical properties of these materials plays a crucial role in the processing stage. Suspensions have been intensively studied in the past. The first study on suspensions of spheres dates back to Einstein (196) who gave the prediction of the suspension viscosity under the hypotheses of Newtonian, inertialess fluid and dilute system. Since then, a great deal of works both experimental and theoretical was focused to the characterization of suspensions with rigid or deformable particles in Newtonian liquids, and even beyond the dilute limit (see e.g., Brenner 1972; Larson 1998). The problem of viscoelastic suspending liquid is in general less studied, and mainly on the experimental side (e.g., Barnes 23). From the theoretical point of view, in order to handle viscoelastic suspensions of practical interest, realistic non-newtonian constitutive equations for the suspending liquid need to be taken into account. Of course, due to the mathematical complexity of the equations involved, use of numerical simulation to tackle specific flow problems is generally required. As it is well known, a typical external flow frequently used in experiments is the oscillatory shear flow. According to the magnitude of the imposed strain, two different regimes can be distinguished. If a small forcing amplitude is applied (Small Amplitude Oscillatory Shear, SAOS), the material exhibits the so called linear viscoelastic behaviour, where the rheological quantities are linear function of the strain. By increasing the forcing amplitude (Large Amplitude Oscillatory Shear, LAOS), a non-linear relationship between stress and strain occurs. For an inertialess emulsion with non-newtonian fluids, Palierne (199) derived an analytical expression for the storage (G') and loss (G'') moduli in SAOS regime. A solid suspension with a viscoelastic matrix can be considered as a limiting case of Palierne results for the inner phase viscosity going to infinity. However, we are not aware of any prediction on the moduli of a solid viscoelastic suspension in LAOS regime. Aim of this work is to analyze dilute, inertialess, non-brownian suspensions of rigid spheres in a viscoelastic medium subjected to oscillatory shear flows, including the response to large forcing amplitudes. Simple but realistic viscoelastic constitutive equations are chosen to properly model non-linear phenomena. The analysis is carried out through a cell model for the dilute suspension, by means of 3D numerical simulations in a cell. Once the velocity and pressure fields within the cell are determined, the bulk shear stress (i.e., the stress obtained by considering the suspension as a whole) is calculated. This latter quantity allows the predictions of both moduli. The analytical prediction from Palierne (199) are used to 2
3 validate the numerical code. Nonlinear results are then obtained when large forcing amplitudes are used. Finally, a Fourier analysis is carried out in order to show the arising of higher harmonic contributions when the forcing amplitude is sufficiently large. 2. GOVERNING EQUATIONS The dilute suspension of spherical particles in a viscoelastic medium is described within the framework of a cell model (Brenner, 1972). Suspension diluteness implies that particle do not feel each other, and thus the suspension can be subdivided in regions (cells) each one containing a single particle. In this view, solving the hydrodynamic problem in a single cell is all what is needed to end up with predictions for the dilute suspension as a whole. The imposed flow field is applied at the border of the cell. Under the assumption of incompressibility, negligible inertia, buoyancy free and isothermal conditions, the hydrodynamic problem of a single sphere in an imposed flow field is governed by the following continuity (mass balance) and momentum balance equations: v = p + τ = (1) where p is the pressure, v is the velocity, and τ is the constitutive extra stress tensor. In this work, the following two constitutive equations for viscoelastic liquids are considered: λ τ+ τ = 2η D (2) αλ λ τ+ τ τ+ τ = 2ηD (3) η where τ indicates the Upper Convected (UC) time derivative of τ, λ is the relaxation time, 2D = v + v T is the rate of deformation tensor and η is the zero-shear viscosity. Equation (2) is the UC Maxwell model and Eq. (3) is the Giesekus model. The latter model contains an additional nondimensional parameter, α, that affects the shear thinning and the elongational thickening. Boundary conditions are assigned at the cell boundaries and on the sphere surface. The farfield condition corresponds to the imposed shear flow: p= p B v = v r B B (4) 3
4 with p B a constant pressure, v B the imposed velocity gradient (constant in space), and r B the position vector from the sphere center to the cell border. SAOS and LAOS regimes are achieved by imposing a co-sinusoidal deformation on the border: γ (t) = γ cos(f t) (5) where γ is the forcing strain amplitude [-] and f is the frequency [t -1 ]. From Eq. (5) it follows that the imposed velocity field is a time-dependent, sinusoidal shear flow. Calling x the flow direction, y the velocity gradient direction, and z the vorticity direction, the imposed velocity gradient is given by: vb xy = γfsin(ft) = γsin(ft) (6) the other components being nil. In Eq. (6), γ = γ f [t -1 ] is the imposed shear rate. No slip boundary conditions are imposed at the sphere surface r = R. We assume that the sphere is torque-free, or "freely rotating" (Brunn, 1976; Leal, 1979), i.e. its rotation is only due to the motion of the surrounding fluid. The torque-free boundary condition at r = R is: r T n da = (7) A S where A s is the sphere surface, n is the normal at the sphere surface, the tensor T is T= pi+ τ, and the integral of the local torque r T n da spans the sphere surface. Due to the symmetry of the imposed shearing flow, only the vorticity component z of Eq. (7) is nonzero. The velocity at the sphere surface is known, and can be written as: ( ) v R,t = ω (t) e z R (8) e z being the unit vector along z. The time-dependent velocity v(r,t) in Eq. (8) is at any point on the sphere surface the tangential velocity of that point. Hence, v(r,t) is everywhere tangential to the surface, and parallel to the shear plane xy. The unknown angular velocity ω(t) in Eq. (8) can then be determined through the z-component of the torque-free condition, Eq. (7). Within the cell description of the dilute system, the bulk tangential stress is obtained by volume averaging of the local tangential stress over the cell: τ xy n = τxydv V (9) Tot Vcell 4
5 with n the number of cells, and τ xy is the shear stress in the cell containing one particle only, with the flow field imposed at the cell boundaries. As usual, the decomposition is made: τ φ φ = τ dv + R µ nµ n ds V V τ (1) xy xy x y p VFcell p SP with V Fcell the volume of fluid in the cell, and V P and S P the volume and surface of the spherical particle of radius R. It can be shown that the first term on the rhs of eq. (1) and the nv surface integral are independent of the volume fraction φ (= V implies the linearity of the bulk stress in φ. P Tot ). The diluteness assumption 3. THE NUMERICAL PROCEDURE By solving eqs. (1) with either eq. (2) or (3) and boundary conditions eqs. (4, 6, 7), the pressure and the velocity, plus the angular velocity of the sphere ω are determined. The 3D flow problem is solved by the finite element method (we used the TFEM package (Hulsen, 27)) on a cubic cell containing the sphere. Due to the symmetry of the problem, only one half (z > ) of the full domain can be considered in order to optimize the computational effort. The flow field is imposed on the boundaries of the computational domain. The cell dimension was so chosen as to give cell invariant pressure and velocity fields. A mesh with tetrahedral elements was chosen, with a higher density of elements close to the sphere, where larger gradients are expected. Convergence of the results was checked by refining the mesh. The angular velocity is monitored by varying the number of elements. A mesh size of about 6 elements suffices to obtain reliable results, whereas deviations of about 1% is observed when the number of elements is decreased to 25. Stabilization techniques are used for the discretization of equations to improve numerical convergence. The momentum equation is discretized using elastic viscous stress splitting (DEVSS-G) formulation (Guénette and Fortin, 1995) combined with the streamline upwind Petrov Galerkin (SUPG) method (Brooks and Hughes, 1982) for the constitutive equation. Moreover, the log-conformation representation is also used, leading to a significant improvement of numerical stability (Fattal and Kupferman, 24; Hulsen et al., 25). The procedure implemented allows to perform transient calculations. In particular, the 5
6 momentum and continuity equations are decoupled from the constitutive equation, and an implicit stress formulation is used. In this formulation the time-discretized constitutive equation is substituted into the momentum balance in order to obtain a Stokes like system. The torque-free condition (eq. (7)) is imposed through constraints on the sphere surface, by means of Lagrange multipliers. In this way, the sphere rotation is automatically calculated by solving the augmented system of equations. At each time step, the solution vector for v, p, and G is computed, with G the projected velocity gradient (see Guénette and Fortin, 1995; Bogaerds et al., 22). Subsequently, the solution for the stress τ (in terms of the conformation tensor) is determined, through a combined Euler forward/backward formula of the constitutive equation with convection taken implicitly. A continuous quadratic interpolation for the velocity, a continuous linear interpolation for pressure and a continuous bilinear interpolation for the stress tensor and auxiliary velocity gradient are used. Finally, at each time step, we need to solve two linear nonsymmetric sparse systems. We use the parallel direct solver PARDISO (Schenk and Gartner, 24). Finally, integrals appearing in eq. (1) were calculated with fourth order quadrature formulas. 4. LINEAR VISCOELASTICITY The dynamic viscoelasticity of the dilute suspension under SAOS is of course affected by the presence of the rigid spherical inclusions. Palierne (199) determined the moduli for an emulsion of non-newtonian liquids in the linear viscoelastic regime. The solid particle suspension is recovered as a special case, and his prediction for the complex modulus, G* ( f) = G' ( f) + ig'' ( f) is: ( ) ( f) * G f 1 3/2 G * + φ = 1 φ (11) with * G the complex modulus of the suspending liquid. Equation (11) states that the complex modulus of the suspension is merely a vertical shift of that of the pure liquid. Once expanded to first order in φ, eq. (11) gives in terms of moduli: ( ) G' f /G' (f) 1 2 ( ) 5 = + φ 5 = + φ G" f /G'' (f) 1 2 (12) 6
7 One can note that the normalized storage and loss moduli are independent of frequency forcing. Both moduli are shifted by the same factor, and the equation for G is a generalization of the Einstein result to the linear viscoelastic case. We now validate our numerical procedure by comparing the results of the simulations with Palierne predictions. The linear viscoelasticity condition is satisfied by choosing a sufficienly small amplitude of the forcing (i.e. γ = γ f, see Eq. (6)). In our simulations, γ was chosen 1-3 s -1 : this value was verified to be sufficiently small so that linear regime was attained. It should be pointed out that in linear viscoelasticity the choice of the viscoelastic model, Eq. (2) or (3), does not affect the moduli because Maxwell behaviour is always recovered in such a regime. Figure 1 shows simulated bulk values for G' and G'' (symbols) for a dilute suspension reported as a function of the frequency, for φ =.5. 1 G''.1 G', G''.1.1 G' unfilled fluid suspension (φ =.5) f Fig. 1. G' and G'' as functions of forcing frequency (f), for a dilute suspension with φ =.5 (symbols) and unfilled fluid (solid lines). The constant vertical shift at all frequencies with respect to the pure liquid behaviour (solid lines in Fig. 1) is clearly observed. A quantitative comparison of the code predictions with prediction by Palierne (199) is reported in Fig. 2, where the relative errors ( ε ', ε " ) are reported as a function of the forcing frequency. In the explored frequency range, quantitative agreement with theory is found. r r 7
8 1e-2 ε' r ε'' r ε' r, ε'' r 1e-3 1e f Fig. 2. Absolute values of the relative error of G' (full circles) and G'' (open circles) as function of forcing frequency (f), for a dilute suspension with φ = LARGE AMPLITUDE OSCILLATIONS In this section, the results of our simulations for LAOS regime are presented and discussed. First, the time-dependent rotation rate of the sphere is investigated. Then, we focus the attention on the bulk properties of the suspension, i.e., its bulk shear stress and storage and loss moduli. Let us first consider how the particle rotation rate is influenced by a large forcing amplitude. As it is well known (Einstein, 196; 1911), an isolated sphere immersed in a Newtonian fluid without inertia subjected to a steady shear flow follows without any lag in time the imposed flow. Indeed, the local fields are independent of time and Newtonianness implies that the stress istantaneously develops (no transient occurs). Einstein's analysis showed a particle rotation rate ω = γ /2, where γ is the externally imposed shear rate. If an oscillatory shear flow is applied, an analogous behaviour is expected. When imposing a sinusoidal shear rate (see eq.(6)), the rotation rate is ω = γ f sin(f t) / 2. D'Avino (27) showed a slowing down of the sphere rotation rate under steady shear flow when a viscoelastic suspending matrix is considered. An analogous result is found here for the case of an imposed sinusoidal shear rate in LAOS. Figure 3 shows the rotation rate of the sphere immersed in viscoelastic liquids (Maxwell and Giesekus respectively) as a function of dimensionless time for a forcing frequency f = 1. and at different strain amplitudes. The dashed line refers to the linear viscoelastic curve whereas the curves with symbols are for different strain amplitudes. 8
9 .5.5 γ =.5 γ = 1. γ = 1. γ = γ = 2. linear viscoelasticity.25 γ = 3. linear viscoleasticity ω/γ. ω/γ t/t t/t Fig. 3. Particle rotation rate in a Maxwell fluid (left) and a Giesekus fluid with α=.2 (right) as a function of time. Linear viscoelastic regime (dashed line), and three strains (solid lines with symbols), at f = 1.. From both figures, it appears that the lowering of angular sphere velocities becomes more pronounced as the strain increases. Passing now to the bulk stress calculation in LAOS, one notes that by increasing the forcing strain, the linear viscoelastic behaviour is lost. For a Maxwell suspending fluid, results for bulk shear stress divided by γ at two different strains for f = 1. are reported in Fig. 4. After a transient, a periodic regime is achieved. In the figure (as in all subsequent similar figures), only a period is reported, after the transient has extinguished. 1. linear viscoelastic γ = 1. γ = 2..5 σ xy /γ t/t Fig. 4. Shear stress divided by γ = γ f for a dilute suspension of spheres in a Maxwell fluid as a function of dimensionless time (with T the oscillation period) for linear viscoelastic 9
10 regime (dashed line) and two strains (solid lines with symbols). For low strains (up to.2), the linear viscoelastic behaviour (dashed line) is found. Well above this value of the forcing amplitude, a slight modification of bulk shear stress is observed. Specifically, maxima and minima reduce in absolute values, and are anticipated with respect to the linear viscoelastic curves, see inset. The effect of the solid spheres on the bulk behaviour of the suspension is even clearer if a Giesekus suspending fluid is considered. The bulk shear stress for such a suspension is reported in Fig. 5, for f = 1. and three different strains, with α = linear viscoelasticity γ = 1..5 γ = 2. γ = 3. σ xy /γ t/t Fig. 5. Shear stress divided by γ = γ f for a dilute suspension of spheres in a Giesekus fluid as a function of time for linear viscoelastic regime (dashed line) and three strains (solid lines with symbols). Again, deviations from linear viscoelastic curves are more and more pronounced as higher strain amplitudes are applied. Furthermore, if the forcing amplitude is sufficiently high, the resulting stress is patently no longer sinusoidal as the response stress appears to be strongly distorted. To better show the impact of the presence of the spherical particle on the bulk shear stress, a comparison with the unfilled fluid for f = 1. and γ = 3. is reported in Fig. 6. In the case of a pure Giesekus liquid, the LAOS response shows a distortion with respect to the linear case. The presence of a filler further deepens the difference. Notice that for a Maxwell fluid, the 1
11 distortion at large amplitudes is not observed at all in the unfilled system, hence it is completely due to the presence of solid inclusions. From the calculated bulk shear stresses, loss and storage moduli of the suspension can be derived. Figure 7 shows the normalized moduli for a Maxwell suspending fluid with φ =.5 reported as function of strain for f =.5 and f = 1.. At vanishing strains, the curves start from the linear viscoelastic value (the Palierne prediction) equal to for the volume fraction considered..4 unfilled φ =.5.2 σ xy /γ t/t Fig. 6. Shear stress divided by γ = γ f as a function of time for γ = 3. and f = 1.. Dashed line is the unfilled Giesekus fluid (α=.2); solid line is the dilute suspension of spheres. A decreasing trend for both moduli is found, more pronounced with increasing frequencies. As just remarked, such a trend is completely due to the presence of the solid particles, as for a pure Maxwell fluid no effect of amplitudes is foreseen. 11
12 f = f = G'/G', G''/G'' G'/G', G''/G'' G'/G' 1.8 G'/G' 1.7 G''/G'' 1.7 G''/G'' γ γ Fig. 7. G'/G' and G''/G'' as function of forcing strain, for a dilute suspension of spheres in a Maxwell fluid with φ = f = f = G'/G', G''/G'' G'/G' G''/G'' G'/G', G''/G'' G'/G' G''/G''.2 unfilled G'/G' unfilled G''/G''.2 unfilled G'/G' unfilled G''/G'' γ γ Fig. 8. G'/G' and G''/G'' (lines with symbols) as function of forcing strain, for a dilute suspension of spheres in a Giesekus fluid with φ =.5. The same quantities for an unfilled fluid are reported with lines without symbols. The normalized loss and storage bulk moduli for a suspending Giesekus fluid are reported in Fig. 8. Again, a decreasing trend from the linear viscoelastic limiting values is found. In this case, however, the effect of the inclusions is superimposed on the characteristic non-linear response of a Giesekus fluid in the LAOS regime. The relative impact of solid particles on the bulk response can be captured by comparing the normalized moduli with those of the unfilled Giesekus fluid (grey lines). As expected, the presence of inclusions leads to higher moduli with respect to the unfilled fluid. However, both quantities are decreasing with strain, and the separation between the filled and unfilled fluid curves becomes progressively smaller at high strains. 12
13 It is now worth showing how the presence of an inclusion affects the loss of linearity with increasing the strain. Figure 9 shows the behaviour of the quantities (G' L - G')/G' L and (G'' L -G'')/G'' L as a function of strain. The pedices L indicate the linear value of the moduli, thus these quantities are nil for vanishing amplitudes of the strain. The curves show that a suspension abandons the linear regime at a lower strain amplitude with respect the pure fluid. Moreover, the curves for a suspension are always higher than those for the unfilled fluid, and the deviations from the linear viscoelastic values are more and more pronounced as γ increases..7.3 (G' L - G')/G' L suspension unfilled fluid (G'' L - G'')/G'' L suspension unfilled fluid γ γ Fig. 9. (G' L - G')/G' L and (G'' L - G'')/G'' L as function of forcing strain γ, for a dilute suspension of spheres in a Giesekus fluid with φ =.5 (black lines). The same quantities for an unfilled fluid are reported with grey lines. The distorsion of the tangential stress waveform from the sinusoidal one, already reported in Fig. 5 for a Giesekus fluid, can be better appreciated by considering the power spectrum as reported in Fig. 1a. Indeed, in the Fourier transformed shear stress signal, higher order odd harmonics appear when the strain amplitude is sufficiently high, as shown in Fig. 1a. Nonlinearities in the response are signalled by the appearance of a nonzero third harmonic (note that even harmonics are nil by symmetry of the tangential stress). Even for dilute systems, one observes a third harmonic around 1% of the fundamental one. In the zoomed inset, also hints of nonnegligible 5 th harmonic are visible. Figure 1b shows the trend of the amplitude of the fundamental and of the third harmonics as a function of the strain. The limiting slope of the fundamental harmonic for vanishing strain indicates, of course, the linear viscoelastic limit, and is correctly captured by the simulations. The curve of the fundamental 13
14 harmonic is a cubic. The third harmonic vanishes at small strain, and increases quadratically with γ γ G* ω 5ω γ G*.8.6 ω 3ω Linear viscoelasticity.4.2 a γ =.5 γ = 1. γ = 2. γ = b. ω 3ω 5ω 7ω γ Fig. 1. Complex modulus (times the strain amplitude) versus frequency for different strain amplitudes (a) and versus the strain amplitude (b)for the first and third harmonic contribution. Data are for a dilute suspension of spheres in a Giesekus (α=.2) fluid with f = 1.. See Fig CONCLUSIONS In this paper, we studied a dilute, non-brownian, inertialess suspension of spheres in a viscoelastic fluid, subjected to an oscillatory shear flow. The motion of the sphere as well as the bulk rheology are investigated. The numerical analysis encompasses both small and large amplitude oscillation regimes. The governing equations are discretized by the finite element method over an unstructured grid surrounding the particle. The rigid-body motion is imposed on the sphere surface through Lagrange multipliers constraints. The bulk stress is obtained by resorting to a cell model. The 3D numerical simulations are validated by comparison with the analytical predictions by Palierne (199). The trends of G' and G'' are investigated for different forcing frequencies and amplitudes. The behaviour is found to be strongly dependent on the constitutive model. However, we found that the presence of the inclusion always increases both moduli. Moreover, for high strains, when the linear viscoelastic regime is abandoned, a decreasing trend of G' and G'' is found in 14
15 agreement with experimental results (Barnes, 23). It is important to point out that such an effect is evident even if a Maxwell suspending fluid is considered, in spite of the fact that the corresponding pure fluid does not exhibit any thinning of both moduli. Finally, a Fourier analysis is performed for the Giesekus model, showing a clear nonzero third harmonic contribution, more and more pronounced with increasing the forcing amplitude. BIBLIOGRAPHY Barnes H. A. in Rheology Reviews, 1, British Society of Rheology (23). Batchelor G.K., The stress system in a suspension of force-free particles, J. Fluid Mech. 41, (197). Bogaerds A.C.B., A.M. Grillet, G.W.M. Peters, F.P.T. Baaijens, Stability analysis of polymer shear flows using the extended pom-pom constitutive equations, J. Non- Newtonian Fluid Mech. 18, (22). Brenner, H. Suspension rheology Proc. Heat and Mass Transfer , (1972). Brooks A.N., T.J.R. Hughes, Streamline upwind/petrov-galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier- Stokes equations, Comp. Meth. Appl. Mech. Eng. 32, (1982). Brunn P.O., Slow motion of a sphere in a second-order fluid, Rheol. Acta 15, (1976). Einstein A., Eine neue bestimmung der moleküldimensionen, Ann. Physik 19, (196); ibidem 34, (1911). Fattal R., R. Kupferman, Constitutive laws for the matrix-logarithm of the conformation tensor, J. Non-Newtonian Fluid Mech., 123, (24). D Avino G., Computational rheology of solid suspensions, PhD Thesis, University of Naples, (27) Guénette R., M. Fortin, A new mixed finite element method for computing viscoelastic flows, J. Non-Newtonian Fluid Mech. 6, (1995). Hulsen M.A., R. Fattal, R. Kupferman, Flow of viscoelastic fluids past a cylinder at high Weissenberg number: stabilized simulations using matrix logarithms, J. Non- Newtonian Fluid Mech., 127, 27-39, (25). Hulsen M.A., TFEM: A toolkit for the finite element method, Userguide. (27). Larson R. G. The Structure and Rheology of Complex Fluids Oxford University Press. (1998) Leal L. G., The motion of small particles in nonnewtonian fluids, J. Non-Newtonian Fluid 15
16 Mech., 5, (1979). Palierne J.F., Linear rheology of viscoelastic emulsions with interfacial tension, Rheol. Acta, 29, (199). Schenk O., K. Gartner, Solving unsymmetric sparse systems of linear equations with PARDISO, J. of Future Generation Computer Systems, 2, (24). 16
Journal of Non-Newtonian Fluid Mechanics
J. Non-Newtonian Fluid Mech. 157 (2009) 101 107 Contents lists available at ScienceDirect Journal of Non-Newtonian Fluid Mechanics journal homepage: www.elsevier.com/locate/jnnfm Effects of confinement
More informationCOMPUTATIONAL RHEOLOGY OF SOLID SUSPENSIONS
Università degli Studi di Napoli Federico II Dipartimento di Ingegneria Chimica Dottorato in Ingegneria Chimica XX Ciclo Tesi di Dottorato COMPUTATIONAL RHEOLOGY OF SOLID SUSPENSIONS Gaetano D Avino Comitato
More informationEffect of Viscoelasticity on the Rotation of a Sphere in Shear Flow
Effect of Viscoelasticity on the Rotation of a Sphere in Shear Flow F. Snijkers, G. D Avino Department of Chemical Engineering, K.U. Leuven, W. de Croylaan 46, B-3001 Leuven, Belgium P.L. Maffettone Dipartimento
More informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More informationRHEOLOGY Principles, Measurements, and Applications. Christopher W. Macosko
RHEOLOGY Principles, Measurements, and Applications I -56081-5'79~5 1994 VCH Publishers. Inc. New York Part I. CONSTITUTIVE RELATIONS 1 1 l Elastic Solid 5 1.1 Introduction 5 1.2 The Stress Tensor 8 1.2.1
More informationNumerical result of complex quick time behavior of viscoelastic fluids in flow domains with traction boundaries
Korea-Australia Rheology Journal Vol. 19, No. 4, December 2007 pp. 211-219 Numerical result of complex quick time behavior of viscoelastic fluids in flow domains with traction boundaries Youngdon Kwon*
More informationViscosity of magmas containing highly deformable bubbles
Journal of Volcanology and Geothermal Research 105 (2001) 19±24 www.elsevier.nl/locate/jvolgeores Viscosity of magmas containing highly deformable bubbles M. Manga a, *, M. Loewenberg b a Department of
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationThe Large Amplitude Oscillatory Strain Response of Aqueous Foam: Strain Localization and Full Stress Fourier Spectrum
The Large Amplitude Oscillatory Strain Response of Aqueous Foam: Strain Localization and Full Stress Fourier Spectrum By F. Rouyer, S. Cohen-Addad, R. Höhler, P. Sollich, and S.M. Fielding The European
More informationCESAME, Louvain, Belgium
Thursday Morning Thursday :2 Johann-Peter-Hebel FM3 PGD-based parametric solutions of double-scale Fokker-Planck equations: application to suspensions involving aggregates of rods Francisco Chinesta, Emmanuelle
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationHEAT TRANSFER OF SIMPLIFIED PHAN-THIEN TANNER FLUIDS IN PIPES AND CHANNELS
HEAT TRANSFER OF SIMPLIFIED PHAN-THIEN TANNER FLUIDS IN PIPES AND CHANNELS Paulo J. Oliveira Departamento de Engenharia Electromecânica, Universidade da Beira Interior Rua Marquês D'Ávila e Bolama, 600
More informationModeling the Rheology and Orientation Distribution of Short Glass Fibers Suspended in Polymeric Fluids: Simple Shear Flow
Modeling the Rheology and Orientation Distribution of Short Glass Fibers Suspended in Polymeric Fluids: Simple Shear Flow Aaron P.R. berle, Donald G. Baird, and Peter Wapperom* Departments of Chemical
More informationApplication Of Optimal Homotopy Asymptotic Method For Non- Newtonian Fluid Flow In A Vertical Annulus
Application Of Optimal Homotopy Asymptotic Method For Non- Newtonian Fluid Flow In A Vertical Annulus T.S.L Radhika, Aditya Vikram Singh Abstract In this paper, the flow of an incompressible non Newtonian
More informationContents. Preface XIII. 1 General Introduction 1 References 6
VII Contents Preface XIII 1 General Introduction 1 References 6 2 Interparticle Interactions and Their Combination 7 2.1 Hard-Sphere Interaction 7 2.2 Soft or Electrostatic Interaction 7 2.3 Steric Interaction
More informationEVALUATION OF NONLINEAR DIFFERENTIAL MODELS FOR THE SIMULATION OF POLYMER MELTS
1 th Fall Rubber Colloquium EVALUATION OF NONLINEAR DIFFERENTIAL MODELS FOR THE SIMULATION OF POLYMER MELTS Jochen Kroll, Stefan Turek, Patrick Westervoß Institute of Applied Mathematics (LS III), TU Dortmund
More informationLecture 2: Constitutive Relations
Lecture 2: Constitutive Relations E. J. Hinch 1 Introduction This lecture discusses equations of motion for non-newtonian fluids. Any fluid must satisfy conservation of momentum ρ Du = p + σ + ρg (1) Dt
More informationSupplementary material to On the rheology of pendular gels and morphological developments in paste- like ternary systems based on capillary attraction
Electronic Supplementary Material (ESI) for Soft Matter. This journal is The Royal Society of Chemistry 214 Supplementary material to On the rheology of pendular gels and morphological developments in
More informationModeling of Suspension Flow in Pipes and Rheometers
Modeling of Suspension Flow in Pipes and Rheometers Nicos S. Martys, Chiara F. Ferraris, William L. George National Institute of Standards and Technology Abstract: Measurement and prediction of the flow
More informationFEM techniques for nonlinear fluids
FEM techniques for nonlinear fluids From non-isothermal, pressure and shear dependent viscosity models to viscoelastic flow A. Ouazzi, H. Damanik, S. Turek Institute of Applied Mathematics, LS III, TU
More informationWall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes
Excerpt from the Proceedings of the COMSOL Conference 9 Boston Wall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes Daoyun Song *1, Rakesh K. Gupta 1 and Rajendra P. Chhabra
More informationELASTIC INSTABILITIES IN CONE{AND{PLATE FLOW: SMALL GAP THEORY. David O. Olagunju. University of Delaware. Newark, DE 19716
ELASTIC INSTABILITIES IN CONE{AND{PLATE FLOW: SMALL GAP THEORY David O. Olagunju Department of Mathematical Sciences University of Delaware Newark, DE 19716 June 15, 1995 Abstract Consider the axisymmetric,
More informationRheology of Soft Materials. Rheology
Τ Thomas G. Mason Department of Chemistry and Biochemistry Department of Physics and Astronomy California NanoSystems Institute Τ γ 26 by Thomas G. Mason All rights reserved. γ (t) τ (t) γ τ Δt 2π t γ
More informationTime dependent finite element analysis of the linear stability of viscoelastic flows with interfaces
J. Non-Newtonian Fluid Mech. 116 (23 33 54 Time dependent finite element analysis of the linear stability of viscoelastic flows with interfaces Arjen C.B. Bogaerds, Martien A. Hulsen, Gerrit W.M. Peters,
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationA numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 6 (2016, 152 161 Research Article A numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows Lei Hou a, Yunqing Feng a,, Lin
More informationA simple method for simulating general viscoelastic fluid flows with an alternate log-conformation formulation
J. Non-Newtonian Fluid Mech. 147 (2007) 189 199 A simple method for simulating general viscoelastic fluid flows with an alternate log-conformation formulation Oscar M. Coronado a, Dhruv Arora a, Marek
More informationViscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior
Viscoelasticity Basic Notions & Examples Formalism for Linear Viscoelasticity Simple Models & Mechanical Analogies Non-linear behavior Viscoelastic Behavior Generic Viscoelasticity: exhibition of both
More informationMonolithic FEM multigrid techniques for the simulation of viscoelastic flow
Monolithic FEM multigrid techniques for the simulation of viscoelastic flow A. Ouazzi, H. Damanik, S. Turek, J. Hron Institute of Applied Mathematics, LS III, TU Dortmund http://www.featflow.de European
More informationTWO-DIMENSIONAL SIMULATIONS OF THE EFFECT OF THE RESERVOIR REGION ON THE PRESSURE OSCILLATIONS OBSERVED IN THE STICK-SLIP INSTABILITY REGIME
1 TWO-DIMENSIONAL SIMULATIONS OF THE EFFECT OF THE RESERVOIR REGION ON THE PRESSURE OSCILLATIONS OBSERVED IN THE STICK-SLIP INSTABILITY REGIME Eleni Taliadorou and Georgios Georgiou * Department of Mathematics
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationCM4655 Polymer Rheology Lab. Torsional Shear Flow: Parallel-plate and Cone-and-plate
CM4655 Polymer heology Lab Torsional Shear Flow: Parallel-plate and Cone-and-plate (Steady and SAOS) Professor Faith A. Morrison Department of Chemical Engineering Michigan Technological University r (-plane
More informationRheometry. II.1 Introduction
II Rheometry II.1 Introduction Structured materials are generally composed of microstructures dispersed in a homogeneous phase [30]. These materials usually have a yield stress, i.e. a threshold stress
More informationUniversity Graz / Austria Institut für Chemie Volker Ribitsch
University Graz / Austria Institut für Chemie Volker Ribitsch 1 Rheology Oscillatory experiments Dynamic experiments Deformation of materials under non-steady conditions in the linear viscoelastic range
More informationSimple constitutive models for linear and branched polymers
J. Non-Newtonian Fluid Mech. 116 (2003) 1 17 Simple constitutive models for linear and branched polymers Roger I. Tanner, Simin Nasseri School of Aerospace, Mechanical and Mechatronic Engineering, University
More informationSANDEEP TYAGI and ANUP K. GHOSH*
Linear Viscoelastic and Transient Behavior of Polypropylene and Ethylene Vinyl Acetate Blends: An Evaluation of the Linear Palierne and a Nonlinear Viscoelastic Model for Dispersive Mixtures SANDEEP TYAGI
More informationThe role of micro-scale inertia in transport processes. - Ganesh Subramanian (JNCASR, Bangalore)
The role of micro-scale inertia in transport processes - Ganesh Subramanian (JNCASR, Bangalore) Micro-scale inertia (suspensions and emulsions) Navier-Stokes equations : Re u t + Re u u = Inertial acceleration
More informationModelling the Rheology of Semi-Concentrated Polymeric Composites
THALES Project No 1188 Modelling the Rheology of Semi-Concentrated Polymeric Composites Research Team Evan Mitsoulis (PI), Professor, NTUA, Greece Costas Papoulias (Research Student), NTUA, Greece Souzanna
More informationLes Houches School of Foam: Rheology of Complex Fluids
Les Houches School of Foam: Rheology of Complex Fluids Andrew Belmonte The W. G. Pritchard Laboratories Department of Mathematics, Penn State University 1 Fluid Dynamics (tossing a coin) Les Houches Winter
More informationRheology of cellulose solutions. Puu Cellulose Chemistry Michael Hummel
Rheology of cellulose solutions Puu-23.6080 - Cellulose Chemistry Michael Hummel Contents Steady shear tests Viscous flow behavior and viscosity Newton s law Shear thinning (and critical concentration)
More informationA new numerical framework to simulate viscoelastic free-surface flows with the finitevolume
Journal of Physics: Conference Series PAPER OPEN ACCESS A new numerical framework to simulate viscoelastic free-surface flows with the finitevolume method Related content - Gravitational collapse and topology
More informationStability analysis of injection molding flows
Stability analysis of injection molding flows Arjen C. B. Bogaerds, Martien A. Hulsen, Gerrit W. M. Peters, a) and Frank P. T. Baaijens Materials Technology, Department of Mechanical Engineering, Dutch
More informationLecture 7: Rheology and milli microfluidic
1 and milli microfluidic Introduction In this chapter, we come back to the notion of viscosity, introduced in its simplest form in the chapter 2. We saw that the deformation of a Newtonian fluid under
More informationTHE 3D VISCOELASTIC SIMULATION OF MULTI-LAYER FLOW INSIDE FILM AND SHEET EXTRUSION DIES
THE 3D VISCOELASTIC SIMULATION OF MULTI-LAYER FLOW INSIDE FILM AND SHEET EXTRUSION DIES Kazuya Yokomizo 1, Makoto Iwamura 2 and Hideki Tomiyama 1 1 The Japan Steel Works, LTD., Hiroshima Research Laboratory,
More informationCHAPTER 2 THERMAL EFFECTS IN STOKES SECOND PROBLEM FOR UNSTEADY MICROPOLAR FLUID FLOW THROUGH A POROUS
CHAPTER THERMAL EFFECTS IN STOKES SECOND PROBLEM FOR UNSTEADY MICROPOLAR FLUID FLOW THROUGH A POROUS MEDIUM. Introduction The theory of micropolar fluids introduced by Eringen [34,35], deals with a class
More informationMODENA. Deliverable 3.2. WP s leader: TU/e. Simulations for foams, dispersion and mixing and developed SW. Principal investigator:
Delivery date: 11-4-2015 MODENA Authors: Patrick Anderson, Martien Hulsen, Christos Mitrias TU/e E-mail : pda@tue.nl Deliverable 3.2 Simulations for foams, dispersion and mixing and developed SW Daniele
More informationInfluence of steady shear flow on dynamic viscoelastic properties of un-reinforced and Kevlar, glass fibre reinforced LLDPE
Bull. Mater. Sci., Vol. 27, No. 5, October 2004, pp. 409 415. Indian Academy of Sciences. Influence of steady shear flow on dynamic viscoelastic properties of un-reinforced and Kevlar, glass fibre reinforced
More informationA Multigrid LCR-FEM solver for viscoelastic fluids with application to problems with free surface
A Multigrid LCR-FEM solver for viscoelastic fluids with application to problems with free surface Damanik, H., Mierka, O., Ouazzi, A., Turek, S. (Lausanne, August 2013) Page 1 Motivation Polymer melts:
More informationMathematical Models and Numerical Simulations for the Blood Flow in Large Vessels
Mathematical Models and Numerical Simulations for the Blood Flow in Large Vessels Balazs ALBERT 1 Titus PETRILA 2a Corresponding author 1 Babes-Bolyai University M. Kogalniceanu nr. 1 400084 Cluj-Napoca
More informationChain-configuration and rate dependent rheological properties in transient networks
Electronic Supplementary Material (ESI) for Soft Matter. This journal is The Royal Society of Chemistry 204 Supplementary Information Chain-configuration and rate dependent rheological properties in transient
More informationState Space Solution to the Unsteady Slip Flow of a Micropolar Fluid between Parallel Plates
International Journal of Scientific and Innovative Mathematical Research (IJSIMR) Volume 2, Issue 10, October 2014, PP 827-836 ISSN 2347-307X (Print) & ISSN 2347-3142 (Online) www.arcjournals.org State
More informationPolymer Dynamics and Rheology
Polymer Dynamics and Rheology 1 Polymer Dynamics and Rheology Brownian motion Harmonic Oscillator Damped harmonic oscillator Elastic dumbbell model Boltzmann superposition principle Rubber elasticity and
More informationInstationary Eulerian viscoelastic flow simulations using time separable Rivlin±Sawyers constitutive equations
J. Non-Newtonian Fluid Mech. 89 (2000) 209±228 Instationary Eulerian viscoelastic flow simulations using time separable Rivlin±Sawyers constitutive equations E.A.J.F. Peters, M.A. Hulsen *, B.H.A.A. van
More informationVISCOELASTIC PROPERTIES OF POLYMERS
VISCOELASTIC PROPERTIES OF POLYMERS John D. Ferry Professor of Chemistry University of Wisconsin THIRD EDITION JOHN WILEY & SONS New York Chichester Brisbane Toronto Singapore Contents 1. The Nature of
More information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More informationColloidal Suspension Rheology Chapter 1 Study Questions
Colloidal Suspension Rheology Chapter 1 Study Questions 1. What forces act on a single colloidal particle suspended in a flowing fluid? Discuss the dependence of these forces on particle radius. 2. What
More informationViscoelastic Flows in Abrupt Contraction-Expansions
Viscoelastic Flows in Abrupt Contraction-Expansions I. Fluid Rheology extension. In this note (I of IV) we summarize the rheological properties of the test fluid in shear and The viscoelastic fluid consists
More information5 The Oldroyd-B fluid
5 The Oldroyd-B fluid Last time we started from a microscopic dumbbell with a linear entropic spring, and derived the Oldroyd-B equations: A u = u ρ + u u = σ 2 pi + η u + u 3 + u A A u u A = τ Note that
More informationA formulation for fast computations of rigid particulate flows
Center for Turbulence Research Annual Research Briefs 2001 185 A formulation for fast computations of rigid particulate flows By N. A. Patankar 1. Introduction A formulation is presented for the direct
More informationNonlinear shape evolution of immiscible two-phase interface
Nonlinear shape evolution of immiscible two-phase interface Francesco Capuano 1,2,*, Gennaro Coppola 1, Luigi de Luca 1 1 Dipartimento di Ingegneria Industriale (DII), Università di Napoli Federico II,
More informationNumerical Simulation of a Bubble Rising in an Environment Consisting of Xanthan Gum
Numerical Simulation of a Bubble Rising in an Environment Consisting of Xanthan Gum Víctor, A, Aguirre 1, Byron, A, Castillo 1 and Christian Narvaez-Muñoz 1 a) 1 Departamento de Ciencias de la Energía
More informationSUMMARY A STUDY OF VISCO-ELASTIC NON-NEWTONIAN FLUID FLOWS. where most of body fluids like blood and mucus are non-newtonian ones.
SUMMARY A STUDY OF VISCO-ELASTIC NON-NEWTONIAN FLUID FLOWS Non-Newtonian fluids abound in many aspects of life. They appear in nature, where most of body fluids like blood and mucus are non-newtonian ones.
More informationJournal of Non-Newtonian Fluid Mechanics
J. Non-Newtonian Fluid Mech. 165 (2010) 1047 1054 Contents lists available at ScienceDirect Journal of Non-Newtonian Fluid Mechanics journal homepage: www.elsevier.com/locate/jnnfm Anisotropy parameter
More informationAblaos Asphalt Binder Large Amplitude Oscillatory Shear
Michael J. Farrar mfarrar@uwyo.edu Jean-Pascal Planche Alec Otto Cookman Steve Salmans Ablaos Asphalt Binder Large Amplitude Oscillatory Shear Binder ETG Baton Rouge, LA Sept. 16, 2014 Outline Linear Viscoelasticity
More informationSuspension-based rheological modeling of crystallizing polymer melts
Rheol Acta (28) 47:643 665 DOI 1.17/s397-8-273-4 ORIGINAL CONTRIBUTION Suspension-based rheological modeling of crystallizing polymer melts Rudi J. A. Steenbakkers Gerrit W. M. Peters Received: 1 July
More informationBoundary Conditions in Fluid Mechanics
Boundary Conditions in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University The governing equations for the velocity and pressure fields are partial
More informationMechanical properties of polymers: an overview. Suryasarathi Bose Dept. of Materials Engineering, IISc, Bangalore
Mechanical properties of polymers: an overview Suryasarathi Bose Dept. of Materials Engineering, IISc, Bangalore UGC-NRCM Summer School on Mechanical Property Characterization- June 2012 Overview of polymer
More informationChapter 3: Newtonian Fluid Mechanics. Molecular Forces (contact) this is the tough one. choose a surface through P
// Molecular Constitutive Modeling Begin with a picture (model) of the kind of material that interests you Derive how stress is produced by deformation of that picture Write the stress as a function of
More informationMadrid, 8-9 julio 2013
VI CURSO DE INTRODUCCION A LA REOLOGÍA Madrid, 8-9 julio 2013 NON-LINEAR VISCOELASTICITY Prof. Dr. Críspulo Gallegos Dpto. Ingeniería Química. Universidad de Huelva & Institute of Non-Newtonian Fluid Mechanics
More informationViscoelastic Mechanical Analysis for High Temperature Process of a Soda-Lime Glass Using COMSOL Multiphysics
Viscoelastic Mechanical Analysis for High Temperature Process of a Soda-Lime Glass Using COMSOL Multiphysics R. Carbone 1* 1 Dipartimento di Ingegneria dei Materiali e della Produzione sez. Tecnologie
More informationStudy of rotation of ellipsoidal particles in combined simple shear flow and magnetic fields
Study of rotation of ellipsoidal particles in combined simple shear flow and magnetic fields Jie Zhang 1, Cheng Wang 1 1. Department of Mechanical and Aerospace Engineering, Missouri University of Science
More informationIn-depth analysis of viscoelastic properties thanks to Microrheology: non-contact rheology
In-depth analysis of viscoelastic properties thanks to Microrheology: non-contact rheology Application All domains dealing with soft materials (emulsions, suspensions, gels, foams, polymers, etc ) Objective
More informationAerodynamic force analysis in high Reynolds number flows by Lamb vector integration
Aerodynamic force analysis in high Reynolds number flows by Lamb vector integration Claudio Marongiu, Renato Tognaccini 2 CIRA, Italian Center for Aerospace Research, Capua (CE), Italy E-mail: c.marongiu@cira.it
More information(2.1) Is often expressed using a dimensionless drag coefficient:
1. Introduction Multiphase materials occur in many fields of natural and engineering science, industry, and daily life. Biological materials such as blood or cell suspensions, pharmaceutical or food products,
More informationRheological And Dielectric Characterization of Thermosetting Polymers. Jeffrey Gotro, Ph.D.
Rheological And Dielectric Characterization of Thermosetting Polymers Outline Introduction Oscillatory parallel plate rheometry Dynamic dielectric measurements Simultaneous dynamic mechanical/dielectric
More informationThe two-dimensional streamline upwind scheme for the convection reaction equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2001; 35: 575 591 The two-dimensional streamline upwind scheme for the convection reaction equation Tony W. H. Sheu*,1
More informationSimulations with Dispersive Mixtures Model
Simulations with Dispersive Mixtures Model P.J.A. Janssen, P.D. Anderson & G.W.M. Peters MT03.05 May 12, 2003 1 Introduction The combination of two different materials in a blend can yield interesting
More informationSIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED
1 SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED WITH LATERAL FRICTIONAL EFFECTS High slope results and comparison with experimental data Nathan Martin1, Ioan Ionescu2, Anne Mangeney1,3 François
More informationThe time-dependent extrudate-swell problem of an Oldroyd-B fluid with slip along the wall
The time-dependent extrudate-swell problem of an Oldroyd-B fluid with slip along the wall Eric Brasseur Unité de Mécanique Appliquée, Université Catholique de Louvain, Bâtiment Euler, 4 6 Avenue Georges
More informationMicro and Macro in the Dynamics of Dilute Polymer Solutions
Micro and Macro in the Dynamics of Dilute Polymer Solutions Ravi Prakash Jagadeeshan Complex Fluid Mechanics Closed form equations Constitutive Equations Stress Calculator Simulations (BDS etc) Homogeneous
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationRHEOLASER LAB MICRORHEOLOGY & END USE PROPERTIES ANALYSIS. MICRORHEOLOGY
RHEOLASER LAB & END USE PROPERTIES ANALYSIS A NEW RHEOLOGY APPROACH TO CHARACTERISE END-USE PROPERTIES THE FIRST READY TO USE & END-USE PROPERTIES ANALYSER Rheolaser Rheolaser is the first Lab ready-to-use
More informationOldroyd Viscoelastic Model Lecture Notes
Oldroyd Viscoelastic Model Lecture Notes Drew Wollman Portland State University Maseeh College of Engineering and Computer Science Department of Mechanical and Materials Engineering ME 510: Non-Newtonian
More informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationLaminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid
Laminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid E. Tejaswini 1*, B. Sreenivasulu 2, B. Srinivas 3 1,2,3 Gayatri Vidya Parishad College of Engineering
More informationResponse of an elastic Bingham fluid to oscillatory shear
Rheologica Acta Rheol Acta 26:428-436 (987) Response of an elastic Bingham fluid to oscillatory shear A. S. Yoshimura and R. K. Prud'homme Department of Chemical Engineering, Princeton University Abstract:
More informationNumerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section
Korea-Australia Rheology Journal Vol. 19, No. 2, August 2007 pp. 67-73 Numerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section Mingkan Zhang, Xinrong Shen, Jianfeng
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationA numerical method for steady and nonisothermal viscoelastic fluid flow for high Deborah and Péclet numbers
A numerical method for steady and nonisothermal viscoelastic fluid flow for high Deborah and Péclet numbers Peter Wapperom Martien A. Hulsen (corresponding author) Jaap P.P.M. van der Zanden Delft University
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationAdvantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations
The Hilltop Review Volume 7 Issue 1 Winter 2014 Article 10 December 2014 Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations Tai-Hsien Wu Western Michigan University
More informationIn the name of Allah the most beneficent the most merciful
In the name of Allah the most beneficent the most merciful Transient flows of Maxwell fluid with slip conditions Authors Dr. T. Hayat & Sahrish Zaib Introduction Non-Newtonian Newtonian fluid Newtonian
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationChristel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011
Notes by: Andy Thaler Christel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011 Many complex fluids are shear-thinning. Such a fluid has a shear
More informationOutline. Motivation Governing equations and numerical methods Results: Discussion:
Bifurcation phenomena in strong extensional flows (in a cross-slot geometry) F. A. Cruz 1,*, R. J. Poole 2, F. T. Pinho 3, P.J. Oliveira 4, M. A. Alves 1 1 Departamento de Engenharia Química, CEFT, Faculdade
More informationEntanglements. M < M e. M > M e. Rouse. Zero-shear viscosity vs. M (note change of slope) Edwards degennes Doi. Berry + Fox, slope 3.4.
Entanglements Zero-shear viscosity vs. M (note change of slope) M < M e Rouse slope 3.4 M > M e Edwards degennes Doi slope 1 Berry + Fox, 1968 Question: Which factors affect the Me: T, P, M, flexibility,
More informationPart III. Polymer Dynamics molecular models
Part III. Polymer Dynamics molecular models I. Unentangled polymer dynamics I.1 Diffusion of a small colloidal particle I.2 Diffusion of an unentangled polymer chain II. Entangled polymer dynamics II.1.
More informationViscoelasticity, Creep and Oscillation Experiment. Basic Seminar Applied Rheology
Viscoelasticity, Creep and Oscillation Experiment Basic Seminar Applied Rheology Overview Repetition of some basic terms Viscoelastic behavior Experimental approach to viscoelasticity Creep- and recovery
More informationNon-linear Stress Response of Non-gap-spanning Magnetic Chains Suspended. in a Newtonian Fluid Under Oscillatory Shear Test: A Direct Numerical
Non-linear Stress Response of Non-gap-spanning Magnetic Chains Suspended in a Newtonian Fluid Under Oscillatory Shear Test: A Direct Numerical Simulation M.R. Hashemi, 1 M.T. Manzari, 1, 2, a) and R. Fatehi
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More information