Lab Project 4a MATLAB Model of Oscillatory Flow
|
|
- Cora Mitchell
- 5 years ago
- Views:
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 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 informationFluid 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)
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 informationPhysics 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 informationDriven 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 informationChapter 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 informationExponential 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 informationPHYSICS 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 informationMechanical 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 informationS.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 informationApplications 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 informationLecture 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 informationSection 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 informationApplied 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 informationFigure 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 information28. 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 informationBasic 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 informationApplied 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 informationPiping 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 informationExperiment 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 informationTopic 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 informationHarmonic 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 informationStudy 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 informationMath 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 informationWater 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 informationCOURSE 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 information150A 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 informationPhysics 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 informationStream 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 informationENG3103 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 informationPhysics 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 informationFrictional 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 informationFACULTY 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 informationLecture 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 informationInternal 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 informationFluid 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 informationIntroduction 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 informationSection 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 informations 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 informationFE 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 informationPrinciples 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 informationChapter 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 information1-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 information2.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 informationMath 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 informationME 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 informationHydraulics 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 informationEngineers 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 informationMOMENTS 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 informationReview: 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 informationELEC9712 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 informationApplied 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 informationChapter 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 informationVORTEX 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 informationTALLINN 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 informationInaugural 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 informationFluid 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 informationPROBLEMS. 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 informationChemical 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 informationExam 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 informationOscillatory 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 informationHarmonic 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 informationOE4625 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 informationI. 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 informationWhat 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 informationLab 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 informationChapter 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 informationAircraft 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 informationChapter 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 informationChapter 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 informationLAB 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 informationThursday, 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 information1.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 informationChapter 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 informationChapter 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)
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 information2. 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 informationa. 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 informationThe... 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 informationCh 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 informationBIOEN 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 informationMEMORANDUM. 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 informationFluid 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 informationUNIVERSITY 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 informationLab 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 informationExternal 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 informationSolutions 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 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 informationChapter 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 informationME8230 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 informationC 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 informationHOW 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 informationMass 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 informationspring 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 information6. 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 informationExperiment 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 informationSteven 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 informationMeasuring 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 informationVibration 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 informationLearning 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