Lab Project 4a MATLAB Model of Oscillatory Flow

Size: px
Start display at page:

Download "Lab Project 4a MATLAB Model of Oscillatory Flow"

Transcription

1 Lab Project 4a MATLAB Model of Oscillatory Flow Goals Prepare analytical and computational solutions for transient flow between the two upright arms of a U-shaped tube Compare the theoretical solutions to experimental data Background: Oscillating systems are commonly described by second order differential equations of the form Ax" + Bx' + Cx = q(t) where x is a state variable such as position, volume, or electrical charge, and q(t) is a forcing function. For a mass-and-spring system, A is the mass, B is a damping coefficient, C is the spring constant, and x is the position of the mass. The equation is typically written as a balance of forces, so every term should have units of force. Liquid oscillations between the two upright tubes in our flow system can be described by this equation IF we assume that flow is laminar, making the pressure drop a linear function of flow rate. If x is a liquid volume, then x' is the flow rate and x" is the acceleration (rate of change of flow). To generate our differential equation, we write a dynamic force balance on the water in the tube. The driving force due to the different heights in the two water columns is Cx, Bx' is the drag force due to viscosity in the tube, and Ax" is the net force that accelerates the fluid. The driving force q(t) is the part of the driving pressure that is not due to the difference in column heights, and is a constant (based on the equilibrium height) for our experiments with transient oscillations. In fluid mechanics it is often convenient to make every term have units of pressure, and the forces can be converted to pressure by dividing x and q(t) by the cross-sectional area of the tube. When flow is turbulent, the pressure drop through a smooth pipe is approximately a function of velocity to the 1.8 power, and the equation becomes non-linear: Ax" + D(x') Cx = q(t). This derivation is not particularly difficult, IF you are familiar with the diagram that relates the Reynolds number to the friction factor for turbulent flow through a tube. Because you might not be familiar with it, it will be derived in lab for any interested students so you may include it with your homework. On the other hand, it is extremely difficult to solve the equation analytically. Instead, we will calculate the time course of the oscillations using a finite difference algorithm in MATLAB. Data collection and MATLAB computation in the turbulent-flow case are just as easy as they are when flow is laminar. print date: 11/11/11 1

2 Previous cohorts of bioengineers have done an experiment to measure the transient behavior of water sloshing back and forth in a U-shaped tube. The basic procedure from 2010 is this: A long tube was mounted as a vertical U, as shown in figure 1. It was filled with enough water so the level, at equilibrium, is about halfway between the curve and the top of the tube. Food coloring and a ruler/meter stick were used to enable video data collection. Using a syringe, water was lifted into one arm of the U, and the water was released and allowed to come to equilibrium. The event was captured on video, and the water height h(t) was extracted from the video. The experiment was repeated for multiple sizes of tubing. Examples of the video and data sets will be provided for you. You do not need to analyze the video yourself; you may use the tabulated data. Procedure Part 1 a) Plot the water heights for one experimental trial from each of the tubing sizes. Do not include the plots here. Instead, categorize each response as over-, under-, or critically damped. b) Estimate the maximum Reynolds number for each of the trials that you plotted. c) Categorize each trial as mostly laminar or mostly turbulent flow. Part 2 Assume that water is contained in a U-shaped tube as shown in Fig. 1. When both ends are open to atmosphere, the water level is the same at some equilibrium height on both sides of the U. A syringe was attached to the left end of the U and was used to draw water up to some level h 1 (0). The water was then released and allowed to come to equilibrium. a) Assuming that flow is laminar, derive a differential equation that describes the transient solution, h(t), for the height of water in one side of the tube after it is released from the initial height h 1 (0). You should get a linear differential equation. b) Find an analytical solution to the linear differential equation in terms of system dimensions and liquid density. c) Starting with the MATLAB code provided for you, create a computational solution to the linear differential equation. Include the relevant lines of code in your MATLAB appendix. d) For each of the trials that you categorized as laminar in part 1, plot the experimental result, analytical solution and computational solution on the same axes. print date: 11/11/11 2

3 h 1 (0) h 1 (t) h 2 (t) Equilibrium height h eq R BEND h 1 (0) is chosen such that the water level never drops below this point. Tube of diameter D Figure 1 U-tube arrangement for flow oscillation experiments. The tube is filled up to the dashed line, before it is pulled up to h 1 (0) on the left side, and released. Part 3 a) Assuming that flow is turbulent, derive a nonlinear differential equation that describes the transient solution, h(t), for the height of water in one side of the tube after it is released from the initial height h 1 (0). Note: One of the steps should be to derive a formula that relates the Fanning friction factor to the velocity of the water in the tube. The steps are shown in the BIOEN 325 course notes, but I expect you will not get exactly the same answer as shown in the notes because you should draw your own tangent line on the Moody diagram. However, your answer might be the same after you round to one decimal place. b) Start with the MATLAB code provided for you. For each of the trials that you categorized as turbulent in part 1, plot the experimental result and computational solution on the same axes. If you did not find any to be turbulent, just compare the one with the highest Reynolds number. Part 4 a) Using the analytical solution from part 2, find the initial height necessary to make the Reynolds number reach 10,000 in the 1 U-tube. Note that this is a linear differential equation, so the starting height is proportional to the maximum Re. print date: 11/11/11 3

4 b) Starting with the MATLAB code provided for you, create and plot the following curves. Assume the starting height h 1 (0) is the one you found in the previous step. o A computational solution to the linear differential equation (assumes laminar flow). o A computational solution to the nonlinear differential equation; Reporting The main challenge in this assignment is to create and solve initial value problems (differential equation plus initial conditions). The experimental results are what they are, and you do not need to match the experimental and simulated data exactly. You should, however, see the same trends over the range of tubing sizes. Discussion: Discuss why some flow arrangements reach equilibrium without oscillating, while others respond with a decaying sinusoid. Some responses look like a combination of monotonically decaying exponential followed by oscillations. Use the equations you derived for laminar and turbulent flow to explain these cases where flow oscillates but does not match a simple decaying sinusoid. Sources of Error: Discuss approximations or simplifications made in the flow-resistance calculations and the initial value problems. State, with justification, how much effect you expect each simplification to have (large or small -- no quantitative error analysis is necessary). Lengthy calculations and new MATLAB code may be stowed in an appendix. You do not need to include MATLAB code if you did not modify it. print date: 11/11/11 4

5 Supplement Creating a numerical solution for a second-order initial value problem Consider the following initial value problem, which is composed of a linear second-order differential equation and initial conditions for position and velocity. Ax" + Bx' + Cx = q(t) x(0) = <initial height>, x'(0)=0. The ODE solvers in MATLAB operate on first-order differential equations, so we first want to reduce our second-order ODE to a system of first-order differential equations. First, we introduce the vector y: x y = x Next, we rewrite the second-order equation above as a system of first-order equations. d x 0 1 x 0 = + dt x C/ A B/ A x q( t)/ A We then create a MATLAB function that calculates the increments in x and x' (i.e., in y) for single time steps. function dy = secondode(t,y) % function to be integrated A=?; B=?; C=?; % Enter values for A, B, C dy = zeros(2,1); dy(1) = y(2); dy(2) = -C/A*y(1)- B/A*y(2) + <Enter q(t) here>/a; Let us suppose that our initial value is 20 cm, 1 our initial velocity is zero, and we want to find the solution from 0 to 9 seconds. These values would be set as follows: y0 = [1; 0]; tspan = [0 9]; Now we can solve the system using ode45(). [t, y] = ode45(@secondode, tspan, y0); The variable t will be a column vector of time values; the output y will be an array in which the first column is x(t) and the second column is x'(t). 1 The units that you use are not going to change your answer when the differential equation is linear, but they can make a big difference if the equation is non-linear. print date: 11/11/11 1

6 If the flow is turbulent, then pressure drop is no longer a linear function of flow rate. In this case the differential equation becomes non-linear in the x' term. Ax" + B(x') β + Cx = q(t) Typically β is some value between 1 and 2, and can be approximated from the Moody friction factor chart and Darcy's pressure drop equation. The MATLAB solution is the same as before except that the y(2) term in the last line includes the exponent β. function dy = secondodenonlin(t,y) % function to be integrated A=?; B=?; C=?; g=? % Enter values for A, B, C dy = zeros(2,1); dy(1) = y(2); dy(2) = -C/A*y(1)- B/A*(y(2))^beta + q/a; This code will return an error when y(2) is negative. Two more changes are needed to complete the y(2) term the first change gets rid of the MATLAB error but introduces a model error (ie the result is wrong). Can you find the new problem and fix it? print date: 11/11/11 2

Three-Tank Experiment

Three-Tank Experiment Three-Tank Experiment Overview The three-tank experiment focuses on application of the mechanical balance equation to a transient flow. Three tanks are interconnected by Schedule 40 pipes of nominal diameter

More information

Fluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 17 Laminar and Turbulent flows

Fluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 17 Laminar and Turbulent flows Fluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 17 Laminar and Turbulent flows Welcome back to the video course on fluid mechanics. In

More information

(Refer Slide Time 1:25)

(Refer Slide Time 1:25) Mechanical Measurements and Metrology Prof. S. P. Venkateshan Department of Mechanical Engineering Indian Institute of Technology, Madras Module - 2 Lecture - 24 Transient Response of Pressure Transducers

More information

Physics 3 Summer 1990 Lab 7 - Hydrodynamics

Physics 3 Summer 1990 Lab 7 - Hydrodynamics Physics 3 Summer 1990 Lab 7 - Hydrodynamics Theory Consider an ideal liquid, one which is incompressible and which has no internal friction, flowing through pipe of varying cross section as shown in figure

More information

Driven Harmonic Oscillator

Driven Harmonic Oscillator Driven Harmonic Oscillator Physics 6B Lab Experiment 1 APPARATUS Computer and interface Mechanical vibrator and spring holder Stands, etc. to hold vibrator Motion sensor C-209 spring Weight holder and

More information

Chapter 10 Flow in Conduits

Chapter 10 Flow in Conduits Chapter 10 Flow in Conduits 10.1 Classifying Flow Laminar Flow and Turbulent Flow Laminar flow Unpredictable Turbulent flow Near entrance: undeveloped developing flow In developing flow, the wall shear

More information

Exponential decay. The deviations in amplitude over 30 periods rise to more than ±20%. Fig 1 a rod and ball pendulum

Exponential decay. The deviations in amplitude over 30 periods rise to more than ±20%. Fig 1 a rod and ball pendulum Exponential decay A counter example There is a common belief that the damping of the motion of a pendulum in air is exponential, or nearly so, in all situations. To explore the limits of that approximation

More information

PHYSICS 211 LAB #8: Periodic Motion

PHYSICS 211 LAB #8: Periodic Motion PHYSICS 211 LAB #8: Periodic Motion A Lab Consisting of 6 Activities Name: Section: TA: Date: Lab Partners: Circle the name of the person to whose report your group printouts will be attached. Individual

More information

Mechanical Energy and Simple Harmonic Oscillator

Mechanical Energy and Simple Harmonic Oscillator Mechanical Energy and Simple Harmonic Oscillator Simple Harmonic Motion Hooke s Law Define system, choose coordinate system. Draw free-body diagram. Hooke s Law! F spring =!kx ˆi! kx = d x m dt Checkpoint

More information

S.E. (Mech.) (First Sem.) EXAMINATION, (Common to Mech/Sandwich) FLUID MECHANICS (2008 PATTERN) Time : Three Hours Maximum Marks : 100

S.E. (Mech.) (First Sem.) EXAMINATION, (Common to Mech/Sandwich) FLUID MECHANICS (2008 PATTERN) Time : Three Hours Maximum Marks : 100 Total No. of Questions 12] [Total No. of Printed Pages 8 Seat No. [4262]-113 S.E. (Mech.) (First Sem.) EXAMINATION, 2012 (Common to Mech/Sandwich) FLUID MECHANICS (2008 PATTERN) Time : Three Hours Maximum

More information

Applications of Second-Order Differential Equations

Applications of Second-Order Differential Equations Applications of Second-Order Differential Equations ymy/013 Building Intuition Even though there are an infinite number of differential equations, they all share common characteristics that allow intuition

More information

Lecture 14: Ordinary Differential Equation I. First Order

Lecture 14: Ordinary Differential Equation I. First Order Lecture 14: Ordinary Differential Equation I. First Order 1. Key points Maple commands dsolve 2. Introduction We consider a function of one variable. An ordinary differential equations (ODE) specifies

More information

Section 11.1 What is a Differential Equation?

Section 11.1 What is a Differential Equation? 1 Section 11.1 What is a Differential Equation? Example 1 Suppose a ball is dropped from the top of a building of height 50 meters. Let h(t) denote the height of the ball after t seconds, then it is known

More information

Applied Fluid Mechanics

Applied Fluid Mechanics Applied Fluid Mechanics 1. The Nature of Fluid and the Study of Fluid Mechanics 2. Viscosity of Fluid 3. Pressure Measurement 4. Forces Due to Static Fluid 5. Buoyancy and Stability 6. Flow of Fluid and

More information

Figure XP3.1 (a) Mass in equilibrium, (b) Freebody diagram, (c) Kinematic constraint relation Example Problem 3.1 Figure XP3.1 illustrates a mass m

Figure XP3.1 (a) Mass in equilibrium, (b) Freebody diagram, (c) Kinematic constraint relation Example Problem 3.1 Figure XP3.1 illustrates a mass m LECTURE 7. MORE VIBRATIONS ` Figure XP3.1 (a) Mass in equilibrium, (b) Freebody diagram, (c) Kinematic constraint relation Example Problem 3.1 Figure XP3.1 illustrates a mass m that is in equilibrium and

More information

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) 28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.

More information

Basic Fluid Mechanics

Basic Fluid Mechanics Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible

More information

Applied Fluid Mechanics

Applied Fluid Mechanics Applied Fluid Mechanics 1. The Nature of Fluid and the Study of Fluid Mechanics 2. Viscosity of Fluid 3. Pressure Measurement 4. Forces Due to Static Fluid 5. Buoyancy and Stability 6. Flow of Fluid and

More information

Piping Systems and Flow Analysis (Chapter 3)

Piping Systems and Flow Analysis (Chapter 3) Piping Systems and Flow Analysis (Chapter 3) 2 Learning Outcomes (Chapter 3) Losses in Piping Systems Major losses Minor losses Pipe Networks Pipes in series Pipes in parallel Manifolds and Distribution

More information

Experiment 5. Simple Harmonic Motion

Experiment 5. Simple Harmonic Motion Reading and Problems: Chapters 7,8 Problems 7., 8. Experiment 5 Simple Harmonic Motion Goals. To understand the properties of an oscillating system governed by Hooke s Law.. To study the effects of friction

More information

Topic 5 Notes Jeremy Orloff. 5 Homogeneous, linear, constant coefficient differential equations

Topic 5 Notes Jeremy Orloff. 5 Homogeneous, linear, constant coefficient differential equations Topic 5 Notes Jeremy Orloff 5 Homogeneous, linear, constant coefficient differential equations 5.1 Goals 1. Be able to solve homogeneous constant coefficient linear differential equations using the method

More information

Harmonic Oscillator. Outline. Oscillatory Motion or Simple Harmonic Motion. Oscillatory Motion or Simple Harmonic Motion

Harmonic Oscillator. Outline. Oscillatory Motion or Simple Harmonic Motion. Oscillatory Motion or Simple Harmonic Motion Harmonic Oscillator Mass-Spring Oscillator Resonance The Pendulum Physics 109, Class Period 13 Experiment Number 11 in the Physics 121 Lab Manual (page 65) Outline Simple harmonic motion The vertical mass-spring

More information

Study fluid dynamics. Understanding Bernoulli s Equation.

Study fluid dynamics. Understanding Bernoulli s Equation. Chapter Objectives Study fluid dynamics. Understanding Bernoulli s Equation. Chapter Outline 1. Fluid Flow. Bernoulli s Equation 3. Viscosity and Turbulence 1. Fluid Flow An ideal fluid is a fluid that

More information

Math 3313: Differential Equations Second-order ordinary differential equations

Math 3313: Differential Equations Second-order ordinary differential equations Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties

More information

Water Circuit Lab. The pressure drop along a straight pipe segment can be calculated using the following set of equations:

Water Circuit Lab. The pressure drop along a straight pipe segment can be calculated using the following set of equations: Water Circuit Lab When a fluid flows in a conduit, there is friction between the flowing fluid and the pipe walls. The result of this friction is a net loss of energy in the flowing fluid. The fluid pressure

More information

COURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour. Basic Equations in fluid Dynamics

COURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour. Basic Equations in fluid Dynamics COURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour Basic Equations in fluid Dynamics Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET 1 Description of Fluid

More information

150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces

150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with

More information

Physics 2310 Lab #3 Driven Harmonic Oscillator

Physics 2310 Lab #3 Driven Harmonic Oscillator Physics 2310 Lab #3 Driven Harmonic Oscillator M. Pierce (adapted from a lab by the UCLA Physics & Astronomy Department) Objective: The objective of this experiment is to characterize the behavior of a

More information

Stream line, turbulent flow and Viscosity of liquids - Poiseuille s Method

Stream line, turbulent flow and Viscosity of liquids - Poiseuille s Method Stream line, turbulent flow and Viscosity of liquids - Poiseuille s Method Dr D. Arun Kumar Assistant Professor Department of Physical Sciences Bannari Amman Institute of Technology Sathyamangalam General

More information

ENG3103 Engineering Problem Solving Computations Semester 2, 2013

ENG3103 Engineering Problem Solving Computations Semester 2, 2013 Assessment: Assignment 2 Due: 16 September 2013 Marks: 100 Value: 10 % Question 1 (70 marks) Introduction You are designing a pipe network system that transfers water from the upper pipe to the lower pipe.

More information

Physics 326 Lab 6 10/18/04 DAMPED SIMPLE HARMONIC MOTION

Physics 326 Lab 6 10/18/04 DAMPED SIMPLE HARMONIC MOTION DAMPED SIMPLE HARMONIC MOTION PURPOSE To understand the relationships between force, acceleration, velocity, position, and period of a mass undergoing simple harmonic motion and to determine the effect

More information

Frictional Losses in Straight Pipe

Frictional Losses in Straight Pipe 2/2/206 CM325 Fundamentals of Chemical Engineering Laboratory Prelab Preparation for Frictional Losses in Straight Pipe Professor Faith Morrison Department of Chemical Engineering Michigan Technological

More information

FACULTY OF CHEMICAL & ENERGY ENGINEERING FLUID MECHANICS LABORATORY TITLE OF EXPERIMENT: MINOR LOSSES IN PIPE (E4)

FACULTY OF CHEMICAL & ENERGY ENGINEERING FLUID MECHANICS LABORATORY TITLE OF EXPERIMENT: MINOR LOSSES IN PIPE (E4) FACULTY OF CHEMICAL & ENERGY ENGINEERING FLUID MECHANICS LABORATORY TITLE OF EXPERIMENT: MINOR LOSSES IN PIPE (E4) 1 1.0 Objectives The objective of this experiment is to calculate loss coefficient (K

More information

Lecture 22. Mechanical Energy Balance

Lecture 22. Mechanical Energy Balance Lecture 22 Mechanical Energy Balance Contents Exercise 1 Exercise 2 Exercise 3 Key Words: Fluid flow, Macroscopic Balance, Frictional Losses, Turbulent Flow Exercise 1 It is proposed to install a fan to

More information

Internal Forced Convection. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.

Internal Forced Convection. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Internal Forced Convection Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Introduction Pipe circular cross section. Duct noncircular cross section. Tubes small-diameter

More information

Fluid Flow Analysis Penn State Chemical Engineering

Fluid Flow Analysis Penn State Chemical Engineering Fluid Flow Analysis Penn State Chemical Engineering Revised Spring 2015 Table of Contents LEARNING OBJECTIVES... 1 EXPERIMENTAL OBJECTIVES AND OVERVIEW... 1 PRE-LAB STUDY... 2 EXPERIMENTS IN THE LAB...

More information

Introduction to First Order Equations Sections

Introduction to First Order Equations Sections A B I L E N E C H R I S T I A N U N I V E R S I T Y Department of Mathematics Introduction to First Order Equations Sections 2.1-2.3 Dr. John Ehrke Department of Mathematics Fall 2012 Course Goals The

More information

Section 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation

Section 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1 General Results

More information

s and FE X. A. Flow measurement B. properties C. statics D. impulse, and momentum equations E. Pipe and other internal flow 7% of FE Morning Session I

s and FE X. A. Flow measurement B. properties C. statics D. impulse, and momentum equations E. Pipe and other internal flow 7% of FE Morning Session I Fundamentals of Engineering (FE) Exam General Section Steven Burian Civil & Environmental Engineering October 26, 2010 s and FE X. A. Flow measurement B. properties C. statics D. impulse, and momentum

More information

FE Exam Fluids Review October 23, Important Concepts

FE Exam Fluids Review October 23, Important Concepts FE Exam Fluids Review October 3, 013 mportant Concepts Density, specific volume, specific weight, specific gravity (Water 1000 kg/m^3, Air 1. kg/m^3) Meaning & Symbols? Stress, Pressure, Viscosity; Meaning

More information

Principles of Convection

Principles of Convection Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid

More information

Chapter 8: Flow in Pipes

Chapter 8: Flow in Pipes Objectives 1. Have a deeper understanding of laminar and turbulent flow in pipes and the analysis of fully developed flow 2. Calculate the major and minor losses associated with pipe flow in piping networks

More information

1-Reynold s Experiment

1-Reynold s Experiment Lect.No.8 2 nd Semester Flow Dynamics in Closed Conduit (Pipe Flow) 1 of 21 The flow in closed conduit ( flow in pipe ) is differ from this occur in open channel where the flow in pipe is at a pressure

More information

2.3 Damping, phases and all that

2.3 Damping, phases and all that 2.3. DAMPING, PHASES AND ALL THAT 107 2.3 Damping, phases and all that If we imagine taking our idealized mass on a spring and dunking it in water or, more dramatically, in molasses), then there will be

More information

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of

More information

ME 230 Kinematics and Dynamics

ME 230 Kinematics and Dynamics ME 230 Kinematics and Dynamics Wei-Chih Wang Department of Mechanical Engineering University of Washington Lecture 6: Particle Kinetics Kinetics of a particle (Chapter 13) - 13.4-13.6 Chapter 13: Objectives

More information

Hydraulics and hydrology

Hydraulics and hydrology Hydraulics and hydrology - project exercises - Class 4 and 5 Pipe flow Discharge (Q) (called also as the volume flow rate) is the volume of fluid that passes through an area per unit time. The discharge

More information

Engineers Edge, LLC PDH & Professional Training

Engineers Edge, LLC PDH & Professional Training 510 N. Crosslane Rd. Monroe, Georgia 30656 (770) 266-6915 fax (678) 643-1758 Engineers Edge, LLC PDH & Professional Training Copyright, All Rights Reserved Engineers Edge, LLC Pipe Flow-Friction Factor

More information

MOMENTS OF INERTIA FOR AREAS, RADIUS OF GYRATION OF AN AREA, & MOMENTS OF INTERTIA BY INTEGRATION

MOMENTS OF INERTIA FOR AREAS, RADIUS OF GYRATION OF AN AREA, & MOMENTS OF INTERTIA BY INTEGRATION MOMENTS OF INERTIA FOR AREAS, RADIUS OF GYRATION OF AN AREA, & MOMENTS OF INTERTIA BY INTEGRATION Today s Objectives: Students will be able to: a) Define the moments of inertia (MoI) for an area. b) Determine

More information

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered

More information

ELEC9712 High Voltage Systems. 1.2 Heat transfer from electrical equipment

ELEC9712 High Voltage Systems. 1.2 Heat transfer from electrical equipment ELEC9712 High Voltage Systems 1.2 Heat transfer from electrical equipment The basic equation governing heat transfer in an item of electrical equipment is the following incremental balance equation, with

More information

Applied Fluid Mechanics

Applied Fluid Mechanics Applied Fluid Mechanics 1. The Nature of Fluid and the Study of Fluid Mechanics 2. Viscosity of Fluid 3. Pressure Measurement 4. Forces Due to Static Fluid 5. Buoyancy and Stability 6. Flow of Fluid and

More information

Chapter 14. Fluid Mechanics

Chapter 14. Fluid Mechanics Chapter 14 Fluid Mechanics States of Matter Solid Has a definite volume and shape Liquid Has a definite volume but not a definite shape Gas unconfined Has neither a definite volume nor shape All of these

More information

VORTEX LEVITATION. Toshiharu Kagawa 1 and Xin Li 2

VORTEX LEVITATION. Toshiharu Kagawa 1 and Xin Li 2 VORTEX LEVITATION Toshiharu Kagawa 1 and Xin Li ABSTRACT In this paper, a new pneumatic levitation method, called vortex levitation, is introduced. Vortex levitation can achieve non-contact handling by

More information

TALLINN UNIVERSITY OF TECHNOLOGY, DIVISION OF PHYSICS 13. STOKES METHOD

TALLINN UNIVERSITY OF TECHNOLOGY, DIVISION OF PHYSICS 13. STOKES METHOD 13. STOKES METHOD 1. Objective To determine the coefficient of viscosity of a known fluid using Stokes method.. Equipment needed A glass vessel with glycerine, micrometer calliper, stopwatch, ruler. 3.

More information

Inaugural University of Michigan Science Olympiad Invitational Tournament. Hovercraft

Inaugural University of Michigan Science Olympiad Invitational Tournament. Hovercraft Inaugural University of Michigan Science Olympiad Invitational Tournament Test length: 50 Minutes Hovercraft Team number: Team name: Student names: Instructions: Do not open this test until told to do

More information

Fluid Mechanics Testbank By David Admiraal

Fluid Mechanics Testbank By David Admiraal Fluid Mechanics Testbank By David Admiraal This testbank was created for an introductory fluid mechanics class. The primary intentions of the testbank are to help students improve their performance on

More information

PROBLEMS. PDF processed with CutePDF evaluation edition

PROBLEMS. PDF processed with CutePDF evaluation edition 114 Chapter 3 Built-In MATLAB Functions PROBLEMS Elementary Math Functions 3.1 Find the cube root of -5, both by using the nthroot function and by raising -5 to the 1/3 power. Explain the difference in

More information

Chemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017

Chemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017 Chemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017 Objective: Text: To introduce the basic concepts of fluid mechanics and heat transfer necessary for solution of engineering

More information

Exam tomorrow on Chapter 15, 16, and 17 (Oscilla;ons and Waves 1 &2)

Exam tomorrow on Chapter 15, 16, and 17 (Oscilla;ons and Waves 1 &2) Exam tomorrow on Chapter 15, 16, and 17 (Oscilla;ons and Waves 1 &2) What to study: Quiz 6 Homework problems for Chapters 15 & 16 Material indicated in the following review slides Other Specific things:

More information

Oscillatory Motion SHM

Oscillatory Motion SHM Chapter 15 Oscillatory Motion SHM Dr. Armen Kocharian Periodic Motion Periodic motion is motion of an object that regularly repeats The object returns to a given position after a fixed time interval A

More information

Harmonic oscillator. U(x) = 1 2 bx2

Harmonic oscillator. U(x) = 1 2 bx2 Harmonic oscillator The harmonic oscillator is a familiar problem from classical mechanics. The situation is described by a force which depends linearly on distance as happens with the restoring force

More information

OE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004

OE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004 OE465 Vaclav Matousek October 13, 004 1 Dredge Vermelding Pumps onderdeel and Slurry organisatie Transport OE465 Vaclav Matousek October 13, 004 Dredge Vermelding Pumps onderdeel and Slurry organisatie

More information

I. Pre-Lab Introduction

I. Pre-Lab Introduction I. Pre-Lab Introduction Please complete the following pages before the lab by filling in the requested items. A. Atomic notation: Atoms are composed of a nucleus containing neutrons and protons surrounded

More information

What we know about Fluid Mechanics. What we know about Fluid Mechanics

What we know about Fluid Mechanics. What we know about Fluid Mechanics What we know about Fluid Mechanics 1. Survey says. 3. Image from: www.axs.com 4. 5. 6. 1 What we know about Fluid Mechanics 1. MEB (single input, single output, steady, incompressible, no rxn, no phase

More information

Lab 11 - Free, Damped, and Forced Oscillations

Lab 11 - Free, Damped, and Forced Oscillations Lab 11 Free, Damped, and Forced Oscillations L11-1 Name Date Partners Lab 11 - Free, Damped, and Forced Oscillations OBJECTIVES To understand the free oscillations of a mass and spring. To understand how

More information

Chapter 8: Flow in Pipes

Chapter 8: Flow in Pipes 8-1 Introduction 8-2 Laminar and Turbulent Flows 8-3 The Entrance Region 8-4 Laminar Flow in Pipes 8-5 Turbulent Flow in Pipes 8-6 Fully Developed Pipe Flow 8-7 Minor Losses 8-8 Piping Networks and Pump

More information

Aircraft Dynamics First order and Second order system

Aircraft Dynamics First order and Second order system Aircraft Dynamics First order and Second order system Prepared by A.Kaviyarasu Assistant Professor Department of Aerospace Engineering Madras Institute Of Technology Chromepet, Chennai Aircraft dynamic

More information

Chapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx

Chapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx Chapter 1 Lecture Notes Chapter 1 Oscillatory Motion Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx When the mass is released, the spring will pull

More information

Chapter 1: Basic Concepts

Chapter 1: Basic Concepts What is a fluid? A fluid is a substance in the gaseous or liquid form Distinction between solid and fluid? Solid: can resist an applied shear by deforming. Stress is proportional to strain Fluid: deforms

More information

LAB 11: FREE, DAMPED, AND FORCED OSCILLATIONS

LAB 11: FREE, DAMPED, AND FORCED OSCILLATIONS Lab 11 ree, amped, and orced Oscillations 135 Name ate Partners OBJECTIVES LAB 11: REE, AMPE, AN ORCE OSCILLATIONS To understand the free oscillations of a mass and spring. To understand how energy is

More information

Thursday, August 4, 2011

Thursday, August 4, 2011 Chapter 16 Thursday, August 4, 2011 16.1 Springs in Motion: Hooke s Law and the Second-Order ODE We have seen alrealdy that differential equations are powerful tools for understanding mechanics and electro-magnetism.

More information

1.1 GRAPHS AND LINEAR FUNCTIONS

1.1 GRAPHS AND LINEAR FUNCTIONS MATHEMATICS EXTENSION 4 UNIT MATHEMATICS TOPIC 1: GRAPHS 1.1 GRAPHS AND LINEAR FUNCTIONS FUNCTIONS The concept of a function is already familiar to you. Since this concept is fundamental to mathematics,

More information

Chapter 13 Lecture. Essential University Physics Richard Wolfson 2 nd Edition. Oscillatory Motion Pearson Education, Inc.

Chapter 13 Lecture. Essential University Physics Richard Wolfson 2 nd Edition. Oscillatory Motion Pearson Education, Inc. Chapter 13 Lecture Essential University Physics Richard Wolfson nd Edition Oscillatory Motion Slide 13-1 In this lecture you ll learn To describe the conditions under which oscillatory motion occurs To

More information

Chapter 10. Solids and Fluids

Chapter 10. Solids and Fluids Chapter 10 Solids and Fluids Surface Tension Net force on molecule A is zero Pulled equally in all directions Net force on B is not zero No molecules above to act on it Pulled toward the center of the

More information

(Refer Slide Time: 1: 19)

(Refer Slide Time: 1: 19) Mechanical Measurements and Metrology Prof. S. P. Venkateshan Department of Mechanical Engineering Indian Institute of Technology, Madras Module - 4 Lecture - 46 Force Measurement So this will be lecture

More information

2. For a S.H.O. determine, (a) the total energy (E), the kinetic and potential energies. of half amplitude:

2. For a S.H.O. determine, (a) the total energy (E), the kinetic and potential energies. of half amplitude: The amplitude of vibration and hence, the energy transferred into the vibrating system is found to depend on the difference between f and, its maximum when the frequency of the external force is equal

More information

a. Follow the Start-Up Procedure in the laboratory manual. Note the safety rules.

a. Follow the Start-Up Procedure in the laboratory manual. Note the safety rules. Lab #1 - Free Vibration Name: Date: Section / Group: Procedure Steps (from lab manual): a. Follow the Start-Up Procedure in the laboratory manual. Note the safety rules. b. Locate the various springs and

More information

The... of a particle is defined as its change in position in some time interval.

The... of a particle is defined as its change in position in some time interval. Distance is the. of a path followed by a particle. Distance is a quantity. The... of a particle is defined as its change in position in some time interval. Displacement is a.. quantity. The... of a particle

More information

Ch 6 Using Newton s Laws. Applications to mass, weight, friction, air resistance, and periodic motion

Ch 6 Using Newton s Laws. Applications to mass, weight, friction, air resistance, and periodic motion Ch 6 Using Newton s Laws Applications to mass, weight, friction, air resistance, and periodic motion Newton s 2 nd Law Applied Galileo hypothesized that all objects gain speed at the same rate (have the

More information

BIOEN LECTURE 18: VISCOELASTIC MODELS

BIOEN LECTURE 18: VISCOELASTIC MODELS BIOEN 326 2013 LECTURE 18: VISCOELASTIC MODELS Definition of Viscoelasticity. Until now, we have looked at time-independent behaviors. This assumed that materials were purely elastic in the conditions

More information

MEMORANDUM. To: Dr. Andrew Kean, American Petroleum Institute (API) From: Josh Baida Abel Estrada

MEMORANDUM. To: Dr. Andrew Kean, American Petroleum Institute (API) From: Josh Baida Abel Estrada MEMORANDUM To: Dr. Andrew Kean, American Petroleum Institute (API) akean@calpoly.edu From: Josh Baida Abel Estrada jbaida@calpoly.edu Cesia Cazares cscazare@calpoly.edu Date: October 27, 2015 RE: Viscosity

More information

Fluid Mechanics. du dy

Fluid Mechanics. du dy FLUID MECHANICS Technical English - I 1 th week Fluid Mechanics FLUID STATICS FLUID DYNAMICS Fluid Statics or Hydrostatics is the study of fluids at rest. The main equation required for this is Newton's

More information

UNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics

UNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics UNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics Physics 117.3 MIDTERM TEST February 11, 009 Time: 90 minutes NAME: (Last) Please Print (Given) STUDENT NO.: LECTURE SECTION (please

More information

Lab 1: Dynamic Simulation Using Simulink and Matlab

Lab 1: Dynamic Simulation Using Simulink and Matlab Lab 1: Dynamic Simulation Using Simulink and Matlab Objectives In this lab you will learn how to use a program called Simulink to simulate dynamic systems. Simulink runs under Matlab and uses block diagrams

More information

External Flow and Boundary Layer Concepts

External Flow and Boundary Layer Concepts 1 2 Lecture (8) on Fayoum University External Flow and Boundary Layer Concepts By Dr. Emad M. Saad Mechanical Engineering Dept. Faculty of Engineering Fayoum University Faculty of Engineering Mechanical

More information

Solutions to Math 53 Math 53 Practice Final

Solutions to Math 53 Math 53 Practice Final Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points

More information

In the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as

In 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 information

Chapter 5 MATHEMATICAL MODELING OF THE EVACATED SOLAR COLLECTOR. 5.1 Thermal Model of Solar Collector System

Chapter 5 MATHEMATICAL MODELING OF THE EVACATED SOLAR COLLECTOR. 5.1 Thermal Model of Solar Collector System Chapter 5 MATHEMATICAL MODELING OF THE EVACATED SOLAR COLLECTOR This chapter deals with analytical method of finding out the collector outlet working fluid temperature. A dynamic model of the solar collector

More information

ME8230 Nonlinear Dynamics

ME8230 Nonlinear Dynamics ME8230 Nonlinear Dynamics Lecture 1, part 1 Introduction, some basic math background, and some random examples Prof. Manoj Srinivasan Mechanical and Aerospace Engineering srinivasan.88@osu.edu Spring mass

More information

C ONTENTS CHAPTER TWO HEAT CONDUCTION EQUATION 61 CHAPTER ONE BASICS OF HEAT TRANSFER 1 CHAPTER THREE STEADY HEAT CONDUCTION 127

C ONTENTS CHAPTER TWO HEAT CONDUCTION EQUATION 61 CHAPTER ONE BASICS OF HEAT TRANSFER 1 CHAPTER THREE STEADY HEAT CONDUCTION 127 C ONTENTS Preface xviii Nomenclature xxvi CHAPTER ONE BASICS OF HEAT TRANSFER 1 1-1 Thermodynamics and Heat Transfer 2 Application Areas of Heat Transfer 3 Historical Background 3 1-2 Engineering Heat

More information

HOW TO GET A GOOD GRADE ON THE MME 2273B FLUID MECHANICS 1 EXAM. Common mistakes made on the final exam and how to avoid them

HOW TO GET A GOOD GRADE ON THE MME 2273B FLUID MECHANICS 1 EXAM. Common mistakes made on the final exam and how to avoid them HOW TO GET A GOOD GRADE ON THE MME 2273B FLUID MECHANICS 1 EXAM Common mistakes made on the final exam and how to avoid them HOW TO GET A GOOD GRADE ON THE MME 2273B EXAM Introduction You now have a lot

More information

Mass on a Spring C2: Simple Harmonic Motion. Simple Harmonic Motion. Announcements Week 12D1

Mass on a Spring C2: Simple Harmonic Motion. Simple Harmonic Motion. Announcements Week 12D1 Simple Harmonic Motion 8.01 Week 1D1 Today s Reading Assignment MIT 8.01 Course Notes Chapter 3 Simple Harmonic Motion Sections 3.1-3.4 1 Announcements Sunday Tutoring in 6-15 from 1-5 pm Problem Set 9

More information

spring magnet Fig. 7.1 One end of the magnet hangs inside a coil of wire. The coil is connected in series with a resistor R.

spring magnet Fig. 7.1 One end of the magnet hangs inside a coil of wire. The coil is connected in series with a resistor R. 1 A magnet is suspended vertically from a fixed point by means of a spring, as shown in Fig. 7.1. spring magnet coil R Fig. 7.1 One end of the magnet hangs inside a coil of wire. The coil is connected

More information

6. Basic basic equations I ( )

6. Basic basic equations I ( ) 6. Basic basic equations I (4.2-4.4) Steady and uniform flows, streamline, streamtube One-, two-, and three-dimensional flow Laminar and turbulent flow Reynolds number System and control volume Continuity

More information

Experiment 7 Simple Harmonic Motion Reading: Bauer&Westfall Chapter 13 as needed (e.g. definitions of oscillation variables)

Experiment 7 Simple Harmonic Motion Reading: Bauer&Westfall Chapter 13 as needed (e.g. definitions of oscillation variables) PHY191 Experiment 7: Simple Harmonic Motion 7/3/9 Page 1 Experiment 7 Simple Harmonic Motion Reading: Bauer&Westfall Chapter 13 as needed (e.g. definitions of oscillation variables) Homework 8: turn in

More information

Steven Burian Civil & Environmental Engineering September 25, 2013

Steven Burian Civil & Environmental Engineering September 25, 2013 Fundamentals of Engineering (FE) Exam Mechanics Steven Burian Civil & Environmental Engineering September 25, 2013 s and FE Morning ( Mechanics) A. Flow measurement 7% of FE Morning B. properties Session

More information

Measuring the Universal Gravitational Constant, G

Measuring the Universal Gravitational Constant, G Measuring the Universal Gravitational Constant, G Introduction: The universal law of gravitation states that everything in the universe is attracted to everything else. It seems reasonable that everything

More information

Vibration Control Prof. Dr. S. P. Harsha Department of Mechanical & Industrial Engineering Indian Institute of Technology, Roorkee

Vibration Control Prof. Dr. S. P. Harsha Department of Mechanical & Industrial Engineering Indian Institute of Technology, Roorkee Vibration Control Prof. Dr. S. P. Harsha Department of Mechanical & Industrial Engineering Indian Institute of Technology, Roorkee Module - 1 Review of Basics of Mechanical Vibrations Lecture - 2 Introduction

More information

Learning Objectives. Lesson 6: Mathematical Models of Fluid Flow Components. ET 438a Automatic Control Systems Technology 8/27/2015

Learning Objectives. Lesson 6: Mathematical Models of Fluid Flow Components. ET 438a Automatic Control Systems Technology 8/27/2015 Lesson 6: Mathematical Models of Fluid Flow Components ET 438a Automatic Control Systems Technology lesson6et438a.pptx 1 Learning Objectives After this presentation you will be able to: Define the characteristics

More information