Graphical Abstract Figure
Graphical Abstract Figure
Close modal

Abstract

The spoiler is vital for optimizing rotary engines, yet its impact on the flow and combustion process within the cylinder remains ambiguous, thereby impeding the optimization efforts for rotary engines. This article, by analyzing the motion of the rotary, discovers a form of high-speed local gas flow within the cylinder due to local pressure differences resulting from interactions among the combustion chamber, cylinder body, and spoiler structure. This phenomenon is named “pressure differential flow” to differentiate it from the forced flow induced by the spoiler. Using mathematical models and three-dimensional simulations, we analyze the intensity of pressure differential flow at various spoiler heights and its regulatory effects on the flow and combustion characteristics within the cylinder. The results indicate that the flow caused by the spoiler in the cylinder is primarily divided into forced flow and pressure differential flow, with the intensity of the latter increasing as the spoiler height increases. When the spoiler height is greater than 75% of the maximum height, the pressure difference flow becomes more apparent, with both forced flow and pressure difference flow coexisting in the cylinder; when the spoiler height is less than 75% of the maximum height, the pressure difference flow is less noticeable, and the forced flow caused by the spoiler dominates. Pressure differential flow can reduce ignition delay and increase the maximum cylinder pressure, but it can also delay ignition timing and reduce combustion stability.

1 Introduction

Internal combustion engines serve as a reliable power source across various industries. Today, cost-effective and dependable small engines are prevalent in motorcycles, all-terrain vehicles, boats, and small aircraft [1,2]. Despite their widespread use, small engines exhibit lower efficiency than larger piston engines. For instance, a specific small-displacement four-stroke reciprocating engine achieves a peak efficiency of only 12–15% [35]. In contrast, in a traditional four-stroke engine, where the piston momentarily halts four times per cycle to adjust its direction of movement, a rotary engine maintains continuous unidirectional rotational motion, resulting in enhanced stability and reduced vibrations [68]. Additionally, the rotary engine boasts high power density, simple design, compact size, and lightweight construction [911], positioning it as a formidable competitor to reciprocating piston engines in applications such as extended-range vehicles and small drones [1214].

However, the rotary engine faces challenges in achieving the fuel efficiency of reciprocating engines due to its narrow and flat combustion chamber shape, high surface area, and low compression ratio [15,16], thereby limiting its application across various fields [17,18]. To enhance fuel economy and optimize rotary engine performance, researchers worldwide have undertaken numerous studies on engine structure [1922]. Kuo et al. [23] investigated the impact of rotary profile on engine compression flow via numerical analysis, revealing that increasing the shape coefficient K can enhance compression efficiency and mixture formation. Tartakovsky et al. [24] focused on optimizing the rotary combustion chamber structure, starting from the in-cylinder flow field, resulting in improved engine performance. Similarly, Jeng et al. [25] emphasized fuel-air mixing as the starting point, enhancing engine performance by optimizing the rotary combustion chamber structure. Additionally, Wei et al. [2628] observed that the position of the rotary combustion chamber influences flame propagation, thus optimizing engine performance by adjusting its position.

Furthermore, scholars have recognized the potential of spoilers within the rotary combustion chamber to enhance engine performance [2931]. Shi et al. [32] investigated the impact of spoiler positioning on rotary engine combustion, demonstrating that closer proximity to the spark plug increases turbulence velocity and dissipation rate in the spark region, which accelerates mixture combustion and enhances combustion characteristics and emission performance. Furthering this line of research, Shi et al. [33] also developed a novel turbulence-induced blade configuration for a hydrogen-enriched rotary engine, revealing that strategic placement within the rotor chamber significantly boosts turbulence, accelerates combustion, and greatly improves thermal efficiency and emissions performance. Complementing these findings, Kordi and Esfahanian [34] conducted a detailed exploration of the axial vane rotary engine using a quasi-dimensional modeling technique, showing that spoiler variations in the cam blade profiles critically impact dynamic behaviors, thereby affecting the engine's time between overhaul.

However, existing studies on spoilers primarily focus on their structural parameters' impact on engine performance, without elucidating the mechanism by which spoilers affect flow and combustion processes within the cylinder [27,3540]. Therefore, by analyzing the motion patterns of the rotary, this study discovers that volume inconsistencies arise between two cylinder areas separated by the spoiler, influenced by rotary movement. If the material exchange between these areas is not promptly completed, pressure differences arise, forcing the gas to flow from higher-pressure areas to lower-pressure areas. To distinguish this flow phenomenon from the forced flow caused by the spoiler, it is termed “pressure differential flow,” resulting from local pressure variations in the cylinder induced by the combustion chamber, cylinder body, and spoiler structure. To investigate the impact of pressure differential flow on in-cylinder flow and combustion processes, this article develops a mathematical model for in-cylinder flow and a three-dimensional simulation model of the engine. Studies on the structural parameter of spoiler height were conducted, calculating the minimum height required to produce pressure differential flow and verifying it through three-dimensional simulation models. The study analyzes spoiler heights' effects on in-cylinder flow patterns, elucidating the impact of pressure differential flow on ignition delay, maximum cylinder pressure, ignition timing, and combustion stability. This article uncovers a unique phenomenon of high-speed local gas flow within rotary engines and studies its causes and effects on in-cylinder flow and combustion characteristics, providing a theoretical basis and guidance for further optimization of rotary engine spoilers.

2 Theoretical Analysis of Rotary Engine Spoiler's Impact on Cylinder Flow

In contrast to conventional engines that convert piston reciprocating motion into crankshaft rotation, the rotary engine's eccentric rotation results in different volumes of the combustion chamber due to varying angles within the cylinder body [41]. Figure 1 illustrates the rotary area, rotating clockwise. In the diagram, region A represents the area formed by the rotary, cylinder, and spoiler in the rotation direction, while region B denotes the area formed in the opposite rotation direction. The space between the spoiler and the rotary cylinder serves as the flow area connecting areas A and B.

Fig. 1
Schematic representation of the rotary area
Fig. 1
Schematic representation of the rotary area
Close modal

When assessing the volume of different cylinder areas, the cylinder's flow process is initially simplified as follows:

  1. The width of the spoiler is disregarded, and the extension line of the rotary's center to the spoiler's center is taken as the virtual boundary between combustion chambers A and B.

  2. The gas within the cylinder is treated as an ideal gas, with constant temperature during rotary selection.

  3. The effects of the spoilers on the gas are neglected.

  4. The effects on the translation distance are excluded, and both cylinder and rotary profiles are considered theoretical.

  5. The impact of leakage on cylinder quality is disregarded.

  6. Chemical reactions do not occur among the substances within the cylinder.

Figure 2 presents a schematic illustrating the geometric relationship between various areas of the combustion chamber, rotary, and cylinder block at a rotation angle α. In the figure, O represents the rotation center, O denotes the rotary center, OA2 and OA1 are parallel to the x-axis, B1 indicates the intersection of the rotary's rotation direction and the cylinder body, B3 represents the intersection of the rotary in the opposite direction of rotation and the cylinder body, B4 signifies the spoiler center, B2 depicts the intersection of the extension line of the rotary center O and the spoiler center B4 with the cylinder body, and C marks the intersection of OB1 and OB2, A1OO=3A2OB3=α.

Fig. 2
Schematic of geometric relationship
Fig. 2
Schematic of geometric relationship
Close modal
The volume F of the combustion chamber in a rotary engine, characterized by a creation radius of R, an eccentricity e, and a cylinder thickness B, can be represented as
(1)
(2)
(3)
(4)
(5)
where F5 denotes the volume of the combustion chamber. If we consider the spoiler center B4, where the vertical distance between the intersection point B1 of the rotary and the cylinder in the direction of rotation is b, then the angle β between B1O and B2O can be formulated as
(6)
Next, the volume FA of area A can be articulated as
(7)
(8)
(9)
(10)
where F5A represents the volume of the combustion chamber within area A. The volume FB of region B can be defined as
(11)
The point P(γ)=P(x,y) on the theoretical cylinder line can be formulated as
(12)
(13)
The coordinates of points B1, B2, and B3 can be represented as
(14)
(15)
(16)
Subsequently, F1A can be articulated as
(17)
The coordinates of point O(xO,yO) can be represented as
(18)
(19)
The coordinates of point C(xC,yC) can be expressed as
(20)
(21)
By integrating Eqs. (14)(21), we can determine the coordinates of points C, B2, and O, and F6 can be formulated as
(22)
In a coordinate system with the rotary center as the origin and OB3 as the direction of the negative y-axis, the point Q(γ)=Q(x,y) on the theoretical rotary line can be represented as
(23)
(24)
The coordinates of point B4 can be expressed as
(25)
Subsequently, F4A can be expressed as
(26)

By integrating Eqs. (1)(26), we can determine the volume FA of area A and the volume FB of area B.

The flow area S between areas A and B can be expressed as
(27)
where S signifies the expanded flow area due to the spark plug chamber, h denotes the maximum height of the spoiler, δ represents the ratio of the current rotary spoiler height to the maximum height of the spoiler, and g(k) is the combustion chamber profile equation of the current rotary at the spoiler position.
Assuming the quantity of gas in the combustion chamber from intake to closure is M, the average molar mass of the gas is m, and the cylinder temperature during the rotary rotation is T, the theoretical average pressure Pt of the combustion chamber at time t can be formulated as
(28)
where R is a constant linked to the units, and Ft represents the volume of the combustion chamber at time t. It is apparent from the equation that the theoretical average pressure of the combustion chamber solely depends on the chamber's volume. Assuming that the average pressure PAt of area A matches the average pressure PBt of area B at time t, then
(29)
(30)
(31)
(32)
(33)
where FAt and FBt represent the volumes of regions A and B at time t, respectively. Similarly, MAt and MBt denote the gas masses of regions A and B at time t, respectively, which can be formulated as
(34)
(35)
If, after Δt has elapsed, the average pressure PAt+Δt in area A stays consistent with the average pressure PBt+Δt in area B, then the gas masses in the two areas should be
(36)
(37)
Subsequently, during the time Δt, the mass of gas exchanged ΔMΔt between areas A and B would be
(38)
At this time, the critical average flow velocity vΔt at the intersection of areas A and B should be
(39)

If the average flow velocity at the junction of areas A and B at time t falls below vΔt, exchange of material between these two areas may not occur within the required timeframe. When Eqs. (15)(20) are integrated, it becomes evident that this discrepancy will cause the average pressure PAt+Δt in area A and the average pressure PBt+Δt in area B at time t+Δt to no longer be equal, thereby resulting in a pressure difference across different areas of the combustion chamber. As a consequence of this pressure differential, the average velocity at the junction of areas A and B will exceed the theoretical average velocity, subsequently influencing the flow and combustion process within the cylinder. This article discusses the local pressure disparity within the cylinder, arising from the combustion chamber, cylinder body, and spoiler structure, along with ensuring high-speed gas flow within the cylinder, termed pressure differential flow.

3 Establishment of a Rotary Engine Simulation Model

To validate the pressure difference flow phenomenon discussed in Sec. 2, this study employs the renormalization group Kε model to simulate airflow dynamics within the cylinder [4244]. Combustion computations were conducted using the primary reference fuel framework mechanism, comprising 41 components and 124 chemical reactions, coupled with the software for analysis of gas-phase chemistry and aerosol formation in the environment (SAGE) detailed chemical kinetic model [45,46]. A three-dimensional simulation model of a rotary engine was developed to investigate in-cylinder flow and combustion processes at various spoiler heights [47]. This study employs a single-cylinder gasoline rotary engine for investigation and validation, with structural parameters detailed in Table 1. The corresponding solid model is illustrated in Fig. 3.

Fig. 3
Three-dimensional model
Fig. 3
Three-dimensional model
Close modal
Table 1

Structural parameters of the rotary engine

ParametersValues
Creation radius (mm)60
Eccentricity (mm)10
Cylinder thickness (mm)41.2
Rated speed (r/min)8000
Single-cylinder volume (cc)133
Compression ratio10
Ignition advance angle (°CA BTD C)30
Ignition sourceSingle spark plug
ParametersValues
Creation radius (mm)60
Eccentricity (mm)10
Cylinder thickness (mm)41.2
Rated speed (r/min)8000
Single-cylinder volume (cc)133
Compression ratio10
Ignition advance angle (°CA BTD C)30
Ignition sourceSingle spark plug

To mitigate the impact of grid count on calculation outcomes and ensure accuracy, it is essential to conduct a grid independence analysis on the simulation model beforehand, selecting an optimal grid count to minimize computation time. This study assesses alterations in in-cylinder pressure across four varying grid counts, maintaining consistent boundary conditions throughout [48]. Details of the boundary condition settings for the rotary engine are outlined in Table 2. Figure 4 depicts the pressure change curves at different mesh counts. Using adaptive mesh refinement based on a 2-mm grid size, the cylinder pressure variation curve aligns closely with the curve at a 1-mm mesh count, indicating stable computational results. Based on the analysis above, the study opts for a 2-mm adaptive mesh division considering computational efficiency, accuracy, and timing [49]. In this article, the rotary engine is modeled with approximately 200,000 hexahedral cells during the simulation process, as shown in Fig. 5.

Fig. 4
Cylinder pressure variation curves under different grid numbers
Fig. 4
Cylinder pressure variation curves under different grid numbers
Close modal
Fig. 5
Computational grid of the rotary engine
Fig. 5
Computational grid of the rotary engine
Close modal
Table 2

Boundary conditions of the rotary engine

BoundaryTemperature (K)Pressure (MPa)
Air inlet3000.101325
Intake300
Exhaust vent5700.101325
Exhaust duct550
rotary400
Spark plug6240.117210
Spark electrode600
BoundaryTemperature (K)Pressure (MPa)
Air inlet3000.101325
Intake300
Exhaust vent5700.101325
Exhaust duct550
rotary400
Spark plug6240.117210
Spark electrode600

4 Impact Analysis of Spoiler Height on Engine Performance

This section delves into the influence of spoiler height, positioned at the rotary's center, on pressure differential flow, as well as its effect on the cylinder flow and the combustion processes. The maximum spoiler height, measured from the bottom edge of the combustion chamber to the rotary profile, stands at 4.19 mm. Simulations were conducted for nine spoilers of varying heights, namely 100% maximum height (4.19 mm), 95% maximum height (3.97 mm), 90% maximum height (3.77 mm), 85% maximum height (3.56 mm), 80% maximum height (3.35 mm), 75% maximum height (3.14 mm), 70% maximum height (2.93 mm), 65% maximum height (2.72 mm), and 60% maximum height (2.51 mm). Figure 6 illustrates cross-sections of the rotary solid model at 100% maximum height, 80% maximum height, and 60% maximum height, with h representing the spoiler height.

Fig. 6
Entity model: (a) 100% maximum height, (b) 80% maximum height, and (c) 60% maximum height
Fig. 6
Entity model: (a) 100% maximum height, (b) 80% maximum height, and (c) 60% maximum height
Close modal

4.1 Impact Analysis of Spoiler Height on Pressure Difference Flow.

As the spoilers of the nine rotary types mentioned above are all situated at the rotary's center, volume changes across different areas remain consistent. Figures 7(a) and 7(b) display the comparison between volume change curves derived from Eqs. (1)(13) and the actual volumes computed by the simulation model for various regions. Figure 7(c) illustrates a comparative diagram between the theoretical flow area computed using Eq. (14) and the actual flow area at 100% maximum height. Table 3 presents the discrepancy between theoretical calculations and actual results. It is evident from the table that the volume calculated using the flow area method exhibits minimal error compared with the actual volume. The volume calculation exhibits an average error of 1.20% in area A and 0.96% in area B. Additionally, the average error in flow area calculation is 7.42%, providing reliable data for in-cylinder flow analysis and computation.

Fig. 7
Comparison between theoretical calculation and actual results: (a) comparison between theoretical volume and actual volume of area A, (b) comparison between theoretical volume and actual volume of area B,and (c) comparison between theoretical flow area and actual flow area at 100% maximum height
Fig. 7
Comparison between theoretical calculation and actual results: (a) comparison between theoretical volume and actual volume of area A, (b) comparison between theoretical volume and actual volume of area B,and (c) comparison between theoretical flow area and actual flow area at 100% maximum height
Close modal
Table 3

Calculation error

ProjectAverage error (%)
Volume of area A1.20
Volume of area B0.96
Flow area7.42
ProjectAverage error (%)
Volume of area A1.20
Volume of area B0.96
Flow area7.42

Figure 8(a) depicts the volume change rates of various regions as computed from theoretical volume adjustments. Observing both Figs. 7 and 8(a), it becomes apparent that within the range of −40°CA to 20°CA, region A initially experiences a decline in volume followed by an increase, whereas region B demonstrates a swift decrease in volume. This varying trend in regional capacity change will result in inconsistent gas densities across the two regions, thereby generating pressure differences between them and ultimately inducing pressure differential flow among different regions. Unequal alterations in regional volume serve as the foundation for pressure differential flow. Figure 8(b) illustrates a curve representing the variation in rotary flow area at different spoiler heights, plotted against the crankshaft angle and calculated using Eq. (14). From the figure, it is evident that with the rotary's fundamental structural parameters, rotation speed, and combustion chamber shape held constant, the spoiler height directly influences the flow area between various regions, consequently resulting in different pressure differential flows. As the spoiler height increases, the flow area between different regions decreases, leading to a more pronounced generation of pressure differential flow.

Fig. 8
Circulation changes: (a) volume change rate in different regions and (b) theoretical flow area at different spoiler heights
Fig. 8
Circulation changes: (a) volume change rate in different regions and (b) theoretical flow area at different spoiler heights
Close modal

Conversely, the lower the spoiler structure, the weaker the effect, potentially resulting in insufficient generation of pressure differential flow. For rotary at the same spoiler positions, flow area trends remain consistent across various spoiler heights. This trend entails a gradual decrease from −40 °CA to 5 °CA, followed by a rapid increase after 5 °CA. Moreover, at the same angles, the flow area decreases with the increasing spoiler height. This variation in flow area constitutes the primary factor influencing pressure difference flow. An excessively large flow area between different regions may expedite a material exchange, resulting in either no or minimal pressure difference due to negligible air density disparities, thereby weakening pressure differential flow. Conversely, a smaller flow area impedes material exchange between different regions, leading to higher-pressure differences attributable to air density variations, thereby enhancing pressure differential flow. Figure 9 depicts a schematic of the flow area at 5°CA and 15°CA, marked with arrows, showing a rapid increase in the flow area post 5°CA due to the spark plug chamber's influence.

Fig. 9
Schematic of flow area at different crankshaft angles: (a) flow area at 5 °CA and (b) flow area at 15 °CA
Fig. 9
Schematic of flow area at different crankshaft angles: (a) flow area at 5 °CA and (b) flow area at 15 °CA
Close modal

Equation (20) was employed to compute the critical average flow velocity vΔt of the rotary with various spoiler heights, and the outcomes are depicted in Fig. 10. Notably, under identical spoiler positions, the critical average flow velocity of rotary with different spoiler heights follows a consistent trend, steadily rising with increasing spoiler height at equivalent angles. Furthermore, the critical flow velocity of rotary with varying spoiler heights exhibits a similar trend of change, gradually increasing from −40 °CA to 5 °CA. Subsequently, post 5 °CA, the flow area increases due to the influence of the spark plug chamber before decreasing. Additionally, at the same angles, the flow area decreases with increasing spoiler height, as illustrated in the figure. The critical velocity serves as a prerequisite for initiating pressure differential flow. As seen from Eq. (20), the critical velocity and the theoretical maximum flow velocity at the regional boundary collectively determine the occurrence and magnitude of pressure differential flow between different regions. When the disparity between the critical velocity and the theoretical maximum velocity at the regional boundary is negative, it signifies that material exchange can occur effectively, resulting in no pressure variance between different regions and thus no pressure differential flow. When the disparity is positive, it indicates that material exchange is influenced by pressure differential among various regions, as exchange cannot be completed promptly. Moreover, a greater disparity suggests a slower-than-anticipated material exchange, resulting in an increased pressure difference between regions. Consequently, this amplifies the strength of pressure differential flow across different areas.

Fig. 10
Critical speed change curve
Fig. 10
Critical speed change curve
Close modal

Figure 11 illustrates the pressure difference between areas A and B at various spoiler heights, as calculated by the three-dimensional simulation model. Positive values denote greater pressure in area A than in area B, resulting in the airflow from A to B, while negative values indicate the reverse. The black solid line represents the pressure difference without a spoiler structure. Figure 12 depicts the maximum speed at the area A and area B junction for different spoiler heights, with the black solid line denoting speeds without a spoiler. The figure indicates that influenced by the flow area, critical velocity, and theoretical maximum flow velocity at the regional interface, both the pressure difference and flow velocity increase as spoiler height increases at the same angle. A spoiler height less than 75% of the maximum height has a minimal effect on pressure, with its pressure variation area similar to that without the spoiler. The figure also reveals that the trend of the maximum speed curve of the spoiler heights ranging from 100% to 80% of the maximum height is consistent. However, for the spoiler heights between 75% and 60% of the maximum height, the pressure difference is smaller, leading to a reduced pressure difference flow. As a result, the trend in this range is similar to that of the structure without spoiler. This indicates that when the spoiler height exceeds 75% of the maximum height, the cylinder's flow is primarily influenced by the pressure difference between different areas and the forced flow generated by the spoiler. Conversely, when the spoiler height is less than 75% of the maximum height, the cylinder's flow is mainly driven by the forced flow created by the spoiler itself.

Fig. 11
Pressure difference change curve
Fig. 11
Pressure difference change curve
Close modal
Fig. 12
Maximum speed change curve
Fig. 12
Maximum speed change curve
Close modal

4.2 Impact Analysis of Spoiler Height on In-Cylinder Flow and Combustion Process.

Figure 13 presents pressure cloud diagrams for structures with spoilers at 100%, 80%, and 75% of the maximum height spoiler, as well as a structure without a spoiler, at −10 °CA, 0 °CA, and 10 °CA. The figure clearly shows a pressure difference between various areas. When the spoiler height exceeds 75% of the maximum height, the regional pressure difference caused by the spoiler is more pronounced. Conversely, when the spoiler height is below 75% of the maximum height, the regional pressure difference is less significant. Compared with the structure without spoilers, structures with spoilers exhibit lower pressure at the time of ignition (−10 °CA to −5 °CA) and in the ignition area (area A), with the pressure in the ignition area decreasing as the spoiler height increases.

Fig. 13
Pressure cloud diagram of different structures
Fig. 13
Pressure cloud diagram of different structures
Close modal

Figure 14 displays velocity cloud diagrams of the structures with spoilers at 100%, 80%, and 75% of the maximum height, as well as for a structure without a spoiler, at −10 °CA, 0 °CA, and 10 °CA. The figure shows that the flow velocity in the cylinder increases with higher spoiler height. When the spoiler height is more than 75% of the maximum height, the pressure difference flow is more pronounced, and the flow field distribution is more uniform. In this case, the cylinder experiences both pressure difference flow and forced flow caused by the spoiler. When the spoiler height is less than 75% of the maximum height, the pressure difference flow is weaker, and the flow field distribution resembles that of a structure without a spoiler. Here, the forced flow in the cylinder is primarily due to the spoiler. The high-speed area in the figure indicates that the pressure difference flow mainly occurs between the spoiler and cylinder body, following the cylinder body's rotation, with the high-speed area peaking at 10 °CA. The pressure difference flow also influences the position and intensity of the vortex within the cylinder. As the spoiler height increases, the vortex intensity in the spark plug chamber decreases. At the crankshaft angle of 10 °CA, the spark plug chamber of the structure with a spoiler at 100% maximum height shows almost no vortex motion. This occurs because the junction between the spark plug chamber and the cylinder is nearly enveloped by the high-speed airflow moving in the rotation direction, creating a closed area that restricts gas flow. As the spoiler height increases, the vortex in the cylinder gradually shifts from the cylinder body to the junction area between the spoiler and the combustion chamber. Structures with higher spoiler heights have less orderly in-cylinder flow in the early stage of combustion (after 10 °CA) than those with lower spoiler heights.

Fig. 14
Velocity cloud diagram of different structures
Fig. 14
Velocity cloud diagram of different structures
Close modal

Figure 15 illustrates the flame morphology for structures at 100%, 80%, and 75% of the maximum height, as well as a structure without a spoiler, at −10 °CA, 0 °CA, and 10 °CA. The figure shows that when the spoiler height exceeds 75% of the maximum height, the pressure difference flow results in a higher flame propagation speed in the early stage of combustion. However, the propagation speed is inconsistent in all directions, potentially not only enhancing fuel combustion and increasing explosion pressure but also leading to poor combustion stability. When the spoiler height is 75% of the maximum height, the pressure difference flow is weak, resulting in a more uniform flow field that resembles the distribution of the structures without spoilers.

Fig. 15
Flame morphology of different structures
Fig. 15
Flame morphology of different structures
Close modal

4.3 Impact Analysis of Spoiler Height on Performance.

Figure 16 illustrates a comparison between the cylinder pressure, cylinder temperature, and heat release rate curves for structures with spoilers at 100%, 80%, and 75% of the maximum height spoiler, as well as for a structure without a spoiler. The figure reveals that the cylinder pressure and temperature curves for the spoiler structures at 100% and 80% of the maximum height follow a similar trend. Conversely, the curves for the spoiler at 75% of the maximum height and the structure without a spoiler exhibit a consistent trend. This aligns with the earlier analysis regarding the impact of spoiler height on cylinder flow. The heat release rate curve shows significant fluctuations for structures with taller spoilers, indicating poor combustion stability attributed to early-stage flow field effects.

Fig. 16
Performance comparison: (a) cylinder pressure, (b) cylinder temperature, and (c) heat release rate
Fig. 16
Performance comparison: (a) cylinder pressure, (b) cylinder temperature, and (c) heat release rate
Close modal

Figures 17 and 18 depict the curves illustrating the maximum cylinder pressure and the angle of the maximum cylinder pressure relative to the top dead center, respectively, as the spoiler height varies. In the figure, the black solid line parallel to the X-axis represents the value without a spoiler structure. The ignition time, defined as the moment when 0.5% of the fuel's heat is released, and the ignition delay, defined as the duration between ignition and when 10% of the fuel's heat is released, are analyzed in Figs. 19 and 20, respectively. In the figure, the black solid line parallel to the X-axis represents the value without a spoiler structure.

Fig. 17
Curve of maximum cylinder pressure variation with spoiler height
Fig. 17
Curve of maximum cylinder pressure variation with spoiler height
Close modal
Fig. 18
Angle between the maximum cylinder pressure and the top dead center of the crankshaft changes with the spoiler height
Fig. 18
Angle between the maximum cylinder pressure and the top dead center of the crankshaft changes with the spoiler height
Close modal
Fig. 19
Ignition time curve with spoiler height
Fig. 19
Ignition time curve with spoiler height
Close modal
Fig. 20
Ignition delay versus spoiler height
Fig. 20
Ignition delay versus spoiler height
Close modal

It can be observed that the maximum cylinder pressure rises with increasing spoiler height, surpassing that of the structure without a spoiler. This trend is attributed to increased flow velocity within the cylinder, mitigating vortex formation and thus enhancing flame propagation speed and combustion acceleration. As the spoiler height increases, ignition time is delayed, primarily due to lower pressure and higher flow velocity in the ignition area of structures with taller spoilers, inhibiting initial fire core formation. Conversely, the ignition delay angle decreases with increasing spoiler height, as taller spoilers induce greater flow velocity, accelerating flame propagation. Consequently, the angle of the maximum cylinder pressure relative to the top dead center decreases with increasing spoiler height, becoming smaller than that of structures without spoilers when spoiler height exceeds 76% of the maximum height.

Structures with different spoiler heights are categorized into two groups based on cylinder flow characteristics: low spoiler height (between 60% and 75% of the maximum height), where the dominant flow in the cylinder is the forced flow induced by the spoiler, and high spoiler height (between 75% and 100% of the maximum height), where both pressure difference flow and forced flow from the spoiler coexist in the cylinder. The change rates of various performance indicators with spoiler height were calculated and are presented in Table 4.

Table 4

Performance change rate

Project60% of maximum height to 75% of maximum height75% of maximum height to 100% of maximum height
Maximum cylinder pressure0.0085 MPa/%0.0379 MPa/%
Angle of cylinder pressure maximum value from top dead center−0.1239 °CA/%−0.1351 °CA/%
Ignition time0.1366 °CA/%0.1696 °CA/%
Ignition delay−0.1805 °CA/%−0.2202 °CA/%
Project60% of maximum height to 75% of maximum height75% of maximum height to 100% of maximum height
Maximum cylinder pressure0.0085 MPa/%0.0379 MPa/%
Angle of cylinder pressure maximum value from top dead center−0.1239 °CA/%−0.1351 °CA/%
Ignition time0.1366 °CA/%0.1696 °CA/%
Ignition delay−0.1805 °CA/%−0.2202 °CA/%

The table reveals that different in-cylinder flow forms resulting from the spoiler have different effects on the performance indicators. Compared with low spoiler height, the change rate of performance at the ignition time for high spoiler height increases from 0.1366 °CA/% to 0.1696 °CA/%, marking a 24.1% increase. The change rate of ignition delay decreases from −0.1805 °CA/% to −0.2202 °CA/%, a decrease of 21.9%, indicating that pressure difference flow has both a negative impact on delaying ignition timing and a positive effect on accelerating combustion speed in the cylinder.

Considering the cylinder flow analysis, the lower pressure and higher flow velocity in area A at the ignition delay moment lead to delayed ignition, while the higher flow velocity and complex flow field distribution postformation of the fire core accelerate the flame propagation and shorten the ignition delay. Moreover, higher spoiler heights significantly increase the rate of maximum cylinder pressure increase, from 0.0085 MPa/% at low spoiler height to 0.0379 MPa/% at high spoiler height, marking a 345.8% increase. This is primarily due to higher flow velocity and more complex flow field distribution in the cylinder, accelerating flame propagation and fuel combustion. However, the impact of different spoiler heights on the change rate of the angle of maximum cylinder pressure from the top dead center is relatively minor, decreasing by only 9.1% from −0.1239 °CA/% at low spoiler height to −0.1351 °CA/% at high spoiler height.

5 Conclusion

This article has investigated the impact mechanism of the rotary engine spoiler on the flow and combustion process within the cylinder. It theoretically examines the pressure difference flow, which arises from local pressure differences in the cylinder caused by the combustion chamber, cylinder body, and spoiler structure, resulting in localized high-speed gas flow within the cylinder. A corresponding mathematical model is established. Various factors influencing the generation and development of differential pressure flow are theoretically analyzed, and the effects of spoiler height on differential pressure flow, in-cylinder flow, and the combustion process are studied. The study leads to the following conclusions:

  1. In rotary engines, the spoiler divides the cylinder space into two regions, influenced by the rotary's movement patterns. If the material exchange between these regions is not completed promptly, it results in a pressure difference, forcing the gas to flow from high-pressure to low-pressure areas. The spoiler height directly determines the flow area between these regions. At the same angle, the flow area decreases with increasing spoiler height. A larger flow area between regions allows rapid material exchange, leading to small or no pressure differences due to minimal air density variations, thus weakening the pressure differential flow. Conversely, a smaller flow area hinders material exchange, causing a higher-pressure difference due to air density variations, thus enhancing the pressure differential flow. The critical speed, necessary to maintain equal pressure differences between regions, along with the theoretical maximum flowrate at the junction, determines whether pressure differential flow occurs and its intensity. If the difference between the critical speed and the theoretical maximum flowrate is negative, material exchange can be completed without creating a pressure difference, eliminating pressure differential flow. If positive, the delayed material exchange due to the pressure difference slows down more than expected, increasing the pressure difference and intensifying the pressure differential flow. Influenced by the flow area, the critical average flow speed at the same angle increases with the spoiler height.

  2. The flow patterns in the cylinder influenced by the spoiler are primarily divided into forced flow directly caused by the spoiler and pressure differential flow induced by local pressure differences within the cylinder. When the spoiler height exceeds 75% of the maximum height, the cylinder flow is primarily influenced by both the pressure difference flow between different areas and the forced flow created by the spoiler. However, when the spoiler height is less than 75% of the maximum height, the flow is mainly driven by the forced flow from the spoiler. Additionally, pressure difference flow primarily occurs in the region between the spoiler and the cylinder body, with high-speed flow area developing along the cylinder body's rotation direction and peaking at 10 °CA. The vortex intensity in the spark plug chamber decreases as the spoiler height increases, with almost no vortex motion at 10 °CA for the 100% maximum height structure. As the spoiler height increases, the in-cylinder vortex shifts from the cylinder body to the junction area of the spoiler and the combustion chamber. Structures with higher spoiler heights exhibit poorer in-cylinder flow orderliness, consistent flame propagation speed, and stable combustion in the early combustion stage compared with those with lower spoiler heights.

  3. Increasing the spoiler height leads to higher maximum cylinder pressure and ignition time, while the angle of the maximum cylinder pressure from the top dead center of the crankshaft and the ignition delay decrease. The pressure difference flow delays the ignition moment but accelerates combustion speed, significantly increasing the rate of the maximum cylinder pressure rise. However, it has minimal effect on the rate of change of the angle between the maximum cylinder pressure and the top dead center.

Acknowledgment

We thank LetPub for its linguistic assistance during the preparation of this manuscript.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

References

1.
Woon
,
K. S.
,
Phuang
,
Z. X.
,
Taler
,
J.
,
Varbanov
,
P. S.
,
Chong
,
C. T.
,
Klemeš
,
J. J.
, and
Lee
,
C. T.
,
2023
, “
Recent Advances in Urban Green Energy Development Towards Carbon Emissions Neutrality
,”
Energy
,
267
(
23
), p.
126502
.
2.
Yang
,
X.
,
Liu
,
J.
, and
Li
,
M.
,
2024
, “
Transportation Energy Consumption Decreases Despite Urban Sprawl: Evidence From Chinese Household Survey
,”
Cities
,
146
(
12
), p.
104738
.
3.
Feng
,
X.
,
Zhao
,
Y.
, and
Yan
,
R.
,
2024
, “
Does Carbon Emission Trading Policy Has Emission Reduction Effect?—An Empirical Study Based on Quasi-Natural Experiment Method
,”
J. Environ. Manage.
,
351
(
29
), p.
119791
.
4.
Kim
,
M.
,
Won
,
W.
, and
Kim
,
J.
,
2017
, “
Integration of Carbon Capture and Sequestration and Renewable Resource Technologies for Sustainable Energy Supply in the Transportation Sector
,”
Energy Convers. Manage.
,
143
(
11
), pp.
227
240
.
5.
Lin
,
Z.
,
Liu
,
S.
,
Qi
,
Y.
,
Chen
,
Q.
, and
Wang
,
Z.
,
2024
, “
Experimental Study on the Performance of a High Compression Ratio SI Engine Using Alcohol/Ammonia Fuel
,”
Energy
,
289
(
25
), p.
129998
.
6.
Zambalov
,
S. D.
,
Yakovlev
,
I. A.
, and
Maznoy
,
A. S.
,
2022
, “
Numerical Assessment of Interchangeability Potential for Syngas-Fueled Rotary Engine
,”
Energy Convers. Manage.
,
266
(
23
), p.
115805
.
7.
Yang
,
Z.
,
Du
,
Y.
,
Geng
,
Q.
, and
He
,
G.
,
2022
, “
Energy Loss and Comprehensive Performance Analysis of a Novel High-Efficiency Hybrid Cycle Hydrogen-Gasoline Rotary Engine Under Off-Design Conditions
,”
Energy Convers. Manage.
,
267
(
23
), p.
115942
.
8.
Yang
,
J.
,
Wang
,
H.
,
Ji
,
C.
,
Chang
,
K.
, and
Wang
,
S.
,
2023
, “
Investigation of Intake Closing Timing on the Flow Field and Combustion Process in a Small-Scaled Wankel Rotary Engine Under Various Engine Speeds Designed for the UAV Application
,”
Energy
,
273
(
22
), p.
127147
.
9.
Ji
,
C.
,
Shi
,
C.
,
Wang
,
S.
,
Yang
,
J.
,
Su
,
T.
, and
Wang
,
D.
,
2019
, “
Effect of Dual-Spark Plug Arrangements on Ignition and Combustion Processes of a Gasoline Rotary Engine With Hydrogen Direct-Injection Enrichment
,”
Energy Convers. Manage.
,
181
(
16
), pp.
372
381
.
10.
Cai
,
T.
,
Zhao
,
D.
, and
Gutmark
,
E.
,
2023
, “
Overview of Fundamental Kinetic Mechanisms and Emission Mitigation in Ammonia Combustion
,”
Chem. Eng. J.
,
458
(
39
), p.
141391
.
11.
Guan
,
Y.
,
Becker
,
S.
,
Zhao
,
D.
,
Xu
,
J.
,
Shahsavari
,
M.
, and
Schluter
,
J.
,
2023
, “
Entropy Generation and CO2 Emission in Presence of Pulsating Oscillations in a Bifurcating Thermoacoustic Combustor With a Helmholtz Resonator at Off-Design Conditions
,”
Aerosp. Sci. Technol.
,
136
(
11
), p.
108204
.
12.
Kotsiopoulos
,
P.
,
Yfantis
,
E.
,
Lois
,
E.
, and
Hountalas
,
D.
,
2001
, “
Diesel and JP-8 Fuel Performance on a Petter AV1 Diesel Engine
,”
39th Aerospace Sciences Meeting and Exhibit
,
Reno, NV
,
Jan. 8–11
, p.
500
.
13.
Dutczak
,
J.
,
2015
, “
Heavy Fuel Engines
,”
Combust. Engines
,
163
(
4
), pp.
34
46
.
14.
Wu
,
H.
,
Wang
,
L.
,
Wu
,
Y.
,
Sun
,
B.
,
Zhao
,
Z.
, and
Liu
,
F.
,
2019
, “
Spray Performance of Air-Assisted Kerosene Injection in a Constant Volume Chamber Under Various in-Cylinder GDI Engine Conditions
,”
Appl. Therm. Eng.
,
150
(
14
), pp.
762
769
.
15.
Attard
,
W. P.
,
Blaxill
,
H.
,
Anderson
,
E. K.
, and
Litke
,
P.
,
2012
, “
Knock Limit Extension With a Gasoline Fueled Pre-Chamber Jet Igniter in a Modern Vehicle Powertrain
,”
SAE Int. J. Engines
,
5
(
3
), pp.
1201
1215
.
16.
Wang
,
C.
,
Zhang
,
F.
,
Wang
,
E.
,
Yu
,
C.
,
Gao
,
H.
,
Liu
,
B.
,
Zhao
,
Z.
, and
Zhao
,
C.
,
2019
, “
Experimental Study on Knock Suppression of Spark-Ignition Engine Fueled With Kerosene via Water Injection
,”
Appl. Energy
,
242
(
22
), pp.
248
259
.
17.
Chen
,
L.
,
Raza
,
M.
, and
Xiao
,
J.
,
2017
, “
Combustion Analysis of an Aviation Compression Ignition Engine Burning Pentanol-Kerosene Blends Under Different Injection Timings
,”
Energy Fuel
,
31
(
9
), pp.
9429
9437
.
18.
Ning
,
L.
,
Duan
,
Q.
,
Wei
,
Y.
,
Zhang
,
X.
,
Yu
,
K.
,
Yang
,
B.
, and
Zeng
,
K.
,
2019
, “
Effects of Injection Timing and Compression Ratio on the Combustion Performance and Emissions of a Two-Stroke DISI Engine Fuelled With Aviation Kerosene
,”
Appl. Therm. Eng.
,
161
(
15
), pp.
114
124
.
19.
Lu
,
Y.
,
Pan
,
J.
,
Fan
,
B.
,
Otchere
,
P.
,
Chen
,
W.
, and
Cheng
,
B.
,
2019
, “
Research on the Application of Aviation Kerosene in a Direct Injection Rotary Engine—Part 1: Fundamental Spray Characteristics and Optimized Injection Strategies
,”
Energy Convers. Manage.
,
195
(
17
), pp.
519
532
.
20.
Yang
,
J.
,
Ji
,
C.
,
Wang
,
S.
,
Wang
,
D.
,
Ma
,
Z.
, and
Zhang
,
B.
,
2018
, “
Numerical Investigation on the Mixture Formation and Combustion Processes of a Gasoline Rotary Engine With Direct Injected Hydrogen Enrichment
,”
Appl. Energy
,
224
(
19
), pp.
34
41
.
21.
Pan
,
J.
,
Ma
,
J.
,
Li
,
J.
,
Liu
,
H.
,
Wei
,
J.
,
Xu
,
J.
,
Zhu
,
T.
,
Zhang
,
H.
,
Li
,
W.
, and
Pan
,
J.
,
2022
, “
Influence of Intake Port Structure on the Performance of a Spark-Ignited Natural Gas Engine
,”
Energies
,
15
(
22
), p.
8545
.
22.
Fan
,
B.
,
Qi
,
X.
,
Pan
,
J.
,
Wu
,
X.
,
Fang
,
J.
,
Lu
,
Q.
, and
Zhang
,
Y.
,
2023
, “
Research on the Hydrogen Injection Strategy of a Turbulent Jet Ignition (TJI) Rotary Engine Fueled With Natural Gas/Hydrogen Blends
,”
Int. J. Hydrogen Energy
,
48
(
52
), pp.
20041
20058
.
23.
Zou
,
R.
,
Yang
,
W.
, and
Zhang
,
Y.
,
2024
, “
Ignition of Position on Combustion Process of X-Type Rotary Engine
,”
Internal Combustion Engine Engineering
,
45
(
3
), pp.
50
58
.
24.
Tartakovsky
,
L.
,
Baibikov
,
V.
,
Gutman
,
M.
,
Veinblat
,
M.
, and
Reif
,
J.
,
2012
, “Simulation of Wankel Engine Performance Using Commercial Software for Piston Engines,” SAE Tech. Paper No. 2012-32-0098.
25.
Jeng
,
D.-Z.
,
Hsieh
,
M.-J.
,
Lee
,
C.-C.
, and
Han
,
Y.
,
2013
, “The Numerical Investigation on the Performance of Rotary Engine With Leakage, Different Fuels and Recess Sizes,” SAE Tech. Paper No. 2013-32-9160.
26.
Wei
,
F. B.
,
Feng
,
P. J.
,
Xian
,
L. Y.
,
Rui
,
C.
, and
Di
,
W.
,
2013
, “
Effect of Pocket Location on Combustion Process in Natural Gas-Fueled Rotary Engine
,”
Appl. Mech. Mater.
,
316–317
(
27
), pp.
73
79
. www.scientific.net/AMM.316-317.73
27.
Fan
,
B. W.
,
Pan
,
J. F.
,
Pan
,
Z. H.
,
Tang
,
A. K.
,
Zhu
,
Y. J.
, and
Xue
,
H.
,
2015
, “
Effects of Pocket Shape and Ignition Slot Locations on the Combustion Processes of a Rotary Engine Fueled With Natural gas
,”
Appl. Therm. Eng.
,
89
(
8
), pp.
11
27
.
28.
Ji
,
C.
,
Su
,
T.
,
Wang
,
S.
,
Zhang
,
B.
,
Yu
,
M.
, and
Cong
,
X.
,
2016
, “
Effect of Hydrogen Addition on Combustion and Emissions Performance of a Gasoline Rotary Engine at Part Load and Stoichiometric Conditions
,”
Energy Convers. Manage.
,
121
(
11
), pp.
272
280
.
29.
Amrouche
,
F.
,
Varnhagen
,
S.
,
Erickson
,
P. A.
, and
Park
,
J. W.
,
2016
, “
An Experimental Study of a Hydrogen-Enriched Ethanol Fueled Wankel Rotary Engine at Ultra Lean and Full Load Conditions
,”
Energy Convers. Manage.
,
123
(
12
), pp.
174
184
.
30.
Su
,
T.
,
Ji
,
C.
,
Wang
,
S.
,
Shi
,
L.
, and
Cong
,
X.
,
2018
, “
Effect of Ignition Timing on Performance of a Hydrogen-Enriched n-Butanol Rotary Engine at Lean Condition
,”
Energy Convers. Manage.
,
161
(
14
), pp.
27
34
.
31.
Chen
,
W.
,
Pan
,
J.
,
Fan
,
B.
,
Liu
,
Y.
, and
Peter
,
O.
,
2017
, “
Effect of Injection Strategy on Fuel-air Mixing and Combustion Process in a Direct Injection Diesel Rotary Engine (DI-DRE)
,”
Energy Convers. Manage.
,
154
(
13
), pp.
68
80
.
32.
Shi
,
C.
,
Zhang
,
P.
,
Ji
,
C.
,
Di
,
L.
,
Zhu
,
Z.
, and
Wang
,
H.
,
2022
, “
Understanding the Role of Turbulence-Induced Blade Configuration in Improving Combustion Process for Hydrogen-Enriched Rotary Engine
,”
Fuel
,
319
(
27
), p.
123807
.
33.
Shi
,
C.
,
Zhang
,
Z.
,
Ji
,
C.
,
Li
,
X.
,
Di
,
L.
, and
Wu
,
Z.
,
2022
, “
Potential Improvement in Combustion and Pollutant Emissions of a Hydrogen-Enriched Rotary Engine by Using Novel Recess Configuration
,”
Chemosphere
,
299
(
25
), p.
134491
.
34.
Zheng
,
Renchuan
,
Gong
,
Erlei
, and
Li
,
Jianzhong
,
2024
, “
Thermal Cycle Analysis of Radial Flow Internal Combustion Wave Rotary Engine Based on Power Ratio Factor
,”
Prop. Technol.
,
45
(
09
), pp.
57
64
.
35.
Wang
,
S.
,
Sun
,
Y.
,
Yang
,
J.
, and
Wang
,
H.
,
2024
, “
Effect of Excess Air Ratio and Ignition Timing on the Combustion and Emission Characteristics of the Ammonia-Hydrogen Wankel Rotary Engine
,”
Energy
,
302
(
25
), p.
131779
.
36.
Feng
,
Z.
,
Pan
,
J.
,
Li
,
P.
,
Fan
,
B.
,
Nauman
,
M.
, and
Yang
,
W.
,
2024
, “
Numeral Study of the Mixture Formation and Combustion Process of an Ammonia/Hydrogen Rotary Engine
,”
Int. J. Hydrogen Energy
,
69
(
6
), pp.
1469
1480
.
37.
Huo
,
S.
,
Fan
,
B.
,
Xu
,
L.
,
Yin
,
G.
,
Pan
,
J.
,
Yang
,
W.
,
Li
,
W.
, and
Wu
,
Y.
,
2024
, “
Combined Effect of Cylinder Shape and Turbulence Blade on the Combustion Performance of a Turbulent Jet Ignition Rotary Engine Using Hydrogen/Natural Gas Blends
,”
Int. J. Hydrogen Energy
,
61
(
5
), pp.
513
527
.
38.
Li
,
W.
,
Fan
,
B.
,
Jiang
,
P.
,
Liu
,
W.
,
Pan
,
J.
,
Huo
,
S.
,
Wu
,
Y.
, and
Lu
,
Q.
,
2024
, “
Exploring the Coupling Mechanism Between Piston Structure and Turbulent Jet Ignition Mode in Hydrogen Mixed Natural Gas Rotary Engines
,”
Int. J. Hydrogen Energy
,
63
(
6
), pp.
905
917
.
39.
Hege
,
J. B.
,
2006
,
The Wankel Rotary Engine: A History
,
McFarland
,
New York
.
40.
Mitianiec
,
W.
,
2015
, “
Modelling and Simulation of Working Processes in Wankel Engine With Direct Hydrogen Injection System
,”
Combust. Engines
,
161
(
2
), pp.
42
52
.
41.
Georgios
,
Z.
,
2005
,
Mathematical and Numerical Modeling of Flow and Combustion Processes in a Spark Ignition Engine
,
Department of Applied Mathematics
,
Hong Kong, China
, p.
32
.
42.
Meng
,
H.
,
Ji
,
C.
,
Wang
,
S.
, and
Yang
,
J.
,
2023
, “
A Review: Centurial Progress and Development of Wankel Rotary Engine
,”
Fuel
,
335
(
29
), p.
127043
.
43.
Fan
,
B.
,
Pan
,
J.
,
Tang
,
A.
,
Pan
,
Z.
,
Zhu
,
Y.
, and
Xue
,
H.
,
2015
, “
Experimental and Numerical Investigation of the Fluid Flow in a Side-Ported Rotary Engine
,”
Energy Convers. Manage.
,
95
(
9
), pp.
385
397
.
44.
Shi
,
C.
,
Ji
,
C.
,
Wang
,
S.
,
Yang
,
J.
,
Ma
,
Z.
, and
Meng
,
H.
,
2020
, “
Potential Improvement in Combustion Behavior of a Downsized Rotary Engine by Intake Oxygen Enrichment
,”
Energy Convers. Manage.
,
205
(
18
), p.
112433
.
45.
Farouk
,
T. I.
,
Xu
,
Y.
,
Avedisian
,
C. T.
, and
Dryer
,
F. L.
,
2017
, “
Combustion Characteristics of Primary Reference Fuel (PRF) Droplets: Single Stage High Temperature Combustion to Multistage ‘Cool Flame’ Behavior
,”
Proc. Combust. Inst.
,
36
(
2
), pp.
2585
2594
.
46.
Liu
,
Y. D.
,
Jia
,
M.
,
Xie
,
M. Z.
, and
Pang
,
B.
,
2012
, “
Enhancement on a Skeletal Kinetic Model for Primary Reference Fuel Oxidation by Using a Semidecoupling Methodology
,”
Energy Fuels
,
26
(
12
), pp.
7069
7083
.
47.
Spreitzer
,
J.
,
Zahradnik
,
F.
, and
Geringer
,
B.
,
2015
, “Implementation of a Rotary Engine (Wankel Engine) in a CFD Simulation Tool With Special Emphasis on Combustion and Flow Phenomena,” SAE Tech. Paper No. 2015-01-0382.
48.
Qi
,
X.
,
Yang
,
Y.
,
Tian
,
L.
,
Wang
,
Z.
, and
Zhao
,
N.
,
2022
, “
A Parallel Methodology of Adaptive Cartesian Grid for Compressible Flow Simulations
,”
Adv. Aerodyn.
,
4
(
1
), p.
21
.
49.
Lee
,
M.
,
Park
,
G.
,
Park
,
C.
, and
Kim
,
C.
,
2020
, “
Improvement of Grid Independence Test for Computational Fluid Dynamics Model of Building Based on Grid Resolution
,”
Adv. Civil Eng.
,
2020
(
30
), pp.
1
11
.