Shear stress and shear rates for µ-slides y-shaped based on Computational Fluid Dynamics (CFD)

Similar documents
2, where dp is the constant, R is the radius of

Laminar Flow. Chapter ZERO PRESSURE GRADIENT

Basic Concepts: Drag. Education Community

CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE

STRUCTURAL ANALYSIS OF A WESTFALL 2800 MIXER, BETA = 0.8 GFS R1. By Kimbal A. Hall, PE. Submitted to: WESTFALL MANUFACTURING COMPANY

Multi-Fidelity Computational Flow Assurance for Design and Development of Subsea Systems and Equipment Simon Lo

Some CFD simulations for the design of the FCC ventilation system. 9/28/2015 A. Rakai EN-CV-PJ 2

An-Najah National University Civil Engineering Department. Fluid Mechanics. Chapter 1. General Introduction

CFD and Thermal Stress Analysis of Helium-Cooled Divertor Concepts

CFD ANALYSIS OF PRESSURE DISTRIBUTION IN SLIDE CONICAL BEARING LUBRICATED WITH NON-NEWTONIAN OIL

Figure 3: Problem 7. (a) 0.9 m (b) 1.8 m (c) 2.7 m (d) 3.6 m

Lecture 2 Flow classifications and continuity

MATHEMATICAL MODELING OF A COMBINED HOT-WATER HEATING SYSTEM BY MEANS OF THE FINITE ELEMENT METHOD

Comparison of Flow and Sedimentation Pattern for three Designs of Storm Water Tanks by Numerical Modelling

Circular Bearing Performance Parameters with Isothermal and Thermo-Hydrodynamic Approach Using Computational Fluid Dynamics

GENERALIZED NEWTONIAN FLUIDS AS LUBRICANTS IN THE HYDRODYNAMIC CONICAL BEARINGS A CFD ANALYSIS

Lecture 5. Labs this week: Please review ME3281 Systems materials! Viscosity and pressure drop analysis Fluid Bulk modulus Fluid Inertance

MOMENTUM PRINCIPLE. Review: Last time, we derived the Reynolds Transport Theorem: Chapter 6. where B is any extensive property (proportional to mass),

Computational model for particle deposition in turbulent gas flows for CFD codes

Studies on flow through and around a porous permeable sphere: II. Heat Transfer

Entropic Evaluation of Dean Flow Micromixers

GENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION

Tutorial for the heated pipe with constant fluid properties in STAR-CCM+

Numerical study of blood fluid rheology in the abdominal aorta

16 Rainfall on a Slope

Differential relations for fluid flow

Supplementary Information for Engineering and Analysis of Surface Interactions in a Microfluidic Herringbone Micromixer

Polymerization Technology Laboratory Course

Practical Analysis Of Turbulent Flow In A Pipe Using Computational Fluid Dynamics

Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel

Simulation of Free Convection with Conjugate Heat Transfer

COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER

Chapter (4) Motion of Fluid Particles and Streams

MODA. Modelling data documenting one simulation. NewSOL energy storage tank

What s important: viscosity Poiseuille's law Stokes' law Demo: dissipation in flow through a tube

PHYSICAL MECHANISM OF CONVECTION

Analysis of Hydrodynamic Journal Bearing Using CFD and FSI Technique

MODELLING PARTICLE DEPOSITION ON GAS TURBINE BLADE SURFACES

Homework 6.1 P = 1000 N. δ δ δ. 4 cm 4 cm 4 cm. 10 cm

NUMERICAL AND EXPERIMENTAL INVESTIGATION OF THE TEMPERATURE DISTRIBUTION INSIDE OIL-COOLED TRANSFORMER WINDINGS

Tutorial for the supercritical pressure pipe with STAR-CCM+

Homework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods

Optimization of a Flocculation-Sedimentation Treatment Plant with Computational Fluid Dynamics

6. Basic basic equations I ( )

FLUID MECHANICS. Lecture 7 Exact solutions

Effect of Shape and Flow Control Devices on the Fluid Flow Characteristics in Three Different Industrial Six Strand Billet Caster Tundish

Applied Aerodynamics - I

Evaluating New Generation Vibrating Tube Sensor for Density Measurement under Process Conditions

report: Computational Fluid Dynamics Modelling of the Vortex Ventilator MK4 Rev 2 Ventrite International

Geomorphology 5. Stream Sediment Stream Sediment

Design and Modeling of Fluid Power Systems ME 597/ABE Lecture 7

Analysis of Fully Developed Turbulent Flow in a AXI-Symmetric Pipe using ANSYS FLUENT Software

Two-scale momentum theory for very large wind farms

CFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate

V/ t = 0 p/ t = 0 ρ/ t = 0. V/ s = 0 p/ s = 0 ρ/ s = 0

FE Exam Fluids Review October 23, Important Concepts

1.060 Engineering Mechanics II Spring Problem Set 2

ME 431A/538A/538B Homework 22 October 2018 Advanced Fluid Mechanics

Mapping of Deformation to Apparent Young s Modulus in Real-Time Deformability Cytometry. Christoph Herold

Numerical analysis of three-lobe journal bearing with CFD and FSI

Chapter 5. The Differential Forms of the Fundamental Laws

Supplementary Information. Continuous Transfer of Liquid Metal Droplets Across a Fluid-Fluid Interface Within an Integrated Microfluidic Chip

1. Introduction, fluid properties (1.1, 2.8, 4.1, and handouts)

HUMBERTO AVILA, M.Sc. Ph.D. student. Factors Affecting Scour of Previously Captured Sediment from Stormwater Catchbasin Sumps. Experimental Design

Direct numerical simulation database for supercritical carbon dioxide

Today s menu. Last lecture. A/D conversion. A/D conversion (cont d...) Sampling

Microfabricated Passive Vapor Preconcentrator/Injector Designed for Microscale Gas Chromatography

International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN

Homework #4 Solution. μ 1. μ 2

Simplified Model of WWER-440 Fuel Assembly for ThermoHydraulic Analysis

FE Fluids Review March 23, 2012 Steve Burian (Civil & Environmental Engineering)

MULTIPLE-CHOICE PROBLEMS:(Two marks per answer) (Circle the Letter Beside the Most Correct Answer in the Questions Below.)

meters, we can re-arrange this expression to give

Drag Coefficient of Tall Building by CFD Method using ANSYS

CFD Analysis of IV Smart - 50 Hz 2280 CMH Impulse Fan

ESO201A: Thermodynamics

Rheometer: Procedure: Part A: Viscosity v Time

CFD Modelling of Hydrodynamic Phenomena and Heat Transfer in Channels of Plate Heat Exchangers

Modeling of Anisotropic Polymers during Extrusion

Investigation of reattachment length for a turbulent flow over a backward facing step for different step angle

CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE

MAHALAKSHMI ENGINEERING COLLEGE TIRUCHIRAPALLI

Computational study on non-asymptotic behavior of developing turbulent pipe flow

Peristaltic Pump. Introduction. Model Definition

Research Article Analysis of Journal Bearing Performance in Two Dimensions

Hydraulic Fundamentals Hydraulics Definition Advent of Oil Hydraulics Pascal s Law Pressure

UNIT I FLUID PROPERTIES AND STATICS

A multiscale framework for lubrication analysis of bearings with textured surface

Boundary Conditions - Inlet

Fluid Flow Fluid Flow and Permeability

MODELING OF THE HYDROPLANING PHENOMENON. Ong G P and Fwa T F Department of Civil Engineering National University of Singapore INTRODUCTION

Natural Convection from a Long Horizontal Cylinder

Digital Simulation for the Behavior of the Flow of Non-Newtonian Fluids in 90 Pipe Bend

IHTC DRAFT MEASUREMENT OF LIQUID FILM THICKNESS IN MICRO TUBE ANNULAR FLOW

Air Drag on a Stratospheric Balloon in Tropical Regions. H.P. Lee, K.M. Lim, B.C. Khoo, Q. Sun

DEVELOPED LAMINAR FLOW IN PIPE USING COMPUTATIONAL FLUID DYNAMICS M.

Introduction to Micro/Nanofluidics. Date: 2015/03/13. Dr. Yi-Chung Tung. Outline

BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW

Applied CFD Project 1. Christopher Light MAE 598

Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition

Transcription:

Shear stress and shear rates for µ-slides y-shaped based on Computational Fluid Dynamics (CFD) 1) Introduction... 1 ) Computational Fluid Dynamics (CFD) Characteristics... 1 3) Using this document: calculating shear stress... 4) Regions of uniform laar shear stress... 3 a) Uniform laar region in one channel section; 100 % flow velocity... 3 b) Uniform laar region in double channel section; 50 % flow velocity... 4 5) Regions of non uniform laar shear stress... 4 6) Shear stress tables... 7 The shear stress is based the dynamical viscosity of water at C, η = 0.01 sec/² (=1 mpa sec = 1cP). Φ flow rate τ shear stress γ shear rate η dynamical viscosity (η = 0.01 sec/²) 1) Introduction The µ-slide y-shaped provides regions of uniform laar flow, and regions of non uniform laar flow. Non uniform laar flow is characterized by flow velocity gradients in small sub millimeter regions. Please keep in d that there is no turbulent flow in µ-slide y-shaped! ) Computational Fluid Dynamics (CFD) Characteristics The CFD calculation was done by ANSYS FLUENT Flow Modeling Software. The parameters are stationary, isothermal, laar, water assumed to be incompressible, density of water 1000 kg/m³, viscosity 0.001 kg/(m sec), and at given flow of mass = 1.001 /. Calculated is the velocity in 1 µm distance from the bottom. ibidi GmbH, Version 1.0, Roman Zantl, October 07 th 009 Page 1 of 8

3) Using this document: calculating shear stress It is assumed that the shear stress is proportional to the velocity calculated in 1 µm distance from the bottom. The value of 100 % velocity is defined as 0.0004 m/sec which is the velocity in the uniform flow region for the single channels. In this region the shear stress is.74 /² at a flow rate of 1/. For calculating the shear stress in another region one needs to take the velocity indicated by the color scale, for example 0.00018 m/sec for the yellow region and use the formula: τ [ = Φ velocity from color [m / sec /.74 velocity 100 % = 0.0004 m / sec 1/ calculation example: flow rate: / region: yellow region 0.00018 m/sec [ 0.00018 m / sec / τ [ =.74 = 3.4 / 0.0004 m / sec 1/ ibidi GmbH, Version 1.0, Roman Zantl, October 07 th 009 Page of 8

4) Regions of uniform laar shear stress Fig. 1: Bottom view of µ-slide y-shaped. The dashed rectangles are indicating zones of uniform laar flow and therefore also uniform laar shear stress. The flow and shear stress in the orange regions is exactly twice as high as in the green regions because the flow is exactly half in the regions between the branching points (yellow), compared to the regions between inlets (big blue spots on both ends) and the branching points. a) Uniform laar region in one channel section; 100 % flow velocity for viscosity of η = 0.01 sec/² Shear stress τ [ =.74 Φ[ Shear rate γ [sec 1 = 7.4 Φ[ for any viscosity η given in [ sec/² sec Shear stress τ [ = η [ 7.4 Φ[ ibidi GmbH, Version 1.0, Roman Zantl, October 07 th 009 Page 3 of 8

b) Uniform laar region in double channel section; 50 % flow velocity for viscosity of η = 0.01 sec/² Shear stress τ [ = 1.137 Φ[ [sec 1 Shear rate γ = 113.7 Φ[ for any viscosity η given in [ sec/² sec Shear stress τ [ = η [ 113.7 Φ[ 5) Regions of non uniform laar shear stress Fig. : Bottom view of µ-slide y-shaped. The dashed rectangles are indicating zones of non uniform laar flow and therefore also non uniform laar shear stress. There are three different regions: on the left side is a branching region with an opening angle of 45. On the right side the branching region provides a 45 opening angle. The both kink regions are symmetrical. ibidi GmbH, Version 1.0, Roman Zantl, October 07 th 009 Page 4 of 8

Fig. 3: Calculated data for the 30 branching region. The arrows are indicating the flow direction. The color scale gives the velocity in [m/sec. Fig. 4: Calculated data for the 45 branching region. The arrows are indicating the flow direction. The color scale gives the velocity in [m/sec. ibidi GmbH, Version 1.0, Roman Zantl, October 07 th 009 Page 5 of 8

Fig. 5: Calculated data for the 45 branching region. The arrows are indicating the flow direction. The color scale gives the velocity in [m/sec. ibidi GmbH, Version 1.0, Roman Zantl, October 07 th 009 Page 6 of 8

6) Shear stress tables Shear stress table - single channel region, uniform flow at η=0.01 s/²: τ [ =.744 Φ[ τ[/² Φ[/ τ[/² Φ[/ τ[/² Φ[/ 0,1 0,04 4 1,76 5 10,99 0, 0,09 5,0 30 13,19 0,3 0,13 6,64 35 15,39 0,4 0,18 7 3,08 40 17,59 0,5 0, 8 3,5 45 19,79 0,6 0,6 9 3,96 50 1,98 0,7 0,31 10 4,40 55 4,18 0,8 0,35 11 4,84 60 6,38 0,9 0,40 1 5,8 65 8,58 1,0 0,44 13 5,7 70 30,78 1, 0,53 14 6,16 75 3,98 1,4 0,6 15 6,60 80 35,17 1,6 0,70 16 7,03 85 37,37 1,8 0,79 17 7,47 90 39,57,0 0,88 18 7,91 95 41,77, 0,97 19 8,35 100 43,97,4 1,06 0 8,79 105 46,17,6 1,14 1 9,3 110 48,36,8 1,3 9,67 115 50,56 3,0 1,3 3 10,11 10 5,76 Shear rate table - single channel region: γ [sec = 7.44 Φ[ γ[1/sec Φ[/ γ[1/sec Φ[/ γ[1/sec Φ[/ 5 0,0 600,64 7000 30,78 10 0,04 700 3,08 7500 3,98 0 0,09 800 3,5 8000 35,17 30 0,13 900 3,96 8500 37,37 40 0,18 1000 4,40 9000 39,57 50 0, 150 5,50 9500 41,77 15 0,55 1500 6,60 10000 43,97 150 0,66 1750 7,69 10500 46,17 175 0,77 1800 7,91 11000 48,36 00 0,88 1850 8,13 11500 50,56 5 0,99 1900 8,35 1000 5,76 50 1,10 000 8,79 1500 54,96 75 1,1 500 10,99 13000 57,16 300 1,3 3000 13,19 13500 59,36 35 1,43 3500 15,39 14000 61,55 350 1,54 4000 17,59 14500 63,75 375 1,65 4500 19,79 15000 65,95 400 1,76 5000 1,98 15500 68,15 450 1,98 5500 4,18 16000 70,35 500,0 6000 6,38 16500 7,55 ibidi GmbH, Version 1.0, Roman Zantl, October 07 th 009 Page 7 of 8

Shear Stress table - double channel region uniform flow at η=0.01 s/²: τ [ = 1.137 Φ[ τ[/² Φ[/ τ[/² Φ[/ τ[/² Φ[/ 0,1 0,09 4 3,5 5 1,98 0, 0,18 5 4,40 30 6,38 0,3 0,6 6 5,8 35 30,78 0,4 0,35 7 6,16 40 35,17 0,5 0,44 8 7,03 45 39,57 0,6 0,53 9 7,91 50 43,97 0,7 0,6 10 8,79 55 48,36 0,8 0,70 11 9,67 60 5,76 0,9 0,79 1 10,55 65 57,16 1,0 0,88 13 11,43 70 61,55 1, 1,06 14 1,31 75 65,95 1,4 1,3 15 13,19 80 70,35 1,6 1,41 16 14,07 85 74,74 1,8 1,58 17 14,95 90 79,14,0 1,76 18 15,83 95 83,54, 1,93 19 16,71 100 87,94,4,11 0 17,59 105 9,33,6,9 1 18,47 110 96,73,8,46 19,35 115 101,13 3,0,64 3 0,3 10 105,5 Shear rate table - double channel region: γ [sec = 113.7 Φ[ γ[1/sec Φ[/ γ[1/sec Φ[/ γ[1/sec Φ[/ 5 0,04 600 5,8 7000 61,55 10 0,09 700 6,16 7500 65,95 0 0,18 800 7,03 8000 70,35 30 0,6 900 7,91 8500 74,74 40 0,35 1000 8,79 9000 79,14 50 0,44 150 10,99 9500 83,54 15 1,10 1500 13,19 10000 87,94 150 1,3 1750 15,39 10500 9,33 175 1,54 1800 15,83 11000 96,73 00 1,76 1850 16,7 11500 101,13 5 1,98 1900 16,71 1000 105,5 50,0 000 17,59 1500 109,9 75,4 500 1,98 13000 114,3 300,64 3000 6,38 13500 118,71 35,86 3500 30,78 14000 13,11 350 3,08 4000 35,17 14500 17,51 375 3,30 4500 39,57 15000 131,90 400 3,5 5000 43,97 15500 136,30 450 3,96 5500 48,36 16000 140,70 500 4,40 6000 5,76 16500 145,09 ibidi GmbH, Version 1.0, Roman Zantl, October 07 th 009 Page 8 of 8