Statistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders

Size: px
Start display at page:

Download "Statistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders"

Transcription

1 Statistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders Toshiki Morita 1, Takashi Watanabe 2 and Yorinobu Toya 3 1. Graduate School of Information Science, Nagoya University, Nagoya, Japan 2. EcoTopia Science Institute, Nagoya University, Nagoya, Japan 3. Nagano National College of Technology, Nagano, Japan Abstract: We performed numerical computations of Taylor vortex flow at small aspect ratio. When the flow develops from its initial state at rest, non-unique flow patterns appear, which depend on the lengths of the cylinders and/or the final Reynolds number. In experiment, it is difficult to keep constant conditions because of the imperfection of the apparatus and the small disturbances of the flow field. Therefore, the experimental result may show its inherit variations. In this study, we focus our attention on the small fluctuations of flows. Two-dimensional numerical simulation is used to estimate the effect on the final modes of Taylor vortex flow. Statistical analysis of the results presented a favorable agreement with our prescribed experimental result. Keywords: Bifurcation, Mode Formation, Taylor Vortex Flow, Computational Fluid Dynamics 1. INTRODUCTION The study of the transition of the flow is important for understanding of the non-linear dynamics. Taylor vortex flow between a rotating inner cylinder and a fixed outer cylinder is one of the models which show non-uniqueness characteristic in the non-linear dynamics. As for Taylor vortex flow, it is experimentally clear that there are some stable states [1], [2]. When a flow develops from an initial state of rest, a non-uniqueness of patterns of various flows appears, which depends on the length of the cylinders and final Reynolds number. In the physical experiment, it is difficult to perform experiments with fixed conditions because of the incompleteness of the apparatus and minute fluctuations of the flow field. In fact, the experimental results showed that formation processes and final modes varied statistically [3], [4]. Numerical studies that focus on minute fluctuations of the flow have been carried out. Shin and Reed [5] used random initial values to investigate decaying properties of a nonlinear unsteady Burgers equation. In the case of Taylor vortex flow, Hirshfeld and Rapaport [6] presented that initial random velocities affected the growth of Taylor vortices. Rigopoulos et. al. [7] investigated new mode vortices caused by the Eckhaus instability. In their study, circular Couette flows with random perturbations were used as initial conditions. Bilson and Bremhorst [8] presented that final modes of turbulent Taylor-Couette flows were affected by the acceleration ratio and the velocity fluctuation field. However, these initial random values do not satisfy physical conditions. This study pays its attention to the minute fluctuations of the fluid. We use the initial velocity field that have minute fluctuations and meet the physical conditions, and investigate the influence on the final modes of Taylor vortex flow by two-dimensional numerical simulation. We compare final modes of the conventional experimental results with the calculation results. 2. COMPUTATIONAL METHOD 2.1. Dominant equations and formulation In this study, we use a model to show in Fig. 1. The dynamic parameters in the numerical analysis are Reynolds number Re and the acceleration of the inner cylinder, and the geometric parameters are the aspect ratio Γ and the radius ratio of the inner cylinder and the outer cylinder η. The Reynolds number is defined as Re=ωR 1 d/ν by the angular velocity ω of the inner cylinder, the gap d of the inner and outer cylinder, the radius R 1 of the inner cylinder and the kinematic viscosity ν of the working fluid. The aspect ratio Γ is given as Γ =L/d by the height of the working fluid L. The radius ratio η is set to be It is assumed that the inner cylinder turns and outer cylinder is fixed. The aspect ratio is 4.0, and the upper and lower ends of the cylinders are fixed. The Reynolds numbers are set from 1000 to 6000 with an interval of Table 1 shows the accelerations of the inner cylinder. The final modes are computed in thirty-six combinations of conditions in total. The basic equations are the Navier-Stokes equations in the cylindrical coordinate system (r,θ,z) and the equation of continuity : u + t 1 Re ( uu) = p + u, u = 0. We use the following function ψ for visualizing: 1 ψ 1 ψ u =, w =. r z r r This function is similar to the stream function of Stokes in the axisymmetry flow. We apply the no-slip condition at the upper and lower walls and the inner and outer walls. The 1101

2 calculation conforms to MAC method. The calculation grid is staggerd grid. The Poisson equation of the pressure is solved in SOR. The time development of the velocity use the fractional step method Definition of the modes The modes of the Taylor vortex flow are decided by the direction of the flow near the end walls, and they are a normal mode and an anomalous mode. The normal mode has a structure where flows through the outer cylinder from the inner cylinder on both end walls in which fluid flows from the inner cylinder to the outer cylinder. The anomalous mode has a structure in which fluid flows from the inner cylinder to the outer cylinder on both end walls or one end wall. The symbol Nn represents normal n cell mode and An is anomalous n cell mode, where integer n stands for the number of cells Generation of initial velocity field First, we rotate the inner cylinder and generate the normal 4-cell mode or anomalous 4-cell mode. Then we completely stop the inner cylinder and take out the flow at the time when the deviations of disturbances relative to the original peripheral velocity of the inner cylinder are 10-5 ~ 10-8 during deceleration. Figure 2 shows contours of the stream function at the initial flow field. The left sides of the contours are the inner cylinders and the right sides are outer cylinders. The starting points of the allows are the contours before deceleration and the ending points of the allows are the contours after deceleration. Fig. 2 shows that original and do not seen. We use each initial flow field and compute well-developed flow and get final modes. Table 1. Definition of the symbol of the angular acceleration Re Symbol Re = 1000 ~ 6000 Angular acceleration of Acceleration AC1 15 rad/s 2 AC2 30 rad/s 2 AC3 45 rad/s 2 AC4 60 rad/s 2 AC5 75 rad/s 2 AC6 90 rad/s 2 Fig.2 Initial flow field at and. Fig.3 Bifurcation with decrease of the Reynolds number (Γ=4.0). Al3 Au3 Fig.1 Cylindrical coordinate system Al5 Au5 Fig.4 Final modes of Taylor vortex flow at Γ=

3 (a) Numerical result (decelerated n4). (a) Numerical result (decelerated n4). Fig.6 The probability of each modes at Γ=4.0, Re=1000. Fig.8 The probability of each modes at Γ=4.0, Re=3000. (a) Numerical result (decelerated ). (a) Numerical result (decelerated ). Fig.7 The probability of each modes at Γ=4.0, Re=2000. Fig.9 The probability of each modes at Γ=4.0, Re=

4 (a) Numerical result (Initial n4). Fig.10 The probability of each modes at Γ=4.0, Re=5000. (a) Numerical result (Initial n4). 3. RESULTS 3.1. Numerical results of flow mode Figure 3 shows the experimental results of the bifurcation of the Taylor vortex flow of experiment at Γ=4.0. The arrow shows the bifurcated direction when the Reynolds number decreases slowly. Six kinds of vortex structures appeared at Γ=4.0:,,,, and. Figure 4 shows the numerical results about contours of the stream function at Γ=4.0. In this figure, the arrows denote the directions of the flow. Figure 4 includes whole modes found in Fig. 3. The minute disturbances of initial velocity field reproduce experimental modes The probability of each modes Figures from Fig. 6 (a) to Fig. 11 (a) show the probability of each mode obtained from numerical computations at Γ = 4.0, Re= Vertical axis shows the probability of each mode in percentage and horizontal axis shows the acceleration of the inner cylinder. The initial flow fields are made by deceleration the flow of. In case the initial flow fields are generated from the flow of, no significant differences is found. Figures from Fig.6 (b) to Fig.11 (b) show the results that are same conditions as those of experiments. The mod of and are classified into to patterns according whether anomalous cell is on upper or lower end walls. However, the figures from Fig.6 (a) to Fig.11 (a) don t classify the modes to adapt the experiment. Figure 6 shows only is found in the results of experiment at Re=2000, but, and are found in those of numerical computation. At Re=2000 in Fig.7, the results of numerical computation generate more probability of than those of experiment. At Re=3000 in Fig. 8, both of the results generate high probability of. Figure 9 shows that the probability of increases. When the Reynolds number is 5000 in Fig. 10,, and are often generated and other modes are hardly gotten. Figure 11 shows, and are often generated in the results of calculation and experiment. In all Reynolds numbers, the final modes appear in the experiment are gotten by the present calculation. 4. CONCLUSIONS We give some pseudo-disturbances to the initial velocity field and find influences on the final modes of the Taylor vortex flow by the two-dimensional numerical computation. We get all final modes shown in bifurcation diagram by the minute difference of the initial flow field in the numerical computation. We reproduce the prior experiment by the present calculation. The calculation results accord with the experiment results qualitatively because the modes appeared in the experiment are made by the numerical computation. Fig.11 The probability of each modes at Γ=4.0, Re=

5 REFERNCE 1. I. Nakamura et al., Japan Society of Mechanical Engineers, , B(1990), p Y. Toya et al., Japan Society of Mechanical Engineers, , B(1990), p I. Nakamura et al., Japan Society of Mechanical Engineers, , B(1994), p J. E. Burkhalter and E. L. Koschmieder, Phys. Fluids, vol. 17 (1974), p Y-C Shin and X B Reed, Jr., Phys. Fluids, vol. 28 (1985), p D. Hirshfeld and D. C. Rapaport, Phys. Review E, vol. 61 (1999), p J. Rigopoulos, J. Sheridan and M. C. Thompson, J. Fluid Mech, vol. 489 (2003), p M. Bilson and K. Bremhorst, J. Fluid Mech., vol. 579 (2006), p

Study on Non-Uniqueness of Taylor Vortex Flow Changing Inner Cylinder Acceleration Time

Study on Non-Uniqueness of Taylor Vortex Flow Changing Inner Cylinder Acceleration Time World Journal of Mechanics, 2018, 8, 301-310 http://www.scirp.org/journal/wjm ISSN Online: 2160-0503 ISSN Print: 2160-049X Study on Non-Uniqueness of Taylor Vortex Flow Changing Inner Cylinder Acceleration

More information

Development of free surface flow between concentric cylinders with vertical axes

Development of free surface flow between concentric cylinders with vertical axes Journal of Physics: Conference Series Development of free surface flow between concentric cylinders with vertical axes To cite this article: T Watanabe et al 005 J. Phys.: Conf. Ser. 14 9 View the article

More information

Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions

Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions Chapter 1 Earth Science Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions Project Representative Yozo Hamano Authors Ataru Sakuraba Yusuke Oishi

More information

Introduction. Statement of Problem. The governing equations for porous materials with Darcy s law can be written in dimensionless form as:

Introduction. Statement of Problem. The governing equations for porous materials with Darcy s law can be written in dimensionless form as: Symbolic Calculation of Free Convection for Porous Material of Quadratic Heat Generation in a Circular Cavity Kamyar Mansour Amirkabir University of technology, Tehran, Iran, 15875-4413 mansour@aut.ac.ir

More information

Evidence for Internal Structures of Spiral Turbulence. Abstract

Evidence for Internal Structures of Spiral Turbulence. Abstract Evidence for Internal Structures of Spiral Turbulence S. Dong Center for Computational & Applied Mathematics, Department of Mathematics, Purdue University, West Lafayette, IN 47907 Abstract We report the

More information

Numerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder

Numerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder Numerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder Ali Kianifar, Edris Yousefi Rad Abstract In many applications the flow that past bluff bodies have frequency nature (oscillated)

More information

Numerical Simulation of Rotating Turbulent Channel Flow by the Vortex in Cell Method

Numerical Simulation of Rotating Turbulent Channel Flow by the Vortex in Cell Method Send Orders for Reprints to reprints@benthamscience.net 30 The Open Transport Phenomena Journal, 2013, 5, 30-41 Open Access Numerical Simulation of Rotating Turbulent Channel Flow by the Vortex in Cell

More information

Self-Excited Vibration in Hydraulic Ball Check Valve

Self-Excited Vibration in Hydraulic Ball Check Valve Self-Excited Vibration in Hydraulic Ball Check Valve L. Grinis, V. Haslavsky, U. Tzadka Abstract This paper describes an experimental, theoretical model and numerical study of concentrated vortex flow

More information

Chuichi Arakawa Graduate School of Interdisciplinary Information Studies, the University of Tokyo. Chuichi Arakawa

Chuichi Arakawa Graduate School of Interdisciplinary Information Studies, the University of Tokyo. Chuichi Arakawa Direct Numerical Simulations of Fundamental Turbulent Flows with the Largest Grid Numbers in the World and its Application of Modeling for Engineering Turbulent Flows Project Representative Chuichi Arakawa

More information

Candidates must show on each answer book the type of calculator used. Only calculators permitted under UEA Regulations may be used.

Candidates must show on each answer book the type of calculator used. Only calculators permitted under UEA Regulations may be used. UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2011 2012 FLUID DYNAMICS MTH-3D41 Time allowed: 3 hours Attempt FIVE questions. Candidates must show on each answer book the type

More information

Validation 3. Laminar Flow Around a Circular Cylinder

Validation 3. Laminar Flow Around a Circular Cylinder Validation 3. Laminar Flow Around a Circular Cylinder 3.1 Introduction Steady and unsteady laminar flow behind a circular cylinder, representing flow around bluff bodies, has been subjected to numerous

More information

1. Introduction, tensors, kinematics

1. Introduction, tensors, kinematics 1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and

More information

Visualization of Traveling Vortices in the Boundary Layer on a Rotating Disk under Orbital Motion

Visualization of Traveling Vortices in the Boundary Layer on a Rotating Disk under Orbital Motion Open Journal of Fluid Dynamics, 2015, 5, 17-25 Published Online March 2015 in SciRes. http://www.scirp.org/journal/ojfd http://dx.doi.org/10.4236/ojfd.2015.51003 Visualization of Traveling Vortices in

More information

Application of a Virtual-Boundary Method for the Numerical Study of Oscillations Developing Behind a Cylinder Near A Plane Wall

Application of a Virtual-Boundary Method for the Numerical Study of Oscillations Developing Behind a Cylinder Near A Plane Wall Fluid Dynamics, Vol. 39, No. 1, 2004, pp. 61 68. Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, 2004, pp. 69 77. Original Russian Text Copyright 2004 by Kit, Nikitin,

More information

FLOW-NORDITA Spring School on Turbulent Boundary Layers1

FLOW-NORDITA Spring School on Turbulent Boundary Layers1 Jonathan F. Morrison, Ati Sharma Department of Aeronautics Imperial College, London & Beverley J. McKeon Graduate Aeronautical Laboratories, California Institute Technology FLOW-NORDITA Spring School on

More information

Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders

Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute

More information

The Reynolds experiment

The Reynolds experiment Chapter 13 The Reynolds experiment 13.1 Laminar and turbulent flows Let us consider a horizontal pipe of circular section of infinite extension subject to a constant pressure gradient (see section [10.4]).

More information

Numerical Investigations on a Curved Pipe Flow

Numerical Investigations on a Curved Pipe Flow Vienna Unversity of Technology Institute of Fluid Mechanics and Heat Transfer Numerical Investigations on a Curved Pipe Flow Bachelor Thesis by Martin Schwab Supervisors Univ.Prof. Dipl.-Phys. Dr.rer.nat.

More information

arxiv:physics/ v1 14 Nov 2006

arxiv:physics/ v1 14 Nov 2006 Controlling the stability transfer between oppositely traveling waves and standing waves by inversion symmetry breaking perturbations A. Pinter, M. Lücke, and Ch. Hoffmann Institut für Theoretische Physik,

More information

The Turbulent Rotational Phase Separator

The Turbulent Rotational Phase Separator The Turbulent Rotational Phase Separator J.G.M. Kuerten and B.P.M. van Esch Dept. of Mechanical Engineering, Technische Universiteit Eindhoven, The Netherlands j.g.m.kuerten@tue.nl Summary. The Rotational

More information

Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5

Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Jingwei Zhu May 14, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of

More information

Physical Properties of Fluids

Physical Properties of Fluids Physical Properties of Fluids Viscosity: Resistance to relative motion between adjacent layers of fluid. Dynamic Viscosity:generally represented as µ. A flat plate moved slowly with a velocity V parallel

More information

Simulation of Aeroelastic System with Aerodynamic Nonlinearity

Simulation of Aeroelastic System with Aerodynamic Nonlinearity Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School

More information

F11AE1 1. C = ρν r r. r u z r

F11AE1 1. C = ρν r r. r u z r F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions

More information

Boundary-Layer Theory

Boundary-Layer Theory Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22

More information

Title. Author(s)Suzuki, Toshiyuki; Iima, Makoto; Hayase, Yumino. CitationPhysics of Fluids, 18(10): Issue Date Doc URL. Rights.

Title. Author(s)Suzuki, Toshiyuki; Iima, Makoto; Hayase, Yumino. CitationPhysics of Fluids, 18(10): Issue Date Doc URL. Rights. Title Surface switching of rotating fluid in a cylinder Author(s)Suzuki, Toshiyuki; Iima, Makoto; Hayase, Yumino CitationPhysics of Fluids, 18(10): 101701 Issue Date 2006-10 Doc URL http://hdl.handle.net/2115/16007

More information

Patterns of Turbulence. Dwight Barkley and Laurette Tuckerman

Patterns of Turbulence. Dwight Barkley and Laurette Tuckerman Patterns of Turbulence Dwight Barkley and Laurette Tuckerman Plane Couette Flow Re = U gap/2 ν Experiments by Prigent and Dauchot Re400 Z (Spanwise) 40 24 o Gap 2 Length 770 X (Streamwise) Examples: Patterns

More information

Numerical Simulation and Its Possibility of Application to Green Environment

Numerical Simulation and Its Possibility of Application to Green Environment Numerical Simulation and Its Possibility of Application to Green Environment Anna Kuwana 1,a, Tetuya Kawamura 2,b 1 Faculty of Science and Technology, Gunma University, 1-5-1 Tenjin-cho, Kiryu City, Gunma

More information

Lagrangian acceleration in confined 2d turbulent flow

Lagrangian acceleration in confined 2d turbulent flow Lagrangian acceleration in confined 2d turbulent flow Kai Schneider 1 1 Benjamin Kadoch, Wouter Bos & Marie Farge 3 1 CMI, Université Aix-Marseille, France 2 LMFA, Ecole Centrale, Lyon, France 3 LMD, Ecole

More information

REE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics

REE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics REE 307 - Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics 1. Is the following flows physically possible, that is, satisfy the continuity equation? Substitute the expressions for

More information

Table of Contents. Foreword... xiii. Preface... xv

Table of Contents. Foreword... xiii. Preface... xv Table of Contents Foreword.... xiii Preface... xv Chapter 1. Fundamental Equations, Dimensionless Numbers... 1 1.1. Fundamental equations... 1 1.1.1. Local equations... 1 1.1.2. Integral conservation equations...

More information

On Taylor-Couette flow: Inferences of near-wall effect

On Taylor-Couette flow: Inferences of near-wall effect On Taylor-Couette flow: Inferences of near-wall effect Md Abdur Razzak, Kim Boon Lua 2, Yin Jen Lee, Boo Cheong Khoo Department of Mechanical Engineering, National University of Singapore,Singapore, 2

More information

Suppression of 3D flow instabilities in tightly packed tube bundles

Suppression of 3D flow instabilities in tightly packed tube bundles Suppression of 3D flow instabilities in tightly packed tube bundles Nicholas Kevlahan kevlahan@mcmaster.ca Department of Mathematics & Statistics CSFD, June 13 15 2004 p.1/33 Collaborators CSFD, June 13

More information

On Decaying Two-Dimensional Turbulence in a Circular Container

On Decaying Two-Dimensional Turbulence in a Circular Container Frontiers of Computational Sciences Y. Kaneda, H. Kawamura and M. Sasai (Eds.) Springer, 2007, pp. 89-95 On Decaying Two-Dimensional Turbulence in a Circular Container Kai Schneider and Marie Farge Univesité

More information

Energy dissipating structures generated by dipole-wall collisions at high Reynolds number

Energy dissipating structures generated by dipole-wall collisions at high Reynolds number Energy dissipating structures generated by dipole-wall collisions at high Reynolds number Duncan Sutherland 1 Charlie Macaskill 1 David Dritschel 2 1. School of Mathematics and Statistics University of

More information

Kinematic and dynamic pair collision statistics of sedimenting inertial particles relevant to warm rain initiation

Kinematic and dynamic pair collision statistics of sedimenting inertial particles relevant to warm rain initiation Kinematic and dynamic pair collision statistics of sedimenting inertial particles relevant to warm rain initiation Bogdan Rosa 1, Hossein Parishani 2, Orlando Ayala 2, Lian-Ping Wang 2 & Wojciech W. Grabowski

More information

Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition

Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow

More information

FEDSM PRESSURE DISTRIBUTION ON THE GROUND BY IMPINGING TWO-DIMENSIONAL JET DUE TO A VORTEX METHOD

FEDSM PRESSURE DISTRIBUTION ON THE GROUND BY IMPINGING TWO-DIMENSIONAL JET DUE TO A VORTEX METHOD Proceedings of the 3rd SME/JSME Joint Fluids Engineering Conferences July 8-23, 999, San Francisco, California, US FEDSM99-685 PRESSURE DISTRIBUTION ON THE GROUND BY IMPINGING TWO-DIMENSIONL JET DUE TO

More information

Side-View Mirror Vibrations Induced Aerodynamically by Separating Vortices

Side-View Mirror Vibrations Induced Aerodynamically by Separating Vortices Open Journal of Fluid Dynamics, 2016, 6, 42-56 Published Online March 2016 in SciRes. http://www.scirp.org/journal/ojfd http://dx.doi.org/10.4236/ojfd.2016.61004 Side-View Mirror Vibrations Induced Aerodynamically

More information

Transition to turbulence in plane Poiseuille flow

Transition to turbulence in plane Poiseuille flow Proceedings of the 55th Israel Annual Conference on Aerospace Sciences, Tel-Aviv & Haifa, Israel, February 25-26, 2015 ThL2T5.1 Transition to turbulence in plane Poiseuille flow F. Roizner, M. Karp and

More information

Experimental setup. Chapter Rheometer

Experimental setup. Chapter Rheometer 21 Chapter 2 Experimental setup The current experiments are designed to examine the effect of volume fraction and Stokes number (and equivalently the Reynolds number) at shear rates sufficiently high enough

More information

Numerical study of multiple periodic flow states in spherical Couette flow

Numerical study of multiple periodic flow states in spherical Couette flow Numerical study of multiple periodic flow states in spherical Couette flow YUAN Li ICMSEC & and LSEC, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100080, CHINA Abstract

More information

SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS

SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS ICAS 2000 CONGRESS SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS H Yao, R K Cooper, and S Raghunathan School of Aeronautical Engineering The Queen s University of Belfast, Belfast BT7 1NN,

More information

Numerical Investigation of Laminar Flow over a Rotating Circular Cylinder

Numerical Investigation of Laminar Flow over a Rotating Circular Cylinder International Journal of Mechanical & Mechatronics Engineering IJMME-IJENS Vol:13 No:3 32 Numerical Investigation of Laminar Flow over a Rotating Circular Cylinder Ressan Faris Al-Maliky Department of

More information

Applied Thermal and Fluid Engineering. Energy Engineering (Thermal Engineering Laboratory)

Applied Thermal and Fluid Engineering. Energy Engineering (Thermal Engineering Laboratory) Applied Thermal and Fluid Engineering Energy Engineering (Thermal Engineering Laboratory) Professor Assoc. Professor Hajime Nakamura Shunsuke Yamada Outline of Research In our laboratory, we have been

More information

Masters in Mechanical Engineering Aerodynamics 1 st Semester 2015/16

Masters in Mechanical Engineering Aerodynamics 1 st Semester 2015/16 Masters in Mechanical Engineering Aerodynamics st Semester 05/6 Exam st season, 8 January 06 Name : Time : 8:30 Number: Duration : 3 hours st Part : No textbooks/notes allowed nd Part : Textbooks allowed

More information

Intermittency in spiral Poiseuille flow

Intermittency in spiral Poiseuille flow Intermittency in spiral Poiseuille flow M. Heise, J. Abshagen, A. Menck, G. Pfister Institute of Experimental and Applied Physics, University of Kiel, 2498 Kiel, Germany E-mail: heise@physik.uni-kiel.de

More information

INFLUENCE OF GEOMETRICAL ASPECT RATIO ON THE OSCILLATORY MARANGONI CONVECTION IN LIQUID BRIDGES R. MONTI, R. SAVINO, M. LAPPA *

INFLUENCE OF GEOMETRICAL ASPECT RATIO ON THE OSCILLATORY MARANGONI CONVECTION IN LIQUID BRIDGES R. MONTI, R. SAVINO, M. LAPPA * AUTHOR S POST PRINT (Romeo Colour: Green) Acta Astronautica (ISSN: 94-5765) Vol. 47, No., pp. 753-76, 2. DOI.6/S94-5765()26- Publisher version available at http://www.sciencedirect.com/science/article/pii/s94576526

More information

Axisymmetric Hopf bifurcation in a free surface rotating cylinder flow

Axisymmetric Hopf bifurcation in a free surface rotating cylinder flow ANZIAM J. 50 (CTAC2008) pp.c251 C265, 2008 C251 Axisymmetric Hopf bifurcation in a free surface rotating cylinder flow S. J. Cogan 1 G. J. Sheard 2 K. Ryan 3 (Received 13 August 2008; revised 24 October

More information

Effect of Liquid Viscosity on Sloshing in A Rectangular Tank

Effect of Liquid Viscosity on Sloshing in A Rectangular Tank International Journal of Research in Engineering and Science (IJRES) ISSN (Online): 2320-9364, ISSN (Print): 2320-9356 Volume 5 Issue 8 ǁ August. 2017 ǁ PP. 32-39 Effect of Liquid Viscosity on Sloshing

More information

Periodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain.

Periodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain. Flow Characterization of Inclined Jet in Cross Flow for Thin Film Cooling via Large Eddy Simulation Naqavi, I.Z. 1, Savory, E. 2 and Martinuzzi, R. J. 3 1,2 The Univ. of Western Ontario, Dept. of Mech.

More information

STEADY, UNSTEADY AND LINEAR STABILITY OF FLOW PAST AN ELLIPTIC CYLINDER S.J.D. D'ALESSIO

STEADY, UNSTEADY AND LINEAR STABILITY OF FLOW PAST AN ELLIPTIC CYLINDER S.J.D. D'ALESSIO CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 4, Number 4, Fall 1096 STEADY, UNSTEADY AND LINEAR STABILITY OF FLOW PAST AN ELLIPTIC CYLINDER S.J.D. D'ALESSIO ABSTRACT. Diecud in thin work is the twedimensional

More information

ME332 FLUID MECHANICS LABORATORY (PART I)

ME332 FLUID MECHANICS LABORATORY (PART I) ME332 FLUID MECHANICS LABORATORY (PART I) Mihir Sen Department of Aerospace and Mechanical Engineering University of Notre Dame Notre Dame, IN 46556 Version: January 14, 2002 Contents Unit 1: Hydrostatics

More information

EULER AND SECOND-ORDER RUNGE-KUTTA METHODS FOR COMPUTATION OF FLOW AROUND A CYLINDER

EULER AND SECOND-ORDER RUNGE-KUTTA METHODS FOR COMPUTATION OF FLOW AROUND A CYLINDER EULER AND SEOND-ORDER RUNGE-KUTTA METHODS FOR OMPUTATION OF FLOW AROUND A YLINDER László Daróczy, László Baranyi MSc student, Professor University of Miskolc, Hungary Department of Fluid and Heat Engineering,

More information

Fluid Dynamics Exercises and questions for the course

Fluid Dynamics Exercises and questions for the course Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r

More information

Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13

Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13 Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 212/13 Exam 2ª época, 2 February 213 Name : Time : 8: Number: Duration : 3 hours 1 st Part : No textbooks/notes allowed 2 nd Part :

More information

Lecture 19. ROLLING WITHOUT SLIPPING. Figure 4.8 Gear rolling in geared horizontal guides. Figure 4.9 Wheel rolling on a horizontal surface.

Lecture 19. ROLLING WITHOUT SLIPPING. Figure 4.8 Gear rolling in geared horizontal guides. Figure 4.9 Wheel rolling on a horizontal surface. Lecture 19. ROLLING WITHOUT SLIPPING Figure 4.8 Gear rolling in geared horizontal guides. Figure 4.9 Wheel rolling on a horizontal surface. The derivation and understanding of velocity and acceleration

More information

NOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

NOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7 Number Pages 3 39 c Institute for Scientific Computing and Information NOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL

More information

Numerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus

Numerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus Nonlinear Analysis: Modelling and Control, 2007, Vol. 12, No. 4, 541 552 Numerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus M. Sankar

More information

DIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING

DIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING DIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING Yukinori Kametani Department of mechanical engineering Keio

More information

Applied Mathematics and Mechanics (English Edition)

Applied Mathematics and Mechanics (English Edition) Appl. Math. Mech. -Engl. Ed., 39(9), 1267 1276 (2018) Applied Mathematics and Mechanics (English Edition) https://doi.org/10.1007/s10483-018-2364-7 Direct numerical simulation of turbulent flows through

More information

Thomas Pierro, Donald Slinn, Kraig Winters

Thomas Pierro, Donald Slinn, Kraig Winters Thomas Pierro, Donald Slinn, Kraig Winters Department of Ocean Engineering, Florida Atlantic University, Boca Raton, Florida Applied Physics Laboratory, University of Washington, Seattle, Washington Supported

More information

Decaying 2D Turbulence in Bounded Domains: Influence of the Geometry

Decaying 2D Turbulence in Bounded Domains: Influence of the Geometry Computational Physics and New Perspectives in Turbulence Y. Kaneda (Ed.) Springer, 2007, pp. 241-246 Decaying 2D Turbulence in Bounded Domains: Influence of the Geometry Kai Schneider 1 and Marie Farge

More information

Exercise: concepts from chapter 10

Exercise: concepts from chapter 10 Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like

More information

Transition processes for junction vortex flow

Transition processes for junction vortex flow J. Fluid Mech. (7), vol. 585, pp. 57 67. c 7 Cambridge University Press doi:.7/s7685 Printed in the United Kingdom 57 Transition processes for junction vortex flow J. J. ALLEN AND J. M. L O P E Z Department

More information

UNIVERSITY OF EAST ANGLIA

UNIVERSITY OF EAST ANGLIA UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must

More information

AC & DC Magnetic Levitation and Semi-Levitation Modelling

AC & DC Magnetic Levitation and Semi-Levitation Modelling International Scientific Colloquium Modelling for Electromagnetic Processing Hannover, March 24-26, 2003 AC & DC Magnetic Levitation and Semi-Levitation Modelling V. Bojarevics, K. Pericleous Abstract

More information

The JHU Turbulence Databases (JHTDB)

The JHU Turbulence Databases (JHTDB) The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW AT Re τ = 5200 DATA SET Data provenance: M. Lee 1 & R. D. Moser 1 Database ingest and Web Services: Z. Wu 2, G. Lemson 2, R. Burns 2, A. Szalay

More information

CVS filtering to study turbulent mixing

CVS filtering to study turbulent mixing CVS filtering to study turbulent mixing Marie Farge, LMD-CNRS, ENS, Paris Kai Schneider, CMI, Université de Provence, Marseille Carsten Beta, LMD-CNRS, ENS, Paris Jori Ruppert-Felsot, LMD-CNRS, ENS, Paris

More information

Vortex Induced Vibrations

Vortex Induced Vibrations Vortex Induced Vibrations By: Abhiroop Jayanthi Indian Institute of Technology, Delhi Some Questions! What is VIV? What are the details of a steady approach flow past a stationary cylinder? How and why

More information

arxiv: v1 [physics.flu-dyn] 21 Jan 2015

arxiv: v1 [physics.flu-dyn] 21 Jan 2015 January 2015 arxiv:1501.05620v1 [physics.flu-dyn] 21 Jan 2015 Vortex solutions of the generalized Beltrami flows to the incompressible Euler equations Minoru Fujimoto 1, Kunihiko Uehara 2 and Shinichiro

More information

Two-dimensional model problem to explain counter-rotating vortex pair formation in a transverse jet

Two-dimensional model problem to explain counter-rotating vortex pair formation in a transverse jet PHYSICS OF FLUIDS 18, 085103 2006 Two-dimensional model problem to explain counter-rotating vortex pair formation in a transverse jet Suman Muppidi and Krishnan Mahesh Aerospace Engineering & Mechanics,

More information

Fine Grid Numerical Solutions of Triangular Cavity Flow

Fine Grid Numerical Solutions of Triangular Cavity Flow Published in : The European Physical Journal - Applied Physics (2007) Eur. Phys. J. Appl. Phys. 2007; Vol 38: pp 97-105 Fine Grid Numerical Solutions of Triangular Cavity Flow Ercan Erturk 1 and Orhan

More information

Fluid Mechanics. Chapter 9 Surface Resistance. Dr. Amer Khalil Ababneh

Fluid Mechanics. Chapter 9 Surface Resistance. Dr. Amer Khalil Ababneh Fluid Mechanics Chapter 9 Surface Resistance Dr. Amer Khalil Ababneh Wind tunnel used for testing flow over models. Introduction Resistances exerted by surfaces are a result of viscous stresses which create

More information

Liquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses

Liquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses Liquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses by W. A. Sirignano Mechanical and Aerospace Engineering University of California, Irvine Collaborators:

More information

Fluid Mechanics II Viscosity and shear stresses

Fluid Mechanics II Viscosity and shear stresses Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small

More information

J. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and

J. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and J. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and turbulent, was discovered by Osborne Reynolds (184 191) in 1883

More information

Intermittency of quasi-static magnetohydrodynamic turbulence: A wavelet viewpoint

Intermittency of quasi-static magnetohydrodynamic turbulence: A wavelet viewpoint Intermittency of quasi-static magnetohydrodynamic turbulence: A wavelet viewpoint Naoya Okamoto 1, Katsunori Yoshimatsu 2, Kai Schneider 3 and Marie Farge 4 1 Center for Computational Science, Nagoya University,

More information

Direct numerical simulation of turbulent Taylor Couette flow

Direct numerical simulation of turbulent Taylor Couette flow J. Fluid Mech. (7), vol. 587, pp. 7 9. c 7 Cambridge University Press doi:.7/s7767 Printed in the United Kingdom 7 Direct numerical simulation of turbulent Taylor Couette flow S. DONG Center for Computational

More information

ROLE OF THE VERTICAL PRESSURE GRADIENT IN WAVE BOUNDARY LAYERS

ROLE OF THE VERTICAL PRESSURE GRADIENT IN WAVE BOUNDARY LAYERS ROLE OF THE VERTICAL PRESSURE GRADIENT IN WAVE BOUNDARY LAYERS Karsten Lindegård Jensen 1, B. Mutlu Sumer 1, Giovanna Vittori 2 and Paolo Blondeaux 2 The pressure field in an oscillatory boundary layer

More information

NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD

NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias

More information

doi: / (

doi: / ( doi: 10.1063/1.1921949(http://dx.doi.org/10.1063/1.1921949) PHYSICS OF FLUIDS 17, 061701 2005 Experimental study of the production of vortex rings using a variable diameter orifice J. J. Allen a Department

More information

Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015

Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 I. Introduction (Chapters 1 and 2) A. What is Fluid Mechanics? 1. What is a fluid? 2. What is mechanics? B. Classification of Fluid Flows 1. Viscous

More information

Contents. Microfluidics - Jens Ducrée Physics: Laminar and Turbulent Flow 1

Contents. Microfluidics - Jens Ducrée Physics: Laminar and Turbulent Flow 1 Contents 1. Introduction 2. Fluids 3. Physics of Microfluidic Systems 4. Microfabrication Technologies 5. Flow Control 6. Micropumps 7. Sensors 8. Ink-Jet Technology 9. Liquid Handling 10.Microarrays 11.Microreactors

More information

Received: 5 September 2017; Accepted: 2 October 2017; Published: 11 October 2017

Received: 5 September 2017; Accepted: 2 October 2017; Published: 11 October 2017 fluids Article A Finite Element Method for Incompressible Fluid Flow in a Noninertial Frame of Reference and a Numerical Study of Circular Couette Flow with Varying Angular Speed C. S. Jog * and Nilesh

More information

AN OVERVIEW OF IMPELLERS, VELOCITY PROFILES AND REACTOR DESIGN

AN OVERVIEW OF IMPELLERS, VELOCITY PROFILES AND REACTOR DESIGN AN OVERVIEW OF IMPELLERS, VELOCITY PROFILES AND REACTOR DESIGN Instructor: Gurmeet Singh (Department Manager of Polymer and Processing Department IOCL R&D center, Faridabad) 1 Presented by: Pranay. K.

More information

Turbulent Rotating Rayleigh-Bénard Convection: DNS and SPIV Measurements

Turbulent Rotating Rayleigh-Bénard Convection: DNS and SPIV Measurements Turbulent Rotating Rayleigh-Bénard Convection: DNS and SPIV Measurements Rudie Kunnen 1 Herman Clercx 1,2 Bernard Geurts 1,2 1 Fluid Dynamics Laboratory, Department of Physics Eindhoven University of Technology

More information

LEAST-SQUARES FINITE ELEMENT MODELS

LEAST-SQUARES FINITE ELEMENT MODELS LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows

More information

M E 320 Supplementary Material Pralav Shetty

M E 320 Supplementary Material Pralav Shetty M E 320 Supplementary Material Pralav Shetty Note: In order to view the demonstrations below, you must first download CDF player to your PC/Mac/Linux. Link for CDF player http://www.wolfram.com/cdf-player/

More information

Flow and heat transfer characteristics of tornado-like vortex flow

Flow and heat transfer characteristics of tornado-like vortex flow Advances in Fluid Mechanics VI 277 Flow and heat transfer characteristics of tornado-like vortex flow Y. Suzuki & T. Inoue Department of Mechanical Sciences and Engineering, Tokyo Institute of Technology

More information

Hajime NAKAMURA and Tamotsu IGARASHI

Hajime NAKAMURA and Tamotsu IGARASHI 622 Heat Transfer in Separated Flow behind a Circular Cylinder for Reynolds Numbers from 120 to 30 000 (2nd Report, Unsteady and Three-Dimensional Characteristics) Hajime NAKAMURA and Tamotsu IGARASHI

More information

Three-dimensional Floquet stability analysis of the wake in cylinder arrays

Three-dimensional Floquet stability analysis of the wake in cylinder arrays J. Fluid Mech. (7), vol. 59, pp. 79 88. c 7 Cambridge University Press doi:.7/s78798 Printed in the United Kingdom 79 Three-dimensional Floquet stability analysis of the wake in cylinder arrays N. K.-R.

More information

The behaviour of high Reynolds flows in a driven cavity

The behaviour of high Reynolds flows in a driven cavity The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,

More information

Intensely swirling turbulent pipe flow downstream of an orifice: the influence of an outlet contraction

Intensely swirling turbulent pipe flow downstream of an orifice: the influence of an outlet contraction 13 th Int. Symp. on Appl. Laser Techniques to Fluid Mechanics, Lisbon, Portugal, June 26-29, 26 Intensely swirling turbulent pipe flow downstream of an orifice: the influence of an outlet contraction Marcel

More information

THE EFFECTS OF NONLINEARITY ON CROSSFLOW RECEPTIVITY

THE EFFECTS OF NONLINEARITY ON CROSSFLOW RECEPTIVITY THE EFFECTS OF NONLINEARITY ON CROSSFLOW RECEPTIVITY Christian Thomas, Philip Hall, Christopher Davies Imperial College London, Cardiff University Keywords: Boundary-layer, crossflow, nonlinearity, disturbance

More information

Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )

Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University

More information

Modeling of turbulence in stirred vessels using large eddy simulation

Modeling of turbulence in stirred vessels using large eddy simulation Modeling of turbulence in stirred vessels using large eddy simulation André Bakker (presenter), Kumar Dhanasekharan, Ahmad Haidari, and Sung-Eun Kim Fluent Inc. Presented at CHISA 2002 August 25-29, Prague,

More information

Lecture 1. Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev

Lecture 1. Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev Lecture Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev Introduction In many cases in nature, like in the Earth s atmosphere, in the interior of stars and planets, one sees the appearance

More information

What is Turbulence? Fabian Waleffe. Depts of Mathematics and Engineering Physics University of Wisconsin, Madison

What is Turbulence? Fabian Waleffe. Depts of Mathematics and Engineering Physics University of Wisconsin, Madison What is Turbulence? Fabian Waleffe Depts of Mathematics and Engineering Physics University of Wisconsin, Madison it s all around,... and inside us! Leonardo da Vinci (c. 1500) River flow, pipe flow, flow

More information

Interfacial waves in steady and oscillatory, two-layer Couette flows

Interfacial waves in steady and oscillatory, two-layer Couette flows Interfacial waves in steady and oscillatory, two-layer Couette flows M. J. McCready Department of Chemical Engineering University of Notre Dame Notre Dame, IN 46556 Page 1 Acknowledgments Students: M.

More information