In the previous section we considered problems where the

Similar documents
An Exact Solution of Navier Stokes Equation

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

15 Solving the Laplace equation by Fourier method

8 Separation of Variables in Other Coordinate Systems

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

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

On the integration of the equations of hydrodynamics

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

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

you of a spring. The potential energy for a spring is given by the parabola U( x)

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

Chapter 3 Optical Systems with Annular Pupils

Math 124B February 02, 2012

1 Similarity Analysis

Applied Aerodynamics

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

Quantum Mechanics II

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

Numerical Integration

One-Dimensional, Steady-State. State Conduction with Thermal Energy Generation

NATURAL CONVECTION HEAT TRANSFER WITHIN VERTICALLY ECCENTRIC DOMED SKYLIGHTS CAVITIES

Review: Electrostatics and Magnetostatics

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.

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

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

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

1D2G - Numerical solution of the neutron diffusion equation

N igerian Journal of M athematics and Applications V olume 24, (2015),

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

FOURIER-BESSEL SERIES AND BOUNDARY VALUE PROBLEMS IN CYLINDRICAL COORDINATES

Substances that are liquids or solids under ordinary conditions may also exist as gases. These are often referred to as vapors.

SCRAPED SURFACE HEAT EXCHANGER (SSHE)

f(k) e p 2 (k) e iax 2 (k a) r 2 e a x a a 2 + k 2 e a2 x 1 2 H(x) ik p (k) 4 r 3 cos Y 2 = 4

4.4 Buoyant plume from a steady heat source

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

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

ASTR415: Problem Set #6

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

To investigate Heat Loss of a Fluid flowing through a Pipeline for Turbulent Flow

2 Governing Equations

ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS

x x2 2 B A ) v(0, t) = 0 and v(l, t) = 0. L 2. This is a familiar heat equation initial/boundary-value problem and has solution

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

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

THERMODYNAMIC OPTIMIZATION OF TUBULAR HEAT EXCHANGERS BASED ON MINIMUM IRREVERSIBILITY CRITERIA

PHYS 301 HOMEWORK #10 (Optional HW)

CBE Transport Phenomena I Final Exam. December 19, 2013

LECTURER: PM DR MAZLAN ABDUL WAHID PM Dr Mazlan Abdul Wahid

COMPARISON OF METHODS FOR SOLVING THE HEAT TRANSFER IN ELECTRICAL MACHINES

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

NUMERICAL SIMULATION OF FLUID FLOW IN ENCLOSED ROTATING FILTER AND DISK

Chemical Engineering 412

CHAPTER 25 ELECTRIC POTENTIAL

=0, (x, y) Ω (10.1) Depending on the nature of these boundary conditions, forced, natural or mixed type, the elliptic problems are classified as

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

Analytical and numerical solutions for torsional flow between coaxial discs with heat transfer

Lecture 8 - Gauss s Law

MAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS

Chapter 7-8 Rotational Motion

T x. T k x. is a constant of integration. We integrate a second time to obtain an expression for the temperature distribution:

3.1 Random variables

Boundary Layers and Singular Perturbation Lectures 16 and 17 Boundary Layers and Singular Perturbation. x% 0 Ω0æ + Kx% 1 Ω0æ + ` : 0. (9.

2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum

Goodness-of-fit for composite hypotheses.

A Relativistic Electron in a Coulomb Potential

Fluid flow in curved geometries: Mathematical Modeling and Applications

Unobserved Correlation in Ascending Auctions: Example And Extensions

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden

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

Introduction to Arrays

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects

Chemical Reaction and Thermal Diffusion Effects on Mass Transfer Flow through an Inclined Plate

Chem 453/544 Fall /08/03. Exam #1 Solutions

A Backward Identification Problem for an Axis-Symmetric Fractional Diffusion Equation

2. Plane Elasticity Problems

The Effect of Density Variation on Heat Transfer in the Critical Region

Dymore User s Manual Two- and three dimensional dynamic inflow models

MATH 417 Homework 3 Instructor: D. Cabrera Due June 30. sin θ v x = v r cos θ v θ r. (b) Then use the Cauchy-Riemann equations in polar coordinates

The second law of thermodynamics - II.

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

Axial Magnetic Field Effect on Taylor-Couette Flow

Lecture 24 Stability of Molecular Clouds

Conservative Averaging Method and its Application for One Heat Conduction Problem

-Δ u = λ u. u(x,y) = u 1. (x) u 2. (y) u(r,θ) = R(r) Θ(θ) Δu = 2 u + 2 u. r = x 2 + y 2. tan(θ) = y/x. r cos(θ) = cos(θ) r.

The Substring Search Problem

Analytical solutions to the Navier Stokes equations

Parallel Plate Heat Exchanger

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms

Liquid gas interface under hydrostatic pressure

Axisymmetric Stokes Flow past a Swarm of Porous Cylindrical Shells

$ i. !((( dv vol. Physics 8.02 Quiz One Equations Fall q 1 q 2 r 2 C = 2 C! V 2 = Q 2 2C F = 4!" or. r ˆ = points from source q to observer

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

Scaling laws for dynamos in rotating spherical shells and applications to planetary magnetic fields

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

Gauss s Law Simulation Activities

Thermodynamic Head Loss in a Channel with Combined Radiation and Convection Heat Transfer

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.

DOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS

Thermal and Solutal Marangoni Mixed Convection Boundary Layer Flow

Physics 2212 GH Quiz #2 Solutions Spring 2016

Transcription:

5.4 Hydodynamically Fully Developed and Themally Developing Lamina Flow In the pevious section we consideed poblems whee the velocity and tempeatue pofile wee fully developed, so that the heat tansfe coefficient was constant with distance along the pipe. In this section we conside poblems in which the velocity is fully developed at the point whee the heat tansfe stats. Futhemoe, as befoe, we conside two cases of constant wall tempeatue and wall heat flux, both by assuming unifom tempeatue at the inlet. Unde these conditions the heat tansfe coefficient is not constant and vaies along the tube. 1

Whiteman and Dake (1980), Lyche and Bid (1956) and Blackwell (1985) studied the case of fully developed flow with themal enty effects fo non-newtonian fluids. Sellas et all (1959) obtained themal enty length solutions fo the case of a Newtonian fluid with constant wall tempeatue and fully developed flow, which ae pesented below.

We make the following assumptions in ode to obtain a closed fom solution fo heat tansfe analysis in a cicula tube. 1. Incompessible Newtonian fluid. Lamina flow 3. Two-dimensional steady state 4. Axial conduction and viscous dissipation is neglected 5. Constant popeties p This does not mean that one cannot obtain analytical solutions when one o moe of the above assumptions is valid, but the solution will be much easie by making the above assumptions. 3

Since the fully developed velocity was aleady obtained in Section 5., we will focus the emphasis on the solution of the enegy equation and bounday conditions fo a developing tempeatue pofile. 4

Constant Wall Tempeatue The descibing dimensionless enegy equation (5.33) and bounday conditions with making the assumptions noted above fo the case of constant wall tempeatue ae u 1 x,0 1 1, x 0 (5.6) 0, x finite o 0, x 0 (5.63) whee T Tw u x/ 0,, u, x T T u Re P o in w m 5

Fo a fully developed lamina flow, the paabolic velocity pofile developed befoe is applicable. u u 1 o 1 m u o Substituting the above equation into the enegy equation (5.6), we get 1 1 (5.64) x Since the above patial diffeential equation is linea and homogeneous, one can apply the method of sepaation of vaiables. The sepaation of vaiables solution is assumed of the fom, x R X x (5.65) 6

The substitution of the above equation into equation (5.64) yields two odinay diffeential equations (5.66) whee X X 0 1 R R R 1 0 dx d X dr d R X, X, R, R dx dx d d and λ is the sepaation constant o eigenvalue. (5.67) The solution fo equation (5.66) is a simple exponential function of the fom x e while the solution of equation (5.67) is of infinite seies efeed to by Stum-Liouville theoy. 7

The final solution can be of the fom, n n exp n n 0 x c R x (5.68) whee λ n ae eigenvalues and R n ae eigenfuncitons it coesponding to equation (5.67) and c n ae constants. 8

The local heat flux, dimensionless mean tempeatue, local Nusselt numbe and mean Nusselt numbe can be obtained fom the following equations, using the tempeatue distibution above T Tw Tin qw k k (5.69) Nu x o k T T G x o o 1 w in n exp n n0 T exp m T n x w 8 m Gn Tin Tw n0 n exp Gn n x hx o qw o n0 k T w Tin k m m 1 G exp x / n0 1 1 G exp x o Nu xdx hx Nu ln 8 k x 0 x 1 whee G n c nr n 1 m x n n n0 n n n n (5.70) (5.71) (5.7) 9

The fist five tems in equations (5.69) and (5.7) ae fully sufficient i to povide accuate solutions to the above infinite seies. The eigenvalues, λ n and G n, to calculate q w, θ m, Nu x and Nu m fo the above poblem ae pesented in Table 5.3. Table 5.4 povides the vaiation of Nu x, Nu m and θ m with distance along the tube. It can be easily obseved fom Table 5.3 that the fully developed tempeatue pofile stats at appoximately x x / 0 Re P (5.73) Theefoe, (L T,T / D) = 0.05ReP whee L T,T is the themal entance length fo constant wall tempeatue. The themal enty length inceases as the Reynolds numbe and Pandtl numbe incease. A vey long themal enty length is needed d fo fluids with a high Pandtl numbe, such as oil. Theefoe, cae should be taken to make a fully developed tempeatue pofile assumption fo fluids with a high Pandtl numbe. 0.1 10

Table 5.3 Eigenvalues and Eigenfunctions of a Cicula Duct; Themal Enty Effect with Fully Developed Lamina Flow and Constant Wall Tempeatue (Blackwell 1985) n λ n / G n 0 3.656 0.749 1.31 0.544 56.9 0.463 3 107.6 0.414 4 174.5 0.38 11

Table 5.4 Nusselt Solution fo Themal Enty Effect of a Cicula Tube fo Fully Developed Lamina Flow and Constant Wall Tempeatue (Blackwell 1985) x + Nu x Nu m θ m 0 1 0.00050005 10.11 15.4 0.940 0.00 8.06 1. 0.907 0.005 6.00 8.94 0.836 0.0 4.17 5.8 0.68 0.04 3.79 4.89 0.457 0.05 3.71 4.64 0.395 0.1 3.658 4.16 0.190 3.657 3.657 0 1

Constant Heat Flux at the Wall The lamina fully developed d themal enty length (developing tempeatue pofile) fo constant wall heat flux is vey simila to constant wall tempeatue, except that dimensionless tempeatue and bounday conditions ae defined as follows Tin T q D / k w (5.74) q k T At w constant (5.75) o 13

Siegel et al. (1958) solved the above poblem fo lamina fully developed flow using sepaation of vaiables and the Stum-Liouville theoy of which the esult is pesented below * 7, x 4x 4 4 *, x c exp nrn n x n1 (5.76) (5.77) The eigenvalues λ n, eigenfunctions R n and constants c n ae pesented in Table 5.5 14

Table 5.5 Eigenvalues and Eigenfunctions fo Themal Enty Effect of a Cicula Tube fo Fully Developed Lamina Flow and Constant Wall Heat Flux (Siegel et al. 1958) n λ n R n (1) C n 1 5.6796-0.49517 0.403483 83.8618 0.395508-0.175111 3 174.167-0.34587 0.105594 4 96.536 0.314047-0.073804 5 450.947-0.915 0.05503570550357 6 637.387 0.73808-0.043483 7 855.850850-0.5985 0.035597035597 15

The local l Nusselt numbe, based on the above solution, is given below and numeical values ae pesented in Table 5.6 48 /11 Nux (5.78) 1 4/11 c exp x R 1 n 0 n n n The themal entance fo constant wall heat flux based on the numeical esults pesented in Table 5.6 is o x 0.05 LT, H 0.05 ReP D (5.79) whee L T,H is the themal enty length fo fully developed flow with constant wall heat flux. 16

Table 5.6 Nusselt Numbe fo Themal Enty Effect with Fully Developed Flow of a Cicula Tube with Constant Wall Heat Flux (Siegel et al. 1958) x + Nu x 0 0.0013 11.5 0.005 9.0 0.005 7.5 0.01 6.1 0.05 5.0 005 0.05 45 4.5 0.1 4.364 4.364 17

The mean tempeatue vaiation can be obtained fom Nusselt numbe, equation (5.78), by using the following equation: T w T w m hx q q D w h Nu k x (6.80) 18