Chauvi ere and Lozinski as deterministic Fokker-Planck (FP) equations. The first option gives rise to stochastic numerical methods [6], [7], which hav
|
|
- Beverly Hood
- 5 years ago
- Views:
Transcription
1 Simulation of dilute polymer solutions using a Fokker-Planck equation Cédric Chauvi ere a and Alexei Lozinski b;λ a Division of Applied Mathematics, Brown University, Box F, Providence, RI 9, USA. b FSTI - ISE - LMF, Ecole Polytechnique Fédérale de Lausanne, CH 5 Lausanne, Switzerland. Abstract The main goal of this work is to rehabilitate the use of the Fokker-Planck equation for simulating complex flows of non-newtonian fluids which do not possess closed-form constitutive equations. Anumerical method based on the Fokker-Planck equation is applied to the simulation of a dilute solution of polymer modeled by FENE dumbbells. The proposed method is tested for both homogeneous planar extensional flow and the flow around a cylinder in a channel. Introduction Non-Newtonian fluids which have a memory of past deformation are called viscoelastic fluids and are present in a wide range of applications e.g. multigrade oils, food processing, biological fluids such as blood. In this paper we shall be concerned with the numerical simulation of dilute polymer solutions which may be regarded as being a low concentration of polymer chains in a Newtonian solvent. Thus, we can assume that the polymer chains do not interact with each other. The solvent is assumed to be incompressible and isothermal and we restrict ourselves to the case of inertialess flows (zero Reynolds number). Therefore, the mass and momentum equations have the form of the Stokes equations in which the extra-stress fi, i.e. the polymeric part of the Cauchy stress tensor, is given as a source term: = rp + s r u + r fi ; () r u =; () where u, p and s denote, respectively, the velocity, pressure and the viscosity of the Newtonian solvent. The extra-stress usually comes from a supplementary equation (closed-form constitutive equation) or a mesoscopic model. Over the past twenty years, numerous numerical methods have been designed for improving the discretization of constitutive equations for viscoelastic flows. It would be too lengthy to enumerate them here and the interested reader is referred to the excellent review in [9]. Mesoscopic models are generally regarded as being (potentially) more realistic than their closedform counterparts but their numerical simulation requires much more work. Mathematically, these models can be written in two equivalent forms (see [7]): either as stochastic differential equations or Λ corresponding author, alexei.lozinski@epfl.ch
2 Chauvi ere and Lozinski as deterministic Fokker-Planck (FP) equations. The first option gives rise to stochastic numerical methods [6], [7], which have become very popular during the last years. In contrast, simulations of dilute polymer solutions through the use of an FP equation are very rare in the literature and have been limited, so far, to very simple flows. Since the pioneering work of Warner [] in 97, followed thirteen years later by Fan [4] for the simulation of simple shear flows for FENE (finitely extensible non-linear elastic) dumbbells, no attempt has been made to use the FP equation for the simulation of more complex flows using this model. It is the purpose of this paper to show that using an FP equation may be a good alternative to stochastic computations for the simulations of dilute polymer solutions, even when complex flows are considered. The paper is organized as follows. In the next section, we recall the equations used to model dilute polymer liquids (the FENE model) and in Section 3 we describe how these equations may be discretized. The computation of the polymeric extra-stress from the solution of our FP equation is then detailed in Section 4. In Section 5 we illustrate our method with two numerical examples: a simple planar extensional flow and the flow past a cylinder in a channel. Finally, in Section 6 we draw some conclusions about the potential of the FP equation to simulate complex viscoelastic flows. Modeling dilute polymer solutions The idea for the modeling of dilute polymeric solutions is to replace the polymer chains by dumbbells which consist of two beads connected by a spring. A dumbbell whose centre of mass is at x at time t is fully described by a configuration end-to-end vector q(t; x). Let ψ(t; x; q) be the probability density function (pdf) of the random process q(t; x). In other words ψ(t; x; q)dq denotes the probability that a dumbbell at position x and at time t is to be found with the configuration vector in the box [q; q + dq]. The pdf satisfies the Fokker-Planck equation, which can be written (see [], [7]) as Dψ Dt +div q (»q F(q))ψ = qψ; (3) where D=Dt is the material derivative,» is the transposed velocity gradient tensor (» ij i =@x j ), is the relaxation time of the fluid and F(q) is the spring force, which takes the following form for FENE dumbbells: q F(q) = jqj b : (4) In the relation above, the parameter b controls the maximum extensibility of the dumbbells which cannot exceed p b. The computation of the velocity field with ()-() is decoupled from the computation of the pdf, so in equation (3), we can assume that the velocity gradient isknown. Since the vector q lies in a disc of radius p b, it is natural to represent it in polar coordinates, i.e q = r cos ; q = r sin ; with r [; p b] and [; ß]: (5) Introducing two functions b ( ) and b ( ) of the angle as follows b ( ) =» cos +» +» b ( ) =» sin +» +» sin ; (6) cos +»» ; (7)
3 Simulation of dilute polymer solutions 3 a detailed expression of equation (3) may be written as Dψ Dt = rb b + br b r + b (b r ) ψ ψ : (8) Equation (8) is completed with the initial condition at time t =, which corresponds to the equilibrium solution (» = ) and takes the following form: ψ eq = b + ßb In the next section, we will show how equation (8) may be discretized. 3 Discretization of the equations b= r : (9) In order to satisfy the boundary conditions at r = and r = p b, we shall use instead of ψ(t; x;r; ) a new unknown ff(t; x; ; ) defined by s ψ(t; x;r; )= ff(t; x; ; ); () where the relation between r and [ ; ] is r = b + and s is a fixed, strictly positive real number. To prevent the spring force (4) from becoming infinite for r = p b,we should have ψ =at r = p b (or = ) and this is achieved by introducing the multiplier (( )=) s into (). We have found experimentally that our numerical method has the best stability properties under s =. We have used this value of s to obtain the results of Section 5. Note that the pdf should be symmetric in q, thatexplains why we have chosen the new variable to be a function of r. In order to decouple the computation in the physical space Ω from the computation in the configuration space, we use the following time-splitting for equation (8) b ff i+= ff i = L FP (ff i+= ); t () ff i+ ff i+= +(u i+ r)ff i+ =; t () where L FP refers to the Fokker-Planck operator appearing on the right-hand side of (8) after the change of variable (). Thus we have reduced the task of solving a four-dimensional problem for ff(t; x; ; ) at a given time t to solving two-dimensional problems for ff i+= ( ; ) and ff i+ (x) respectively. Each step of the splitting ()-() requires a specific type of discretization to account for the type of the equation and the space in which the solution is sought. For equation (), we search for an approximate solution ff( ; ) of the form ff( ; ) = XNR NF i= X j= ff ijh i ( ) cos(j )+ XNR NF i= X j= ff ijh i ( ) sin(j ); (3) where fh i ( )g»i»nr are Lagrange interpolating polynomials based on the Gauss-Legendre points i (see the book of Canuto et al. [] for more details). Note that the set f i g is chosen so that it does not include the points = and = since the boundary conditions there are already
4 4 Chauvi ere and Lozinski taken into account by (). Only the Fourier modes of even order are kept in (3) because of the symmetry of ff( ; ). Inserting (3) into (), we form the product of () with a test function fh k ( ) cos(l )g»k»n R»l»NF or fh k( ) sin(l )g»k»n R and integrate over configuration space. The integrals with respect to are evaluated using Gauss-Lobatto quadrature rule whereas the integrals with»l»nf respect to can be computed analytically. The resulting linear system for the unknown coefficients ffij and ffij is LU decomposed and a back substitution gives the numerical solution of equation (). We note that this spectral method conserves the integral of the pdf over the configuration space provided s in () is an integer from to N R. An SUPG element-by-element spectral element method with constant upwinding factor, as detailed by Chauvi ere and Owens [3], is used to discretize the hyperbolic equation () in Ω. Periodic boundary conditions are applied at inflow andoutflow for the flow around a cylinder in a channel. The Stokes equations ()-() are discretized using a spectral element method, as described in the paper by Owens and Phillips [8], for example. 4 Computation of the extra-stress Having computed the pdf at x at time t, the extra-stress which will serve as a source term for the Stokes equations ()-() is given by what is known as the Kramers expression [7]. For the FENE model and a two dimensional problem it is fi (x;t) = p b + ( I+ < q Ω F(q) >) (4) b where the symbol Ω denotes the tensor product of two vectors and p is the polymeric viscosity of the fluid. The brackets < > in (4) denote the statistical average over configuration space, i.e. < q Ω F(q) >= Z D(; p b) q Ω F(q)ψ(t; x; q)dq: (5) Using (4),(5) and (), the evaluation of (4) gives the following expressions for the components of the extra-stress! fi xx = p b + XNR ψ + ßb! i ( + i ) s i (ff i + ff b 6 i) ; (6) fi xy = p fi yy = p b + ßb b 6 b + where! i are the quadrature weights. 5 Numerical results b XNR i= ψ + ßb 6 5. Planar extensional flow i= s! i ( + i ) i ff i; (7) XNR i= s! i ( + i ) i (ff i ff i)! ; (8) As a first test, we consider a homogeneous planar extensional flow for which the velocity field is given by (u x ;u y )=( "x; "y), and " is the extension rate. One can find an analytical expression of the pdf
5 Simulation of dilute polymer solutions 5 for the steady-state solution of (3) in this case (see []). It has the form b= ψ = C jqj exp( "(q b q )) (9) where C is some normalization constant. We report the results of two computations for different values of the extension rate (" = and " = 5), the other parameters being p = = and b =. Table gives the computed values of the xx-component of the extra-stress and the relative error for different resolutions in the case " =. We can see that the error is very small even for low resolutions. At this extension rate, the solution is very smooth (see Figure, which represents the pdf for the computational mesh (N F ;N R )=(6; 3)) and only afewpoints are needed for our spectral method to properly capture the solution. When the extension rate is increased up to " = 5, getting accurate numerical results requires much more refined meshes, as is shown in Table. This is because the exact solution at high extension rates possesses huge gradients which are very localized. Capturing such a solution requires highly refined meshes. Figure shows the numerical oscillatory solution for a low level of discretization (N F ;N R )=(6; 3) and Figure 3 shows how the solution can be properly captured when the mesh is refined up to (N F ;N R )=(; 4): In the next section we will apply the same numerical method to a complex flow. 5. Flow in complex geometry The problem of planar viscoelastic flow around a cylinder (see Figure 4) is used to test our method for more complex geometries. This benchmark problem is rheologically rich since it offers both shear flows (near solid walls) and extensional flows in the wake of the cylinder. We choose the aspect ratio Λ=R=H ==, where H is the width of the channel and R is the radius of the cylinder. We set the ratio s =( s + p ) equal to.59, and b equal to. The time steps are chosen equal to t =:5. The flow domain is divided into 6 spectral elements and polynomials of degree N =6are used in the two spatial direction for the representation of the variables (with exception of the pressure). The pressure basis polynomials are chosen to be of degree less in each direction so as to ensure that the discrete problem is well posed (see []). The iterations are stopped when the following convergence criterion is fulfilled for all collocation points x Ω : ku(x;t i+ ) u(x;t i )k» 4 : () t For this problem, we may define the Deborah number as De = U R ; () where U is the average velocity at entry or exit. For the FENE model, this benchmark problem has only been investigated with stochastic simulations of Monte-Carlo type [5]. We fix the Deborah number equal to.8, obtained by taking ( ; U;R)=(:8; ; ) and we increase the resolution of the configuration space from (N F ;N R )=(4; 9) up to (N F ;N R )=(6; 3). The convergence of our numerical scheme with mesh refinement is then established by plotting the component fi xx of the extra-stress on the axis of symmetry (jxj > ) and on the cylinder surface (jxj < ). This is shown in Figure 5, where for the two finest meshes the results match very well. We also show the contours of fi xx over the whole physical domain in Figure 6, to demonstrate the smoothness of the plot. Our results are obtained on a personal computer with the CPU Pentium III at 8MHz. The CPU time per time step was.3, 6.4 and 4.5 seconds for the meshes mentioned at Figure 5. time steps were needed to reach the steady state.
6 6 Chauvi ere and Lozinski 6 Conclusion In this paper, we have shown that computing complex viscoelastic flows using an FP equation is now possible. The key ingredient to make such a simulation possible is the time-splitting which reduces the computation of a four-dimensional problem to the computation of two-dimensional problems. The advantage of this approach over stochastic simulations is that the resulting solution is not polluted by statistical noise. The method introduced in this paper is very general and can also be applied to models describing concentrated polymer solutions such as reptation models []. Acknowledgment. The authors would like to thank Robert Owens for helpful discussions and for proofreading the manuscript. The work of the second author has been supported by the Swiss National Science Foundation, grant number References [] Bird RB, Armstrong RC, Hassager O. Dynamics of polymeric liquids, vol., Kinetic theory, nd ed. New York, Chichester, Brisbane, Toronto, Singapore: Wiley-Interscience, 987. [] Canuto C, Hussaini MY, Quarteroni A, Zang TA. Spectral Methods in Fluid Dynamics. Berlin, Heidelberg: Springer Verlag, 988. [3] Chauvi ere C, Owens RG. A new spectral element method for the reliable computation of viscoelastic flow. Comput. Methods Appl. Mech. Engrg. ; 9: [4] Fan XJ. Viscosity, first normal-stress coefficient, and molecular stretching in dilute polymer solutions. J. Non-Newtonian Fluid Mech. 985; 7:5-44. [5] van Heel APG, Hulsen MA, van den Brule BHAA. On the selection of parameters in the FENE-P model. J. Non-Newtonian Fluid Mech. 998; 75:53-7. [6] Hulsen MA, van Heel APG, van den Brule BHAA. Simulation of viscoelastic flows using Brownian configuration fields. J. Non-Newtonian Fluid Mech. 997; 7:79-. [7] Öttinger HC. Stochastic Processes in Polymeric Fluids. Berlin: Springer Verlag, 996. [8] Owens RG, Phillips TN. Steady viscoelastic flow past a sphere using spectral elements. Int. J. Numer. Meth. Engrg. 996; 39: [9] Owens RG, Phillips TN. Computational Rheology. Singapore: Imperial College Press/World Scientific,. [] Warner HR. Kinetic theory and rheology of dilute suspensions of finitely extensible dumbbells. Ind. Eng. Chem. Fundam. 97; :
7 Simulation of dilute polymer solutions 7 List of Figures Plot of the pdf ψ for " = and (N F ;N R )=(6; 3): Plot of the pdf ψ for " = 5 and (N F ;N R )=(6; 3): Plot of the pdf ψ for " = 5 and (N F ;N R )=(; 4): Cylinder radius R placed symmetrically in a D channel of half width H Profile of the component fi xx on the cylinder surface and in the wake of the cylinder for three levels of discretization Contour of fi xx at De =:8 for FENE model List of Tables Numerical results for planar extensional flow at " = (the exact value of fi xx is ) Numerical results for planar extensional flow at " = 5 (the exact value of fi xx is 4: )
8 8 Chauvi ere and Lozinski Figure : Plot of the pdf ψ for " = and (N F ;N R )=(6; 3): Figure : Plot of the pdf ψ for " = 5 and (N F ;N R )=(6; 3):
9 Simulation of dilute polymer solutions Figure 3: Plot of the pdf ψ for " =5and (N F ;N R ) = (; 4): H R Figure 4: Cylinder radius R placed symmetrically in a D channel of half width H.
10 Chauvi ere and Lozinski 9 8 (N F,N R )=(6,3) (N F,N R )=(5,) (N F,N R )=(4,9) 7 Txx x Figure 5: Profile of the component fi xx on the cylinder surface and in the wake of the cylinder for three levels of discretization. Figure 6: Contour of fi xx at De =:8 for FENE model.
11 Simulation of dilute polymer solutions N R N F Computed Relative error 5 8: : : : : : : : : :6 8: :5 Table : Numerical results for planar extensional flow at " = (the exact value of fi xx is ). N R N F Computed Relative error 5 diverging : :35 : : 3 5 4: : : : 7 Table : Numerical results for planar extensional flow at " = 5 (the exact value of fi xx is 4: )
Chapter 1 Introduction
Chapter 1 Introduction This thesis is concerned with the behaviour of polymers in flow. Both polymers in solutions and polymer melts will be discussed. The field of research that studies the flow behaviour
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 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 informationExistence of global weak solutions to implicitly constituted kinetic models of incompressible homogeneous dilute polymers
1 / 31 Existence of global weak solutions to implicitly constituted kinetic models of incompressible homogeneous dilute polymers Endre Süli Mathematical Institute, University of Oxford joint work with
More informationExplaining and modelling the rheology of polymeric fluids with the kinetic theory
Explaining and modelling the rheology of polymeric fluids with the kinetic theory Dmitry Shogin University of Stavanger The National IOR Centre of Norway IOR Norway 2016 Workshop April 25, 2016 Overview
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 informationEvents of High Polymer Activity in Drag Reducing Flows
Flow Turbulence Combust (27) 79:23 32 DOI.7/s494-7-975-5 Events of High Polymer Activity in Drag Reducing Flows C. Tesauro B. J. Boersma M. A. Hulsen P. K. Ptasinski F. T. M. Nieuwstadt Received: May 26
More informationA Numerical Study of Several Viscoelastic Fluid Models
A Numerical Study of Several Viscoelastic Fluid Models Corina Putinar January 3, 6 Abstract Viscoelastic fluids are a type of fluid with a solvent and immersed elastic filaments which create additional
More informationA MAXIMUM ENTROPY PRINCIPLE BASED CLOSURE METHOD FOR MACRO-MICRO MODELS OF POLYMERIC MATERIALS. Yunkyong Hyon. Jose A. Carrillo. Qiang Du.
Kinetic and Related Models Website: http://aimsciences.org c American Institute of Mathematical Sciences Volume 1, Number, June 8 pp. 171 18 A MAXIMUM ENTROPY PRINCIPLE BASED CLOSURE METHOD FOR MACRO-MICRO
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 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 informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationThe Polymers Tug Back
Tugging at Polymers in Turbulent Flow The Polymers Tug Back Jean-Luc Thiffeault http://plasma.ap.columbia.edu/ jeanluc Department of Applied Physics and Applied Mathematics Columbia University Tugging
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 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 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 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 informationModeling of a Floc-Forming Fluid Using a Reversible Network Model
ORIGINAL PAPER Journal of Textile Engineering (2014), Vol.60, No.4, 61-67 2014 The Textile Machinery Society of Japan Modeling of a Floc-Forming Fluid Using a Reversible Network Model Ya m a m o t o Takehiro
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 informationSystematic Closure Approximations for Multiscale Simulations
Systematic Closure Approximations for Multiscale Simulations Qiang Du Department of Mathematics/Materials Sciences Penn State University http://www.math.psu.edu/qdu Joint work with C. Liu, Y. Hyon and
More informationarxiv: v1 [math.ap] 10 Jul 2017
Existence of global weak solutions to the kinetic Peterlin model P. Gwiazda, M. Lukáčová-Medviďová, H. Mizerová, A. Świerczewska-Gwiazda arxiv:177.2783v1 [math.ap 1 Jul 217 May 13, 218 Abstract We consider
More informationThe FENE-L and FENE-LS closure approximations to the kinetic theory of finitely extensible dumbbells
J. Non-Newtonian Fluid Mech. 87 (1999) 179 196 The FENE-L and FENE-LS closure approximations to the kinetic theory of finitely extensible dumbbells G. Lielens, R. Keunings,V.Legat CESAME, Division of Applied
More informationOn the Boundary Layer Flow of a Shear Thinning Liquid over a 2-Dimensional Stretching Surface
Advanced Studies in Theoretical Physics Vol. 12, 2018, no. 1, 25-36 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/astp.2018.71259 On the Boundary Layer Flow of a Shear Thinning Liquid over a 2-Dimensional
More informationCoil-stretch transition and the breakdown of computations for viscoelastic fluid flow around a confined cylinder
Coil-stretch transition and the breakdown of computations for viscoelastic fluid flow around a confined cylinder Mohit Bajaj Department of Chemical Engineering, Monash University, Melbourne, Australia
More informationvs. Chapter 4: Standard Flows Chapter 4: Standard Flows for Rheology shear elongation 2/1/2016 CM4650 Lectures 1-3: Intro, Mathematical Review
CM465 Lectures -3: Intro, Mathematical //6 Chapter 4: Standard Flows CM465 Polymer Rheology Michigan Tech Newtonian fluids: vs. non-newtonian fluids: How can we investigate non-newtonian behavior? CONSTANT
More informationNavier-Stokes Flow in Cylindrical Elastic Tubes
Navier-Stokes Flow in Cylindrical Elastic Tubes Taha Sochi University College London, Department of Physics & stronomy, Gower Street, London, WC1E 6BT Email: t.sochi@ucl.ac.uk. bstract nalytical expressions
More informationDeficiencies of FENE dumbbell models in describing the rapid stretching of dilute polymer solutions
Deficiencies of FENE dumbbell models in describing the rapid stretching of dilute polymer solutions Indranil Ghosh, Gareth H. McKinley, a) Robert A. Brown, and Robert C. Armstrong b) Department of Chemical
More informationA spectral vanishing viscosity method for stabilizing viscoelastic flows
J. Non-Newtonian Fluid Mech. 115 (3) 15 155 A spectral vanishing viscosity method for stabilizing viscoelastic flows X. Ma, V. Symeonidis, G.E. Karniadakis Division of Applied Mathematics, Center for Fluid
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 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 informationA MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationSmoothed Dissipative Particle Dynamics: theory and applications to complex fluids
2015 DPD Workshop September 21-23, 2015, Shanghai University Smoothed Dissipative Particle Dynamics: Dynamics theory and applications to complex fluids Marco Ellero Zienkiewicz Centre for Computational
More informationFinal Report: DE-FG02-95ER25239 Spectral Representations of Uncertainty: Algorithms and Applications
Final Report: DE-FG02-95ER25239 Spectral Representations of Uncertainty: Algorithms and Applications PI: George Em Karniadakis Division of Applied Mathematics, Brown University April 25, 2005 1 Objectives
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
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 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 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 informationDIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER
Conference on Turbulence and Interactions TI26, May 29 - June 2, 26, Porquerolles, France DIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER E. Leriche, Laboratoire d Ingénierie
More informationKinetic Theory for Rigid Dumbbells
Kinetic Theory for Rigid Dumbbells Hector D. Ceniceros hdc@math.ucsb.edu University of California @ Santa Barbara Mathematics Department University of California, Santa Barbara Universidade de São Paulo,
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationMICRO-MACRO METHODS FOR THE MULTISCALE SIMULATION OF VISCOELASTIC FLOW USING MOLECULAR MODELS OF KINETIC THEORY
to appear in Rheology Reviews 2004, D.M. Binding and K. Walters (Eds.), British Society of Rheology, November 2003 MICRO-MACRO METHODS FOR THE MULTISCALE SIMULATION OF VISCOELASTIC FLOW USING MOLECULAR
More informationNumerical study of the flow of FENE-CR viscoelastic fluids in a planar cross-slot geometry
Numerical study of the flow of FENE-CR viscoelastic fluids in a planar cross-slot geometry G.N. Rocha 1, R.J. Poole 2, M.A. Alves 3 and P.J. Oliveira 1 1 Universidade da Beira Interior, Departamento de
More informationNumerical modelling of shear-thinning non-newtonian flows in compliant vessels
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2007; 00:1 [Version: 2002/09/18 v1.01] Numerical modelling of shear-thinning non-newtonian flows in compliant vessels M.
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
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 informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationON THE EFFECT OF CORNER CURVATURE IN LAMINAR FLOW THROUGH A CROSS-SLOT CHANNEL
XXIX CILAMCE Iberian Latin American Congress on Computational Methods in Engineering November 4-7, 2008, Maceió-Alagoas, Brazil ON THE EFFECT OF CORNER CURVATURE IN LAMINAR FLOW THROUGH A CROSS-SLOT CHANNEL
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationTWO-DIMENSIONAL MAGMA FLOW *
Iranian Journal of Science & Technology, Transaction A, Vol. 34, No. A2 Printed in the Islamic Republic of Iran, 2010 Shiraz University TWO-DIMENSIONAL MAGMA FLOW * A. MEHMOOD 1** AND A. ALI 2 1 Department
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES A conservative scheme for the Fokker-Planck equation with applications to viscoelastic polymeric fluids Hana Mizerová Bangwei She Preprint No. 34-2018
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 informationSOME THOUGHTS ON DIFFERENTIAL VISCOELASTIC MODELS AND THEIR NUMERICAL SOLUTION
CONFERENCE ON COMPLEX FLOWS OF COMPLEX FLUIDS University of Liverpool, UK, March 17-19, 2008 SOME THOUGHTS ON DIFFERENTIAL VISCOELASTIC MODELS AND THEIR NUMERICAL SOLUTION Paulo J. Oliveira Universidade
More informationAnalysis of the Oseen-Viscoelastic Fluid Flow Problem
Analysis of the Oseen-Viscoelastic Fluid Flow Problem Vincent J. Ervin Hyesuk K. Lee Louis N. Ntasin Department of Mathematical Sciences Clemson University Clemson, South Carolina 9634-0975 December 30,
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 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 informationA new numerical approach for Soret effect on mixed convective boundary layer flow of a nanofluid over vertical frustum of a cone
Inter national Journal of Pure and Applied Mathematics Volume 113 No. 8 2017, 73 81 ISSN: 1311-8080 printed version); ISSN: 1314-3395 on-line version) url: http://www.ijpam.eu ijpam.eu A new numerical
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 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 informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationComputer Fluid Dynamics E181107
Computer Fluid Dynamics E181107 2181106 Transport equations, Navier Stokes equations Remark: foils with black background could be skipped, they are aimed to the more advanced courses Rudolf Žitný, Ústav
More informationOn the oscillatory tube flow of healthy human blood
On the oscillatory tube flow of healthy human blood Miguel A. Moyers-Gonzalez Department of Mathematics and Statistics, University of Canterbury, New Zealand Robert G. Owens Département de mathématiques
More informationMultiscale Modelling of Complex Fluids: A Mathematical Initiation
Multiscale Modelling of Complex Fluids: A Mathematical Initiation Claude Le Bris 1,2 1,2 and Tony Lelièvre 1 CERMICS, École Nationale des Ponts et Chaussées, 6 & 8, avenue Blaise Pascal, F-77455 Marne-La-Vallée
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 information[Supplementary Figures]
[Supplementary Figures] Supplementary Figure 1 Fabrication of epoxy microchannels. (a) PDMS replica is generated from SU-8 master via soft lithography. (b) PDMS master is peeled away from PDMS replica
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationPaper CIT UNSTEADY VISCOELASTIC FLOW PAST A CYLINDER IN A CHANNEL
Proceedings of the 10o Brazilian Congress of Thermal Sciences and Engineering -- ENCIT 004 Braz. Soc. of Mechanical Sciences and Engineering -- ABCM, Rio de Janeiro, Brazil, Nov. 9 -- ec. 03, 004 Paper
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
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 informationSimulations with adaptive modeling
Simulations with adaptive modeling Alexandre Ern CERMICS Ecole des Ponts, ParisTech INdAM Workshop, Cortona, 21/09/06 Introduction Engineering problems often involve models of different complexity, scales
More information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
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 information2009 Best Paper Understanding and Quantification of Die Drool Phenomenon During Polypropylene Extrusion Process
2009 Best Paper Understanding and Quantification of Die Drool Phenomenon During Polypropylene Extrusion Process Print (10)» 2010 Best Paper An Engineering Approach to the Correction of Rotational Flow
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
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 information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationQUIZ 2 OPEN QUIZ WHEN TOLD THERE ARE TWO PROBLEMS OF EQUAL WEIGHT. Please answer each question in a SEPARATE book
2.341J MACROMOLECULAR HYDRODYNAMICS Spring 2012 QUIZ 2 OPEN QUIZ WHEN TOLD THERE ARE TWO PROBLEMS OF EQUAL WEIGHT Please answer each question in a SEPARATE book You may use the course textbook (DPL) and
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 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 informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationhttp://www.springer.com/3-540-30725-7 Erratum Spectral Methods Fundamentals in Single Domains C. Canuto M.Y. Hussaini A. Quarteroni T.A. Zang Springer-Verlag Berlin Heidelberg 2006 Due to a technical error
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 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 informationNumerical analysis for the BCF method in complex fluids simulations
for the BCF method in complex fluids simulations Tiejun Li School of Mathematical Sciences, Peking University, Beijing 100871 tieli@pku.edu.cn Joint work with Weinan E and Pingwen Zhang CSCAMM conference,
More informationVisualization of polymer relaxation in viscoelastic turbulent micro-channel flow
Supplementary Information for Visualization of polymer relaxation in viscoelastic turbulent micro-channel flow Authors: J. Tai, C. P. Lim, Y. C. Lam Correspondence to: MYClam@ntu.edu.sg This document includes:
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 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 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 informationRheology and Constitutive Equations. Rheology = Greek verb to flow. Rheology is the study of the flow and deformation of materials.
Rheology and Constitutive Equations Rheology = Greek verb to flow Rheology is the study of the flow and deformation of materials. The focus of rheology is primarily on the study of fundamental, or constitutive,
More informationUnsteady Boundary Layer Flow and Symmetry Analysis of a Carreau Fluid
Mathematics and Statistics: Open Access Received: Dec, 5, Accepted: Jan 5, 6, Published: Jan 8, 6 Math Stat, Volume, Issue http://crescopublications.org/pdf/msoa/msoa--.pdf Article Number: MSOA-- Research
More informationRheology of a dilute suspension of spheres in a viscoelastic fluid under LAOS
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,
More informationOn the Hysteretic Behaviour of Dilute Polymer Solutions in Relaxation Following Extensional Flow
On the Hysteretic Behaviour of Dilute Polymer Solutions in Relaxation Following Extensional Flow R. Sizaire, G. Lielens, I. Jaumain, R. Keunings and V. Legat CESAME, Division of Applied Mechanics, Université
More informationMODELING THREE-DIMENSIONAL NON-NEWTONIAN FLOWS IN SINGLE- SCREW EXTRUDERS
MODELING THREE-DIMENSIONAL NON-NEWTONIAN FLOWS IN SINGLE- SCREW EXTRUDERS Christian Marschik Wolfgang Roland Bernhard Löw-Baselli Jürgen Miethlinger Johannes Kepler University Institute of Polymer Extrusion
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More informationRevisit to Grad s Closure and Development of Physically Motivated Closure for Phenomenological High-Order Moment Model
Revisit to Grad s Closure and Development of Physically Motivated Closure for Phenomenological High-Order Moment Model R. S. Myong a and S. P. Nagdewe a a Dept. of Mechanical and Aerospace Engineering
More informationCESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 CESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP Yiolanda Damianou, George Kaoullas,
More information