A spherical vessel has insulation extending from π1=0.50Β m
to π2=0.55Β m. Take the insulation conductivity as
π=0.040Β Wβmβ1βKβ1. Its inner and outer insulation
surface temperatures are 80βΒ°C and
20βΒ°C. Assume steady radial conduction and neglect
heat flow through supports and pipe connections. What is the heat loss
through the insulation?
Worked solution
The given temperatures are at the insulation surfaces, so the required
resistance is the insulation resistance alone. From Eq. 7,
If the given temperatures were fluid temperatures, the fluid-film
resistances would also be needed. For a vessel containing a cold fluid,
the same resistance method gives the inward heat gain.
Consider a 0.8-m-high and 1.5-m-wide glass window with a thickness of 8
mm and a thermal conductivity of π=0.78Β W/(mβ Β°C).
Determine the steady rate of heat transfer through this glass window and
the temperature of its inner surface for a day during which the room is
maintained at 20Β Β°C while the temperature of the outdoors
is β10Β Β°C. Take the heat transfer coefficients on the
inner and outer surfaces of the window to be
β1=10Β W/(m2β Β°C) and
β2=40Β W/(m2β Β°C), which includes the effects of
radiation.
Worked solution
Heat passes from the room air to the inner glass surface, through the
glass, and then from the outer glass surface to the surroundings. These
three steps form a series resistance network.
Figure 8: Heat flows from room air at 20 degrees Celsius through the inside film, glass and outside film to air at minus 10 degrees Celsius. The two glass surface temperatures are the unknown intermediate temperatures.
Data
Quantity
Symbol
Value
Window height
π»
0.8Β m
Window width
π
1.5Β m
Glass thickness
πΏ
8Β mm=0.008Β m
Glass thermal conductivity
π
0.78Β Wβmβ1βKβ1
Room air temperature
πβ,1
20βΒ°C
Outside air temperature
πβ,2
β10βΒ°C
Inside heat transfer coefficient
β1
10Β Wβmβ2βKβ1
Outside heat transfer coefficient, including radiation
β2
40Β Wβmβ2βKβ1
A temperature difference of 1βΒ°C equals 1Β K.
We can therefore express the supplied conductivity and heat transfer
coefficients using kelvin-based units in the calculation.
We need the heat loss Λπ and the inside glass surface temperature
ππ ,in.
Assumptions
Steady state.
Unidirectional heat flow, normal to the glass.
Heat transfer through the edges and frame is neglected.
No heat generation within the glass.
Constant thermal conductivity and heat transfer coefficients.
Uniform temperature on each glass surface.
The supplied outside heat transfer coefficient includes radiation, so no
separate radiation resistance is added.
Calculations
First calculate the area normal to heat flow:
π΄=π»π
Substituting the window dimensions gives
π΄=0.8Γ1.5=1.2Β m2
The inside film resistance is
π in=1β1π΄
Using the inside coefficient and window area gives
π in=110Γ1.2=0.083333Β KβWβ1
For conduction through the glass,
π π=πΏππ΄
Substituting the thickness, conductivity and area gives
π π=0.0080.78Γ1.2=0.008547Β KβWβ1
The outside film resistance is
π out=1β2π΄
Substituting the outside coefficient gives
π out=140Γ1.2=0.020833Β KβWβ1
The same heat rate passes through all three resistances, so they add in
series:
π total=π in+π π+π out
Substituting the three resistances gives
π totalβ0.112714Β KβWβ1
The temperature difference between the air on either side is
Ξπ=πβ,1βπβ,2
Substituting the air temperatures gives
Ξπ=20β(β10)=30Β K
The heat loss can now be calculated using Eq. 2. Writing that
relation for this window gives
Λπ=πβ,1βπβ,2π total
Substituting the temperature difference and total resistance gives
Λπ=300.112714β266.16Β W
Heat flows from the warmer room to the colder outdoors. To find the
inside glass temperature, consider the temperature drop across the
inside film:
πβ,1βππ ,in=Λππ in
The calculated heat rate gives an inside-film temperature drop of
Ξπin=266.16Γ0.083333β22.180Β K
Rearranging the inside-film relation gives
ππ ,in=πβ,1βΛππ in
Subtracting the temperature drop from the room temperature gives
ππ ,in=20β22.180=β2.180βΒ°C
We can check this result by working from the outside air towards the
room. The outer glass surface must be warmer than the outside air by
Λππ out:
The temperature drop through the glass is therefore approximately
β2.180β(β4.455)=2.275Β K.
Dividing this drop by π π
recovers the heat rate of approximately 266Β W.
The surface temperatures also have the expected order: temperature
falls from the room air to the inner glass surface, then to the outer
surface and finally to the outside air.
Answer
The heat loss from the room through the window is
Λπβ266Β W
The inside glass surface temperature is
ππ ,inββ2.18βΒ°C
Most of the resistance is in the two fluid films, with the largest drop
across the inside film. This explains why the glass surface is much
colder than the room air. Adding a trapped air layer, as in double
glazing, increases the resistance and reduces heat loss at the same air
temperatures. The calculated inside surface is below freezing, so frost
may form if sufficient moisture reaches it.
A wire has radius π =0.002Β m, length πΏ=0.50Β m,
thermal conductivity π=15Β Wβmβ1βKβ1 and electrical
conductivity ππ=106Β Sβmβ1. Its current density is
πΌ=106Β Aβmβ2, and its surface is held at
π0=25Β Β°C. Use the wire model above to calculate the
centre temperature, mean temperature and heat loss through the curved
surface. Verify that the heat loss equals the total heat generation.
Calculation and answer
The electrical current produces heat at the volumetric rate
ππ=πΌ2/ππ=106Β Wβmβ3.
From Eq. 11, the centre-to-surface temperature difference is
106(0.002)24Γ15=0.066667Β K.
The mean rise from
Eq. 12 is half this value, 0.033333Β K. Adding
these rises to π0 gives πmaxβ‘=25.0667Β Β°C and
Β―Β―Β―Β―π=25.0333Β Β°C.
The surface flux is
ππ(π )=106Γ0.0022=1000Β Wβmβ2.
Multiplying by the curved area
2ππ πΏ=0.00628319Β m2
gives
Λππ =6.28319Β W.
Independently, the total heat generation is
ππππ 2πΏ=6.28319Β W. This equals the surface heat loss,
confirming the energy balance. The wire loses approximately
6.28Β W at steady state.
A furnace wall separates hot gas at 600Β Β°C from air
at 25Β Β°C. From the hot side, the layers are
0.12Β m of firebrick, 0.060Β m of insulation
and 0.0020Β m of steel. Their thermal conductivities are
1.2, 0.050 and 16Β Wβmβ1βKβ1, respectively.
The hot-side and cold-side heat transfer coefficients are 25 and
10Β Wβmβ2βKβ1, respectively.
Assume steady, one-dimensional heat transfer, constant properties,
perfect contact and negligible radiation. Calculate the heat flux and
all surface and interface temperatures. Determine which solid layer has
the largest temperature drop. Then calculate the heat flux if the
insulation thickness is doubled, with all other conditions unchanged.
Worked solution
Heat passes through the hot-side film, firebrick, insulation, steel and cold-side film in series. We need the heat flux and the temperature at every solid interface.
Figure 10: The furnace wall has three solid layers between hot gas and outside air. A common heat flux passes through all five resistance elements.
Data
Quantity
Symbol
Value
Hot-gas and cold-air temperatures
πβ,ππ
600,25βΒ°C
Firebrick thickness and conductivity
πΏπ,ππ
0.12Β m,1.2Β Wβmβ1βKβ1
Insulation thickness and conductivity
πΏπ,ππ
0.060Β m,0.050Β Wβmβ1βKβ1
Steel thickness and conductivity
πΏπ ,ππ
0.0020Β m,16Β Wβmβ1βKβ1
Hot-side and cold-side coefficients
ββ,βπ
25,10Β Wβmβ2βKβ1
Assumptions
Steady one-dimensional heat flow with constant properties.
Perfect contact between the layers.
No generation, edge loss or separate radiative transfer.
Calculations
We can calculate the heat flux directly by using π β²β²=π΄π , with units
of m2βKβWβ1. The temperature drop is then
Ξπ=πβ²β²π β²β². For each solid layer, π β²β²π=πΏπ/ππ; for each
film, π β²β²β=1/β. In order from hot gas to cold air, these values are
Element
Resistance per unit area, m2βKβWβ1
Hot gas film, 1/ββ
0.040
Firebrick, πΏπ/ππ
0.100
Insulation, πΏπ/ππ
1.200
Steel, πΏπ /ππ
0.000125
Cold air film, 1/βπ
0.100
The total is π β²β²total=1.440125Β m2βKβWβ1.
Dividing the gas-to-air temperature difference by this resistance gives
πβ²β²=600β251.440125=399.3Β Wβmβ2
The temperature drop across each element is πβ²β²π β²β²π. Subtracting
these drops successively from the hot-gas temperature gives
Hot gas: 600Β Β°C.
Hot firebrick surface: 584.0Β Β°C.
Firebrick-insulation interface: 544.1Β Β°C.
Insulation-steel interface: 64.98Β Β°C.
Outer steel surface: 64.93Β Β°C.
Cold air: 25Β Β°C.
The insulation has the largest temperature drop, 479.1Β K.
The drops through the firebrick and steel are 39.93Β K and
0.0499Β K, respectively. Since the heat flux is the same
through all three layers, these differences are determined by πΏπ/ππ.
The insulation has the largest resistance because of its combination
of thickness and low conductivity.
Doubling the insulation thickness adds another
1.200Β m2βKβWβ1 to the total resistance. The new heat
flux is
The heat loss decreases by 45.5%. It does not halve because the films,
firebrick and steel still contribute the same resistances. For the
original wall, adding all five temperature drops recovers the full
575Β K difference between the gas and air.
Answer
The outward heat flux is 399Β Wβmβ2. The solid temperatures from the hot side are 584.0, 544.1, 64.98 and 64.93βΒ°C. Doubling the insulation thickness reduces the flux to 218Β Wβmβ2, a reduction of 45.5 percent.
The two cold-side interface temperatures retain guard digits to show the small contact-temperature drop.
A thin spherical vessel of radius 0.50Β m holds boiling
liquid nitrogen at 77Β K. Its shell is at the liquid
temperature, and the surrounding air is at 293Β K. Insulation
with conductivity 0.040Β Wβmβ1βKβ1 extends from
the vessel surface to an outer radius of 0.54Β m.
Take the outside heat transfer coefficient as
10Β Wβmβ2βKβ1 for both insulated and bare vessels.
Neglect shell resistance, radiation and heat transfer through supports.
The liquid remains at its boiling temperature.
Calculate the heat leak into the insulated vessel, the nitrogen boil-off
rate and the outer insulation temperature. Compare the heat leak and
boil-off rate with those of the bare vessel. The latent heat is
198Β kJβkgβ1. Use the supplied outside coefficient;
no convection correlation is required.
Worked solution
Warm air supplies heat through the outside film and spherical insulation. That heat evaporates nitrogen at its maintained boiling temperature.
Figure 11: The cold spherical vessel has inner insulation radius ri and outer radius ro. Heat flows inward from ambient air through the outside film and insulation.
Data
Quantity
Symbol
Value
Vessel and outer insulation radii
ππ,ππ
0.50,0.54Β m
Liquid and air temperatures
ππ,πβ
77,293Β K
Insulation conductivity
π
0.040Β Wβmβ1βKβ1
Outside coefficient
β
10Β Wβmβ2βKβ1
Nitrogen latent heat
βππ
198000Β Jβkgβ1
Assumptions
Steady heat transfer with constant properties.
Negligible shell resistance, support losses and radiation.
The liquid remains at its boiling temperature.
The same supplied outside coefficient is used for the insulated and bare comparisons.
Calculations
Heat flows from the warmer air through the outside film and then through
the insulation to the nitrogen. The film and spherical conduction
resistances are in series. Using Eq. 7 gives
Λπππ=293β77π πππ +π ππππ=670.8Β W
The outside film uses the area at π=0.54Β m. The insulation
resistance uses both radii to account for the changing area through the
layer. A plane-wall expression with a single area would only approximate
the spherical conduction resistance.
The heat entering the vessel evaporates nitrogen while its temperature
remains constant. Dividing the heat rate by the latent heat gives
between the air and
liquid temperatures as expected for inward heat flow.
For the bare vessel, there is no insulation resistance. The outside film
now acts over the spherical area at π=0.50Β m, giving
Λπππππ=10(4π)(0.50)2(293β77)=6786Β W
The resulting boil-off rate is
Λπππππ=123.4Β kgβhβ1
The insulation reduces both heat leak and boil-off by 90.1% for the
specified conditions. This comparison uses the same β in both cases.
If β were calculated from a convection correlation, changes in diameter
and surface temperature could also affect its value.
Answer
With insulation, the heat leak is 671Β W, nitrogen boil-off is 12.2Β kgβhβ1 and the outside insulation surface is at 274.7Β K. The bare vessel receives 6.79Γ103Β W and boils off 123Β kgβhβ1.
A long solid cylinder has radius 0.010Β m, length
1.0Β m and thermal conductivity
5.0Β Wβmβ1βKβ1. Uniform heat generation occurs
at ππ=2.0Γ106Β Wβmβ3. Its curved surface is
cooled by fluid at 300Β K, with a heat transfer coefficient
β=200Β Wβmβ2βKβ1.
Assume steady radial conduction, constant properties and negligible heat
loss through the ends. Derive the temperature profile and calculate the
surface and centre temperatures. Check the overall heat balance.
Compare two modifications: doubling β, and replacing the cylinder by
four cylinders of half the radius. Each smaller cylinder has the same
length, conductivity, heat generation per unit volume and original heat
transfer coefficient. The coolant temperature is unchanged.
Worked solution
The radial geometry is the same as the wire model in Figure 9. Here the surface temperature must be found from convection before the centre temperature can be evaluated.
Data
Quantity
Symbol
Value
Cylinder radius and length
π ,πΏ
0.010,1.0Β m
Thermal conductivity
π
5.0Β Wβmβ1βKβ1
Volumetric heat generation
ππ
2.0Γ106Β Wβmβ3
Coolant temperature
πβ
300Β K
Surface coefficient
β
200Β Wβmβ2βKβ1
Assumptions
Steady radial conduction with uniform heat generation.
Constant conductivity and convection coefficient.
Negligible heat flow through the ends.
Each comparison cylinder sees the same specified coolant temperature.
Calculations
The governing equation is the same as for the heated wire. At the axis,
symmetry requires a zero radial temperature gradient:
The surface temperature is now unknown, so we need a convection boundary
condition at π=π . At steady state, the conductive flux reaching the
surface equals the convective flux leaving it:
The two temperature rises come from different parts of the model. The
surface rise follows from the overall generation-to-convection balance:
ππ βπβ=πππ 2β=2.0Γ106(0.010)2(200)=50.0Β K
The conduction profile then gives the additional centre rise:
ππππ₯βππ =πππ 24π=2.0Γ106(0.010)24(5.0)=10.0Β K
The surface-to-coolant temperature difference is 50.0Β K.
The centre-to-surface difference is a further 10.0Β K.
The temperatures are therefore ππ =350.0Β K and
πmaxβ‘=π(0)=360.0Β K.
The total heat generated is ππππ 2πΏ=628.3Β W. The surface
heat loss is
β(2ππ πΏ)(ππ βπβ)=628.3Β W,
so the
overall energy balance is satisfied. This check also confirms the surface
boundary condition used to obtain the temperature profile.
Doubling β halves the surface-to-coolant difference to
25.0Β K. The centre-to-surface difference stays at
10.0Β K because ππ, π and π are unchanged. The
new centre temperature is
300+25+10=335.0Β K.
Four cylinders of half the radius have the same total volume and heat
generation as the original cylinder, but twice its total curved surface
area. In each smaller cylinder, the surface-to-coolant difference is
25.0Β K and the centre-to-surface difference is
2.50Β K. The centre temperature is therefore
300+25+2.50=327.5Β K.
The smaller radius reduces both temperature differences. At fixed ππ,
β and π, the surface-to-coolant difference varies as π , whereas
the centre-to-surface difference varies as π 2. These results assume
that every smaller cylinder is exposed to the same coolant temperature
and specified heat transfer coefficient.
Answer
The original cylinder has ππ =350.0Β K, ππππ₯=360.0Β K and total heat loss 628Β W. Doubling β lowers the centre to 335.0Β K. Four cylinders of half the radius lower each centre to 327.5Β K while preserving the original total heat generation.
A hot spherical particle made out of carbon steel (1.25 cm diameter,
500Β Β°C, π=40Β W/(mβ K)) is cooled by a
fan, with convection coefficient of 110Β W/(m2β K) and
air temperature of 25Β Β°C. How long would it take for the
sphere to reach a temperature of 100Β Β°C?
Worked solutionFigure 1: The steel sphere cools in air. Its radius is half its diameter. After checking the Biot number, the lumped model treats its centre and surface as having the same temperature.
Data
Quantity
Value
Diameter
π·=0.0125Β m
Initial temperature
ππ=500Β Β°C
Target temperature
π=100Β Β°C
Air temperature
πβ=25Β Β°C
Heat transfer coefficient
β=110Β Wβmβ2βKβ1
Steel conductivity
π=40Β Wβmβ1βKβ1
Steel density
π=8131Β kgβmβ3
Steel specific heat capacity
πΆπ=434Β Jβkgβ1βKβ1
The original slide does not supply density or specific heat capacity. These two values come from the lecturerβs notes.
Assumptions
Uniform initial temperature.
Constant material properties and heat transfer coefficient.
Constant air temperature.
Radiation and heat transfer through supports are neglected.
No internal heat generation.
Calculations
First test whether conduction within the sphere is fast enough to keep its
internal temperature nearly uniform. For a sphere,
πΏπ=ππ΄π =π·6=0.01256=0.0020833Β m
The Biot number compares internal conduction resistance with the surface
convection resistance:
Biπ=βπΏππ=110(0.0020833)40=0.00573
Since Biπ<0.1, the lumped model is appropriate. The energy balance
is
π=πππΆπβπ΄π =ππΆππΏπβ=8131(434)(0.0020833)110=66.84Β s
Integrating the balance, as in Eq. 3, gives the time to the
requested temperature:
π‘=πlnβ‘ππβπβπβπβ
Substituting the temperatures and time constant gives
π‘=66.84lnβ‘500β25100β25=123.4Β s(9)
The required excess-temperature ratio is 75/475=0.1579, between
πβ1 and πβ2. A time between one and two time constants is therefore
consistent with the result. Radiation would increase heat loss, particularly
at the initial high temperature, so the result applies to convection alone.
Answer
The sphere reaches 100βΒ°C after approximately
123Β s, or 2.06Β min, in the convection-only model.
In a manufacturing facility, 5 cm diameter brass spheres
(π=111Β W/(mβ K), π=8500Β kg/m3,
πΆπ=385.2Β J/(kgβ K)) initially at 121Β Β°C
are quenched in a water bath at 20Β Β°C for a period of
2 minutes at a rate of 120 spheres per minute. If the convection heat
transfer coefficient 240Β W/(m2β K), determine:
The temperature of the spheres after quenching
The rate at which heat needs to be removed from the water in order to keep its temperature constant at 20Β Β°C.
Worked solutionFigure 2: Radial section of a sphere in a fluid bath, showing radius, centre temperature, surface temperature and the external convection condition. In the lumped model the centre and surface temperatures are approximately equal.
Data
Quantity
Value
Sphere diameter
π·=0.050Β m
Brass conductivity
π=111Β Wβmβ1βKβ1
Brass density
π=8500Β kgβmβ3
Specific heat capacity
πΆπ=385.2Β Jβkgβ1βKβ1
Initial temperature
ππ=121Β Β°C
Bath temperature
πβ=20Β Β°C
Residence time
π‘=120Β s
Heat transfer coefficient
β=240Β Wβmβ2βKβ1
Sphere throughput
Λπ=120/60=2Β sβ1
Assumptions
Uniform initial temperature in each sphere.
Constant material properties and heat transfer coefficient.
Uniform bath temperature, maintained at 20βΒ°C.
No internal heat generation.
Heat transfer from the bath to other surroundings is neglected.
Steady throughput for the bath heat-duty calculation.
Calculations
The characteristic length for a sphere is
πΏπ=π·6=0.0506=0.0083333Β m
Test the internal temperature gradients using
Biπ=βπΏππ=240(0.0083333)111=0.01802
This is below 0.1, so each sphere can be treated as a lumped thermal
capacitance. Its time constant is
π=ππΆππΏπβ=8500(385.2)(0.0083333)240=113.69Β s
At steady throughput, the cooling equipment must remove heat at the same
average rate as the spheres release it to the bath:
Λπbath=ΛππΈ1=12060(14.11)=28.22Β kW(10)
The exit temperature lies between the bath temperature and the initial
sphere temperature, as required. There are 240 spheres in the bath at
steady operation, but heat duty depends on the throughput of two spheres
per second. Inventory determines how many spheres are cooling at once.
Answer
The spheres leave at approximately 55.2βΒ°C. The required
heat-removal rate is 28.2Β kW.
A flat wall of fire-clay brick (0.5 m thick) initially has a uniform
temperature of 200 K. One of its faces is suddenly exposed to hot gas at
1200 K. If the heat transfer coefficient on the hot side is
7.38Β W/(m2β K) and the other face of the wall is insulated,
determine:
The time it takes for the centre of the wall to reach 600 K
The temperature of the insulated wall at the time evaluated in (a)
Properties of fire-clay brick:
π=1.125Β W/(mβK)
πΆπ=919Β J/(kgβK)
π=2310Β kg/m3
Worked solution
The insulated face has zero normal temperature gradient. Reflecting the
brick about that face produces a symmetric wall heated from both sides.
Figure 5: A brick wall heated by gas on one face and insulated on the other. Reflection makes the insulated face the symmetry plane of a wall twice as thick.
Data
Quantity
Value
Real brick thickness
π=0.50Β m
Initial temperature
ππ=200Β K
Gas temperature
πβ=1200Β K
Heat transfer coefficient
β=7.38Β Wβmβ2βKβ1
Brick conductivity
π=1.125Β Wβmβ1βKβ1
Brick density
π=2310Β kgβmβ3
Specific heat capacity
πΆπ=919Β Jβkgβ1βKβ1
Reflected-wall half-thickness
πΏ=0.50Β m
Physical brick-centre position
π₯=0.25Β m, π=π₯/πΏ=0.5
Target temperature at the brick centre
π=600Β K
Measure position from the insulated face. Reflection makes that face the midplane of the equivalent wall.
Assumptions
Uniform initial temperature.
Unidirectional conduction normal to the wall.
Perfect insulation at the rear face.
Constant material properties, gas temperature and heat transfer coefficient.
No internal heat generation.
Radiation is neglected.
Calculations
The insulated face obeys βπ/βπ₯=0, which is the symmetry
condition of the finite-wall solution. First calculate the diffusivity:
For this one-sided wall, π/π΄π =0.50Β m also. The Biot number is
well above 0.1, so a spatial temperature profile is needed. The one-term
wall solution from Eq. 18 is
At the real brick centre, the requested temperature corresponds to
Ξ(0.5,π‘)=600β1200200β1200=0.600
The value 0.600 applies at the physical centre of the brick,
π₯=0.25Β m or π=0.5. The reference position for the symmetric
wall solution is π₯=0, at the insulated face. This is the centre of the
equivalent reflected wall, so it is a different location from the physical
centre of the brick.
At π=0, the cosine factor is one and the one-term solution gives
Ξ0=Ξ(0,π‘)=π΄1πβπ21Fo. At the physical
brick centre, it therefore gives
Ξ(0.5,π‘)=Ξ0cosβ‘(π1/2). Divide by this cosine factor
to find the dimensionless temperature at the insulated face:
Thus 0.600 and 0.73096 are dimensionless temperatures at different
positions at the same time. The value 0.73096 is the insulated-face
reference value needed to find the elapsed time.
Since Fo>0.2, the one-term approximation is appropriate. Convert
this dimensionless time to seconds:
π‘=FoπΏ2πΌ=0.34452(0.50)25.2992Γ10β7=1.6253Γ105Β s=45.1Β h
At the insulated face, π=0 and the cosine factor is one. Its temperature is
π0=πβ+Ξ0(ππβπβ)=1200+0.73096(200β1200)=469.0Β K
The insulated face is colder than the brick centre because heat enters from
the opposite face. The ordering 200<π0<600<1200Β K is consistent
with one-sided heating. Constant properties over this temperature range are
an idealisation of the classroom model.
Answer
The brick centre reaches 600Β K after approximately
45.1Β h. At that time, the insulated face is at
469Β K.
A steel cylinder (π·=30.5Β cm, πΏ=61Β cm) is initially
at 645 K. It is immersed in an oil bath with a constant temperature of 310K.
If the convection heat transfer coefficient at the steel surface is 34Β W/(m2β K),
determine the temperatures at the following locations:
at the centre of the cylinder after 1 hour?
at the centre of the top surface of the cylinder?
Properties of steel:
π=42.9Β W/(mβK)
πΆπ=473Β J/(kgβK)
π=7820Β kg/m3
Worked solution
An axial section identifies the two requested points and the two directions
of conduction.
Data
Quantity
Value
Cylinder diameter
π·=0.305Β m
Full cylinder length
πΏ=0.610Β m
Half-length for the plane-wall factor
πΏπ€=0.305Β m
Cylinder radius
ππ=0.1525Β m
Initial temperature
ππ=645Β K
Oil temperature
πβ=310Β K
Heat transfer coefficient
β=34Β Wβmβ2βKβ1
Steel conductivity
π=42.9Β Wβmβ1βKβ1
Steel density
π=7820Β kgβmβ3
Specific heat capacity
πΆπ=473Β Jβkgβ1βKβ1
Time for both requested temperatures
π‘=3600Β s
Assumptions
Uniform initial temperature.
Axisymmetric conduction in both radial and axial directions.
Constant material properties.
Constant oil temperature and uniform heat transfer coefficient on all surfaces.
No internal heat generation.
Radiation is neglected.
Calculations
The ends cannot be neglected because the length is only twice the diameter.
For a constant-property body with these separable boundary conditions, the
short-cylinder excess-temperature ratio is the product of a wall factor
and a long-cylinder factor:
ππ=πβ+Ξπ(ππβπβ)=310+0.56437(645β310)=499.1Β K
The face centre is cooler than the body centre and both temperatures remain
between the initial and bath temperatures. A position chart supplies a
ratio relative to its own centre value. Multiplying a position ratio alone
by the other centre factor would omit part of the cooling history.
Answer
After one hour, the body centre is at approximately 522Β K and
the centre of the top face is at approximately 499Β K.
An egg, initially at a uniform temperature of 5Β Β°C, is
dropped into boiling water at 95Β Β°C. Taking the convection
heat transfer coefficient to be 1200Β W/(m2β Β°C),
and approximating the egg to be a 5-cm diameter sphere, determine how long
it will take for the centre of the egg to reach 70Β Β°C.
Worked solutionFigure 6: The egg is modelled as a sphere heated by surrounding water. The requested temperature is at the centre, which can lag behind the surface.
Data
Quantity
Value
Model sphere diameter
π·=0.050Β m
Initial temperature
ππ=5Β Β°C
Target centre temperature
π0=70Β Β°C
Water temperature
πβ=95Β Β°C
Heat transfer coefficient
β=1200Β Wβmβ2βKβ1
Conductivity
π=0.627Β Wβmβ1βKβ1
Thermal diffusivity
πΌ=0.151Γ10β6Β m2βsβ1
The conductivity and diffusivity come from the lecturerβs notes.
Assumptions
Homogeneous spherical body with a uniform initial temperature.
Radial conduction.
Constant properties and heat transfer coefficient.
Constant water temperature.
No internal heat generation.
Shell resistance and energy changes due to cooking are neglected.
Calculations
Test lumped capacitance using πΏπ=π·/6:
Biπ=βπ·6π=1200(0.050)6(0.627)=15.95
The large Biot number rules out a uniform internal temperature. The sphere
solution uses the radius ππ=π·/2=0.025Β m instead:
This just exceeds 0.2, so the usual one-term criterion is satisfied.
The time is
π‘=Foπ2ππΌ=0.20842(0.025)20.151Γ10β6=863Β s=14.4Β min
A chart reading of Foβ0.20 gives about 13.8 minutes,
consistent with this estimate. The large Biot number means that internal
conduction limits the response. Increasing the external coefficient has
less effect than it would in a lumped body. Real eggs have a nonspherical
shape and different internal regions, so the calculation is an idealised
heating estimate.
Answer
The centre reaches 70βΒ°C after approximately
14Β min in the spherical, constant-property model.
What minimum burial depth would you recommend to avoid freezing under
conditions for which the soil, initially at a uniform temperature of
15Β Β°C, is subjected to a constant surface temperature of
β10Β Β°C for 3 months?
Properties of soil:
π=0.52Β W/(mβ K)
πΆπ=1840Β J/(kgβ K)
π=2050Β kg/m3
Worked solution
Measure depth downwards from the cold ground surface. The pipe must lie
below the calculated zero-degree isotherm.
Data
Quantity
Value
Initial soil temperature
ππ=15Β Β°C
Prescribed surface temperature
ππ =β10Β Β°C
Temperature at the required depth
π=0Β Β°C
Elapsed time
π‘=90Β days
Soil conductivity
π=0.52Β Wβmβ1βKβ1
Soil density
π=2050Β kgβmβ3
Specific heat capacity
πΆπ=1840Β Jβkgβ1βKβ1
For this calculation, take three months as 90 days.
Assumptions
Uniform initial soil temperature.
Unidirectional conduction normal to a plane ground surface.
Semi-infinite soil depth over the time considered.
Constant thermal properties.
Constant prescribed surface temperature.
No internal heat generation.
Soil moisture phase change is neglected.
Calculations
The surface temperature is supplied directly, so use the prescribed-temperature
semi-infinite solution rather than a convection boundary condition. The thermal
diffusivity is
Inverting the complementary error function gives π=0.370807.
The corresponding depth is
π₯=2πβπΌπ‘=2(0.370807)(1.0354)=0.7679Β m
Temperature increases with depth during cooling, so a pipe deeper than
this isotherm remains above 0βΒ°C within the model.
The calculated depth is consistent with the metre-scale diffusion distance.
The far boundary must be well below the disturbed region for the
semi-infinite approximation to hold.
Soil moisture releases latent heat on freezing, which this calculation
omits. Snow cover and a changing surface temperature would also alter the
temperature profile. The result is a classroom estimate based on the
specified temperature history.
Answer
The 0βΒ°C isotherm reaches approximately
0.77Β m. A burial depth of at least 0.8Β m lies below
it in this model; 1Β m provides a larger margin.
A large slab of full thickness 0.100 m starts uniformly at 200 degrees
Celsius. Both broad faces cool in fluid at 20 degrees Celsius with
β=20 W mβ»Β² Kβ»ΒΉ. Use π=1.0 W mβ»ΒΉ Kβ»ΒΉ, π=2000 kg mβ»Β³ and
πΆπ=1000 J kgβ»ΒΉ Kβ»ΒΉ. Neglect edges, radiation and property variation.
After 30 minutes, find the centre, surface and volume-average temperatures,
and the energy released per square metre of slab plan area. For the plane-wall
one-term solution at Bi=1, use π1=0.860334 and
π΄1=1.11913. Check that this approximation is appropriate.
Worked solutionFigure 8: Both slab faces cool under the same conditions. The solution uses half-thickness L, measured from the symmetry plane to either exposed surface.
Data
Quantity
Value
Full slab thickness
2πΏ=0.100Β m
Half-thickness
πΏ=0.050Β m
Initial temperature
ππ=200Β Β°C
Fluid temperature
πβ=20Β Β°C
Cooling time
π‘=1800Β s
Heat transfer coefficient
β=20Β Wβmβ2βKβ1
Conductivity
π=1.0Β Wβmβ1βKβ1
Density
π=2000Β kgβmβ3
Specific heat capacity
πΆπ=1000Β Jβkgβ1βKβ1
Supplied first eigenvalue
π1=0.860334
Supplied first coefficient
π΄1=1.11913
Assumptions
Uniform initial temperature.
Unidirectional conduction through the slab thickness.
Identical convection conditions on both broad faces.
Constant properties, fluid temperature and heat transfer coefficient.
Edge heat transfer and radiation are neglected.
No internal heat generation.
Calculations
Both broad faces transfer heat. If π΄plan is the area of one
face, π=2πΏπ΄plan and π΄π =2π΄plan. Therefore
πΏπ=ππ΄π =πΏ=0.050Β m.
The Biot number is
Biπ=βπΏππ=20(0.050)1.0=1.0
Lumped capacitance is unsuitable. Use the spatial wall solution, with
π₯=0 at the midplane. First calculate the diffusivity:
This includes heat leaving both faces. The ordering
20<ππ <Β―Β―Β―Β―π<ππ<200βΒ°C is consistent with cooling
from the surfaces. Using the centre temperature for the energy balance
would underestimate the heat released, because most of the slab is cooler
than its centre.
Answer
After 30 minutes, ππ=174.3βΒ°C,
ππ =120.6βΒ°C and Β―Β―Β―Β―π=156.0βΒ°C.
The energy released is 8.80Β MJβmβ2 of plan area.
Air at a pressure of 6Β kN/m2 and temperature of
300Β Β°C flows with a velocity of 10Β m/s over
a flat plate 0.5 m long. Estimate the cooling rate per unit width of the
plate needed to maintain the surface at a temperature of 27Β Β°C.
π@1Β atm=30.84Γ10β6Β m2/s
πairΒ atΒ 165βΒ°C=3.6427Γ10β2Β W/(mβ K)
PrairΒ atΒ 165βΒ°Cβ0.684
SolutionFigure 2: Air flows along the plate, and the velocity and thermal boundary layers grow from its leading edge. Use the full plate length for the average coefficient.
The plate geometry is the boundary-layer geometry in Figure 1. The stream runs along the 0.50 m length, and the cooling rate is requested per metre of width.
Data
Air temperature πβ=300Β Β°C
Plate temperature ππ =27Β Β°C
Plate length πΏ=0.50Β m
Air speed π’β=10Β mβsβ1
Absolute pressure π=6Β kPa
The supplied conductivity, Prandtl number and one-atmosphere viscosity are evaluated near the film temperature.
Assumptions
Steady state.
Smooth plate aligned with the free stream.
Constant fluid properties at the film temperature.
Ideal-gas density and pressure-independent dynamic viscosity.
Radiation and buoyancy effects are neglected.
Calculations
The hot air heats the colder plate. The required cooling rate is the heat
received by one exposed plate face. First evaluate the film temperature:
ππ=ππ +πβ2=27+3002=163.5Β Β°C
The supplied properties at 165Β Β°C approximate this temperature.
At fixed temperature, ideal-gas density is proportional to absolute pressure.
Since π=π/π and the dynamic viscosity is nearly pressure independent,
the kinematic viscosity increases as pressure falls:
As a geometry check, the laminar thickness estimate gives
πΏ(πΏ)β5πΏβReπΏ=0.0255Β m,
about five percent of the plate length. Ignoring the low-pressure viscosity
correction would overestimate β by a factor of about 4.11. The estimated
molecular mean free path is much smaller than this boundary-layer thickness,
so the low pressure does not by itself invalidate the continuum model.
Answer
The cooling system must remove approximately
570Β W per metre of plate width from the exposed face.
A metallic cylinder (12.7 mm in diameter and 94 mm long) is heated
internally by an electrical heater and is subjected to a cross flow of air.
The air velocity and temperature are maintained at 10Β m/s and
26.2Β Β°C. The heater power dissipation is measured to be
39.1 W, while the average cylinder surface temperature is 128.4Β Β°C.
Compare the heat transfer coefficient from experimental observations and
correlation predictions.
SolutionFigure 4: Air flows across the cylinder, perpendicular to its axis. The convecting lateral area is pi D L when the end faces are neglected.
Data
Quantity
Value
Cylinder diameter
π·=0.0127Β m
Cylinder length
πΏ=0.094Β m
Air speed
π’β=10Β mβsβ1
Air temperature
πβ=26.2Β Β°C
Surface temperature
ππ =128.4Β Β°C
Measured heat dissipation
Λπ=39.1Β W
Air density
π=0.995Β kgβmβ3
Dynamic viscosity
π=2.082Γ10β5Β Paβs
Air conductivity
π=0.030Β Wβmβ1βKβ1
Prandtl number
Pr=0.700
The film temperature is the mean of the surface and air temperatures, 77.3 degrees Celsius or 350.45 K. Use the lecturerβs property set near 350 K.
Assumptions
Steady state.
Uniform mean surface temperature.
Constant fluid properties at the film temperature.
The stated 39.1Β W leaves through the curved cylinder surface.
Radiation and heat transfer through the ends are neglected in the comparison.
The long-cylinder cross-flow correlations approximate the finite cylinder.
Calculations
The measured coefficient follows from Newtonβs relation. The curved area is
The power law is 13.6 percent below the measured coefficient; Churchill
and Bernstein is 5.9 percent below it. The experimental value is an effective
coefficient under the assumed heat balance. The original slide diagram also
shows 46Β W heater input and 6.9Β W losses,
whose difference is 39.1Β W. Using the full 46Β W
with the curved area would include heat leaving by another route.
Radiation would reduce the heat assigned to convection. For illustration,
π=0.5 and surroundings at the air temperature give
Λπrad=πππ΄(π4π βπ4surr)=0.5(5.6704Γ10β8)(0.00375043)[(401.55)4β(299.35)4]=1.91Β W
Subtracting this illustrative loss gives a convection-only measured
coefficient of about 97.0Β Wβmβ2βKβ1.
The emissivity is not supplied, so this is a sensitivity check. The finite
aspect ratio πΏ/π·=7.40 and uncertain support losses also limit the comparison.
Answer
The measured effective coefficient is approximately
102Β Wβmβ2βKβ1. The cylinder power law predicts
88.2Β Wβmβ2βKβ1 and Churchill-Bernstein predicts
96.0Β Wβmβ2βKβ1 under the stated assumptions.
Air flows along a smooth flat plate at 15 m sβ»ΒΉ. The measured average
skin-friction coefficient is Β―Β―Β―Β―πΆπ=0.0040. At the film temperature,
use π=1.18 kg mβ»Β³, πΆπ=1005 J kgβ»ΒΉ Kβ»ΒΉ and Pr=0.72.
Estimate the average heat transfer coefficient with the Colburn analogy.
Explain why the uncorrected Reynolds analogy gives a different answer.
Assume matching flow and thermal conditions for the analogy and negligible
pressure drag.
Solution
Data
Air speed π’β=15Β mβsβ1
Average skin-friction coefficient Β―Β―Β―Β―πΆπ=0.0040
Film density π=1.18Β kgβmβ3
Specific heat capacity πΆπ=1005Β Jβkgβ1βKβ1
Prandtl number Pr=0.72
Assumptions
Matching flow and thermal conditions for the Colburn analogy.
Negligible pressure drag.
Constant fluid properties.
Calculations
Use average quantities throughout. Equation Eq. 13 gives
The plain Reynolds analogy would set Β―Β―Β―Β―Β―Β―St=0.0020
and give 35.6 W mβ»Β² Kβ»ΒΉ. It assumes Pr equal to one. At the supplied Pr,
the Colburn correction increases the estimate by a factor of 1.245.
No viscosity or plate length is required because skin friction is measured.
Using a local friction coefficient with an average heat transfer coefficient
would mix two different quantities. Total drag on a sphere or cylinder also
cannot replace the plate skin-friction coefficient in this calculation.
Answer
The Colburn estimate is Β―Β―Β―Β―β=44.3Β Wβmβ2βKβ1. The uncorrected Reynolds analogy gives 35.6Β Wβmβ2βKβ1, because it omits the Prandtl-number correction.
Air flows at 0.80 m sβ»ΒΉ over an isothermal flat plate. In the laminar region,
the local coefficient at π₯1=0.25 m is βπ₯1=12.0 W mβ»Β² Kβ»ΒΉ.
Use π=1.60Γ10β5 mΒ² sβ»ΒΉ and constant properties. At what distance
will the local coefficient be 6.0 W mβ»Β² Kβ»ΒΉ? Find the average coefficient
between the leading edge and that distance, and check the laminar assumption.
Solution
Data
Air speed π’β=0.80Β mβsβ1
First position π₯1=0.25Β m
Coefficient at the first position βπ₯1=12.0Β Wβmβ2βKβ1
For fixed speed and properties, the local relation has
Nuπ₯βRe1/2π₯βπ₯1/2.
Because βπ₯=Nuπ₯π/π₯, it follows that βπ₯βπ₯β1/2:
βπ₯2βπ₯1=(π₯1π₯2)1/2π₯2=π₯1(12.06.0)2=1.00Β m
Halving a local coefficient requires four times the distance, not twice.
The layer thickens downstream and the wall temperature gradient decreases.
The length average is a different quantity. Integrating the local law gives
Thus the average over the first 1.00 m is 12.0 W mβ»Β² Kβ»ΒΉ even though the
local coefficient at its trailing edge is 6.0 W mβ»Β² Kβ»ΒΉ. The ideal leading-edge
singularity is integrable.
At π₯2,
Reπ₯2=0.80(1.00)π=50000,
below the
transition range discussed above. The whole calculation stays in the
laminar region. A measured halving distance beyond transition would need
the appropriate changing flow model.
Answer
The local coefficient falls to half at π₯2=1.00Β m. The average coefficient over that length is 12.0Β Wβmβ2βKβ1. The Reynolds number is 50000, consistent with laminar flow.
Air (1 atm, 77Β Β°C) enters a 5 mm I.D. tube with a bulk
average velocity of 2.5Β m/s. The velocity profile is developed
and the thermal profile is βdevelopingβ. The tube length is 1.0 m, and a
constant (uniform) heat flux is imposed at the tube surface on the air over
the entire length. The bulk average exit temperature of the air stream is
127Β Β°C. Determine:
The heat transfer coefficient at the exit (βπΏ)
The heat flux imposed at the surface of the tube (ππ )
The surface temperature (ππ ) of the tube at the exit
Use ideal-gas air with π air=287.05Β J/(kgβK) and
π=101325Β Pa. At the mean bulk temperature of 375.15 K, use
π=2.1805Γ10β5Β Paβs,
π=0.03184Β W/(mβK) and πΆπ=1011.5Β J/(kgβK).
SolutionFigure 3: Heat enters through the tube wall while air flows from inlet Ti to outlet Te. In this example the inward wall flux is uniform.
Data
Quantity
Value
Inside diameter
π·=0.0050Β m
Heated length
πΏ=1.0Β m
Inlet mean speed
π’π=2.5Β mβsβ1
Inlet bulk temperature
ππ=77Β Β°C
Outlet bulk temperature
ππ=127Β Β°C
Inlet air density
ππ=1.00810Β kgβmβ3
Dynamic viscosity
π=2.1805Γ10β5Β Paβs
Air conductivity
π=0.03184Β Wβmβ1βKβ1
Specific heat capacity
πΆπ=1011.5Β Jβkgβ1βKβ1
The wall heat flux is uniform, and the inlet velocity profile is developed. The mean bulk temperature is 102 degrees Celsius or 375.15 K. Use transport properties near 375 K, but pair inlet velocity with inlet
density:
Constant viscosity, conductivity and heat capacity; ideal-gas density varies with temperature at approximately constant pressure.
Heat supplied at the wall enters the air.
Axial conduction, viscous heating and kinetic and potential energy changes are neglected.
The tube is circular and the flow remains laminar.
Calculations
First check the flow regime and whether the temperature profile has developed
by the exit. Using the diameter in the Reynolds number, as defined in Eq. 18, gives
The flow is laminar. The thermal-entry estimate from Eq. 19 is
πΏπ‘β0.0575Reπ·Prπ·=0.0575(577.91)(0.69271)(0.0050)=0.1151Β m
The temperature profile develops after heating starts, but the exit at
1.0Β m lies well beyond this estimated entrance length.
Use the developed local constant-flux result from Eq. 12, Nuπ·=48/11:
The wall is 5.73Β K hotter than the exit air, as required for
inward heat flow. Near the inlet the coefficient is larger, so the same
heat flux needs a smaller wall-to-bulk difference. The developed coefficient
must not be applied over that entrance region. The large axial Peclet product
ReβPrπΏπ·β8.01Γ104
also supports neglecting
axial conduction. At the same volume flow and temperature rise, replacing
air by water would require a much larger heat input because water has a much
larger volumetric heat capacity.
Answer
At the exit, βπΏβ27.8Β Wβmβ2βKβ1.
The imposed heat flux is ππ β159Β Wβmβ2, and the
exit wall temperature is ππ (πΏ)β132.7Β Β°C.
A horizontal tubular condenser is used to condense saturated steam at
80Β Β°C. The condenser is a shell and tube one with brass
tubes (π=110Β W/(mβ K), 1.59 cm OD, 1.34 cm ID). Steam
is outside of the tubes (β=10760Β W/(m2β K)) and cooling
water enters the tube at 20Β Β°C with a velocity of
1.4Β m/s. Water leaves the tubes at 40Β Β°C.
If the rate of cooling water supply is 55000 kg/h and the latent heat of
condensation of steam at 80Β Β°C is 2304 kJ/kg, calculate:
The number of tubes required
The length of each tube
Use inlet water density ππ=998.2Β kg/m3 at
20 degrees Celsius. At the bulk mean of 30 degrees Celsius, use
π=7.98Γ10β4Β Paβs,
π=0.615Β W/(mβK) and πΆπ=4180Β J/(kgβK).
Treat 1.4 m/s as the nominal inlet speed; after rounding the number of tubes
up, use the exact total mass flow to calculate the duty and Reynolds number
of each tube.
SolutionFigure 4: Water is heated as it passes through a condenser tube. The heat flux is inward; the wall temperature and local flux need not be uniform. The steam film and brass wall are included in the overall coefficient.
Data
Quantity
Value
Steam temperature
πsteam=80Β Β°C
Outer diameter
π·π=0.0159Β m
Inner diameter
π·π=0.0134Β m
Brass conductivity
ππ€=110Β Wβmβ1βKβ1
Outside coefficient
βπ=10760Β Wβmβ2βKβ1
Water inlet temperature
ππ=20Β Β°C
Water outlet temperature
ππ=40Β Β°C
Nominal inlet water speed
π’π=1.4Β mβsβ1
Total water mass flow
Λπ=55000Β kgβhβ1
Inlet water density
ππ=998.2Β kgβmβ3
Dynamic viscosity
π=7.98Γ10β4Β Paβs
Water conductivity
π=0.615Β Wβmβ1βKβ1
Specific heat capacity
πΆπ=4180Β Jβkgβ1βKβ1
The rounded water properties are consistent with the inlet and bulk mean
states. See the IAPWS liquid-water property release.
The inlet density determines flow area; mean-temperature viscosity,
conductivity and heat capacity determine heat transfer.
Assumptions
Steady state.
Equal flow through parallel tubes.
Constant properties and heat transfer coefficients.
Developed turbulent flow in smooth tubes.
Clean walls, with fouling resistance neglected.
Negligible heat loss to the surroundings.
Calculations
The nominal inlet speed first determines a trial mass flow per tube:
Heat passes through the steam film, the brass wall and the water film in
series. The steam temperature is fixed, but the inner wall is cooler.
Using the cylindrical conduction resistance, the overall coefficient on
inside area satisfies
Thus ππ=4113Β Wβmβ2βKβ1.
The water film supplies about 63.5 percent of the total resistance, the steam
film 32.2 percent and the brass wall 4.3 percent. Changing the wall conductivity
alone would therefore have little effect. Area bases must remain consistent,
with πππ΄π=πππ΄π.
The driving differences are 80β20=60Β K at the water inlet
and 80β40=40Β K at the outlet. Their logarithmic mean, defined in Eq. 15, is
Ξπlnβ‘=ΞππβΞππlnβ‘(ΞππΞππ)=60β40lnβ‘(6040)=49.326Β K
The per-tube heat duty follows from the water enthalpy rise:
Λπtube=ΛπtubeπΆπ(ππβππ)=0.195869(4180)(40β20)=16375Β W
Set this equal to ππππ·ππΏΞπlnβ‘ and solve for length:
The ratio πΏ/π·πβ143 supports the developed turbulent approximation.
The ideal isothermal steam source needs no flow-arrangement correction to
the logarithmic mean. The arithmetic mean 50Β K would overestimate
the driving force by about 1.37 percent.
As a steam-side check, the total water duty and the stated latent heat give
The latent heat checks the steam requirement; it does not determine tube
count. These are clean-tube estimates using the stated correlation. A
practical design also needs allowances for fouling, pressure drop and
available tube sizes.
Answer
Using properties at the stated water temperatures, provide 78 parallel tubes
of approximately 1.92Β m heated length each.
Oil enters a 16 mm diameter pipe with Re=900 and
Pr=60. Estimate the hydrodynamic and thermal entrance lengths.
Is a 5 m pipe long enough for both profiles to develop? Assume a uniform
velocity and temperature at the inlet and heating beginning there.
Solution
Data
Diameter π·=0.016Β m
Pipe length πΏ=5.0Β m
Reynolds number Re=900
Prandtl number Pr=60
Assumptions
Uniform inlet velocity and temperature.
Heating begins at the pipe inlet.
Constant properties.
Calculations
The Reynolds number is below 2300, so use the laminar estimates in
Eq. 19, converting the diameter to 0.016 m:
πΏβ=0.0575(900)(0.016)=0.828Β mπΏπ‘=0.0575(900)(60)(0.016)=49.7Β m
The velocity profile develops within the pipe, but the temperature profile
does not. The ratio πΏπ‘/πΏβ=Pr=60 expresses the slower thermal
diffusion relative to momentum diffusion in this oil. A heat transfer
calculation for the 5 m pipe needs a thermal-entry correlation. A developed
velocity profile alone does not justify a developed Nusselt number.
If heating starts after a long unheated approach, measure thermal entry from
the start of heating. The two entrance lengths then have different origins.
Answer
The hydrodynamic entrance length is πΏβ=0.828Β m and the thermal entrance length is πΏπ‘=49.7Β m. The 5.0Β m pipe permits velocity development but remains thermally developing.
Water enters a 2.0Β m heated tube at 20Β Β°C.
The inside diameter is 0.020Β m, mass flow is
0.020Β kgβsβ1 and πΆπ=4180Β Jβkgβ1βKβ1.
Adopt a constant effective coefficient β=400Β Wβmβ2βKβ1 to isolate the effect of the wall
condition. Neglect axial conduction and heat loss.
Find the outlet temperature and total heat input for two separate cases:
a wall maintained at 80Β Β°C, and a uniform inward wall heat flux
of 10000Β Wβmβ2. Explain how the wall-to-bulk temperature difference changes
with distance in each case.
SolutionFigure 5: The same flowing stream is heated in two separate cases. One fixes the wall temperature; the other fixes inward heat flux. The energy balance relates inlet and outlet bulk temperatures in both cases.
Data
Quantity
Value
Tube length
πΏ=2.0Β m
Diameter
π·=0.020Β m
Mass flow
Λπ=0.020Β kgβsβ1
Inlet temperature
ππ=20Β Β°C
Specific heat capacity
πΆπ=4180Β Jβkgβ1βKβ1
Supplied effective coefficient
β=400Β Wβmβ2βKβ1
Wall temperature in the first case
ππ =80Β Β°C
Wall heat flux in the second case
ππ =10000Β Wβmβ2
Assumptions
Steady state.
Constant fluid properties.
Axial conduction and heat loss are neglected.
Calculations
The heated area is π΄π =ππ·πΏ=0.12566 mΒ² and the stream heat-capacity
rate is ΛππΆπ=83.6 W Kβ»ΒΉ.
βπ΄π ΛππΆπ=400[π(0.020)(2.0)]0.020(4180)=0.6013ππ=ππ β(ππ βππ)πββπ΄π ΛππΆπ=80β(80β20)πβ0.6013=47.11Β Β°CΛπ=83.6(47.11β20)=2267Β W
The driving difference decreases from 60 K to 32.89 K. Even with a constant
coefficient, the heat flux falls from 24000 to 13155 W mβ»Β². Fully developed
flow does not require a constant temperature difference.
For constant wall heat flux, the overall balance instead gives
Λπ=ππ π΄π =10000[π(0.020)(2.0)]=1257Β Wππ=ππ+ΛπΛππΆπ=20+125783.6=35.03Β Β°Cππ βππ=ππ β=25.0Β K
The bulk temperature rises linearly, and the wall temperature rises with it,
from 45.0Β Β°C to 60.03Β Β°C. Their difference stays constant under
the stipulated constant coefficient. In a real entrance region, β(π₯)
varies and a constant imposed flux need not give a constant difference.
The two cases prescribe different heat inputs, so their duties need not agree.
Answer
For fixed wall temperature, ππ=47.1Β Β°C and Λπ=2.27Γ103Β W, with a decreasing wall-to-bulk difference. For fixed heat flux, ππ=35.0Β Β°C and Λπ=1.26Γ103Β W, with a constant 25.0Β K wall-to-bulk difference.
A flat wall of fire-clay brick (0.5 m thick) initially has a uniform
temperature of 200K. One of its faces is suddenly exposed to hot gas at
1200K. If the heat transfer coefficient on the hot side is
7.38Β W/m2.K and the other face of the wall is insulated, determine:
The time it takes for the centre of the wall to reach 600K
Properties of fire-clay brick:
π=1.125Β W/m.K;
πΆπ=919Β J/kg.K;
π=2310Β kg/m3
Estimate the heat transfer coefficient, assuming that it stays constant
throughout the transient stage! Does the number quoted above make sense?
Solution
Reflecting the wall across its insulated face allows us to use the symmetric
plane-wall solution. The half-thickness labelled πΏ in the schematic is
denoted π below, to distinguish it from the vertical convection length.
Data
Quantity
Value
Actual wall thickness
π=0.50Β m
Initial brick temperature
ππ=200Β K
Hot-gas temperature
πβ=1200Β K
Target temperature
π=600Β K
Brick conductivity
ππ€=1.125Β Wβmβ1βKβ1
Brick density
ππ€=2310Β kgβmβ3
Brick specific heat capacity
πΆπ,π€=919Β Jβkgβ1βKβ1
Specified transient coefficient
β=7.38Β Wβmβ2βKβ1
Air conductivity at 700 K
π=0.0530Β Wβmβ1βKβ1
Air kinematic viscosity at 700 K
π=6.58Γ10β5Β m2βsβ1
Air Prandtl number at 700 K
Pr=0.697
The lecturerβs initial-film air properties allow an independent check of the specified coefficient. The slide supplies no wall height for this convection estimate.
Assumptions
One-dimensional transient conduction through the wall.
Constant solid properties and heat transfer coefficient.
No internal heat generation.
Heat transfer through the edges is neglected.
Radiation is neglected in the classroom model.
A vertical panel of height πΏ=1.0Β m is used only to estimate the convection coefficient.
Calculations
First check the supplied coefficient. Initially the exposed wall is at
200Β K, giving
ππ=1200+2002=700Β K
and
π½=1/ππ for ideal-gas air. The vertical-plate Rayleigh number is
Churchill-Chu gives about 8.36Β Wβmβ2βKβ1 for the
same conditions. The supplied value 7.38Β Wβmβ2βKβ1
is therefore a plausible convection estimate. The very large temperature
range makes constant-property natural-convection calculations approximate.
Use the supplied β=7.38Β Wβmβ2βKβ1 for the requested
transient calculation.
Measure π₯ from the insulated face, with the exposed face at π₯=π.
The physical centre of the brick is at π₯=π/2. Reflection across the
insulated face produces an equivalent symmetric slab of thickness 2π=1.0Β m.
Its half-thickness is π, so it is 0.50Β m, rather than
0.25Β m, that belongs in the slab Biot and Fourier numbers.
The dimensional time follows from Fo=πΌπ€π‘/π2:
π‘=Foπ2πΌπ€β0.3445(0.50)25.29938Γ10β7β1.625Γ105Β sβ45.15Β h
The calculated Fourier number exceeds 0.2, supporting the first-term
approximation. The insulated face is the centre of the reflected slab,
but it is not the physical centre of the original brick. At π₯=0 the
cosine is one, giving about 62.65Β h to reach 600Β K.
These times describe different positions. If the estimated
β=6.95Β Wβmβ2βKβ1 replaces the stipulated value,
the physical-centre result is about 46.51Β h.
The diffusion scale π2/πΌπ€β131Β h is consistent
with these heating times.
Freezing β is approximate because both the driving temperature difference
and fluid properties change during heating. Radiation can matter much more.
For example, with emissivity 0.8, a 600Β K surface facing
1200Β K surroundings has an effective radiation coefficient
about 147Β Wβmβ2βKβ1. The convection-only model would
therefore overestimate heating time in a furnace with those radiative conditions.
Answer
Using the supplied constant coefficient, the physical centre of the wall
reaches 600Β K after approximately 45.2Β h.
The separate natural-convection estimate gives
ββ6.95Β Wβmβ2βKβ1, supporting the order of
magnitude of the stated 7.38Β Wβmβ2βKβ1.
Determine the surface temperature of a cylindrical tank (0.75 m Diameter,
1.2m high), which contains a transformer immersed in an oil bath, producing
a uniform surface temperature condition. All heat loss from the surface
may be assumed to be due to natural convection to surrounding air at 295 K.
The heat dissipation rate from the transformer is constant at 1.5 kW.
The tank side and top have different natural-convection characteristic lengths.
The base is taken to rest on the floor, as in the lecturerβs calculation.
Figure 3: Cylindrical transformer tank with diameter, height, ambient temperature and separate characteristic lengths for its side and top.
Data
Quantity
Value
Tank diameter
π·=0.75Β m
Tank height
π»=1.2Β m
Heat dissipation
Λπ=1500Β W
Ambient air temperature
πβ=295Β K
Air conductivity at 347.5 K
π=0.0294Β Wβmβ1βKβ1
Air kinematic viscosity at 347.5 K
π=2.05Γ10β5Β m2βsβ1
Air Prandtl number at 347.5 K
Pr=0.718
Air conductivity at 340 K
π=0.0289Β Wβmβ1βKβ1
Air kinematic viscosity at 340 K
π=1.98Γ10β5Β m2βsβ1
Air Prandtl number at 340 K
Pr=0.720
Ideal-gas expansion coefficient
π½=1/ππ
The oil bath gives an approximately uniform surface temperature. Use the lecturerβs air properties at the two trial film temperatures.
Assumptions
Steady state.
The floor-supported base is adiabatic.
The side behaves as a vertical plate.
Fluid properties are constant at each trial film temperature.
Calculations
The side area and vertical characteristic length are
π΄side=ππ·π»=π(0.75)(1.2)=2.82743Β m2
The side characteristic length is
πΏside=π»=1.2Β m
For the upward-facing top, use area divided by perimeter:
π΄top=ππ·24=π(0.75)24=0.441786Β m2
The top characteristic length is
πΏtop=π΄topππ·=π·4=0.1875Β m
The two surface heat rates add. Define their combined conductance πΊ:
The coefficients depend on the unknown surface temperature through Rayleigh
number and film properties. Guess ππ , calculate the heat loss and adjust
the guess until the heat loss equals 1500Β W.
For the first trial, let ππ =400Β K. Then
Ξπ=105Β K and
For the side, use the turbulent vertical-plate correlation
Β―Β―Β―Β―Β―Β―Β―Β―Nuπ»=0.10Ra1/3π», applicable here above
Raπ»β109. For the hot upward-facing top, use
Β―Β―Β―Β―Β―Β―Β―Β―Nuπ·/4=0.15Ra1/3π·/4 in the
range 107β²Raπ·/4β²1011.
Substitute the coefficients and areas into Eq. 27:
Λπ400=[5.0488(2.82743)+7.5732(0.441786)](400β295)=1850.2Β W
The first guess loses more than 1500Β W, so reduce the surface
temperature. At ππ =385Β K, the film is at 340Β K
and Ξπ=90Β K. Repeating the same relations gives
For the small final correction, hold the 340Β K property set
fixed. Both correlations have exponent 1/3, so ββΞπ1/3
and the heat balance gives ΛπβΞπ4/3. Therefore,
Ξπ=90(15001528.0)3/4=88.760Β K(30)
The surface temperature is
ππ =295+88.760=383.760Β K.
Its film temperature is 339.38Β K, close to the property-table
value used for the final correction. Further precision would require updating
those properties. The side supplies about 81 percent of the heat loss.
Check the vertical-cylinder approximation using Eq. 20.
At the corrected temperature, Grπ»β1.13Γ1010,
so the required diameter is
π·minβ‘=35π»Gr1/4π»β35(1.2)(1.13Γ1010)1/4β0.129Β m
The actual diameter 0.75Β m exceeds this value. Radiation,
excluded by the question, could be significant in an actual tank. For
π=0.8 and surroundings at 295Β K, the radiation
coefficient near the calculated temperature is about
7.2Β Wβmβ2βKβ1, comparable to convection.
Including radiation would reduce the surface temperature. A raised tank
would also lose heat from its downward-facing base, requiring that surfaceβs
own natural-convection correlation.
Answer
The convection-only estimate is ππ β384Β K, or
111Β Β°C, for an adiabatic floor-supported base.
A long horizontal metal cylinder of diameter 0.180Β m starts
uniformly at 180Β Β°C in still air at 30Β Β°C.
Use solid properties ππ€=45Β Wβmβ1βKβ1,
π=7800Β kgβmβ3 and πΆπ=500Β Jβkgβ1βKβ1.
At the initial film temperature, use air properties
π=0.032Β Wβmβ1βKβ1,
π=2.40Γ10β5Β m2βsβ1 and Pr=0.70.
Estimate the natural-convection coefficient and the centre temperature after
20 minutes. Neglect end effects and radiation, and hold the estimated
coefficient and properties constant over the cooling interval. If a fan
instead produces a coefficient of 30Β Wβmβ2βKβ1, can the same transient model
still be used?
SolutionFigure 4: The horizontal cylinder loses heat from its curved surface. In the still-air case there is no imposed air speed; the fan case supplies an effective coefficient. End losses are neglected.
Data
Quantity
Value
Cylinder diameter
π·=0.180Β m
Initial temperature
ππ=180Β Β°C
Air temperature
πβ=30Β Β°C
Cooling time
π‘=1200Β s
Solid conductivity
ππ€=45Β Wβmβ1βKβ1
Solid density
ππ€=7800Β kgβmβ3
Solid specific heat capacity
πΆπ,π€=500Β Jβkgβ1βKβ1
Air conductivity
π=0.032Β Wβmβ1βKβ1
Air kinematic viscosity
π=2.40Γ10β5Β m2βsβ1
Air Prandtl number
Pr=0.70
Coefficient with the fan
β=30Β Wβmβ2βKβ1
Assumptions
Uniform initial solid temperature.
Constant properties and estimated heat transfer coefficient.
End effects and radiation are neglected.
Ideal-gas air for the expansion coefficient.
Calculations
Still air permits buoyancy-driven motion, so use natural convection around
a horizontal cylinder. The initial film temperature is
ππ=180+302+273.15=378.15βK.
For ideal-gas air,
π½=1/ππ=0.002644 Kβ»ΒΉ.
The convection correlation uses the cylinder diameter:
Diameter in Rayleigh number and π·/4 in the lumped Biot number answer
different questions. Rayleigh number describes buoyant flow around the body;
the Biot number compares internal and external thermal resistances.
With the fan, recalculate
Biπ=30(0.045)45=0.0300.
Lumped capacitance remains valid, and the same exponential gives 152.2Β Β°C after 20 minutes. If the forced coefficient had not been
supplied, it would require a suitable cross-flow correlation and a new Biot
check. One cannot assume that every increase in β preserves lumped behaviour.
Natural-convection β changes as the cylinder cools. The small 6.8 K
temperature change supports this initial-coefficient estimate over 20 minutes.
For a larger change, update properties and β(π) in the transient balance.
Radiation is excluded by the problem; this is a convection-only estimate.
Answer
The initial natural-convection coefficient is Β―Β―Β―Β―β=6.79Β Wβmβ2βKβ1, giving a centre temperature of approximately 173.2Β Β°C after 20 minutes. With the supplied fan coefficient, lumped capacitance remains valid and gives 152.2Β Β°C.
Water in the bottom of a narrow metal tube is held at a constant
temperature of 293 K. The total pressure of air (assumed dry) is
1.01325Γ105Β Pa and the temperature is 293 K. Water
evaporates and diffuses through the air in the tube, and the diffusion path
π§2βπ§1 is 0.1524 m. Calculate the rate of evaporation at steady state
in kgmol/(sβ m2). The diffusivity of water vapour at 293 K
and 1 atm pressure is 0.250Γ10β4Β m2/s. Assume that
the system is isothermal. The vapour pressure of water at 1 atm and 293 K
is 17.535 mmHg.
Additional note:
Moving boundaries or pseudo-steady state diffusion. How long would it take
for the water level to decrease by 0.1 m?
SolutionFigure 3: Water vapour diffuses through the gas above the liquid to the tube mouth. As the liquid surface falls, the gas diffusion distance increases.
The evaporation path is the gas space between the water surface and the tube mouth, as shown in Figure 2. Water vapour moves upwards. Air has zero net flux across the liquid interface.
Data
Quantity
Value
Initial gas diffusion length
πΏ0=0.1524Β m
Temperature
π=293Β K
Pressure
π=101325Β Pa
Diffusivity
π·π΄π΅=2.50Γ10β5Β m2βsβ1
Surface vapour pressure
ππ΄,π =17.535Β mmHg
Liquid density for the moving-boundary calculation
ππΏ=1000Β kgβmβ3
Water molar mass
ππ΄=18Β kgβkmolβ1
Final gas diffusion length
πΏπ=0.1524+0.1000=0.2524Β m
The source unit kgmol is equivalent to kmol. Liquid density and molar mass are supplementary inputs needed to convert liquid-depth loss to an amount evaporated.
Assumptions
Unidirectional diffusion through the gas space.
Steady gas profile for the initial flux calculation.
Pseudo-steady gas profile as the liquid surface recedes.
Ideal gas behaviour.
Constant temperature, pressure and diffusivity.
Vapour-liquid equilibrium at the water surface.
Zero net air flux through the liquid surface.
Negligible natural convection.
Calculations
At the liquid surface, equilibrium fixes the vapour partial pressure. At the mouth, the dry air gives ππ΄2=0. Converting the surface vapour pressure to Pa gives
ππ΄1=17.535760(101325)=2337.8Β Pa
The surface mole fraction and total gas concentration are
The initial rate is about 0.248Β kgβmβ2 per day. The stagnant-carrier flux is only 1.17 percent above the equimolar estimate because the water vapour is dilute.
For the falling liquid level, the gas length increases and the evaporation flux decreases. Define the constant logarithmic driving force,
As a check, removing 0.1Β m of water requires 5.556Β kmolβmβ2. Even if the initial flux stayed constant, evaporation would take
5.5561.593Γ10β7=3.49Γ107Β s.
The falling level lengthens the diffusion path, so the calculated time must be greater than that lower bound.
The gas diffusion time is πΏ20/π·π΄π΅β929Β s, much shorter than the liquid retreat time. This supports the pseudo-steady approximation. Maintaining the stipulated isothermal, convection-free conditions for 536 days would be difficult in an experiment.
Answer
The initial evaporation flux is ππ΄=1.59Γ10β7Β kgmolβmβ2βsβ1, directed towards the tube mouth. The corresponding mass flux is 2.87Γ10β6Β kgβmβ2βsβ1.
Under the pseudo-steady moving-boundary model, the water level takes 4.63Γ107Β s, or about 536 days, to fall by 0.1Β m.
Calculate the rate of diffusion of butanol at 20Β Β°C under
unidirectional steady state conditions through a 0.1 cm thick film of
water when the concentrations of butanol at the opposite sides of the
film are, 10% and 4% butanol by weight, respectively.
π·π΄π΅ (butanol in water solution)=5.9Γ10β6Β cm2/s.
The densities of 10% and 4% butanol solutions at 20Β Β°C may
be taken as 0.971 and 0.992 g/cc, respectively. Molecular weight of Butanol
is 74, and that of water is 18.
Solution
Data
Quantity
Value
Film thickness
πΏ=0.1Β cm
Diffusivity
π·π΄π΅=5.9Γ10β6Β cm2βsβ1
Butanol mass fraction at side 1
π€π΄1=0.10
Density at side 1
π1=0.971Β gβcmβ3
Butanol mass fraction at side 2
π€π΄2=0.04
Density at side 2
π2=0.992Β gβcmβ3
Butanol molar mass
ππ΄=74Β gβmolβ1
Water molar mass
ππ΅=18Β gβmolβ1
Assumptions
Steady state.
Unidirectional diffusion across the film.
Isothermal conditions at 20Β Β°C.
No chemical reaction.
Zero net water flux.
Constant diffusivity.
Total molar concentration approximated by its arithmetic mean across the film.
Calculations
Let A be butanol and B be water, with positive flux from the 10 percent side to the 4 percent side. The stagnant-water model uses mole fractions, so first convert the supplied mass fractions. On a 100Β g basis, side 1 contains 10/74 mol butanol and 90/18 mol water. Hence
π₯π΄1=10741074+9018=0.0263158
For side 2, the same conversion gives
π₯π΄2=474474+9618=0.0100334
The total molar concentration equals the total moles divided by the solution volume. Since a 100Β g sample occupies 100/ππΒ cm3,
The total molar concentration changes by about 7 percent across the film. Treating it as uniform is an approximation for this liquid mixture. Constant temperature and pressure alone do not make a liquidβs molar concentration constant.
With zero net water flux, the governing total-flux relation is
The full result is 1.85 percent larger, consistent with dilute butanol.
The flux is about 3.23Β kgβmβ2 per day. Comparing it with the water-evaporation example requires the concentration difference and film thickness as well as the diffusivity.
Answer
Butanol diffuses from the 10 percent side to the 4 percent side at ππ΄=5.05Γ10β7Β kmolβmβ2βsβ1, equivalent to a mass flux of 3.74Γ10β5Β kgβmβ2βsβ1.
A binary ideal gas is at 300 K and 100 kPa. Across a 5.0 mm plane film,
the mole fraction of A falls from 0.30 to zero. Use
π·π΄π΅=2.0Γ10β5 mΒ² sβ»ΒΉ and π =8.314 J molβ»ΒΉ Kβ»ΒΉ.
Compare the steady molar flux of A for equimolar counter-diffusion and for
zero net flux of B. Take positive flux from the high-A face to the low-A
face. Explain the difference and determine the flux ratio if the high-face
mole fraction is reduced to 0.010, with all other data unchanged.
Solution
Data
Quantity
Value
Temperature
π=300Β K
Pressure
π=100000Β Pa
Film thickness
πΏ=0.0050Β m
Diffusivity
π·π΄π΅=2.0Γ10β5Β m2βsβ1
Gas constant
π =8.314Β Jβmolβ1βKβ1
Side 1 mole fraction in the first comparison
π¦π΄1=0.30
Side 2 mole fraction
π¦π΄2=0
Side 1 mole fraction in the second comparison
π¦π΄1=0.010
Assumptions
Steady state.
One-dimensional diffusion across a plane film.
Ideal gas behaviour.
Constant temperature, pressure and diffusivity.
No reaction in the film.
Calculations
Positive flux points from the high-A face towards the low-A face. The total concentration is
A enters the film without a compensating net flux of B, producing molar flow in Aβs direction. That bulk flow transports additional A. B still diffuses towards the high-A side, but its bulk contribution cancels its diffusion.
so the ratio approaches one as π¦β0. The two models then give nearly the same flux.
Answer
At π¦π΄1=0.30, the equimolar flux is 0.0481Β molβmβ2βsβ1 and the stagnant-carrier flux is 0.0572Β molβmβ2βsβ1. The latter is 18.9 percent larger. At π¦π΄1=0.010, their ratio falls to 1.005.
A given gas A can diffuse through the solid wall of a tube in which it
flows. Find the rate of gas diffusion in terms of the tube dimensions
(inside wall radius π π, outside wall radius π π), the gas-solid
diffusivity, and the gas concentrations at the walls.
Assumptions:
steady state
flux in π direction only
no chemical reaction
A concentration in gas is quite dilute βπ₯π΄βΌ0
Constant π·π΄π΅
Solution
The radial shell in Figure 3 shows why a cylindrical balance must account for the changing area.
Figure 3: A stationary cylindrical wall and an enlarged radial shell. The molar rate equals the local flux multiplied by the cylindrical area 2ΟrL.
Data
Quantity
Value
Inner wall radius
π π
Outer wall radius
π π
Axial length
πΏ
Constant gas-solid diffusivity
π·π΄π΅
Dissolved concentration at the inner face
ππ΄π in kmolβmβ3
Dissolved concentration at the outer face
ππ΄π in kmolβmβ3
Assumptions
Steady state.
Radial diffusion only.
No chemical reaction.
Dilute dissolved A in a stationary solid matrix.
Constant diffusivity.
Negligible end effects.
Calculations
Choose positive flux radially outwards. A steady shell has no accumulation or reaction, so its inward and outward molar rates are equal. From Eq. 10, the governing radial balance is
πππ(πππ΄π)=0
Consequently, the total outward molar rate is independent of radius,
Λππ΄=2πππΏππ΄π
The source denotes this rate by ππ΄. In a stationary solid, transport of dilute dissolved A follows Fickβs law,
ππ΄π=βπ·π΄π΅πππ΄ππ
Substitute the flux into the rate expression, separate the variables and integrate between the two wall faces,
The concentration profile follows from the same differential equation. Since πβπππ΄/ππ is constant, integration from π π to an arbitrary radius π gives
The profile recovers both specified wall concentrations. When ππ΄π>ππ΄π, its gradient is negative and the rate is positive outwards. The units of Eq. 14 are
just as cylindrical conduction can be written using a log-mean area.
The specified concentrations must be concentrations dissolved in the solid. A bulk gas concentration requires a solubility relation before it can be used as a wall boundary condition. For example, ππ΄=ππ΄ππ΄, with ππ΄ in kmolβmβ3βPaβ1, gives the permeability π·π΄π΅ππ΄.
A spherical ball of solid, nonporous naphthalene is suspended in still air.
The naphthalene ball slowly sublimes, releasing naphthalene into the
surrounding air. It may be assumed that the sublimation occurs as a
molecular diffusion-limited process. Through consideration of the
diffusion equation, estimate the time required to reduce the diameter of
the ball from 2 to 0.5 cm when the surrounding air is at 347 K and 1 atm.
Data for naphthalene:
MW = 128 g/mol
Density = 1.145 g/cm3
Diffusivity in air = 8.19Γ10β6Β m2/s
Vapour pressure = 666 Pa at 347 K
Solution
The spherical geometry in Figure 2 applies at each instantaneous ball radius. Vapour diffuses outwards through air while the solid surface slowly recedes.
Data
Quantity
Value
Temperature
π=347Β K
Pressure
π=1Β atm=101325Β Pa
Diffusivity
π·π΄π΅=8.19Γ10β6Β m2βsβ1
Saturation pressure
πsatπ΄=666Β Pa
Solid density
ππ =1145Β kgβmβ3
Molar mass
ππ΄=128Β kgβkmolβ1
Initial radius
π π=0.010Β m
Final radius
π π=0.0025Β m
The supplied dimensions are diameters. The table gives the corresponding radii.
Assumptions
Spherically symmetric diffusion.
Pseudo-steady gas profile as the surface recedes.
Isothermal sublimation at 347Β K.
Ideal gas behaviour.
Constant pressure, diffusivity and solid density.
Vapour-solid equilibrium at the surface.
Negligible naphthalene concentration far from the ball.
Zero net air flux through the solid surface.
Negligible natural convection and gas-phase reaction.
Calculations
From Eq. 11 at a fixed instantaneous radius, the steady radial balance gives
π(π2ππ΄π)ππ=0.
The same total rate crosses each concentric sphere,
Λππ΄=4ππ2ππ΄π
Air cannot cross the solid surface, so ππ΅π=0. Writing Eq. 12 in radial form retains the resulting molar flow,
The instantaneous sublimation rate is then Λππ΄=4ππ ππ·π΄π΅πΉ. The sphere loses moles as its radius decreases. Differentiating its solid inventory gives
much shorter than the sublimation time. This supports a pseudo-steady gas profile. Replacing πΉ by the dilute approximation π¦π΄π gives 61.60 h, only 0.33 percent higher.
The integrated balance can also be expressed as a diameter-square law,
Here πΎ has units m2βsβ1. The composition factor must remain in πΎ because the vapour pressure determines the driving force. A different thermal model would be needed if sublimation cooling changes the surface temperature appreciably.
Answer
The diffusion-limited model predicts π‘=2.21Γ105Β s, or about 61.4 h, for the diameter to decrease from 2 cm to 0.5 cm.
Consider a catalytic tubular reactor, where a dilute solution of solute A
in solvent S is in fully developed, laminar flow in the region π§<0.
It encounters the catalytic wall in the region 0β€π§β€πΏ, where
solute A is instantaneously and irreversibly rearranged to an isomer B.
Write the appropriate diffusion equation (in differential form) for this
problem, along with the boundary conditions. Assume that the flow is
isothermal and neglect the presence of B.
SolutionFigure 7: Axial liquid flow in a tube with a catalytic wall. For the instantaneous wall reaction, the labelled wall concentration is zero. Species A is consumed at the wall as the liquid flows downstream.
In this diagram, π=2π and π=Β―Β―Β―Β―π£ππ 2. Its concentration labels
denote the molar concentration ππ΄ used below.
Data
Tube radius π
Catalyst-coated length 0β€π§β€πΏ
Mean axial velocity Β―Β―Β―Β―π£
Uniform inlet concentration ππ΄0
Wall reaction A converts instantaneously to its isomer B
The dilute liquid enters in fully developed laminar flow. The question asks for a model and boundary conditions, so no numerical solution is required.
Assumptions
Steady state.
Axisymmetric flow and concentration field.
Constant properties.
No radial velocity.
No homogeneous reaction in the liquid.
Uniform inlet concentration.
Axial diffusion neglected only when axial convection dominates over the length of interest.
The specified isothermal flow and instantaneous wall reaction are supplied features of the model.
Calculations
Reaction occurs on the wall, so it enters through a boundary condition. It does not supply a volumetric π π΄ term in the liquid. The steady, axisymmetric form of Eq. 10 is
For illustration, use Β―Β―Β―Β―π£=0.1Β mβsβ1, πΏ=0.1Β m and π·π΄π΅=10β9Β m2βsβ1. These are supplementary scale values rather than data supplied in the question. They give
PeπΏ=0.1(0.1)10β9=107
Axial convection then dominates axial diffusion over this length, away from axial end regions. The resulting working equation is
The reduced equation needs inlet data and two radial conditions. The inlet supplies uniform A. Instantaneous reaction holds the catalytic-wall concentration at zero. Centreline symmetry gives zero radial gradient. These conditions form the system
The inlet and wall data are interpreted as one-sided conditions at the inlet-wall corner, where the idealised concentration jumps. Every term in the governing equation has units kmolβmβ3βsβ1.
at π=π , because diffusive arrival must equal surface consumption.
This formulation is the mass-transfer Graetz problem. In the fully developed concentration region, a constant wall concentration gives Sh=ππππ‘/π·π΄π΅=3.66, where ππ‘=2π and ππ uses the difference between mixing-cup and wall concentrations. The entrance-region coefficient is larger. The thermal problem with the same geometry and boundary conditions gives Nu=3.66.
when axial diffusion can be neglected. Apply uniform inlet concentration, zero A concentration at the catalytic wall, and zero radial concentration gradient at the centreline, as stated above.
If axial diffusion is significant, use Eq. 30 and provide an additional axial boundary condition.
A dilute gas species A crosses a plane film to a catalyst, where it undergoes
an isomerisation π΄βΆπ΅. The gas is at 500 K and 100 kPa.
The bulk mole fraction is 0.020, film thickness is 1.5 mm, diffusivity is
6.0Γ10β5 mΒ² sβ»ΒΉ and surface rate constant is ππ =0.015 m sβ»ΒΉ.
The catalyst area is 0.040 mΒ². Assume steady one-dimensional isothermal
transport with constant pressure and diffusivity, no gas-phase reaction,
and negligible normal bulk velocity. Let π§=0 be the bulk edge and π§=πΏ
the catalyst.
Derive the concentration profile from a shell balance. Find the surface mole
fraction and molar removal rate. Which resistance is greater, and what
happens in the limits of very slow and very fast surface reaction?
SolutionFigure 8: A diffuses from the bulk through a nonreacting film to the catalyst. A finite reaction rate leaves a nonzero surface concentration, determined by matching film transport and surface consumption.
Data
Quantity
Value
Temperature
π=500Β K
Pressure
π=100000Β Pa
Bulk mole fraction
π¦π΄π=0.020
Film thickness
πΏ=0.0015Β m
Diffusivity
π·π΄π΅=6.0Γ10β5Β m2βsβ1
Surface reaction coefficient
ππ =0.015Β mβsβ1
Surface area
π΄π =0.040Β m2
Assumptions
Steady state.
One-dimensional transport across a plane film.
Ideal gas behaviour.
Constant temperature, pressure and diffusivity.
No gas-phase reaction.
Negligible normal bulk velocity.
Calculations
For a film slice of area π and thickness Ξπ§, no accumulation or
volumetric reaction occurs. The shell balance is
Both concentration-based resistances have units s mβ»ΒΉ:
πΏ/π·π΄π΅=25.0 and 1/ππ =66.7. Their sum gives
ππ΄=ππ΄ππΏπ·π΄π΅+1ππ .
Surface reaction supplies the larger
resistance, but neither resistance is negligible.
As ππ β0, ππ΄π βππ΄π and removal tends to zero.
As ππ ββ, ππ΄π β0 and removal approaches the
diffusion-limited value π·π΄π΅ππ΄π/πΏ.
The isomerisation exchanges one gas mole for one gas mole, consistent with
zero net molar flow. Other stoichiometries require checking that assumption.
For a cylindrical reactor with axial flow, the earlier catalytic-tube example
retains axial convection. Its finite-rate wall condition uses the same
principle: arriving diffusive flux equals surface consumption. Its radial
profile must be found with the axial concentration field, rather than assumed
linear as in this plane, nonconvecting film.
Answer
The surface mole fraction is π¦π΄π =0.0146 and the molar removal rate is Λππ΄=2.10Γ10β4Β molβsβ1. Surface reaction supplies the larger resistance, 66.7Β sβmβ1, compared with 25.0Β sβmβ1 for the film. The slow-reaction limit gives zero removal; the fast-reaction limit gives the diffusion-limited flux π·π΄π΅ππ΄π/πΏ.
At a catalyst, gas R reacts by π βΆ2π. Take the coordinate
direction towards the catalyst and suppose ππ =0.010 mol mβ»Β² sβ»ΒΉ.
Find ππ and the total molar flux. Is equimolar counter-diffusion
appropriate? Compare with a reaction π βΆπ at the same
reactant flux.
Solution
Data
Surface reaction π βΆ2π
Reactant flux toward the catalyst ππ =0.010Β molβmβ2βsβ1
Comparison reaction π βΆπ
Positive flux points toward the catalyst.
Assumptions
Steady surface reaction.
No surface accumulation.
Only the stated reaction produces or consumes R and P.
Calculations
Products move away from the surface. Two moles of P leave for each mole of
R arriving, so
the bulk term is negative.
It opposes Rβs inward diffusive transport.
The flux relation can also be written ππ =βππ/2. The negative sign
expresses opposing directions; the factor comes from stoichiometry.
For π βΆπ, ππ=βππ and the total molar flux is zero.
That reaction permits the equimolar model under the other stated binary-film
assumptions. Constant temperature and pressure alone do not determine the
flux ratio.
Answer
For π βΆ2π, ππ=β0.020Β molβmβ2βsβ1 and the total molar flux is β0.010Β molβmβ2βsβ1, directed away from the catalyst. Equimolar counter-diffusion does not apply. For π βΆπ, ππ=β0.010Β molβmβ2βsβ1 and the total molar flux is zero.
An idealised stationary water layer of depth 0.010 m absorbs oxygen at its
upper face. That face is held at ππ΄0=0.25 mol mβ»Β³. Oxygen consumption
throughout the liquid follows π π΄=βπβ²β²β²1ππ΄, with
πβ²β²β²1=0.0020 sβ»ΒΉ and π·π΄π΅=2.0Γ10β9 mΒ² sβ»ΒΉ. The bottom is
impermeable. Assume a steady one-dimensional profile and constant properties.
Formulate the problem and estimate the oxygen influx. Explain how the model
would change if consumption occurred only on a catalytic bottom instead.
Solution
Data
Liquid depth πΏ=0.010Β m
Fixed top concentration ππ΄0=0.25Β molβmβ3
Measure π₯ downwards from the upper surface. Consumption occurs throughout
the liquid, so it is a homogeneous reaction and belongs in the local balance:
for a first-order surface reaction.
The two rate constants have different units, sβ»ΒΉ in the liquid and m sβ»ΒΉ at
the surface.
A pond can require an unsteady absorption model or a steady reaction model.
Its size alone does not decide. The stated consumption, flow, initial state
and boundary conditions determine which model is appropriate.
Answer
The oxygen influx is 5.00Γ10β7Β molβmβ2βsβ1 into the liquid. If reaction occurs only at the bottom, use π·π΄π΅πβ²β²π΄=0 in the interior and the bottom condition
A catalyst converts gas A instantaneously by π΄βΆ3π΅.
A plane gas film of thickness πΏ=2.0 mm separates the surface from a
well-mixed bulk gas with π¦π΄0=0.20. The binary gas is at 500 K and
1.5 bar, with π·π΄π΅=2.0Γ10β5 mΒ² sβ»ΒΉ. Assume steady,
one-dimensional transport, constant temperature, pressure and diffusivity,
and no reaction inside the film. The catalyst is impermeable except for the
reactant and product exchange specified by the reaction.
Taking π§ from the bulk edge towards the catalyst, derive π¦π΄(π§) and
the inward flux of A. Compare with equimolar counter-diffusion across the
same film. Does a constant species flux imply a linear mole-fraction profile?
SolutionFigure 9: A crosses the gas film toward the catalyst. The surface reaction produces B, whose outward flux changes the mixture molar flow.
Data
Quantity
Value
Film thickness
πΏ=0.0020Β m
Bulk mole fraction
π¦π΄0=0.20
Temperature
π=500Β K
Pressure
π=150000Β Pa
Diffusivity
π·π΄π΅=2.0Γ10β5Β m2βsβ1
Instantaneous catalyst reaction
π΄βΆ3π΅
Assumptions
Steady state.
One-dimensional transport across a plane film.
Binary ideal gas mixture.
Constant temperature, pressure and diffusivity.
No gas-phase reaction.
Surface exchange occurs only through the stated reaction.
Calculations
A film shell balance gives πππ΄/ππ§=0 because there is no accumulation or
gas-phase reaction. The boundary conditions are π¦π΄(0)=π¦π΄0 and
π¦π΄(πΏ)=0. Surface stoichiometry fixes
ππ΅=β3ππ΄, so the total flux is ππ΄+ππ΅=β2ππ΄.
Product expansion creates outward molar flow even though there is no
externally imposed flow through the film.
The actual flux in this film model is 15.9 percent below that value.
Outward bulk flow opposes the inward transport of A. This reverses the
effect of the inward molar flow in the contracting 2π΄βΆπ΅
case above.
At π§=πΏ/2,
π¦π΄=β1.4β12=0.09161,
compared with
0.100 for a straight line. The species flux is constant, but its diffusive
and bulk contributions vary with composition. Constant flux therefore does
not require a linear profile. Both imposed endpoints are recovered exactly.
As π¦π΄0β0, the logarithm approaches 2π¦π΄0 and the dilute
equimolar estimate is recovered.
Answer
The inward flux is ππ΄=0.0607Β molβmβ2βsβ1, 15.9 percent below the equimolar estimate. The profile is
For the reaction and air composition above, find the particle radius and burnout time in terms of π 0, π, π·, ππΆ and ππΆ. Quantify the flux correction relative to pure diffusion.
Calculation and answer
Data
Quantity
Value
Initial particle radius
π 0
Gas total molar concentration
π
Effective oxygen diffusivity
π·
Carbon density
ππΆ
Carbon molar mass
ππΆ
Ambient oxygen mole fraction
π¦π2,β=0.21
Reaction stoichiometry
2.5 mol oxygen consumed per 3 mol carbon
Assumptions
Spherical particle with uniform solid density.
Pseudo steady diffusion in the surrounding gas.
Instantaneous surface reaction and no gas-phase reaction.
Constant temperature, total concentration and effective diffusivity.
Calculations
The steady gas calculation supplies the oxygen flux at each current radius. Writing Eq. 5 for the inward magnitude gives
The radius expression applies for 0β€π‘β€π‘π. The oxygen flux is 2.04 percent smaller than the pure-diffusion estimate. No numerical burnout time is determined without the carbon density and operating gas properties.
A spherical drug capsule (0.652 cm diameter) has an initial concentration
of 68.9 mg/cm3 Dramamine.
What is the residual concentration of Dramamine at the centre of the spherical bead after 48 hours?
The capsule is intended to release the drug in the gastrointestinal system.
Its nontoxic gelatinous material allows the drug to diffuse out gradually.
The drug within the bead is the source of A; the surrounding fluid is the sink.
Data: π·π΄π΅ (Dramamine in the gel matrix of the capsule) =
3Γ10β7Β cm2/s at body temperature of 37Β Β°C.
Calculation and answerFigure 5: Drug initially fills the spherical bead uniformly. It diffuses radially to a sink at the surface. Total release requires the concentration averaged over the whole bead.
The spherical diffusion model is the one shown in Figure 4. The maintained sink fixes the surface concentration, while the centre concentration is the requested unknown.
Data
Quantity
Value
Bead diameter
π=0.652Β cm
Bead radius
π π=0.326Β cm
Initial drug concentration
πΆππ=68.9Β mgβcmβ3
Drug diffusivity
π·=3Γ10β7Β cm2βsβ1
Duration
π‘=48Β h=172800Β s
Body temperature
π=37Β Β°C
Assumptions
Spherical, stationary gel with constant radius.
Uniform initial drug concentration.
Well-stirred sink maintains zero surface concentration.
No drug reaction within the gel.
Constant diffusivity.
Calculations
The sink removes drug from the bead surface while internal diffusion supplies it from the centre. A prescribed surface concentration corresponds to infinite mass Biot number. First calculate the Fourier number using the radius:
The centre retains 1.62 percent of its initial concentration while only 0.493 percent of the initial drug remains in the entire bead. The centre has the largest residual concentration, so these two results have the expected order. Finite external resistance or swelling-dependent diffusivity requires a different model.
Answer
After 48 hours, the concentration at the centre is approximately
πΆππ=1.12Β mgβcmβ3
The corresponding additional estimate is 99.5 percent of the drug released.
Chemical vapour deposition takes place at 1100Β Β°C, where
phosphorus diffuses into silicon. The surface concentration is maintained at
πΆπ΄π=2.5Γ1020 phosphorus atoms per cm3 of solid silicon.
For comparison, pure silicon contains about 5Γ1022 silicon atoms
per cm3, so the phosphorus concentration is dilute.
Treat the phosphorus-rich coating as an effectively infinite source relative
to the amount transferred. Predict the depth at which the phosphorus
concentration reaches one percent of the surface value after one hour,
2.5Γ1018 phosphorus atoms per cm3. Also predict the
concentration profile after one hour.
Calculation and answer
The coating and diffusion direction are shown in Figure 7. The calculated profile in Figure 8 can be used to check the recovered one-percent depth.
Diffusivity at 1100 degrees Celsius π·=6.5Γ10β13Β cm2βsβ1
Duration π‘=3600Β s
Target concentration
0.01πΆππ =2.5Γ1018Β atomsβcmβ3
Assumptions
One-dimensional diffusion into stationary silicon.
Uniform temperature and constant diffusivity.
No reaction in the silicon.
An effectively infinite coating maintains the surface concentration.
The wafer is sufficiently thick for the semi-infinite approximation.
Calculations
Phosphorus enters from the maintained coating and diffuses into the initially phosphorus-free silicon. The semi-infinite diffusion solution on an atomic-concentration basis is
Therefore erfβ‘(π)=0.99, and the inverse error function gives π=1.821386. The source table reading of 1.82 gives the same depth to three significant figures.
The one-hour diffusion length is
2βπ·π‘=2β6.5Γ10β13(3600)=9.67471Γ10β5Β cm
The requested depth follows from the definition of π:
At π₯=0.967471Β πm, the profile gives 3.93248Γ1019Β atomsβcmβ3. At π₯=2Β πm it gives 8.6523Γ1017Β atomsβcmβ3. Substitution of the calculated target depth recovers one percent of the surface concentration.
For an illustrative 300Β πm wafer, the disturbance-depth estimate is 4βπ·π‘=1.93494Β πm and the thickness-based Fourier number is 2.6Γ10β6. The far boundary is therefore well beyond the diffusion region. The one-percent criterion defines a concentration depth; an electrical junction depth also requires the waferβs background dopant concentration.
Answer
After one hour the one-percent concentration depth is 1.76Β πm. The one-hour phosphorus concentration profile is given by Eq. 27.
A stationary polymer initially contains a uniform solute concentration of
0.20 mol mβ»Β³. At π‘=0 one face is held at 2.0 mol mβ»Β³ by a maintained
reservoir. Use constant diffusivity π·π΄π΅=1.0Γ10β10 mΒ² sβ»ΒΉ,
negligible external resistance and no reaction. The polymer is 0.020 m thick
and transport is normal to the exposed face; neglect edges.
After one hour, find the concentration 0.60 mm below the surface, the
instantaneous surface flux and the cumulative uptake per unit exposed area.
Check the semi-infinite approximation, and describe the effect of increasing
the duration to four hours.
SolutionFigure 9: Solute enters an initially uniform body from a maintained surface concentration. Position is measured inward from the exposed face.
Data
Quantity
Value
Initial concentration
ππ΄π=0.20Β molβmβ3
Surface concentration
ππ΄π =2.0Β molβmβ3
Diffusivity
π·π΄π΅=10β10Β m2βsβ1
Evaluation depth
π₯=0.000600Β m
Evaluation time
π‘=3600Β s
Body thickness
π»=0.020Β m
Assumptions
One-dimensional diffusion into a stationary polymer.
Uniform initial concentration and maintained surface concentration.
Negligible external resistance and edge effects.
No reaction and constant diffusivity.
Calculations
At 3600 s, βπ·π΄π΅π‘=0.000600 m. The estimated disturbance depth
is 4βπ·π΄π΅π‘=0.00240 m, well short of the 0.020 m far boundary.
The semi-infinite solution is appropriate on this scale.
Multiplying the final flux by the duration would give half this uptake.
The surface flux was greater earlier and falls as π‘β1/2.
After four hours, the disturbance depth and cumulative uptake both double,
while instantaneous surface flux halves. The disturbance depth is then
0.00480 m, still below the far boundary. At a fixed depth the concentration
also changes; only points at equal π₯/βπ‘ have the same concentration
ratio. These results require the maintained surface concentration and
constant diffusivity specified in the problem.
Answer
After one hour, ππ΄=1.06Β molβmβ3 at 0.60 mm, ππ΄(0,π‘)=1.69Γ10β7Β molβmβ2βsβ1 and ππ΄=1.22Γ10β3Β molβmβ2. At four hours the uptake and disturbance depth double and the surface flux halves.
A 0.10 m deep layer of still water at 298 K initially contains dissolved
molecular oxygen uniformly at 10 g mβ»Β³. At π‘=0 a regulated pure-oxygen
headspace maintains an oxygen partial pressure of 2.0 atm above the water.
Use Henryβs law in the form ππ΄=π»ππ΄π with
π»=0.80 atm mΒ³ molβ»ΒΉ, and π·π΄π΅=2.1Γ10β9 mΒ² sβ»ΒΉ.
The molar mass of molecular oxygen is 32.0 g molβ»ΒΉ.
How long does the concentration 5.0 mm below the interface take to reach
20 g mβ»Β³? Assume negligible gas-side resistance, rapid interfacial
equilibration, constant properties, no reaction and no liquid motion.
Check whether a semi-infinite approximation is consistent with the result.
For the calculation,
erfβ‘(1.03608)=0.857143.SolutionFigure 10: Oxygen enters water through the top surface. Henry's law sets the dissolved surface concentration; the target is at a specified depth below that surface.
Data
Quantity
Value
Water depth
πΏ=0.10Β m
Oxygen diffusivity
π·π΄π΅=2.1Γ10β9Β m2βsβ1
Oxygen partial pressure
ππ΄=2.0Β atm
Henryβs constant
π»=0.80Β atmβm3βmolβ1
Initial mass concentration
10Β gβmβ3
Target mass concentration
20Β gβmβ3
Evaluation depth
π₯=0.0050Β m
Oxygen molar mass
ππ΄=32.0Β gβmolβ1
Assumptions
One-dimensional diffusion in stationary water.
Uniform initial concentration.
Maintained gas partial pressure and rapid interfacial equilibrium.
Negligible gas-side resistance.
Constant properties and no reaction.
Calculations
The initial concentration and the interfacial equilibrium concentration are
different. Henryβs law fixes the latter on a molar basis:
The regulated headspace makes this a maintained surface boundary, as
required by Eq. 22. The surface value would change if gas were
depleted from an unregulated finite headspace.
Using mass concentration πΆπ, which is proportional to molar concentration,
the semi-infinite solution is
Thus erfβ‘(π)=6/7 and π=1.03608. Solving for time,
π‘=π₯24π·π΄π΅π2=(0.0050)24(2.1Γ10β9)(1.03608)2=2773Β s=46.2Β min
The disturbance-depth estimate is 4βπ·π΄π΅π‘=0.00965 m, much
less than the 0.10 m water depth. The far boundary is therefore unimportant
on this time scale within the stated diffusion-only model.
At twice the target depth, reaching the same concentration requires four
times the duration, about 3.08 hours. Using 16 g molβ»ΒΉ would describe atomic
oxygen and impose the wrong surface concentration for this problem.
Natural convection, mixing or biological consumption would require revisiting
the stated model rather than changing only the diffusion time.
Answer
The target concentration is reached after approximately 2.77Γ103Β s=46.2Β min. The disturbance depth is about 9.65 mm, so the semi-infinite approximation is consistent with the 100 mm water depth.
Air flows at 2 m sβ»ΒΉ over a hypothetical volatile solid A with an exposed
area of 10β3Β m2. The measured sublimation rate is
2.29Γ10β4Β molβsβ1. Both the air and the surface are
maintained at 293 K; external heating supplies the latent heat. The pressure
is 1.013Γ105Β Pa and incoming air contains negligible A.
For this model solid, use the stipulated equilibrium vapour pressure
ππ΄,π =4740Β Pa and diffusivity in air
π·π΄π΅=1.5Γ10β5Β m2βsβ1. Use air properties
π=1.55Γ10β5Β m2βsβ1, Pr=0.71,
πΆπ=1007Β Jβkgβ1βKβ1 and molar mass
ππ΅=0.029Β kgβmolβ1.
Determine the concentration-basis mass transfer coefficient and estimate
the heat transfer coefficient using the Chilton-Colburn analogy. This
adaptation replaces the inconsistent CO2 interface state in slide 17.
Solution
Data
Quantity
Value
Exposed area
π΄π =10β3Β m2
Measured molar rate
Λππ΄=2.29Γ10β4Β molβsβ1
Gas temperature
π=293Β K
Gas pressure
π=101300Β Pa
Species diffusivity
π·π΄π΅=1.5Γ10β5Β m2βsβ1
Kinematic viscosity
π=1.55Γ10β5Β m2βsβ1
Surface equilibrium pressure
ππ΄,π =4740Β Pa
Prandtl number
Pr=0.71
Specific heat capacity
πΆπ=1007Β Jβkgβ1βKβ1
Air molar mass
ππ΅=0.029Β kgβmolβ1
These values are stipulated for the hypothetical solid and air.
Assumptions
Steady sublimation over the exposed area.
Negligible A in the bulk air.
Ideal gas behaviour and constant supplied properties.
Low normal mass flux and matching heat and mass transfer conditions.
Calculations
The gas-side surface concentration follows from ideal-gas equilibrium:
The surface fraction is about 0.0468, so the low-flux approximation is plausible but not exact.
External heating maintains the specified surface temperature. The estimated heat coefficient describes how this flow would respond to a small temperature difference; it does not supply latent heat when the air and surface are isothermal.
Answer
The concentration-basis mass transfer coefficient is ππβ0.118Β mβsβ1. The corresponding heat transfer coefficient is ββ184Β Wβmβ2βKβ1.
A stream of air at 100 kPa pressure and 290 K is flowing on the top
surface of a thin flat sheet of solid naphthalene of length 1 m with a
velocity of 20 m/s and surface temperature of 310K. Calculate:
The average mass transfer coefficient over the flat plate
The rate of loss of naphthalene from the surface per unit width
Data @ ππ=300Β K:
π·π΄π΅=6Γ10β6Β m2/s (naphthalene vapour in air)
π=1.5Γ10β5Β m2/s (kinematic viscosity of air)
πΆπ΄π =1Γ10β5Β kmol/m3 (conc. of naphthalene at the air-solid interface)
SolutionFigure 3: Air flows along the naphthalene sheet. The velocity boundary layer is the same type as for a heated plate; the concentration boundary layer replaces the thermal boundary layer in the mass transfer analogy.
Use the concentration boundary layer in Figure 2, with the leading edge at the start of the sheet. The one-metre dimension is along the flow; width remains unspecified.
Data
Sheet length πΏ=1Β m
Air speed π=20Β mβsβ1
Film diffusivity π·π΄π΅=6Γ10β6Β m2βsβ1
Film kinematic viscosity π=1.5Γ10β5Β m2βsβ1
The table uses 0.037, which changes this estimate by about one percent.
For an untripped plate, the laminar starting length is
π₯π‘π=5Γ105ππ=0.375βm.
The mixed average relation
Β―Β―Β―Β―Β―Β―Β―ShπΏ=(0.037Re0.8πΏβ871)Sc1/3
gives 2806.4, Β―Β―Β―Β―ππ=0.01684 m sβ»ΒΉ and Λππ΄/width=1.684Γ10β4 mol sβ»ΒΉ mβ»ΒΉ.
This is the preferred estimate when transition occurs naturally at the stated Reynolds number.
Both answers are conditional on the specified flow model; they are not interchangeable rounding choices.
With ππ΄=0.128 kg molβ»ΒΉ the source fully turbulent estimate corresponds to 2.62 kg dayβ»ΒΉ per metre of width.
Answer
With the sourceβs fully turbulent model, Β―Β―Β―Β―ππ=0.0237Β mβsβ1 and the loss per unit width is 2.37Γ10β4Β molβsβ1βmβ1. For an untripped plate with natural transition, the preferred estimates are 0.0168Β mβsβ1 and 1.68Γ10β4Β molβsβ1βmβ1.
Dry air at 1.013Γ105Β Pa blows across a thermometer
whose bulb has been covered with a dampened wick. The classical βwet bulbβ
thermometer indicates a steady-state temperature is reached by a small
amount of liquid water evaporating into a large reservoir of unsaturated
vapour-gas mixture. The thermometer reads 290 K. At this temperature,
the following properties were evaluated:
Property
Value
Vapour pressure of water
1.94Γ103Β Pa
Density of air
1.219 kg/m3
Latent heat of vaporisation of water
2461 kJ/kg
Prandtl number
0.71
Schmidt number
0.61
Specific heat of air
1.006 kJ/kg.K
What is the temperature of the dry air?
Solution
Data
Quantity
Value
Wet-bulb temperature
ππ =290Β K
Surface vapour pressure
ππ΄,π =1940Β Pa
Air density
π=1.219Β kgβmβ3
Specific heat capacity
πΆπ=1006Β Jβkgβ1βKβ1
Prandtl number
Pr=0.71
Schmidt number
Sc=0.61
Latent heat
π=2.461Γ106Β Jβkgβ1
Water molar mass
ππ΄=0.018Β kgβmolβ1
Assumptions
Steady wick temperature and evaporation rate.
Dry incoming air and saturated gas at the wet surface.
Radiation and conduction through the thermometer stem are neglected.
Sensible energy transport by injected vapour and blowing are neglected.
Constant properties and matching heat and mass transfer conditions.
Calculations
At steady state convection supplies evaporationβs latent heat.
Substituting the latent heat, concentration and coefficient ratio gives
πββππ =2.461Γ106(0.018)(0.80465)(9.023Γ10β4)=32.161Β K
Adding this temperature difference to the wet-bulb reading gives
πβ=290+32.161=322.16Β K
Velocity cancels because the same flow affects both coefficients.
This cancellation holds under the stated analogy and energy assumptions; it does not remove ventilation requirements or radiative errors of a real thermometer.
Answer
The estimated dry-air temperature is πβ=322.2Β K, approximately 49.0βΒ°C.
Water evaporates from a wetted surface of rectangular shape, with an air
stream (1 atm, 35Β Β°C) blown over the surface at a velocity
of 15 m/s. Heat transfer measurements indicate that for air at 1 atm and
35Β Β°C, the average heat transfer coefficient is given
by the following empirical relation:
The transport properties are rounded model inputs. The surface saturation
pressure is consistent with 35 degrees Celsius using the
NIST Antoine correlation.
Assumptions
Steady transfer from a uniformly wet surface.
Dry incoming air.
The surface is maintained at the stated air temperature.
Constant supplied properties.
Heat and mass transfer share the same flow and surface conditions.
Normal evaporation flow is small.
Calculations
The measurement gives
Β―Β―Β―Β―β=21Γ150.6=106.63βWβmβ2βKβ1.
Compute Sc=0.61115.
The matching-flow Chilton-Colburn relation is
Β―Β―Β―Β―ππ=Β―Β―Β―Β―βππΆπ(PrSc)2/3
Substituting the coefficient and property values gives
Show how this equation can be used to predict the mass-transfer
coefficient from a naphthalene rod that is used to form the inner tube
of an annular duct. π1 and π2 are the diameters of the inner and
outer tubes, respectively, whereas π and ππ are the viscosity
of the fluid at the bulk and surface temperatures.
Solution
Data
Inner diameter π1
Outer diameter π2
Reynolds number based on inner diameter Reπ1
Empirical correlation constant πΆ
No numerical geometry, velocity or value of the correlation constant is supplied. Here the symbol for the correlation constant does not denote concentration.
Assumptions
The heat and mass problems have matching geometry and flow regime.
Uniform surface concentration corresponds to the heat correlationβs wall condition.
Low mass flux leaves the velocity field essentially unchanged.
The gas is dilute and isothermal, so wall and bulk viscosities are approximately equal.
The translated correlation is used within the applicable Reynolds and Schmidt ranges.
Calculations
Keep the same geometry, Reynolds definition and experimental range.
Replace Nusselt by Sherwood and Prandtl by Schmidt.
At isothermal dilute conditions ππ βπ, giving
Shπ1=πΆRe0.8π1Sc1/3(π2π1)0.15
ππ=π·π΄π΅π1Shπ1(17)
The viscosity factor can be retained if wall temperature or composition changes viscosity appreciably, but the heat correlation alone does not validate a mass-composition correction.
The Reynolds and diameter-ratio factors describe flow and geometry, so they remain unchanged. The viscosity correction becomes unity under the stated isothermal dilute approximation.
Dry air (310 K, 1.013Γ105Β Pa, 1.5 m/s) is flowing
through a pipe lined with a felt material that is continuously saturated
with water at 290 K. The pipe is 0.15 m in diameter and 6 m long. The
saturated felt material has a roughness index, π/π·, of 0.00001. Assuming
that the temperature of the air and pipe wall is constant, determine the
rate at which water must be added to keep the felt continuously saturated.
Solution
The inlet, outlet and wetted wall are shown in Figure 4. Concentration rises along the pipe, so the downstream driving force is smaller than the inlet driving force.
Data
Quantity
Value
Pipe diameter
π=0.15Β m
Pipe length
πΏ=6Β m
Bulk air speed
π=1.5Β mβsβ1
Dry inlet temperature
ππ=310Β K
Wet wall temperature
ππ =290Β K
Surface vapour pressure
ππ΄,π =1917Β Pa
Kinematic viscosity
π=1.569Γ10β5Β m2βsβ1
Diffusivity
π·π΄π΅=2.626Γ10β5Β m2βsβ1
Fanning friction factor
πΆπ=0.007
Water molar mass
ππ΄=0.018Β kgβmolβ1
The vapour pressure, transport properties and friction factor come from the lecturerβs notes.
Assumptions
Steady flow and evaporation.
Constant bulk velocity and supplied properties.
Uniform saturated wall concentration.
Each cross-section has a representative mixed bulk concentration.
Axial diffusion and normal-flow effects are neglected.
A developed-flow average coefficient approximates the pipe transfer.
The pipe spans 40 diameters; the roughly ten-diameter turbulent development estimate motivates a developed average approximation, not exact development at the inlet.
The concentration model treats the maintained wall and gas temperatures through the supplied surface concentration and constant properties. It does not solve the axial thermal field. The friction analogy gives the coefficient as
The source tableβs Dittus-Boelter-type mass row gives about 0.00688 m sβ»ΒΉ, but Sc=0.5975 falls below its stated 0.7 limit.
That comparison is an extrapolation and is not the governing answer.
The primary friction estimate also sits at the approximate lower edge of the sourceβs Chilton-Colburn range; treat the result as an engineering estimate.
For a slice ππ₯, steady balance is inflow minus outflow plus wall input equals zero.
Inflow is ππ΄πππ΄(π₯), outflow is ππ΄πππ΄(π₯+ππ₯), and wall input is
The required water replenishment rate is approximately 2.07Γ10β4Β kgβsβ1, or 0.745Β kgβhβ1. The exit bulk concentration is 0.434Β molβmβ3.
A 225 mm diameter pan of water is maintained at 37Β Β°C
is exposed to a dry, still atmosphere at 17Β Β°C. Estimate
the heat and mass transfer coefficients, as well as the evaporation rate
from the pan,
a. Assuming the convection is driven by species transfer;
b. Assuming that the convection is driven by temperature gradients.
Solution
Data
Quantity
Value
Pan diameter
π=0.225Β m
Rounded surface temperature
ππ =310Β K
Rounded air temperature
πβ=290Β K
Pressure
π=101325Β Pa
Gravitational acceleration
π=9.81Β mβsβ2
Kinematic viscosity
π=15.89Γ10β6Β m2βsβ1
Diffusivity
π·π΄π΅=2.6Γ10β5Β m2βsβ1
Prandtl number
Pr=0.707
Specific heat capacity
πΆπ=1006Β Jβkgβ1βKβ1
Air conductivity
π=0.0263Β Wβmβ1βKβ1
Film density
π=1.1769Β kgβmβ3
Expansion coefficient
π½=1/300Β Kβ1
Surface saturation pressure
πsatπ΄=6220Β Pa
Water molar mass
ππ΄=0.018Β kgβmolβ1
Air molar mass
ππ΅=0.029Β kgβmolβ1
Use the sourceβs rounded temperatures and the supplementary properties from the lecturerβs notes.
Assumptions
Steady natural convection above an upward-facing wet surface.
Dry bulk air and saturated gas at the water surface.
Uniform maintained surface temperature.
Constant supplied film properties and ideal gas behaviour.
Low normal evaporation flux.
Composition buoyancy and thermal buoyancy are considered separately, as requested.
Calculations
Area is
π΄π =π(0.225)24=0.039761βm2
and horizontal characteristic length is
πΏ=π΄π ππ =π4=0.05625βm.
Using diameter instead gives a Grashof number 64 times too large.
These are the two requested isolated-mechanism answers.
For comparison, using dry-air density at 290 K and saturated-air density at 310 K gives a total density difference that includes both mechanisms.
The mean density 1.16616 and total difference 0.10515 kg mβ»Β³ give Gr=6.2366Γ105, Sh=13.419 and ππ=0.006202 m sβ»ΒΉ.
The corresponding total-density estimate is Λππ΄=1.0712Γ10β5 kg sβ»ΒΉ, about 38.6 g hβ»ΒΉ.
This approximately 39 g hβ»ΒΉ estimate already includes thermal expansion. Adding a thermal Grashof contribution would count that expansion twice.
The real pan has reinforcing thermal and composition buoyancy.
The total-density route is a possible single-group engineering approximation; an exact coupled answer requires a validated heat-and-species natural-convection model.
The two isolated evaporation rates cannot be added to predict the combined rate, because the same velocity field transports both heat and species.
The saturation fraction is 0.0614, so low-flux and dry-air-property assumptions remain approximate.
The maintained water temperature implies external heat supply; no latent-heat constraint has been omitted from an adiabatic balance because this pan is not adiabatic.
Answer
For composition-driven convection alone, ππ=0.00443Β mβsβ1, ββ4.76Β Wβmβ2βKβ1 and Λππ΄β27.6Β gβhβ1.
For temperature-driven convection alone, β=6.03Β Wβmβ2βKβ1, ππ=0.00562Β mβsβ1 and Λππ΄β34.9Β gβhβ1.
A water tray is flush with the floor of a rectangular duct. The tray is
0.60 m long along the air flow and 0.25 m wide. The duct is 0.25 m wide
and 0.12 m high. Initially dry air enters at 2.0 m sβ»ΒΉ, 30βΒ°C
and 101325 Pa. A heater maintains the water at 30βΒ°C.
Assume a uniform air stream over the tray and a concentration boundary layer
beginning at its leading edge. Neglect upstream velocity-profile effects.
Use π=1.60Γ10β5 mΒ² sβ»ΒΉ, π·π΄π΅=2.60Γ10β5 mΒ² sβ»ΒΉ,
water saturation pressure 4240 Pa, ππ΄=0.0180 kg molβ»ΒΉ and
π =8.314 J molβ»ΒΉ Kβ»ΒΉ. Calculate the average mass transfer coefficient and
evaporation rate. Estimate the outlet relative humidity and use it to check
the assumption of dry bulk air over the entire tray. Use the dilute-transfer
flat-plate model without a Stefan-flow correction.
Solution
Data
Quantity
Value
Tray length
πΏ=0.60Β m
Tray width
π=0.25Β m
Duct height
π»=0.12Β m
Air speed
π=2.0Β mβsβ1
Temperature
π=303.15Β K
Pressure
π=101325Β Pa
Kinematic viscosity
π=1.60Γ10β5Β m2βsβ1
Diffusivity
π·π΄π΅=2.60Γ10β5Β m2βsβ1
Surface saturation pressure
πsatπ΄=4240Β Pa
Water molar mass
ππ΄=0.0180Β kgβmolβ1
Assumptions
Steady external transfer with a boundary layer starting at the tray leading edge.
Uniform incoming velocity and negligible upstream profile effects.
Maintained saturated surface at the air temperature.
Constant properties and dilute transfer without a Stefan-flow correction.
Dry bulk air is used for the initial rate estimate and then checked.
Calculations
The concentration boundary layer grows along the tray, so use its length
πΏ=0.60 m for external transfer. The duct hydraulic diameter is not the
length in this flat-plate correlation:
Under the assumed mixed-bulk description, the bulk concentration reaches
only 1.67 percent of the surface saturation concentration. The driving force
therefore falls by less than 2 percent from inlet to outlet. This supports
using dry bulk air in the first evaporation estimate. A much wetter outlet
would require coupling local transfer to the changing bulk humidity.
This humidity check does not remove all model approximations. The surface
mole fraction is
4240101325=0.04185,
so a Stefan-flow correction can
also change the dilute estimate by a few percent. The imposed heater avoids
an evaporative-cooling calculation; without it, the surface temperature and
vapour pressure would need a coupled energy balance.
Answer
The average mass transfer coefficient is Β―Β―Β―Β―ππ=0.00670Β mβsβ1 and the evaporation rate is 0.110Β kgβhβ1. The estimated outlet relative humidity is 1.67 percent, supporting the dry-bulk approximation.
A rigid spherical pellet of diameter 0.010 m has a thin wet coating containing
0.0020 g of volatile solvent per cmΒ² of pellet surface. Solvent-free air flows
past it at 0.10 m sβ»ΒΉ. Both the air and pellet surface are maintained at 300 K
and 101325 Pa. External temperature control supplies the evaporation heat.
Use air properties π=1.18 kg mβ»Β³, π=1.85Γ10β5 Pa s,
π=0.0263 W mβ»ΒΉ Kβ»ΒΉ and Pr=0.707. The solvent has
π·π΄π΅=9.62Γ10β6 mΒ² sβ»ΒΉ, saturation pressure 1270 Pa at 300 K,
molar mass 0.0780 kg molβ»ΒΉ and liquid density 876.5 kg mβ»Β³.
For this low-Re sphere, adopt the matching correlations
Estimate the heat and mass transfer coefficients and the drying time.
Assume the exposed area stays saturated at the stated vapour pressure until
the solvent is exhausted, coating resistance is negligible, and normal
evaporation flow is small. Check the thin-coating assumption.
SolutionFigure 5: A thin solvent coating surrounds the pellet. Its surface concentration is fixed by vapour equilibrium, and solvent transfers into the surrounding air. The diagram shows the external gas film; the model tracks coating loss rather than diffusion inside the pellet.
Data
Quantity
Value
Pellet diameter
π·=0.010Β m
Air speed
π’=0.10Β mβsβ1
Initial solvent loading
ππ΄/π΄π =0.0200Β kgβmβ2
Solvent diffusivity
π·π΄π΅=9.62Γ10β6Β m2βsβ1
Solvent saturation pressure
πsatπ΄=1270Β Pa
Solvent molar mass
ππ΄=0.0780Β kgβmolβ1
Liquid density
ππΏ=876.5Β kgβmβ3
Air and surface temperature
π=ππ =300Β K
Use the supplied air properties.
Assumptions
Quasi steady external transfer and constant supplied properties.
Negligible solvent in the bulk air.
The exposed coating remains saturated until the solvent is exhausted.
Negligible coating resistance and normal evaporation flow.
Surface temperature is maintained by external heat supply.
Fixed pellet diameter is used and the coating thickness is checked.
Calculations
Use pellet diameter for both the flow Reynolds number and the transfer groups:
This Reynolds number is within the low-Re use described for the sphere
correlation in the external forced-convection notes. The two transfer groups
differ because Schmidt and Prandtl numbers differ. Their additive two is
the stagnant-sphere limit and must be retained.
Recover the dimensional coefficients with the corresponding transport
property:
The loading converts as 0.0020 g cmβ»Β² =0.0200 kg mβ»Β², so the
initial solvent mass is ππ΄0=6.283Γ10β6 kg. The constant-rate
estimate is
π‘πππ¦=ππ΄0ππ΄ππ΄π΄π =0.02000.0780(0.003742)=68.5Β s
The area cancels because both inventory and evaporation rate scale with it.
Diameter still affects drying time through the transfer coefficient.
The solvent-equivalent coating thickness is
0.0200876.5=2.28Γ10β5Β m,
only 0.46 percent of the pellet
radius. Treating the diameter as fixed is consistent with a thin coating
provided the nonvolatile coating is also thin.
The surface mole fraction is 0.01253, supporting the dilute approximation.
Actual paint can dry more slowly as solvent activity falls or diffusion
through the coating becomes limiting; the calculated time uses the stated
constant-activity model.
Answer
The heat transfer coefficient is β=16.5Β Wβmβ2βKβ1, the mass transfer coefficient is ππ=0.00735Β mβsβ1 and the constant-rate drying time is approximately 68.5 s. The solvent-equivalent thickness is 22.8Β πm, only 0.46 percent of the pellet radius.
Saturated steam at 130ββC is flowing inside a steel pipe having an ID
of 2.09 cm and an OD of 2.67 cm. The pipe is insulated with 3.8 cm of insulation on the
outside. The convective coefficient for the inside steam surface of the pipe is
estimated to be βπ = 5680 W/m2.K, while on the outside of the insulation
it is estimated to be βπ = 11.4 W/m2.K. The mean thermal conductivities
of the metal and the insulation material are 45 W/m.K and 0.064
W/m.K, respectively. Calculate the heat loss for 1 m of pipe using the method
of resistances if the surrounding air is at 26ββC.
Worked solution
Heat flows radially from the steam, through the steel and insulation, and
then into the air. The same heat rate passes through four resistances.
Figure 1: Cross-section of an insulated pipe. The inner and outer steel radii are r1 and r2; insulation extends to r3. Heat crosses the inner film, steel, insulation and outer film in series.
Data
Quantity
Symbol
Value
Pipe length
πΏ
1.0Β m
Inside diameter
π·1
0.0209Β m
Outside steel diameter
π·2
0.0267Β m
Radial insulation thickness
π‘
0.0380Β m
Steel conductivity
ππ
45Β Wβmβ1βKβ1
Insulation conductivity
ππ
0.064Β Wβmβ1βKβ1
Steam and air temperatures
ππ π‘πππ,ππππ
130,26ββC
Inside and outside coefficients
βπ,βπ
5680,11.4Β Wβmβ2βKβ1
We need the heat loss from the specified one-metre length.
Assumptions
Steady radial heat flow with negligible end losses.
Constant conductivities and heat transfer coefficients.
Perfect thermal contact between steel and insulation.
No separate radiation or fouling resistance.
Calculations
First convert diameters to radii. Insulation thickness is already a radial
length, so it is added directly to the outside steel radius:
π1=π·12=0.01045Β m,π2=π·22=0.01335Β mπ3=π2+π‘=0.01335+0.0380=0.05135Β m
Each film uses its own cylindrical area, 2πππΏ. Conduction through a
cylindrical layer uses both bounding radii. The four resistances are
Element
Resistance relation
Value, KβWβ1
Inside steam film
π ππ=1/(βπ2ππ1πΏ)
0.00268
Steel wall
π π =lnβ‘(π2/π1)/(2πππ πΏ)
0.000866
Insulation
π π=lnβ‘(π3/π2)/(2ππππΏ)
3.35
Outside air film
π ππ=1/(βπ2ππ3πΏ)
0.272
For example, the largest resistance is obtained from
It lies above the air temperature and below the steam temperature. The
insulation supplies about 92 percent of the total resistance. Its
large temperature drop is consistent with the small heat loss.
Answer
The heat loss is approximately 28.7Β W from one metre of pipe,
or 28.7Β Wβmβ1. This is a heat rate per pipe length,
not a heat flux per surface area.
In a laboratory experiment, an electric current is passed through a slab of thickness
2πΏ. As a result, heat generation ππ takes place in the slab. The slab is
insulated on both sides (thickness = t) to satisfy the safety considerations related to
the danger stemming from the fact that an electric current is passing. In addition, the
insulated slab is cooled on both sides by a fluid of temperature ππ, with convective
coefficient of h. Determine the temperature at the surface of the insulation, assuming
that the thermal conductivities for the slab and the insulation material are constants,
denoted by ππ and ππ, respectively.
Worked solution
Choose π₯=0 at the slab centre and π₯=πΏ at its right face.
Symmetry gives ππ/ππ₯=0 at the centre.
The slab energy balance is
ππ π2πππ₯2+ππ=0.
Integrating with the symmetry condition gives
ππππ₯=βπππ₯ππ .
The heat flux leaving either slab face, in Wβmβ2, is therefore
A pipe designed to transport liquid metal is to be embedded in a 1.2 m thick wall. The
inside surface of the wall has a constant surface temperature of 925 K, whereas its
external surface is exposed to air with a convective heat transfer coefficient of 23
W/m2.K and 300 K. The thermal conductivity of the wall varies with
temperature according to π=0.073+3.942Γ10β4π, where T is temperature in
Kelvin, and k is in W/m.K. It is desired that the pipe is located at a point
in the wall where the temperature is < 650 K. Calculate how far from the hot surface
should the pipe be located and the heat flux for the wall.
Worked solution
Data
Quantity
Symbol
Value
Wall thickness
πΏ
1.2Β m
Hot-face temperature
πβ
925Β K
Outside air
πβ
300Β K
Outside coefficient
β
23Β Wβmβ2βKβ1
Conductivity
π(π)
0.073+3.942Γ10β4πΒ Wβmβ1βKβ1
Pipe temperature limit
ππππππ‘
650Β K
We need the outward heat flux and the depth at which the wall reaches the pipe-temperature limit.
Assumptions
Steady one-dimensional conduction with no heat generation.
Negligible disturbance of the wall temperature field by the embedded pipe.
The given conductivity relation applies throughout the temperature range.
Calculations
Write π(π)=π+ππ, where π=0.073βWβmβ1βKβ1 and π=3.942Γ10β4βWβmβ1βKβ2.
Let π₯ increase from the hot face towards the air, πΏ=1.2βm, and ππ=π(πΏ).
With no generation the positive heat flux is constant:
πβ²β²=βπ(π)ππππ₯.
Integrating across the wall and applying convection gives
In a production facility, 4-cm-thick large plates (k = 107 W/mβ βC, π = 8520 kg/m3, πΆπ = 381.2 J/kgβ βC) that
are initially at a uniform temperature of 20ββC are heated by passing
them through an oven maintained at 500ββC. The plates remain in the
oven for a period of 7 min. Taking the convection heat transfer coefficient to be h =
120 W/m2β βC, determine the surface temperature of the plates
when they come out of the oven. Assume L = W = 1 m.
Worked solutionFigure 1: Both faces of the plate are heated by the oven air. Use half the plate thickness for the plane-wall solution, not its one-metre plan dimension.
Data
Quantity
Symbol
Value
Plate dimensions
π»,π,2πΏ
1.0,1.0,0.040Β m
Initial and oven temperatures
ππ,πβ
20,500ββC
Heating duration
π‘
420Β s
Convection coefficient
β
120Β Wβmβ2βKβ1
Conductivity
π
107Β Wβmβ1βKβ1
Density
π
8520Β kgβmβ3
Specific heat
ππ
381.2Β Jβkgβ1βKβ1
Here πΏ=0.020Β m denotes half the plate thickness in the
plane-wall model. The one-metre dimensions in the question set the face area.
Assumptions
Both broad faces see the same oven temperature and coefficient.
Edge heat transfer and radiation are neglected in the plane-wall comparison.
Properties remain constant and the initial temperature is uniform.
Calculations
The worked solution compares lumped capacitance with the one-term plane-wall solution.
For the comparison, solve π1tanβ‘π1=π΅π with
π΅π=βπΏ/π and use
The finite-wall result is about 0.38 degrees above the lumped estimate.
The large Fourier number makes higher eigenmodes negligible.
The temperature is not exactly uniform, but the small Biot number explains the close agreement.
Answer
The one-term plane-wall estimate gives a surface temperature of 279.4ββC. The simpler lumped estimate is 279.1ββC, consistent with the small Biot number.
During the quenching process of stainless steel ball bearings, it is known that the
temperature of the bearings should not reach below 850ββC prior to
entering the water to minimise the developments of surface defects. Ball bearings with
diameter of 1.2 cm (π = 8085 kg/m3, πΆπ = 0.480
kJ/kgβ βC, k = 15.1 W/mβ βC and πΌ =
3.91 Γ 10β6m2/s) are taken out of the oven at a uniform temperature of
900ββC, and are exposed to 30ββC air for a while
before they are dropped into the water. Assuming that the heat transfer coefficient of
the air is 125 W/m2.K, how long can the spheres be left in air before
quenching?
Worked solutionFigure 2: A small spherical bearing cools in air. Its radius is half the supplied diameter, and the exposed area is the full spherical surface.
For a sphere the volume-to-area length is one third of its radius.
With π =0.006βm, πΏπ=π /3=0.002βm and
π΅ππ=βπΏππ=0.01656<0.1.
Use the lumped cooling relation to solve for time:
In the Betty Crocker Cookbook, it is stated that it takes 2 h 45 minutes to roast a 3.2
kg ribs that is initially at 4.5ββC in an oven that is maintained at
163ββC. It is recommended that a meat thermometer be used to monitor
the cooking, and the rib is considered medium-rare when the thermometer, inserted into
the centre of the thickest part of the meat, registers 60ββC. If the
ribs is approximated as a homogeneous spherical object with π = 1200
kg/m3, πΆπ = 4.1 kJ/kgβ βC, k = 0.45
W/mβ βC and πΌ = 0.91 Γ 10β7m2/s,
determine:
a. The heat transfer coefficient at the surface of the ribs
b. The temperature at the outer surface of the ribs when it is done (medium rare)
Worked solutionFigure 3: A spherical body has distinct centre and surface temperatures during heating. The measured centre temperature is used to infer the coefficient at the outer boundary.
Data
Quantity
Symbol
Value
Mass
π
3.2Β kg
Density
π
1200Β kgβmβ3
Initial and oven temperatures
ππ,πβ
4.5,163ββC
Observed centre temperature
ππ
60ββC
Elapsed time
π‘
9900Β s
Conductivity
π
0.45Β Wβmβ1βKβ1
Supplied diffusivity
πΌ
0.91Γ10β7Β m2βsβ1
Assumptions
Use the given rounded diffusivity for the transient calculation. Assume radial
conduction, constant properties, no evaporation, no generation, and no separate
radiative heat flux. The effective coefficient describes the stated idealised
boundary, not a measured material property.
Calculations
The first step is to infer the radius from
ππ=4ππ 33.
For π΅ππ =βπ /π, the positive eigenvalues satisfy
1βππcotβ‘ππ=π΅ππ . The full spherical solution is
Use ππ(0)=1. Retain enough terms for convergence; the usual one-term
criterion πΉπ>0.2 may not hold here. The first eigenvalue lies between zero
and π, with each subsequent one between adjacent multiples of π.
The mass and density fix the radius without a geometric guess.
π=3.21200=0.0026667βm3
and π =0.0860254βm.
With π‘=9900βs, πΉπ=0.121737.
The observed centre temperature requires
Ξ0=60β1634.5β163=0.649842.
For a trial π΅ππ , solve the eigenvalue equation in each interval and evaluate Eq. 1 at the centre.
Adjust π΅ππ until the sum equals 0.649842.
The following trials show how the inverse calculation brackets the measured
centre temperature. Each row uses its own eigenvalues in the same series.
Trial π΅ππ
First eigenvalue π1
Centre ratio Ξ0
Predicted centre temperature
10.0
2.83630
0.690587
53.54ββC
20.0
2.98572
0.643649
60.98ββC
17.8639
2.96741
0.649842
60.00ββC
The first trial heats the centre too slowly; the second heats it too quickly.
Solving between them gives the fitted Biot number. Once it is fitted, use
that same boundary coefficient to predict the surface temperature.
Five terms already agree with ten terms to better than 10β10 in the temperature ratio at the fitted coefficient.
The surface is much hotter than the centre during heating.
Lumped capacitance is unsuitable because π΅ππ=π΅ππ /3=5.95.
A one-term inverse calculation also introduces avoidable error at this Fourier number.
The handwritten sourceβs trial value π΅ππ β30 does not reproduce the stated centre temperature with a converged series; its coefficient is not used as the answer.
Answer
The fitted effective coefficient is β=93.5Β Wβmβ2βKβ1 and the surface temperature is 156.5ββC. These values use the converged spherical series and the stipulated constant-property model.
Water at 16ββC enters a 2.54 cm ID tube that is used to cool a
nuclear reactor. The water flow rate is 0.114 m3/min. Using the Colburn
analogy, determine the total rate of heat transfer, the exiting water temperature and
the wall temperature at the exit, if the length of the tube is 4.5 m and the wall
condition is one of uniform heat flux of 5.7 Γ 103kJ/hr.m2.
Worked solutionFigure 1: Water flows through the heated tube. The specified inward heat flux supplies the bulk enthalpy rise; the local outlet coefficient determines the exit wall temperature.
Data
Quantity
Symbol
Value
Tube diameter and length
π·,πΏ
0.0254,4.5Β m
Volume flow
Λπ
0.114Β m3βminβ1
Inlet temperature
πππ
16ββC
Uniform wall flux
πβ²β²
5.7Γ103Β kJβhβ1βmβ2
Water property
Constant value
Density π
1000βkgβmβ3
Specific heat ππ
4180βJβkgβ1βKβ1
Dynamic viscosity π
0.001βPaβs
Prandtl number ππ
7.0
Assumptions
Use the source solutionβs rounded water properties throughout. Assume steady
single-phase flow, fully developed turbulent heat transfer, negligible axial
conduction and no heat loss elsewhere. Ignore entrance effects. The prescribed
heat flux acts on the inner wetted surface.
Here ππΉ is the Fanning friction factor, one quarter of the Darcy factor.
The correlation applies to smooth turbulent flow in the approximate range
4000<π π<105. Use
Consider the flow of oil in a 30 cm diameter pipeline at an average velocity of 2 m/s
and temperature of 20ββC. A section of the pipeline (200 m long)
passes through icy waters of a lake at 0ββC. Measurements indicate
that the surface temperature of the pipe is very nearly 0ββC.
Disregarding the thermal resistance of the pipe material, determine:
a. The temperature of the oil when the pipe leaves the lake;
Assume steady hydrodynamically developed flow at the start of the cold section,
a uniform inlet bulk temperature, constant properties and negligible viscous
heating and axial conduction. Neglect pipe-wall resistance, so the inner wall
is also isothermal at ππ =0ββC. The thermal boundary layer
starts at the lake entrance. Check π π=π’π·/π and ππ=ππππ/π.
Calculations
For laminar flow, estimate the thermal entry length with
The correlation models thermal development with a developed velocity profile;
ignore wall-to-bulk viscosity corrections under the stated constant-property
assumption. Do not substitute the fully developed value ππ’=3.66 prematurely.
The oil velocity profile is laminar but its thermal boundary layer is still developing.
The supplied properties give π π=665.93, ππ=10373.55 and πΏπ‘=1.1916Γ105βm, far longer than the cold section.
The large Prandtl number explains the long thermal entry region.
Equation Eq. 2 gives πΊπ§=10362.04, Β―Β―Β―Β―Β―Β―Β―Β―ππ’π·=37.3168 and Β―Β―Β―Β―β=18.0364βWβmβ2βKβ1.
The mass flow rate is
Λπ=ππ’ππ·24=125.538βkgβsβ1.
The dimensionless exponent compares the wall conductance with the heat-capacity flow rate:
The oil cools only 0.286 degrees because the heat-capacity flow rate is large.
This positive heat-loss convention is the magnitude of the negative heat input to the oil.
Retaining intermediate precision avoids the handwritten sourceβs rounded outlet temperature changing the reported duty to about 68.4 kW.
Answer
The oil exits at 19.7ββC and loses 67.5Β kW over the cold section. The unrounded energy balance uses a bulk temperature drop of 0.286 K. The flow is laminar but thermally developing.
A thick-walled steel pipe (k = 60 W/m.K) carrying hot water is cooled
externally by a cross-flow airstream at a velocity of 20 m/s and a temperature of
25ββC. The inner and outer diameters of the pipe are 20 mm and 25 mm,
respectively. At a certain location along the pipe, the mean temperature of the water is
80ββC. Assuming the flow inside the tube is fully developed with a
Reynolds number of 20,000, find the heat transfer rate to the air stream per unit pipe
length. What will be the heat transfer rate if the surrounding air is stagnant?
Worked solutionFigure 2: Heat passes from the water through the inside film, steel wall and outside air film. The inside and outside coefficients act on different cylindrical areas.
Data
Quantity
Symbol
Value
Inside and outside diameters
π·π,π·π
0.020,0.025Β m
Steel conductivity
ππ π‘πππ
60Β Wβmβ1βKβ1
Water and ambient air temperatures
ππ,πβ
80,25ββC
Internal Reynolds number
Reπ
20000
Cross-flow speed
π’β
20Β mβsβ1
The property approximations are
Property
Water
Moving air
Stagnant air film
Conductivity π, Wβmβ1βKβ1
0.67
0.0263
0.0282
Prandtl number ππ
2.2
0.707
0.707
Kinematic viscosity π, m2βsβ1
Not needed
15.89Γ10β6
18.2Γ10β6
Assumptions
Steady radial wall conduction with negligible fouling and axial wall conduction.
An isolated long horizontal pipe with radiation excluded.
The local bulk water temperature is fixed in both air cases.
Water and moving-air values follow the source solution. The stagnant-air
values approximate air at 325 K and stay fixed while the film temperature
in π½=1/ππ is iterated. Use π=9.81βmβsβ2.
Calculations
For internal water use the sourceβs Colburn form; for cylinder cross-flow use
the stated Zukauskas approximation with the surface-Prandtl correction set to one:
The outside correlation applies to 1000<π ππ<2Γ105 and πππ<10.
This explicitly uses the sourceβs 1/3 internal exponent rather than the
Dittus-Boelter cooling exponent. Resistances per unit length have units
mβKβWβ1:
The handwritten resistance labels transpose the inner and outer coefficients; their sum and approximate 489.5 W per metre result are nevertheless consistent.
Here each resistance uses its correct coefficient and diameter.
For stagnant air, trial an outer surface temperature between 25 and 80 degrees and equate the conductive/internal flux to the natural-convection flux.
The natural-convection Rayleigh number lies within the specified range.
The final film is close to the supplied 325 K property approximation.
Stagnant air removes much less heat and leaves the wall closer to the water temperature.
Radiation was excluded, so this is a convection-only comparison.
Answer
The heat loss is approximately 490Β Wβmβ1 in moving air and 32.3Β Wβmβ1 in stagnant air. The latter is the convection-only result for the specified horizontal pipe.
A tank (diameter = 1m, height = 3 m) containing liquid methanol (MW = 32 g/mol) has its
top open to the atmosphere. The level of methanol in the tank is maintained at 0.5 m.
The gas space inside the tank is stagnant and the methanol vapours are immediately
dispersed once they exit the tank. Engineers onβsite are concerned that this open tank
may be emitting a considerable amount of methanol vapour.
a. What is the emission rate of methanol vapour from the tank (in kgΒ methanol/day) when the temperature of the tank is 30ββC? State all
assumptions and boundary conditions.
b. If the temperature of the tank is raised to 40ββC, what is the new
emission rate?
π·π΄π΅ (methanol/air) @ 298 K = 0.162 Γ 10β4m2/s.
Worked solutionFigure 1: The open tank has the same one-dimensional gas-space model as this evaporation tube. The diffusion distance is the tank height minus the maintained liquid depth.
Data
Quantity
Symbol
Value
Tank diameter
π
1.0Β m
Gas depth above maintained liquid level
πΏ
3.0β0.5=2.50Β m
Atmospheric pressure
π
101325Β Pa=760Β mmHg
Methanol molar mass
ππ΄
0.032Β kgβmolβ1
Reference diffusivity
π·π΄π΅(298)
1.62Γ10β5Β m2βsβ1
Saturation pressure at 30 and 40 degrees Celsius
ππ ππ‘
163,265Β mmHg
We need the total daily emission from the open tank at each temperature.
Assumptions
Steady one-dimensional diffusion through the gas space.
Uniform temperature in each calculation and ideal-gas behaviour.
Constant pressure, maintained liquid level and saturated vapour at the surface.
Zero net air flux through the liquid; dry gas at the tank mouth.
Natural convection, wind penetration and transport through the walls are neglected.
Calculations
Let A denote methanol and B denote air. Coordinate π points upwards from
the liquid. Zero net air flux permits Stefan flow of the mixture, so use
Here π =8.314Β Paβm3βmolβ1βKβ1. The boundary values are
π¦π΄(0)=ππ ππ‘/π and π¦π΄(πΏ)=0. At the same pressure, estimate the
change in diffusivity from
π·π΄π΅(π)=π·π΄π΅(298)(π298)3/2
Integrate Eq. 1 across the gas space and multiply the flux by the
open cross-sectional area, molar mass and seconds per day.
The maintained liquid level fixes the gas diffusion distance throughout each calculation.
Since ππ΄ is constant over the planar gas space, integration gives
Consider the evaporation of water droplets (diameter ~ 0.1 mm) in still air at
50ββC, 1 atm pressure, with 20% relative humidity. As the droplets
are very small, it can be considered that their terminal velocity is negligible, i.e.
they are suspended in air. The droplet temperature can be assumed to be at
50ββC and its vapour pressure is 12 kPa. If the diffusion coefficient
of water vapour in air is 0.288 cm2/s,
a. Calculate the initial rate of evaporation of a single droplet, assuming steady state
condition.
b. Determine the time required for the water droplet to evaporate completely.
Worked solutionFigure 2: The droplet evaporates into air with a fixed far-field humidity. Its surface concentration is saturated, while its radius decreases as liquid is lost. The dashed outer boundary is taken to infinity in this problem.
Data
Property
Value
Total pressure π
101325βPa
Temperature π
323.15βK
Water saturation vapour pressure ππ ππ‘
12000βPa
Water-air diffusivity π·π΄π΅
0.288βcm2βsβ1
Liquid water density ππΏ
1000βkgβmβ3
Water molar mass ππ΄
0.018βkgβmolβ1
Assumptions
Assume spherical symmetry, negligible droplet velocity, constant properties
and no air dissolving in the droplet. The gas profile adjusts much faster
than the droplet shrinks, so treat it as steady at each instantaneous radius
π(π‘). An external heat supply maintains the liquid and gas at 50 degrees;
ignore evaporative cooling. Far-field humidity stays fixed in an effectively
infinite gas domain. At π=π, π¦π =ππ ππ‘/π. At πββ,
π¦β=0.20ππ ππ‘/π. Relative humidity is a ratio of water partial
pressure to saturation pressure, not a gas mole fraction.
Calculations
Derive π2ππ΄=constant from a spherical gas shell balance. Integrate
Eq. 1 radially using ππ΅=0. Then use the shrinking liquid balance
where Λππ΄=4ππ2ππ΄ is the molar evaporation rate.
Use the instantaneous rate as the radius decreases; do not divide the initial
mass by the initial rate to obtain the lifetime.
The spherical area increases outward while the total molar evaporation rate stays constant.
A steady shell balance gives
π(π2ππ΄)ππ=0.
Write π2ππ΄=πΆ and integrate the flux law between π and infinity:
Dividing it by the initial mass rate gives 0.4175 s and underestimates the lifetime because the rate decreases with radius.
The gas diffusion time
π20π·π΄π΅=8.68Γ10β5βs
is much shorter than π‘π, which supports the quasi-steady approximation.
Complete evaporation is the πβ0 limit of this continuum model.
Cooling and molecular-scale effects are outside the assumptions.
Answer
The initial mass evaporation rate is 1.25Γ10β9Β kgβsβ1 per droplet. The shrinking-sphere model gives a lifetime of 0.626Β s, with the droplet temperature maintained at 50ββC.
A cylindrical graphite (pure carbon) rod (length = 25 cm, initial diameter = 2 cm) is
placed in a combustion chamber with flowing air stream at 1100 K and 2 atm total
pressure. A stagnant gas boundary layer 5 mm thick is formed around the surface of the
rod. At 1100 K, solid carbon oxidises to carbon dioxide gas,
CΒ (s)+O2Β (g)βCO2Β (g)
Assuming that the oxidation is limited by molecular diffusion of O2 through
the stagnant gas film, and outside of the gas film, the bulk composition of air stream
prevails,
a) Estimate the initial rate of CO2 production from the rod;
b) Write an equation for the rate of change in the radius of the rod as a function of
time.
Data: Density of solid graphite = 2.25 g/cm3; π·π΄π΅ (CO2-air)
@ 273 K and 1 atm = 0.136 cm2/s
Worked solutionFigure 1: Radial diffusion through a gas film around the graphite rod. The radius a decreases during oxidation, while the film thickness remains fixed. The drawing is a cylindrical cross-section; oxygen diffuses inward and carbon dioxide outward.
Data
Quantity
Symbol
Value
Initial radius and rod length
π0,β
0.010,0.25Β m
Gas-film thickness
πΏ
0.0050Β m
Gas temperature and pressure
π,π
1100Β K,2Β atm
Property
Value
Graphite density ππΆ
2250βkgβmβ3
Carbon molar mass ππΆ
0.012βkgβmolβ1
Carbon-dioxide molar mass ππΆπ2
0.044βkgβmolβ1
Bulk oxygen mole fraction π¦π
0.21
Reference effective gas diffusivity π·0
0.136βcm2βsβ1 at 273 K and 1 atm
Assumptions
Treat the sourceβs CO2-air diffusivity as a common effective diffusivity for
the oxygen and carbon-dioxide counterdiffusion. This is an approximation to
multicomponent air. Nitrogen has zero net flux, and the reaction exchanges
one gas mole for one gas mole, so total gas molar flux is zero.
Estimate the diffusivity by
π·=π·01βatmπ(π273βK)3/2.
Assume isothermal quasi-steady radial diffusion, no gas-phase reaction,
constant density and length, and no transport through the ends. The maintained
bulk gas has π¦π2(π+πΏ)=π¦π, while the fast reaction gives
π¦π2(π)=0. Keep πΏ fixed as π decreases. Derive the constant
πππ2 from a cylindrical shell balance, integrate Eq. 1, and
relate oxygen consumption to the changing carbon mass. Do not use the initial
area as a constant area for the shrinking rod.
Calculations
Use molar flux ππ΄ in molβmβ2βsβ1 and ideal-gas
concentration π=π/(π π), with
π =8.314βPaβm3βmolβ1βKβ1.
With negligible total molar flux, diffusion obeys
ππ΄=βππ·ππ¦π΄ππ .(1)
The coordinate π increases away from the reacting surface. Reactant flux
is therefore negative. Report removal and production rates as positive magnitudes.
The oxygen consumption balances carbon-dioxide production mole for mole at the rod surface.
A steady cylindrical shell balance gives
π(πππ2)ππ=0.
The positive oxygen consumption rate π satisfies π=β2ππβππ2.
Integrating Eq. 1 between π and π+πΏ gives
Carbon consumption uses ππΆ; carbon-dioxide mass production uses ππΆπ2.
The gas-film diffusion time πΏ2/π·=0.455βs is much shorter than the initial shrinkage scale
π0β£ππππ‘β£=2.97Γ104βs,
supporting quasi-steady transport initially.
The model excludes combustion heating, carbon-monoxide formation, end effects and changes in the flow as the rod shrinks.
Answer
The initial carbon-dioxide production rate is 9.91Γ10β4Β molβsβ1, or 4.36Γ10β5Β kgβsβ1. Radius evolves according to Eq. 3, starting from π0=0.010Β m; initially
Nitric oxide (NO) emission from automobile exhaust can be reduced by using a catalytic
converter. The following reaction occurs at the catalytic surface:
NO+COβ0.5Β N2+CO2
The concentration of NO is reduced by passing the exhaust gases over the surface of the
catalyst, and the rate of reduction at the catalyst surface is as follows:
πππ=βπβ²β²πΆππ(0),
where πΆππ(0) is the concentration of NO at the catalyst surface.
As a first approximation, it may be assumed that NO reaches the surface by
one-dimensional diffusion through a thin gas film of thickness L that adjoins the
surface, and that the effects of bulk motion is negligible. Determine the mole fraction
of NO at the catalytic surface and the NO removal rate for a catalytic surface area of
200 cm2.
Data: The exhaust gas is at 500ββC and 1.2 bar, the mole fraction of
NO at L is 0.15, π·π΄π΅ = 10β4m2/s, πβ²β² = 0.05 m/s and the film
thickness is L = 1 mm.
Worked solutionFigure 2: Nitric oxide diffuses through a nonreacting gas film to the catalyst. Its surface concentration is not zero when surface reaction is finite. The solution measures x from the catalyst toward the bulk.
Data
Property
Value
Bulk NO mole fraction π¦π
0.15
Film thickness πΏ
0.0010βm
Effective NO diffusivity π·
1.0Γ10β4βm2βsβ1
Surface rate constant ππ
0.050βmβsβ1
Catalytic area π΄π
200βcm2
NO molar mass πππ
0.030βkgβmolβ1
Assumptions
Assume steady isothermal one-dimensional diffusion with constant pressure and
π·, no gas-phase reaction, sufficient CO, and the sourceβs first approximation
of negligible bulk-motion contribution to NO transport. The reaction changes
total gas mole number, so this approximation is not exact equimolar diffusion.
Do not interpret the 0.15 mole fraction as realistic dilute exhaust data.
It is the supplied idealised problem value. At the surface,
The film resistance πΏ/π·=10βsβmβ1 and reaction resistance 1/ππ =20βsβmβ1 add.
Reaction supplies twice the film resistance, but neither can be ignored.
The calculated surface fraction in Eq. 4 confirms finite reaction, rather than a perfectly absorbing surface.
These results apply only to the stated film approximation.
Answer
The surface NO mole fraction is 0.100. The positive removal rate is 1.87Γ10β3Β molβsβ1, equivalent to 5.60Γ10β5Β kgβsβ1, within the stipulated film approximation.
The surface of a mild steel component is commonly hardened by packing the component in a
carbonaceous material in a furnace at a high temperature for a predetermined time.
Consider such a component with a uniform initial carbon concentration of 0.15% by mass.
The component is packed in a carbonaceous material and is placed in a highβtemperature
furnace. The diffusion coefficient of carbon in steel at the furnace temperature is
given to be 4.8 Γ 10β10m2/s, and the equilibrium concentration of
carbon in the iron at the interface is determined from equilibrium data to be 1.2% by
mass. Determine how long the component should be kept in the furnace for the mass
concentration of carbon 0.5 mm below the surface to reach 1%.
Worked solutionFigure 1: Carbon enters the steel through a maintained surface concentration. The target is at a depth measured inward from that surface; the far boundary remains unaffected.
Data
Quantity
Symbol
Value
Initial, surface and target mass fractions
π€π,π€π ,π€
0.0015,0.012,0.010
Target depth
π₯
0.00050Β m
Diffusivity at furnace temperature
π·
4.8Γ10β10Β m2βsβ1
Assumptions
Assume constant diffusivity, no bulk motion or reaction inside each solid,
uniform initial concentration and a maintained surface concentration.
The surface condition changes at π‘=0. For constant density, carbon mass
fraction obeys the same diffusion equation as concentration.
Calculations
Define the
normalised concentration
Ξ=πΆβπΆπ πΆπβπΆπ .(1)
Use the same ratio with mass fraction π€ replacing πΆ in problem 1.
The error function is
For interpolation, erfβ‘(0.170)=0.189992 and
erfβ‘(0.171)=0.191088. A calculator with an inverse error
function can also be used.
Treat the solid as a semi-infinite medium with a planar surface at π₯=0,
where π₯ increases into the steel. Other surfaces are far beyond the diffusion
penetration depth during the treatment. Assume no external transport resistance
and constant steel density. Use fractions 0.0015, 0.012 and 0.010, or use
percentages consistently. The semi-infinite solution is
Interpolation of the supplied table gives about 4482βs.
The sourceβs rounded π§=0.17 gives about 4505βs.
The characteristic penetration distance 2βπ·π‘=2.93βmm must be small relative to the component dimensions for this planar semi-infinite model to apply.
The answer excludes furnace heating time.
Answer
The required diffusion time at furnace temperature is 4.48Γ103Β s, approximately 74.7Β min. Furnace warm-up time is additional.
A large sheet of material 40 mm thick contains dissolved hydrogen at a uniform
concentration of 3 kmol/m3. The sheet is exposed to a fluid stream, causing
the concentration of the dissolved hydrogen to be reduced suddenly to zero at both
surfaces. This surface condition is maintained constant thereafter. If the mass
diffusivity of hydrogen is 9 Γ 10β7m2/s, how much time is required to
bring the density of dissolved hydrogen to 1.2 kg/m3 at the centre of the
sheet?
Worked solutionFigure 2: Hydrogen leaves both faces of the sheet under equal surface conditions. Use half-thickness L between the centre plane and either face.
Data
Quantity
Symbol
Value
Sheet thickness and half-thickness
2πΏ,πΏ
0.040,0.020Β m
Initial molar concentration
πΆπ
3.0Β kmolβmβ3
Target centre mass concentration
ππ»2,0
1.2Β kgβmβ3
Surface concentration
πΆπ
0
Diffusivity
π·
9.0Γ10β7Β m2βsβ1
Molecular hydrogen molar mass
ππ»2
2.0Β kgβkmolβ1
The worked solution uses molecular hydrogen molar mass ππ»2=2.0βkgβkmolβ1.
Assumptions
Assume one-dimensional diffusion normal to the faces, negligible edge effects
and no external transport resistance. Set π₯=0 at the centre, with faces at
π₯=Β±πΏ and half-thickness πΏ=20βmm. Define πΉπ=π·π‘/πΏ2.
For πΉπβ³0.2, use the first-term approximation
Ξ0β4πexpβ‘(βπ2πΉπ4).(3)
Use Eq. 1 with πΆπ =0. Convert the target mass concentration to molar
concentration before forming the ratio. Check that the calculated πΉπ is
consistent with the approximation.
The sheet target is one fifth of its initial hydrogen concentration after unit conversion.
The time uses the half-thickness, not the total thickness:
π‘=πΉππΏ2π·=333.415βs=5.55691βmin.
Since πΉπ>0.2, the first-term approximation is appropriate.
The next term at this time is about β2.47Γ10β8 in Ξ0, so retaining it does not change the reported time.
The source chart gives πΉπβ0.75 and π‘β333βs.
Using the full 40 mm thickness in place of πΏ would overestimate time by a factor of four.
Answer
The centre reaches 1.2Β kgβmβ3 of molecular hydrogen after 333Β s, or 5.56Β min. The resulting Fourier number supports the first-term approximation.
Consider a small spherical bead (D = 0.1 cm) that initially contains 0.2 mmol/L of the
drug griseofulvin (species A). Determine the time it will take for the concentration of
griseofulvin to reach 10% of its initial value, if its diffusivity within the bead
material is 1.5 Γ 10β7cm2/s. Upon release from the bead, the drug is
immediately consumed so that the surface concentration can be assumed to be zero.
Worked solutionFigure 3: Drug diffuses radially from an initially uniform spherical bead to a sink. The centre concentration and total amount remaining are different measures.
Data
Quantity
Symbol
Value
Diameter and radius
π,π
0.10,0.050Β cm
Initial concentration
πΆπ
0.20Β mmolβLβ1
Target centre concentration
πΆ0
0.020Β mmolβLβ1
Surface concentration
πΆπ
0
Diffusivity
π·
1.5Γ10β7Β cm2βsβ1
The source question does not specify the location of the target concentration. This worked solution uses the centre, following the supplied solution.
Assumptions
Assume spherical symmetry, a fixed bead radius π=π/2, constant temperature
and no reaction inside the bead. The centre has zero radial concentration
gradient; the surface at π=π has πΆπ =0. Define πΉπ=π·π‘/π2.
Check the approximation after computing πΉπ. Centre concentration differs
from the volume-averaged concentration. This question uses the centre,
matching the coordinatorβs worked solution.
The bead centre reaches the specified fraction when its Fourier number is about 0.304.
The target is πΆ0=0.020βmmolβLβ1, giving Ξ0=0.10.
At this time the second term of the exact centre series is β1.25Γ10β5, so the first-term ratio is accurate to about 0.013%.
The exact series gives about 5058.64βs.
The source rounds πΉπ to 0.3 and obtains 5000βs or 1.39 hours.
The diameter must first be halved.
Initial concentration cancels from the ratio because this constant-property diffusion model is linear.
Answer
The bead centre reaches ten percent of its initial concentration after approximately 5.06Γ103Β s, or 84.3Β min. This is a centre-concentration result, not a volume-average result.
In a manufacturing process, an organic solvent (methyl ethyl ketone, MEK) is used to
dissolve a thin coating of polymer film away from a nonporous flat surface (length 20
cm, width 10 cm, see Figure 1). The thickness of the polymer film is initially uniform
at 0.2 mm. In the process, MEK is added at a flow rate of 30 cm3/s to an
open flat pan of length 30 cm and width 10 cm. The depth of the liquid MEK in the pan is
maintained at 2.0 cm. Assuming that the concentration of the dissolved polymer in the
bulk solvent is essentially zero (ππ΄β=0), determine the rate of removal of
polymer coating.
Data: diffusion coefficient of dissolved polymer in MEK = 3 x 10β6cm2/s; kinematic viscosity of MEK = 6 x 10β3cm2/s; density
of polymer material (solid) = 1.05 g/cm3; density of MEK = 0.8
g/cm3; maximum solubility of polymer in MEK = 0.04 g/cm3.
Figure 1: Pan and polymer coating. The pan is 30 cm long and 10 cm wide, the coating is 20 cm long and initially 0.2 mm thick, and the liquid depth is 2 cm. Original redrawing of Figure 1 from the source question sheet.Worked solution
Data
Quantity
Symbol
Value
Coated length and width
πΏ,π
20,10Β cm
Liquid depth and initial coating thickness
π»,πΏ
2.0,0.020Β cm
Solvent volume flow
π
30Β cm3βsβ1
Diffusivity
π·π΄π΅
3.0Γ10β6Β cm2βsβ1
Kinematic viscosity
π
0.0060Β cm2βsβ1
Surface and bulk mass concentrations
πΆπ΄,π ,πΆπ΄,π
0.040,0Β gβcmβ3
Polymer density
ππ
1.05Β gβcmβ3
Assumptions
Use mass concentration πΆπ΄, so ππ(πΆπ΄,π βπΆπ΄,π) is a mass flux.
In centimetre units, πΆπ΄ in g cmβ»Β³ and ππ in cm sβ»ΒΉ give g cmβ»Β² sβ»ΒΉ.
Assume isothermal, dilute transport with constant properties, negligible
reaction in the liquid and rapid dissolution maintaining saturation at the
solid-liquid interface. The coating is thin enough to neglect changes in
flow geometry during an instantaneous rate calculation.
Model the exposed upper coating face as a laminar flat plate with a fresh
concentration boundary layer at its leading edge. Neglect pan end effects,
sidewall transfer and the velocity-profile corrections of a shallow channel.
This is the sourceβs external-flow approximation, not a resolved pan flow.
An experiment on the dissolution of fatty deposits in blood vessels involves the setup
of a tube (0.8 cm inner diameter, 50 cm in length) that is coated on the inner surface
with a thin layer (0.01 cm thick) of lipophilic material that represents the fatty
deposits, called species A (π = 1.10 g/cm3). A synthetic solution
(π = 0.02 cm2/s, π = 1.04 g/cm3) is pumped through the
tube at a rate of 20 cm3/s. The maximum solubility of A in the synthetic
solution is 20 mgΒ A/cm3 solution.
a. Develop a material balance to predict the concentration of A in the solution exiting
the tube
b. Experimental measurements indicate that the concentration of A in the synthetic
solution exiting the tube is 0.1 mgΒ A/cm3. Determine the diffusion
coefficient of A in the solution.
Worked solutionFigure 2: Liquid flows through a tube with a saturated coated wall. The wall adds dissolved material, so the outlet bulk concentration exceeds the inlet concentration.
Data
Quantity
Symbol
Value
Liquid bore diameter
π
0.80Β cm
Coated length
πΏ
50Β cm
Volume flow
π
20Β cm3βsβ1
Kinematic viscosity
π
0.020Β cm2βsβ1
Saturated wall concentration
πΆπ΄,π
20Β mgβcmβ3
Inlet and measured outlet concentrations
πΆπ΄,π,πΆπ΄,π
0,0.10Β mgβcmβ3
The supplied solution uses Shπ=πππ/π·π΄π΅=3.66 for fully
developed laminar mass transfer. The calculation below checks that assumption.
It treats the quoted diameter as the initial liquid bore and takes the inlet
concentration to be zero.
Assumptions
Assume steady liquid transport during the measurement, negligible axial
diffusion, no bulk reaction, constant flow rate and a saturated coated wall.
The cross-section balance uses a mixed-mean concentration, not a uniform
local radial profile. Assume an upstream uncoated length establishes the
laminar velocity profile before the coated section. The source correlation
3.66 additionally requires a developed concentration profile at fixed wall
concentration. Reynolds number alone does not establish this condition.
An upstream uncoated tube can develop the velocity profile, but cannot develop a solute profile before dissolution starts.
The supplied 3.66 model therefore does not provide a justified physical diffusivity for this entrance-region experiment.
The data establish the effective ππ and mass removal rate under the stated balance.
Inferring a reliable diffusivity additionally needs an appropriate developing-flow correlation or a resolved transport model, with the velocity entrance condition specified.
The conditional numerical result must not be reported as a validated measured diffusivity.
Coating thickness and solid density describe the available inventory; liquid density is unnecessary when kinematic viscosity is supplied.
Answer
The measured outlet concentration gives an effective ππ=7.98Γ10β6Β mβsβ1 and a mass pickup of 2.0Β mgβsβ1. Assuming Sh=3.66 gives π·π΄π΅=1.74Γ10β8Β m2βsβ1, but the concentration-entry check fails. This diffusivity is conditional and is not validated by the experiment.