notes viscosity
play

Notes Viscosity Thursday: class will begin late, 11:30am In - PowerPoint PPT Presentation

Notes Viscosity Thursday: class will begin late, 11:30am In reality, nearby molecules travelling at different velocities occasionally bump into each other, transferring energy Differences in velocity reduced (damping) Measure this


  1. Notes Viscosity � Thursday: class will begin late, 11:30am � In reality, nearby molecules travelling at different velocities occasionally bump into each other, transferring energy • Differences in velocity reduced (damping) • Measure this by strain rate (time derivative of strain, or how far velocity field is from rigid motion) • Add terms to our constitutive law cs533d-term1-2005 1 cs533d-term1-2005 2 Strain rate Viscous stress � At any instant in time, measure how fast chunk � As with linear elasticity, end up with two parameters if we want isotropy: of material is deforming from its current state viscous = 2 µ ˙ • Not from its original state � ij � ij + � ˙ � kk � ij • So we � re looking at infinitesimal, incremental strain • µ and � are coefficients of viscosity (first and second) updates • These are not the Lame coefficients! Just use the same symbols • Can use linear Cauchy strain! • � damps only compression/expansion � So the strain rate tensor is � Usually � � -2/3 µ (exact for monatomic gases) � So end up with � � + � u j � u i � ij = 1 � � ˙ � � � � + � u j viscous = µ � u i � u k � x j � x i 2 � 2 � � � ij � � ij � � � � x j � x i � x k � 3 � cs533d-term1-2005 3 cs533d-term1-2005 4

  2. Navier-Stokes Nondimensionalization � Navier-Stokes equations include the viscous � Actually go even further stress � Select a characteristic length L � Incompressible version: ( ) • e.g. the width of the domain, u t + u � � u + 1 � � p = g + 1 � � � µ � u + � u T � And a typical velocity U � � u = 0 • e.g. the speed of the incoming flow � Often (but not always) viscosity µ is constant, � Rescale terms and this reduces to • x � =x/L, u � =u/U, t � =tU/L, p � =p/ � U 2 µ u t + u � � u + 1 � � p = g + � � 2 u so they all are dimensionless U 2 + � • Call � = µ / � the “kinematic viscosity” u ' t + u ' �� u ' + � p ' = Lg UL � 2 u ' cs533d-term1-2005 5 cs533d-term1-2005 6 Nondimensional parameters Why nondimensionalize? � Re=UL/ � is the Reynold � s number � Think of it as a user-interface issue � It lets you focus on what parameters matter • The smaller it is, the more viscosity plays a • If you scale your problem so nondimensional role in the flow parameters stay constant, solution scales • High Reynold � s numbers are hard to simulate � Code rot --- you may start off with code which has true dimensions, but as you hack around � Fr= is the Froude number U g L they lose meaning • The smaller it is, the more gravity plays a role • If you � re nondimensionalized, there should be only in the flow one or two parameters to play with • Note: often can ignore gravity (pressure � Not always the way to go --- you can look up material constants, but not e.g. Reynolds gradient cancels it out) numbers u t + u � � u + � p = (0, � 1,0) + 1 Re � 2 u Fr 2 cs533d-term1-2005 7 cs533d-term1-2005 8

  3. Other quantities Boundary conditions � Inviscid flow: � We may want to carry around auxiliary • Solid wall: u•n=0 quantities • Free surface: p=0 (or atmospheric pressure) • E.g. temperature, the type of fluid (if we have • Moving solid wall: u•n=u wall •n a mix), concentration of smoke, etc. � Also enforced in-flow/out-flow • Between two fluids: u 1 •n=u 2 •n and p 1 =p 2 + �� � Use material derivative as before � Viscous flow: � E.g. if quantity doesn � t change, just is • No-slip wall: u=0 transported (“advected”) around: • Other boundaries can involve coupling tangential Dq components of stress tensor… Dt = q t + u � � q � = 0 � � � Pressure/velocity coupling at boundary: • u•n modified by � p/ � n advection cs533d-term1-2005 9 cs533d-term1-2005 10 What now? Vorticity � Can numerically solve the full equations � How do we measure rotation? � = � � u • Will do this later • Vorticity of a vector field (velocity) is: • Not so simple, could be expensive (3D) • Proportional (but not equal) to angular velocity of a rigid body - off by a factor of 2 � Or make assumptions and simplify them, � Vorticity is what makes smoke look then solve numerically interesting • Simplify flow (e.g. irrotational) • Turbulence • Simplify dimensionality (e.g. go to 2D) cs533d-term1-2005 11 cs533d-term1-2005 12

  4. Vorticity equation Vorticity equation continued � Start with N-S, constant viscosity and density � Simplify with more vector identities, and assume u t + u � � u + 1 � � p = g + � � 2 u incompressible to get: D � � Take curl of whole equation Dt = � � � u + � � 2 � ( ) ( ) + � � 1 � � u t + � � u � � u � � p = � � g + � � � � 2 u � Important result: Kelvin Circulation Theorem � Lots of terms are zero: • Roughly speaking: if � =0 initially, and there � s no • g is constant (or the potential of some field) viscosisty, � =0 forever after (following a chunk of • With constant density, pressure term too fluid) ( ) = � � � � 2 u � � u t + � � u � � u � If fluid starts off irrotational, it will stay that way � Then use vector identities to simplify… (in many circumstances) ) = � � 2 � � u ( ( ) � u + 1 ( ) � � u t + � � � � u 2 � u 2 ( ) = � � 2 � � t + � � � � u cs533d-term1-2005 13 cs533d-term1-2005 14 Potential flow Potential in time � If velocity is irrotational: � What if we have a free surface? � � u = 0 � Use vector identity u• � u=( �� u) � u+ � |u| 2 /2 � Assume • Which it often is in simple laminar flow • incompressible ( � •u=0), inviscid, irrotational ( �� u=0) � Then there must be a stream function (potential) such that: • constant density u = � � • thus potential flow (u= �� , � 2 � =0) � Solve for incompressibility: � Then momentum equation simplifies (using G=-gy for gravitational potential) � � � � = 0 2 + 1 ( ) � u + 1 u t + � � u 2 � u � � p = g � If u•n is known at boundary, we can solve it 2 + 1 � � t + 1 2 � u � � p = �� G cs533d-term1-2005 15 cs533d-term1-2005 16

  5. Bernoulli’s equation Water waves � Every term in the simplified momentum � For small waves (no breaking), can reduce geometry of water to 2D heightfield equation is a gradient: integrate to get 2 u 2 + p � Can reduce the physics to 2D also � t + 1 � = � G • How do surface waves propagate? � Plan of attack • (Remember Bernoulli � s law for pressure?) • Start with potential flow, Bernoulli � s equation � This tells us how the potential should • Write down boundary conditions at water surface evolve in time • Simplify 3D structure to 2D cs533d-term1-2005 17 cs533d-term1-2005 18 Set up Boundaries � y = 0 � We � ll take y=0 as the height of the water at � At sea floor (y=-H), v=0 rest � At sea surface (y=h � 0), v=h t � H is the depth (y=-H is the sea bottom) • Note again - assuming very small horizontal motion � h is the current height of the water at (x,z) � y = h t � Simplification: velocities very small (small amplitude waves) � At sea surface (y=h � 0), p=0 • Or atmospheric pressure, but we only care about pressure differences • Use Bernoulli � s equation, throw out small velocity squared term, use p=0, � t = � gh cs533d-term1-2005 19 cs533d-term1-2005 20

  6. Finding a wave solution Dispersion relation � Plug in � =f(y)sin(K•(x,z)- � t) � So the wave speed c is c = � • In other words, do a Fourier analysis in g k = k tanh kH horizontal component, assume nothing much happens in vertical � Notice that waves of different wave- • Solving � 2 � =0 with boundary conditions on � y numbers k have different speeds ( ) gives cosh K ( y + H ) � = A � • Separate or disperse in time ( ) sin K � ( x , z ) � � t ( ) K sinh K H � For deep water (H big, k reasonable -- not tsunamis) tanh(kH) � 1 • Pressure boundary condition then gives (with k=|K|, the magnitude of K) g c = � = gk tanh kH k cs533d-term1-2005 21 cs533d-term1-2005 22 Simulating the ocean Energy spectrum � Far from land, a reasonable thing to do is � Fourier decomposition of height field: � • Do Fourier decomposition of initial surface ( ) e ( ) � x , z ( ) ˆ h ( x , z , t ) = � 1 i , j h i , j , t height i , j 2 • Evolve each wave according to given wave S ( K ) = ˆ � “Energy” in K=(i,j) is h ( K ) speed (dispersion relation) � Update phase, use FFT to evaluate � Oceanographic measurements have found � How do we get the initial spectrum? models for expected value of S(K) • Measure it! (oceanography) (statistical description) cs533d-term1-2005 23 cs533d-term1-2005 24

Recommend


More recommend