1. Introduction
The decreasing crude oil sources, the increasing concern about greenhouse gas pollution, and the strict emission standards force researchers into neverending effort to design less polluting and more fuel efficient vehicles. Freepiston linear alternator (FPLA) which is a freepiston engine coupled to a linear alternator, converting piston’s kinetic energy into electricity directly, has attracted considerable research interests recently by a number of research groups worldwide due to its potential advantages listed bellow:
High efficiency. 1) Direct conversion of piston motion into electric energy with shorter energy transferring path. 2) Reduction of friction between piston ring and cylinder wall due to the elimination of the dynamic piston side forces.
High power capacity. With a freepiston engine, heavy mechanical components as crankshaft, flywheel and camshaft are eliminated and the generator is integrated into the engine directly. Therefore, the powertrain will have fewer components, lower weight and also occupies less space.
Controllable compression ratio. The stroke and thus the compression ratio can be dynamically controlled through the electrical load.
Multifuel feasibility. The variable compression ratio feature provides the opportunity for the engine to be operated with wide range of fuels [1].
Few moving parts. There is only one major moving part: the piston, connecting rod and translator assembly. As a consequence, the engine requires less lubrication and has less friction wear.
Good transient response. The freepiston engine reaches its operating point almost immediately since it does not have any energy accumulators, such as flywheel [2].
In the recent two decades, some literatures have already been published about the design and numerical simulation of FPLA. Martin Goertz and Lixin Peng evaluated several feasible hybrid powertrain concepts and the results indicated that freepiston linear generator was one of the most promising candidates of the future powertrain configuration [3]. West Virginia University had already demonstrated stable operation of a spark ignited FPLA prototype with bore of 36.5mm, maximum possible stroke of 50mm and 316W output power was produced at 79V working at full load [4]. Ehab Shoukry did some research about parametric study of a twostroke direct injection linear engine using zero dimensional single zone models [5].The European Union has been doing a subject of Freepiston Energy Converter (FPEC) since 2002 aimed to develop an efficient new technology suitable for vehicle propulsion, auxiliary power unit and distributed power generation. The design and optimization of linear alternator and crankless engine, control strategy, numerical simulation and power management have been carried out systemically [6]. Dr. Peter Van Blarigan at Sandia National Laboratory presented the design of a dual piston freepiston engine generator with 40kW electric power output and the engine employed homogeneous charge compression ignition (HCCI) operating on a variety of hydrogencontaining fuels [7, 8]. Jakob Fredriksson and Ingemar Denbratt at Chalmers University of Technology investigated a twostroke freepiston engine with different fuels using BOOST and SENKIN [9]. Czech Technical University had successfully developed a direct injection FPLA prototype recently, and steady operation had been realized based on precise motion control. When the prototype was running with frequency of 27Hz and compression of 9, the average power output was approximately 350W, but the efficiency had not been reported [10]. Mikalsen and Roskilly proposed a design of a single cylinder freepiston engine generator with gasfilled bounce chamber, simulated its working process and discussed the effects of parameters in a wide operating range to the engine’s performance [11]. A novel approach of modeling the freepiston engines through the introduction of solutiondependent mesh motion in an engine CFD toolkit OpenFOAM was also presented [12].
Dimensionless analysis using Buckingham’s theory was performed base on the parameters of a spark ignited FPLA prototype (as can be seen in Fig.1) in the following contents. The objectives of this work are:
To decrease the number of variables and still use the equations involved in the modeling of FPLA.
To generalize the study of FPLA to a wide range of geometrical designs in order to find out the best geometric dimensions with most favorable performance.
To establish a scale for the similitude concept so the experiments done on the existing prototype could be used to predict performances for FPLA with various sizes and dimensions.
2. System modeling
2.1. Dynamic modeling
The Modeling of FPLA and numerical simulation of its performance have already been done by a lot of researchers. Usually the model consists of a dynamic calculation of the piston motion and a thermodynamic calculation of the engine’s thermodynamic events. The piston motion is governed by the interaction of forces that act on the piston simultaneously, which are incylinder pressure force from each cylinder, electromagnetic force, friction force and inertia force, as is shown in Fig.2.
The Newton’s second law is applied:
Since the FPLA doesn’t have a crankshaft, it is free of the dynamic friction caused by the crankshaft mechanism and the friction force is caused mainly by the interactions between cylinder and piston rings, cylinder and piston skirt, which is small compared with the load of the linear alternator. Therefore, the friction force is considered to be constant during the whole cycle in the calculation [13].
2.2. Modeling of the linear alternator
The linear alternator consists of two main components, a stator and a translator. The permanent magnets are mounted on the stator and the translator is the moving portion of the machine which is made up of coils. A schematic of a threephase, “U” shaped linear alternator with permanent magnet (PM) excitation is shown in Fig.3.
The FPLA operates on the same basic physical principles as conventional rotary alternators. The principle that governs the voltage generating operation of the alternator is Faraday’s law expressed as [14]:
The permanent magnets create a magneto motive force (MMF) in the air gap between the stator and the winding coils as shown in Fig.4, and it can be described by the following mathematic equation:
Where
The mean value of the MMF can be obtained from the singleorder truncated Fourier series [15]:
Where,
So the flux density in the air gap due to PM is:
Where
Both experimental measurements and numerical calculation by finite element method showed that the flux in the air gap of the PM exited linear alternator in Fig.4 could be assumed to be sinusoidal supporting the above result [16].
Therefore, the flux contained in the differential element
Then the total flux contained in the coil of one phase at random position
Thus, the induced electromotive force produced in the coil of one phase is:
The induced current from the load circuit can be derived in the following equations:
The magnetic force has the opposite direction to the direction of the translator’s movement. According to Ampere’s law, it is described in the following equation:
When it comes to threephase linear alternator, the third phase is derived from another two according to the following equation [17]:
So the total electromagnetic force produced by a threephase linear alternator is:
Where
In order to develop a first approach analysis, it’s assumed that the alternator circuit works in resonance condition, which means that the current and voltage have the same phase, so that it is possible to consider the circuit as only a resistance [18]. Thus, the electromagnetic force is proportional to the speed of the translator:
2.3. Thermodynamic modeling
The thermodynamic model is derived based on the first law of thermodynamics and ideal gas law. It consists of the calculation of the process of scavenging, compression, combustion, expansion and exhaust. The zerodimensional, single zone model is used to describe the thermodynamic process.
Appling the first law of thermodynamics and ideal gas law on the cylinder as an open thermodynamic system, shown in Fig.5, and assuming that the specific heat
For the case of compression and expansion process neglecting the crevice flow and the leakage, the first law of thermodynamics applied to the cylinder content becomes:
Considering the cylinder content is ideal gas, and then at every instant the ideal gas law is satisfied:
Substitution and mathematical manipulation yield the following equation which is used to calculate the incylinder pressure at each time step.
In the combustion model, since the engine is crankless, a timebased Wiebe functions (as opposed to a conventional crank anglebased approach) is used to express the mass fraction burned in combustion process [4]:
and the heat rate released during combustion is:
The incylinder heat transfer effect is modeled according to Hohenberg [19]:
Since most heat transfer models, like the ones proposed by Woschni and Hohenberg, are made for Diesel engines. This means that they take radiative heat transfer effect into account, which is hardly present in premixed combustion. Hence, in the numerical simulations a factor of 0.5 is introduced to reduce the heat transfer coefficient [9].
So the total energy that is used to increase the incylinder pressure in equation (18) can be expressed in the following equation:
Exhaust blown down is modeled to be a polytrophic expansion process while the exhaust port is opening and the scavenging ports are still covered by the piston [11].
For twostroke spark ignition engines with under piston or crankcase scavenging, the scavenging efficiency is about 0.7~0.9 [20],, a scavenging efficiency of 0.8 is introduced to evaluate the effects of incomplete scavenging effect. The moment the scavenging ports are open, the pressure and temperature are assumed to be the same with the scavenging conditions and the incoming gases mix entirely with the burned gases.
3. Dimensionless analysis
3.1. Dimensionless modeling of FPLA
The dimension of a physical quantity is associated with combinations of mass, length, time, temperature and heat quantity, represented by symbols [M], [L], [T], [θ] and [H] respectively, each raised to rational powers. Since heat quantity is also a kind of energy, its dimension is replaced with [ML^{2}T^{2}].
The basic variables that involved were defined based on the dynamic and thermodynamic modeling of FPLA. The variables and their dimensions are listed in Tab.1.
Variables  Symbol  Dimensions 
Bore  D  [L] 
Piston area  A  [L^{2}] 
Volume of the cylinder  V  [L^{3}] 
Effective stroke length  L_{eff}  [L] 
Total stroke length  L_{tot}  [L] 
Translator ignition position  x_{ign}  [L] 
Mass of the translator  m  [M] 
Load coefficient  M  [MT^{1}] 
Friction force  F_{f}  [MLT^{2}] 
Energy  Q  [ML^{2}T^{2}] 
Scavenge pressure  p_{0}  [ML^{1}T^{2}] 
Scavenge temperature  T_{0}  [θ] 
Pressure  p  [ML^{1}T^{2}] 
Gas constant  R  [θL^{1}T^{2}] 
Temperature  T  [θ] 
Wall temperature  T_{w}  [θ] 
Combustion duration  t_{c}  [T] 
Compression ratio  ε  [1] 
indicated efficiency  η_{i}  [1] 
Effective efficiency  η_{e}  [1] 
Time  t  [T] 
Frequency  f  [T^{1}] 
The effective stroke length is the distance between the upper edge of the exhaust port and the cylinder and the total stroke length is distance that the translator can travel from cylinder head to cylinder head [21].
Since the FPLA is controlled based on the displacement feedback of displacement sensor, the ignition timing is defined by the position of the translator where the controller gives out the ignition signals.
According to Buckingham’s theory [22], four fundamental reference variables which are independent from each other were chosen, and the other variables are described with these reference variables in their index form and every redefined variable is dimensionless. The four reference variables chosen were bore, mass of the translator, scavenging pressure and scavenging temperature.
Dividing the basic variables by the reference variables results in the dimensionless variables:
Dimensionless bore
Dimensionless piston area
Dimensionless cylinder volume
Dimensionless total stroke length
Dimensionless mass of the translator
Dimensionless scavenge air pressure
Dimensionless scavenging temperature
Dimensionless temperature
Dimensionless wall temperature
Dimensionless load coefficient
Since the dimensionless translator ignition position is changing with different stroke length and fixed translator ignition position, in the dimensionless calculation, its actual dimensionless value taken is defined by the compression ratio the engine has already achieved when the spark plug ignites, as is described in Fig.6.
Then the dimensionless translator ignition position can be deduced:
Substituting all the dimensionless variables into equations (1) and (18), the following equations can be derived:
Where
The dimensionless gas constant was deduced based on the ideal gas law:
Since the incylinder temperature is strongly transient, the dimensionless incylinder temperature is acquired using the dimensionless ideal gas law:
So the dimensionless variables defined in equations (28) (30) and (31) are:
Dimensionless time
Dimensionless combustion duration
Dimensionless frequency
Dimensionless friction force
Dimensionless pressure
Dimensionless energy
Dimensionless gas constant
The combustion model and heat transfer equations can also be transferred to their dimensionless form using the dimensionless parameters we have already obtained:
All the equations were solved with a numerical simulating program in Matlab. The program starts calculation using a group of mathematic equations describing the dimensionless dynamic and dimensionless thermodynamic processes of FPLA, including the heat transfer rate, the incylinder pressure, the incylinder temperature, the fraction of fuel that burnt, the load and the work done, etc. At each time step, the program calls the dynamic subroutine and updates the displacement, the velocity, the averaged velocity and the acceleration. After the engine stabilizes, the program calculates the engine frequency, the compression ratio, the indicated power, the frictional power, the effective power output and the effective efficiency.
The dimensionless form of the output parameters are:
Dimensionless compression ratio
Dimensionless frictional work
Dimensionless averaged speed
Dimensionless indicated work
Dimensionless effective work
Dimensionless indicated efficiency
Dimensionless effective efficiency
Dimensionless effective power output
Dimensionless frictional power
The dimensionless compression ratio and efficiency have just the same value with their dimensional forms, which are used to validate the correctness of the dimensionless process.
3.2. Operating ranges
The variables of the dimensionless analysis were chosen to cover the normal operating ranges of twostroke freepiston engine. The base case represented the parameters of the experimental spark ignited, twostroke FPLA prototype that was built in Beijing Institute of Technology. The basic parameters are listed in Tab.2.
The effective stroke length to bore ratio, dimensionless load coefficient, dimensionless translator ignition position, dimensionless combustion duration and dimensionless input energy (opening proportion of the throttle) were the variable input factors of the parametric study. The dimensionless operating matrixes of each parameter are listed in Tab.3.
Basic parameters  Value  Dimensionless value 
D/mm  34  1 
Leff/mm  20  0.5882 
Ltot/mm  36  1.0588 
Rl/Ω  2.5  × 
Rs/Ω  2  × 
M/N·(m·s1)1  55.3  0.7746 
m/kg  1.74  1 
p0/bar  1  1 
T0/K  293  1 
Tw/K  453  1.5461 
tc/ms  5  0.2210 
εign  4  × 
α  30%  × 
Parameters  Value  
case1  case2  case3  case4  case5  case6  case7  
Leff*  0.5  0.6  0.7  0.8  0.9  1.0  1.1 
RL/Ω  2  2.5  3  4  ×  ×  × 
M/N·(m·s1)1  62.2320  55.3174  49.7856  41.4880  ×  ×  × 
M*  0.8714  0.7746  0.6971  0.5809  ×  ×  × 
εign  3  4  5  6  ×  ×  × 
tc/ms  3  4  5  6  ×  ×  × 
tc*  0.1326  0.1768  0.2210  0.2652  ×  ×  × 
α  25%  30%  35%  ×  ×  ×  × 
The dimensionless effective stroke length was defined as the first independent variable while the dimensionless load coefficient, dimensionless translator ignition position, dimensionless combustion duration and dimensionless input energy were taken as second independent variables. When analyzing the effects of each second variable, the other second variables would be set to equal to the base case which was decided by the FPLA prototype.
4. CFD calculation of combustion process for validation
4.1. Combustion modeling of FPLA
Since the combustion model used in the numerical program is zero dimensional Wiebe function and some parameters like combustion duration are of great uncertainty, the accuracy of the numerical calculated results is suspectable. Nowadays, multidimensional CFD computational tools have become an integral part of the engine design process due mainly to advances in computing capabilities and improvements in the modeling techniques utilized. In this study, in order to validate the results of dimensionless analysis, a multidimensional commercial CFD software AVL_Fire was used to evaluate the effects of translator ignition position with different effective stroke length to bore ratio.
As the FPLA doesn’t have a crankshaft mechanism, the dynamics of the piston is totally different from conventional engine. The dynamics were defined based on the results of a zero dimensional FPLA modeling mentioned in the former paragraphs. The dynamics and thermodynamics equations of FPLA in section 2 were solved using a numerical simulating program in Matlab, and some of the parameters were defined according to the experimental data measured. Then the dynamics of the FPLA were incorporated into AVL_Fire to define the movement of the piston. The piston motion profile was described with two arrays of numbers, one of which represented the ECA (Here ECA is equivalent crank angle which is used to note the port timings. However, it is only a time notation since the freepiston engine does not have a crankshaft to define the piston’s motion and
Only compression, combustion and expansion processes of the freepiston engine were calculated in order to minimize the number of computational cells (intake port, scavenging ports, exhaust port and scavenging case were not included in the combustion process). The computational model of the cylinder is shown in Fig.7, and the basic geometry is defined based on the FPLA prototype. Due to the symmetry of the cylinder ports layout, it is only necessary to model half of the geometry in order to minimize the computational cost.
4.2. Boundary conditions and combustion model
The boundary conditions were chosen to reflect the physical conditions which exist in the validation model and the prototype engine. Constant wall temperatures were also used. The standard
The initial conditions in the cylinder, such as pressure, temperature, EGR ratio and kinetic energy were defined based on the results of multidimensional scavenging calculation using CFD tools. In order to investigate the effects of ignition timings, the initial conditions were the same with base case. The input energy in the cylinder had just the same value as it is in the numerical simulation program.
Probability Density Function (PDF) model which takes into account the simultaneous effects of both finite rate chemistry and turbulence was chosen to describe the combustion process. The benefits of the PDF approach lie in the fact that it provides a complete statistical description of the scalar quantities under consideration. Thus, it allows first (mean values), second (variance), and even higher (skewness) order moments to be easily extracted, and that the term expressing the rate of chemical reaction appears in closed form.
4.3. CFD calculation cases
Based on the basic geometry of the FPLA, two kinds of effective stroke length to bore ratio and four ignition compression ratios were chosen in the CFD calculation. The other parameters are the same with the base case mentioned before.
As the piston dynamics is changing with different operating conditions, the piston motion profiles have to be defined first in the numerical simulation program, and then the required data in the CFD calculation can be derived, which are listed in Tab.4. The other parameters are based on the FPLA prototype. The revolution of the engine doesn’t has real meanings as the free piston engine does not have a crankshaft, and it is bring forward just to complete the combustion process required by the CFD software.
Practically, the maximum compression ratio is confined by the geometry of the chamber since the roof of cylinder and piston is not flat, as is shown in Fig.7. Thus, we have to make sure that the compression ratios of the typical effective stroke length chosen do not exceed their maximum value.
Leff/mm  D/mm  Leff /D  εign  tign/ECA  ε  f/Hz  EGR  Qin/J  Revolution/rpm 
23  34  0.6765  3  319.5  21.0  46.7  20%  17.84  2803.7 
23  34  0.6765  4  328.5  22.2  45.0  20%  17.84  2697.8 
23  34  0.6765  5  333.6  21.6  43.4  20%  17.84  2606.4 
23  34  0.6765  6  337.4  20.7  42.2  20%  17.84  2531.6 
35  34  1.0294  3  320.7  11.4  34.1  20%  27.15  2046.4 
35  34  1.0294  4  328.6  13.8  33.7  20%  27.15  2025.0 
35  34  1.0294  5  333.7  14.6  33.0  20%  27.15  1980.2 
35  34  1.0294  6  335.4  14.5  32.1  20%  27.15  1928.0 
The piston motion profiles with different operating conditions listed in Tab.4 are shown in Fig.8.
5. Results and discussion
5.1. Zero dimensional models validation
In order to make sure most part of the numerical simulation program is correct, the numerical simulated incylinder pressure is compared with experimental derived pressure in Fig.9. The thermodynamic model used in the numerical analysis is a single zone model. In single zone models the cylinder mixture composition, pressure and temperature are considered to be uniform and the energy release rate is modeled using experiential model. Matching the experimentally and numerically derived pressure profiles was a complicated task due to the simplification and unknown variables, such as the actual heat addition, the combustion duration and the actual load when the engine was running. Seen in Fig.9, the numerical simulation model used to simulate the operation of FPLA proved to be in agreement with the experimental data with a certain combination of the variables.
5.2. Effects of dimensionless parameters
5.2.1. Effects of dimensionless effective stroke length
The stroke to bore ratio is one of the core design variables in internal combustion engines, relating combustion chamber surface area to its volume and piston area to stroke length [11]. In order to predict the performance of FPLA with various sizes and dimensions, seven cases of dimensionless effective stroke length were chosen as the basic variable to investigate the performance of the FPLA in wide operating ranges.
Fig.10 illustrates the change of dimensionless velocity versus dimensionless displacement for different effective stroke length to bore ratios while the other parameters remain the same with the base case. It can be seen that as the dimensionless effective stroke length increases, the dimensionless velocity increases. However, the profile of all the curves is similar to each other, which means this kind of freepiston engine has its own specific characteristics.
As is shown in Figs.11~13, the dimensionless compression ratio, dimensionless frequency and dimensionless incylinder peak pressure keep decreasing as the effective stroke length to bore ratio increases. These are because as the dimensionless effective stroke length increases, the translator has to travel longer strokes and more energy is wasted overcoming friction. And as the effective stroke length to bore ratio increases, the cylinder contains more charge and the charge would contain more energy during compression stoke; therefore the dimensionless compression ratio would decrease.
The highest dimensionless effective efficiency is achieved while L_{ eff } ^{*} is 0.8 under the basic working conditions, and the peak point is mainly affected by the other four dimensionless parameters, which will be discussed in the following sections.
As the effective stroke length to bore ratio increases, the dimensionless energy input to the engine every cycle increases as a result of increasing the volume of the cylinder. Since the dimensionless effective power output is also strongly determined by the dimensionless frequency and dimensionless effective efficiency, the highest effective power output is achieved while L_{ eff } ^{*} is 0.9 under the basic working conditions.
5.2.2. Effects of dimensionless load coefficient
Increasing the dimensionless load coefficient means the load demand of the linear alternator is increasing and the electromagnetic force produced by the linear alternator is increasing. Four different dimensionless load coefficients (M^{*}1>M^{*}2>M^{*}3>M^{*}4) were chosen to investigate the effects of changing the load of the linear alternator. The load coefficient was varied by changing the value of the load resistance. According to the results calculated, the dimensionless load coefficient has large impact on different parameters studied and can affect the operating condition of FPLA.
According to Figs.14~15, as the dimensionless load coefficient increases, the dimensionless compression ratio and dimensionless frequency decrease since bigger electromagnetic force is acting on the translator. The highest dimensionless effective efficiency is changing with different dimensionless load coefficient and effective stroke length to bore ratio. As is shown in Fig.14, when the effective stroke length to bore ratio is less than 0.67, smaller dimensionless load coefficient would lead to a higher dimensionless effective efficiency and when the effective stroke length to bore ratio is more than 1.0, the larger the load coefficient the higher the dimensionless effective efficiency. The reason behind these is believed to be caused by the percentage of heat released before top dead center (TDC), which is strongly determined by the frequency of the translator.
As is shown in Fig.15, smallest dimensionless load coefficient lead to the highest dimensionless power output although the dimensionless effective efficiency is the lowest since the dimensionless frequency with smaller load coefficient is higher. It is more obvious when the effective stroke length to bore ratio is more than 1.0 since smaller load coefficient lead to higher dimensionless effective efficiency and higher dimensionless frequency. Therefore, we can conclude that the main factor that controls the power output of FPLA is its frequency.
5.2.3. Effects of dimensionless translator ignition position
Ignition timing is one of the major parameters that control the engine's operating conditions, such as frequency and compression ratio. Since the dimensionless ignition timing is changing with different dimensionless stroke length, the ignition timing is defined by the compression ratio the engine has already achieved when the spark plug ignites in the calculation, and it means that the lower the ignition compression ratio is the bigger the ignition advance is.
According to some literatures [3][5], it’s held that an earlier combustion in diesel freepiston engines would lead to more waste of energy to reverse the translator, thus the efficiency and frequency would drop. However, according to the results of spark ignited FPLA obtained in this paper, with different effective stroke length to bore ratio the best ignition advance differs with each other, since an early ignition is associated with negative work in the compression stroke and a late ignition is associated with low peak incylinder pressure, as is shown in Fig.16.
As is described in Figs.17~18, with smaller effective stroke length to bore ratio (closer to 0.5), a bigger ignition advance would lead to higher dimensionless compression ratio, higher dimensionless effective efficiency, higher dimensionless frequency and higher dimensionless effective power output. The reason is that with small dimensionless effective
stroke length, the dimensionless frequency of FPLA is high and most of the energy is released after TDC. Thus, the incylinder peak pressure is higher with a bigger ignition advance, which will help improve the performance of the engine. With a high effective stroke length to bore ratio (closer to 1.1), the frequency of the engine decreases a lot since the translator has to travel a longer stroke and a bigger proportion of energy will be released before TDC, which is associated with negative work in the compression stroke. According to the results derived, when the dimensionless effective stroke length is longer than 1.0, the ignition compression ratio of 5 would leads to the best engine performance.
The dimensionless effective power output is determined by dimensionless effective efficiency and dimensionless frequency, as has been discussed before. As is shown in Fig.18, the biggest dimensionless power output is achieved when effective stroke length to bore ratio is 0.9 and ignition compression ratio is 4. Since the dimensionless frequency has little deviation with different ignition compression ratios, the dimensionless effective power output has similar trends with the dimensionless effective efficiency.
In order to analysis the effects of different ignition timings, the combustion duration was assumed to be invariant. However, the combustion duration is strongly depend on the working conditions of the engine, thus CFD tools were taken to analysis the effects of different ignition timings to verify the dimensionless results later.
5.2.4. Effect of dimensionless combustion duration
The modeling of the heat release in freepiston engine is one of the factors with the highest degree of uncertainty in the simulation model [11]. The piston motion of freepiston engines differs significantly from that of conventional engines and very little research exists on how this influences the combustion process. In the dimensionless calculation, the heat release rate is defined by the combustion duration and shorter combustion duration will lead to a faster heat release rate. Based on the base case, four cases of combustion duration were chosen to instigate its effects to the engine’s performances.
Seen in Fig.19, a shorter combustion duration which means a faster heat release rate would lead to a higher compression ratio and higher effective efficiency when the dimensionless effective stroke length is less than 0.68 and 0.75. However, as the dimensionless effective stroke length increases, the dimensionless frequency will decrease and more energy will be released before TDC. For shorter combustion duration a lot more percentage of energy is released before TDC, which is associated with more negative work in the compression stroke. Thus, shorter combustion duration would lead to a lower dimensionless compression ratio and lower dimensionless effective efficiency with a longer dimensionless effective stroke length and fixed ignition compression ratio.
As is shown in Fig.20, shorter combustion duration leads to a higher frequency with smaller dimensionless effective stroke length and as dimensionless effective stroke length grows, shorter combustion duration leads to faster decreasing of dimensionless frequency as more energy is released before TDC to stop the translator. The dimensionless effective power output is determined by the dimensionless frequency and dimensionless effective efficiency and it has a similar trend with dimensionless efficiency.
Therefore, with a longer effective stroke length to bore ratio it is recommended to postpone the ignition timing to achieve a good performance of the freepiston engine.
5.2.5. Effects of dimensionless input energy
The freepiston engine investigated in this paper is a sparkignited engine and the input energy is varied by changing the opening proportion of the throttle. For FPLA, a much narrow range of operating speeds is expected to be utilized, which is due to the electrical generating scheme employed by the device [23]. Therefore, the opening proportion of the throttle is confined to low speed range. According to the load of FPLA, efficient generation will be achieved by operation at a fixed oscillating rate.
The effects of different dimensionless input energy while other parameters remain the same with the base case are shown in Figs.21~22. As expected, with more input energy, the dimensionless frequency, dimensionless compression ratio and dimensionless effective power output of the engine are increasing since more energy is released in the combustion process.
The amount of energy input to the engine is strictly determined by the load of FPLA. If we keep increasing the amount of input energy, the current load coefficient is not suitable for the current load coefficient and the speed of the translator will keep increasing since extra energy cannot be extracted, and at last the piston will crush with the cylinder head, which is strictly forbidden. However, if we decrease the amount of input energy, the translator will stop since the amount of energy is not enough to sustain the stable operation of the engine. Therefore, the operation range of the engine is confined by the load of the linear alternator, and the amount of the input energy has to be adjusted with the load coefficient to obtain a higher efficiency or higher power output.
5.3. CFD calculated results
In order to verify the results of dimensionless translator ignition position of spark ignited freepiston engines, multidimensional CFD tools were used to calculate the combustion process of the FPLA with four different ignition timings and two kinds of effective stroke length to bore ratio..
Nomenclature
a  combustion constant  R  air gas constant 
A  top area of the piston  RL  load resistance 
Acyl  heat transfer area  Rs  internal resistance of coils 
b  combustion form factor  t  time 
B  magnetic induction intensity  t0  time combustion begins 
cV  constant volume specific heat  tc  combustion duration 
D  cylinder diameter  tign  ignition timing 
f  frequency  T  temperature 
Fe  electromagnetic force  T0  scavenge temperature 
Ff  friction force  Tw  wall temperature 
g  air gap length  U  internal energy 
h  heat transfer coefficient 

mean piston speed 
hm  thickness of the permanent magnet  V  displaced volume of the cylinder 
H  length of the coils cutting magnetic lines  Veff  effectively compressed volume of the cylinder 
Hc  magnetic field strength  Vign  volume of the cylinder when ignite 
He  enthalpy output  W  work done 
Hi  enthalpy input  We  effective work 
iL  current in the load circuit  Wf  frictional work 
L  induction  Wi  indicated work 
Ltot  total stroke length  x  displacement of the translator 
Leff  effective stroke length  xign  translator ignition position 
m  translator mass  xs  half of maximum stroke length 
min  mass of the charge  α  opening proportion of throttle 
M  load coefficient  γ  specific heat ratio 
MF  mean magneto motive force  ε  compression ratio 
n  polytrophic exponent  εign  ignition compression ratio 
Ncoil  number of turns in the coil  εind  induced voltage 
p  incylinder absolute pressure  Φ  flux passing through the coil 
p0  scavenge pressure  λ  total flux pass through the coil 
pL  pressures in the left cylinder  μ0  vacuum permeability 
pR  pressures in the right cylinder  τ  pole pitch 
Pe  effective power output  τp  width of PM 
Pf  frictional power  ηe  effective efficiency 
Q  energy  ηi  indicated efficiency 
Qc  heat released in combustion 

velocity 
Qht  heat transfer 

acceleration 
Qin  total input energy 
The incylinder pressure curves with different ignition compression ratio while L_{ eff } ^{*}=0.6765 are shown in Fig.23. It is clear that smaller ignition compression ratio or bigger ignition advance leads to higher peak pressure which is in agreement with the dimensionless results.
The incylinder pressure curves with different ignition compression ratio while L_{
eff
}
^{*}=1.0294 are shown in Fig.24. The sequence of the peak pressure achieved with different ignition compression ratio is
The combustion duration calculated via CFD is about 4.4~5.6ms with different ignition timings and effective stroke length, which has some deviation with the value in numerical simulating program which is defined based on the heat release rate of FPLA prototype. The deviations can be eliminated by using an iterative procedure between the numerical simulating program and CFD calculation when calculating a specific freepiston engine.
6. Conclusion
A detailed dimensionless modeling and dimensionless parametric study of spark ignited FPLA was presented to build up a guideline for the design of FPLA prototype with desired operating performances. The parameters of the numerical simulation program were amended by comparing the simulated incylinder pressure with experimentally derived data. At last CFD calculation of the combustion process was carried out to verify the effects of translator ignition position with two kinds of typical effective stroke length to bore ratios.
According to the dimensionless results, it can be concluded that:
For FPLA, a much narrow range of low operating speeds is expected to be utilized, which is due to the electrical generating scheme employed by the device. Therefore, a bigger stroke to bore ratio is favorable to decrease the to and fro frequency of the translator.
According to the load of FPLA, efficient power generation will be achieved by operating at a fixed oscillating rate. With smaller effective stroke length to bore ratio, bigger load coefficient is advantageous to achieve a higher effective efficiency while smaller load coefficient would lead to higher effective efficiency with bigger effective stroke length to bore ratio. Smaller load coefficient would lead to higher effective power output.
It has been found that an optimum ignition advance is available for the freepiston engine to achieve its best performance since earlier ignition is associated with more negative work in the compression stroke and a later ignition is associated with low peak incylinder pressures.
The efficiency of the engine is mainly associated with the proportion of the energy released before TDC which is associated with negative work to stop the translator. With a longer effective stroke length to bore ratio it is recommended to postpone the ignition timing to achieve a good performance of the engine.
According to the CFD calculated results with typical effective stroke length to bore ratio and ignition timings, the dimensionless results were reasonable.
References
 1.
Hannson J. Leksell M. Carlesson F.  2.
Mikalsen [. R. Roskilly A. P. The control of a freepiston engine generator. Part 2: Engine dynamics and piston motion control . Appl Energy2009 doi: 10.1016/ j.apenergy.2009.06.035  3.
Goertz M. Peng L. X.  4.
Atkinson C. Petreanu S. Clark N. N. Atkinson R. J.  5.
Shoukry E. Taylor S. Clark N. Numerical simulation for parametric . SAE paper 2002 011739, 2002  6.
Max ,E.  7.
Blarigan P. V. Paradiso N. Goldsborough S. S.  8.
Blarigan P. V. 2002 U.S. hydrogen program review, NREL/CP61032405, 2002  9.
Fredrisksson J. Denbratt I. 2004 011871, 2004  10.
Nemecek P. Vysoky O.  11.
Mikalsen R. Roskilly A. P. a twostroke free piston engine for electric power generation . Appl. Therm. Eng.2007 doi: 10.1061/j.applthermaleng. 2007.04.009  12.
Mikalsen ,R .Roskilly ,A.P .  13.
Xiao J. et al. Motion characteristic of a free piston linear engine . Appl Energy2009 doi:10.1016/j.apenergy.2009.07.005  14.
Cawthorne W. R. Famouri P. Chen J. D. Development of a linear alternatorengine for hybrid electric vehicle application.  15.
Wang ,J.B .,Howe ,H .  16.
Deng Z. Bold I. Nasar S. A. Fields in permanent magnet linear synchronous machines . IEEE Transactions on magnets.MAG22 22 2 1986  17.
Němeček P. Vysoký O. 2004EO  18.
Caresana F. Comodi G. Pelagalli L. 2004 320037, 2004  19.
Hohenberg G. F. Advanced approaches for heat transfer calculations. SAE Special Publications.449 61 79 , 1979  20.
Stone R. Introduction to internal combustion engine. 0768004950 of Automotive Engineers, Inc. Warrendale,1999  21.
Nagy C. T. Linear engine development for series hybrid electrical vehicl es. Dissertation, West Virginia: West Virginia University,2004  22.
Buckingham ,Edgar 1914 On Physically Similar Systems: Illustrations of the Use of Dimensional Analysis. Phys. Rev. 4: 345.doi:10.1103/PhysRev.4.345  23.
Goldsborough ,S.S. ,Blarigan ,P.V.  24.
Goldsborough S. S. Blarigan P. V. 2003  25.
Bergman M. Fredriksson J. Golovitchev V. I. 2008 012423, 2008