Optimal Stroke Path for Reciprocating Heat Engines

By testing piston motion in reciprocating heat engines as a control variable, one could find piston trajectories, different from the conventional near sinusoidal motion that should increase power production. *is results from minimizing frictional losses. *e purpose of this study is to determine piston trajectories that are optimal for noncombustion strokes in reciprocating engines, in the sense of minimizing frictional dissipation and hence maximizing efficiency and power. *e optimal piston traces for noncombustion strokes are determined by using a combination of optimal control theory and models for the thermodynamic irreversibilities. Hence, the results are germane to external combustion engines and to the noncombustion strokes of internal combustion engines. *e optimal piston traces or trajectories obtained here can be viewed as some of the building blocks from which optimal overall cycles can be constructed.

One possible way to increase the power delivery of reciprocating heat engines is to vary the piston trajectory relative to its conventional near sinusoidal motion [17].e idea of treating piston motion as a control variable in an attempt to minimize thermodynamic losses and hence maximize power production was raised in [17,18], in which attempts were made to find the optimal piston trajectory for Otto and Diesel cycle engines.e exercise was performed by using optimal control theory in concerts with models for the key thermodynamic irreversibilities [17,18].An important finding in these works was that the potential improvement in the efficiency of the reciprocating engines is not negligible and could be as high as 10-15% of current engine efficiencies.e improvements stem from lowering frictional dissipation, and heat leaks in the internal combustion engine cycles considered.ese earlier analyses, however, adopted oversimplified models for the influence of combustion processes on engine dynamics [20,21].Two other limitations were that (1) the heat generated due to friction was modeled as being totally transferred to the engine cooling system, namely, none of that resulted in the heating of the engine working fluid and (2) the only type of friction considered was the rubbing friction of a piston well-lubricated cylinder surface, where power dissipation goes as the square of piston speed.
e aim of this study is to determine piston trajectories that are optimal for noncombustion strokes in reciprocating engines, in the sense of minimizing frictional dissipation and hence maximizing efficiency and power.In this study, we determine the optimal piston trajectories for the strokes of reciprocating heat engines, in the sense of maximizing power production.However, we restrict our analyses to noncombustion strokes due to the complex nature of modeling combustion and its influence on engine performance.e complex models dictate pure numerical solutions.Among the benefits of doing this is the capability of the analytical solution to enable more explicit and transparent results.Hence, our results are directly applicable to external combustion engines and to the noncombustion strokes of internal combustion engines.
Although we restrict the study to noncombustion strokes, the same methods of optimal control theory could be used to optimize combustion strokes as was done in [17,18,19].
e two key irreversibilities modeled are friction and heat leak, and a range of friction sources are considered: mechanical and/or fluid friction.In addition, our results account for frictional dissipation heating the engine working fluid.
e piston trajectories we determine could hence be viewed as some of the building blocks from which one can calculate the optimal piston motion for various engine cycles.In addition, one would have to calculate the fraction of the total cycle time allotted to each stroke to tailor the results to the particulars of any engine cycle under consideration.
ese are the cycle-specific calculations.e achievable maximum power with these optimal piston trajectories will be compared to the power that can be attained with conventional near sinusoidal piston motion.Furthermore, sensitivity studies on important engine parameters such as compression ratio, maximum piston acceleration, type of friction (i.e., functional dependence on piston speed), and the degree to which frictional dissipation heats up the engine working fluid are performed.e potential improvements in engine power for the strokes analyzed here shown to be of the order of few percent.

Formulation of the Problem and Modeling Assumptions
Figure 1 is a schematic of the system analyzed here: a piston moving inside a cylinder of fixed, given cross-sectional area A. e piston moves along the vertical axis, and the volume V of the working fluid can change from a maximum value of V max to a minimum value of V min , with the engine compression ratio r defined as r � V max /V min .Piston motion as a function of time during the stroke is described by the time evolution of the gas volume V(t), or, equivalently, piston position x(t) � V(t)/A.Piston velocity is described by _ x(t) � dx(t)/dt.For the convenience of analysis in what follows, piston velocity is defined as positive for a given stroke.e steady-state operation is only analyzed: namely, start-up and turnoff transients are not considered.e general picture could also include a regenerator inside the cylinder, as well as two-piston cylinders, as illustrated in Figure 2. Stirling and related cycles are pertinent examples.In that case, we adopt the approximation suitable for current high-quality regenerators of a regeneration effectiveness of close to 100%.Namely, the ratio of the heat exchanged in the regenerator to the heat stored in the regenerator is essentially unity [22][23][24].However, as the gas passes through the regenerator, fluid friction losses will occur and are accounted for in our model.e ideal-gas approximation is adequate for most engine operating conditions and permits closed analytical solutions to be derived.e formalism described below can also be used for real gas behavior, but results would then have to be generated numerically [17,18,[20][21][22][23][24].
e types of noncombustion engine strokes included in our analysis are compression, expansion, and constant volume.e two former types of strokes are found in all engine cycles, and the latter is relevant for Stirling cycles, among others.e heat leak is assumed to occur with constant thermal conductance κ and according to a linear heat transfer law.Viewing κ as a constant is usually a good approximation due to the relatively large heat transfer areas in most cylinders.Hence, the changes in κ that are due to piston motion have a negligible effect.e rate of heat leak Q leak is then given by where T wall is the temperature of the cylinder wall and T gas is the temperature of the working fluid (ideal gas here), both of which are taken to be spatially uniform.e rate at which work is dissipated as friction is taken to be comprised of a static term and a dynamic term and is given by where α 0 and α 1 are constants for the static and dynamic terms, respectively; the exponent m characterizes the dominant type of dynamic friction; and recall that piston velocity _ x is defined as positive.It is important to note that the units of α 1 change with m.For example, m � 2 would describe rubbing friction for the piston against a well-lubricated cylindrical surface; and m � 3 would pertain to fluid friction when turbulence is induced [25], which might be relevant for modeling losses in engines with internal exchangers, such as Stirling cycles [22][23][24].e results will be derived for a general m, and illustrative examples will be presented for m � 2 and 3.
In its most general form, the optimization problem can be stated as follows.We aim to maximize the useful work w over a stroke of given, fixed time τ: where we can express _ w � dw/dt as a function f 1 of three variables: piston position x, piston velocity _ x, and the internal energy E of the gas and is given by e expression for _ w will be the difference between the pressure-volume (mechanical p − V) work, if any, and the frictional dissipation noted in equation ( 2) is given by where δ � −1 for a compression branch, δ � 0 for a constant volume branch, and δ � 1 for an expansion branch.R is the universal gas constant; C v is the specific heat at constant volume; and the ideal gas law has been used for the first term on the right-hand-side of equation ( 5) (for example, in the analyses of [17,18], α 0 � 0 and m � 2).It is important to note that, in equations ( 2)-( 5), the absolute value of the velocity is used in the calculations.
A constant volume stroke need not trivially mean a stationary piston.For example, in the Stirling and related cycles, which are double-piston engines for each cylinder, both pistons can move together at the same speed such that no net change occurs in gas volume, yet there are frictional dissipation and heat leak on these strokes (Figure 2).It is for this reason that constant-volume strokes are included in the analysis and refer to strokes where piston motion occurs.e dynamic constraint for the time evolution of E, as a function f 2 of x, _ x, and E, is given by which in this case is given by where β is a dimensionless number between zero and unity that signifies the fraction of heat dissipated as friction that heats up the engine working fluid (for example, in the analyses of [17,18], β � 0) and n is the number of gas moles in the cylinder.We will also need to incorporate a realistic upper bound, a max , on piston acceleration € x ≡ d 2 x/dt 2 , which is given by e problem can be solved with optimal control theory [26] once the specific modeling assumptions are expressed.We will present solutions for fully externally dissipative friction (β � 0) and fully internally dissipative friction (β � 1) (all realistic systems span the intermediate cases).

Method of Solution
Following the procedures of optimal control theory [26], we define a modified Lagrangian function L as given by with all the terms on the right-hand side of equation ( 9) being given in equations ( 4)- (7).λ is the Lagrange multiplier and, due to the dynamic constraints, is not a constant but depends on time [26].e two control variables are taken as x and E, for which the Euler-Lagrange equations are, respectively, given by ese equations can be arranged in the two coupled equations and are given by Although the initial state of our system is known, the energy at the end of the stroke is not known.erefore, an additional boundary condition is required [26] and is given by which is equivalent to the following boundary condition as given by λ(t � τ) � 0.
In summary, we know the initial conditions and the final conditions are given by Modelling and Simulation in Engineering Now, we have a set of ordinary differential equations and boundary conditions that can be solved to yield the optimal motion x(t) for producing the power of given stroke, or, equivalently, minimizing dissipated work for a particular stroke.In general, the set of differential equations is not linear and must be solved numerically.However, as will now be argued, realistic engine operating conditions correspond to parameter regimes for which the equations can be solved analytically and in closed form, which makes the solution more transparent and amenable to sensitivity studies.

Negligible Role of Heat Leak
For most practical reciprocating engines, on noncombustion strokes, heat leaks are negligible compared to the friction losses and/or p − V work [17,18,[20][21][22][23][24].Just as one specific illustrative case, consider a Stirling engine with the compression ratio of 3, κ � 100 W/K, hydrogen ideal working fluid with R/C v � 0.4, 0.015 kg of gas in the cylinder, a stroke time of 0.008 s, and frictional losses characterized by m � 2 in equation (2).It then turns out that the difference in calculated engine power, between the cases of ignoring heat leak and taking it fully into account, is negligible.Even if κ is as much as ten times greater, the error introduced by ignoring it altogether is less than 0.5% (not percentage points).We have confirmed this by generating the exact solution numerically for the above realistic κ value, and then comparing it against the approximate analytic solution, the latter taken at κ � 0. erefore, we proceed by generating solutions for the case of κ � 0; namely, the only dissipation is due to friction.Similar arguments supporting this approximation could be found in [17,18].
Static friction does not affect the optimal solution; only the dynamic friction term can influence the optimal path.To see this in the governing equations, integrate equations ( 5) and (7), with negligible heat leak (κ � 0) to obtain (16) where E f denotes the energy at the end of the stroke, and where r denotes the compression ratio.e path dependence for frictional losses is contained in the dynamic friction terms (α 1 terms) only.

Parameters for Conventional Piston Motion
As a basis for later comparison, consider the conventional piston motion, which is approximated here as having a sinusoidal velocity profile (see Figure 3 in the next section): where x i � x(t � 0) is the piston position at the beginning of the stroke and τ is the stroke time.e maximum acceleration for the sinusoidal motion is given by which will later serve as a lower bound for the maximum acceleration in the optimal trajectories.In typical reciprocating engines, frictional dissipation results in an efficiency or power loss of 20% [17,18,[20][21][22][23][24].Since the optimal trajectories determined below are based on minimizing frictional losses only and should reduce them by 5-75% for realistic engine parameters, and improvements in engine efficiency should be in the range of 1-15%.

Solution with Externally Dissipative Friction
Consider the case of externally dissipative friction; namely, β � 0 in equation (7). is problem was solved in [17] for the case of m � 2 (equation ( 2)), and the optimal solution being that piston velocity should be constant.e solution is unaffected by different values of the exponent m.However, the magnitude of the frictional losses, and the reduction therein relative to conventional piston motion, does depend on the exponent m.
Strictly constant piston velocity, however, would require infinite acceleration at the start of the stroke and infinite deceleration at the end of the stroke.As shown in [17], incorporating the constraint of a given maximum acceleration yields a solution in which the piston is accelerated at its maximum permissible value at the start of (3) optimal motion, externally dissipative friction, a * � 2; (4) optimal motion, externally dissipative friction, a * � 1. e abscissa is time divided by stroke time (t/τ), and the ordinate is piston velocity relative to its maximum value for the conventional trajectory.a * � maximum piston acceleration relative to maximum piston acceleration for sinusoidal motion.4 Modelling and Simulation in Engineering the stroke, followed by constant velocity, and completed by deceleration, at the maximum permissible value and at the end of the stroke (Figure 3) [17].In the language of optimal control theory, an additional constraint is accounted for (8) with which the Lagrangian is linear in piston acceleration.Hence, piston acceleration has a "bang-bang" solution [26], whereby it assumes its maximum value, and the sign changing (acceleration or deceleration) depends on whether one is at the beginning or at the end of the stroke.e optimal piston motion for externally dissipative friction is described by the velocity profile [17]: where t 1 is the switching time and is given by

Reduced Losses with Externally Dissipative Friction
e work dissipated due to friction can now be evaluated for conventional piston motion and optimal piston motion.It is convenient to define nondimensional acceleration, a * , as the maximum piston acceleration for the optimal path, relative to maximum piston acceleration for sinusoidal motion, and its value is given by e optimal switching time t 1 (from the constant acceleration to the constant velocity branch) can then be expressed as e optimal switching time for this case depends on maximum permissible acceleration only (Figure 4).
As noted above in Section 4, in equations ( 16) and ( 17), differences in frictional dissipation due to different piston motions arise solely from dynamic friction, i.e., α 1 term in equation (2).Hence, the problem of evaluating potential improvements in engine efficiency that stem from optimal piston motion reduces to evaluating the reduction in dynamic frictional dissipation only.
e ratio of dynamic frictionally dissipated work for conventional piston motion to that for the optimal piston motion can easily be shown as where Γ denotes the gamma function (gamma(•) � exp(gammaln(•))).Note that, with equation ( 23), the righthand side of equation ( 24) depends on a * only.
Figure 5 shows this relative improvement as a function of nondimensional maximum piston acceleration, for which two values of the frictional exponent m are considered (m � 2, 3, equation ( 2)).In the limit of unbounded acceleration, this improvement is 23% for m � 2 and 64% for m � 3, with most of these gains already being reached at only twice the maximum acceleration of conventional piston motion.In Figures 3, 5, and 6, our choice of nondimensional variables results in plots that do not have an explicit dependence on the compression ratio.
Figure 6 shows the relative improvement as a function of the frictional exponent (externally dissipative friction) m in the range 2 ≤ m ≤ 3 and for three values of nondimensional maximum piston acceleration: unbounded, 2, and 1.

Solution with Internally Dissipative Friction
is case corresponds to solving the governing equations with β � 1 (and κ � 0) so that the optimal piston trajectory is described by the following differential equation: On constant volume strokes, of the type often encountered in Stirling cycles, with δ � 0, the optimal piston trajectory calls for constant velocity.When bounded acceleration is accounted for, the solution of the optimal piston motion is identical to the case with externally dissipative friction, which is analyzed in Section 6.
In the more general case of expansion or compression strokes, the solution to equation ( 25) for piston velocity is where a * is the maximum piston acceleration relative to maximum acceleration for conventional trajectory (equation (22)).e optimal switching times for ending the initial constant acceleration branch, t 1 , and commencing the final constant deceleration branch, t 2 , are not symmetrically located about midstroke, as they were for the case of externally dissipative friction.ese optimal switching times are determined by the boundary conditions that (1) piston velocity vanishes at the end of the stroke ( _ ey can be obtained numerically as the solutions to the following two coupled nonlinear equations: ) . ( Illustrations of optimal motion and conventional piston motion, when friction is dissipated internally, are presented in Figures 7 and 8. ese include sensitivity study to key system parameters, such as compression ratio, maximum piston acceleration, and frictional exponent in equation (2).

Reduced Losses with Internally Dissipative Friction
Recall that to determine the potential improvement due to optimal piston motion, one needs to evaluate only the reduction in dynamic frictional dissipation.Dynamic frictional dissipation for the conventional and the optimal trajectories can be evaluated via numerical integration or analytically with the evaluation of hypergeometric functions.Graphical illustrations of the magnitude of this improvement and how it depends on key system parameters are presented in Figures 9 and 10.Although the shapes of the ordinates approximately appear constant, their numerical values change between 1.2 and 1.3 for m � 2 and between 1.65 and 1.75 for m � 3. e changes are nonlinear in behavior.2)).Maximum piston acceleration is expressed relative to its maximum value for conventional sinusoidal motion.Maximum piston acceleration is expressed relative to its maximum value for conventional sinusoidal motion. 6 Modelling and Simulation in Engineering e potential improvement in engine efficiency is slightly larger than that for the externally dissipative case.
e influence of compression ratio is modest, whereas the functional form of frictional losses, and piston maximum acceleration, can have a marked effect on the relative improvement in engine power.

Summary and Conclusions
One prospective method to increase the efficiency and power of reciprocating engines is to reduce frictional losses by modifying piston motion.e purpose of this paper (study) is to determine piston trajectories that are optimal for noncombustion strokes in reciprocating engines, in the sense of minimizing frictional dissipation and hence  (3) optimal motion, internally dissipative friction, a * � 2; (4) optimal motion, internally dissipative friction, a * � 1. e abscissa is time divided by stroke time (t/τ), and the ordinate is piston velocity relative to its maximum value for the conventional trajectory.Modelling and Simulation in Engineering maximizing efficiency and power.e potential improvements in overall engine efficiency have also been evaluated and should lie in the range of around 2-10%, depending on how severe frictional losses were prior to the modification of piston motion.e optimized engine strokes that we have calculated could serve as building blocks in the construction of optimal cycles in external combustion of optimal cycles in the external combustion engines and for noncombustion strokes in internal combustion engines.
Several qualifications are in order.First, no experimental verification has yet been attempted.Second, specific functional forms for frictional losses have been assumed.However, these forms appear to cover the behavior of real engines [20][21][22][23][24], with friction stemming from static, surface rubbing friction, and fluid turbulence contributions.ird, altering piston motion engenders changes in standard crankshafts, which in turn can change frictional losses external to the engine.ese could improve or worsen overall engine efficiency, depending on how such crankshafts would be designed and built.is issue is not analyzed here.
Optimal piston motion turns out to be sensitive to where the heat generated by frictionally dissipated work goes: either to an external cooling system (externally dissipative), or to heating the engine working fluid (internally dissipative), or, in reality, somewhere between these two extremes.We have derived closed-form analytic solutions for optimal piston motion for both cases and compared them to conventional sinusoidal piston paths.
e improvement in engine efficiency for the optimal paths, relative to conventional piston motion, turns out to be mildly sensitive to compression ratio, with a significant effect arising from the functional form of frictional losses and from maximum piston acceleration.Only when piston accelerations in excess of those achievable with conventional sinusoidal motion can be reached will the improvements in engine efficiency be substantial.Fortunately, these improvements in engine efficiency increase rapidly with maximum attainable acceleration for values just above those of conventional motion and then asymptote (Figures 5 and  10).
e mechanisms to achieve the optimal piston motion and the experimental verification of the optimal piston motion are two good research ideas for future research.

Figure 3 :
Figure 3: Piston velocity plotted against time, for one engine (expansion) stroke for (1) conventional sinusoidal motion;(2) optimal motion, externally dissipative friction, unbounded acceleration; (3) optimal motion, externally dissipative friction, a * � 2; (4) optimal motion, externally dissipative friction, a * � 1. e abscissa is time divided by stroke time (t/τ), and the ordinate is piston velocity relative to its maximum value for the conventional trajectory.a * � maximum piston acceleration relative to maximum piston acceleration for sinusoidal motion.

Figure 4 :
Figure 4: Switching time as a function of a * in the case of externally dissipative friction.

Figure 5 :
Figure5: Ratio of work dissipated as friction with conventional piston motion relative to optimal piston motion, plotted against maximum piston acceleration, for two values of frictional exponent m (equation (2)).Maximum piston acceleration is expressed relative to its maximum value for conventional sinusoidal motion.

Figure 6 :
Figure6: Ratio of work dissipated as friction with conventional piston motion relative to optimal piston motion, plotted against frictional exponent m (equation (2)), for three values of maximum piston acceleration: unbounded acceleration, a * � 2, and a * � 1. Maximum piston acceleration is expressed relative to its maximum value for conventional sinusoidal motion.

3 Figure 8 :
Figure 8: Example of sensitivity of optimal piston motion to compression ratio r, for internally dissipative friction.Maximum nondimensional piston parameter a * � 1. Frictional exponent m � 2. Plots are for (1) conventional motion;(2) optimal motion with r � 3; (3) optimal motion with r � 8. e abscissa is time divided by stroke time (t/τ), and the ordinate is piston velocity relative to its maximum value for the conventional trajectory.

Figure 9 :
Figure9: to compression ratio of reduction in frictional losses.Calculations are for unbounded acceleration.Ordinate is dynamic frictional losses for conventional piston motion relative to that for optimal piston motion, for internally dissipative friction.Abscissa is the compression ratio.e two curves are for different frictional exponents m and show the marked dependence on m.

Figure 10 :
Figure10: Sensitivity of reduction in frictional losses to maximum piston acceleration.Ordinate is dynamic frictional losses for conventional piston motion relative to that for optimal piston motion, for internally dissipative friction.Abscissa is maximum piston acceleration relative to its maximum value for the conventional trajectory.e two curves are for different frictional exponents m and show the marked dependence on m.
Ratio of frictional losses of sinusoidal motion and frictional losses of optimal motion vs. m for different values of a *