Chemical Engineering Tutorials

Wednesday, 17 September 2025

VAPOR COMPRESSION REFRIGERATION CYCLES

Carnot Refrigeration Cycle (Carnot Heat Pump)

A refrigeration cycle is just a reversed heat engine cycle.

Therefore, heat is transferred from a low temperature level to a high temperature level. However, according to the second law of thermodynamics this cannot be accomplished without the use of external energy. The working mediums that are used in compression refrigeration systems are called refrigerants.

The following figure shows a typical Carnot refrigeration cycle:


  1. The refrigerant evaporates at a constant temperature and pressure as heat is absorbed from the low temperature region in the evaporator.
  2. The vapor leaving the evaporator is then compressed to a higher pressure and external work is required in this process.
  3. Heat is rejected to the high temperature region in the condenser as the refrigerant condenses at constant temperature and pressure.
  4. To complete the cycle, the liquid from the condenser is returned to its original state by an expansion process.

The overall process is represented on a T-S diagram in the following figure:


The coefficient of performance, COP is used to measure the performance of a refrigeration cycle. It is the ratio of the refrigeration obtained to the work required, i.e.,


The COP of a Carnot refrigeration cycle is given by:


A ton is a common unit used in practice to describe the refrigeration effect. One ton of refrigeration is the term used to refer to 12,000 Btu/h.

Therefore, a chiller or condensing unit with a cooling capacity of 60,000 Btu/h is said to have a capacity of 5 tons. One ton of refrigeration approximates to 3.5 kW of cooling.

As seen from the T-S graph of a Carnot refrigeration cycle, both the compression and expansion steps in the Carnot refrigeration occur within the two-phase region. The compression of a two-phase mixture, however, is very difficult in practice and is generally avoided. This problem can be eliminated by simply allowing the refrigerant to evaporate completely in the evaporator resulting in the production of a saturated vapor.

On the other hand, the fluid passing through the expander is mostly liquid and its specific volume is relatively low. Thus, the amount of work produced by the expander is not appreciable. For this reason, much less expensive and almost maintenance-free throttling expansion devices are preferred over expanders in practice.

The modified cycle and its representation on a T-S diagram are shown in the following two images:



The vapor compression refrigeration cycles can also be represented on a P-H diagram as shown:


Refrigerants

The design of a vapor compression refrigeration system is greatly influenced by the physical, thermodynamic and chemical properties of the refrigerant used.

The desirable properties for a refrigerant can be summarized as follows:

  • Positive evaporating pressures: This prevents leakage of atmospheric air into the system during operation.
  • Moderately low condensing pressures: This allows the use of light weight equipment on the high-pressure side of the system.
  • Low freezing point.
  • High latent heat of vaporization and relatively high critical temperatures: A high latent heat means a high refrigeration effect per kg of refrigerant circulated and low power cost for circulation.
  • Low cost.
  • Inertness and stability
  • Not be toxic, irritating or flammable.

 

However, no single compound meets all these requirements.

Ammonia and sulfur dioxide were the early refrigerants for commercial use due to their high latent heats of vaporization. However, they have the obvious drawbacks of being highly toxic and corrosive, and with NH3 being flammable.

Research efforts in the 1920’s led to the conclusion that small molecules having C-F bonds (fluorocarbons) were suitable choices, with dichlorodifluoromethane (CCl2F2) having the best properties. However, these compounds accumulated in the atmosphere and resulted in destruction of O3 molecules in the ozone layer.

Hydrofluorocarbons (HFCs) were developed as alternatives to the ozone-depleting refrigerants. HFC-134a (1,1,1,2-tetrafluoroethane - CF3CH2F) has become the refrigerant of choice to replace CFC-12 in most refrigeration and auto air conditioning systems.


Monday, 1 September 2025

The Gaussian Integral

The Gaussian or Probability Integral is an essential concept in mathematics particularly in the fields of probability theory, statistics and quantum mechanics.

The Gaussian integral, closely related to the erf function, is the integral of the one-dimensional Gaussian function over (-∞, ∞).

The Gaussian Integral can also be defined as the integral of the exponential of -x2 over the entire real line.

It can be calculated using the trick of combining two one-dimensional Gaussians:


In this case, Since the variable in the integral is a dummy variable i.e, it integrates out in the end, we can rename from x to y.

When switching to polar coordinates we get:

Example

Prove the following:

Solution

Convert the integral into the polar coordinates (r, θ) where x2 + y2 = r2

and dxdy = rdrdθ:


Evaluating the integrals:


Therefore:





Thursday, 21 August 2025

CONTROL VOLUME APPROACH

After studying the various type of flows, we need to address how to determine the velocity in the flow field. Two approaches as used for this:

The Lagrangian Approach

  • This is used in solid mechanics and involves describing particle’s motion by position as a function of time.
  • Can be used to describe the motion of an object falling due to gravity: s = ½gt2. At any time, the distance from the body’s original position is known.
  • This approach is difficult to use in fluid mechanics because a fluid is a continuous medium i.e., a single fluid volume changes shape, and different fluid particles within the fluid volume are traveling at different velocities. Thus, due to the nature of fluids, the Lagrangian Approach is generally not a desirable analysis method.

The Eulerian (or control volume) Approach

  • This is preferred in fluid mechanics.
  • In this method, a region in the flow field is chosen for study. For example, consider flow draining from a sink, as shown below. A control volume is chosen around the region of study and is bounded by the dashed line called the control surface while everything outside is called the surroundings. The control volume or shape is chosen for convenience in solving the problem.

  • Generally, the control volume shape is selected so that fluid and flow properties can be evaluated at locations where mass crosses the control surface or, if no mass enters or exits, where energy crosses the control surface. Furthermore, the control surface can move or change shape with time as illustrated above.
  • The control volume is to fluids as the free-body diagram is to solids. 


The Eulerian approach is suitable in solving fluid mechanical problems. The aim is to develop equations of fluid dynamics that are conservation equations each developed from a general conservation equation:

  • the continuity equation (conservation of mass),
  • the momentum equation (conservation of linear momentum),
  • the energy equation (conservation of energy).


Let N be defined as a flow quantity (mass, momentum, or energy) associated with a fluid volume or system of particles and n represents the flow quantity per unit mass. Thus:

Using the following image consider a system of particles at two different times: V1 at t1 and V2 at t2 where V1 is bounded by the solid line and V2 is bounded by the double line. V1 consists of VA and VB while V2 consists of VB and VC. The control volume is bounded by the dashed line.



The amount of flow quantity N contained in V1 is the amount in VA and the amount in VB at t1 which is NA1 + NB1. The amount of flow quantity contained in V2 is the amount in VB and VC at t2, which is NB2 + NC2. During the time interval, the change in N is, therefore:



To obtain a specific limiting expression for this last term, let us consider a differential area dA through which fluid particles flow:


The fluid velocity at the differential area, dA is Ṽ, which has components normal and tangential to dA: Vn and Vt. The tangential velocity carries no fluid out of the control volume with it as all fluid leaving dA is in the Vn direction. During the time interval Δt, the mass of fluid crossing dA is:


By substituting this into Equation 3, we get the general conservation equation:


Equation 4 gives us a relationship between the various quantities associated with a system of particles i.e., equation 4 simply means:




Wednesday, 16 July 2025

KINEMATICS OF FLOW

Types of flow and how to characterize them:

  • Closed-conduit flows: Type of flow that is completely enclosed by restraining solid surfaces. Examples include flow through a pipe.
  • Open-channel flows: In this flow type there is one surface exposed to atmospheric pressure. Examples include flow in a river and flow in a spillway.
  • Unbounded flows: In this flow the fluid is not in contact with any solid surface. Examples are the jet that issues from a tap and the jet from a can of spray paint.

KINEMATICS OF FLOW

In fluid flow situations velocity often needs to be determined. However, velocity generally varies in the flow field. The flow can thus be classified according to how the velocity varies.

One-dimensional flow occurs if the parameters of both the fluid and the flow are constant at any cross section normal to the flow. It can also occur when represented by average values over the cross section. Although the flow velocity can change from point to point it remains constant at each location.

The figure below illustrates velocity distributions for a one-dimensional flow:

Two-dimensional flow occurs when the fluid or flow parameters have a gradient in two directions. For flow in a pipe, the velocity at any cross section is parabolic thus the velocity is a function of the radial coordinate. Additionally, a pressure gradient exists in the axial direction that maintains the flow i.e., a pressure difference between inlet and outlet that causes the fluid to flow. Even though the fluid flows in one direction, due to the pressure and velocity gradient, the flow is classified as two-dimensional. This is illustrated as follows:

Another example of two-dimensional flow occurs when at the constant area section velocity is a function of one variable but at the convergent section, velocity is a function of two space variables. Furthermore, a pressure gradient exists that maintain the flow. This is illustrated as follows:

Another example of a one-directional, two-dimensional flow where gradients exist in two dimensions is illustrated below:


Three-dimensional flow occurs when the fluid velocity or flow parameters vary with respect to all three space variables. Thus, a gradient exists in three directions.

Steady Flow occurs when conditions do not vary with time or when fluctuating variations are small with respect to mean flow values, and the mean flow values do not vary with time.

Unsteady Flow occurs when flow conditions change with time.

Quasi-steady Flow occurs in some unsteady flows, where it is permissible or even necessary to assume that steady flow exists to obtain a solution.



Wednesday, 18 June 2025

Ideal Reactors - Solved Examples

 Example 1

The following reaction is irreversible and first order:

The reaction is carried out in a PFR with 80 tubes. Each tube has a diameter of 5cm and a height of 1m. The feed consists of reactant A and 30% inerts flowing at 200kg/h with a pressure of 50bar and temperature of 600°C (873K). The output conversion is 90%. What is the average residence time? (The molecular weight of A is 60g/mol & R = 0.0821 L • atm/ mol • K)

Solution

The space time for a first order reaction in a PFR is given by the following equation:

The mass flow in each tube is: (200/80) = 1.25kg/h

In molar basis:

Thus, the initial concentration of A is given by:

We can now calculate the initial volumetric flow:

The volume of each tube can be calculated as follows:

The residence time can be calculated as follows:

The expansion parameter, εA can be determined as follows:

The rate constant, k, can now be solved using equation 1:

Thus, the average residence time can now be solved using:


Example 2

The following second order irreversible reaction is carried out in gas phase in a PFR reactor:


The PFR feed consists of 50% by weight reactant A (MW = 40) and the rest is inert (MWinert = 20). The reaction occurs at a constant temperature of 70C and pressure of 5.25atm. The rate constant is 400m3/(kmol.ks).

For a production rate of R set at 30kmol/h with a 40% conversion, what should be the volume of the reactor?

Solution

For a second-order reaction, the rate can be expressed as:

Substituting the PFR design equation into the above equation and integrating we obtain:


With a mass basis of 1.0g, we can calculate εA. The molar values for each feed components are 0.5/40 = 0.0125 mol of A and 0.5/20 = 0.025 mol of inert. Thus, the volume balance is as shown below:

The expansion parameter, εis = (0.05 - 0.0373)/0.0373 = 0.340. The initial concentration of A is:

Equation 1 can now be solved: 

Given that the output molar flow FR = 30kmol/h, the volumetric flow rate, vo is:

Thus, the reactor volume is:










Friday, 2 May 2025

Solved Example for Thermodynamics of Energy Conversion

For the following examples, Appendix B for steam values that I have referred to in these questions was obtained from:

M.D Koretsky, Engineering and Chemical Thermodynamics, Wiley, 2004.

Example 

Consider a power plant operating on a Rankine cycle using steam as the working fluid. The boiler pressure is 2.5MPa and the steam leaving the boiler is superheated to a temperature 126°C above its saturation temperature. The condenser temperature is 50°C and it discharges saturated liquid. The efficiency of the turbine is 0.90 and of the pump 0.80 as compared to reversible and adiabatic machines operating at the same pressure ranges.

a) Sketch the cycle on a T-S diagram.

b) Calculate the thermal efficiency of the cycle.

c) Compare the thermal efficiency of this Rankine cycle with the thermal efficiency of a Carnot cycle receiving heat at the saturation temperature of steam at 2.5MPa and discarding it at 50°C.

Solution

a) 

b) In order to calculate the thermal efficiency, we need to calculate enthalpy values at states 1, 2, 3 and 4.





c)








VAPOR COMPRESSION REFRIGERATION CYCLES

Carnot Refrigeration Cycle (Carnot Heat Pump) A refrigeration cycle is just a reversed heat engine cycle. Therefore, heat is transferre...