Three-Tank Experiment

Similar documents
Lab Project 4a MATLAB Model of Oscillatory Flow

Dynamic Modeling. For the mechanical translational system shown in Figure 1, determine a set of first order

Modeling Basics: 4. Numerical ODE Solving In Excel 5. Solving ODEs in Mathematica

3.2 Systems of Two First Order Linear DE s

Study fluid dynamics. Understanding Bernoulli s Equation.

Fluid Flow Analysis Penn State Chemical Engineering

Introduction to Chemical Engineering Computing, 2 nd edition, Bruce A. Finlayson, Wiley (2012)

Basic Fluid Mechanics

LOSSES DUE TO PIPE FITTINGS

Linear Transport Relations (LTR)

M E 320 Professor John M. Cimbala Lecture 38

Lesson 6 Review of fundamentals: Fluid flow

ENGI 3424 First Order ODEs Page 1-01

ξ Da nc n = c Pe ξ = c o ( θ ) (95) c ξ = 0 (96) θ = 0 ; c= c i ξ ( ) (97) (98) Pe z = u L E z

Project TOUCAN. A Study of a Two-Can System. Prof. R.G. Longoria Update Fall ME 144L Prof. R.G. Longoria Dynamic Systems and Controls Laboratory

Mass of fluid leaving per unit time

Module 2: External Flows Lecture 6: Exact Solution of Boundary Layer Equation. The Lecture Contains: Blasius Solution. Numerical Approach

Process Control and Instrumentation Prof. A. K. Jana Department of Chemical Engineering Indian Institute of Technology, Kharagpur

Nodalization. The student should be able to develop, with justification, a node-link diagram given a thermalhydraulic system.

MTE 119 STATICS FINAL HELP SESSION REVIEW PROBLEMS PAGE 1 9 NAME & ID DATE. Example Problem P.1

6 Two-layer shallow water theory.

Lab 3A: Modeling and Experimentation: Two-Can System

Modeling and Experimentation: Compound Pendulum

The Enigma of The Transition to Turbulence in a Pipe

Lab-Report Control Engineering. Proportional Control of a Liquid Level System

Applied Fluid Mechanics

CHAPTER 80 NUMERICAL METHODS FOR FIRST-ORDER DIFFERENTIAL EQUATIONS

BERNOULLI EQUATION. The motion of a fluid is usually extremely complex.

ERRATA CRANE FLOW OF FLUIDS THROUGH VALVES, FITTINGS AND PIPE TECHNICAL PAPER NO. 410 METRIC VERSION

Solution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0

Objectives. Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation

M E 320 Supplementary Material Pralav Shetty

7.6 Example von Kármán s Laminar Boundary Layer Problem

Hydraulic (Fluid) Systems

Chemical Engineering 140. Chemical Process Analysis C.J. Radke Tentative Schedule Fall 2013

Enhancing Computer-Based Problem Solving Skills with a Combination of Software Packages

Calibrate Rotameter and Orifice Meter and Explore Reynolds #

Lecture 22. Mechanical Energy Balance

Chapter Four fluid flow mass, energy, Bernoulli and momentum

Q1 Give answers to all of the following questions (5 marks each):

LEAKLESS COOLING SYSTEM V.2 PRESSURE DROP CALCULATIONS AND ASSUMPTIONS

Ordinary Differential Equations. Monday, October 10, 11

M E 320 Professor John M. Cimbala Lecture 05

Chapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers

FLUID MECHANICS PROF. DR. METİN GÜNER COMPILER

Extended MHD simulation of Rayleigh-Taylor/Kelvin-Helmholtz instability

Use of Differential Equations In Modeling and Simulation of CSTR

LECTURE 1 THE CONTENTS OF THIS LECTURE ARE AS FOLLOWS:

PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Abstract

Math 307 E - Summer 2011 Pactice Mid-Term Exam June 18, Total 60

1 Systems of First Order IVP

Introduction to Fluid Flow

Simulation of Supercritical Fluid Extraction Process Using COMSOL Multiphysics

MTH 103 Group Activity Problems (W1B) Name: Types of Functions and Their Rates of Change Section 1.4 part 1 (due April 6)

MASS, MOMENTUM, AND ENERGY EQUATIONS

Math 128A Spring 2003 Week 12 Solutions

Atmospheric pressure. 9 ft. 6 ft

Physics 3 Summer 1990 Lab 7 - Hydrodynamics

Numerical modelling of phase change processes in clouds. Challenges and Approaches. Martin Reitzle Bernard Weigand

Conservation of Angular Momentum

EGN 3353C Fluid Mechanics

Lecture 24. Design of flow meters

FRICTION LOSS ALONG A PIPE

Math 122 Fall Handout 11: Summary of Euler s Method, Slope Fields and Symbolic Solutions of Differential Equations

Power System Transient Stability Simulation

Computational Physics (6810): Session 8

FLUID MECHANICS D203 SAE SOLUTIONS TUTORIAL 2 APPLICATIONS OF BERNOULLI SELF ASSESSMENT EXERCISE 1

y = sin(x) y = x x = 0 x = 1.

Numerical method for approximating the solution of an IVP. Euler Algorithm (the simplest approximation method)

Chapter 1: Introduction

MATHCAD SOLUTIONS TO THE CHEMICAL ENGINEERING PROBLEM SET

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016

Mixture of ideal gas and radiation: pressure effects

In steady flow the velocity of the fluid particles at any point is constant as time passes.

Chapter (6) Energy Equation and Its Applications

FORMULA SHEET. General formulas:

External Flows. Dye streak. turbulent. laminar transition

Differential Equations

Section 2.1 Differential Equation and Solutions

Math 113 (Calculus II) Final Exam KEY

Solution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0

Exploring Nonlinear Oscillator Models for the Auditory Periphery

Dynamic Loads CE 543. Examples. Harmonic Loads

BAE 820 Physical Principles of Environmental Systems

1. Consider the initial value problem: find y(t) such that. y = y 2 t, y(0) = 1.

ME332 FLUID MECHANICS LABORATORY (PART I)

A population is modeled by the differential equation

Chapter 15: Fluid Mechanics Dynamics Using Pascal s Law = F 1 = F 2 2 = F 2 A 2

Physics 123 Unit #1 Review


Lab 1: Damped, Driven Harmonic Oscillator

Empirical versus Direct sensitivity computation

Application of Reconstruction of Variational Iteration Method on the Laminar Flow in a Porous Cylinder with Regressing Walls

CLASS SCHEDULE 2013 FALL

Chapter (4) Motion of Fluid Particles and Streams

FDM for parabolic equations

Lab 1: damped, driven harmonic oscillator

Project #1 Internal flow with thermal convection

Chemical Engineering 374

Bernoulli and Pipe Flow

Transcription:

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 3/4" leading from the bottom of each tank to the other two tanks (see Figure 1). Two of the tanks are identical and have the diameter of 11". The diameter of the third tank is 13". The tanks will be filled with water to different heights, which will drive flow of water between the tanks. The purpose of this experiment is to investigate time dependence of water levels in the tanks as the system reaches equilibrium. The experimental observations should be compared with theoretical predictions which will be obtained by solving a system of nonlinear differential equations. Therefore, students should review numerical methods (e.g., Euler, Runge-Kutta) for integration of differential equations. Figure 1. Three-tank system. 1

Theory The governing equations in this experiment are obtained from applying force and mass balances on the fluid flow between tanks. Mass balance for the flow of water leaving one of the i-th tanks is: dh i dt = Q i A i, (1) where h i is the level of water in the i-th tank, Q i is the flow rate of water leaving or entering the i- th tank, and A i is the cross-sectional area of the i-th tank. A more useful relationship can be obtained by rewriting the flow rate in terms of the average fluid velocity v i of water in the pipe leading to the i-th tank, dh i dt = A p v A i, (2) i where A p is the pipe area. To obtain the average velocity, we need to apply the mechanical balance to the pipe system. For simplicity, let us first perform this analysis for a 2-tank system shown in Figure 2. After that, we will extend the analysis to the 3-tank system. Figure 2. Schematics of the two-tank system. Two-Tank System The mechanical balance between points 1 and 2 of the 2-tank system is P 1 P 2 ρ = 1 2 K 12v 2 12, (3) where K 12 is the friction coefficient of the pipe connecting the two tanks, v 12 is the fluid velocity in this pipe, P i is the hydrostatic pressure at the bottom of the i-th tank (i = 1, 2), and ρ is the water density. The hydrostatic pressure difference is Substituting this into Eq. (3) and solving for v 12 we get: P 1 P 2 = ρg(h 1 h 2 ) (4) 2

v 12 = 2g(h 1 h 2 ) K 12. (5) This equation assumes that h 1 h 2. Note that if the h 1 < h 2, Eq. (5) will yield an imaginary number for v 12, whereas physically it should be a negative value. To account for this, Eq. (5) should be rewritten as follows: where v ij = Sign(h i h j ) 2g h i h j K ij, (6) 1, x > 0 Sign(x) = { 0, x = 0 1, x < 0 is the sign function. The indices ij in Eq. (6) denote the location of the pipe network (for example, K 12, is the friction coefficient of the pipes between tanks 1 and 2). Eqs. (5) and (6) adhere to the following sign convention: v ij > 0 if water flows from tank i to tank to tank j. Now that we have an expression for the velocity between two tanks, we can plug it into the mass balance Eq. (2) to get a system of first-order differential equations describing the water level in each of the tanks. For simplicity, let us consider flow between two tanks of the same diameter. In this case, the mass balance simplifies to Integrating this equation, we obtain where (7) dh 1 dt = dh 2 dt. (8) h 1 (t) = h 2 (t) + C (9) C = h 1 (0) + h 2 (0) (10) is the integration constant related to the initial levels h 1 (0) and h 2 (0). Substituting Eq. (9) into Eq. (6) and then substituting the resulting expression into the mass balance Eq. (2) for tank 2 yields: dh 2 dt = A p A 2 v 12 = A p A 2 Sign(h 1 h 2 ) 2g C 2h 2 K 12 (11) Assuming that K 12 is independent of the fluid velocity, the differential equation (11) can be solved analytically to obtain an expression for dependence of h 2 on time. Plotting a linearized model of the result is useful to obtain the friction coefficient from the slope of the graph. This 3

friction coefficient can be used as an approximation in the three-tank system to simplify calculations. Three-Tank System In the three-tank system, water leaving one tank flows into two other tanks, depending on the levels of water in the tanks. Performing a mass balance on the tee below tank 1, we get: v = v 12 + v 13 (12) Since the area of the pipe is the same throughout the system, the flow rate reduces to just the average velocities. v 12 is the velocity of water going from tank 1 to tank 2 and similarly, v 13 is the velocity of water going from tank 1 to tank 3. Plugging in Eq. (6) for both v 12 and v 13 and then substituting Eq. (12) into our mass balance Eq. (2) for each of the tanks will yield three first-order differential equations. The friction coefficients in the pipe network are given by: K ij = K ij,f + K ij,s (13) K f is a constant and it is the summation of the friction contributions from the fittings, K s is the skin friction of the pipe. The skin friction coefficient for flow through a pipe is: K s = 4 L f, (14) D where L is the length of the pipe, D is the diameter of the pipe, and f is the Fanning friction factor. The flow between the tanks can be assumed to be turbulent for the entire duration of the experiment (see prelab questions). Therefore, we can use Blasius equation for the Fanning friction factor: f = 0.079Re.25 (15) Substituting Eq. (15) into Eq. (14) and then plugging in the friction coefficients into the system of equations will yield three non-linear first-order differential equations. Solution of this system requires a numerical solution using, e.g., Euler or Runge-Kutta method. The friction coefficients depend on the velocity of the system, which means a solver is needed to first obtain the velocity to obtain the friction coefficients. This must be done at all water levels and seeded in during the numerical solver process which is difficult to manage. Instead, friction coefficients obtained from the two-tank analysis can be used as an approximation. Since the pipe networks between tanks 1 and 2 and tanks 1 and 3 are different, two different coefficients must be obtained. 4

Objectives 1. Perform 2-tank flow experiments and obtain friction coefficients in pipes connecting each pair of tanks. 2. Perform 3-tank flow experiments at various initial water levels in the tanks. 3. Make a theoretical prediction for transient flows in the 3-tank system. 4. Compare the experimental results with the theoretical prediction. 5