ME 5302 Computational fluid mechanics
|
|
- Alban Porter
- 5 years ago
- Views:
Transcription
1 ME 5302 Computational fluid mechanics Assignment Part II Thursday, 23 April 2015 National University of Singapore ARAVIND BASKAR A
2 Question 1 Solution: Problem Statement:
3 Analysis: The conservative form of the navier stokes equation has been employed for the analysis. Pressure correction is used for updating the values of velocity and temperature. The algorithm for solving the problem is as shown below: Input parameters Boundary Conditions & Stability Criteria Temporary Velocity using NS equation - Conservative form Pressure Correction and error check Update Velocity Energy equation solution Plot Results Conservative form of Navier Stokes Equation is as follows: u t + u2 x + uv x = Pr u ( 2 x u P y2) x v t + v2 y + uv y = Pr u ( 2 x u P y2) + Pr. Ra. T x T t + ut x + vt y = Pr T ( 2 x T y 2)
4 The plots at various resolutions are as follows: Fig (a) Grid Size:30x30, Velocity profiles at Ra=10 3 Fig (b) Grid Size: 30x30, Velocity contour at Ra=10 3
5 X 10 3 X Fig (c) Grid Size: 30x30, Velocity contour at Ra=10 5
6 Fig (c) Grid Size: 30x30, Velocity contour at Ra=10 5 X 10-1
7 Fig (c) Grid Size: 30x30, Velocity contour at Ra=10 3 Fig (c) Grid Size: 30x30, Velocity contour at Ra=10 3
8 Fig (c) Grid Size: 20x20, Velocity contour at Ra=10 3 X 10-1
9 Fig (d) Grid Size: 20x20, Velocity Contour
10 Matlab Code % ME Computational Fluid Dynamics % Assignment - Cavity Flow with natural convection % Date - 23/4/2015 % Aravind Baskar - A J % Input parameters nx=10;ny=10; % Grid Points in x & y hx=1/nx;hy=1/ny; % Grid Size in x & y T_left=1; % Left Wall Temperature u_top=0; % Initial Velocity Ra=1e5; % Rayleigh Number Pr=0.7; % Prandtl Number tf=0.15; % Total time Beta=1.5; % Correction factor MaxIt=100; % Iterations MaxErr=0.001; % Error dt=1e-4; % Time step size Maxstep=tf/dt; % Time steps % Initialize vectors u=zeros(nx+1,ny+2); u_temp=zeros(nx+1,ny+2); v=zeros(nx+2,ny+1); v_temp=zeros(nx+2,ny+1); P=zeros(nx+2,ny+2); T=ones(nx+2,ny+2); t=0; %% Time Loop % Boundary conditions for n=1:maxstep %% Left Wall u(1,2:ny+1)=0; v(1,1:ny+1)=2*0-v(2,1:ny+1); T(1,2:ny+1)=T_left; %% Right Wall u(nx+1,2:ny+1)=0; v(nx+2,1:ny+1)=2*0-v(nx+1,1:ny+1); T(nx+1,2:ny+1)=T(nx+1,2:ny+1); %% Top Wall v(2:nx+1,1)=0;
11 u(1:nx+1,ny+2)=2*u_top-u(1:nx+1,ny+1); T(2:nx+1,ny+2)=T(2:nx+1,ny+1); %% Bottom Wall v(2:nx+1,ny+1)=0; u(1:nx+1,1)=2*0-u(1:nx+1,2); T(2:nx+1,1)=T(2:nx+1,2); % Temporary velocity and temperature for i=2:nx for j=2:ny+1 usqr_1=(0.5*(u(i+1,j)+u(i,j)))^2; usqr_2=(0.5*(u(i-1,j)+u(i,j)))^2; uv_x1=(0.5*(u(i,j)+u(i,j+1)))*(0.5*(v(i,j)+v(i+1,j))); uv_x2=(0.5*(u(i,j)+u(i,j-1)))*(0.5*(v(i,j-1)+v(i+1,j-1))); A=((usqr_1-usqr_2)/hx)+((uv_x1-uv_x2)/hy); D2_u=((u(i+1,j)-(2*u(i,j))+u(i-1,j))/hx^2)+((u(i,j+1)-(2*u(i,j))+u(i,j- 1))/hy^2); u_temp(i,j)=u(i,j)+(dt*((pr*d2_u)-a)); for i=2:nx+1 for j=2:ny vsqr_1=(0.5*(v(i,j+1)+v(i,j)))^2; vsqr_2=(0.5*(v(i,j)+v(i,j-1)))^2; uv_y1=(0.5*(u(i,j)+u(i,j+1)))*(0.5*(v(i,j)+v(i+1,j))); uv_y2=(0.5*(u(i-1,j+1)+u(i-1,j)))*(0.5*(v(i,j)+v(i-1,j))); B=((vsqr_1-vsqr_2)/hy)+((uv_y1-uv_y2)/hx); D2_v=((v(i+1,j)-(2*(v(i,j)))+v(i-1,j))/hx^2)+((v(i,j+1)-(2*(v(i,j)))+v(i,j- 1))/hy^2); v_temp(i,j)=v(i,j)+(dt*((pr*d2_v)+(pr*ra*t_temp)-b)); % Pressure Correction for it=1:maxit P_chk=P; for i=2:nx+1 for j=2:ny+1 if(i==2&&j==2) (i==nx+1&&j==2) (i==2&&j==ny+1) (i==nx+1&&j==ny+1) C=1; elseif(i>=3&&i<=nx&&j==2) (i==2&&j>=3&&j<=ny) (i==nx+1&&j>=3&&j<=ny ) (i>=3&&i<=nx&&j==ny+1) C=2/3;
12 else C=1/2; P(i,j)=((C*Beta/(hx^2+hy^2))*((hy^2*(P(i+1,j)+P(i-1,j)))+... (hx^2*(p(i,j+1)+p(i,j-1)))-(((hx*hy^2)/dt)*(u_temp(i,j)-... u_temp(i-1,j)))-(((hx^2*hy)/dt)*(v_temp(i,j)-v_temp(i,j-1)))))+ ((1-Beta)*(P(i,j))); Err=0; for i=2:nx+1 for j=2:ny+1 Err=Err+(abs(P_chk(i,j)-P(i,j-1))); ; ; if Err<=MaxErr break % Velocity Updation for i=2:nx for j=2:ny+1 u(i,j)=u_temp(i,j)-((dt/hx)*(p(i+1,j)-p(i,j))); for i=2:nx+1 for j=2:ny v(i,j)=v_temp(i,j)-((dt/hy)*(p(i,j+1)-p(i,j))); for i=2:nx+1 for j=2:ny+1 A_T=(0.5*(u(i-1,j)+u(i,j)))*((T(i+1,j)-T(i-1,j))/(2*hx))+((0.5*(v(i,j- 1)+v(i,j)))*((T(i,j+1)-T(i,j-1))/(2*hy))); D2_T=((T(i+1,j)-(2*T(i,j))+T(i-1,j))/hx^2)+(T(i,j+1)-(2*T(i,j))+T(i,j- 1)/hy^2); T(i,j)=T(i,j)+(dt*(D2_T-A_T)); t=t+dt;
13 % Plotting results P(1:nx+2,1)=0;P(1:nx+2,ny+2)=0; P(1,1:ny+2)=0; P(nx+2,1:ny+2)=0; u_cnt(1:nx+1,1:ny+1)=0.5*(u(1:nx+1,1:ny+1)+u(1:nx+1,2:ny+2)); v_cnt(1:nx+1,1:ny+1)=0.5*(v(1:nx+1,1:ny+1)+v(2:nx+2,1:ny+1)); P_cnt(1:nx+1,1:ny+1)=0.25*(P(1:nx+1,1:ny+1)+P(2:nx+2,1:ny+1)+P(1:nx+1,2:ny+2)+P(2:nx+2, 2:ny+2)); wt(1:nx+1,1:ny+1)=((v(2:nx+2,1:ny+1)-v(1:nx+1,1:ny+1))/hx)-((u(1:nx+1,2:ny+2)- u(1:nx+1,1:ny+1))/hy); T_cnt(1:nx+1,1:ny+1)=0.25*(T(1:nx+1,1:ny+1)+T(2:nx+2,1:ny+1)+T(1:nx+1,2:ny+2)+T(2:nx+2, 2:ny+2)); x(1:nx+1)=(0:nx); y(1:ny+1)=(0:ny); figure(1), quiver(x,y,(rot90(fliplr(u_cnt))),(rot90(fliplr(v_cnt)))),... xlabel('nx'),ylabel('ny'),title('velocity Vector Plot');axis('square','tight'); figure(2), subplot(211),contourf(x,y,rot90(fliplr(u_cnt)),30),... xlabel('nx'),ylabel('ny'),title('u-velocity');axis('square','tight');colorbar; subplot(212),contourf(x,y,rot90(fliplr(v_cnt)),30),... xlabel('nx'),ylabel('ny'),title('v-velocity');axis('square','tight');colorbar; figure(3),contourf(x,y,rot90(fliplr(p_cnt)),30),... xlabel('nx'),ylabel('ny'),title('pressure');axis('square','tight');colorbar; figure(4),contourf(x,y,rot90(fliplr(wt)),30),... xlabel('nx'),ylabel('ny'),title('vorticity');axis('square','tight');colorbar; figure(5),contourf(x,y,rot90(fliplr(t_cnt)),30),... xlabel('nx'),ylabel('ny'),title('temperature');axis('square','tight');colorbar; Inferences: The velocity of flow is increased as the Rayleigh number is increased from 10 3 to Convergence time reduces as the Rayleigh number is increased. Velocity contours are almost zero at the centre and the gradient increases as it nears the walls. Temperature can undergo an inversion at the centre of the cavity Temperature gradient is increased as it nears the walls as Ra increases Ra Numin Numax Nuavg Umax Vmax Numin Numax Nuavg Umax Vmax 30X30 20X
14 References: 1. Numerical simulation of natural convection in a square cavity by simple generalized differential quadrature method by C.Shu,K.H.A Wee, Computer & Fluids Journal (Pg ). 2. D.C. Wan, B.S.V. Patnaik and G.W. Wei, A new benchmark quality solution for the buoyancydriven cavity by discrete singular convolution. Taylor & Francis, Numerical Heat Transfer, Part B, Vol. 40, pages , D.A. Mayne, A.S. Usmani and M. Crapper, h-adaptive finite element solution of high Rayleigh number thermally driven cavity problem. International Journal of Numerical Methods for Heat & Fluid Flow,Vol. 10, No. 6, pages , 2000.
A Front-tracking/Finite-Volume Navier-Stokes Solver for Direct Numerical Simulations of Multiphase Flows. Grétar Tryggvason
A Front-tracking/Finite-Volume Navier-Stokes Solver for Direct Numerical Simulations of Multiphase Flows Grétar Tryggvason April 3, 20 2 Contents Introduction 5. Problem Specification........................
More informationNatural Convection in Parabolic Enclosure Heated from Below
www.ccsenet.org/mas Modern Applied Science Vol. 5, No. 3; June 011 Natural Convection in Parabolic Enclosure Heated from Below Dr. Ahmed W. Mustafa (Corresponding auther) University of Tikrit, College
More informationME469. Computational Methods in Fluid Mechanics. Handout #11 GI - W2015
ME469 Computational Methods in Fluid Mechanics Handout #11 GI - W2015 Thermally Driven Cavity The Problem δt/δn=0# g# T h# T c# δt/δn=0# Water in a heated cavity. What is the effect of gravity? What is
More informationProject 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions
Project 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions R. S. Sellers MAE 5440, Computational Fluid Dynamics Utah State University, Department of Mechanical and Aerospace Engineering
More informationTwo-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 informationLAB 5. INTRODUCTION Finite Difference Method
LAB 5 In previous two computer labs, you have seen how the analytical techniques for solving electrostatic problems can be approximately solved using numerical methods. For example, you have approximated
More informationHandout 8 MATLAB Code for the Projection Method and BC Details for the Lid-Driven Cavity Problem
Handout 8 MATLAB Code for the Projection Method and BC Details for the Lid-Driven Cavity Problem Let s simplify the equations of the last page of Handout 7 for the case h = Δx = Δy because the code given
More informationDOING PHYSICS WITH MATLAB ELECTRIC FIELD AND ELECTRIC POTENTIAL: POISSON S EQUATION
DOING PHYSICS WITH MATLAB ELECTRIC FIELD AND ELECTRIC POTENTIAL: POISSON S EQUATION Ian Cooper School of Physics, University of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS For
More informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
More informationA Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer
Copyright c 2008 ICCES ICCES, vol.6, no.1, pp.13-18 A Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer K agamani Devi 1,D.W.Pepper 2 Summary Over the past few years, efforts have
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
More informationA STUDY ON NATURAL CONVECTION HEAT TRANSFER IN COMPLEX BOUNDARIES
http://doi.org/10.4867/jpe-017-01-11 JPE (017) Vol.0 (1) Mohapatra, C. R. Preliminary Note A STUDY ON NATURAL CONVECTION HEAT TRANSFER IN COMPLEX BOUNDARIES Received: 3 February 017 / Accepted: 01 April
More informationA NEW BENCHM ARK QUALITY SOLUTION FOR THE BUOYANCY-DRIVEN CAVITY BY DISCRETE SINGULAR CONVOLUTION. D. C. Wan, B. S. V. Patnaik, and G. W.
Numerical Heat Transfer, Part B, 40: 199± 228, 2001 Copyright # 2001 Taylor & Francis 1040-7790 /01 $12.00 +.00 A NEW BENCHM ARK QUALITY SOLUTION FOR THE BUOYANCY-DRIVEN CAVITY BY DISCRETE SINGULAR CONVOLUTION
More informationEffect of an adiabatic fin on natural convection heat transfer in a triangular enclosure
American Journal of Applied Mathematics 2013; 1(4): 78-83 Published online November 10, 2013 (http://www.sciencepublishinggroup.com/j/ajam) doi: 10.11648/j.ajam.20130104.16 Effect of an adiabatic fin on
More informationNATURAL CONVECTION OF AIR IN TILTED SQUARE CAVITIES WITH DIFFERENTIALLY HEATED OPPOSITE WALLS
Proceedings of the International onference on Mechanical Engineering 0 (IME0 8-0 December 0, Dhaka, Bangladesh IME- NATURAL ONVETION OF AIR IN TILTED SQUARE AVITIES WIT DIFFERENTIALLY EATED OPPOSITE WALLS
More informationº Õߪí À À Ë ËÕπ â«ß Õ µ «â«µ Àπàß ÿ æ Ë Ëß øíß å π µ «Õ ªìπøíß å π π ß»
««æ π. ªï Ë 31 Ë 4 µÿ - π«2551 643 º Õߪí À À Ë ËÕπ â«ß Õ µ «â«µ Àπàß ÿ æ Ë Ëß øíß å π µ «Õ ªìπøíß å π π ß» µ π» «1 À µ å»ÿπ å ß µ ÕßÀ «ß ª ÿ π 12121 ËÕ 1 惻 π 2550 µõ ËÕ 20 ÿπ π 2551 àõ «µ Àπàß ÿ æ Ë
More informationAnalysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle
Nonlinear Analysis: Modelling and Control, 2011, Vol. 16, No. 1, 89 99 89 Analysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle S. Parvin,
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (24 ) 55 556 th International Conference on Mechanical Engineering, ICME 23 Analysis of heat transfer and flow due to natural
More informationVisualization of Natural Convection in Enclosure. Filled with Porous Medium by Sinusoidally. Temperature on the One Side
Applied Mathematical Sciences, Vol., 2012, no. 97, 801-812 Visualization of Natural Convection in Enclosure Filled with Porous Medium by Sinusoidally Temperature on the One Side Paweena Khansila Department
More informationNumerical computation of three-dimensional incompressible Navier Stokes equations in primitive variable form by DQ method
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2003; 43:345 368 (DOI: 10.1002/d.566) Numerical computation of three-dimensional incompressible Navier Stokes equations
More informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
More informationFREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER
Free Convective Heat Transfer From an Object at Low Rayleigh Number FREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER Md. Golam Kader and Khandkar Aftab Hossain * Department of Mechanical
More informationComputation Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computation Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. # 40 What is there in implementing a CFD Code So far, we
More informationNatural convection adjacent to a sidewall with three fins in a differentially heated cavity
ANZIAM J. 48 (CTAC2006) pp.c806 C819, 2007 C806 Natural convection adjacent to a sidewall with three fins in a differentially heated cavity F. Xu 1 J. C. Patterson 2 C. Lei 3 (Received 31 August 2006;
More informationConvection Workshop. Academic Resource Center
Convection Workshop Academic Resource Center Presentation Outline Understanding the concepts Correlations External Convection (Chapter 7) Internal Convection (Chapter 8) Free Convection (Chapter 9) Solving
More informationLattice Boltzmann Method for Fluid Simulations
1 / 16 Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas Simon Fraser University April 7, 2011 2 / 16 Ludwig Boltzmann and His Kinetic Theory of Gases The Boltzmann Transport
More informationAdditive Manufacturing Module 8
Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p
More informationSpatial discretization scheme for incompressible viscous flows
Spatial discretization scheme for incompressible viscous flows N. Kumar Supervisors: J.H.M. ten Thije Boonkkamp and B. Koren CASA-day 2015 1/29 Challenges in CFD Accuracy a primary concern with all CFD
More informationDNS of the Taylor-Green vortex at Re=1600
DNS of the Taylor-Green vortex at Re=1600 Koen Hillewaert, Cenaero Corentin Carton de Wiart, NASA Ames koen.hillewaert@cenaero.be, corentin.carton@cenaero.be Introduction This problem is aimed at testing
More informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
More informationTable 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 informationA Mixed Finite Element Formulation for Solving Phase Change Problems with Convection
A Mixed Finite Element Formulation for Solving Phase Change Problems with Convection Youssef Belhamadia 1, Abdoulaye S. Kane 2, and André Fortin 3 1 University of Alberta, Campus Saint-Jean and Department
More informationA numerical study of heat transfer and fluid flow over an in-line tube bank
Fluid Structure Interaction VI 295 A numerical study of heat transfer and fluid flow over an in-line tube bank Z. S. Abdel-Rehim Mechanical Engineering Department, National Research Center, Egypt Abstract
More informationTHERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW
THERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW Luigi De Giorgi, Carlo Cima, Emilio Cafaro Dipartimento di Energetica, Politecnico di Torino, Torino, Italy Volfango Bertola School
More informationNumerical Investigation into Natural Convection and Entropy Generation in a Nanofluid-Filled U-Shaped Cavity
Entropy 015, 17, 5980-5994; doi:10.3390/e17095980 Article OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Numerical Investigation into Natural Convection and Entropy Generation in a Nanofluid-Filled
More informationEffect of Hartmann Number on Free Convective Flow in a Square Cavity with Different Positions of Heated Square Block
Effect of Hartmann Number on Free Convective Flow in a Square Cavity with Different Positions of Heated Square Block Abdul Halim Bhuiyan, M. A. Alim, Md. Nasir Uddin Abstract This paper is concerned with
More informationSELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE
Computational Thermal Sciences, 3 (1): 63 72 (2011) SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE M. Zhao, 1, M. Yang, 1 M. Lu, 1 & Y. W. Zhang 2 1
More informationLecture 2: Fundamentals 2/3: Navier-Stokes Equation + Basics of Discretization
Lecture 2: Fundamentals 2/3: Navier-Stokes Equation Basics of Discretization V.Vuorinen Aalto University School of Engineering CFD Course, Spring 2018 January 15th 2018, Otaniemi ville.vuorinen@aalto.fi
More informationRAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8
HEFAT01 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 01 Malta RAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8 Leong S.S. School of Mechanical
More information10. Buoyancy-driven flow
10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation
More informationBenchmark solutions for the natural convective heat transfer problem in a square cavity
Benchmark solutions for the natural convective heat transfer problem in a square cavity J. Vierendeels', B.Merci' &L E. Dick' 'Department of Flow, Heat and Combustion Mechanics, Ghent University, Belgium.
More informationHeat transfer increase with thin fins in three dimensional enclosures
157 Heat transfer increase with thin fins in three dimensional enclosures R. L. Frederick & S. Samper Universidad de Chile, Departamento de Ingeniería Mecánica, Santiago, Chile Abstract Heat transfer enhancement
More informationMeysam ATASHAFROOZ, Seyyed Abdolreza GANDJALIKHAN NASSAB, and Amir Babak ANSARI
THERMAL SCIENCE: Year 014, Vol. 18, No., pp. 479-49 479 NUMERICAL INVESTIGATION OF ENTROPY GENERATION IN LAMINAR FORCED CONVECTION FLOW OVER INCLINED BACKWARD AND FORWARD FACING STEPS IN A DUCT UNDER BLEEDING
More informationNumerical Investigation of Fluid and Thermal Flow in a Differentially Heated Side Enclosure walls at Various Inclination Angles
Numerical Investigation of Fluid and Thermal Flow in a Differentially Heated Side Enclosure walls at Various Inclination Angles O. A. SHAHRUL Faculty of Mechanical & Manufacturing Engineering Universiti
More informationA Finite Element Analysis on MHD Free Convection Flow in Open Square Cavity Containing Heated Circular Cylinder
American Journal of Computational Mathematics, 2015, 5, 41-54 Published Online March 2015 in SciRes. http://www.scirp.org/journal/ajcm http://dx.doi.org/10.4236/ajcm.2015.51003 A Finite Element Analysis
More informationNUMERICAL SIMULATION OF THERMAL CONVECTION IN A CLOSED CAVITY IN THE PRESENCE OF A THIN HORIZONTAL HEATED PLATE
Proceedings of CHT-1 ICHMT International Symposium on Advances in Computational Heat Transfer July 1-6, 1, Bath, England CHT-1-NC17 NUMERICAL SIMULATION OF THERMAL CONVECTION IN A CLOSED CAVITY IN THE
More informationNatural Convection and Entropy Generation in a Porous Enclosure with Sinusoidal Temperature Variation on the Side Walls
Avestia Publishing Journal of Fluid Flow, Heat and Mass Transfer Volume 1, Year 14 Journal ISSN: 368-6111 DOI: 1.11159/jffhmt.14.4 Natural Convection and Entropy Generation in a Porous Enclosure with Sinusoidal
More informationMIXED CONVECTION IN A SQUARE CAVITY WITH A HEAT-CONDUCTING HORIZONTAL SQUARE CYLINDER
Suranaree J. Sci. Technol. Vol. 17 No. 2; April - June 2010 139 MIXED CONVECTION IN A SQUARE CAVITY WITH A HEAT-CONDUCTING HORIZONTAL SQUARE CYLINDER Md. Mustafizur Rahman 1 *, M. A. Alim 1 and Sumon Saha
More informationAbstract. Introduction
Combined forced and natural convection in a square cavity - numerical solution and scale analysis A.T. Franco/ M.M. Ganzarolli'' "DAMEC, CEFET, PR 80230-901, Curitiba, PR Brasil >>DE, FEM, UNICAMP 13081-970,
More informationDynamics in the Earth s core. Philippe Cardin, ISTerre, Université Grenoble Alpes et CNRS
Dynamics in the Earth s core Philippe Cardin, ISTerre, Université Grenoble Alpes et CNRS Doctoral training on internal Earth, Barcelonnette, oct 2016 Sources of motions inside the core Core cooling and
More informationDOING PHYSICS WITH MATLAB. ELECTRIC FIELD AND ELECTRIC POTENTIAL: POISSON S and LAPLACES S EQUATIONS
DOING PHYSICS WITH MATLAB ELECTRIC FIELD AND ELECTRIC POTENTIAL: POISSON S and LAPLACES S EQUATIONS Ian Cooper School of Physics, University of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECTORY FOR MATLAB
More informationModule 6: Free Convections Lecture 26: Evaluation of Nusselt Number. The Lecture Contains: Heat transfer coefficient. Objectives_template
The Lecture Contains: Heat transfer coefficient file:///d /Web%20Course%20(Ganesh%20Rana)/Dr.%20gautam%20biswas/Final/convective_heat_and_mass_transfer/lecture26/26_1.html[12/24/2014 6:08:23 PM] Heat transfer
More informationSimulation 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 informationSIMULATION OF MIXED CONVECTIVE HEAT TRANSFER USING LATTICE BOLTZMANN METHOD
International Journal of Automotive and Mechanical Engineering (IJAME) ISSN: 2229-8649 (Print); ISSN: 2180-1606 (Online); Volume 2, pp. 130-143, July-December 2010 Universiti Malaysia Pahang DOI: http://dx.doi.org/10.15282/ijame.2.2010.3.0011
More informationFluid 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 informationNumerical 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 informationNumerical solutions of 2-D incompressible driven cavity flow with wavy bottom surface
American Journal of Applied Mathematics 015; (1-1): 0-4 Published online December 5, 014 (http://www.sciencepublishinggroup.com/j/ajam) doi: 10.11648/j.ajam.s.01500101.14 ISSN: 0-004 (Print); ISSN: 0-006X
More informationAbstract: This paper numerically investigates the physical mechanism of flow instability
Hua-Shu Dou*, Gang Jiang, Numerical simulation of flow instability and heat transfer of natural convection in a differentially heated cavity, International Journal of Heat and Mass Transfer, 103 (2016)
More informationSIMPLE Algorithm for Two-Dimensional Flow - Fluid Flow and Heat Transfer
SIMPLE Algorithm for Two-Dimensional Flow - Fluid Flow and Heat Transfer by Professor Jung-Yang San Mechanical Engineering Department National Chung Hsing University Two-dimensional, transient, incompressible
More informationThe Effect of Undulated Wall on Natural Convection Heat Transfer in Vertical Channels
Proceedings of the 3 rd International Conference on Fluid Flow, Heat and Mass Transfer (FFHMT 16) Ottawa, Canada May 2 3, 2016 Paper No. 123 The Effect of Undulated Wall on Natural Convection Heat Transfer
More informationInfluence of rheological behavior of nanofluid on heat transfer
Influence of rheological behavior of nanofluid on heat transfer ADNAN RAJKOTWALA AND JYOTIRMAY BANERJEE * Department of Mechanical Engineering National Institute of Technology Surat (Gujarat) - 395007
More informationWALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID
WALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID Zia Ghiasi sghias@uic.edu Dongru Li dli@uic.edu Jonathan Komperda jonk@uic.edu Farzad
More informationThe investigation on SIMPLE and SIMPLER algorithm through lid driven cavity
29, Issue 1 (2017) 10-23 Journal of Advanced Research in Fluid Mechanics and Thermal Sciences Journal homepage: www.akademiabaru.com/arfmts.html ISSN: 2289-7879 The investigation on SIMPLE and SIMPLER
More informationNATURAL CONVECTION WITHIN TRAPEZOIDAL ENCLOSURE WITH TWO BAFFLES: EFFECT OF VARIOUS ANGLES OF INCLINATION
NATURAL CONVECTION WITHIN TRAPEZOIDAL ENCLOSURE WITH TWO BAFFLES: EFFECT OF VARIOUS ANGLES OF INCLINATION Éliton Fontana, eliton_fontana@hotmail.com Universidade Federal de Santa Catarina - UFSC Adriano
More informationLEAST-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 informationNumerical Study of Free Convection Heat Transfer in a Square Cavity with a Fin Attached to Its Cold Wall
Heat Transfer Research, 2011, Vol. 42, No. 3 Numerical Study of Free Convection Heat Transfer in a Square Cavity with a Fin Attached to Its Cold Wall SAEID JANI, 1* MEYSAM AMINI, 2 and MOSTAFA MAHMOODI
More informationCombined Natural Convection and Thermal Radiation in an Inclined Cubical Cavity with a Rectangular Pins Attached to Its Active Wall
Periodicals of Engineering and Natural Sciences ISSN 2303-4521 Vol.5, No.3, November 2017, pp. 347~354 Available online at:http://pen.ius.edu.ba Combined Natural Convection and Thermal Radiation in an
More informationOPTIMAL DESIGN OF CLUTCH PLATE BASED ON HEAT AND STRUCTURAL PARAMETERS USING CFD AND FEA
International Journal of Mechanical Engineering and Technology (IJMET) Volume 9, Issue 5, May 2018, pp. 717 724, Article ID: IJMET_09_05_079 Available online at http://www.iaeme.com/ijmet/issues.asp?jtype=ijmet&vtype=9&itype=5
More informationEntropy 2011, 13, ; doi: /e OPEN ACCESS. Entropy Generation at Natural Convection in an Inclined Rectangular Cavity
Entropy 011, 13, 100-1033; doi:10.3390/e1305100 OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Article Entropy Generation at Natural Convection in an Inclined Rectangular Cavity Mounir
More informationNUMERICAL STUDY OF MELTING OF TIN WITHIN A RECTANGULAR CAVITY INCLUDING CONVECTIVE EFFECTS
NUMERICAL STUDY OF MELTING OF TIN WITHIN A RECTANGULAR CAVITY INCLUDING CONVECTIVE EFFECTS Christiano Garcia da Silva Santim, chrisoff22@yahoo.com.br Luiz Fernando Milanez, milanez@fem.unicamp.br Universidade
More informationc. The Grashof number is the ratio of buoyant forces to viscous forces acting on a fluid.
QUESTION 1. (0 pts) With respect to free convection: a. What is an extensive, quiescent fluid? (4 points) b. What are the two major physical considerations or forces for free convection? (4 points) c.
More informationTransactions on Engineering Sciences vol 18, 1998 WIT Press, ISSN
Simulation of natural convection in a reservoir P. Jelmek*, V. Havlik\ R. Cerny\ P. Pfikryl" * Czech Technical University, Faculty of Civil Engineering, Department of Physics, Thdkurova 7, 166 29 Prague
More informationSolidification of a binary alloy: Finite-element, single-domain simulation and new benchmark solutions
Journal of Computational Physics 216 (2006) 247 263 www.elsevier.com/locate/jcp Solidification of a binary alloy: Finite-element, single-domain simulation and new benchmark solutions Michael Le Bars *,
More informationEFFECT OF VARYING THE HEATED LOWER REGION ON FLOW WITHIN A HORIZONTAL CYLINDER
ISTP-1, 5, PRAGUE 1 TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA EFFECT OF VARYING THE HEATED LOWER REGION ON FLOW WITHIN A HORIZONTAL CYLINDER S. S. Leong School of Mechanical and Manufacturing Engineering
More informationPAPER 331 HYDRODYNAMIC STABILITY
MATHEMATICAL TRIPOS Part III Thursday, 6 May, 016 1:30 pm to 4:30 pm PAPER 331 HYDRODYNAMIC STABILITY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal
More informationBenchmark solutions for twodimensional. transfer problems in irregular regions using multigrid method. Jingfa Li, Bo Yu and Min Wang.
Special Issue Article Benchmark solutions for twodimensional fluid flow and heat transfer problems in irregular regions using multigrid method Advances in Mechanical Engineering 2015, Vol. 7(11) 1 17 Ó
More informationTurbulent Natural Convection in an Enclosure with Colliding Boundary Layers
Turbulent Natural Convection in an Enclosure with Colliding Boundary Layers Abstract Mutuguta John Wanau 1* 1. School of Pure and Applied Sciences, Murang a University of Technology, P.O box 75-10200,
More informationComputational Techniques to Solve Electromagnetic Problems by Using FEM
Computational Techniques to Solve Electromagnetic Problems by Using FEM N.FÜSUN SERTELLER (*), DURSUN ÜSTÜNDAĞ (**) * Electrical Education Department, Faculty of Technical Education, **Mathematics Department,
More informationLeast Squares Finite Element Methods for Large Scale Incompressible Flows
Least Squares Finite Element Methods for Large Scale Incompressible Flows by Tate T. H. Tsang Department of Chemical & Materials Engineering University of Kentucky Lexington, KY 40506 tsang@engr.uky.edu
More informationNumber of pages in the question paper : 06 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationNovel Computational Fluid Dynamics technique for Incompressible flow and flow path design of a novel centrifugal compressor
Novel Computational Fluid Dynamics technique for Incompressible flow and flow path design of a novel centrifugal compressor A thesis submitted to the Graduate School of the University of Cincinnati in
More informationParabolic Flow in Parallel Plate Channel ME 412 Project 4
Parabolic Flow in Parallel Plate Channel ME 412 Project 4 Jingwei Zhu April 12, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of this project is to develop and apply a computer
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationRésonance et contrôle en cavité ouverte
Résonance et contrôle en cavité ouverte Jérôme Hœpffner KTH, Sweden Avec Espen Åkervik, Uwe Ehrenstein, Dan Henningson Outline The flow case Investigation tools resonance Reduced dynamic model for feedback
More informationNumerical Study for Free Convection Heat Transfer inside an Inclined Cavity with Cylindrical Obstacles
International ournal of Engineering and Technology Volume 3 No. 5, May, 13 Numerical Study for Free Convection eat Transfer inside an Inclined Cavity with Cylindrical Obstacles Khudheyer S. Mushatat College
More informationThe solution of the discretized incompressible Navier-Stokes equations with iterative methods
The solution of the discretized incompressible Navier-Stokes equations with iterative methods Report 93-54 C. Vuik Technische Universiteit Delft Delft University of Technology Faculteit der Technische
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationEnhancement of Natural Convection Heat Transfer in a Square Enclosure with Localized Heating from Below
Enhancement of Natural Convection Heat Transfer in a Square Enclosure with Localized Heating from Below Onyango. O. M 1 ; Sigey. J. K 2 ; Okelo. J. A 3, Okwoyo. J. M 4 1, 2, 3 Department of Pure and Applied
More informationBEM for turbulentfluidflow. L. Skerget & M. Hribersek Faculty of Mechanical Engineering, University ofmaribor, Slovenia.
BEM for turbulentfluidflow L. Skerget & M. Hribersek Faculty of Mechanical Engineering, University ofmaribor, Slovenia. Abstract Boundary-Domain Integral Method is applied for the solution of incompressible
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 90 (2014 ) 418 424 10th International Conference on Mechanical Engineering, ICME 2013 Effect of Lewis Number on Unsteady Double
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 5 NATURAL CONVECTION HEAT TRANSFER BASIC CONCEPTS MECHANISM OF NATURAL
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationNatural Convection Heat Transfer from a Rectangular Block Embedded in a Vertical Enclosure
Chapter 5 Natural Convection Heat Transfer from a Rectangular Block Embedded in a Vertical Enclosure Xiaohui Zhang Additional information is available at the end of the chapter http://dx.doi.org/10.5772/52666
More informationMaximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection
Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection Ahmed Waheed Mustafa 1 Mays Munir Ismael 2 AL-Nahrain University College of Engineering Mechanical Engineering Department ahmedwah@eng.nahrainuniv.edu.iq
More informationTowards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below
Copyright 2008 Tech Science Press FDMP, vol.4, no.4, pp.263-269, 2008 Towards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below G. Accary 1, S. Meradji
More informationNATURAL CONVECTION FLOW IN A SQUARE CAVITY WITH INTERNAL HEAT GENERATION AND A FLUSH MOUNTED HEATER ON A SIDE WALL
Journal of Naval Architecture and Marine Engineering December, 2010 DOI: 10.3329/jname.v7i2.3292 http://www.banglajol.info NATURAL CONVECTION FLOW IN A SQUARE CAVITY WITH INTERNAL HEAT GENERATION AND A
More informationEnhancement of the momentum interpolation method on non-staggered grids
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 33: 1 22 Enhancement of the momentum interpolation method on non-staggered grids J. Papageorgakopoulos, G. Arampatzis,
More information