Show Summary Details
More options …

# Nonlinear Engineering

### Modeling and Application

Editor-in-Chief: Nakahie Jazar, Gholamreza

Managing Editor: Skoryna, Juliusz

CiteScore 2018: 1.33

SCImago Journal Rank (SJR) 2018: 0.313
Source Normalized Impact per Paper (SNIP) 2018: 0.439

Open Access
Online
ISSN
2192-8029
See all formats and pricing
More options …
Volume 8, Issue 1

# Mathematical Modelling Comparison of a Reciprocating, a Szorenyi Rotary, and a Wankel Rotary Engine

Luis F. Espinosa
/ Petros Lappas
Published Online: 2018-10-15 | DOI: https://doi.org/10.1515/nleng-2017-0082

## Abstract

This paper provides an explanation of the geometry, design, and operational principles for the three engines; having special emphasis in the Szorenyi rotary engine which has a deforming rhombus revolving inside a mathematically defined stator. A basic ideal mathematical simulation of those engines were performed, assuming the Otto cycle for the three engines. Also, it assumes the volumetric efficiency of 100%, a wide-open throttle (WOT), no knock nor any mechanical or thermal losses. This simulation focuses on how the fuel burns during combustion, creating pressure and thus, net work. A comparison in pressure traces and cycle performance is made. The study concludes analysing and comparing the ignition advance; finding the best advance for each engine thus the net work between the three engines during one working cycle. Finally, this paper analyses how the different volume change ratio for the combustion chamber of the Szorenyi, Wankel and the reciprocating engine have an effect in the pressure, net work and thermal efficiency generated inside the chamber during combustion for every working cycle.

## 1 Introduction

Rotary engines have seen many different designs over the years. The most successful of those, the Wankel engine, achieved production, but it had different problems such as high hydrocarbon emissions, increased oil and fuel consumption, more expensive manufacturing costs, and limited fuel flexibility [1]. For all this reason, they have not been developed as much as reciprocating engines. In 2004, Peter Szorenyi got a patent approved with his partner Peter King of a hinged rotor internal combustion engine. The Szorenyi engine may offer some advantages over the Wankel rotary engine and over the reciprocating engine [2].

This paper explores the potential of the Szorenyi engine compared to the exiting Wankel rotary engine and the reciprocating engine. In consequence, the advantages, disadvantages of the Szorenyi engine against the reciprocating engine and Wankel rotary engine will be determined.

By running a series of basic fuel combustion simulations at different advance angles, and same angular speed. Their pressure, volumes, net-work and thermal efficiency will be compared under the same working conditions. Also, it is assumed optimal conditions: not taking into account autoignition, thermal, mechanical and friction losses. In consequence, the differences on how the volume change, will affect the pressure change ratio and net work achieved from them during the cycle.

## 2 Background

Reciprocating engines consist of a piston inside a cylinder to create variable pressure and volume [3]. But there are not the only internal combustion engines. Rotary engines are a middle step between reciprocating engines and turbines. Rotary engines form ‘a plurality of variable volume chambers’ between a rotor and the outer part known as the stator [4]. Unfortunately, Rotary engines have several disadvantages as higher fuel consumption and emissions. Also, they are a less fuel flexible engine, with sealing, lifespan, compression ratio disadvantages and longer flame travel path and other combustion effects that makes the combustion unstable [58].

## 2.1.1 The Wankel Engine

The Wankel engine is a rotary engine where its most successful variety has a triangular shaped rotor. This rotor ‘performs a planetary rotary movement relative to the outer body’ [4]. This movement is generated using a planetary gearing, at the side plate to maintain the correct phase between the rotor and the eccentric shaft rotations but, it generates unbalanced forces. This forces can be cancelled by just adding balancing weights to balance the rotor [7, 8]. It also forces the rotor to revolve at two thirds of the crank’s angular velocity, meaning that , it produces three power stroke every two crankshaft revolution [5, 79]. Finally, breathing is through ports which can be in the stator face or at the sides of the housing [5, 9, 10].

## 2.1.1.1 Wankel Engine Geometry

The Wankel engine has an epitrochoid stator in which the rotor rotates. (1) and (2) define the coordinates in X and Y axis of the stator shown in Fig. 1. Also, they were used to determine the volume of one chamber at any given time [7, 11, 12].

Fig. 1

Wankel Engine Geometry

$X=e⋅cos3β+Rw⋅cosβ$(1)

$Y=e⋅sin3β+Rw⋅sinβ$(2)

## 2.1.2 The Szorenyi Engine

The Szorenyi engine has a four-segment hinged rotor assembly which deforms and adapts continuously to the stator profile during its rotation, changing from a square to a rhombus and back. Therefore, the Szorenyi is a rotary engine with four combustion chambers [2, 13]. Since, there is no gearing, it has four power strokes per crankshaft revolution. Thus, one rotor of the Szorenyi engine generates the same number of power strokes as a reciprocating 4 stroke engine with eight cylinders [2, 13]. Therefore, the Szorenyi engine should have higher power output at lower RPM having a lower mechanical complexity [13].

## 2.1.2.1 Szorenyi Geometry

The stator curve, patented as the Szorenyi curve, is determined by the right isosceles triangle, as shown in Fig. 2 [6]. It is ‘generated by the base extremes A and B of an isosceles right angle, translating and simultaneously rotating triangle’ [2]. In the present paper in Fig. 1 is the angle used to generate the x and y co-ordinates of the stator profile. The crankshaft angle of the Szorenyi engine, as used in this paper, is defined as the angle between the positive vertical axis and the centre of the rotor segment. Using the notation in Fig. 1, (3) and (4) define all points on the stator curve in Cartesian coordinates [13]. Then is possible to find the volume at any given time applying the same basic concept than in the Wankel engine volume [7].

Fig. 2

Szorenyi engine curve construction diagram

$XA=12{c2−sin22θ−sin2θ}{sinθ−cosθ}$(3)

$YA=12{c2−sin22θ−sin2θ}{1−2sin2θsinθ+cosθ}$(4)

## 3.1 Thermal Simulation

For the purposes of the present simulation, the total volume chamber for the three engines was defined on 125cc. The volumes on several stages were calculated and verified using CAD software. The simulation was executed with a spreadsheet to calculate the chamber volume, fuel fraction burned, pressure and the end-gas temperature at any time. The calculations used general and specific parameters for each engine in order to have a net work output per chamber cycle of the three engines under equivalent conditions. These parameters are described in Table 1, for the modelling of the engines [9, 10, 1420].

Table 1

Due to the limited scope of this paper, it should be noted that the present simulation does not consider pumping, heat, mechanical losses, knock nor their different efficiencies.

Note: The same compression ratio can be achieved using a different ω value by removing or adding material to the rotor face. The aim of the present work is not to determine the optimum value ‘ω’. Although it will affect the volume history of the combustion chamber. The scope of this paper does not extend to optimising the ‘ω’ value. By varying the ‘ω‘ value, the volume history of the chamber could be optimised and this might lead to different results for the Szorenyi engine than those presented in the paper.

## 3.2 Engines Combustion Chamber

For the Szorenyi and Wankel rotary engines, only one chamber was modelled for each engine. The engine modelling determined the volume of the combustion chamber at any time during the power stroke cycle. Then the points of the rotor segment in contact with the stator at any time were determined. After that, the volume between the rotor and stator at any time was calculated. The resulting area was multiplied by the rotor depth and a scale factor to determine the working volume [7].

For the Szorenyi rotary engine, in order to keep the simulation as simple as possible, the value ’ω’ shown in (3) and (4) was modified until the desired compression ratio was achieved. In order to validate the geometry, the coordinates generated in the mathematical simulation were exported to a CAD software as seen in Fig. 3. Also, Fig. 3 illustrates the concept of calculating the projected area of the working chamber of the Szorenyi rotary engine.

Fig. 3

Szorenyi engine chamber cross-sectional area at 0, 30, 60 and 90 degrees of the crankshaft angle.

For the Wankel rotary engine, in order to have the desired compression ratio, the rotor surface is further modified by a formation of a calculated rotor pocket, which reduces the volume in the chambers. Thus, increasing the compression ratio until the desired value [17].

As in the Szorenyi engine, careful consideration was taken in the Wankel rotary engine order to determine the area between the rotor and stator at any time as shown in Fig. 4 [7].

Fig. 4

Wankel engine chamber cross-sectional area at different rotor positions corresponding at every 45 degrees’ rotation of the crankshaft

Likewise, the cylinder volume at any crank angle position was calculated for the reciprocating engine [10].

## 3.2.1 Fuel Fraction Burned – Wiebe Function

To create an accurate model of the engine which determines the temperature and pressure through the power stroke, it is important to know how much fuel is being burn inside the combustion chamber. The Wiebe function is used to determine this basic parameter as shown in (5) [9, 10, 2023].

$Xb=1−exp[−C(θ−θ0Δθ)m+1];C=ln(1−xd)$(5)

All three engines rotate at the same speed but, their design forces them to have combustion cycles with different duration. In order to compare all three engines combustion in an equivalent way, the fuel mass fraction burned ‘mfb’, which is the fuel burned since the spark was generated, was compared to the percentage of the combustion period as it is shown in Fig. 5. It can be seen that the reciprocating and Wankel engines have a very similar mfb when it is compared to the percentage of combustion. And the Szorenyi engine has a slightly lower curve, which means that the fuel is burned slower at the beginning and faster at the end. This will have an effect on how the pressure is generated in the combustion chamber.

Fig. 5

Fuel Fraction burned During Combustion

## 3.3 Pressure Change Determination

The incremental change in combustion chamber pressure is obtained based on the ideal gas law, the first law of thermodynamics, Mayer law, ratio of specific heats, and heat release equation [10, 15, 24]. The resultant pressure increment can be expressed in the form of (6) [15, 23].

$dP=QLHV*mf*dXb*(k−1)−k*Po*dVV$(6)

## 3.4 Work Determination and Thermal Efficiency

The resulting data from the simulation, based on the P-V diagram, and assuming the friction value to be zero, makes it possible to determine the net work done per cycle [10].

Heat cannot be fully converted into work [24]. Therefore, thermal efficiency was calculated, in order to indicate how effective, the engine is in converting heat input into a mechanical work [10].

## 4 Results and Discussion

The present research and modelling of combustion processes is based on the parameters presented previously at Table 1. Several simulations, using a Wiebe based thermodynamic model, were run to determine the best advance of ignition. The best advance of ignition is the one which provides the maximum net work for each engine.

The model only accounts for the positive net work generated by the engine. It does not take into account the negative work of intake and exhaust. It assumes a wide-open throttle.

## 4.1 Volume Displaced

In relation to the crankshaft position, the volume in the combustion chamber of the Szorenyi, the Wankel rotary engine and the reciprocating engine are different throughout the working cycle. The volume history of all of them was compared during the compression phase as shown in Fig. 6.

Fig. 6

Volume vs % of compression phase

Fig. 6 indicates that the Szorenyi engine has the slowest and better rate of volume change near TDC and the fastest near BDC. A slower rate of change near TDC is closer to the ideal Otto cycle because it brings a steadier volume for combustion. Therefore, a higher fuel mass fraction is burned.

## 4.2 Net Work

The net peak work output with respect to the ignition timing is shown in Fig. 7.

Fig. 7

Variation of net work produced by the Szorenyi engine, Wankel engine and the reciprocating engine at different advance of ignition

Consequently, it is observed in Fig. 7 the net work of the Szorenyi rotary engine is the greatest per working cycle, followed by the Wankel rotary engine and finally the reciprocating engine as summarised in Table 2.

Table 2

Net work comparison between engines

P-V diagrams of the three engines at their optimum ignition advance are shown in Fig. 8. This results in a highest maximum pressure for the Szorenyi engine, followed by the Wankel engine. Both maximum pressure and work are higher than the reciprocating engine.

The enclosed areas of the P-V diagram were computed with the result that the Wankel engine ideally produces 0.10% more net work per cycle than the reciprocating engine. Also, the Szorenyi engine ideally produces 0.462% more work per cycle than the equivalent reciprocating engine.

Fig. 8

P-V diagram comparison of Szorenyi, Wankel and Reciprocating engine

Furthermore, in Fig. 9, the top pressure of the Wankel engine is 4.58% higher than the reciprocating engine and the Szorenyi engine is 5.515% higher than the same reciprocating engine. This leads to a higher ideal net work per cycle for the Szorenyi engine followed by the Wankel engine and the least net work per cycle is produced by the reciprocating engine.

Fig. 9

Pressure comparison in a complete cycle of the Szorenyi, Wankel and reciprocating engines

## 4.3 Thermal Efficiency

Thermal efficiency was calculated for the Szorenyi, Wankel and the reciprocating engine. The results are shown in Table 3.

Table 3

Fuel conversion efficiency comparison between Szorenyi, Wankel and Reciprocating engine

It can be seen in Table 3 that under the same working conditions, the Wankel engine has a thermal efficiency 0.079% higher than the reciprocating engine. Also, the Szorenyi engine has a thermal efficiency 0.462% higher than the same reciprocating engine.

## 4.4 Power

Power is a function of time and torque, therefore, for one chamber of each engine, and the three engines rotating at the same speed at the crankshaft, 3000 RPM, will have different chamber variation. Therefore, the ideal power output is stated in Table 4 for each engine at the same crankshaft speed.

Table 4

Power output comparison of one chamber of the Szorenyi, Wankel and a Piston engine

## 5 Conclusions

As a result of mathematical modelling, it is established that:

• The Szorenyi engine is more sensitive to changes in the advance of ignition respect to the crankshaft rotation angle than the reciprocating engine and it has a greater effect on the net work produced.

• The Szorenyi engine has a shorter rotation of the crankshaft during combustion than the reciprocating and Wankel engines. Therefore, it must have a smaller advance of ignition than the reciprocating and Wankel engine with respect to the crankshaft.

• The net work per cycle produced by the Wankel engine is 0.079% higher than the reciprocating engine. Also, the Szorenyi rotary engine is 0.462% higher than the reciprocating engine for the same displaced volume.

• The Szorenyi engine has the highest net work per cycle, and therefore the highest thermal efficiency, as shown in Table 2 and Table 3. The Szorenyi engine’s net work and thermal efficiency is followed by the Wankel engine and, finally, the reciprocating engine has the lowest net work per cycle.

• The Szorenyi engine higher power output is 100.90% higher than the reciprocating engine. Also, the Wankel engine is 33.43% higher than the reciprocating engine. These differences are created by the net work, thermal efficiency and number of power strokes per crankshaft revolution of each engine when all engines are working at the same crankshaft speed.

## References

• [1]

S. Pehan and B. Kegl. (2001). Rotary Engine Design [SAE Technical Paper 2001-01-3194]. Available: http://dx.doi.org/10. 4271/2001-01-3194 Google Scholar

• [2]

P. Szorenyi, “Hinged rotor internal combustion engine,” U.S. Patent US 6.718.938 B2, Apr. 13, 2004. Google Scholar

• [3]

Z. Spakovszky, M. Greitzer, and I. Waitz. (2008, October 10). Thermodynamics & Propulsion. Available: http://web.mit.edu/16.unified/www/SPRING/propulsion/notes/node25.html Google Scholar

• [4]

F. Wankel and E. Hoeppner, “Rotary internal combustion engine,” US Patent 2988065, June 13, 1961. Google Scholar

• [5]

L. Peng, “Fundamental research on a unique rotary machine,” PhD, School of Graduate Studies, McMaster University, Canada, 1994. Google Scholar

• [6]

R. W. Richardson, “Automotive Engines for the 1980’s,” in Eaton’s Worldwide Analysis of Future Automotive Power Plants, Eaton Corporation, ed, 1973, 40-41. Google Scholar

• [7]

K. Yamamoto, Rotary Engine. Japan: Toyo Koyo co. Ltd., 1971. Google Scholar

• [8]

K. Yamamoto, Rotary engine. Hiroshima, Japan: Sankaido, 1981. Google Scholar

• [9]

T. J. Norman, “A performance model of a spark ignition Wankel engine : including the effects of crevice volumes, gas leakage, and heat transfer,” Master´s Degree, Dept. of Mechanical Engineering, Massachusetts Institute of Technology, 1983. Google Scholar

• [10]

J. B. Heywood, Internal combustion engine fundamentals vol. 930: Mcgraw-hill New York, 1988.

• [11]

R. Sierens, R. Baert, D. E. Winterbone, and P. C. Baruah, “A Comprehensive Study of Wankel Engine Performance,” 1983.

• [12]

M. Alexandru-Catalin and B. N. S. Liviu-Constantin, “ROTARY INTERNAL COMBUSTION ENGINES,” Universitatii Maritime Constanta. Analele, 2014, 15, 93-99

• [13]

P. King. (2013, August 2013) The Szorenyi Rotary-the four chamber rotary engine. Autoengineer [Technical Talk]. 4.

• [14]

J. Heywood and P. Watts, “Parametric studies of Fuel Consumption and NO Emissions of Dilute Spark-Ignition Engine Operation Using a Cycle Simulation,” in Institute of Mechanical Engineers Conference Publication C, 1979, 117-127.

• [15]

G. H. Abd Alla, “Computer simulation of a four stroke spark ignition engine,” Energy Conversion and Management, 2002, 43, 1044-1047.

• [16]

A. Boretti, S. Jiang, and J. Scalzo, “A Novel Wankel Engine Featuring Jet Ignition and Port or Direct Injection for Faster and More Complete Combustion Especially Designed for Gaseous Fuels,” SAE Technical Paper 0148-7191, 2015.

• [17]

M. S. Raju and E. A. Willis, “Three-Dimensional Analysis and Modeling of a Wankel Engine,” SAE Technical Paper 0148-7191, 1991.

• [18]

A. Boretti, “CAD/CFD/CAE Modelling of Wankel Engines for UAV,” SAE Technical Paper 0148-7191, 2015.

• [19]

S. D. Hires, R. J. Tabaczynski, and J. M. Novak, “The Prediction of Ignition Delay and Combustion Intervals for a Homogeneous Charge, Spark Ignition Engine,” SAE transactions, 1978, 1053-1067.

• [20]

J. I. Ghojel, “Review of the development and applications of the wiebe function: A tribute to the contribution of ivan wiebe to engine research,” International Journal of Engine Research, 2010, 11, 297-312.

• [21]

G. Abbaszadehmosayebi and L. Ganippa, “Determination of specific heat ratio and error analysis for engine heat release calculations,” Applied Energy, 6/1/ 2014, 122, 144.

• [22]

G. Abbaszadehmosayebi and L. Ganippa, “Characterising Wiebe Equation for Heat Release Analysis based on Combustion Burn Factor (Ci),” Fuel, 3/1/ 2014, 119, 301

• [23]

S. Klein and G. Nellis. (2011). Thermodynamics. Available: http://RMIT.eblib.com.au/patron/FullRecord.aspx?p=807134 Google Scholar

• [24]

S. Klein and G. Nellis, Thermodynamics: Cambridge University Press, 2012.

Revised: 2018-03-28

Accepted: 2018-07-06

Published Online: 2018-10-15

Published in Print: 2019-01-28

Citation Information: Nonlinear Engineering, Volume 8, Issue 1, Pages 389–396, ISSN (Online) 2192-8029, ISSN (Print) 2192-8010,

Export Citation