Next Article in Journal
Do Investment Strategies Matter for Trading Global Clean Energy and Global Energy ETFs?
Next Article in Special Issue
Inkjet 3D Printed MEMS Electromagnetic Multi-Frequency Energy Harvester
Previous Article in Journal
Steady State Modeling and Performance Analysis of a Wind Turbine-Based Doubly Fed Induction Generator System with Rotor Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design Optimization of Tubular Heat Exchangers for a Free-Piston Stirling Engine Based on Improved Quasi-Steady Flow Thermodynamic Model Predictions

1
Department of Mechanical System Design Engineering, Seoul National University of Science and Technology, Seoul 01811, Korea
2
Department of Mechanical and Automotive Engineering, Seoul National University of Science and Technology, Seoul 01811, Korea
*
Authors to whom correspondence should be addressed.
Energies 2022, 15(9), 3326; https://doi.org/10.3390/en15093326
Submission received: 18 March 2022 / Revised: 27 April 2022 / Accepted: 27 April 2022 / Published: 3 May 2022
(This article belongs to the Special Issue Vibration-Based Energy Harvesters)

Abstract

:
This paper presents the design optimization of a heat exchanger for a free-piston Stirling engine (FPSE) through an improved quasi-steady flow (iQSF) model and a central composite design. To optimize the tubular hot heat exchanger (HHX) design, a design set of central composite designs for the design factors of the HHX was constructed and the brake power and efficiency were predicted through the iQSF model. The iQSF model is improved because it adds various heat and power losses based on the QSF model and applies a heat transfer model that simulates the oscillating flow condition of an actual Stirling engine. Based on experimental results from the RE-1000, an FPSE developed by Sunpower, the iQSF model significantly improves the prediction error of the indicated power from 66.9 to 24.9% compared to the existing QSF model. For design optimization of the HHX, the inner diameter and the number of tubes with the highest brake power and efficiency were determined using a regression model, and the tube length was determined using the iQSF model. Finally, the brake output and efficiency of FPSE with the optimized HHX were predicted to be 7.4 kW and 36.4%, respectively, through the iQSF analysis results.

1. Introduction

Due to increasing concerns about the excessive use of fossil fuels, air pollution, and global warming, environmental regulations imposed on conventional heat engines are becoming increasingly stringent. Increasing population trends and more modern lifestyles worldwide are, however, demanding higher levels of energy. To meet these conflicting demands, the European Union (EU) took the initiative to raise the ratio of renewable energy usage beyond 20% out of the total energy consumption by 2020 by setting up the EU 2009 renewable energy action plan [1].
The Stirling engine is an external combustion engine that operates on a closed cycle. It consists of two isochoric heating/cooling and two isothermal compression/expansion processes based on the Stirling cycle.
The most important advantage of the Stirling engine is that it can use any type of heat source for fuel. Owing to this advantage, the Stirling engine can produce electricity using eco-friendly thermal energy sources such as biogas, biomass, and waste heat. In particular, the free-piston Stirling engine (FPSE), first proposed by Beale [2,3], replaces the slider-crank mechanism of the mechanical Stirling engine with a spring connected to the pistons. By reducing both the engine weight and the mechanical friction, it has a semi-permanent lifespan, leading to technological advances in the Stirling engine.
The Stirling engine consists of two types of parts; the first type is the moving parts, the displacer and power piston linked to the spring/load, and the second consists of the thermal energy transfer parts, specifically the hot heat exchanger, the regenerator, and the cold heat exchanger. The flow of the working fluid is induced between two enclosed spaces by the reciprocating motion of the displacer. While passing through the heat exchangers and regenerator, the temperature of the working fluid changes. The compression and expansion of the working fluid led by the temperature change result in the reciprocating motion of the power piston and, consequently, mechanical power.
Stirling engine analysis models for predicting performance outcomes have been developed by various researchers, and Dyson [4] classified Stirling engine analysis models into five major categories.
The zeroth-order model, first proposed by Beale [5], is an analysis model that introduces Beale’s number, which predicts engine performance outcomes through correlations between operating conditions and performance based on experimental data of commercial engines. This model is the simplest, but the performance prediction accuracy is quite low.
The first-order model is an isothermal model that was initially proposed by Schmidt [6]. The isothermal model assumes that each space of the engine is isothermal; it calculates the volume and pressure of the engine space based on the ideal gas state equation and then predicts the pressure-volume (PV) power. Performance predictions by this model remain limited due to the isothermal assumption of the idealized engine space, but this type is widely used by those who design initial engines due to its simplicity.
The second-order model is an adiabatic model first proposed by Finkelstein [7] and Urieli [8]. The adiabatic model divides the engine space into two working spaces (compression space and expansion space) and three heat exchanger spaces (cooler, regenerator, and heater) and assumes that the working and heat exchanger spaces are adiabatic and isothermal, respectively. Each space defines a thermodynamic state quantity based on the energy equation and predicts the engine performance by numerical integration according to the time change. In particular, the non-ideal adiabatic model proposed by Urieli [8] is a quasi-steady flow (QSF) model that predicts heat transfer characteristics by assuming the flow in the heat exchanger as a steady-state flow. It is combined with various loss models [9,10,11] for reasonable performance predictions. It is widely used in optimization research due to its advantages of accuracy and fast analysis speeds. Based on these advantages, studies of various second-order models have recently been conducted [12,13].
Third-order models use the finite difference method to predict engine performance outcomes numerically by dividing the engine space into several connected one-dimensional nodes, defining each node using a differential equation that follows the laws of conservation of energy and conservation of mass. This model was performed by Finkelstein [14], Urieli [15], and Berchowitz [16]. In particular, Gedeon developed Sage [17,18,19], a program dedicated to the Stirling cycle based on the third-order model.
A fourth-order analysis refers to a model that divides and analyzes the engine space into multidimensional microspaces using computational fluid dynamics (CFD). It was studied by Mahkamov [20], Tew [21], and Wilson [22]. This model provides detailed analysis results of factors such as the flow velocity, flow pattern, temperature, and pressure distribution of the working fluid inside the engine. Despite these advantages, a disadvantage is that it requires a new model design each time depending on the engine shape. It also requires long calculation times for the analysis.
This study presents the design optimization of tubular heat exchangers of the FPSE developed by Nexergy lab at Seoultech [23] using the iQSF model and the development and verification of the iQSF model. The iQSF model is based on the second-order QSF model, and the third-order model is applied to the regenerator space, which is divided into several one-dimensional nodes. The model is verified through the experimental results of the RE-1000 engine developed by NASA’s Lewis Research Center [24]. The optimized design of the heat exchanger also targets the engine developed by this institution. The optimized design of the heat exchanger aims to increase performance outcomes by reducing losses. The design parameters of the heat exchanger are the tube diameter, length, and number of tubes. Each simulation case was determined from the design of experiment (DOE), which is a rapid and efficient method for determining the effects of variables [25]. The simulation results were fitted with a regression model [26], and the effects of the design variables were studied. The optimized design variables of the heat exchanger, which exhibited good performance and low losses, were obtained.

2. Description of the Developed FPSE

Figure 1 shows a cross-sectional view of a 3 kW FPSE design developed by Nexergy lab at Seoultech [27]. This engine was developed to use the exhaust gas waste heat below 500 °C from a 100 kW diesel generator as a heat source. The hot heat exchanger (HHX) and cold heat exchanger (CHX) are the shell-tube types, and 600 tubes with an inner diameter of 1.7 mm are applied to maximize the heat exchange area. However, using too many tubes will likely lead to welding defects, and in fact, many welding defects occurred here. Therefore, in this study, instead of the existing shell-tube type HHX, a tube-type HHX with dozens of tubes is used, and the design of the HHX tubes is optimized. In addition, the heat source temperature is increased to 600 °C to increase the engine performance. Table 1 presents the design parameters of the developed FPSE.

3. Range of the Design Parameters

3.1. Tubular Heat Exchanger Literature Review

To determine the range of the HHX design parameters, Table 2 shows the tube design parameters of several kilowatt Stirling engines employing tubular HHX developed by various institutions around the world. The design parameters of the tubular HHX are the tube inner diameter, tube length, and the number of tubes.
RE-1000 [24], a 1 kW-class FPSE developed in NASA’s Lewis facility in 1986, has an HHX consisting of 34 tubes whose inner diameter and length are ~2.4 mm and 183 mm, respectively. The SHARP (Stirling helium advanced reliable prototype) engine [28], 2 kW-class FPSE, was developed by Sunpower Inc. in 1989. The HHX of this engine is a tube type, and the inner diameter, length, and the number of tubes are 2.2 mm, 140 mm, and 60, respectively. The engine system was built by connecting two engine heads in parallel. These two engines produce 4 kW of electrical output, but the system electrical output is 3 kW, excluding the power consumption in the power requirements of the packaging (coolant pump, burner, battery charging). Thorsen [29] and Carlsen’s engines [30] are both crank-driven engines developed by the Technical University in Denmark, with an electrical output of 3 kW and 10 kW, respectively. The number of tubes in both engines is 24, which is relatively few compared to other engines. In particular, the tube inner diameter of Carlsen’s engine is 8 mm, which is more than twice that of other engines. The GPU-3 engine [31] is well known as a representative crank-driven Stirling engine, and its engine test results have been used as a reference to verify various crank-driven engine prediction models [32,33,34]. The 7 kW class Stirling engine developed by Qnergy Inc. is commercialized [35]; it is characterized by well-optimized configurations and component designs compared to the other prototype-level engines. Qnergy’s Stirling engine was presented by estimating the tube design parameters based on the released patent [31] because there was no data on which the design parameters were disclosed.
To summarize the range of tube parameters for all engines, as shown in Table 2, the tube inner diameter, number of tubes, and length are 2.2–8 mm, 24–80, and 140–365 mm, respectively. The tube parameter range in this study for optimizing the HHX design of the developed FPSE was determined by considering the shape of the developed FPSE while generally including the parameter ranges above. Each parameter was categorized into three levels, as presented in Table 3.
The characteristics of both objective functions are such that they offer greater brake power and efficiency of the engine. To find the optimal point-of-response values from the selected design factors, three-level design points were chosen for the three design variables, as shown in Table 4, using the central composite design (CCD). In addition to CCD, DOE that can be applied includes full factorial and Taguchi design. The efficient CCD method was used since the full factorial increases the analysis cases. Taguchi design is mainly used for the robust design that can minimize the influence of noise factors. Since the analysis model is used for the optimization of this paper, there is no need to consider the influence of noise factors, so the CCD method is used.

3.2. Design Constraints

Figure 2 shows the shape structure of the cylinder head and the heat exchanger tube. The cylinder head maintains the shape used in the previously developed engine. The HHX is U-shaped and is divided into an inner tube, an outer tube, and a bent tube. A number of tubes are arranged circumferentially in the cylinder head. The diameter of the outer tube array is fixed at 95 mm because it is a flow passage connected to the regenerator, and the diameter of the inner tube array changes according to the radius r b of the bent tube. The tube length l t is defined by Equation (1):
l t = r b π + 2 b + a  
where a = 0.0084 r b 2 + 1.4684 r b + 1.4434 .
Here, b is the straight length of the inner tube and a is the length difference between the inner and outer tubes, which can be expressed as a second-order polynomial that increases with the length of r b . For the tube bend radius r b , a minimum value of 15 mm is required for production. The minimum value of the tube length determined as a design variable is 100 mm, meaning that the maximum value of r b according to Equation (1) is approximately 25 mm when b is 0 mm.
In addition, the maximum number of applicable tubes varies according to the tube bend radius; this is determined by Equation (2) below:
n t . m a x = 2 π 95 2 r b / d o + l g
where d o is the outer diameter of the tube and l g is the gap between the tubes. The tube thickness is fixed at 0.5 mm in consideration of the charge pressure of the engine, and the gap between the tubes remains at least 1 mm to simplify the machining of the holes.
Figure 3 shows the maximum applicable number of tubes according to the inner tube diameter and tube bending radius. As r b increases, the number of tubes decreases for the same tube inner diameter. The maximum number of tubes that can be arranged in a given cylinder head shape is approximately 80, but the number of tubes here increases to 130 in order to identify the performance trend of the engine according to the number of tubes as well as the optimal design of the HHX.

4. Improved Quasi-Steady Flow (iQSF) Model

4.1. Modeling Procedure

The iQSF model is an improved model of a quasi-steady flow (QSF) model that is based on a second-order analysis [8], except for the regenerator part, where a third-order analysis is applied in iQSF to minimize errors by considering the effects of various losses in the analyses. In the iQSF model, the current FPSE can be largely divided into seven space cells: the wall of the cooler, the matrix of the regenerator, the wall of the heater, the expansion space, the bounce space for the displacer, the compression space, and the bounce space for the power piston. For the calculations of the thermodynamic characteristics of the system, appropriate assumptions were made for the heat transfer analyses. The working spaces including expansion and compression spaces were assumed to be adiabatic. The bounce spaces were modeled as closed systems with the assumption of an adiabatic process. In the heater and cooler regions, the heat addition to and rejection from the working fluid occur by forced convection under the isothermal wall temperature conditions. The regenerator is modeled as consecutive finite elements where convection heat transfer occurs between the matrix of the regenerator and the working fluid. The additional assumptions needed to solve the energy equations and to determine the performance of the FPSE are as follows:
  • The working fluid is considered an ideal gas.
  • Changes in the kinetic and potential energies are negligible.
  • The incompressible effect of the working fluid is negligible.
  • The mass of the working fluid inside the engine is constant.
  • All flows of the working fluid inside the engine are assumed to be periodic steady-state flows.
  • The motion of DP and PP follows a sinusoidal waveform with a constant velocity, amplitude, and phase.
  • The temperature distribution of the working fluid along the radial direction of the engine space is negligible.
Based on these assumptions, the governing equations and losses for each space cell of the engine were derived as shown in Table 5 and Table 6, respectively, and were used for analyzing the thermodynamic characteristics of the FPSE by propagating the dynamics of one space cell to another.
The friction factor of the heat exchanger is based on Moody’s chart, and friction factor correlation equations ignoring the surface roughness are used [36]. In addition, the convective heat transfer coefficient was defined using the Reynolds analogy and the J-factor analogy, which are heat transfer correlation equations that depend on the Reynolds number in the pipe flow. The correlation equations of the heat exchanger are shown in Table 7 [36,37,38]. The friction coefficient and convective heat transfer coefficient of the regenerator were referenced from a study that experimentally measured the flow velocity, temperature, and pressure drop of the flow inside the regenerator [39]. The correlation equations of the regenerator are shown in Table 8.
Structural conduction loss implies heat conduction loss through the external structure of the engine, specifically the engine case and the outer surface of the cylinder. The regenerator efficiency loss refers to the loss of thermal energy that cannot be reheated and is wasted during the reheating process of the regenerator. The appendix gap loss refers to the heat conduction loss through the gap and structure of the displacer and leakage loss implies the heat loss caused by the leakage of operating fluids. The seal friction loss and flow friction loss refer to losses caused by the viscous friction of the working fluid in the seal structure and caused while passing through the heat exchangers, respectively. The gas spring hysteresis loss refers to the loss due to the formation of hysteresis caused by the fluid viscosity distribution.
By integrating the governing equations and the losses, the pressure and temperature distribution in each space cell can be obtained, as shown in Figure 4, which also defines the variables. Descriptions of the temperature and pressure distribution are given in detail in the references [8].
Based on the quantity of state determined through the equations in Table 5 and Table 6, the expansion work, compression work, and indicated power generated through the Stirling cycle are defined as follows:
W ˙ e = p e d V e d θ ,   W ˙ c = p c d V c d θ
P i n d = f W ˙ e + W ˙ c
Here, f is the operating frequency. The brake power applying the power losses, consisting of the gas spring hysteresis loss, heat exchanger friction loss, and seal friction loss, can be defined as the follows:
P b r a k e = P i n d f W ˙ g s h , j + W ˙ h x f , i + W ˙ s f , j d θ
Based on this, the indicated efficiency and brake efficiency of the engine reflecting the heat conduction loss through the structure can be defined as follows:
η i n d = P i n d f Q ˙ h + Q ˙ c o n d . c a s e + Q ˙ c o n d . c y l d θ
η b r a k e = P b r a k e f Q ˙ h + Q ˙ c o n d . c a s e + Q ˙ c o n d . c y l d θ

4.2. Validity of the Model

For verification of the iQSF model, test result #1011 of the RE-1000 engine developed by NASA was used [24]. Test result #1011 is the standard experimental condition suggested in the literature, and the analytical models developed by many researchers were compared with the results of experiment #1011 [8,40]. The RE-1000 engine is a beta-type structure in which two pistons are located on the same axis, as shown in Figure 5. The DP and PP bounce spaces are connected to DP and PP, respectively, and act as a gas spring of the piston. The dashpot load acts as a mechanical load on PP. The heat exchanger consists of a tube-type heater, a fin-type cooler, and a woven mesh-type regenerator. Table 9 shows the design parameters of RE-1000 and the operating conditions of #1011.
Figure 6 shows a flow chart of the calculation process of the iQSF model. First, the engine design parameters and operating conditions serve as input. At this time, the step size of the iQSF model numerical analysis was determined to be one cycle per hundred cycles. The initial values of the iQSF model variables are determined, and numerical integration is conducted on the ordinary differential equations of the iQSF model using the fourth-order Runge–Kutta method. At this time, the analysis is repeated based on the convergence standard when the net heat transfer value of the regenerator during one cycle is less than 0.1% of the indicated output. The net heat transfer amount of the regenerator to zero during one cycle means that the Stirling cycle has reached a periodic steady state. At the same time, convergence is checked every one cycle, and in the absence of convergence, the initial value is reset and calculations are repeated to derive the result.
Figure 7 shows the predicted heat input and indicated power of the RE-1000 engine according to the number of finite elements when the regenerator space is divided into one-dimensional finite elements. As the number of the regenerator elements increases from 10 to 100, the heat input and the indicated power constantly converge from 7523 W to 4615 W and 1085 W to 1286 W, respectively. When the number of regenerator elements changed from 80 to 100, it was judged that the analysis results converged because the indicated power changed by less than 0.5%. Based on these results, all analyses were performed by dividing the regenerator space into 100 finite elements.
Figure 8 shows an energy flow diagram of the RE-1000 engine through the iQSF model analysis. The heat input is 4615 W, the heat including the loss of 3329 W is removed, and the indicated power of 1286 W is generated. The indicated power generates 1126 W of brake power, and the remaining 160 W is dissipated as power loss. The seal leakage loss is the largest, accounting for 55% of the total loss. This is presumed to be the cause of the large pressure difference between the engine space and the bounce space due to the high charge pressure. The heat exchanger friction loss occurs mostly in the regenerator with 60.5 W for the regenerator, 15.8 W for the heater, and 2.5 W for the cooler. The analysis found that a large amount of viscous friction arises due to the narrow flow passage of the regenerator.
Table 10 shows the prediction results of the Formosa model, the QSF model, and the iQSF model under identical RE-1000 test conditions [8,24,40]. First, the iQSF model predicts the test results with high accuracy compared to the existing QSF model. In particular, the predicted indicated power and efficiency of the iQSF model show high accuracy with less than 25% error compared to the test. In the Formosa results, the predicted indicated power showed a low error of 15.6% compared to the test results. However, the heat in and rejection outcomes are predicted with smaller values than those in the test results. The analysis shows that the Formosa model predicts the performance based on the idealized first-order analysis, which explains the lower predicted values compared to the test results.

5. Regression Analysis

The results of the iQSF model analysis for the CCD case are listed in Table 11. Figure 9 shows effects of terms for (a) the brake power and (b) efficiency. The parameters with the strongest effects with regard to brake power and efficiency were the product of the tube’s inner diameter ( d i ) and the number of tubes ( n t ); the p-value is less than 0.05, indicating a statistically significant effect on the brake power and efficiency. This is because n t · d i term is the factor that increases the dead volume the most. As a single term, performance is greatly affected in the order of n t , d i , and l t . In particular, l t has a relatively low effect on performance compared to n t and d i . The rest of the terms except for the n t · d i term also affect the performance, but the p-value is greater than 0.05, meaning that statistically, it is not considered to have a significant effect on the performance. However, in order to identify the performance trend according to the determined tube design parameters, the following regression models for the brake power and efficiency were established by considering all terms:
P b r a k e W = 2074 + 1776 d i + 2.6 l t + 89.2 n t 102.5 d i 2 0.0498 l t 2 0.177 n t 2 + 1.76 d i · l t 9.04 d i · n t + 0.046 l t · n t
η b r a k e % = 10.54 + 9.59 d i 0.009 l t + 0.445 n t 0.504 d i 2 0.000159 l t 2 0.0009 n t 2 + 0.00555 d i · l t 0.0451 d i · n t + 0.000167 l t · n t
R-square represents the fitness of the regression model for the interpreted data. R-square for the indicated power and the efficiency calculated a high fit above 0.85.

6. Analysis of the Influence of Variables

We analyzed the effects of the design variables on (a) the brake power and (b) the efficiency, as shown in Figure 10. Comprehensively examining the effects of tube design parameters on the brake power and efficiency, the optimal value for the brake power and efficiency depends on the tube’s inner diameter. It is shown that when the tube’s inner diameter increases, the heat transfer area, which increases the performance, and the dead volume, which decreases the performance, increase simultaneously, forming an optimum value. An increase in the number of tubes significantly increases the heat transfer area compared to the dead volume, thus increasing both the brake power and the efficiency. As the tube length increases, the brake power finds its optimum value, but the efficiency decreases. This occurs because when the tube length increases, the heat transfer area increases, leading to the advantage of absorbing a considerable amount of heat. Thus, the brake power can increase, but the efficiency decreases because the friction loss increases.
Next, in order to identify the performance for the n t   · d i term, which has the greatest influence on the brake power and efficiency, Figure 11 shows the predicted brake power and efficiency according to the inner tube diameter and the number of tubes using the regression model. Figure 11a shows the brake power for tube lengths of 100 mm, 200 mm, and 300 mm. Overall, the brake power increases when the tube inner diameter is small and with a high number of tubes, as this design can significantly increase the heat transfer area compared to the dead volume. Additionally, when both the tube inner diameter and the number of tubes are increased, the heat exchanger dead volume becomes very large and thus the brake power is reduced. The brake power increases in the order of tube lengths of 200 mm, 100 mm, and 300 mm; a similar trend is shown in Figure 10.
However, the realizable area in the cylinder structure is r b 25 mm. In all graphs in Figure 11, when r b is 15 mm, 20 mm, and 25 mm, the realizable area is indicated by a line. In the graph showing a case with a tube length of 200 mm, the highest brake power of 7707 W (case 17) was predicted when the inner diameter of the tube was 8.5 mm and the number of tubes was 39.
Figure 11b shows the brake efficiency for tube lengths of 100 mm, 200 mm, and 300 mm. The brake efficiency trend according to the tube inner diameter and the number of tubes is similar to that of the brake power, but for the tube length, the maximum brake efficiency is shown at 100 mm, in contrast to the brake power result. With a limited r b value, the highest possible efficiency is 38.55% (case 16), and the tube design parameters at this time are the tube inner diameter of 8.2 mm, number of tubes of 40, and tube length of 100 mm.
The tube inner diameter and the number of tubes that generate maximum brake power and efficiency are very similar at about ~8.5 mm and ~39, respectively, while the tube length is different at 100 mm and 200 mm, respectively. Because the tube length can change freely under the limited r b design condition, design optimization for the tube length was performed through an iQSF analysis.
Figure 12 shows (a) the brake output and (b) efficiency with respect to changes in the tube length when the tube inner diameter is 8.5 mm and the number of tubes is 39 through the iQSF analysis. The brake efficiency has slight fluctuations in value due to numerical error; thus, a clear trend was confirmed through a cubic polynomial fitting line. First, looking at the brake power, when the tube length is 200 mm, the result of the iQSF analysis was approximately 7400 W, but in the regression model, the brake power is 7707 W, which causes a discrepancy of about 4%, stemming from the error of the regression model. As mentioned above, because the R-squared of the regression model is about 0.9, there may be an error of 10% or less. The brake efficiency also has a regression error of less than 10%. The brake power calculated through the iQSF analysis is the highest at a tube length of 260 mm, and the brake efficiency rate is the highest at a tube length of 180 mm. The brake power variation was nearly identical within 0.3% in the tube length range of 220–320 mm, and the brake efficiency variation was similar within 0.5% in the tube length range of 140–220 mm, meaning that the tube length was set to 220 mm, which can increase both the output and the efficiency. Finally, the tube inner diameter, tube length, and the number of tubes were determined to be 8.5 mm, 220 mm, and 39 tubes, respectively. The optimal design condition is defined as case 19.

7. Loss Analysis for the Tube Design Parameters

Figure 13 shows the portion of the total loss for the optimal HHX design conditions (inner tube diameter 8.5 mm, tube length 220 mm, number of tubes 39) in the optimal case 19 determined above. The regenerator efficiency loss accounted for the largest portion at 42.3%, followed by the appendix gap loss and DP seal leakage loss at 28.1% and 14.6%, respectively. This generally shows values similar to those of the parasitic loss of a well-designed Stirling engine [8].
Figure 14 shows both the major loss for the design condition that generates the maximum brake power for cases 16, 17, and 18 in Figure 11 and the major loss in the worst case 2 with the lowest brake power in Table 11. The major losses are the DP seal leakage loss, heat exchanger flow friction loss, regenerator efficiency loss, and appendix gap loss, and the remaining losses are not indicated due to their small values. The tube inner diameter and the number of tubes for the worst case are 3 mm and 30, respectively, lower than in other cases, meaning that the flow friction loss generated inside the HHX is very high at approximately 3500 W. Among the input parameters of the iQSF analysis, two pistons motions operate under a fixed condition such that the volume flow rate of all cases is identical, whereas the flow area inside the HHX is very narrow, resulting in a large flow friction loss due to the high gas flow rate. Moreover, as the flow resistance that arises when the working gas moves from the expansion space to the HHX space increases, the flow rate of the working gas delivered to the compression space through the DP gap also increases. Accordingly, the DP seal leakage loss is also very high at approximately 7500 W. As a result, the flow rate of the working gas through the regenerator is reduced such that the regenerator inefficiency loss is relatively small compared to the other cases. The other three design conditions are similar in terms of the magnitude of the major losses due to the well-designed HHX.

8. Conclusions

This paper presents the development and verification of an improved quasi-steady flow (iQSF) model and undertakes design optimization of the heat exchanger of the FPSE developed by Nexergy lab at Seoultech using the iQSF model.
The iQSF model additionally applied heat losses and power losses based on the QSF model. In addition, the heat transfer model of the oscillating flow condition, which is the flow of the working fluid in an actual Stirling engine, was applied, and FDM was applied to the regenerator space to improve the prediction of the internal flow of the regenerator. Model verification was performed using the experimental results of the RE-1000 engine, an FPSE developed by Sunpower. The iQSF model significantly improved the model prediction error from 66.9 to 24.9% compared to the existing QSF model.
In order to optimize the design of the HHX of the engine developed by our institution using the iQSF model, the selected design factors of the HHX were the tube inner diameter, the tube length, and the number of tubes, and the design range was determined through a literature review of kilowatt Stirling engines. To find the optimal point-of-response values from the selected design factor, the HHX design set of CCD was constructed and the brake power and efficiency were predicted through the iQSF model. The parameters with the greatest effects on the brake power and efficiency were the product of the tube inner diameter and the number of tubes, and single factors that affected the performance were the number of tubes, tube inner diameter, and tube length in that order.
To optimize the HHX design, using the regression model derived from response surface regression, the inner diameter and number of tubes with the highest brake power and efficiency were determined to be 8.5 mm and 39, respectively. The tube length with high design freedom was determined to be 220 mm using the iQSF model. As a result, the brake power and efficiency of the FPSE with the optimized HHX were predicted to be 7.4 kW and 36.4%, respectively, through the iQSF model.
Finally, the major losses of three cases with good performance and the worst cases were analyzed through the iQSF model. It was confirmed that as the inner diameter and number of tubes decreased, the flow resistance increased such that not only the flow friction loss but also the DP seal leakage loss increased significantly. Through this, it was confirmed that the flow friction loss and the DP seal leakage loss are very important design considerations in the design of the HHX of the FPSE.
A limitation of the iQSF model used in this paper is that the piston motion must be constrained. The iQSF model is suitable for use in a mechanical type Stirling engine in which the piston motion is fixed, but there is a limit to using it in the FPSE in which the piston motion is changed. Therefore, in the future, we plan to develop a model that can predict the piston motion by combining the iQSF model and the equation of motion of the piston of the FPSE.

Author Contributions

Conceptualization, K.S.; methodology, K.S. and T.Y.K.; software, D.-J.K. and Y.P.; validation, D.-J.K. and Y.P.; formal analysis, D.-J.K.; investigation, D.-J.K. and T.Y.K.; writing—original draft preparation, D.-J.K. and T.Y.K.; writing—review and editing, D.-J.K.; supervision, K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT (NRF-2018R1D1A1B07047868).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Symbols
Aarea [m2]
Awtwetted area [m2]
alength difference between the inner and outer tubes [m]
bstraight length of the inner tube [m]
cpspecific heat at constant pressure [J/kg·K]
cvspecific heat at constant volume [J/kg·K]
dhydraulic diameter [m]
ditube inner diameter [m]
dotube outer diameter [m]
foperating frequency [Hz]
fffanning friction factor
gmass flux [kg/m2·s]
Henthalpy [J]
hconvective heat transfer coefficient [W/m2·K]
kthermal conductivity [W/m·K]
Llength [m]
lttube length [m]
mmass [kg]
NuNusselt number
ntNumber of tubes
Ppower [W]
ppressure [Pa]
Qheat [J]
Rgas constant [J/kg·K]
ReReynolds number
rbradius of the bent tube [m]
StStanton number
Ttemperature [K]
Twwall temperature [K]
ttime [s]
ugas velocity [m/s]
Vvolume [m3]
Wwork [J]
Xamplitude [m]
Greek letters
γspecific heat ratio
εeffectiveness
ηefficiency
θcrank angle [°]
κthermal conductivity [W/m·K]
μviscosity [Pa·s]
Subscripts
appappendix gap
bbounce space
ccompression space
ckcompression space to cooler space
condconduction
cylcylinder
ddisplacer
eexpansion space
frfree flow
gshgas spring hysteresis
hheater
heheater space to expansion space
hxheat exchanger friction
indindicated
kcooler
krcooler space to regenerator space
ppower piston
rregenerator
rehreheat
shshuttle
slseal leakage
whheater wall
wkcooler wall

References

  1. Jackson, J. Ensuring emergency power for critical municipal services with natural gas-fired combined heat and power (CHP) systems: A cost–benefit analysis of a preemptive strategy. Energy Policy 2007, 35, 5931–5937. [Google Scholar] [CrossRef]
  2. Beale, W.T. Free Piston Stirling Engines-Some Model Tests and Simulations; 0148-7191, SAE Technical Paper; SAE: Warrendale, PA, USA, 1969. [Google Scholar]
  3. Beale, W.T. Stirling Cycle Type Thermal Device. U.S. Patent US3552120A, 5 January 1971. [Google Scholar]
  4. Dyson, R.; Wilson, S.; Tew, R. Review of computational stirling analysis methods. In Proceedings of the 2nd International Energy Conversion Engineering Conference, Providence, RI, USA, 16–19 August 2004; p. 5582. [Google Scholar]
  5. Walker, G.; Senft, J.R. Free Piston Stirling Engines; Springer: Berlin/Heidelberg, Germany, 1985. [Google Scholar]
  6. Schmidt, G. The theory of Lehmann’s calorimetric machine. Z. Des Ver. Dtsch. Ing. 1871, 15, 98–112. [Google Scholar]
  7. Finkelstein, T. Generalized Thermodynamic Analysis of Stirling Engines; 0148-7191, SAE Technical Paper; SAE: Warrendale, PA, USA, 1960. [Google Scholar]
  8. Urieli, I.; Berchowitz, D.M. Stirling Cycle Engine Analysis; Taylor & Francis: Boca Raton, FL, USA, 1984. [Google Scholar]
  9. Parlak, N.; Wagner, A.; Elsner, M.; Soyhan, H.S. Thermodynamic analysis of a gamma type Stirling engine in non-ideal adiabatic conditions. Renew. Energy 2009, 34, 266–273. [Google Scholar] [CrossRef]
  10. Babaelahi, M.; Sayyaadi, H. Simple-II: A new numerical thermal model for predicting thermal performance of Stirling engines. Energy 2014, 69, 873–890. [Google Scholar] [CrossRef]
  11. Hosseinzade, H.; Sayyaadi, H. CAFS: The Combined Adiabatic–Finite Speed thermal model for simulation and optimization of Stirling engines. Energy Convers. Manag. 2015, 91, 32–53. [Google Scholar] [CrossRef]
  12. Sowale, A.; Anthony, E.J.; Kolios, A.J. Optimisation of a quasi-steady model of a free-piston stirling engine. Energies 2019, 12, 72. [Google Scholar] [CrossRef] [Green Version]
  13. Torres García, M.; Carvajal Trujillo, E.; Vélez Godiño, J.A.; Sánchez Martínez, D. Thermodynamic model for performance analysis of a Stirling engine prototype. Energies 2018, 11, 2655. [Google Scholar] [CrossRef] [Green Version]
  14. Finkelstein, T. Computer analysis of Stirling engines. Advances in Cryogenic Engineering. Adv. Cryogenic Eng. 1975, 20, 269–282. [Google Scholar]
  15. Urieli, I.; Rallis, C.J.; Berchowitz, D.M. Computer simulation of Stirling cycle machines. In Proceedings of the 12th Intersociety Energy Conversion Engineering Conference, Washington, DC, USA, 28 August–2 September 1977; pp. 1512–1521. [Google Scholar]
  16. Berchowitz, D.M.; Rallis, C. A computer and experimental simulation of stirling cycle machines. In Proceedings of the 13th Intersociety Energy Conversion Engineering Conference, San Diego, CA, USA, 20–25 August 1978; pp. 1730–1738. [Google Scholar]
  17. Gedeon, D. Sage-Object-oriented software for Stirling machine design. In Proceedings of the Intersociety Energy Conversion Engineering Conference, Monterey, CA, USA, 7–12 August 1994; p. 4106. [Google Scholar]
  18. Gedeon, D. Sage Stirling-Cycle Model-Class Reference Guide; Gedeon Associates: Athens, OH, USA, 1999; p. 16922. [Google Scholar]
  19. Gedeon, D. Sage User’s Guide; Gedeon Associates: Athens, OH, USA, 2009; p. 16922. [Google Scholar]
  20. Mahkamov, K.; Djumanov, D. Three-dimensional CFD modeling of a Stirling engine. In Proceedings of the 11th International Stirling Engine Conference, Rome, Italy, 19–21 November 2003; pp. 97–107. [Google Scholar]
  21. Ibrahim, M.B.; Zhang, Z.-G.; Tew, R.C., Jr.; Gedeon, D.; Simon, T.W. CFD modeling of free-piston Stirling engines. In Proceedings of the 36th Intersociety Energy Conversion Engineering Conference, Savannah, GA, USA, 29 July–2 August 2001. [Google Scholar]
  22. Wilson, S.; Dyson, R.; Tew, R.; Ibrahim, M. Multi-D CFD Modeling of Free-Piston Stirling Convertor at NASA GRC. In Proceedings of the 2nd International Energy Conversion Engineering Conference, Providence, RI, USA, 16–19 August 2004; p. 5673. [Google Scholar]
  23. Lee, K.-S.; Lee, S.-H.; Park, J.-H.; Choi, J.-Y.; Sim, K.-H. Design and experimental analysis of a 3 kW single-phase linear permanent magnet generator for stirling engines. IEEE Trans. Magn. 2018, 54, 1–5. [Google Scholar] [CrossRef]
  24. Schreiber, J.G.; Geng, S.M.; Lorenz, G.V. Re-1000 Free-Piston Stirling Engine Sensitivity Test Results; National Aeronautics and Space Administration Report; National Aeronautics and Space Administration: Washinton, DC, USA, 1986. [Google Scholar]
  25. Khamneh, M.E.; Askari-Paykani, M.; Shahverdi, H.; Hadavi, S.M.M.; Emami, M. Optimization of spring-back in creep age forming process of 7075 Al-Alclad alloy using D-optimal design of experiment method. Measurement 2016, 88, 278–286. [Google Scholar] [CrossRef]
  26. Kanlayasiri, K.; Boonmung, S. Effects of wire-EDM machining variables on surface roughness of newly developed DC 53 die steel: Design of experiments and regression model. J. Mater. Process. Technol. 2007, 192, 459–464. [Google Scholar] [CrossRef]
  27. Lee, C.-W.; Kim, D.-J.; Kim, S.-K.; Sim, K.-H. Design Optimization of Flexure Springs for Free-Piston Stirling Engines and Experimental Evaluations with Fatigue Testing. Energies 2021, 14, 5156. [Google Scholar] [CrossRef]
  28. Lane, N.; Berchowitz, D.; Shade, D.; Karandikar, A. Development of a high frequency Stirling engine-powered 3 kW (e) generator set. In Proceedings of the 24th Intersociety Energy Conversion Engineering Conference, Washington, DC, USA, 6–11 August 1989; pp. 2213–2218. [Google Scholar]
  29. Thorsen, J.E.; Bovin, J.; Carlsen, H. 3 kW Stirling engine for power and heat production. In Proceedings of the IECEC 96—31st Intersociety Energy Conversion Engineering Conference, Washington, DC, USA, 11–16 August 1996; Volume 1282, pp. 1289–1294. [Google Scholar]
  30. Carlsen, H. Operating experience with a 10 kW hermetic Stirling Engine. In Proceedings of the Intersociety Energy Conversion Engineering Conference, Monterey, CA, USA, 7–12 August 1994; p. 3802. [Google Scholar]
  31. Timoumi, Y.; Tlili, I.; Nasrallah, S.B. Design and performance optimization of GPU-3 Stirling engines. Energy 2008, 33, 1100–1114. [Google Scholar] [CrossRef]
  32. Formosa, F.; Despesse, G. Analytical model for Stirling cycle machine design. Energy Convers. Manag. 2010, 51, 1855–1863. [Google Scholar] [CrossRef] [Green Version]
  33. Babaelahi, M.; Sayyaadi, H. Analytical closed-form model for predicting the power and efficiency of Stirling engines based on a comprehensive numerical model and the genetic programming. Energy 2016, 98, 324–339. [Google Scholar] [CrossRef]
  34. Ahmed, F.; Hulin, H.; Khan, A.M. Numerical modeling and optimization of beta-type Stirling engine. Appl. Therm. Eng. 2019, 149, 385–400. [Google Scholar] [CrossRef]
  35. Augenblick, J.E.; McFadden, G.D.; Peterson, A.A.; Rink, K.K. Tubular Heat Exchange. U.S. Patent US206355A, 15 August 2013. [Google Scholar]
  36. Rogers, G.F.C.; Mayhew, Y.R. Engineering Thermodynamics, Work and Heat Transfer; Longman Scientific: Harlow, UK, 1992. [Google Scholar]
  37. Incropera, F.P.; DeWitt, D.P.; Bergman, T.L.; Lavine, A.S. Fundamentals of Heat and Mass Transfer; Wiley: New York, NY, USA, 1996; Volume 6. [Google Scholar]
  38. Geankoplis, C.J. Transport Processes and Separation Process Principles, 4th ed.; Person Education: London, UK, 2003. [Google Scholar]
  39. Tanaka, M.; Yamashita, I.; Chisaka, F. Flow and heat transfer characteristics of the Stirling engine regenerator in an oscillating flow. Ser. 2, Fluids engineering, heat transfer, power, combustion, thermophysical properties. JSME Int. J. 1990, 33, 283–289. [Google Scholar]
  40. Formosa, F. Coupled thermodynamic–dynamic semi-analytical model of free piston Stirling engines. Energy Convers. Manag. 2011, 52, 2098–2109. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Cross-sectional view of a 3 kW FPSE design developed by Nexergy lab at Seoultech.
Figure 1. Cross-sectional view of a 3 kW FPSE design developed by Nexergy lab at Seoultech.
Energies 15 03326 g001
Figure 2. Shape structure of the cylinder head and the HHX tube.
Figure 2. Shape structure of the cylinder head and the HHX tube.
Energies 15 03326 g002
Figure 3. Maximum number of tubes applicable according to applicable tube inner diameter and tube bending radius.
Figure 3. Maximum number of tubes applicable according to applicable tube inner diameter and tube bending radius.
Energies 15 03326 g003
Figure 4. The thermodynamic system, temperature, and pressure levels of the FPSE.
Figure 4. The thermodynamic system, temperature, and pressure levels of the FPSE.
Energies 15 03326 g004
Figure 5. Schematic of the RE-1000 engine (built by NASA).
Figure 5. Schematic of the RE-1000 engine (built by NASA).
Energies 15 03326 g005
Figure 6. Flow chart of the calculation procedure of the iQSF model.
Figure 6. Flow chart of the calculation procedure of the iQSF model.
Energies 15 03326 g006
Figure 7. Predicted heat in and indicated power for the RE-1000 engine using the finite difference method (FDM) for the regenerator.
Figure 7. Predicted heat in and indicated power for the RE-1000 engine using the finite difference method (FDM) for the regenerator.
Energies 15 03326 g007
Figure 8. Energy flow of iQSF model for the RE-1000 engine.
Figure 8. Energy flow of iQSF model for the RE-1000 engine.
Energies 15 03326 g008
Figure 9. Effects of terms for (a) the brake power and (b) efficiency.
Figure 9. Effects of terms for (a) the brake power and (b) efficiency.
Energies 15 03326 g009
Figure 10. Variation of (a) the brake power and (b) brake efficiency according to the design factors of the tube inner diameter, tube length, and number of tubes.
Figure 10. Variation of (a) the brake power and (b) brake efficiency according to the design factors of the tube inner diameter, tube length, and number of tubes.
Energies 15 03326 g010
Figure 11. Predicted (a) brake power and (b) efficiency according to the inner tube diameter and the number of tubes using the regression model.
Figure 11. Predicted (a) brake power and (b) efficiency according to the inner tube diameter and the number of tubes using the regression model.
Energies 15 03326 g011
Figure 12. Predicted (a) brake power and (b) efficiency with respect to tube length changes when the tube inner diameter is 8.5 mm and the number of tubes is 39 through the iQSF model.
Figure 12. Predicted (a) brake power and (b) efficiency with respect to tube length changes when the tube inner diameter is 8.5 mm and the number of tubes is 39 through the iQSF model.
Energies 15 03326 g012
Figure 13. Portion of the total loss for the optimal HHX design conditions in the optimal case 19 (inner tube diameter 8.5 mm, tube length 220 mm, number of tubes 39).
Figure 13. Portion of the total loss for the optimal HHX design conditions in the optimal case 19 (inner tube diameter 8.5 mm, tube length 220 mm, number of tubes 39).
Energies 15 03326 g013
Figure 14. Major loss for the design condition that generates the maximum brake power for cases 16, 17, and 18 in Figure 11 and the major loss in the worst case 2 with the lowest brake power in Table 11.
Figure 14. Major loss for the design condition that generates the maximum brake power for cases 16, 17, and 18 in Figure 11 and the major loss in the worst case 2 with the lowest brake power in Table 11.
Energies 15 03326 g014
Table 1. Design parameters of the FPSE.
Table 1. Design parameters of the FPSE.
ParameterValueUnit
Operating conditionsWorking gasHelium
Charge pressure22bar
Heater wall temp.600°C
Cooler reject temp.25°C
Operating frequency52Hz
DP and PP amplitude11.6mm
Heat exchanger geometryNumber of CHX tubes600ea.
Inner diameter of CHX tubes1.75mm
Length of the CHX tube75mm
Regenerator geometryRegenerator typeWoven mesh
Fiber diameter150μM
Porosity88%
Working space geometryExpansion space length15mm
Compression space length25mm
Driving part geometryPP and DP diameter170mm
DP rod diameter40mm
DP-liner clearance140μM
PP-liner clearance19μM
Table 2. Design parameters of HHX for several kilowatt-class Stirling engines from the literature.
Table 2. Design parameters of HHX for several kilowatt-class Stirling engines from the literature.
Model or Author#1. RE-1000 [24]#2. SHARP [24]#3. Thorsen [29]#4. GPU-3 [31]#5. QB-80 [35]#6. Carlsen [30]
Figure Energies 15 03326 i001 Energies 15 03326 i002 Energies 15 03326 i003 Energies 15 03326 i004 Energies 15 03326 i005 Energies 15 03326 i006
Power output, kW1234710
Typefree-pistonfree-pistonkinematickinematicfree-pistonkinematic
HHX
para
meter
Number of tubes346024408024
Tube inner diameter, mm2.362.233.0238
Tube length, mm183140140245.3300365
Table 3. Values of the design factors of the HHX.
Table 3. Values of the design factors of the HHX.
FactorCodesLevel Value
−101
Tube inner diameter (mm)di369
Tube length (mm)lt100200300
Number of tubesnt3080130
Table 4. HHX tube design set of central composite design (CCD).
Table 4. HHX tube design set of central composite design (CCD).
CaseTube Inner
Diameter
(mm)
Tube Length
(mm)
Tube
Number
CaseTube Inner
Diameter
(mm)
Tube Length
(mm)
Tube
Number
case 1310030case 9630080
case 2330030case 106200130
case 3320080case 11910030
case 43100130case 12930030
case 53300130case 13920080
case 6620030case 149100130
case 7610080case 159300130
case 8620080
Table 5. Governing equations of space cells for the iQSF model (here, i = c, k, r, h, e, and j = p, d, rod).
Table 5. Governing equations of space cells for the iQSF model (here, i = c, k, r, h, e, and j = p, d, rod).
Space CellGoverning Equations
Compression space H ˙ s l , d + H ˙ s l , p + H ˙ s l , r o d c p m ˙ c k T c k = c v R p ˙ c V c + c p R p c V ˙ c
Expansion space Q ˙ s h Q ˙ c o n d , d H ˙ s l , d + c p m ˙ h e T h e = c v R p ˙ e V e + c p R p e V ˙ e
Heater Q ˙ h + c p m ˙ r n , h T r n , h m ˙ h e T h e = c v R p ˙ h V h
Regenerator Q ˙ r 1 Q ˙ r e h 1 + c p m ˙ k r 1 T k r 1 m ˙ r 1 r 2 T r 1 r 2 = c v R p ˙ r 1 V r 1

Q ˙ r n Q ˙ r e h , n + c p m ˙ r n 1 , r n T r n 1 , r n m ˙ r n , h T r n , h = c v R p ˙ r n V r n
Cooler Q ˙ k + c p m ˙ c k T c k m ˙ k r 1 T k r 1 = c v R p ˙ k V k
Table 6. Definitions of losses and corresponding equations (here, i = c, k, r, h, e, and j = p, d, rod).
Table 6. Definitions of losses and corresponding equations (here, i = c, k, r, h, e, and j = p, d, rod).
LossesEquations
External heat lossesStructure conduction loss Q ˙ c o n d , c a s e + Q ˙ c o n d , c y l = k c a s e A c a s e L c a s e T w h T w k + k c y l A c y l L c y l T w h T w k
Internal heat lossesRegenerator efficiency loss Q ˙ r = 1 ε h r A w t , r T w r T r
Appendix gap loss Q ˙ · s h + Q ˙ c o n d , d = π d d 2 h g a p k g a s X d 2 T e T c L a p p + k d A d L d T e T c
leakage loss Q ˙ s l , j = m ˙ s l , j c p T i n . j T o u t . j
Power lossesSeal friction loss W ˙ s f , j = u j μ L s e a l , j π d j u j H j + μ d j h j p o u t , j p i n , j 2
Flow friction loss Q ˙ f r , i = Δ p i A f r , i g · i ρ · i
Gas spring hysteresis loss W ˙ · g s h = 1 32 ω γ 3 γ 1 T w p · k g a s ( V ^ b V b ) 2 A w t
Table 7. Correlation equations of the heat exchanger.
Table 7. Correlation equations of the heat exchanger.
ModelCorrelation
Rogers and Mayhew [36]Re ≤ 2000,
2000 < Re ≤ 4000,
Re > 4000,
f f = 16 Re 1  
f f = 7.343 × 10 4 Re 0.314
f f = 0.0791 Re 0.25
Reynolds analogy [37,38]
(Re ≤ 4000)
f f = Re 2 = Nu ,   h = 1 2 f f Pr 1 g c p
Chilton-Colburn J-factor analogy [37,38]
(Re > 4000)
f f 2 = St P r 2 / 3 ,   h = 1 2 f f Pr 2 / 3 g c p
Table 8. Correlation equations of the regenerator [39].
Table 8. Correlation equations of the regenerator [39].
Flow TypeOscillating Flow
Regenerator typeWoven wire (Stainless steel)
Wire diameter0.05–0.23 mm
Porosity0.645–0.754
Friction factor f f = 175 Re + 1.6
Heat transfer coefficient Nu = 0.33 Re 0.67
Table 9. Design parameters and operating conditions of the RE-1000 engine for the #1011 test.
Table 9. Design parameters and operating conditions of the RE-1000 engine for the #1011 test.
Design ParameterValueUnitDesign ParameterValueUnit
Heater
(tube type)
Number of tubes34 DPDiameter56.4mm
Tube inner diameter2.36mmAppendix gap clearance0.381mm
Heater volume27.33cm3Seal clearance0.033mm
Cooler
(fin type)
Fin number135 DP rodDiameter16.7mm
Cooler volume20.43cm3Seal clearance0.03mm
Regenerator
(woven mesh type)
Porosity75.9%PPDiameter57.2mm
Wire diameter0.0889mmSeal clearance0.033mm
Void volume56.37cm3Operating conditionHeater temperature600°C
SpaceExpansion space volume27.74cm3Cooler temperature25°C
Compression space volume54.8cm3Charge pressure70.6bar
DP bounce volume30.79cm3Frequency30Hz
PP bounce volume2615cm3DP stroke24.5mm
DisplacerDiameter56.4mmPP stroke28mm
Appendix gap clearance0.381mmDP-PP phase angle57.5°
Seal clearance0.033mm
Table 10. Comparison of various model results.
Table 10. Comparison of various model results.
Formosa Model [40]QSF Model [8]iQSF ModelExperiment (#1011) [24]
Heat in (W)3777 (−6.46%)4949.0 (22.6%)4614.7 (14.3%)4038
Heat rejection (W)−2107 (−30.5%)−2742.4 (−9.5%)−3329.2 (9.8%)−3032
Indicated power (W)1100 (15.6%)1717.4 (66.9%)1286.8 (24.9%)1030
Indicated efficiency (%)29.2 (23.3%)34.5 (35.3%)27.9 (9.4%)25.5
Table 11. Results of the CCD case.
Table 11. Results of the CCD case.
CaseBrake Power (W)Brake Efficiency (%)CaseBrake Power (W)Brake Efficiency (%)
case 1463922.71case 9755935.92
case 213507.53case 10742435.4
case 3715733.4case 11667836.04
case 4796738.7case 12744336.18
case 5754335.52case 13699934.87
case 6769137.44case 14652733.62
case 7744638.24case 15626928.45
case 8765736.86
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, D.-J.; Park, Y.; Kim, T.Y.; Sim, K. Design Optimization of Tubular Heat Exchangers for a Free-Piston Stirling Engine Based on Improved Quasi-Steady Flow Thermodynamic Model Predictions. Energies 2022, 15, 3326. https://doi.org/10.3390/en15093326

AMA Style

Kim D-J, Park Y, Kim TY, Sim K. Design Optimization of Tubular Heat Exchangers for a Free-Piston Stirling Engine Based on Improved Quasi-Steady Flow Thermodynamic Model Predictions. Energies. 2022; 15(9):3326. https://doi.org/10.3390/en15093326

Chicago/Turabian Style

Kim, Dong-Jun, Yeongchae Park, Tae Young Kim, and Kyuho Sim. 2022. "Design Optimization of Tubular Heat Exchangers for a Free-Piston Stirling Engine Based on Improved Quasi-Steady Flow Thermodynamic Model Predictions" Energies 15, no. 9: 3326. https://doi.org/10.3390/en15093326

APA Style

Kim, D. -J., Park, Y., Kim, T. Y., & Sim, K. (2022). Design Optimization of Tubular Heat Exchangers for a Free-Piston Stirling Engine Based on Improved Quasi-Steady Flow Thermodynamic Model Predictions. Energies, 15(9), 3326. https://doi.org/10.3390/en15093326

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop