HMT Manual v2 1 1
HMT Manual v2 1 1
CH2092E
Name: __________________________
CH2092E
LABORATORY RECORD
Name: _____________________________________________________________
Roll No.:____________________________________________________________
Certified that this is the bonafide record of the work done in CH2092E Heat and Mass
Transfer Laboratory of this department at National Institute of Technology Calicut, by
Mr/Miss: ____________________________________________
Place: _____________________
Faculty Incharge
(Signature with Date)
INDEX
10
OBJECTIVE
To determine the total resistance of a composite wall and plot the temperature gradient
along the composite wall
INTRODUCTION
If the temperature varies with location in some stationary medium be it a solid or a fluid,
then a “temperature gradient” is said to exist in that medium. Energy in transit due to
such a temperature gradient is called conduction heat transfer. The degree of
temperature variation within a conducting medium is quantified using the mathematical
tools of vector calculus; in particular, the magnitude of the temperature gradient is found
from the derivative of temperature with respect to distance, and the gradient points in the
direction of steepest temperature increase. In this laboratory exercise, you will study
conduction heat transfer in detail, exploring the relationship between heat flux and
temperature gradient, geometrical effects, composite structures, and interfacial
resistance.
This mode of heat transfer is extremely important in industry, with applications in
automotive systems, aerospace systems, chemical and materials processing, electronics
cooling, space conditioning and environmental control, and medicine. In order to
concentrate on the fundamentals of conduction, your experiments will be restricted to
one-dimensional, steady state conditions, where the heat transfer occurs in only a
single coordinate direction and is independent of time.
Heat is always transferred when the temperature difference exists between two bodies.
There are three basic modes of heat transfer:
THEORY
In conduction heat transfer problems, the object being studied is usually a solid.
Convection problems involve a fluid medium. Radiation heat transfer problems involve
either solid or fluid surfaces, separated by a gas, vapor, or vacuum. There are several
ways to correlate the geometry, physical properties, and temperature difference of an
object with the rate of heat transfer through the object. In conduction heat transfer, the
most common means of correlation is through Fourier’s Law of Conduction. The law, in
its equation form, is used most often in its rectangular or cyclic. Let's assume that we
have a combination of different materials put together to form a composite structure like
the composite wall. Let's also assume that the cross- sectional area normal to the flow of
heat transfer is constant and the periphery is insulated so that heat flow is
one-dimensional. Taking only one of the slabs for now, the heat transfer is governed by
Fourier's Law: cyclic form (pipes and cylinders), both of which are presented below.
Fourier’s Law
The fundamental relation between the heat flux, q , and the temperature gradient, ∇T,
q = − kT (1)
Thoughtful inspection of Fourier’s Law can tell you a lot about conduction heat transfer.
Namely, the conductive heat flux is a vector—it has a magnitude and direction. The
magnitude of the conductive heat flux is proportional to the temperature gradient, and its
direction is coincident with the temperature gradient. Since the temperature gradient is
in the direction of maximum temperature change, the heat flux is in the direction of
maximum temperature change; i.e., it is perpendicular to the isotherms. The minus sign
means the conduction heat flux is in the direction of decreasing temperature. Many
engineers regard Fourier’s law as defining the thermal conductivity; this is a healthy
point of view that you may find useful.
In Cartesian (x, y, z) and cylindrical (r, Φ, z) coordinates, Fourier's Law is expressed as,
q q y ˆ qz T ˆ T ˆ T
q = iˆ x + ˆj +k = −kT = −k iˆ + j +k (2)
Ax Ay Az x y z
where qx, is the heat transfer rate in the x-direction, , is the thermal conductivity of the
T
material, Ax is the cross-sectional area of the material normal to the x-direction, and
x
is the temperature gradient in the x-direction.
One-Dimensional Conduction in a Cartesian System
Consider the situation shown in Figure 1, where a cylinder of constant cross-sectional
area has a high temperature T1 imposed on its left end and a lower temperature T2 on its
right end.
Assuming negligible heat loss from the cylinder periphery and that the x-coordinate is
aligned with the cylinder axis, the heat transfer and temperature gradient will occur in
the x-direction only (i.e., heat transfer is by one-dimensional conduction). For steady
state conditions (∂T/∂t = 0) with no heat generation, the heat diffusion equation becomes
dq x
=0 (3)
dx
d dT
=0 (4)
dx dx
The first form states that under these conditions the heat flux, q x is a constant,
independent of x. If the thermal conductivity can be assumed to be constant, the second
form reduces to
d 2T
=0 (5)
dx 2
This can be integrated twice to obtain the general form of the temperature distribution.
T ( x) = C1 x + C2 (6)
Applying boundary conditions T(0) = T1, at the left and T(L) = T2 at the right end gives
T2 − T1
C2 = T1 and C1 = (7)
L
T ( x) − T1 x
= (8)
T2 − T1 L
This result demonstrates the well-known finding that for steady one dimensional heat
conduction through a plane wall with constant thermal conductivity and no volumetric
heat generation, the temperature varies linearly with x.
Now using Fourier’s law to relate the heat flux to the temperature gradient
Q dT − k (T1 − T2 )
q x = = −k = (9)
Ax dx L
qx, is the heat transfer rate in the x-direction , and A, is the cross-sectional area normal to
the x-direction.
Substituting this equation into the solution for T(x), the temperature distribution may
also be expressed in the form
qx
T ( x) = T1 − x (10)
kA
In analogy with Ohm’s law for electrical resistance R elec =∆E/I, the thermal resistance is
defined as the ratio of the driving potential ΔT to the corresponding transfer rate qx. For
steady heat transfer through a plane wall, the thermal resistance for conduction is
therefore given by
T1 −T 2 L
Rcond = = (11)
qx kA
DATA:
Wall thickness of composite slab = 0.045 m
Bakelite:
Diameter = 25 cm
Thickness = 10 mm
Brass:
Diameter = 25 cm
Thickness = 10 mm
Mild Steel:
Diameter = 25 cm
Thickness = 25 mm
OBSERVATION TABLE:
LA L L
Rcond = + B + C
k A A k B A kC A
APPARATUS DESCRIPTION:
The apparatus consists of a central heater sandwiched between two sheets. Three types
of slabs are provided both sides of heater, which forms a composite structure. A small
hand press frame is provided to ensure the perfect contact between the slabs. A
dimmerstat is provided for varying the input to the heater and input measurements are
carried out by a voltmeter and an ammeter. Thermocouples are embedded between
interfaces of the slabs, to read the temperature at the surface.
APPARATUS:
SPECIFICATIONS:
Electricity Supply: 1 Phase, 220 V AC, 2 Amp.
Table for set-up support.
TECHNICAL DETAILS:
Slab assembly arranged symmetrically provided on both sides of heater.
Slab Material and Slab Size:
Cast Iron : 250 mm diameter & 20 mm thick.
Bakelite : 250 mm diameter & 15 mm thick.
Press Wood : 250 mm diameter & 12 mm thick.
Heater : Nichrome wire.
CONTROL PANEL COMPRISING OF
Digital Voltmeter : 0-300 Volts.
Digital Ammeter : 0-2 Amp.
Variac : 0-230 V, 2 Amp.
Digital Temp. Indicator : 0-200 °C, with multi-channel switch
Temperature Sensors : RTD PT-100 type.
EXPERIMENTAL PROCEDURE:
1) See that plates are symmetrically arranged on both sides of the heater plates.
2) Operate the hand press properly to ensure perfect contact between the plates.
3) Close the box by cover sheet to achieve steady environmental
4) Start the supply of heater by varying the dimmerstat. Adjust the input to desired value
5) Take readings of all the thermocouples at an interval of 10 minutes until fairly steady
temperatures are achieved and rate of rise is negligible.
6) Note down the reading in observation table.
PRECAUTIONS:
1) Keep the dimmer stat at zero before switching on the power supply.
2) Do not disturb the thermocouples while testing.
3) Do not increase the power supply above 150 V.
4) Increase the voltage slowly.
MODEL CALCULATIONS:
S. No.: 1
A = π × 0.252 / 4 = 0.04906 m2
= (45-36)/713.4121
LA L L
Equivalent resistance method: Rcond = + B + C
k A A k B A kC A
QL A
For mild steel, k A = W/mK = (35*25*10-3 )/(0.04906) *(45-44)
A(TA − TB )
kA = 17.8353 W/mK
QL B
For Bakelite, k B = W/mK = (35*10*10-3 )/(0.04906) *(44-37)
A(TB − TC )
kB = 1.01916 W/mK
QLC
For Brass, k C = W/mK = (35*10*10-3 )/(0.04906) *(37-36)
A(TC − TD )
kC = 7.1341W/mK
GRAPH:
Plot temperature in y-axis and thickness or length of the wall along x-axis
50
40
TEMPERATURE IN OC
30
20
10
0
0 25 35 45
DISATNCE IN MM
RESULTS:
The heat transfer through a composite wall through conduction has been investigated.
The temperature vs length of the walls has been plotted. The average total resistance of
the given composite wall has been found to be 0.25715 Km2/W.
The thermal conductivity of the three materials namely Bakelite, Brass and Mild Steel
are 1.01916 W/mK, 7.1341 W/mK and 17.8353 W/mK respectively for calculations
number 1.
DISCUSSION:
1. This experiment helps us to determine which material has a higher or lower rate
of heat dissipation.
SAFETY PRECAUTIONS:
1. Avoid direct contact with hot surfaces – The composite wall can reach high
temperatures; use heat-resistant gloves if handling components.
2. Allow components to cool before touching – After the experiment, let the system
cool down before disassembling.
3. Use heat insulation – To minimize heat loss and accidental burns, ensure proper
insulation around the heated sections.
The rate of heat Q0 transferred to the surrounding fluid is proportional to the object's
exposed area A, and the difference between the object temperature Tw and the fluid
free-stream temperature T∞.
The constant of proportionality h is termed the convection heat-transfer coefficient.
Other terms describing h include film coefficient and film conductance. The Dittus –
Boelter correlation (1930) is a common and particularly simple correlation useful for
many applications. This correlation is applicable when forced convection is the only
mode of heat transfer; i.e., there is no boiling, condensation, significant radiation, etc.
The accuracy of this correlation is anticipated to be ±15%.
The Dittus- Boelter equation, has two forms as follows:
Nu = 0.023 Re0.8 Pr0.4, for ‘heating’ (temperature of wall > temperature of fluid)
Nu = 0.026 Re0.8 Pr0.3, for ‘cooling’ (temperature of wall < temperature of fluid)
SPECIFICATIONS:
Test pipe = 33 mm internal diameter, 500 mm long
Diameter of pipe = 22 mm
Cd = 0.64
Cross sectional area of pipe = 8.3 × 10-4 m2
OBSERVATIONS:
S. Hw Voltage Current T1 T2 T3 T4 T5 T6 T7
No (cm) (V) (A) (℃) (℃) (℃) (℃) (℃) (℃) (℃)
PROCEDURE:
1. Switch ON the mains system
2. Switch ON the blower.
3. Adjust the flow using gate valve to a desired difference in the manometer level.
4. Start heating the test section with the help of dimmerstat and adjust the desired
heat input with the help of Voltmeter (125 V) and Ammeter.
5. Take readings of all the seven thermocouples after the steady state is reached.
6. Note down the heater input.
7. Repeat the experiment for different air flow rate or power rates
PRECAUTIONS:
1. Keep the dimmer stat at zero position before switching ON the power supply.
2. Increase the voltmeter gradually.
3. Do not stop the blower in between the testing period.
4. Do not disturb thermocouples while testing. Operate selector switch of the
thermocouple gently. Don’t exceed 200 Watts
5. Operate selector switch of the temperature indicator gently.
ρa0 =
ρa =
Ta = (T1 + T7)/2 =
Ts = (T2 + T3 + T4 + T5 + T6)/5 =
ha = hw × (ρw/ρa) =
ao = π do / 4 =
Q = Cd × ao × √(2gha) (m3/s) =
ma = Q × ρa =
ap = π dp2 / 4 =
v = Q / ap (m/s) =
A = π×D×L (m2) =
ε = 0.4
NOMENCLATURE:
ha Head of air
hw Head of water
do Diameter of orifice = 22 mm
v Velocity of air
Ts Surface temperature
Ta Ambient temperature
Nu Nusselt number
Nu = __________NRe______ NPr______
Dv a
N Re = =
Cp
N Pr = =
k
Nu =
htheo = Nu k / D =
RESULTS:
The heat transfer coefficient has been found from experiments and theoretical
expression of Dittus Boelter equation. The error is about _________ %.
The experimental heat transfer coefficient is about ________________ W/m2K.
DISCUSSION:
INTRODUCTION:
Thermal radiation is electromagnetic radiation generated by the thermal motion of
charged particles in matter. All matter with a temperature greater than absolute zero
emits thermal radiation. The energy transfer by conduction and convection requires
some material carrier, but the transfer of energy by radiation can take place in a vacuum
without a material carrier between the heat source and the receiver.
A black body is an idealized physical body that absorbs all incident electromagnetic
radiation. Because of this perfect absorptivity at all wavelengths, a black body is also the
best possible emitter of thermal radiation, which it radiates incandescently in a
characteristic, continuous spectrum that depends on the body's temperature. At
Earth-ambient temperatures this emission is in the infrared region of the electromagnetic
spectrum and is not visible. The object appears black, since it does not reflect or emit
any visible light.
Absorptivity: If the amounts of radiation energy absorbed, reflected, and transmitted
when radiation strikes a surface are measured in percentage of the total energy in the
incident electromagnetic waves. The total energy would be divided into three groups,
they are called Absorptivity (α), Reflectivity (ρ), and Transmissivity (τ), with α+ρ+τ = 1.
Absorptivity: It is the fraction of irradiation absorbed by a surface.
Reflectivity: It is the fraction reflected by the surface.
Transmissivity: It is the fraction transmitted by the surface.
A body is considered transmit some of the radiation waves falling on its surface. If
electromagnetic waves are not transmitted through the substance it is therefore called
opaque. When radiation waves hit the surface of an opaque body, some of the waves are
reflected back while the other waves are absorbed by a thin layer of the material close to
the surface. For engineering purposes all materials are thick enough to be considered as
opaque, reducing α + ρ = 1.
CH2092E Heat and Mass Transfer Laboratory Page No.
National Institute of Technology Calicut Department of Chemical Engineering
The emissivity of a material (usually written ε or e) is the relative ability of its surface to
emit energy by radiation. It is the ratio of energy radiated by a particular material to
energy radiated by a black body at the same temperature. A true black body would have
ε = 1 while any real object would have ε < 1. Emissivity is a dimensionless quantity.
Stefan Boltzmann law:
The Stefan–Boltzmann law, also known as Stefan's law, states that the total energy
radiated per unit surface area of a black body per unit time (known variously as the
black-body irradiance, energy flux density, radiant flux, or the emissive power) q , is
directly proportional to the fourth power of the black body's thermodynamic
temperature T (also called absolute temperature):
q = AT 4
where,
q = rate of heat transfer (W)
σ = 5.67 × 10-8 W/m2K4, Stefan Boltzmann constant
A = surface area of the emitting body (m2)
The Planck law gives the intensity radiated by a blackbody as a function of frequency (or
wavelength). Let a blackbody have temperature T. Let be the energy density per unit
2h 3
u v () = 3
c exp( h /( kT )) − 1
EXPERIMENTAL APPARATUS:
The apparatus consist of a water heater jacket of hemisphere shape, where a hemisphere
is surrounded by hot water. A copper test disc is fitted at the center of jacket. The hot
water is obtained from hot water tank, fitted to the panel in which water is heated by an
electric heater. The hot water is taken around the hemisphere, so that hemisphere
temperature rises. The test disc is then inserted at the center. Thermocouples are fitted
inside hemisphere to find the average hemisphere temperature. Another thermocouple is
fitted at the center of the test disc to measure the temperature of the test disc. A timer
with a small buzzer is provided to note down the disc temperature at the time interval of
5 sec.
When the blackened disc is inserted at the centre of the hemisphere, heat is transferred to
the disc from hemisphere by radiation and its temperature begins to rise and from
temperature rise rate (it is be the intervals of 5 seconds) Using this data, Stefan –
Boltzmann constant is determined
0 70 140
5 75 145
10 80 150
15 85 155
20 90 160
25 95 165
30 100 170
35 105 175
40 110 180
45 115 185
50 120 190
55 125 195
60 130 200
65 135
PROCEDURE:
1. Ensure that all the valves are closed except the overflow valve on the test setup
2. Ensure that copper disc is pulled down and water is present in the top water tank
3. Heat the water in top water tank to about 80℃ as measured by thermocouple T5
4. Open the gate valve very slowly so that water flows from hot water tank and fills
around the hemisphere. Close the gate valve immediately if hot water overflows
5. Measure the temperatures on hemisphere using thermocouples T1, T2 and T3
6. Change reading to T4 and Push the copper test disc up to reach the centre
7. Note temperature T4 every five seconds until about 200 seconds
8. Pull the copper test disc down to its original position
9. Once the water around the hemisphere has cooled, drain it and refill the top tank
EQUATIONS:
According to Stefan –Boltzmann law, energy radiated by a body per unit area per unit
time is given by
R = εσT4, where R is energy radiated per unit area per time
ε = emissivity of the material
σ = Stefan –Boltzmann constant = 5.67 x 10-8 W/m2K4
T = temperature in Kelvin scale.
For a black body emissivity =1, hence, R = σ T4
In the given experiment set up the net heat transfer to the disc per second is
∆Q/∆t = m CpdT/ dt = σA[Th4 –Td4]
From the above equation,
mC p dT / dt
=
A(Th4 − Td4 )
CALCULATIONS:
T1 + T2 + T3
Hemisphere temperature = + 273.15 = = _______ K
3
mC p dT / dt
= =
A(Th4 − Td4 )
GRAPH:
Plot the temperature T (in ℃) in y-axis and time t (in seconds) in x-axis.
RESULTS:
Heat transfer in a black body radiation setup was studied and the Stefan Boltzmann
DISCUSSION:
SAFETY PRECAUTIONS:
To determine the heat transferred between the hot and cold fluid flowing through a shell
and tube heat exchanger. Determine the internal and external overall heat transfer
coefficients of heat transfer for this setup. Determine the temperature efficiencies.
INTRODUCTION:
A heat exchanger is a piece of equipment built for efficient heat transfer from one
medium to another. The media may be separated by a solid wall, so that they never mix,
or they may be in direct contact. Shell and tube heat exchanger is a class of heat
exchanger designs. It is the most common type of heat exchanger in oil refineries and
other large chemical processes, and is suited for higher-pressure applications. This type
of heat exchanger consists of a shell (a large pressure vessel) with a bundle of tubes
inside it. One fluid runs through the tubes, and another fluid flows over the tubes
(through the shell) to transfer heat between the two fluids. The set of tubes is called a
tube bundle, and may be composed by several types of tubes: plain, longitudinally
finned, etc. In this way, waste heat can be put to use. This is an efficient way to conserve
energy. Heat exchangers with only one phase (liquid or gas) on each side can be called
one-phase or single-phase heat exchangers. Two-phase heat exchangers can be used to
heat a liquid to boil it into a gas (vapor), sometimes called boilers, or cool a vapor to
condense it into a liquid (called condensers), with the phase change usually occurring on
the shell side. There are also single pass, and multi pass exchangers, based on the
number of passes of the shell side or tube side liquids.
We apply the same formula of Q = m Cp ΔT to their shell and tube side fluids.
Qe = mH C pH (THi − THo )
Qa = mC C pC (TCo − TCi )
Here,
Qe = heat emitted by the hot fluid (W)
mH = mass flow rate of the hot fluid (kg/s)
CpH = specific heat of the hot fluid at its average temperature (J/kg K)
THi = inlet temperature of the hot fluid (K)
THo = outlet temperature of the hot fluid (K)
Qa = heat absorbed by the cold fluid (W)
mC = mass flow rate of the cold fluid (kg/sec)
CpC = specific heat of the cold fluid at its average temperature (J/kg K)
TCo = outlet temperature of the cold fluid (K)
TCi = inlet temperature of the cold fluid (K)
The average temperature is:
H + C
m =
2
Qe
U=
ATln
where, A is the area of heat transfer (m2) and Q is heat emitted from the hot fluid (W)
APPARATUS DESCRIPTION:
The apparatus is a shell and tube heat exchanger itself. It has an external heater
attached to it, to heat the hot flow water. The liquids used in both the sides of the
apparatus id water. Thus, water is heated and sent into the hot water inlet at point 2,
shown in Figure 1. Cold water is sent into the cold water inlet at point 1. Similarly, the
hot water is removed from the hot water outlet at point 4. And the cold water is
removed from the cold water outlet at point 3.
The hot water flows in the tubes, and the cold water flows in the shell. The baffles are
provided to promote the counter flow of the shell side liquid. The hot fluid is heated to
a constant temperature by an external heater. The cold fluid is let directly from the tap.
To measure the temperatures, there are four thermo couples at the points 1, 2, 3, 4.
Their readings respectively give the temperatures TCi, THi, TCo, THo. To measure the
flow rate of the hot and cold fluids, two rotameters are fixed to the apparatus. One
measures the flow rate of the cold fluid as it enters, and the other measures the flow
rate of the hot fluid as it leaves. There is also a control valve and a bypass valve to
control the flow of the hot water. There is an indication panel where the readings of the
thermocouples can be read.
PROCEDURE:
1. Fill water in the water bath and switch on the heaters.
2. Adjust the required temperature of hot water.
3. Adjust the control valve. Allow hot water to recycle in bath through by pass
switching on the magnetic pump.
4. Start flow of cold water through annulus and set the exchanger to counter flow.
5. Measure the flow rate of cold water and hot water by using rotameter or beakers.
6. Keeping the flow same, wait till steady state condition is reached (30 minutes).
7. Record the temperature of the hot water and cold water side and also the flow
rate accurately.
8. Now change the flow of the cold water using the control valve.
9. Repeat the procedure again
OBSERVATIONS:
S. Hot Water Cold Water
No V mL t (s) THi ℃ THo ℃ V mL t (s) TCi ℃ TCo ℃
SAMPLE CALCULATIONS:
S. No._____________
ΔTLMTD = = _____________ K
Uo = Qavg/(Ao ΔTLMTD) =
= =_____________ W/m2K
Ui = Qavg/(Ai ΔTLMTD) =
= =_____________ W/m2K
ηH =
ηC =
ηm =
Overall Results
RESULTS:
The heat transfer between the hot and cold fluid flowing through a _______ current
shell and tube heat exchanger has been investigated. Using log mean temperature, the
average internal and external overall heat transfer coefficients were determined to be
DISCUSSION:
SAFETY PRECAUTIONS:
THEORY:
The general boiling is a convection process involving a change is phase from liquid the
vapour. Boiling may occur when a liquid is to contact with a surface maintained at a
temperature higher than the saturation temperature of the liquid of heat is added to a
liquid from a submerged solid surface. The boiling process is referred to as pool boiling.
In this process the vapour produced may form bubbles, which glows and subsequently
detach them from the surface, rising to the free surface due to buoyancy effects.
PROCEDURE:
1. Distilled water is to be taken in the container.
2. The heaters are to be completely submerged.
3. The heater coil and best heater wire are connected across the studs and electrical
connections are made.
4. Heater coil is switched ON and left ON till the required bulk temperature of water
reaches the desired value.
5. Then bulk heating coil is switched OFF and the test heater is switched ON.
6. The boiling phenomenon on wire is observed, by varying the voltage across it using
the dimmerstat.
7. As the voltage is increased further and a point is reached when wire breaks (keep
rigorous note of the current) and at this point the voltage and current are to be noted
8. Change the nichrome wire and repeat with different bulk temperature
Figure
OBSERVATION TABLE:
S. No TWater (℃) Voltage V (V) Current I (A) Q (W) qE (W/m2)
APPARATUS DESCRIPTION:
SPECIFICATIONS:
FORMULA:
( l − v )g
1/ 4
GRAPH:
Plot qE in y-axis and temperature in x-axis
SAMPLE CALCULATION:
S. No.: 3
= 5.7491*10-5 m2
= 27*3.57/(5.7491*10-5 )
RESULTS:
We have successfully studied the boiling heat transfer phenomena and determined the
critical heat flux of pool boiling at 60 ℃, 80 ℃ and 100 ℃. The variation is also
plotted as a graph.
DISCUSSION:
3. As the bult temperature increases, the current required to break the wire
SAFETY PRECAUTIONS:
1. When changing nichrome wire, avoid touching the heating element and hot
water.
2. Turn off main switch before changing the wire.
OBJECTIVE
INTRODUCTION
THEORY
If two gases are inter-diffusing with a continual supply of fresh gas and removal of the
products of diffusion, this diffusion reaches an equilibrium state with constant
concentration gradients. This is known as steady state diffusion. If also there is no total
flow in either direction the rates of diffusion of A and B, N A and NB are equal but have
opposite sign.
As per Dalton’s law, the total concentration of the two components CA and CB is constant
dC A dC
=− B
dx dx
NOMENCLATURE
In the case where gas A is diffusing through stagnant gas, B (non-diffusing) the flow
carries both components in proportions to their partial pressure
The total transfer of A is the sum of this proportion of the flow and the transfer by
diffusion
This is the expression used for the experimental determination of vapour diffusion
coefficients in gases by evaporation from a liquid surface in a narrow bore tube and
measuring the fall of the level of this surface. The distance of the liquid surface below
the open end of the tube is measured before and after evaporation over a definite period.
If the variation in level is small, then arithmetic mean of these two readings is taken as
the value of x. In case there is an appreciable change of level, the value of x is
determined by the integration between the initial and final readings of the level.
DATA
P = ____________________ kN/m2
R = ____________________ J/mol K
MA = ____________________ kg/kmol
ρl = ____________________ kg/m3 at T = ______________
ρl = ____________________ kg/m3 at T = ______________
OBSERVATIONS
T = ___________________ K x0 = ______________ cm
x, cm
S. No θ, s (x0 – x), cm θ/(x0 – x), s/cm
3 decimals
1 300
2 600
3 900
4 1200
5 1500
6 1800
T = ___________________ K x0 = ______________ cm
x, cm
S. No θ, s (x0 – x), cm θ/(x0 – x), s/cm
3 decimals
1 300
2 600
3 900
4 1200
5 1500
6 1800
Usually, x0 will not be measured accurately nor is the effective distance for diffusion, x,
at time θ. However, accurate values of (x – x0 ) are available, hence we rewrite Equation
(2) as:
CALCULATIONS:
C2
ln( VP) = C1 + + C3 ln T + C4T C5 (T in K)
T
where,
C1 = 69.006 C2 = −5599.6
C3 = −7.0985 C4 = 6.2237 x 10-6
C5 = 2
= _____________ kN/m2
= _____________ kN/m2
GRAPHS:
Find the slope S from the plot of θ /(x0 − x) in y-axis vs (x0 − x) in x-axis
D = Const. T1.5/P
D can be determined by drawing a curve of DP vs. T.
DESCRIPTION
The equipment consists of a T tube made of glass, placed in a constant temperature water
bath. The temperature of the bath is controlled by the DTC. Air pump is used to supply
the air, passed through the T tube. Volatile component is filled in the T tube and air
passed over it by the pump and change in the level is seen by the sliding microscope.
PROCEDURE
For T =
P
CT = =
RT
(VP) A
CA = CT =
P
From graph, S =
C B1 = CT =
P − (VP) A
CB2 = C B1 =
P
C B1 − C B 2
C BM = =
C
ln B1
CB2
AC BM
D AB = =
2M A SC ACT
For T =
P
CT = =
RT
(VP) A
CA = CT =
P
From graph, S =
C B1 = CT =
P − (VP) A
CB2 = C B1 =
P
C B1 − C B 2
C BM = =
C
ln B1
CB2
AC BM
D AB = =
2M A SC ACT
RESULTS:
The effect of temperature on diffusion coefficient has been studied. The values of
diffusion coefficient at different temperatures are:
DISCUSSIONS:
SAFETY PRECAUTIONS:
SIMPLE DISTILLATION
OBJECTIVE
To verify Rayleigh’s equation for simple distillation of methanol-water mixture.
To verify material balance equation, FxF = DyD,avg + WxW
APPARATUS
CHEMICALS
1. 50 mole% methanol-water mixture.
2. Pure methanol
THEORY
OBSERVATIONS:
Weight (g)
In the figure, the mole fraction of more volatile component (here A) is used to denote the
compositions of the mixtures. Consider the mixture to be initially in the liquid phase
denoted by point G. On addition of heat, the temperature of the mixture rises. Let us
consider that two components, the state of the system have reached point H when the
vapour pressure of the solution equals the constant pressure exerted. The liquid mixtureis
at the point of vaporization, or the liquid mixture has reached its bubble point. When
more heat is added, the temperature further rises to point P. Now the system exists in two
phases; a liquid phase represented by point L and a vapour phase represented by point K.
The composition of the liquid phase suitably changed so that the vapour pressure of the
solution is equal to the pressure maintained. This is the property made use of in
distillation to effect a separation between A and B. On further addition of heatthe state of
the system reaches point M where the mixture gets completely vaporized. The system at
this state said to beat at its dew point. The bottom line is saturated liquidline and the top
curve is the saturated vapour curve. Vapour phase K and liquid phase L coexist and
hence they said to be in equilibrium. While the mole fraction of A in the liquid phase is
denoted by x, the mole fraction of A in the vapour phase in equilibrium which is denoted
by y*. If there is a separation by distillation, y* has to be different from x.
With reference to the figure, let us consider that the state of the system has reached the
point H. A drop of vapour forms which has compositions represented by J. By some
method this vapour is removed from the vessel in which the equilibrium distillation is
carried out and is condense elsewhere. The composition of liquid phase moves up along
the saturated liquid line. Addition of more heat forms more vapour, but all vapour is
removed from the region of distillation immediately on the formation. The process of
distillation is stopped when the composition of the liquid phase has changed to the extent
required. In practice, this can only be approximated since,
• Removal of differential amounts if vapour makes the process infinitely slow.
• Entrainment cannot be fully avoided.
• There will be cooling and condensing of the vapour before removal from the
region of distillation.
OBSERVATIONS
Volume of Volume Weight of mix Mixture Mole %
Mix
methanol of water + specific Density of
No
(ml) (ml) gravity bottle(g) (g/ml) Methanol
1 0 30 0.000
2 3 27
3 6 24
4 9 21
5 12 18
6 15 15
7 18 12
8 21 9
9 24 6
10 27 3
11 30 0 1.000
VA A
MA
xA = =
V A A VB B
+
MA MB
The amount of liquid that is distilled can be related to the initial and final liquid phase
composition using an equation called Rayleigh’s equation. For a binary mixture of A and
B assume that at any time during the course of distillation there are L moles of liquid in
the still of composition x mole fraction of A and the amount of distillate ( =dL) of mole
fraction y* in equilibrium with the liquid, is formed. Let us denote the corresponding
change in the liquid phase composition by dx.
Let us determine the corresponding material balance for A
Lx = y * L + ( L − dL)( x − dx)
where F and W are the total number of moles of liquid at any instant during distillation
having composition xF and xW respectively. This equation can be used to determine F, W,
xF, and xW for simple distillation when any three of these are known.
PROCEDURE
1. Add 300 ml of given methanol-water mixture in a simple distillation flask after
exactly measuring the volume using a measuring cylinder. Assemble the distillation
apparatus, start circulation of water through condenser and switch on the heating.
2. Take distilled water in one burette and methanol in another. Determine the specific
gravity of pure methanol. Clean and dry a 100 mL conical flask. Prepare a 10 % by
volume mixture of methanol and water in this by mixing 3 mL methanol with 27 mL
water and find out the specific gravity of this mixture.
3. Similarly, find the specific gravity values of several known compositions of
methanol-water mixtures. It is better to do the experiment in increasing percentages
of methanol and to shake the dry bottle between successive densitydeterminations so
that little amount of solution remain in the bottle and will not affect much.
4. When required amount of distillate (50 % of feed volume) has been collected, stop
heating and disassemble the apparatus. Measure the volume of distillate collected in
that of residue remaining using a measuring cylinder. Determine the specific gravity
of this solution also
CALCULATIONS:
Volume of Methanol-Water mixture as feed =
Volume of distillate collected =
Volume of residue after distillation =
Molecular wt. of methanol =
Molecular wt. of water =
Density of water =( - )/25 = ____________ g/mL
Density of methanol =( - )/25 = ____________ g/mL
Density of feed mixture =( - )/25 = ____________ g/mL
Density of distillate =( - )/25 = ____________ g/mL
Density of residue =( - )/25 = ____________ g/mL
The equilibrium relation between Methanol (A) – Water (B) system at 1 atm :
𝟏
xA yA*
𝒚∗ −𝒙
0 0 -
0.02 0.1 12.5
0.04 0.2 6.25
0.06 0.3 4.17
0.09 0.4 3.23
0.14 0.5 2.78
0.22 0.6 2.61
0.3 0.66 2.78
0.4 0.73 3.03
0.5 0.78 3.57
0.6 0.82 4.54
0.7 0.87 5.88
0.8 0.92 8.33
0.9 0.96 16.67
1.0 1.0 -
GRAPHS:
Plot the mixture density in y and methanol mole fraction in x-axis. Find xW, xF and yD,avg.
Plot 1/(y*-x) in y-axis and methanol mole fraction in x-axis. Find area under curve
From graph:
For feed, F
From graph, xF =
= =
For Distillate, D
= =
= =
For Residue, W
From graph, xw =
xw ×Molecular Weight of methanol
mass fraction of methanol in Residue =
xw ×MW of methanol+(1−xw )×MW of water
= =
Material Balance:
xF dx
The area under the curve between xF and xW gives the value of ∫x
W y∗ −x
F
ln ( ) = =
W
Hence, the Rayleigh’s equation and material balance for the system has been verified.
RESULTS:
The Rayleigh’s equation and material balances for the given methanol - water system
has been verified with errors present within allowable limits. The error in material
DISCUSSION:
SAFETY PRECAUTIONS:
ADSORPTION ISOTHERM
OBJECTIVE
To verify the applicability of the Freundlich equation for adsorption of acetic acid on
activated carbon and to determine the value of constants k and n in the equation for
adsorption at room temperature.
APPARATUS
Reagent bottles, conical flasks, filtering funnel, burette, pipette, and volumetric flask
CHEMICALS REQUIRED
Acetic acid, activated carbon, 0.1N NaOH solution, Oxalic acid
THEORY
Adsorption isotherm exploits the ability of certain solids preferentially to concentrate
specific substances from solution to their surfaces. In this, components of either liquid or
gaseous solution may be separated from each other. Some applications of adsorption in
gaseous separations are de-humidification of air and other gases, removal of
objectionable odors of impurities from industrial gases, recovery of valuable solvent
vapours from dilute mixtures with air and other gases and fractionation of mixtures of
hydrocarbon gases containing such substances like methane, ethane, and propylene.
Typical liquid separation is the removal of moisture dissolved in gasoline,
decolourisation of petroleum products and aqueous sugar solution, removal of
objectionable taste and odour from water and fractionation of a mixture of aromatics and
paraffinic hydrocarbons.
Adsorption can be physical or chemical. Physical adsorption or Vander Waal’s
adsorption can be readily reserved and is the result of intermolecular forces of attraction
between solid and the substance adsorbed. The equilibrium partial pressure of the
substance adsorbed equals that of contacting gas phase. This is lower than the vapour
pressure and so adsorption is not simple condensation. By lowering the pressure of gas
phase by raising the temperature, the adsorbate is readily removed or desorbed in
unchanged form, reversible adsorption is seen in liquids as well.
OBSERVATIONS:
Room Temperature = 25 oC
Standardization of NaOH
Indicator: Phenolphthalein
Average V1 = 12.1 ml
Preparation of mixtures
2 20 200 5
3 30 200 5
4 40 200 5
5 50 200 5
Indicator: Phenolphthalein
2 20 180 10 9 8.8
CALCULATIONS
= 0.0825 N
solutions where the fraction of the original solvent, which may be adsorbed, is small. For
dilute solutions, the adsorption isotherm may be described by an empirical expression,
V (C0 − Ce )
qe = = KCe1 / n
m
1
In linear form of Freundlich isotherm, can be written as ln qe = ln k + (n) ln Cewhere
PROCEDURE
1. Record the room temperature.
6. Repeat the same procedure with bottle numbers B2, B3, B4, B5 with the
difference that the volumes of solution pipette out from the flask F1 are 20 ml,
30 ml, 40 ml and 50 ml respectively.
7. Shake the bottles well for one hour using bottle shaker for adsorption to come
into equilibrium.
8. Prepare a standard solution of Oxalic acid (0.1 N) and use this to standardize
the supplied NaOH solution.
9. After 1 h shaking, filter the solution and titrate in the five bottles. Take the
solutions in the burette and titrate it against 10 ml NaOH.
CH2092E Heat and Mass Transfer Laboratory Page No.
National Institute of Technology Calicut Department of Chemical Engineering
= 0.002744 g/ml
= 0.004 g/ml
V2B (C0 − Ce )
qe =
m
=200(0.004-.002744)/5
= 0.05024
ln Ce = ln(0.002744)
= -5.8983
ln qe = ln(0.05024)
= -2.9909
C0 Ce
3
Bottle V/m (cm of (g solute/ (g solute/
solution/ g qe ln Ce ln qe
No ml of ml of
adsorbent) solution) solution)
1 50 0.004 0.002744 0.05024 -5.8983 -2.9909
GRAPH:
CALCULATIONS:
From graph,
Thus, n = 1.70588
Intercept = ln k = -0.85
Thus, k = 0.4274
RESULTS:
n = 1.70588
k = 0.4274
DISCUSSION:
The experiment aims to determine the relationship between the concentration of solute
in the liquid phase and amount absorbed into a solid at equilibrium. Data obtained is
fitted to Freundlich isotherm and has been verified within experimental accuracy.
SAFETY PRECAUTIONS:
1. Alway wear lab coat, safety goggles, and gloves to protect against chemical
3. Review the MSDS Sheet for Acetic acid and use appropriate PPF such as gloves
to prevent skin and eye irritation. In case contact with skin, immediately flush
The rate of drying can be determined for a sample of a substance by suspending it over
an electronic balance in the duct. The weight of the drying sample can then be
measured as a function of time.
NOMENCLATURE:
A = Drying surface area, m2
a0 = Area of orifice, m2
ap = Area of pipe, m2
D = Diameter, m
do = Diameter of orifice, m
dp = Diameter of pipe, m
G = Mass flow rate of air, kg/s
h = Head loss, m
N = Drying rate, kg/m2s
Nc = Constant drying rate, kg/m2s
W = Weight of moist solid, kg
R = Manometric difference, m
S = Weight of dry solid, kg
X = Moisture content of a solid, kg water / kg of dry solid
Xc = Critical moisture content, kg water / kg of dry solid
θ = Time, s
ρa = Density of air, kg/m3
ρw = Density of water, kg/m3
DESCRIPTION
The set up consists of an insulated double wall chamber. Inside the chamber, a tray is
attached directly to an electronic weighing balance, which is fitted on bottom outside
the chamber. Material for drying is placed in the tray and regular loss of weight is
monitored. Air from a blower passes through a heating chamber is used as a drying
agent. Flow control and bypass valve are fitted to regulate the airflow.
UTILITIES REQUIRED
1. Electricity Supply: Single Phase, 220V AC, 50 Hz, 5-15 amp. socket with earth
connection.
2. Bench Area Required: 0.75 m x 2 m.
3. Wet Solid.
EXPERIMENTAL PROCEDURE
1. Record the weight of tray alone. Then dry sand and tray.
2. Place the tray with solid inside the drying chamber.
3. Start the blower and fix the air flow rate and also its inlet temperature by adjusting
the energy regulator.
4. When the desired conditions of temperature and air velocity are reached (in about
10-15 min), weight the sample tray again and put known amount of water (6-8 mL) in
it to give desired initial moisture content.
5. Keep the tray gently in the drying chamber and start the stopwatch
6. Record the balance reading at every 3 minute interval for about 40 minutes.
7. Drying is assumed to be complete when at least 3 consecutive readings are
unchanged. The wet/dry bulb temperatures at the inlet and outlet of the drying chamber
and the air flow rate (manometer reading fixed across the orifice) are recorded at least
three times during the course of a run to give average operating conditions.
8. The same steps are repeated for other runs at different operating conditions like air
flow rate.
do =
dp =
d =
S =
Cd =
TDry,in =
TDry,out =
TWet,out =
ρw (at T = ________) =
ρa (at T = ________) =
TABLE:
R = _____ cm
No. θ, s W, g X ΔX Δθ, s N, kg/m2s
1
10
11
12
13
14
R = _____ cm
No. θ, s W, g X ΔX Δθ, s N, kg/m2s
1
10
11
12
13
14
GRAPHS:
Plot the drying rate N (kg/m2s) vs moisture content X (kg of water/kg of dry solid) i.e.
drying rate curve. Mark the critical moisture content, Xc where drying rate starts to fall
CALCULATIONS:
ao = πdo2/4 =
ap = πdp2/4 =
A = πd2/4 =
H = R w − 1 =
a
C d ao a p a 2 gH
G= =
a 2p − ao2
W −S
X = =
S
S X
N =− =
A
RESULTS:
The drying rate curve has been plotted where the different regions of constant drying
region, varying drying regions, critical moisture content and equilibrium moisture
content have been marked.
The critical moisture content is ___________________ at _____________________
DISCUSSION:
SAFETY PRECAUTIONS:
(a) (b)
Figure 2: (a) Type I system and (b) Type II system
R lineME x E − x M
= = (1)
E lineRM xM − x R
where xE, xM, and xR represent the compositions of any one component in mixture E,
M, and R respectively. Similarly, if a mixture at M has been removed from it a mixture
of composition E, the new mixture is on the straight-line EM extended in the direction
away from E and the location of this at R is given by Equation 1.
CH2092E Heat and Mass Transfer Laboratory Page No.
National Institute of Technology Calicut Department of Chemical Engineering
OBSERVATIONS:
No. Benzene Acetic Water No. Benzene Acetic Water
(ml) Acid (ml) (ml) (ml) Acid (ml) (ml)
1 6.6 9 1 10 1 9 0.1
2 3.6 8 2 11 2 8 0.2
3 1.7 7 3 12 3 7 0.3
4 1.4 6 4 13 4 6 0.4
5 1.2 5 5 14 5 5 0.6
6 0.9 4 6 15 6 4 0.9
7 0.6 3 7 16 7 3 1.5
8 0.4 2 8 17 8 2 2.6
9 0.3 1 9 18 9 1 6.1
CALCULATION:
= 1.0156 g/mL
= 1.0632 g/mL
= 0.8796 g/mL
PROCEDURE
TABLE:
No. Weight Fraction No. Weight Fraction
Sample Calculations:
For Mixture 1 ( S.No. : 9):
GRAPH:
Plot the weight fractions on a ternary graph with different markers for 1-9 and 10-18.
RESULTS:
The binodal solubility curve for benzene - water - glacial acetic acid system has been
drawn at room temperature and pressure. From the shape of the binodal curve, we
observe that the system is of Type I, where 2 components are only partially miscible,
and 3rd is miscible with both. The binodal solubility curve exhibits a teardrop shape,
with the immiscibility region lying inside the curve.
DISCUSSION:
1. The experiment identified the immiscibility region, plotted the binodal curve,
petrochemicals.
SAFETY PRECAUTIONS:
1. Wear gloves, goggles, and lab coats to prevent exposure to toxic or corrosive
solvents. Work in a well-ventilated area or fume hood.
2. Symbols on Benzene and Acetic Acid bottles indicate that both are flammable,
and they are toxic and cause skin irritation respectively.
4. Keep flammable solvents away from heat sources. Use spill trays and have a fire
extinguisher ready.