Download
1 / 81
810 likes | 821 Views
A comprehensive overview of viscoelastic flows in polymeric liquids, exploring mathematical models, constitutive laws, and molecular theories. Covers open questions, challenges, and approximations used in research.
E N D
Turning polymeric liquidsinto theorems Michael Renardy Department of Mathematics Virginia Tech Blacksburg, VA 24061-0123, USA
“I came to New York because I heard the streets were paved with gold. • When I got here, I learned three things: • The streets are not paved with gold. • They are not paved at all. • I am expected to pave them.” • Objective: • Why should mathematicians be interested in viscoelastic flows? • What are some of the major open questions? • What challenges arise that are qualitatively different from what everybody • is familiar with? • Mathematical issues in viscoelastic flows: • Existence • Flow stability • The high Weissenberg number limit • Controllability
A quick review of models of non-Newtonian fluids Balance Laws Conservation of mass and momentum: v: velocity T: extra stress tensor p: pressure : density
Constitutive Laws How is the stress tensor related to the motion? Linear models Newtonian fluid: Linear viscoelasticity: G: Stress relaxation modulus (usually assumed completely monotone) Viscosity:
Maxwell model: Take G(s)= exp(- s). Then we have the alternative formulation Similarly, if G is a linear combination of exponentials, we can formulate a differential system. Jeffreys model: T is a linear combination of a Maxwell and Newtonian term. 1/ is called a relaxation time. Weissenberg or Deborah number: Ratio of the relaxation time to a time scale typical of the flow.
Nonlinear models: Generalized Newtonian fluid Viscosity of polymeric fluids in shear flows generally decreases strongly with increasing shear rate. But: Behavior in shear is different from behavior in elongation, and model misses out on normal stresses.
Integral models: Nonlinear generalization of Boltzmann theory Relative deformation gradient: Let y(x,t,s) be the position at time s of the particle that occupies position x at time t. We have Define In general, the stress is a functional of the history of the relative deformation gradient:
Material frame indifference: Stresses result only from deformations, they are unaffected by rotations of the medium. For any rotation Q, we must have Relative Cauchy strain: We have where the functional is isotropic:
K-BKZ integral model: Motivated by analogy with elasticity: Here C is defined analogously as above, but with y(x,t,s) replaced by y(x,t), which is the position of a particle in the equilibrium configuration of the medium.
Differential models A system of ODEs relates the stress to the velocity gradient. Easiest models to analyze or simulate numerically. Examples: Upper convected Maxwell (UCM) model: Phan-Thien Tanner: Giesekus: Johnson-Segalman:
The upper convected Maxwell model has the alternative integral form In general, differential models cannot be transformed to integral form and vice versa.
Molecular theories “The time has come,” the penguin said, “To speak of many things: Of flowing macromolecules, And little beads and springs That join together into ‘chains’ Or even ‘stars’ or ‘rings’.” Robert Byron Bird (with help from Lewis Carroll) • Dilute solution theories • Network theories • Reptation theories.
Polymer molecule is modeled as two spherical beads connected by a spring. • Forces acting on beads: • Spring force. • Drag from surrounding fluid. • Random force due to Brownian motion.
Force balance on each bead: With R=r2-r1, this leads to If the stochastic forces are described by a Wiener process with magnitude kT, we have an equivalent Fokker-Planck equation for a probability density (R,x,t):
Stress tensor: With n denoting the number density of dumbbells, the contribution of the tensions in the springs to the stress is Hookean dumbbells: F(R)=HR. Use the notation Let C is known as a conformation tensor. For Hookean dumbbells, we find
For S=Tp-nkTI, this yields the upper convected Maxwell model: Nonlinear dumbbells: F(R)=(|R|2)R. Peterlin approximation: F(R)=(<|R|2>)R. With the Peterlin approximation, we get
Approximations: • Type A: There is an identified small parameter . If the equations are correct up to • some power of , we can expect that the solutions are also correct up to some power • of . • Type B: The equations are too difficult. But our advisor want us to finish a thesis, • our personnel committee wants us to publish a paper, our dean wants us to get a • grant, etc. We cannot just say we are stuck. So we replace the equations with • simpler ones that might still reproduce essential aspects of the problem. • The Peterlin approximation is Type B. But it can become Type A in two limiting cases: • If is approximately constant. This is a reasonable assumption for R small. • If R is approximately constant. This is a reasonabe assumption for steady • homogeneous flows with large strain rates.
More issues: Dumbbells → chains. Hydrodynamic interaction. Effect of the other polymers in the surrounding fluid.
Types of constitutive theories: • Differential models • Integral models • Microstructure models • Partly solved: • Local existence for initial value problems • Stability of the rest state and global existence for small data • Existence of steady flows for small data • Remaining issues: • Microstructure and integral models are less completely studied than differential • models. • Behavior of distribution function at the limiting extension. • Characterization of decay rates (when they are not exponential). • Conditions at infinity for unbounded flows.
But what about global existence? • Existence proofs for initial value problems have two parts: • An argument for local existence, typically based on proving convergence • of some approximation scheme. • 2. A priori estimates showing solutions do not blow up and can be continued. • The Newtonian case • Assume, for simplicity, periodic boundary conditions. • If we multiply by v and integrate we find
This is enough to guarantee global existence (but not uniqueness) of a weak solution. In two dimensions, we can do more. Take the curl of the equation of motion, and let be the vorticity. We find and hence This suffices to prove global existence of smooth solutions.
Do non-Newtonian fluids help? • It was once widely believed that non-Newtonian fluids might have “better” existence • results. This is the case only in very limited instances, e.g. • Shear thickening generalized Newtonian fluids. • The “thermodynamic” second order fluid: • Great theorems if 1>0 and 2=21. But 1<0 in real fluids. • 3. Equations based on “cutting corners” (e.g. keeping the inertial nonlinearities • and neglecting the constitutive ones). • In general, however, the situation for non-Newtonian fluids is far less understood than • the Newtonian case. The reason for this is the lack of a priori estimates.
The corotational Jeffreys fluid (Lions and Masmoudi) Energy estimate: Multiply momentum equation by v, and constitutive equation by T/(2) and integrate. The result is
Lions and Masmoudi prove that this a priori bound can be used to prove global existence of weak solutions. But this only works because of a miracle: Compare the upper convected Maxwell model: Note that The correct energy estimate is but this is much too weak to infer existence of weak solutions.
Why is this an a priori estimate at all? Note that It follows that positive definiteness of T+I is preserved. Only solutions with positive definite stresses are physically relevant. The Challenge: The only global existence results of any kind for realistic models of polymeric fluids are for one-dimensional shear flows and they are based on treating the viscoelastic terms as a “controlled” perturbations of a Newtonian problem.
A remark on the side If you feel nervous about drinking a glass of water while global well-posedness of Navier-Stokes is unresolved … you might consider that the water must form a free surface jet before it gets into the glass. Free surface jets break into droplets. Asymptotic solutions for breakup have been studied for Newtonian as well as viscoelastic fluids. Heuristic arguments and numerical evidence suggest there should be general theorems linking the initial shape of the jet to the eventual asymptotics of the breakup (or, if no breakup occurs, the asymptotics for t). No such theorems have been proved, even for one-dimensional approximations of the Newtonian case.
Stability problems in viscoelastic flows • Viscoelastic flows have been found to show many instabilities, even at zero • Reynolds number. • Examples: • Shear flows with curved streamlines. • Entry flows. • Extrudate instabilities. • Elastic mechanisms for jet breakup. • Interfacial instabilities.
“Usual” analysis of instabilities • Linearize about a known base state. • Look for eigenvalues crossing the imaginary axis at critical values of the • relevant parameter (e.g. Reynolds or Weissenberg number) • Do a bifurcation analysis. • All of this has been justified rigorously for the Navier-Stokes equations. • However, there are major unresolved problems for viscoelastic flows.
Abstract framework: C0 semigroups Consider an evolution problem with solution For matrices, we know that the eigenvalues of exp(At) correspond to those of A, but in infinite dimensions the situation is more complicated. In general, we can split the spectrum of A and of exp(At) into isolated eigenvalues of finite multiplicity and the “rest,” known as the essential spectrum. The isolated eigenvalues cause no problem, but we cannot always infer the essential spectrum of exp(At) from that of A. It is actually possible for exp(At) to have an essential spectrum even if A does not. But does this happen in “real” problems?
Example: with 2-periodic boundary conditions. The natural setup is to set and look for (u,v) in the function space H1 L2. The eigenvalues of A are purely imaginary, and there is no essential spectrum. However, the essential spectrum of exp(At) has radius exp(t/2). Newtonian flows: Essential spectra exist only in unbounded geometries. Even then, linear stability of the flow can be inferred from the spectrum of A. Viscoelastic flows: Essential spectra always exist, they are notoriously difficult to compute, and we have no rigorous connection between stability and spectrum.
Advective equations (R. Shvydkoy) These are problems of the form where q satisfies periodic boundary conditions and A is a pseudodifferential operator of order zero (something like, for instance the inverse of a second order differential operator times another second order differential operator). Shvydkoy’s result does not say stability is determined by the spectrum of but instead it characterizes what else we need to determine stability.
Pseudodifferential operators of order zero have the property that where A0 is homogeneous of order zero: For instance, we could have
To determine the stability of an advective equation, we have to do two things: • Find the discrete eigenvalues. • Study the stability of the amplitude b in the following ODE system: • The good news: Creeping flows of fluids with differential constitutive laws fit into • this framework. • The bad news: At this point, this only works for periodic problems.
The high Weissenberg number limit Analogy: high Reynolds number limit for Newtonian flow Problem: pour a pint of beer We ignore the difficult free surface flow and concentrate on the flow inside the tap, which we idealize as a straight pipe.
Attempt 1: Poiseuille flow solution. This is a solution at any Reynolds number, but is not observed at high Reynolds number. Attempt 2: Euler equations. Allow for any parallel flow profile, and numerous other solutions. Too many solutions to predict anything. Attempt 3: Assume approximately uniform flow away from the walls, and try to match the boundary conditions with boundary layers. Textbook example of boundary layers: Boundary layers in fluids are “sort of” like this example, but only “sort of.” There is the well known Blasius solution. But how does a boundary layer that looks like this fit into a tap that looks like that?
The actual description of this seemingly simple problem is a patchwork of heuristics, formal approximation, truncated models etc. We really cannot expect the situation to be any better in the high Weissenberg number limit of viscoelastic flows. Theorems will necessarily be limited in scope, addressing only specific aspects of the overall problem. Excerpt from careerplanner.com: Skill Requirements for: "Bartender" Active Listening -- Giving full attention to what other people are saying, taking time to understand the points being made, asking questions as appropriate, and not interrupting at inappropriate times.Speaking -- Talking to others to convey information effectively.Social Perceptiveness -- Being aware of others' reactions and understanding why they react as they do.Mathematics -- Using mathematics to solve problems.
High Weissenberg Number Asymptotics The upper convected Maxwell model and the Euler equations Steady creeping flow of the UCM: If we formally set W=, we get
If we assume T is of rank 1, i.e. T=qqT, where is normalized in such a way that then we obtain the equations The first two equations are like the compressible Euler equations with q replacing the velocity and replacing the density. However, there is no “equation of state” linking p to .
Analogue of potential flow Assume a two-dimensional flow where q is parallel to v, and is constant along streamlines. Then we obtain precisely the incompressible Euler equations. All solutions of these equations, in particular potential flows, reappear as solutions of the high Weissenberg number viscoelastic problem.
Boundary layers: Model problem Suppose f is periodic in x. If we formally let W → , u becomes equal to the average of f, but for y=0, we must have u=f(x). There is a boundary layer when y is of order 1/W. In fact, the solution is In viscoelastic fluids, similar boundary layers arise in the stress, due to the presence of the term
Boundary layer for the UCM fluid Governing equations:
The boundary is at y=0. If we assume a shear rate of order 1 at the wall, then T12 is of order 1, and T11 is of order W there. If we balance, for instance, the term u( T11/ x) against T11/W, we find that the term becomes significant when u (i.e. y) is of order 1/W. This motivates the scaling
