THERMO Spoken Here! ~ J. Pohl © ~ 2017 (D5600-207)

Atmospheric Engine

Thomas Newcomen in­vented this piston and cylinder engine in 1712. The power stroke of the atmopheric steam engine is an ingenous thermodynamic event. A fine antimation of the engine's operation has been coded. Animation URL!.

"How it works," is explained using three idealized instances: (i) Initial Condition, (ii) Water Injection, and lastly, (iii) Equilibration of the water and steam.

(i) Initial Condition:  Briefly described, from a fired boiler, steam is piped into the cylinder bottom and out its top until all air is ejected and the cylinder wall and piston are "steam" hot. Next, as the piston rises from least to max­imum height, steam at 100°C (g1) fills the maximum cylinder volume. Fitted to the large piston/cylinder is a smaller "injection" cylinder full of water (f1). Our system is the steam and water. Fig 1. depicts the situation.

It is convenient to calculate properties of the Initial Condition now. Table 1. contains necessary assumptions and properties of water.

Table 1:  Assumptions - Properties @ 1 atmosphere
Steam (g1)    Vcyl = 1.000m³    T = 100°C
  vg = 1.667m³/kg,  ug = 2506kJ/kg
Water (f1)  Vinj = 0.001m³    T = 24°C
  vf = 0.001m³/kg  uf = 100.8kJ/kg

Necomen's Engine requires the three equations, two physical principles: volume is constant and energy constant (Eqn-6), along with tabular properties for water be solved simultaneously. The system volume is that of the steam and injection water.


The discrete masses of the system are determined via specific volumes and volumes.


The cylinder of steam saturated at one atmosphere and the injector water saturated at 24°C have internal energies. The system internal energy is:



(ii) Injection: Into the steam in the sealed piston/cylinder is violently jetted the injection water. We assume this to occur instantly (hence no heat). The system mass is constant and while the piston can move, there is too little time for volume to change (hence no work). Also the internal energy remains what it was. For the injection event, mass, volume, and internal energy are constant.


(iii) Equilibration: When water is sprayed into the hot steam, it will condense and the resulting pressure will be less than the original, 101.3kPa. Our interest is to prove the "result of quench" will be lower pressure. The thermodynamic event is called "equilibration," whereupon the initially separate steam and water transform from initial temperatures and one atmosphere to a same temperature and lower pressure. Equations must be written then used with water properties to prove this fact.

For equilibration, system volume is constant:


For equilibration, system internal energy is constant:



Below are "final forms" of Equation (4) and (5) algebraically arranged such that each has its second quality left of equality. These equations (6) apply to the "equilibated" or quenched state." Below Eqn-6, Table 2 is a brief listing of properties of water.


Table 2:  Saturation Properties for Water
p(kPa) vf(m³/kg) vfg(m³/kg) uf(kJ/kg) ufg(kJ/kg)
90 0.001 1.87 404.9 2097
60 0.001 2.73 359.7 2489

Necomen's Engine requires the three equations, two physical principles: volume is constant and energy constant (Eqn-6), along with tabular properties for water be solved simultaneously.

The two equations and table constitute a set that can be solved for the second temperature and pressure. But since the table is involved, we must use trial and error. So we...
   (i) Guess the equilibrated pressure, p*.
   (ii) Use Table 3. to obtain water properties:  vf(p*), vfg(p*), uf(p*), and ufg(p*).
   (iii) When the properties of water at p* are substituted into Eqns-6 two values of xfg(p*)    are obtained. If the two values xfg(p*) are not equal - try again.

Table 3 presents trial pressures and the qualities (xfg,2) associated for tries at 60kPa and 90kPa. We observe the "constant volume" curve and the "constant internal energy" to curve cross each other - indicating a solution of about 75kPa.

Table 3:  Trial and Error Data
p(kPa) [xfg,2] V = const [xfg,2] U = const
90 0.33 0.28
60 0.23 0.26

75kPa seems a pretty high pressure. We would like it to be lower. But this is our first, guessed case of the apparatus. This is not Necomen's solution. Further design might yield better performance. However, I'm tired of this now... Maybe I'll do more later!

Tags: None