Preliminary Development of a Free Piston Expander–linear Generator for Small-scale Organic Rankine Cycle (orc) Waste Heat Recovery System

A novel free piston expander-linear generator (FPE-LG) integrated unit was proposed to recover waste heat efficiently from vehicle engine. This integrated unit can be used in a small-scale Organic Rankine Cycle (ORC) system and can directly convert the thermodynamic energy of working fluid into electric energy. The conceptual design of the free piston expander (FPE) was introduced and discussed. A cam plate and the corresponding valve train were used to control the inlet and outlet valve timing of the FPE. The working principle of the FPE-LG was proven to be feasible using an air test rig. The indicated efficiency of the FPE was obtained from the p–V indicator diagram. The dynamic characteristics of the in-cylinder flow field during the intake and exhaust processes of the FPE were analyzed based on Fluent software and 3D numerical simulation models using a computation fluid dynamics method. Results show that the indicated efficiency of the FPE can reach 66.2% and the maximal electric power output of the FPE-LG can reach 22.7 W when the working frequency is 3 Hz and intake pressure is 0.2 MPa. Two large-scale vortices are formed during the intake process because of the non-uniform distribution of velocity and pressure. The vortex flow will convert pressure energy and kinetic energy into thermodynamic energy for the working fluid, which weakens the power capacity of the working fluid.


Introduction
The internal combustion (IC) engine, which is the power system of an automobile, consumes a tremendous amount of energy.As energy shortages are continuously worsening, the energy-efficient utilization of IC engines has attracted increasing attention from researchers for a long time now.Only one-third of the energy from fuel combustion is utilized effectively in IC engines, whereas the rest is wasted in various forms of cooling, exhaust, and mechanical losses; exhaust energy accounts for the largest proportion of wasted energy in IC engines [1][2][3][4].Therefore, recovering exhaust energy is an effective means to save vehicle energy.In the case of waste heat recovery (WHR), the Organic Rankine Cycle (ORC) system is an effective technical solution and a promising means for industrialization of energy savings [5][6][7].To date, research on ORC applications has mostly focused on large-scale heat sources, such as solar energy, geothermal energy, and industrial waste heat, etc. [8][9][10][11][12].The lack of a suitable expander has caused difficulties in building corresponding ORC systems for recovering 2 of 18 small-scale waste heat, such as vehicle exhaust energy, etc.At present, developing an ideal small-scale expander or power output system that can convert waste heat into mechanical energy or electric energy to recover exhaust energy from vehicle IC engines is a crucial objective of ORC system research [13][14][15].
An expander is a key component in an ORC system.Compared with a centrifugal expander, a positive displacement expander is characterized by low flow rate, large expansion ratio, and low rotational speed [15][16][17].A positive displacement expander should therefore be more suitable for a small-scale ORC system [18,19].However, low cost and low rotational speed scroll expanders are frequently used in 0.1-2 kW experimental systems, whereas their performance in situations involving higher power has not been verified yet [20][21][22].A piston expander has higher thermal efficiency for small steam flow rate conditions and a good power output/size ratio.This type of expander is generally built or renovated based on existing commercial IC engines.If the pressure ratio on the expander is high, then piston expanders are suitable.Moreover, these expanders are more robust than scroll expanders.Gao et al. [23] constructed a mathematical model to evaluate the performance of a Rankine cycle system with a reciprocating piston expander and conducted a preliminary experiment.The results showed that the introducing a heat recovery system could increase engine power output by 12% when a diesel engine was operated at 80 kW/2590 rpm.Glavatskaya et al. [24] investigated the performance of a reciprocating expander for an automotive WHR application.The results showed that a maximum power output of 7 kW could be obtained at a high operating point, and the isentropic efficiency of the reciprocating expander varied from 55% to 70%.Chiong et al. [25] presented a concept for a new piston expander utilizing a nozzle as part of a secondary steam cycle to recover exhaust energy.Through simulation, the nozzle piston expander was found to increase the output power from a minimum of 0.73 kW to a maximum of 4.75 kW.The aforementioned survey of Rankine cycle systems applied to passenger vehicles indicates that the piston expander is a promising technology for WHR applications compared with other types of expander.This technology is innovative in the automotive industry [23][24][25][26].However, a comparatively large weight and a complicated structure restrict the extensive application of piston expanders in small-scale vehicle ORC systems.
To overcome the aforementioned shortcomings of piston expanders, the concept of a free piston expander (FPE), which is characterized by cancelling the crank-link mechanism, has been presented [27][28][29][30][31][32].Zhang et al. [27] developed a FPE that could be used in the transcritical CO 2 refrigeration cycle to replace the throttling valve.An expansion machine was integrated into a compressor in this device, and thus, the crankshaft became unnecessary.Moreover, a slide valve was designed to replace the traditional inlet/outlet valves to control the charge/discharge process.The design efficiencies of the expander and the compressor were assumed to be 60% and 70%, respectively; the experimental studies using p-V diagrams showed that device isentropic efficiency could reach 62%.Han et al. [28,29] constructed a simulation model and a test bench for a similar Rankine cycle WHR system using a FPE.Subsequently, a preliminary experiment was conducted.The power output of the free piston mechanism was eventually consumed through hydraulic buffers and reached 2.08 kJ, whereas piston displacement was 52 mm when the evaporator outlet pressure of the working fluid was 0.7 MPa.Weiss et al. [30] designed a FPE to produce power using a low-temperature energy source, and the FPE model was studied under various conditions.The results showed that decreasing piston mass reduced piston stroke length, but increased operation frequency and power output.Power output and energy recovery efficiency reached 25.6 mW and 80%-90%, respectively.Champagne et al. [31] conducted a preliminary experimental analysis of a small-scale FPE.The results showed that lubricants significantly affected the seal and leakage of a FPE.A thick lubricant seal was better in static configurations than in dynamic testing.Nickl et al. [32] developed a free piston compressor-expander unit with three expansion stages, and the isentropic efficiency of this device could reach 70% during the tests.
In the present research, a novel free piston expander-linear generator (FPE-LG) is proposed based on previous design experiences and the considerable literature on FPE prototypes mentioned earlier [27][28][29][30][31][32].Compared with the ORC power generation system that uses a conventional expander, Energies 2016, 9, 300 3 of 18 an ORC power generation system based on FPE-LG exhibits certain advantages such as compact structure, operation flexibility, and effective power output under the action of a two-phase flow in a partial load condition [30,31].In addition, FPE-LG is more suitable for ethanol and the water Rankine cycle because of its large expansion ratio.A small-scale ORC system with FPE-LG can transform waste heat from IC engine exhaust into electricity and subsequently supply power to the powertrain or auxiliary load of the vehicle.
In this research, a novel FPE-LG was developed, and several preliminary experiments were performed using an air test rig.For the intake and exhaust processes of the FPE, the dynamic characteristics of the in-cylinder flow field were analyzed using a computational fluid dynamics (CFD) method and Fluent software.

Working Principle and Function of the Free Piston Expander-Linear Generator
The FPE-LG mainly consists of two FPEs and a linear generator (LG).The LG mover is connected to the piston rods through flexible joints.This LG mover executes a reciprocating linear movement because of the gas pressure difference between the top of the two free pistons.The kinetic energy of the mover can be converted into electric energy by the LG.The valve train of the FPE is driven by a servo motor.In a formal engine application, the operation power of the servo motor will be provided by the engine-generated electrical output.When the FPE is in an intake-expansion stroke, the piston moves from top dead center (TDC) to bottom dead center (BDC).The working fluid flows first into the cylinder (i.e., intake process, the inlet valve is open and the outlet valve is closed), and then expands in the cylinder (i.e., expansion process, the inlet valve and outlet valves are both closed).When the FPE is in an exhaust stroke, the piston moves from BDC to TDC, and then the working fluid flows out of the cylinder (i.e., exhaust process, the inlet valve is closed and the outlet valve is open).When FPE A is in intake-expansion stroke, FPE B is in exhaust stroke, and vice versa.The LG is the sole energy output device of the FPE-LG, and thus, selecting an appropriate LG in the conceptual design stage is crucial.
For the same geometry size, a flat permanent magnet LG is easily constructed and exhibits good excitation; meanwhile, its output voltage, current effective value, efficiency, and specific power are slightly higher than those of the tubular permanent magnet LG with the same load [33,34].Therefore, a flat permanent magnet LG is selected to match with the FPEs in this research.The manufactured FPE-LG prototype is shown in Figure 1.
Energies 2016, 9, 300 3 of 17 an ORC power generation system based on FPE-LG exhibits certain advantages such as compact structure, operation flexibility, and effective power output under the action of a two-phase flow in a partial load condition [30,31].In addition, FPE-LG is more suitable for ethanol and the water Rankine cycle because of its large expansion ratio.A small-scale ORC system with FPE-LG can transform waste heat from IC engine exhaust into electricity and subsequently supply power to the powertrain or auxiliary load of the vehicle.
In this research, a novel FPE-LG was developed, and several preliminary experiments were performed using an air test rig.For the intake and exhaust processes of the FPE, the dynamic characteristics of the in-cylinder flow field were analyzed using a computational fluid dynamics (CFD) method and Fluent software.

Working Principle and Function of the Free Piston Expander-Linear Generator
The FPE-LG mainly consists of two FPEs and a linear generator (LG).The LG mover is connected to the piston rods through flexible joints.This LG mover executes a reciprocating linear movement because of the gas pressure difference between the top of the two free pistons.The kinetic energy of the mover can be converted into electric energy by the LG.The valve train of the FPE is driven by a servo motor.In a formal engine application, the operation power of the servo motor will be provided by the engine-generated electrical output.When the FPE is in an intake-expansion stroke, the piston moves from top dead center (TDC) to bottom dead center (BDC).The working fluid flows first into the cylinder (i.e., intake process, the inlet valve is open and the outlet valve is closed), and then expands in the cylinder (i.e., expansion process, the inlet valve and outlet valves are both closed).When the FPE is in an exhaust stroke, the piston moves from BDC to TDC, and then the working fluid flows out of the cylinder (i.e., exhaust process, the inlet valve is closed and the outlet valve is open).When FPE A is in intake-expansion stroke, FPE B is in exhaust stroke, and vice versa.The LG is the sole energy output device of the FPE-LG, and thus, selecting an appropriate LG in the conceptual design stage is crucial.
For the same geometry size, a flat permanent magnet LG is easily constructed and exhibits good excitation; meanwhile, its output voltage, current effective value, efficiency, and specific power are slightly higher than those of the tubular permanent magnet LG with the same load [33,34].Therefore, a flat permanent magnet LG is selected to match with the FPEs in this research.The manufactured FPE-LG prototype is shown in Figure 1.Energies 2016, 9, 300 4 of 18 Figure 2 presents the schematic diagram of a small-scale ORC system using the FPE-LG for recovering exhaust energy from vehicle engines.The key components are the evaporator, the FPE-LG, the condenser, the working fluid tank, and the working fluid pump.During the ORC WHR system operation process, the organic working fluid is pressurized by the pump and absorbs heat from low-medium grade waste heat source (i.e., engine exhaust) in the evaporator.Then, this fluid turns into high pressure and high temperature vapor and flows into the FPE-LG.This unit can convert the thermodynamic energy of the working fluid into electric energy (electric power can be consumed by the load resistance).After expanding in the FPE, the low pressure and low temperature vapor flows into the condenser.Then, this vapor turns into liquid and flows into the tank.The overall operation performance of the FPE-LG changes with working conditions, such as evaporation pressure, condensation pressure, and load resistance, etc.As shown in Figure 2, on the one hand, FPE converts the energy stored in the working fluid that enters the FPE into the kinetic energy of the free piston.The piston-mover group is driven to perform an oscillating linear movement according to the energy conservation law (various losses are neglected).On the other hand, the electromotive force is in proportion to the motion speed of the piston-mover group according to Faraday's law of electromagnetic induction.Hence, a reasonable conceptual design of the FPE is beneficial for reducing losses in the energy conversion process.Such a design also significantly influences the overall energy conversion efficiency of the FPE-LG.We primarily focus on the FPE in this research because of the aforementioned reasons.
Energies 2016, 9, 300 4 of 17 Figure 2 presents the schematic diagram of a small-scale ORC system using the FPE-LG for recovering exhaust energy from vehicle engines.The key components are the evaporator, the FPE-LG, the condenser, the working fluid tank, and the working fluid pump.During the ORC WHR system operation process, the organic working fluid is pressurized by the pump and absorbs heat from lowmedium grade waste heat source (i.e., engine exhaust) in the evaporator.Then, this fluid turns into high pressure and high temperature vapor and flows into the FPE-LG.This unit can convert the thermodynamic energy of the working fluid into electric energy (electric power can be consumed by the load resistance).After expanding in the FPE, the low pressure and low temperature vapor flows into the condenser.Then, this vapor turns into liquid and flows into the tank.The overall operation performance of the FPE-LG changes with working conditions, such as evaporation pressure, condensation pressure, and load resistance, etc.As shown in Figure 2, on the one hand, FPE converts the energy stored in the working fluid that enters the FPE into the kinetic energy of the free piston.The piston-mover group is driven to perform an oscillating linear movement according to the energy conservation law (various losses are neglected).On the other hand, the electromotive force is in proportion to the motion speed of the piston-mover group according to Faraday's law of electromagnetic induction.Hence, a reasonable conceptual design of the FPE is beneficial for reducing losses in the energy conversion process.Such a design also significantly influences the overall energy conversion efficiency of the FPE-LG.We primarily focus on the FPE in this research because of the aforementioned reasons.

Configuration and Main Parameters
According to a large number of measured experimental data, the restricted conditions for the FPE design are listed as follows: (1) Intake temperature of the working fluid is 300-350 K; (2) Intake pressure of the working fluid is 0.8-1 MPa; (3) Maximal expansion ratio of the FPE is 9; (4) Intake volume flow rate of the working fluid is 10-20 m 3 /h; (5) Exhaust back pressure is 0.11-0.13MPa.
Figure 3 shows the configuration diagram of the manufactured FPE.The main components include: valve train, piston, piston rod, cylinder cover, cylinder block, cylinder head, and rear cover, etc.The inlet and outlet valves are normally closed because of the precompression force of the valve spring.According to the precompression and stiffness of the valve spring as well as the maximal working frequency (8 Hz) of the FPE, the torque output of the servo motor should be higher than 10 N•m to open the inlet (or outlet) valve reliably during the intake process (or exhaust process).Thus, the specific servo motor (operating voltage is 200 V, rated torque output is 15 N•m) is used to drive

Configuration and Main Parameters
According to a large number of measured experimental data, the restricted conditions for the FPE design are listed as follows: (1) Intake temperature of the working fluid is 300-350 K; (2) Intake pressure of the working fluid is 0.8-1 MPa; (3) Maximal expansion ratio of the FPE is 9; (4) Intake volume flow rate of the working fluid is 10-20 m 3 /h; (5) Exhaust back pressure is 0.11-0.13MPa.
Figure 3 shows the configuration diagram of the manufactured FPE.The main components include: valve train, piston, piston rod, cylinder cover, cylinder block, cylinder head, and rear cover, etc.The inlet and outlet valves are normally closed because of the precompression force of the valve spring.According to the precompression and stiffness of the valve spring as well as the maximal working Energies 2016, 9, 300 5 of 18 frequency (8 Hz) of the FPE, the torque output of the servo motor should be higher than 10 N¨m to open the inlet (or outlet) valve reliably during the intake process (or exhaust process).Thus, the specific servo motor (operating voltage is 200 V, rated torque output is 15 N¨m) is used to drive the valve train.The servo motor is connected with the drive shaft of the valve train through a flexible shaft coupling.The piston rod of the FPE is connected with the mover of the LG through flexible joints.The specific design objectives of prototype performance are determined by considering the requirements for related projects and current experimental research status, and are listed as follows: output power of the FPE-LG is 1 kW and its maximal working frequency is 8 Hz.Hence, the main structure parameters of the FPE-LG prototype were determined according to the design method introduced in [35], as listed in Table 1: Energies 2016, 9, 300 5 of 17 the valve train.The servo motor is connected with the drive shaft of the valve train through a flexible shaft coupling.The piston rod of the FPE is connected with the mover of the LG through flexible joints.The specific design objectives of prototype performance are determined by considering the requirements for related projects and current experimental research status, and are listed as follows: output power of the FPE-LG is 1 kW and its maximal working frequency is 8 Hz.Hence, the main structure parameters of the FPE-LG prototype were determined according to the design method introduced in [35], as listed in Table 1:  As mentioned earlier, piston expanders have several drawbacks that will restrict their extensive application in small-scale ORC systems.In comparison, the FPE-LG is an electrically and mechanically integrated unit that does not require a crankshaft.The output power of a single reciprocating piston expander may be larger than that of the FPE-LG.In its final form, however, the FPE-LG will be designed to be modular, and thus, additional installation space is unnecessary.The number of these units can be adjusted flexibly according to the output power requirement for smallscale vehicle ORC systems.In addition, this FPE-LG is the first generation prototype.We will increase the working frequency gradually in a further study based on the preliminary results of motion characteristic in the current research.The final design objective of the working frequency ranges from 30 Hz to 40 Hz.As mentioned earlier, piston expanders have several drawbacks that will restrict their extensive application in small-scale ORC systems.In comparison, the FPE-LG is an electrically and mechanically integrated unit that does not require a crankshaft.The output power of a single reciprocating piston expander may be larger than that of the FPE-LG.In its final form, however, the FPE-LG will be designed to be modular, and thus, additional installation space is unnecessary.The number of these units can be adjusted flexibly according to the output power requirement for small-scale vehicle ORC systems.In addition, this FPE-LG is the first generation prototype.We will increase the working frequency gradually in a further study based on the preliminary results of motion characteristic in the current research.The final design objective of the working frequency ranges from 30 Hz to 40 Hz.

Energy Loss Mechanisms
The FPE-LG is an ideal energy conversion device for a small-scale ORC WHR system because of its characteristics of less mechanical loss and compact configuration.Compared with the ideal working process of the FPE, several primary loss mechanisms are involved in the actual operation process of the FPE.Such mechanisms can be summarized as follows: (1) Incomplete expansion.During the expansion process, the expansion of the working fluid in the cylinder is incomplete because of the limitation in actual piston stroke and valve timing.In-cylinder pressure is higher than exhaust back pressure at the end of expansion process; (2) Heat transfer.The working process of the FPE is not adiabatic because of the temperature difference between the FPE cylinder (or working fluid) and the ambient environment.Thus, heat exchange occurs between the FPE cylinder (or working fluid) and the ambient environment; (3) Flow loss during intake and exhaust processes.The inlet and outlet valves are used to control inlet and outlet valve timing, respectively; hence, flow loss occurs in the actual intake and exhaust processes; (4) Mechanical friction.The primary mechanical friction occurs in the following contact segments: piston and cylinder wall, piston rod and guide sleeve, drive shaft and bearing/cover, cam plate and valve slider, etc.A wear-resistant material with high strength (e.g., aluminum alloy and aluminum bronze) is used in designing the FPE process to reduce wear and improve the fatigue life of relevant components; (5) Leakage.When the cylinder is full of working fluid, there exists more or less leakage in the gap between the piston and cylinder wall.

Lubrication and Sealing
The maximal in-cylinder pressure of the FPE is limited to 0.8-1 MPa.The working fluid may leak through various gaps during the FPE operation process, including the gaps between the piston and cylinder wall, between the cylinder block and cylinder cover, and between the cylinder block and cylinder head.Hence, ring grooves on both ends of the cylinder are designed to install sealing rings constituted by polytetrafluoroethylene (PTFE).The piston ring is also constituted by PTFE.PTFE can reduce the friction between the sealing ring and the corresponding contact segments as well as perform a sealing function because of its self-lubricating property.Mechanical wear mainly occurs in the valve train and in the piston-mover group during the operation process of the FPE.The lubricating oil is injected through the oil holes processed on the rear cover and into the enclosed chamber between the rear cover and the cylinder head, where the cam plate of the valve train and some moving parts are located.Therefore, the corresponding contact segments can be lubricated through splash lubrication when the cam plate rotates during the operation process of the FPE.The piston rod performs a rapid reciprocating movement when passing through the guide sleeve.Meanwhile, installation errors exist in the alignment, which can lead to friction and wear.Hence, aluminum bronze is selected as the constituent material of the guide sleeve.Heat treatment technology is also adopted to satisfy the requirements of strength and reliability.Meanwhile, lubricating grease is inserted into the inner hole of the guide sleeve.This design is beneficial for reducing friction and wear.

Valve Train
The valve train signifcantly influences the stable operation of the FPE and flow loss during its intake and exhaust processes [36,37].Two kinds of valve trains can be selected for the FPE, namely, mechanical control valve manner and electromagnetic control valve manner.The former is used in this research after considering the manufacturing cost and working frequency of the FPE [37].As shown in Figure 4, a mechanical control valve train using a rotational cam plate has been developed.The specific servo motor (operating voltage is 200 V, rated torque output is 15 N¨m) is utilized to drive the valve train.The output end of the servo motor is connected to the drive shaft of the valve train through a flexible shaft coupling.Figure 4 illustrates the working process of the FPE and the corresponding behavior of the valve train.The entire working process (cycle) of the FPE includes two strokes: intake-expansion and exhaust.An intake-expansion stroke includes the intake and expansion processes, whereas an exhaust stroke includes only the exhaust process.For the initial position, as shown in Figure 4a, under the action of return force of the valve spring, the inlet and outlet valves are both closed, and the piston-mover group is located at TDC.During the intake process, as shown in Figure 4b, the cam plate axially lifts up the inlet valve slider and then opens the inlet valve by overcoming the return force of the valve spring.The working fluid flows into the cylinder through the inlet port and the inlet valve; the piston-mover group moves from TDC to BDC.During the expansion process, as shown in Figure 4c, the inlet and outlet valves are both closed under the action of the return force of the valve spring.The working fluid expands in the cylinder and pushes the piston-mover group to BDC.During the exhaust process, as shown in Figure 4d, the cam plate axially lifts up the outlet valve slider and then opens the outlet valve by overcoming the return force of the valve spring.The pistonmover group moves from BDC to TDC, and working fluid flows out of the cylinder and through the outlet port and the outlet valve.

Air Test Rig of the Free Piston Expander-Linear Generator
As the manufactured FPE-LG is a novel energy conversion device for a small-scale ORC WHR system, the air test rig is first established to verify the working principle and feasibility of the FPE-LG.According to the test results, FPE-LG may be subsequently integrated into the ORC WHR system after implementing several improvements (or adjustments) to the FPE-LG.Figure 5 shows the air test rig for the FPE-LG.This rig mainly includes an air compressor, a compressed air tank, the FPE-LG, an electric load, different kinds of sensors, data acquisition system, etc. Compressed air is used as the working fluid in the experiment because of two main reasons [38].First, supplying and consuming compressed air in the experiment are convenient.Moreover, the properties of compressed air are close to those of ideal gas, and compressed air does not exhibit any phase change in the experiment.Second, although the properties of compressed air and the vapor state of an organic working fluid are significantly different, this difference only has a slight influence on the test results of the feasibility Figure 4 illustrates the working process of the FPE and the corresponding behavior of the valve train.The entire working process (cycle) of the FPE includes two strokes: intake-expansion and exhaust.An intake-expansion stroke includes the intake and expansion processes, whereas an exhaust stroke includes only the exhaust process.For the initial position, as shown in Figure 4a, under the action of return force of the valve spring, the inlet and outlet valves are both closed, and the piston-mover group is located at TDC.During the intake process, as shown in Figure 4b, the cam plate axially lifts up the inlet valve slider and then opens the inlet valve by overcoming the return force of the valve spring.The working fluid flows into the cylinder through the inlet port and the inlet valve; the piston-mover group moves from TDC to BDC.During the expansion process, as shown in Figure 4c, the inlet and outlet valves are both closed under the action of the return force of the valve spring.The working fluid expands in the cylinder and pushes the piston-mover group to BDC.During the exhaust process, as shown in Figure 4d, the cam plate axially lifts up the outlet valve slider and then opens the outlet valve by overcoming the return force of the valve spring.The piston-mover group moves from BDC to TDC, and working fluid flows out of the cylinder and through the outlet port and the outlet valve.

Air Test Rig of the Free Piston Expander-Linear Generator
As the manufactured FPE-LG is a novel energy conversion device for a small-scale ORC WHR system, the air test rig is first established to verify the working principle and feasibility of the FPE-LG.According to the test results, FPE-LG may be subsequently integrated into the ORC WHR system after implementing several improvements (or adjustments) to the FPE-LG.Figure 5 shows the air test rig for the FPE-LG.This rig mainly includes an air compressor, a compressed air tank, the FPE-LG, an electric load, different kinds of sensors, data acquisition system, etc. Compressed air is used as the working fluid in the experiment because of two main reasons [38].First, supplying and consuming compressed air in the experiment are convenient.Moreover, the properties of compressed air are close to those of ideal gas, and compressed air does not exhibit any phase change in the experiment.Second, although the properties of compressed air and the vapor state of an organic working fluid are significantly different, this difference only has a slight influence on the test results of the feasibility of the FPE-LG as well as the valve train.

Experimental Procedure
The flow rate of the working fluid can be regulated by adjusting the opening degree of the ball valves.The temperature sensors and pressure sensors are installed into the two FPEs.A displacement sensor is also installed and used to indicate the displacement variation of the piston-mover group during the working process of the FPE-LG.A phase difference of 180° is observed between the cam plates of the two FPEs under the control of the two servo motors, which ensures that the working process and the piston stroke of the two FPEs are coordinated.A flowmeter is installed to measure the total flow rate of the working fluid for the two FPEs.The electrical output of the LG can be converted to direct current electricity by using the three-phase bridge rectifier and the smoothing capacitor.The main test instrumentations of the experimental procedure are listed in Table 2.

Three-Dimensional Numerical Simulation of the Free Piston Expander
To explore the dynamic characteristics of the in-cylinder flow field during the intake and exhaust processes of the FPE, three-dimensional (3D) numerical simulation of the FPE is implemented by using Fluent V14.0 based on the control volume methodology of the combined dynamic meshing technique [39].

Grid Generation
The 3D grid is built using Gambit 6.0 because of the complex geometric model established in Catia V5R21.The mesh model of the FPE is shown in Figure 6.The computational grid of the FPE consists of 2,031,293 elements and 3,937,918 nodes.The internal calculation fluid region of the FPE is composed of three sub-fluid regions, which are separated by the inlet and outlet valves.The unstructured tetrahedral grid is used to mesh the inlet and outlet fluid regions (fluid regions B and C, respectively, in Figure 6) with the smallest unit size of 0.2 mm.The submap structured hexahedral

Experimental Procedure
The flow rate of the working fluid can be regulated by adjusting the opening degree of the ball valves.The temperature sensors and pressure sensors are installed into the two FPEs.A displacement sensor is also installed and used to indicate the displacement variation of the piston-mover group during the working process of the FPE-LG.A phase difference of 180 ˝is observed between the cam plates of the two FPEs under the control of the two servo motors, which ensures that the working process and the piston stroke of the two FPEs are coordinated.A flowmeter is installed to measure the total flow rate of the working fluid for the two FPEs.The electrical output of the LG can be converted to direct current electricity by using the three-phase bridge rectifier and the smoothing capacitor.The main test instrumentations of the experimental procedure are listed in Table 2.

Three-Dimensional Numerical Simulation of the Free Piston Expander
To explore the dynamic characteristics of the in-cylinder flow field during the intake and exhaust processes of the FPE, three-dimensional (3D) numerical simulation of the FPE is implemented by using Fluent V14.0 based on the control volume methodology of the combined dynamic meshing technique [39].

Grid Generation
The 3D grid is built using Gambit 6.0 because of the complex geometric model established in Catia V5R21.The mesh model of the FPE is shown in Figure 6.The computational grid of the FPE consists of 2,031,293 elements and 3,937,918 nodes.The internal calculation fluid region of the FPE is composed of three sub-fluid regions, which are separated by the inlet and outlet valves.The unstructured tetrahedral grid is used to mesh the inlet and outlet fluid regions (fluid regions B and C, respectively, in Figure 6) with the smallest unit size of 0.2 mm.The submap structured hexahedral mesh with the minimum unit size of 0.4 mm is adapted to the in-cylinder fluid region (fluid region A in Figure 6) because of its regular cylindrical shape.The valves and the piston are meshed by the unstructured tetrahedral grid with the smallest unit sizes of 0.2 mm and 0.4 mm, respectively.The gap between the valve and the valve seat is 0.1 mm.The grid of the sub-fluid regions near the valves is refined with the minimum unit size of 0.1 mm after considering the dramatic change in the flow field.Interface technology is applied to ensure the variable interpolation and delivery of the overlapping segments in different fluid regions.

Boundary Conditions and Dynamic Mesh
The numerical simulation of unsteady flow in the FPE is implemented using the renormalization group-based k-ε turbulence model in Fluent V14.0 [39].The thermodynamic properties of the working fluid can be determined using REFPROP software developed by the National Institute of Standards and Technology (NIST).The calculation difficulties are attributed to the grid varying constantly as a result of the movement of the valves and the piston during the working process of the FPE.Thus, dynamic mesh definition is necessary for the aforementioned moving components.The motion characteristics of the valves can be converted into the equivalent relationship between valve displacement and elapsed time according to cam plate profile and working frequency.Consequently, a control program is compiled to regulate the dynamic grid for the valves.This control program can be derived from the experimental results of the piston movement and the force analysis results during the working process of the FPE through MATLAB Simulink.
The working conditions of the FPE for the numerical simulation are set as follows: intake pressure is 0.2 MPa, intake temperature is 300 K, and working frequency of the FPE-LG is 3 Hz.Air is used as working fluid in the CFD simulation model.The property parameters of air are obtained from the average temperature and average pressure of the inlet and outlet of the FPE using NIST-REFPROP.The outlet boundary condition is set to atmospheric pressure, whereas heat transfer between the cylinder wall and the ambient environment is disregarded.Simultaneously, the average time step, which can be adjusted to suit the moving situations of the valves and the piston, is set to 1 × 10 −4 s.In this manner, negative volume will not be generated when the dynamic mesh updates the position.Residual, pressure, and turbulence intensity variations can be monitored in real time to observe the convergence and results of the numerical calculation.

Unsteady Flow Behavior during the Intake and Exhaust Processes
Figure 7 shows the complex vortex flow in the cylinder when the rotational angle of the cam plate is 75° after the beginning of the curved surface for lifting the inlet valve during the intake process.At the beginning of the intake process, the working fluid with high pressure is injected into the cylinder through the narrow flow area between the inlet valve and the inlet valve seat.Thus, given the small opening degree of the inlet valve, the Laval nozzle effect may occur in the narrow flow area.Simultaneously, the narrow flow area between the inlet valve and the inlet valve seat exhibits evident pressure drop, throttling loss, and turbulence intensity.A local high-pressure zone

Boundary Conditions and Dynamic Mesh
The numerical simulation of unsteady flow in the FPE is implemented using the renormalization group-based k-ε turbulence model in Fluent V14.0 [39].The thermodynamic properties of the working fluid can be determined using REFPROP software developed by the National Institute of Standards and Technology (NIST).The calculation difficulties are attributed to the grid varying constantly as a result of the movement of the valves and the piston during the working process of the FPE.Thus, dynamic mesh definition is necessary for the aforementioned moving components.The motion characteristics of the valves can be converted into the equivalent relationship between valve displacement and elapsed time according to cam plate profile and working frequency.Consequently, a control program is compiled to regulate the dynamic grid for the valves.This control program can be derived from the experimental results of the piston movement and the force analysis results during the working process of the FPE through MATLAB Simulink.
The working conditions of the FPE for the numerical simulation are set as follows: intake pressure is 0.2 MPa, intake temperature is 300 K, and working frequency of the FPE-LG is 3 Hz.Air is used as working fluid in the CFD simulation model.The property parameters of air are obtained from the average temperature and average pressure of the inlet and outlet of the FPE using NIST-REFPROP.The outlet boundary condition is set to atmospheric pressure, whereas heat transfer between the cylinder wall and the ambient environment is disregarded.Simultaneously, the average time step, which can be adjusted to suit the moving situations of the valves and the piston, is set to 1 ˆ10 ´4 s.In this manner, negative volume will not be generated when the dynamic mesh updates the position.Residual, pressure, and turbulence intensity variations can be monitored in real time to observe the convergence and results of the numerical calculation.

Unsteady Flow Behavior during the Intake and Exhaust Processes
Figure 7 shows the complex vortex flow in the cylinder when the rotational angle of the cam plate is 75 ˝after the beginning of the curved surface for lifting the inlet valve during the intake process.At the beginning of the intake process, the working fluid with high pressure is injected into the cylinder through the narrow flow area between the inlet valve and the inlet valve seat.Thus, given the small opening degree of the inlet valve, the Laval nozzle effect may occur in the narrow flow area.Simultaneously, the narrow flow area between the inlet valve and the inlet valve seat exhibits evident pressure drop, throttling loss, and turbulence intensity.A local high-pressure zone A is formed in the half cavity around the inlet valve because of the inlet port wall.Simultaneously, gas flows along the valve axes with a tremendously high velocity in the cylinder, which results in low pressure in the gas flow core region.Thus, a local high-pressure zone B is formed in the small cylindrical cavity between the cylinder head and cylinder wall (below the inlet valve).The mainstream continues flowing along the tangential direction of the cylinder wall and the piston top.A local high-pressure zone C is formed in a short channel (between the outlet valve and the cylinder wall) below the outlet valve.Subsequently, gas flows toward the main intake stream along the tangential direction of the cylinder head surface because of the pressure gradient between the core region of the mainstream and zone C. Gas flow is merged into the newly injected intake flow because of viscous force.Thecore region of the vortices is another high-pressure zone (marked with the letter D in Figure 7).Two large-scale vortices V 1 and V 2 are formed during the intake process.
Energies 2016, 9, 300 10 of 17 tangential direction of the cylinder head surface because of the pressure gradient between the core region of the mainstream and zone C. Gas flow is merged into the newly injected intake flow because of viscous force.The core region of the vortices is another high-pressure zone (marked with the letter D in Figure 7).Two large-scale vortices V1 and V2 are formed during the intake process.Figure 8 shows the gas velocity distribution and streamlines during the intake process at different rotational angles of the cam plate after the beginning of the curved surface for lifting up the inlet valve (ranges from 2° to 88°).When the opening degree of the inlet valve is relatively small, a large flow resistance occurs, which results in the low pressure in the cylinder.Hence, gas flow around the inlet valve head is injected into the cylinder along the tangential direction of the inlet valve head through the narrow flow area between the inlet valve and the inlet valve seat because of the considerable pressure difference.Moreover, gas flow is separated into two branches after it comes in contact with the piston top and subsequently compresses the original gas in the cylinder respectively.Consequently, the distribution of velocity and pressure becomes non-uniform.As shown in Figure 8a,b, no obvious large-scale vortex is formed in the cylinder because of the relatively small cylinder volume.Subsequently, the volume of intake flow gradually increases with piston movement.Two large-scale vortices (V1 and V2) are formed in the cylinder, as shown in Figure 8c, because of the restriction imposed by the solid wall (which is constructed by the cylinder wall, piston top, and cylinder head surface).Vortices V1 and V2 are gradually enhanced with the increase in inlet valve opening degree because the volume of intake flow rapidly increases with the increase in the flow area between the inlet valve and the inlet valve seat, as shown in Figure 8d.When the inlet valve is fully opened, the gas continues flowing into the cylinder.However, the pressure difference between the inlet port and the cylinder becomes small.At this point, the distribution of gas velocity and pressure generally becomes relatively uniform in the cylinder, and two stable vortices are formed in the cylinder.The in-cylinder velocity field distribution is layered, and the gas flow near the solid wall has a high velocity and gradually spreads out.Irregular small-scale vortices are formed in the local high-pressure zone C below the outlet valve, as shown in Figure 8e, because of the decrease in overall vortex velocity.The overall velocity of gas flow around the intake valve increases again with the decrease in inlet valve opening degree at the end of intake process.Furthermore, vortices V1 and V2 Figure 8 shows the gas velocity distribution and streamlines during the intake process at different rotational angles of the cam plate after the beginning of the curved surface for lifting up the inlet valve (ranges from 2 ˝to 88 ˝).When the opening degree of the inlet valve is relatively small, a large flow resistance occurs, which results in the low pressure in the cylinder.Hence, gas flow around the inlet valve head is injected into the cylinder along the tangential direction of the inlet valve head through the narrow flow area between the inlet valve and the inlet valve seat because of the considerable pressure difference.Moreover, gas flow is separated into two branches after it comes in contact with the piston top and subsequently compresses the original gas in the cylinder respectively.Consequently, the distribution of velocity and pressure becomes non-uniform.As shown in Figure 8a,b, no obvious large-scale vortex is formed in the cylinder because of the relatively small cylinder volume.Subsequently, the volume of intake flow gradually increases with piston movement.Two large-scale vortices (V 1 and V 2 ) are formed in the cylinder, as shown in Figure 8c, because of the restriction imposed by the solid wall (which is constructed by the cylinder wall, piston top, and cylinder head surface).Vortices V 1 and V 2 are gradually enhanced with the increase in inlet valve opening degree because the volume of intake flow rapidly increases with the increase in the flow area between the inlet valve and the inlet valve seat, as shown in Figure 8d.When the inlet valve is fully opened, the gas continues flowing into the cylinder.However, the pressure difference between the inlet port and the Energies 2016, 9, 300 11 of 18 cylinder becomes small.At this point, the distribution of gas velocity and pressure generally becomes relatively uniform in the cylinder, and two stable vortices are formed in the cylinder.The in-cylinder velocity field distribution is layered, and the gas flow near the solid wall has a high velocity and gradually spreads out.Irregular small-scale vortices are formed in the local high-pressure zone C below the outlet valve, as shown in Figure 8e, because of the decrease in overall vortex velocity.The overall velocity of gas flow around the intake valve increases again with the decrease in inlet valve opening degree at the end of intake process.Furthermore, vortices V 1 and V 2 enter the dissipation stage as the overall distributions of velocity and pressure become uniform in the cylinder.Irregular small-scale vortices are formed below the inlet valve (i.e., the narrow channel between the inlet valve and the cylinder), as shown in Figure 8f.Vortices form in the cylinder, and they change with inlet valve opening degree and piston displacement during the intake process.Vortex flow converts pressure energy and kinetic energy into thermodynamic energy for the working fluid, which weakens the power capacity of the working fluid within a certain period [40,41].Several practical approaches should be adopted to adjust the intake flow and reduce the energy losses caused by large-scale vortex flow, and thus, improve the energy conversion efficiency of the FPE.Examples of such approaches include optimizing the inner structure of the cylinder, modifying inlet valve shape and size, and reducing the length of the inlet port.
Figure 9 shows the gas velocity distribution during the exhaust process at different cam plate rotational angles after the beginning of the curved surface for lifting the inlet valve (ranges from 185° to 275°).Compared with that in the intake process, flow velocity is lower during the exhaust process.Gas flow with high velocity mainly occurs in the flow area between the outlet valve and the outlet valve seat.For in-cylinder gas flow, compared with exhaust velocity, velocity of vortex flow can be ignored.Vortices form in the cylinder, and they change with inlet valve opening degree and piston displacement during the intake process.Vortex flow converts pressure energy and kinetic energy into thermodynamic energy for the working fluid, which weakens the power capacity of the working fluid within a certain period [40,41].Several practical approaches should be adopted to adjust the intake flow and reduce the energy losses caused by large-scale vortex flow, and thus, improve the energy conversion efficiency of the FPE.Examples of such approaches include optimizing the inner structure of the cylinder, modifying inlet valve shape and size, and reducing the length of the inlet port.
Figure 9 shows the gas velocity distribution during the exhaust process at different cam plate rotational angles after the beginning of the curved surface for lifting the inlet valve (ranges from 185 ˝to 275 ˝).Compared with that in the intake process, flow velocity is lower during the exhaust process.Gas flow with high velocity mainly occurs in the flow area between the outlet valve and the outlet valve seat.For in-cylinder gas flow, compared with exhaust velocity, velocity of vortex flow can be ignored.
The aforementioned work can provide a theoretical foundation to improve and optimize the conceptual design of next-generation FPEs.
Figure 9 shows the gas velocity distribution during the exhaust process at different cam plate rotational angles after the beginning of the curved surface for lifting the inlet valve (ranges from 185° to 275°).Compared with that in the intake process, flow velocity is lower during the exhaust process.Gas flow with high velocity mainly occurs in the flow area between the outlet valve and the outlet valve seat.For in-cylinder gas flow, compared with exhaust velocity, velocity of vortex flow can be ignored.

Analysis of the Actual Working Process of the Free Piston Expander
Figure 10 shows the variation of in-cylinder pressure and piston displacement at various working frequencies of FPE A when the intake pressure is 0.2 MPa and the intake temperature is 300 K.As shown in Figure 10, the actual stroke length of the piston decreases with the increase in working frequency.When working frequency is low, the piston-mover group moves rapidly from TDC to BDC and reaches BDC during the intake process because intake time is long and the electromagnetic force is insufficiently large.Then, the piston-mover group stops moving for a certain period until the exhaust stroke (process) begins.No effective expansion process actually occurs in the cylinder of the FPE, which results in an extremely low energy conversion efficiency.In addition, collision between the piston and the cylinder head may occur in the other FPE.This becomes the largest damage risk for such a system [42].Two possible solutions are proposed to solve the aforementioned problem.On the one hand, the working frequency of the FPE should be increased (higher than 2.5 Hz).On the other hand, the LG should be efficiently matched with the FPEs, and a reasonable valve timing and valve duration should be determined.

Analysis of the Actual Working Process of the Free Piston Expander
Figure 10 shows the variation of in-cylinder pressure and piston displacement at various working frequencies of FPE A when the intake pressure is 0.2 MPa and the intake temperature is 300 K.As shown in Figure 10, the actual stroke length of the piston decreases with the increase in working frequency.When working frequency is low, the piston-mover group moves rapidly from TDC to BDC and reaches BDC during the intake process because intake time is long and the electromagnetic force is insufficiently large.Then, the piston-mover group stops moving for a certain period until the exhaust stroke (process) begins.No effective expansion process actually occurs in the cylinder of the FPE, which results in an extremely low energy conversion efficiency.In addition, collision between the piston and the cylinder head may occur in the other FPE.This becomes the largest damage risk for such a system [42].Two possible solutions are proposed to solve the aforementioned problem.On the one hand, the working frequency of the FPE should be increased (higher than 2.5 Hz).On the other hand, the LG should be efficiently matched with the FPEs, and a reasonable valve timing and valve duration should be determined.Furthermore, a natural frequency of operation exists if the piston working frequency is increased from 3 Hz to a different value that eliminates the piston dwell time at top or bottom dead center.In reference [42], a simulation model on the FPLA (which is similar to the FPE-LG except for combustion process) is formed in a MATLAB/Simulink environment.According to the theoretical analysis, the system can be equivalent to a linear single-degree-of-freedom discrete system with variable damping Furthermore, a natural frequency of operation exists if the piston working frequency is increased from 3 Hz to a different value that eliminates the piston dwell time at top or bottom dead center.In reference [42], a simulation model on the FPLA (which is similar to the FPE-LG except for combustion process) is formed in a MATLAB/Simulink environment.According to the theoretical analysis, the system can be equivalent to a linear single-degree-of-freedom discrete system with variable damping and stiffness.It can be considered that the natural frequency is only decided by the mass and the equivalent stiffness.It can be written as: where P 0 is the cylinder pressure at the equilibrium position, L is the oscillation amplitude, γ is the polytropic exponent, A is the piston cross sectional area.For the FPLA system, the natural frequency has some relationship with P 0 and L. They are variables varying with the external excitation.All the variables can be set to the proper ranges.This means that the FPLA has a compression ratio range which can cover the full working range of this engine in the studied case.
As mentioned in reference [42], the piston motion frequency depends on intake pressure, moving mass, and loads, etc.It is not constrained.Thus, each combination of these parameters decides one working condition.Theoretically these parameters can be gradually adjusted to achieve certain frequency.However, at present, some mechanical wears exist between the cam plate and the valve sliders.Meanwhile, the value of external load resistance is small during the tests.Hence, from the view of safety, the working frequency (more than 3 Hz) has not been continued increasing in preliminary tests.This problem will be tried to resolve and the natural frequency will be estimated in the future work based on the method in reference [42] and more experimental studies.

Indicated Efficiency of the Free Piston Expander
In this research, the indicated efficiency of FPE A is analyzed based on the p-V indicator diagram of FPE A. The in-cylinder pressure variation of FPE A, as well as the comparison between the ideal working process and actual working process, is shown in Figure 11 when the working frequency is 3 Hz, the intake pressure is 0.2 MPa, and the load resistance is 5 Ω.
Energies 2016, 9, 300 13 of 17 variables can be set to the proper ranges.This means that the FPLA has a compression ratio range which can cover the full working range of this engine in the studied case.
As mentioned in reference [42], the piston motion frequency depends on intake pressure, moving mass, and loads, etc.It is not constrained.Thus, each combination of these parameters decides one working condition.Theoretically these parameters can be gradually adjusted to achieve certain frequency.However, at present, some mechanical wears exist between the cam plate and the valve sliders.Meanwhile, the value of external load resistance is small during the tests.Hence, from the view of safety, the working frequency (more than 3 Hz) has not been continued increasing in preliminary tests.This problem will be tried to resolve and the natural frequency will be estimated in the future work based on the method in reference [42] and more experimental studies.

Indicated Efficiency of the Free Piston Expander
In this research, the indicated efficiency of FPE A is analyzed based on the p-V indicator diagram of FPE A. The in-cylinder pressure variation of FPE A, as well as the comparison between the ideal working process and actual working process, is shown in Figure 11 when the working frequency is 3 Hz, the intake pressure is 0.2 MPa, and the load resistance is 5 Ω.In Figure 11, Process e1-a1-b1 is the actual intake process, Process e-a-b is the ideal intake process, Process b1-c is the actual expansion process, Process b-c is the ideal expansion process, Process c-e1 is the actual exhaust process, and Process c-d-e is the ideal exhaust process.We define the indicated efficiency of the FPE as the ratio of the area surrounded by the curves representing the actual working process to the area surrounded by the curves representing the ideal working process.Thereby, according to test data shown in Figure 11, the indicated efficiency of FPE A is 66.2%.Based on the p-V indicator diagram of the FPE, we define the indicated energy loss ratio of the FPE as the ratio of the area between the curves representing the actual working process and the curves representing the ideal working process to the area surrounded by the curves representing the ideal working process.In Figure 11, Process e 1 -a 1 -b 1 is the actual intake process, Process e-a-b is the ideal intake process, Process b 1 -c is the actual expansion process, Process b-c is the ideal expansion process, Process c-e 1 is the actual exhaust process, and Process c-d-e is the ideal exhaust process.We define the Energies 2016, 9, 300 14 of 18 indicated efficiency of the FPE as the ratio of the area surrounded by the curves representing the actual working process to the area surrounded by the curves representing the ideal working process.Thereby, according to test data shown in Figure 11, the indicated efficiency of FPE A is 66.2%.Based on the p-V indicator diagram of the FPE, we define the indicated energy loss ratio of the FPE as the ratio of the area between the curves representing the actual working process and the curves representing the ideal working process to the area surrounded by the curves representing the ideal working process.According to the test data shown in Figure 11, the total energy loss ratio of FPE A is 33.8%.The energy loss ratio that corresponds to intake process is 12.9% (i.e., the ratio of the area between the curves representing the actual intake process and the curves representing the ideal intake process to the area surrounded by the curves representing the ideal working process), the energy loss ratio that corresponds to the expansion process is 4.9% (i.e., the ratio of the area between the curves representing the actual expansion process and the curves representing the ideal expansion process to the area by the curves representing the ideal working process), and the energy loss ratio that corresponds to the exhaust process is 16% (i.e., the ratio of the area between the curves representing the actual exhaust process and the curves representing the ideal exhaust process to the area surrounded by the curves representing the ideal working process).
Figure 12a shows the p-V indicator diagram for FPE A at different intake pressures when the working frequency, the load resistance, and the intake temperature are 2 Hz, 5 Ω, and 300 K, respectively.We define variable V on the X-axis as sweep volume of the expander (i.e., total cylinder volume less clearance volume).The clearance volume includes two parts: the first part is formed with the cylinder head, the cylinder block, and the piston when the piston reached TDC; and the second part is formed with the cylinder cover, the cylinder block, and the piston when the piston reached BDC.With the increase of intake pressure, the power output and the actual piston stroke length both increase.The difference between the intake pressure and the peak value of the in-cylinder pressure also increases.As shown in Figure 12b, we define variable P in on the X-axis as the intake pressure of the expander.The indicated efficiency of FPE A decreases with the increase in intake pressure.The energy loss ratios that correspond to the intake and exhaust processes evidently increase.
Energies 2016, 9, 300 14 of 17 Figure 12a shows the p-V indicator diagram for FPE A at different intake pressures when the working frequency, the load resistance, and the intake temperature are 2 Hz, 5 Ω, and 300 K, respectively.We define variable V on the X-axis as sweep volume of the expander (i.e., total cylinder volume less clearance volume).The clearance volume includes two parts: the first part is formed with the cylinder head, the cylinder block, and the piston when the piston reached TDC; and the second part is formed with the cylinder cover, the cylinder block, and the piston when the piston reached BDC.With the increase of intake pressure, the power output and the actual piston stroke length both increase.The difference between the intake pressure and the peak value of the in-cylinder pressure also increases.As shown in Figure 12b, we define variable Pin on the X-axis as the intake pressure of the expander.The indicated efficiency of FPE A decreases with the increase in intake pressure.The energy loss ratios that correspond to the intake and exhaust processes evidently increase.

Electric Power Output of Free Piston Expander-Linear Generator
The variation tendency of the electric power output with the intake pressure is shown in Figure 13 when the load resistance and the working frequency of the FPE-LG are 5 Ω and 3 Hz, respectively.As can be seen in Figure 13, the electric power output of the FPE-LG generally increases substantially with the improvement of the intake pressure, but the increasing extent differs.The variation of the electric power output is small at low working frequencies, whereas the increase in the electric power output is comparatively obvious at high working frequencies.The maximal electric power output of the FPE-LG can reach 22.7 W when the intake pressure is 0.2 MPa.

Electric Power Output of Free Piston Expander-Linear Generator
The variation tendency of the electric power output with the intake pressure is shown in Figure 13 when the load resistance and the working frequency of the FPE-LG are 5 Ω and 3 Hz, respectively.As can be seen in Figure 13, the electric power output of the FPE-LG generally increases substantially with the improvement of the intake pressure, but the increasing extent differs.The variation of the electric power output is small at low working frequencies, whereas the increase in the electric power output is comparatively obvious at high working frequencies.The maximal electric power output of the FPE-LG can reach 22.7 W when the intake pressure is 0.2 MPa.

Electric Power Output of Free Piston Expander-Linear Generator
The variation tendency of the electric power output with the intake pressure is shown in Figure 13 when the load resistance and the working frequency of the FPE-LG are 5 Ω and 3 Hz, respectively.As can be seen in Figure 13, the electric power output of the FPE-LG generally increases substantially with the improvement of the intake pressure, but the increasing extent differs.The variation of the electric power output is small at low working frequencies, whereas the increase in the electric power output is comparatively obvious at high working frequencies.The maximal electric power output of the FPE-LG can reach 22.7 W when the intake pressure is 0.2 MPa.

Conclusions
Based on this initial work, the FPE-LG shows promising potential for continued development as a simple and efficient thermo-electric conversion unit.This research focuses on the development of

Conclusions
Based on this initial work, the FPE-LG shows promising potential for continued development as a simple and efficient thermo-electric conversion unit.This research focuses on the development of the FPE-LG, as well as on the conceptual design and 3D numerical simulation of the FPE.In addition, several preliminary experiments based on an air test rig have been conducted.The main research findings can be summarized as follows: (1) The working principle of the FPE-LG is proven feasible through the air test rig.However, further test of the FPE-LG in a small-scale ORC system should be performed in the future; (2) The energy-conversion efficiency of the expander is obviously affected by the intake parameters.
The indicated efficiency of the FPE can reach 66.2% and the maximal electric power output of the FPE-LG can reach 22.7 W when the working frequency is 3 Hz and the intake pressure is 0.2 MPa; (3) Two large-scale vortices are formed during the intake process.To improve the energy conversion efficiency of the FPE, several practical approaches should be adopted to adjust the intake flow and to reduce the energy losses caused by large-scale vortex flow.
Several conclusions may be drawn from this design study.First, a linear generator with a larger thrust force are desirable to increase output power of FPE-LG and avoid collisions between the piston and the cylinder head at low working frequencies.Second, the mechanical strength of related parts in the valve train needs to be improved to further increase the working frequency.Meanwhile, the cam plate intake duration needs to be decreased, whereas the duration of the expansion process needs to be increased.Third, in its final form, installation space is a premium consideration for the FPE-LG.This unit will be designed to be modular, and thus, the number of these units can be flexibly adjusted according to the output power requirement for a wide variety of small-scale vehicle ORC systems.
These results are greatly beneficial to next generation FPE-LG development.More experimental testing of the optimized FPE-LG performance will be conducted using the guidelines established in this work as a foundation.Future CFD simulation with organic working fluid will be conducted in an attempt to further reveal the in-cylinder flow characteristics.

Figure 3 .
Figure 3. Exploded view of the FPE configuration (small fasteners are omitted for clarity).

Figure 3 .
Figure 3. Exploded view of the FPE configuration (small fasteners are omitted for clarity).

Figure 4 .
Figure 4. Working process of FPE and corresponding behaviors of the valve train: (a) initial position; (b) intake process; (c) expansion process; (d) exhaust process.

Figure 4 .
Figure 4. Working process of FPE and corresponding behaviors of the valve train: (a) initial position; (b) intake process; (c) expansion process; (d) exhaust process.

Energies 2016, 9 , 300 8 of 17 Figure 5 .
Figure 5. Air test rig of the FPE-LG: (a) the schematic diagram of experimental system; (b) the picture of experimental system.

Figure 5 .
Figure 5. Air test rig of the FPE-LG: (a) the schematic diagram of experimental system; (b) the picture of experimental system.
valve and the valve seat is 0.1 mm.The grid of the sub-fluid regions near the valves is refined with the minimum unit size of 0.1 mm after considering the dramatic change in the flow field.Interface technology is applied to ensure the variable interpolation and delivery of the overlapping segments in different fluid regions.

Figure 6 .
Figure 6.Computational domains of the FPE.

Figure 6 .
Figure 6.Computational domains of the FPE.

Figure 7 .
Figure 7. Vortex flow in cylinder during the intake process: (a) Main view; (b) View 1: cross-section vertical to the inlet valve axis, 10 mm below the inlet valve; (c) View 2: cross-section vertical to the inlet valve axis, 50 mm below the inlet valve; (d) View 3: longitudinal-section along the valve axes.

Figure 7 .
Figure 7. Vortex flow in cylinder during the intake process: (a) Main view; (b) View 1: cross-section vertical to the inlet valve axis, 10 mm below the inlet valve; (c) View 2: cross-section vertical to the inlet valve axis, 50 mm below the inlet valve; (d) View 3: longitudinal-section along the valve axes.

Figure 11 .
Figure 11.p-V indicator diagram of the FPE A.

Figure 11 .
Figure 11.p-V indicator diagram of the FPE A.

Figure 12 .
Figure 12. p-V diagram and indicated efficiency for different intake pressures: (a) p-V diagram for different intake pressures; (b) indicated efficiency and energy losses.

Figure 12 .
Figure 12. p-V diagram and indicated efficiency for different intake pressures: (a) p-V diagram for different intake pressures; (b) indicated efficiency and energy losses.

Table 1 .
The main construction parameters of FPE.

Table 1 .
The main construction parameters of FPE.

Table 2 .
The main test instrumentations of the experimental procedure.

Table 2 .
The main test instrumentations of the experimental procedure.