Countercurrent Double-Pipe - Maple Programming Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Applications and Example Worksheets : Science and Engineering : applications/CountercurrentHeatExchanger

Countercurrent Double-Pipe

Heat Exchanger

Introduction

This application models the temperature dynamics of a countercurrent double-pipe heat exchanger. Three partial differential equations describe the following:

 

• 

The heat balances across the tube- and shell-side liquids

• 

A heat balance across the tube-wall (taking into account the heat flow from the shell- and tube-side liquids, and conduction along the length of the tube)

 

The equations are solved numerically, and the temperature profiles are plotted.  The heat exchanger is assumed to be perfectly insulated. Densities, specific heat capacities, heat transfer coefficients, and thermal conductivities are assumed to be constant.

 

Parameters

Specific heat capacities for the tube-side fluid, shell-side fluid, and the tube wall:

Cpt4085:Cps4186:Cpw380:

Densities:

rhot800:rhos1200:rhow8000:

Flow rates:

Ft1:Fs1:

Heat transfer coefficients:

Ut40000:Us40000:

Thermal conductivity of tube-wall:

kw109:

Length of heat exchanger:

L1:

Internal and external diameter of inner tube, and internal diameter of outer tube:

Di0.05:Do0.06:Dis0.1:

Partial Differential Equations Derived from Energy Balances

Tube wall heat balance:

pde1π4Do2Di2CpwrhowtTwx,t=Ut π Di Ttx,tTwx,tUs π DoTwx,tTsx,t+kw π4Do2Di22x2Twx,t:

Tube-side heat balance:

pde2rhot Cpt π4Di2tTtx,t=CptFtxTtx,tπDiUtTtx,tTwx,t:

Shell-side heat balance:

pde3rhos Cps π4Dis2Do2tTsx,t=Cps Fs xTsx,t+πDoUsTwx,tTsx,t:

Initial and boundary conditions:

ibcTsx,0=300,TsL,t=300,Ttx,0=360,Tt0,t=360,Twx,0=330,D1Tw0,t=0,D1TwL,t=0:

Solution and Results

solpdsolvepde1,pde2,pde3,ibc,numeric,time=t,range=0.. L:

p1sol:-plotTt,x= 0.5 L,t=0..5,axes=boxed,color=black,legend=Tube-side liquid:p2sol:-plotTs,x=0.5 L,t=0..5,axes=boxed,color=blue,legend=Shell-side liquid:p3sol:-plotTw,x=0.5 L,t=0..5,axes=boxed,color=red,legend=Tube wall: plotsdisplayp1,p2,p3,labels=Time (s),Temperature (C),title=Temperature Halfway along Heat Exchanger,titlefont=Helvetica, 18,Bold

p1sol:-plotTt,t=5,x=0..L,axes=boxed,color=black,legend=Tube-side liquid:p2sol:-plotTs,t=5,x=0..L,axes=boxed,color=blue,legend=Shell-side liquid:p3sol:-plotTw,t=5,x=0..L,axes=boxed,color=red,legend=Tube wall:plotsdisplayp1,p2,p3,labels=Distance (m),Temperature (C),title=Temperature Profile along Heat Exchanger,titlefont=Helvetica, 18,Bold