NUMERICAL SIMULATION OF FLUID FLOW IN ENCLOSED ROTATING FILTER AND DISK

Similar documents
A new class of exact solutions of the Navier Stokes equations for swirling flows in porous and rotating pipes

NATURAL CONVECTION HEAT TRANSFER WITHIN VERTICALLY ECCENTRIC DOMED SKYLIGHTS CAVITIES

Effect of Downstream Flapping Plate on the Flow Field Characteristics behind a Circular Cylinder

Axisymmetric Stokes Flow past a Swarm of Porous Cylindrical Shells

SCRAPED SURFACE HEAT EXCHANGER (SSHE)

An Exact Solution of Navier Stokes Equation

UNSTEADY MIXED CONVECTION WITH HEATING EFFECTS IN A HORIZONTAL ANNULUS WITH CONSTANT ROTATING OUTER CYLINDER AND CONSTANT MAGNETIC FIELD

Stress, Cauchy s equation and the Navier-Stokes equations

EFFECTS OF SWIRL IMPARTED AT INLET ON THE RECIRCULATION BUBBLE GENERATED IN THE TURBULENT FLOW IN AN AXI-SYMMETRIC SUDDEN EXPANSION

Do not turn over until you are told to do so by the Invigilator.

In the previous section we considered problems where the

Spatial Development of Viscous Flow Induced by Wave Propagation over Vortex Ripples

PROBLEM SET #3A. A = Ω 2r 2 2 Ω 1r 2 1 r2 2 r2 1

Exercise sheet 8 (Modeling large- and small-scale flows) 8.1 Volcanic ash from the Eyjafjallajökull

Institute of Fundamental Sciences-Mathematics,Massey University, Palmerston North,NZ

7.2.1 Basic relations for Torsion of Circular Members

Analysis of vortex breakdown in an enclosed cylinder based on the energy gradient theory

Analysis of high speed machining center spindle dynamic unit structure performance Yuan guowei

The Combined Effect of Chemical reaction, Radiation, MHD on Mixed Convection Heat and Mass Transfer Along a Vertical Moving Surface

CBE Transport Phenomena I Final Exam. December 19, 2013

PARTICLE TRANSPORT AND DEPOSITION IN POROUS STRUCTURES: EFFECTS OF PARTICLE PROPERTIES, POROSITY AND REYNOLDS NUMBER

Chapter Introduction to Finite Element Methods

INFLUENCE OF GROUND INHOMOGENEITY ON WIND INDUCED GROUND VIBRATIONS. Abstract

A Most Useful Device of Studying Electrode Processes: The Rotating Disk Electrode

Application of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates

, and the curve BC is symmetrical. Find also the horizontal force in x-direction on one side of the body. h C

Diffusion and Transport. 10. Friction and the Langevin Equation. Langevin Equation. f d. f ext. f () t f () t. Then Newton s second law is ma f f f t.

12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,

Liquid gas interface under hydrostatic pressure

Radial Inflow Experiment:GFD III

Hydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods

THERMODYNAMIC OPTIMIZATION OF TUBULAR HEAT EXCHANGERS BASED ON MINIMUM IRREVERSIBILITY CRITERIA

Design, Simulation and Experimental of the Very Low Head Turbine with Minimum Pressure and Free Vortex Criterions

Instability of Taylor-Couette Flow between Concentric Rotating Cylinders

Effect of drag on the performance for an efficient wind turbine blade design

Chapter 2 Modelling of Fluid Flow and Heat Transfer in Rotating-Disk Systems

ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS

B. Spherical Wave Propagation

Research Article EXPERIMENTAL STUDY OF HEAT TRANSFER CHARACTERISTICS OF STAINLESS STEEL FIBROUS FLOW INSULATOR

International Journal of Mathematical Archive-9(5), 2018, Available online through ISSN

arxiv: v1 [physics.flu-dyn] 21 Dec 2018

Analysis and Optimization of a Special Type of Dielectric Loaded Resonant Cavity for Mobile Communication Filters

STUDY ON 2-D SHOCK WAVE PRESSURE MODEL IN MICRO SCALE LASER SHOCK PEENING

PROBLEM SET 5. SOLUTIONS March 16, 2004

Applied Aerodynamics

Incompressible Viscous Flows

Direct numerical simulation of Open Von Kármán Swirling Flow *

PRESSURE DRAWDOWN EQUATIONS FOR MULTIPLE-WELL SYSTEMS IN CIRCULAR-CYLINDRICAL RESERVOIRS

Natural Convection Heat Transfer in Horizontal Concentric Annulus between Outer Cylinder and Inner Flat Tube

A 1. EN2210: Continuum Mechanics. Homework 7: Fluid Mechanics Solutions

Right-handed screw dislocation in an isotropic solid

Construction Figure 10.1: Jaw clutches

Modeling of the fermentation in an internal loop airlift reactor

ASTR415: Problem Set #6

J. Electrical Systems 1-3 (2005): Regular paper

is the instantaneous position vector of any grid point or fluid

1) Consider an object of a parabolic shape with rotational symmetry z

A scaling-up methodology for co-rotating twin-screw extruders

T. Raja Rani. Military Technological College, Muscat, Oman. Abstract

LINEAR AND NONLINEAR ANALYSES OF A WIND-TUNNEL BALANCE

CHAPTER 25 ELECTRIC POTENTIAL

On the Sun s Electric-Field

Averaging method for nonlinear laminar Ekman layers

Numerical solution of diffusion mass transfer model in adsorption systems. Prof. Nina Paula Gonçalves Salau, D.Sc.

Numerical simulation of velocity and temperature fields in natural circulation loop

Raju Ananth Structural Integrity Associates, Inc Hellyer Avenue, Suite 210 San Jose, CA

Asymptotic Model of the Flow in Metallic and Elastomeric Progressive Cavity Pump

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere

Identification of the Hardening Curve Using a Finite Element Simulation of the Bulge Test

DEMONSTRATION OF INADEQUACY OF FFOWCS WILLIAMS AND HAWKINGS EQUATION OF AEROACOUSTICS BY THOUGHT EXPERIMENTS. Alex Zinoviev 1

MECHANICAL PULPING REFINER MECHANICAL PULPS

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s

The influence of a micropolar fluid on peristaltic transport in an annulus: application of the clot model

A Three-Dimensional Magnetic Force Solution Between Axially-Polarized Permanent-Magnet Cylinders for Different Magnetic Arrangements

DESIGN OF INTERMEDIATE RING STIFFENERS FOR COLUMN-SUPPORTED CYLINDRICAL STEEL SHELLS

General Railgun Function

Performance Prediction of a Cavitating Rim Driven Tunnel Thruster

Stress Intensity Factor

Two Dimensional Inertial Flow of a Viscous Fluid in a Corner

Three dimensional flow analysis in Axial Flow Compressors

Supplementary Figure 1. Circular parallel lamellae grain size as a function of annealing time at 250 C. Error bars represent the 2σ uncertainty in

1 Similarity Analysis

1 Dark Cloud Hanging over Twentieth Century Physics

Magnetic Field and Inductance Calculations in Theta-Pinch and Z-Pinch Geometries

MODELING OF CYLINDRICAL COUETTE FLOW OF RAREFIED GAS. THE CASE OF ROTATING INNER CYLINDER *

PES 3950/PHYS 6950: Homework Assignment 6

Review: Electrostatics and Magnetostatics

Fluid flow in curved geometries: Mathematical Modeling and Applications

2 Governing Equations

Designing a Sine-Coil for Measurement of Plasma Displacements in IR-T1 Tokamak

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

Computational Methods of Solid Mechanics. Project report

BIFURCATION ANALYSIS FOR A QUASI-LINEAR CENTRIFUGAL PENDULUM ABSORBER

Chapter 7-8 Rotational Motion

Dynamics of Rotating Discs

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM

Duality between Statical and Kinematical Engineering Systems

A dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids

Chapter 13 Gravitation

ELECTRONICS COOLING FAN NOISE PREDICTION

Transcription:

VOL., NO. 9, SEPTEMBER 010 ISSN 1819-6608 ARPN Jounal of Engineeing and Applied Sciences 006-010 Asian Reseach Publishing Netwok (ARPN). All ights eseved. www.apnjounals.com NUMERICAL SIMULATION OF FLUID FLOW IN ENCLOSED ROTATING FILTER AND DISK Yaan Taamneh Depatment of Mechanical Engineeing, Tafila Technical Univesity Tafila, Jodan E-Mail: ytaamneh@yahoo.com ABSTRACT Computational fluid dynamics is used to investigate the behaviou of the fluid flow nea otating filte and disk suface in closed housing. The study focused on the lamina flow and wide gap egion. A commecial finite volume package FLUENT was used to visualie the natue of the flow suounding the otating disk. The fluid is assumed to be Newtonian, incompessible, non-fouling and isothemal. It was found that the viscous dag coefficient on the otating filte and disk suface is stongly govened by the Reynolds numbe. The otating filte in a closed housing was found to poduce high viscous dag (shea stess) compaed to otating disk. Inceasing the otation speed of the filte incease the static pessue inside the filte due to centifugal foce induced inside it. Theefoe, a back pessue phenomenon may be aised. Viscous dag coefficient fo the case of disk otating in closed housing was found to be in good ageement with pevious expeimental esults. Keywods: CFD simulation, otating filte, lamina flow, wide gap. 1. INTRODUCTION The estimation of the toque due to skin fiction on otating disks is impotant because of its application in the calculation of losses in many diffeent industial machiney [1-4]. The pioneeing wok on the flow due to otating disk has been done by von Kaman [-6]. He was pehaps the fist who deive a theoetical expession fo the toque on a smooth disk otating in an infinite fluid when the flow is tubulent. Ealy chaacteiation of the hydodynamics of a otating disk in a closed cylindical enclosue was epoted by Daily and Nece in 1960 [7-10]. They showed that the hydodynamic of a otating disk with adius and axial gap s in housing ae govened by two paametes. Gap Reynolds numbe whee the width of the gap is the chaacteistic length ( Re s = ω s / υ ) o the adial Reynolds numbe whee the adius of the disk is the chaacteistic ( Re = ω / υ ), the second paamete is the atio between the width of the gap and the disk adius ( s/ ). By measuing the toque acting on the shaft and coelating the fictional toque coefficient with the Reynolds numbe and axial cleaance, they wee able to diffeentiate fou flow egimes. Region Ι: Lamina flow and naow gap ( Re 4 Re.10 ), the bounday s, layes ae meged and poduce a shea ate vaying invesely with spacing s. Region ΙΙ: lamina flow and sepaate bounday layes ( Res 4, Re.10 ), sepaate lamina bounday layes ae fomed and the shea ate is independent on s. between the two layes will develop a one that moves like a solid body with a otation speed of Kω whee 0 < K >1. Region ΙΙΙ: Tubulent flow and naow gap ( Re.10 and s/ 0.0), the flow is chaacteised by meged tubulent bounday layes. Region ΙV: Tubulent flow and wide gap ( Re.10 and s/ 0.0), two tubulent bounday layes with the angula velocity Kω, whee 0 < K > 1. Thee seems to be lack of computational woks on estimating the skin fiction fo a poous disk otating in a closed cylinde flow in axisymmetic lamina flow egime. Theefoe, this pape aims to povide a CFD simulation study that focused on lamina flow egime and wide gap fo a smooth poous disk otating in a closed cylinde ( Re.10 and Res 4). Anothe sub goal of the pesent study is to test whethe FLUENT, a commecial Computational Fluid Dynamics softwae package, is capable of poviding the solutions fo the poblem unde consideation.. THEORY.1. Govening equation The steady lamina flow of viscous incompessible fluid and axisymmatic nea an enclosed otating disk with a constant angula velocity is consideed. Unde these conditions the incompessible Navie- Stokes equations in cylindical pola coodinates educe to. v v v 1 p v v v v θ + = + + + v ν (1) ρ v θ vθ v vθ vθ vθ vθ v + = + v ν + + () v v v + v v + v v + 1 p v 1 v v = + ν + + ρ = 0 whee v,, v ae the adial, tangential and axial v θ components of the velocity, ρ is the fluid density and ν the (3) (4) 48

VOL., NO. 9, SEPTEMBER 010 ISSN 1819-6608 ARPN Jounal of Engineeing and Applied Sciences 006-010 Asian Reseach Publishing Netwok (ARPN). All ights eseved. www.apnjounals.com dynamic viscosity of the fluid. P denotes the fluid pessue... Bounday conditions Since the geomety of an axisymmetic geomety is consideed a half geomety pat otated about an axis paallel to the velocity inlet is modeled. The bottom bounday of the domain is modeled as an axis bounday. The suface of the housing is modeled as a wall. The voids of the otating filte disk is modeled as poous jump, whee as the solid pat as a wall. The exit of the pemeable filte disk is modeled as pessue outlet. The flow inlet to the closed housing is modeled as velocity inlet. Figue-1. Schematic of the enclosed otating disk filte showing the solution domain and Computational gid with bounday conditions... Computational fluid dynamics A finite volume method is employed using commecial softwae FLUENT 6. to solve the govening equations subject to specified bounday conditions. The model has been solved in two dimensions and the entie egion is meshed (see Figue 1). The CFD tool can also accommodate numeous epesentation of flow in poous media. Fo a Newtonian isothemal and incompessible fluid, flow though a poous medium can be descibed using the poous jump model [8]. No slip bounday conditions ae set at all solid wall boundaies. The flow inside the poous media is assumed to be lamina. Moe cells ae constucted nea the suface of the disk to compensate the high velocity gadient in the bounday laye egion of the viscous flow. A commecial softwae GAMBIT is used fo gid geneation. The coupling between the pessue and velocity fields is achieved using PISO. A second ode upwind scheme is used fo the convection. Hee in this study, the dimensionless toque coefficient (viscous dag coefficient) fo a disk wetted on both sides and a Reynolds numbe can be defined as: C M M = 1 ρω o M = o τ θ 0 ( πd) Re ω o = () ν Whee 0 is the oute adius of the disk, M is the total toque equied to tun a disk and, τ the tangential wall θ shea stess on the disk. The gid independence is achieved by compaing the esults of the diffeent gid cell sie. It was found that 4000 cells ae satisfactoy, and any incease beyond this sie would lead to an insignificant change in the esulting solution. 3. RESULTS AND DISCUSSIONS CFD Simulation esults fo smooth otating disk in closed housing ae fistly compaed to expeimental data to veify the validity of the CFD simulation solution. Figue shows the total viscous dag coefficient as a function of Reynolds numbe fo the case of otating disk. As can be seen fom Figue, thee is an excellent ageement in the Reynolds numbe dependence of C M between CFD simulations in this study and the expeimental measued dependence by []. The effects of Reynolds numbe on the total viscous dag coefficient fo otating filte and disk ae shown in Figue 3. It is clea that C M values gadually decease with incease in Reynolds numbe fo otating filte and disk. It can also be seen that the otating filte exhibit lage viscous dag coefficient compaed to otating disk due to the suction 49

VOL., NO. 9, SEPTEMBER 010 ISSN 1819-6608 ARPN Jounal of Engineeing and Applied Sciences 006-010 Asian Reseach Publishing Netwok (ARPN). All ights eseved. www.apnjounals.com effect on the bounday layes. This means, that moe enegy input is needed to tun the filte. The diffeence in the value of the viscous dag between the otating filte and disk decease significantly by inceasing the Reynolds numbe especially at high Reynolds numbe. The eason behind that efes to the fluid flow tubulences. Figues 4 and Figue show the axial distibution of swil and adial velocities component espectively fo otating filte and disk. It is clea that fo otating filte the swil velocity component decay shaply away fom the filte suface. This in tuns decease the bounday laye thickness nea the otating filte suface and enhance the shea ate. The axial distibution adial velocity components close to the otating filte suface inceasing due to the effect of suction (see Figue ). Figue-. Compaison of computed viscous dag of disk otating in a housing with the expeimental coelation Schlichting []. Figue-3. Vaiation of the viscous dag as a function of Reynolds numbe fo otating filte and disk. 0

VOL., NO. 9, SEPTEMBER 010 ISSN 1819-6608 ARPN Jounal of Engineeing and Applied Sciences 006-010 Asian Reseach Publishing Netwok (ARPN). All ights eseved. www.apnjounals.com Figue-4. The axial distibution of swil velocity fo otating filte and disk at a adial distance /R = 0. fo Re =.3 10 4 Figue-. The axial distibution of adial velocity fo otating filte and disk at a adial distance /R = 0. fo Re =.3 10 4 The effect of filte otation on the flow patten is illustated clealy on steamline contou Figue 6 (a). The otation speed of the filte disk leads to fluid ciculation. Figue 6 (b) shows the contous of static pessue fo otating filte in closed housing. It can be seen that the local pessue suounding the otating filte is almost constant. Wheeas, the static pessue inside the filte vaies accoding to the centifugal foce induced inside. Since the centifugal foce depends only on the otation speed of the filte and not on the local pessue inside the housing, we expect that by inceasing the otation speed of the filte a back pessue is poduced and thus low pemeate ate is obtained. 1

VOL., NO. 9, SEPTEMBER 010 ISSN 1819-6608 ARPN Jounal of Engineeing and Applied Sciences 006-010 Asian Reseach Publishing Netwok (ARPN). All ights eseved. www.apnjounals.com (a) (b) Figue-6. (a) Contous of static pessue in Pa (b) vaiation of contous of steam function fo otating disk at Re =.3 10 4 4. CONCLUSIONS The steady lamina flow of viscous incompessible fluid induced by a smooth otating filte and disk in closed housing using computational fluid dynamics ae caied out. The natue of the fluid flow nea a otating disk and filte inside a closed housing was visualised. The dependency of the viscous dag coefficient on the Reynolds numbe was shown. It was found that the otating filte exhibit lage viscous dag coefficient due to the suction effect on the bounday layes. A high otation speed of filte may poduced a back pessue which could be educed by inceasing the inlet feed pessue. Compaison the simulation esults with the expeimental data validate the commecially-available softwae FLUENT in poviding a easonable good solution of complicated flow stuctues, including flow with ciculation. It emains to be veified the effect of filte pemeability, feed pessue, angula velocity of the filte on the otating filte pefomance. Such simulation ae cuently unde way and esults will be epoted in fothcoming manuscipt REFERENCES [1] Daily JW, Nece RE. 1960. Chambe dimension effects on induced flow and fictional esistance of enclosed otating disks. J. Basic Eng. 8: 648-63. [] J. Mukes and C.G. Calsson. 1988. Cossflow Filtation. Wiley, New Yok. [3] Y. Taamneh. 009. Numeical and Expeimental Investigation of Dynamic Filtes. PhD Thesis. TU Kaiseslauten, Gemany. [4] R. Booeae, M. Y. Jaffin, L. Dingand P. Paullie. 000. Influence of geomety and angula velocity on pefomance of a otating disk filte. AICHE j. 46: 7-6. [] Schlichting H. Bounday laye Theoy, New Yok: McGaw-Hill; 1969. [6] J.A. Engle. 1997. Investigation of membane filtation in a otating disc geomety: using of computational fluid dynamics and laboatoy evaluation. PhD thesis. Rice Univesity, Houston. [7] C. A. Sea, M. R. Wiesne. 000. A compaison of otating and stationay membane disk filtes using computational fluid dynamics. J. Memb. Sci. 16: 19-9. [8] Y. Taamneh, S. Rippege. 010. Appoximate Solution Using Momentum Integal in Estimating the Shea Stess ove a Poous Rotating Disk. Intenational Jounal of Dynamic of Fluids. 6() : 1-. [9] Y. Taamneh, S. Rippege. 008. Pefomance of single and double shaft disk sepaato. Physical Sepaation in Science and Engineeing. 008 (1): 1-

VOL., NO. 9, SEPTEMBER 010 ISSN 1819-6608 ARPN Jounal of Engineeing and Applied Sciences 006-010 Asian Reseach Publishing Netwok (ARPN). All ights eseved. [10] Y. Taamneh, L. Steinke, S. Rippege. 008. Investigation of Dynamic Filtes Using CFD. 10 th Wold Filtation Congess Poceedings. : 47-476. [11] Y. Taamneh, L. Steinke and S. Rippege. 007. Einfluss de Geometievehältnisse auf die Sepaation mittels otieende Membanen. Chemie Ingenieu Technik. 79(9): 147-147. www.apnjounals.com 3