Microstructure and Particle-Phase Stress in a Dense Suspension
|
|
- Patience Sims
- 5 years ago
- Views:
Transcription
1 Microstructure and Particle-Phase Stress in a Dense Suspension Jim Jenkins, Cornell University Luigi La Ragione, Politecnico di Bari Predict microstructure and particle stresses in slow, steady, pure straining of a dense, single layer of spheres in a viscous fluid Singh & Nott, J Fluid Mech. 412, 279 (2) Jenkins & La Ragione, J. Fluid Mech. 763, 218 (215)
2 Pure Straining θ B d (BA) A BA) Area fraction: ν Average number of near contacts: k Short-range repulsive force: F
3 Elements of the Analysis Pair translations and rotations Pair forces and moments Force and moment balances Averages over directions and neighbors for equilibrium Rate equations for the separation and orientation Distribution of particle orientation and particle flux Averaging over particle orientations for the stress
4 Kinematics Treat BA exactly Relative velocity of centers of pair BA v α (BA) = d dt s(ba) ˆdα (BA) + 2a d dt θ(ba)ˆt α (BA) Relative velocity of points of near contact of BA v (BA) α + a( ω (A) + ω (B) (BA) )ˆt α Treat na on average Relative velocity of centers of pair na v i (na) = 2aD αβ ˆdβ (na) Relative velocity of points of near contact of na v α (na) + aω (A)ˆt α (na)
5 Interaction Forces Jeffrey & Onishi (1984), Jeffrey (1992) Pair BA F (BA) α = 6πµaK (BA) αβ v (BA) β F ˆd (BA) s (BA) α 9.54πµa 2 ( ˆt (BA) (BA) (BA) β D βγ ˆd γ )ˆt α +πµa 2 ln K (BA) = 1 αβ 4 a s (BA).96 ω(a)ˆt (BA) α + πµa 2 ln a ˆd (BA) s (BA) α ˆd (BA) β ln a s (BA) a s (BA) ˆt (BA) α (BA)ˆt β ω(b)ˆt (BA) α Pairs na F (na) α = 12πµa 2 K (na) αβ D βγ ˆd (na) γ + πµa 2 F s ˆd α (na) + 2πµa 2 K (na) αβ = 1 a 4 ln a s.96 ω(a)ˆt α (na) ln 1 s.96 ( ˆt (na) (na) (na) β D βγ ˆd γ )ˆt α s ˆd (na) α ˆd (na) β ln a s ˆt (na) α (na)ˆt β
6 Force Balance, Particle A N ( A ) n B F α (BA) + F α (na) = F α (AB) = F α (BA), ˆd α (AB) = ˆd α (BA), ˆt α (AB) = ˆt α (BA) = 4πµa 2 ln +πµa 2 Difference of force balances = 3πµa a s (BA) 2ln +6πµa 2 a s a s (BA) Normal component d dt s(ba) 2 F s (BA) ˆd (BA) α J αβγ D βγ 2 F s Y (BA) ˆd α α Tangential component a s (BA) d dt θ(ba) 19πµa 2ˆt (BA) (BA) α D ˆdβ αβ.96 S + 6πµa a 2 s ˆt (BA) α J αβγ D βγ (S = ω A + ω B )
7 Moment Balance, Particle A N ( A ) n B ε αβ F α (BA) ˆdβ (BA) + ε αβ F α (na) ˆdβ (na) = (ε 12 = - ε 12 = 1, ε 11 = ε 22 = ) = 2 ln + ln a s (BA) a s (BA) Sum of Moment Balances d dt θ(ba) 9.54ˆt (BA) (BA) α D ˆdβ αβ + k 1 2 ln a s S + 2 ln a s.96 ε A D βγ αβ αγ
8 Neighborhood Averages Isotropic representations Y α N ( A ) n B ˆdα (na) = ξˆd α (BA) A αβ N ( A ) n B ˆdα (na) ˆd β (na) = β 1 δ αβ + β 2 ˆdα (BA) ˆdβ (BA) J αβγ N ( A ) n B ˆdα (na) ˆdβ (na) ˆd γ (na) (BA) (BA) = α ˆdβ 1ˆdα ˆd (BA) γ + α ( 2 δ ˆdk (BA) αβ + δ ˆdβ (BA) αγ + δ ) (BA) ˆdα βγ
9 Neighborhood Averages s a = π 8νg g = 16 7ν 16(1 ν) 2 Isotropic distribution function π Ψ α..γ = 2 I( θ )ˆd α..ˆd γ d θ I( θ ) =, < θ < π / 3 3(k 1) / 4π, π / 3 < θ < π α 1 =, α 2 = 3 3(k 1) 16π α β 1 = (k 1) (k 1) 16π, β 2 = 3 3(k 1) 8π = 2α ξ = 3 3(k 1) 4π = 4α
10 Trajectory Equations!γ D αβ =!γ ˆd (BA) α D ˆdβ (BA) αβ =!γ cos2θ, ˆd α (BA) = (sinθ,cosθ) ˆt α (BA) D αβ ˆdβ (BA) =!γ sin2θ Lengths dimensionless by a and F F / ( πµa 3!γ ) 1 s ds dγ = 2 3 F 1 s + 4 α s + 4 α s cos2θ ln 1 s dθ dγ = α 1 s sin2θ 1 2 ln 1 s S Tangential force difference - Sum moments S = 12α (4α k +1) 1/ s ln(1/ s).96 sin2θ
11 Contact Distribution A(θ)dθ: number of contacts within dθ A(θ) dθ dγ = constant and A(θ) π/2 dθ = k 4 Knowledge of k determines the constant of integration. Implement as differential equations I(θ) A( θ ) d θ θ di dθ = A, I() =, I(π / 2) = k / 4 da dθ = A d θ! θ! dθ d θ! dθ =! θ θ +! θ s ds dθ
12 Differential Equations ds dθ = ds / dγ dθ / dγ f[s,θ; s(ν),k] da dθ = A dθ / dγ θ dθ dγ + s dθ dγ ds dθ di dθ = A dγ dθ = 1 dθ / dγ Boundary Conditions x 1 θ = θ + (θ 1 θ )x s(θ ) =.3 s(θ 1 ) =.3 I(θ ) = I(θ 1 ) = k / 4 γ(θ ) = γ(θ 1 ) = 1.1
13 Separation \ k = 3., γ = 1.1, F = 1-4 ν =.55 ν =.6 ν =.65 2 s θ The repulsive force creates asymmetry about π/4.
14 Inverse Separation k = 3., γ = 1.1, F = 1-4 ν =.55 ν =.6 ν = /s θ The asymmetry of the separation clarified.
15 Angular Velocity 5 k = 3., γ = 1.1, F = dθ/dγ ν =.55 ν =.6 ν = θ A(θ) dθ dγ = constant
16 Circumferential Distribution 4 35 k = 3., γ = 1.1, s = s 1 =.3 A, F = A+1, F = θ The repulsive force shifts the anisotropy due to the flow, symmetric about π/4, to a distribution that is asymmetric.
17 Circumferential Distribution k = 3., γ = 1.1, F = 1-4 ν =.55 ν =.6 ν =.65 A θ Asymmetry about π/4 in both the distribution and interval.
18 T αβ = na Stress (dimensional) 2π A(θ)F α (BA) ˆdβ (BA) dθ n: number/area F α (BA) = N ( A ) (na) F α n B t αβ = T αβ aµ!γ, n = ν πa 2, F = F πa 3 µ!γ t αβ = ν α 2π A(θ) ( 4F + 6cos2θ) s ˆd ˆdβ α dθ +ν s α 2π A(θ)sin2θ ˆt ˆdβ α dθ +ν 12α (3α k +1) 1/ s ln(1/ s) 2π A(θ)ln 1 s sin2θ ˆt ˆdβ α dθ ˆd α (AB) = (sinθ,cosθ), ˆt α (AB) = (cosθ, sinθ)
19 Pressure p 1 2 (t xx + t yy ) p = 4ν α s π/2 ( ) A(θ) 2F + 3cos2θ dθ Shear stress τ 1 2 (t xx t yy ) π/2 ( ) τ = 4ν α A(θ) 2F + 3cos2θ cos2θdθ s +ν α 2π s A(θ)sin 2 2θdθ +ν 6α (3α k +1) 1/ s ln(1/ s) 2π A(θ)ln 1 s sin2 2θdθ
20 Stresses 25 k = 3., γ = 1.1, F = 1-4 Pressure (o), Shear Stress (x) ν Asymmetry about π/4 creates particle pressure.
21 Initial/Final Separation 2 k = 3., γ = 1.1, F = 1-4 Pressure (o), Shear Stress (x) s, s 1 The initial/final separation of the typical trajectory is where the shear stress exhibits a maximum (.3).
22 Conclusions A rough implementation of particle equilibrium for particles assumed to interact as pairs through points of near contact produces expressions for the normal and tangential components of their relative velocities as functions of their orientation with respect to the flow. From these velocities, the separation of neighboring particles and their circumferential distribution can be determined as functions of orientation. The existence of asymmetries in the circumferential distribution of particles influences the shear stress and produces a pressure associated with viscous interactions between particle pairs at points of near contact.
Chapter 6: Momentum Analysis
6-1 Introduction 6-2Newton s Law and Conservation of Momentum 6-3 Choosing a Control Volume 6-4 Forces Acting on a Control Volume 6-5Linear Momentum Equation 6-6 Angular Momentum 6-7 The Second Law of
More informationThe response of an idealized granular material to an incremental shear strain
The response of an idealized granular material to an incremental shear strain Luigi La Ragione and Vanessa Magnanimo Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Bari, Italy E-mail: l.laragione@poliba.it;
More informationChapter 6: Momentum Analysis of Flow Systems
Chapter 6: Momentum Analysis of Flow Systems Introduction Fluid flow problems can be analyzed using one of three basic approaches: differential, experimental, and integral (or control volume). In Chap.
More information(r i F i ) F i = 0. C O = i=1
Notes on Side #3 ThemomentaboutapointObyaforceF that acts at a point P is defined by M O (r P r O F, where r P r O is the vector pointing from point O to point P. If forces F, F, F 3,..., F N act on particles
More informationdt Now we will look at the E&M force on moving charges to explore the momentum conservation law in E&M.
. Momentum Conservation.. Momentum in mechanics In classical mechanics p = m v and nd Newton s law d p F = dt If m is constant with time d v F = m = m a dt Now we will look at the &M force on moving charges
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 informationClusters in granular flows : elements of a non-local rheology
CEA-Saclay/SPEC; Group Instabilities and Turbulence Clusters in granular flows : elements of a non-local rheology Olivier Dauchot together with many contributors: D. Bonamy, E. Bertin, S. Deboeuf, B. Andreotti,
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 informationMHA042 - Material mechanics: Duggafrågor
MHA042 - Material mechanics: Duggafrågor 1) For a static uniaxial bar problem at isothermal (Θ const.) conditions, state principle of energy conservation (first law of thermodynamics). On the basis of
More informationOscillatory Motion. Solutions of Selected Problems
Chapter 15 Oscillatory Motion. Solutions of Selected Problems 15.1 Problem 15.18 (In the text book) A block-spring system oscillates with an amplitude of 3.50 cm. If the spring constant is 250 N/m and
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationI. Elastic collisions of 2 particles II. Relate ψ and θ III. Relate ψ and ζ IV. Kinematics of elastic collisions
I. Elastic collisions of particles II. Relate ψ and θ III. Relate ψ and ζ IV. Kinematics of elastic collisions 49 I. Elastic collisions of particles "Elastic": KE is conserved (as well as E tot and momentum
More informationChapter 1. Viscosity and the stress (momentum flux) tensor
Chapter 1. Viscosity and the stress (momentum flux) tensor Viscosity and the Mechanisms of Momentum Transport 1.1 Newton s law of viscosity ( molecular momentum transport) 1.2 Generalization of Newton
More informationName (please print) π cos(θ) + sin(θ)dθ
Mathematics 2443-3 Final Eamination Form B December 2, 27 Instructions: Give brief, clear answers. I. Evaluate by changing to polar coordinates: 2 + y 2 3 and above the -ais. + y d 23 3 )/3. π 3 Name please
More informationand in each case give the range of values of x for which the expansion is valid.
α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω Mathematics is indeed dangerous in that it absorbs students to such a degree that it dulls their senses to everything else P Kraft Further Maths A (MFPD)
More informationName (please print) π cos(θ) + sin(θ)dθ
Mathematics 2443-3 Final Eamination Form A December 2, 27 Instructions: Give brief, clear answers. I. Evaluate by changing to polar coordinates: 2 + y 2 2 and above the -ais. + y d 2(2 2 )/3. π 2 (r cos(θ)
More informationCalculus 152 Take Home Test 2 (50 points)
Calculus 5 Take Home Test (5 points) Due Tuesday th November. The following test will be done at home in order to ensure that it is a fair and representative reflection of your own ability in mathematics
More informationWater is sloshing back and forth between two infinite vertical walls separated by a distance L: h(x,t) Water L
ME9a. SOLUTIONS. Nov., 29. Due Nov. 7 PROBLEM 2 Water is sloshing back and forth between two infinite vertical walls separated by a distance L: y Surface Water L h(x,t x Tank The flow is assumed to be
More informationUniformity of the Universe
Outline Universe is homogenous and isotropic Spacetime metrics Friedmann-Walker-Robertson metric Number of numbers needed to specify a physical quantity. Energy-momentum tensor Energy-momentum tensor of
More informationFinal Examination F.5 Mathematics M2 Suggested Answers
Final Eamination F.5 Mathematics M Suggested Answers. The (r + )-th term C 9 r ( ) 9 r r 9 C r r 7 7r For the 8 term, set 7 7r 8 r 5 Coefficient of 8 C 9 5 5. d 8 ( ) set d if > slightly, d we have
More informationVector Geometry Final Exam Review
Vector Geometry Final Exam Review Problem 1. Find the center and the radius for the sphere x + 4x 3 + y + z 4y 3 that the center and the radius of a sphere z 7 = 0. Note: Recall x + ax + y + by + z = d
More informationDEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD R. A. Amaro
More informationHandout 6: Rotational motion and moment of inertia. Angular velocity and angular acceleration
1 Handout 6: Rotational motion and moment of inertia Angular velocity and angular acceleration In Figure 1, a particle b is rotating about an axis along a circular path with radius r. The radius sweeps
More informationInterfacial dynamics
Interfacial dynamics Interfacial dynamics = dynamic processes at fluid interfaces upon their deformation Interfacial rheological properties: elasticity, viscosity, yield stress, Relation between macroscopic
More informationCOLLECTIVE HYDRODYNAMICS OF SWIMMING MICROORGANISMS
COLLECTIVE HYDRODYNAMICS OF SWIMMING MICROORGANISMS Timothy J Pedley DAMTP, University of Cambridge APS-DFD, Salt Lake City, 2007 OUTLINE 1 Bioconvection Mechanisms Model Instability of a uniform suspension
More information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
More informationDynamics 12e. Copyright 2010 Pearson Education South Asia Pte Ltd. Chapter 20 3D Kinematics of a Rigid Body
Engineering Mechanics: Dynamics 12e Chapter 20 3D Kinematics of a Rigid Body Chapter Objectives Kinematics of a body subjected to rotation about a fixed axis and general plane motion. Relative-motion analysis
More informationComposites Design and Analysis. Stress Strain Relationship
Composites Design and Analysis Stress Strain Relationship Composite design and analysis Laminate Theory Manufacturing Methods Materials Composite Materials Design / Analysis Engineer Design Guidelines
More informationStress, Strain, Mohr s Circle
Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected
More informationDiscontinuous Shear Thickening
Discontinuous Shear Thickening dynamic jamming transition Ryohei Seto, Romain Mari, Jeffrey F. Morris, Morton M. Denn Levich Institute, City College of New York First experimental data Williamson and Hecker
More informationVortex Motion and Soliton
International Meeting on Perspectives of Soliton Physics 16-17 Feb., 2007, University of Tokyo Vortex Motion and Soliton Yoshi Kimura Graduate School of Mathematics Nagoya University collaboration with
More informationThe curious properties of spin
The curious properties of spin February, 07 The Stern-Gerlach experiment The Schrödinger equation predicts degenerate energy levels for atoms electron states that differ only in the z-component of their
More informationPROBLEM Copyright McGraw-Hill Education. Permission required for reproduction or display. SOLUTION. ω = 29.6 rad/s. ω = = 36 3.
PROLEM 15.1 The brake drum is attached to a larger flywheel that is not shown. The motion of the brake drum is defined by the relation θ = 36t 1.6 t, where θ is expressed in radians and t in seconds. Determine
More informationEnergy during a burst of deceleration
Problem 1. Energy during a burst of deceleration A particle of charge e moves at constant velocity, βc, for t < 0. During the short time interval, 0 < t < t its velocity remains in the same direction but
More informationIII.H Zeroth Order Hydrodynamics
III.H Zeroth Order Hydrodynamics As a first approximation, we shall assume that in local equilibrium, the density f 1 at each point in space can be represented as in eq.(iii.56), i.e. [ ] p m q, t)) f
More informationPhysicsAndMathsTutor.com
1. A uniform circular disc has mass 4m, centre O and radius 4a. The line POQ is a diameter of the disc. A circular hole of radius a is made in the disc with the centre of the hole at the point R on PQ
More informationDifferential Equations: Homework 8
Differential Equations: Homework 8 Alvin Lin January 08 - May 08 Section.6 Exercise Find a general solution to the differential equation using the method of variation of parameters. y + y = tan(t) r +
More informationASSIGNMENT COVER SHEET omicron
ASSIGNMENT COVER SHEET omicron Name Question Done Backpack Ready for test Drill A differentiation Drill B sketches Drill C Partial fractions Drill D integration Drill E differentiation Section A integration
More informationA Two-Phase Continuum Theory for Windblown Sand
A Two-Phase Continuum Theory for Windblown Sand Jim Jenkins Civil and Environmental Engineering Cornell University Alexandre Valance Institut de Physique de Rennes Université de Rennes 1 Suspension and
More informationHOMEWORK SOLUTIONS MATH 1910 Sections 8.2, 8.3, 8.5 Fall 2016
HOMEWORK SOLUTIONS MATH 191 Sections 8., 8., 8.5 Fall 16 Problem 8..19 Evaluate using methods similar to those that apply to integral tan m xsec n x. cot x SOLUTION. Using the reduction formula for cot
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationS12.1 SOLUTIONS TO PROBLEMS 12 (ODD NUMBERS)
OLUTION TO PROBLEM 2 (ODD NUMBER) 2. The electric field is E = φ = 2xi + 2y j and at (2, ) E = 4i + 2j. Thus E = 2 5 and its direction is 2i + j. At ( 3, 2), φ = 6i + 4 j. Thus the direction of most rapid
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationChapter 22. Dr. Armen Kocharian. Gauss s Law Lecture 4
Chapter 22 Dr. Armen Kocharian Gauss s Law Lecture 4 Field Due to a Plane of Charge E must be perpendicular to the plane and must have the same magnitude at all points equidistant from the plane Choose
More informationMA3D1 Fluid Dynamics Support Class 5 - Shear Flows and Blunt Bodies
MA3D1 Fluid Dynamics Support Class 5 - Shear Flows and Blunt Bodies 13th February 2015 Jorge Lindley email: J.V.M.Lindley@warwick.ac.uk 1 2D Flows - Shear flows Example 1. Flow over an inclined plane A
More informationPart 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA
Part 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA Review of Fundamentals displacement-strain relation stress-strain relation balance of momentum (deformation) (constitutive equation) (Newton's Law)
More informationExaminer: D. Burbulla. Aids permitted: Formula Sheet, and Casio FX-991 or Sharp EL-520 calculator.
University of Toronto Faculty of Applied Science and Engineering Solutions to Final Examination, June 017 Duration: and 1/ hrs First Year - CHE, CIV, CPE, ELE, ENG, IND, LME, MEC, MMS MAT187H1F - Calculus
More informationStrain Transformation equations
Strain Transformation equations R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur-613 401 Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 8 Table of Contents 1. Stress transformation
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationThe viscosity of a dilute suspension of rough spheres
J. Fluid Mech. (000), vol. 4, pp. 339 367. Printed in the United Kingdom c 000 Cambridge University Press 339 The viscosity of a dilute suspension of rough spheres By HELEN J. WILSON AND ROBERT H. DAVIS
More informationSurface stress-induced self-assembly
Surface stress-induced self-assembly Zhigang Suo Work with Wei Lu Solid surface Liquid surface Stretched solid surface Energy per surface atom depends on stretch Stretched liquid surface Energy per surface
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 information1.1. Rotational Kinematics Description Of Motion Of A Rotating Body
PHY 19- PHYSICS III 1. Moment Of Inertia 1.1. Rotational Kinematics Description Of Motion Of A Rotating Body 1.1.1. Linear Kinematics Consider the case of linear kinematics; it concerns the description
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationShear stress of a monolayer of rough spheres
J. Fluid Mech. (2002), vol. 452, pp. 425 441. c 2002 Cambridge University Press DOI: 10.1017/S002211200100688 Printed in the United Kingdom 425 Shear stress of a monolayer of rough spheres By HELEN J.
More informationMath Calculus II Homework # Due Date Solutions
Math 35 - Calculus II Homework # - 007.08.3 Due Date - 007.09.07 Solutions Part : Problems from sections 7.3 and 7.4. Section 7.3: 9. + d We will use the substitution cot(θ, d csc (θ. This gives + + cot
More information2 Basic Equations in Generalized Plane Strain
Boundary integral equations for plane orthotropic bodies and exterior regions G. Szeidl and J. Dudra University of Miskolc, Department of Mechanics 3515 Miskolc-Egyetemváros, Hungary Abstract Assuming
More information1 4 (1 cos(4θ))dθ = θ 4 sin(4θ)
M48M Final Exam Solutions, December 9, 5 ) A polar curve Let C be the portion of the cloverleaf curve r = sin(θ) that lies in the first quadrant a) Draw a rough sketch of C This looks like one quarter
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationSolutions Ph 236b Week 1
Solutions Ph 236b Week 1 Page 1 of 7 Solutions Ph 236b Week 1 Kevin Barkett and Mark Scheel January 19, 216 Contents Problem 1................................... 2 Part (a...................................
More informationx + ye z2 + ze y2, y + xe z2 + ze x2, z and where T is the
1.(8pts) Find F ds where F = x + ye z + ze y, y + xe z + ze x, z and where T is the T surface in the pictures. (The two pictures are two views of the same surface.) The boundary of T is the unit circle
More information2.6 Oseen s improvement for slow flow past a cylinder
Lecture Notes on Fluid Dynamics.63J/.J) by Chiang C. Mei, MIT -6oseen.tex [ef] Lamb : Hydrodynamics.6 Oseen s improvement for slow flow past a cylinder.6. Oseen s criticism of Stokes approximation Is Stokes
More informationMonte Carlo Radiation Transfer I
Monte Carlo Radiation Transfer I Monte Carlo Photons and interactions Sampling from probability distributions Optical depths, isotropic emission, scattering Monte Carlo Basics Emit energy packet, hereafter
More informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
More informationAssignment 6 Solution. Please do not copy and paste my answer. You will get similar questions but with different numbers!
Assignment 6 Solution Please do not copy and paste my answer. You will get similar questions but with different numbers! This question tests you the following points: Integration by Parts: Let u = x, dv
More informationNonlinear stability of steady flow of Giesekus viscoelastic fluid
Nonlinear stability of steady flow of Giesekus viscoelastic fluid Mark Dostalík, V. Průša, K. Tůma August 9, 2018 Faculty of Mathematics and Physics, Charles University Table of contents 1. Introduction
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More information1. Evaluate the integrals. a. (9 pts) x e x/2 dx. Solution: Using integration by parts, let u = x du = dx and dv = e x/2 dx v = 2e x/2.
MATH 8 Test -SOLUTIONS Spring 4. Evaluate the integrals. a. (9 pts) e / Solution: Using integration y parts, let u = du = and dv = e / v = e /. Then e / = e / e / e / = e / + e / = e / 4e / + c MATH 8
More information7a3 2. (c) πa 3 (d) πa 3 (e) πa3
1.(6pts) Find the integral x, y, z d S where H is the part of the upper hemisphere of H x 2 + y 2 + z 2 = a 2 above the plane z = a and the normal points up. ( 2 π ) Useful Facts: cos = 1 and ds = ±a sin
More informationHW6. 1. Book problems 8.5, 8.6, 8.9, 8.23, 8.31
HW6 1. Book problems 8.5, 8.6, 8.9, 8.3, 8.31. Add an equal strength sink and a source separated by a small distance, dx, and take the limit of dx approaching zero to obtain the following equations for
More informationSuggested Homework Solutions
Suggested Homework Solutions Chapter Fourteen Section #9: Real and Imaginary parts of /z: z = x + iy = x + iy x iy ( ) x iy = x #9: Real and Imaginary parts of ln z: + i ( y ) ln z = ln(re iθ ) = ln r
More informationMassachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004
Massachusetts Institute of Technology Department of Physics Course: 8.09 Classical Mechanics Term: Fall 004 Final Examination December 17, 004 Instructions Do not start until you are told to do so. Solve
More informationROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring
ROTATIONAL STRESSES INTRODUCTION High centrifugal forces are developed in machine components rotating at a high angular speed of the order of 100 to 500 revolutions per second (rps). High centrifugal force
More informationPhysics 141 Rotational Motion 2 Page 1. Rotational Motion 2
Physics 141 Rotational Motion 2 Page 1 Rotational Motion 2 Right handers, go over there, left handers over here. The rest of you, come with me.! Yogi Berra Torque Motion of a rigid body, like motion of
More informationPotential flow of a second-order fluid over a sphere or an ellipse
J. Fluid Mech. (4), vol. 5, pp. 5. c 4 Cambridge niversity Press DOI:.7/S4954 Printed in the nited Kingdom Potential flow of a second-order fluid over a sphere or an ellipse By J. W A N G AND D. D. J O
More informationTime-Dependent Conduction :
Time-Dependent Conduction : The Lumped Capacitance Method Chapter Five Sections 5.1 thru 5.3 Transient Conduction A heat transfer process for which the temperature varies with time, as well as location
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationa x Questions on Classical Solutions 1. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress
Questions on Classical Solutions. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress σ xy = T is applied at infinity. Determine the value of the stress σ θθ on the edge
More informationThermocapillary Migration of a Drop
Thermocapillary Migration of a Drop An Exact Solution with Newtonian Interfacial Rheology and Stretching/Shrinkage of Interfacial Area Elements for Small Marangoni Numbers R. BALASUBRAMANIAM a AND R. SHANKAR
More informationElastic Scattering. R = m 1r 1 + m 2 r 2 m 1 + m 2. is the center of mass which is known to move with a constant velocity (see previous lectures):
Elastic Scattering In this section we will consider a problem of scattering of two particles obeying Newtonian mechanics. The problem of scattering can be viewed as a truncated version of dynamic problem
More informationUnit 13 Review of Simple Beam Theory
MIT - 16.0 Fall, 00 Unit 13 Review of Simple Beam Theory Readings: Review Unified Engineering notes on Beam Theory BMP 3.8, 3.9, 3.10 T & G 10-15 Paul A. Lagace, Ph.D. Professor of Aeronautics & Astronautics
More informationHölder regularity estimation by Hart Smith and Curvelet transforms
Hölder regularity estimation by Hart Smith and Curvelet transforms Jouni Sampo Lappeenranta University Of Technology Department of Mathematics and Physics Finland 18th September 2007 This research is done
More informationPLANAR KINETIC EQUATIONS OF MOTION (Section 17.2)
PLANAR KINETIC EQUATIONS OF MOTION (Section 17.2) We will limit our study of planar kinetics to rigid bodies that are symmetric with respect to a fixed reference plane. As discussed in Chapter 16, when
More informationGG612 Lecture 3. Outline
GG61 Lecture 3 Strain and Stress Should complete infinitesimal strain by adding rota>on. Outline Matrix Opera+ons Strain 1 General concepts Homogeneous strain 3 Matrix representa>ons 4 Squares of line
More informationMoving layers and interfaces : a thermodynamical approach of Wear Contact
Moving layers and interfaces : a thermodynamical approach of Wear Contact C. Stolz M. Peigney, M. Dragon-Louiset LMS, CNRS-UMR7649, Ecole polytechnique, Palaiseau LaMSID, CNRS-UMR2832, EdF R&D, Clamart
More informationRatcheting deformation in thin film structures
Ratcheting deformation in thin film structures Z. SUO Princeton University Work with MIN HUANG, Rui Huang, Jim Liang, Jean Prevost Princeton University Q. MA, H. Fujimoto, J. He Intel Corporation Interconnect
More informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
More information0 = p. 2 x + 2 w. z +ν w
Solution (Elliptical pipe flow (a Using the Navier Stokes equations in three dimensional cartesian coordinates, given that u =, v = and w = w(x,y only, and assuming no body force, we are left with = p
More informationC:\Users\whit\Desktop\Active\304_2012_ver_2\_Notes\4_Torsion\1_torsion.docx 6
C:\Users\whit\Desktop\Active\304_2012_ver_2\_Notes\4_Torsion\1_torsion.doc 6 p. 1 of Torsion of circular bar The cross-sections rotate without deformation. The deformation that does occur results from
More information( ) = x( u, v) i + y( u, v) j + z( u, v) k
Math 8 ection 16.6 urface Integrals The relationship between surface integrals and surface area is much the same as the relationship between line integrals and arc length. uppose f is a function of three
More informationContact Mechanics and Elements of Tribology
Contact Mechanics and Elements of Tribology Lectures 2-3. Mechanical Contact Vladislav A. Yastrebov MINES ParisTech, PSL Research University, Centre des Matériaux, CNRS UMR 7633, Evry, France @ Centre
More informationProblem Set 5 Math 213, Fall 2016
Problem Set 5 Math 213, Fall 216 Directions: Name: Show all your work. You are welcome and encouraged to use Mathematica, or similar software, to check your answers and aid in your understanding of the
More informationAn analytic form for the pair distribution function and rheology of a dilute suspension of rough spheres in plane strain flow
Under consideration for publication in J. Fluid Mech. 1 An analytic form for the pair distribution function and rheology of a dilute suspension of rough spheres in plane strain flow By H E L E N J. W I
More informationMAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design
MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design Cornell University, Fall 2009 Nicholas Zabaras Materials Process Design and Control Laboratory Sibley School of Mechanical and Aerospace
More informationMicrostructural studies for rheology. Chapter 7: Microstructural studies for rheology. Separation of length scales. Micro & macro views
Chapter 7: Microstructural studies for rheology Microstructural studies for rheology To calculate the flow of complex fluids, need governing equations, in particular, the constitutive equation relating
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More informationUniversity of Alberta. Math 214 Sample Exam Math 214 Solutions
University of Alberta Math 14 Sample Exam Math 14 Solutions 1. Test the following series for convergence or divergence (a) ( n +n+1 3n +n+1 )n, (b) 3 n (n +1) (c) SOL: n!, arccos( n n +1 ), (a) ( n +n+1
More informationMathematical Relativity, Spring 2017/18 Instituto Superior Técnico
Mathematical Relativity, Spring 2017/18 Instituto Superior Técnico 1. Starting from R αβµν Z ν = 2 [α β] Z µ, deduce the components of the Riemann curvature tensor in terms of the Christoffel symbols.
More information