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

: 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 pre ‐ dicted 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 de ‐ termined 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.


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 semipermanent 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.

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.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.4mm and 183 mm, respectively.The SHARP (Stirling helium advanced reliable prototype) engine [28], 2 kWclass 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.

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  of the bent tube.The tube length  is defined by Equation (1): where  0.0084 1.4684 1.4434.Here,  is the straight length of the inner tube and  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  .For the tube bend radius  , 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  according to Equation ( 1) is approximately 25 mm when  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: where  is the outer diameter of the tube and  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  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.

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: 1.The working fluid is considered an ideal gas.2. Changes in the kinetic and potential energies are negligible.3. The incompressible effect of the working fluid is negligible.4. The mass of the working fluid inside the engine is constant.5.All flows of the working fluid inside the engine are assumed to be periodic steadystate flows.6.The motion of DP and PP follows a sinusoidal waveform with a constant velocity, amplitude, and phase.7. 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 Tables 5 and 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.Reynolds analogy [37,38] (Re ≤ 4000)  Nu, ℎ  Pr  Chilton-Colburn J-factor analogy [37,38] (Re > 4000) 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 Tables 5 and 6, the expansion work, compression work, and indicated power generated through the Stirling cycle are defined as follows: Here,  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: 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:

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 onedimensional 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.

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 ( ) and the number of tubes ( ); the p-value is less than 0.05, indicating a statistically significant effect on the brake power and efficiency.This is because  •  term is the factor that increases the dead volume the most.As a single term, performance is greatly affected in the order of  ,  , and  .In particular,  has a relatively low effect on performance compared to  and  .The rest of the terms except for the  •  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: R-square represents the fitness of the regression model for the interpreted data.Rsquare for the indicated power and the efficiency calculated a high fit above 0.85.

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  •  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  25 mm.In all graphs in Figure 11, when  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  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  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.

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 welldesigned 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.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.

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.

Figure 1 .
Figure 1.Cross-sectional view of a 3 kW FPSE design developed by Nexergy lab at Seoultech.

Figure 2 .
Figure 2. Shape structure of the cylinder head and the HHX tube.

Figure 3 .
Figure 3. Maximum number of tubes applicable according to applicable tube inner diameter and tube bending radius.

Figure 4 .
Figure 4.The thermodynamic system, temperature, and pressure levels of the FPSE.

Figure 6 .
Figure 6.Flow chart of the calculation procedure of the iQSF model.

Figure 7 .
Figure 7. Predicted heat in and indicated power for the RE-1000 engine using the finite difference method (FDM) for the regenerator.

Figure 8 .
Figure 8. Energy flow of iQSF model for the RE-1000 engine.

Figure 9 .
Figure 9. Effects of terms for (a) the brake power and (b) efficiency

Figure 10 .
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 11 .
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 12 .
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 13 .
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 14 .
Figure 14.Major loss for the design condition that generates the maximum brake power for cases 16, 17, and 18 in Figure11and the major loss in the worst case 2 with the lowest brake power in Table11.

Table 1 .
Design parameters of the FPSE.

Table 3 .
Values of the design factors of the HHX.

Table 4 .
HHX tube design set of central composite design (CCD).

Table 5 .
Governing equations of space cells for the iQSF model (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).

Table 7 .
Correlation equations of the heat exchanger.

Table 9 .
Design parameters and operating conditions of the RE-1000 engine for the #1011 test.

Table 10 .
Comparison of various model results.

Table 11 .
Results of the CCD case.