Overall Performance Evaluation of Small Scale LNG Production Processes

The liquefied natural gas (LNG) is considered a viable solution to replace oil-based engines (common in heavy-duty truck and naval industry) reducing the environmental impact in the transport sector. Since liquefaction plants represent energy intensive processes, the best configurations/operation assessment is of primary importance. In this paper, a novel general procedure for the thermodynamic design and optimization, engineering design and off-design evaluation for small-scale LNG production systems is presented. The procedure can be used for the complete design and performance evaluation of plug & play facilities at filling stations for vehicles/boats, with the contemporary benefits of reducing pollutant emission in the city/port area and operating as electrical storage, coupled with renewable generators. Furthermore, the procedure has been applied to a case study (ferry boat operating at the main canal in the port of Ravenna, Italy), evaluating the optimal size for the integrated wind plant by minimizing the electricity introduction into the grid. The obtained results show 78 kW as optimal wind size, allowing the LNG plant to operate 187 h/year in design and 4720 h/year in off-design conditions, with electricity surplus around 33 MWh/year. A prototype will be installed to reduce pollutant emissions and test this technology as a storage option for renewable sources.


Introduction
It is well known that in decades to come, natural gas (NG) is expected to cover an increasing portion of the supply demand in all the energy fields. Among all fossil fuels, it is the one with the best combination of natural reserves and acceptable emissions levels [1]. Focusing on the transport sector, at present, the fossil fuels utilization is responsible of about 26% of the total greenhouse gases emissions in the European Union, with the largest impact determined by road transport (72.9%), followed by aviation (13.3%) and the naval sector (12.8%) [2]. The naval industry is also responsible for local air pollution, since it causes consistent concentration of NO x , SO 2 and PM in port areas [3]. These emissions have risen substantially since the 1970s with the growth of international trade [4]. Ships' propulsion systems are currently based on internal combustion engines fueled with diesel fuel or MDO (Marine Diesel Oil), the latter being often adopted due to its lower price compared to more refined fuels, despite the consistent content of sulfur [5]. The EU recently set stringent limitations on fuels for ships: starting from 2020, in the naval sector, the current fuels must be replaced with products with a maximum concentration of sulfur equal to 0.5% (at present it is equal to 3.5%) [6]. In this context, the conversion of oil engines for working with liquefied natural gas (LNG) can become a convenient 1.
the thermodynamic design and optimization (presented in Section 2) on a small-scale LNG production system, varying the values of the supply pressure, in order to identify the optimum layout and parameters (e.g., compressor pressure ratio, temperature and pressure in each plant section, etc.) that minimizes the total specific electric energy consumption (e); 2.
the development of a model (Section 3) which, based on the parameters obtained in the thermodynamic optimization, allows to develop the engineering design and dimensioning of each component of the plant. The obtained results (i.e., compressor volumetric efficiency and displacement, etc.) are used as starting point for the off-design analysis; 3.
the development of an off-design model (Section 3) that estimate the LNG system behavior when it is coupled with a renewable non programmable generator; In this study, this general procedure, which can be used for the design and the performance evaluation of a small scale LNG production plant, coupled with renewable generators, has been applied to a case study (Section 4) consisting in the installation of an LNG production plant at the port of Ravenna. In more detail, the optimal size of the renewable plant to be coupled with the LNG production system has been defined and the yearly operation of the integrated system has been evaluated.

Process Description
A schematic representation of the LNG production layout is presented in Figure 1. It includes the basic processes for the liquefaction of the natural gas, i.e., a compression unit, a cooling section, an expansion device, a separation tank (flash tank) and a storage tank. The most generic case shall comprehend a secondary compression line to compress the recirculated vapor coming from the flash tank for the mixing with the fresh natural gas (NG) stream. As it will be better clarified in the following, the presence of a compression in the secondary line depends on the hypotheses and parameters assumed in the analysis, and in particular, on pressure values selected for the gas grid and for the storage tank. The cooling section is composed of an ammonia compression chiller and of two regenerative heat exchangers, where the low-pressure recirculated vapor is used to cool down the stream coming from the primary compression line.
To limit the consumption of the compressors, both compression lines are cooled by inter and after-coolers, whose optimized thermodynamic design was assessed in [24]. As an expansion device, a Joule-Thompson valve and a turbo-expander are alternatively considered in the simulations. The expander mechanical power is given to the compressors, leading to a reduction of the energy externally provided to the process; on the other side, the installation of an expander increases the investment costs, thus the energy savings must be consistent to support this option.

Thermodynamic Modeling
The thermodynamic optimization was realized through a physical-mathematical model, developed in Visual Basic for Application (VBA) environment. In more detail, through an iterative procedure, the model allows to determine the optimal plant layout and the physical states (pressure, temperature and mass flow rate) of each sections in order to achieve the optimal values of the performance indicators usually considered for these power plants (as will be explained in the Section 2.3). Furthermore, by means of the iterative resolution of the developed thermodynamic model, a complete thermodynamic evaluation of the power plant is carried out by varying the key parameters affecting the performance of the plant, namely (i) the layout, (ii) the pressure at the exit of the primary compressor line (p 6 ) and (iii) the storage tank pressure (p 10 ). In particular, these last two parameters are varied between a minimum and a maximum value, and all their possible combinations are considered.
During the calculation, the developed model also defines the configuration of the secondary compressor line on the basis of the value of natural gas grid pressure. The second compressor line options are (with reference to Figure 1): (1) compressor C2-1 and inter-cooler C2 + compressor C2-2 and after-cooler C2; (2) compressor C2-1 and inter-cooler C2; (3) no secondary compressor line.
Furthermore, depending on the considered value of the gas network supply pressure and for given values of the pressures at the exit of the primary compressor line and at the storage tank, the total compression ratio for each compression train is completely defined, if the pressure losses through the various heat exchanger are assumed. The developed model opportunely splits the total compression ratio between the two compressor stages by applying the minimum compressor specific work criterion. This criterion has been modified by the Authors [24] in order to take into account also the pressure losses through the compression train.
As it regards the cooling section, the components HE3 and HE1 are simply calculated by the method of logarithmic mean temperature difference, while HE2 performance is estimated by assuming a value of the temperature at the chiller outlet (T 8 ) as input.
Finally, it is possible to choose between the use of Joule-Thomson (JT) valve or a turbo-expander (TE) in the expansion section. In this last case, the expander performance is estimated by considering a constant value of isentropic efficiency.
For further details about the mathematical models adopted in this thermodynamic optimization tool, please refer to [24].
Four scenarios have been taken into account relating to the value of natural gas supply pressure: high pressure supply network, 70 bar (HPN), -medium-high pressure supply network, 30 bar (MHPN), -medium-low pressure supply network, 15 bar (MLPN), -low pressure supply network, 3 bar (LPN) According to the most common pressure values in the various branches of Italian national gas network [25].
Furthermore, the pressure at the exit of primary compressor line (i.e., the maximum pressure of the cycle) has been assumed varying between 150 bar and 300 bar, while the storage pressure ranges between 3 bar and 15 bar. These values have been chosen considering the current technology regarding the gas compressors and the available technology for LNG storage tanks on vehicles.
The main flow chart of the calculation model, which accomplishes the thermodynamic design and optimization phase (level 1) is presented in Section 3, as part of the complete procedure for the design and performance evaluation. In addition, further details on the thermodynamic optimization procedure can be found in [17].

Performance Indicators
As mentioned in Section 2.2, the aim of the thermodynamic optimization analysis is to find the layout and the thermodynamic parameters that allow to minimize the energy consumption, while maximizing the LNG production. To this purpose, the evaluation is made by means of the total specific electric energy consumption, e, already introduced in [24]. It is expressed in kJ/kg and defined by Equation (1): where P tot [kW] is the total electric power required by the whole liquefaction process and . m LNG [kg/s] is the LNG produced mass flow rate. The amount of the consumed electric power corresponds to the power required by the two compression lines minus the power supplied by the expander, if installed: For a fair comparison between the various cases, all the simulations are referred to a production of LNG equal to 1 kg/s, considering that all the results are scalable to different plant sizes. The quality of the NG in the flash tank, (x f t = m vap m tot ), is also taken into consideration as an indicator of the process performance, since it strongly influences the total consumption of the compressors. In fact, the higher the quality in the flash tank, the higher the fraction of vapor that needs to be compressed and recirculated to the compressor C2, which rises its required power, increasing also the flow rate and consumption of the compressor C1 due to the higher flow rate coming from the mixer. From a different viewpoint, at given mass flow rate of natural gas from the grid, the mass flow rate of produced LNG decreases at higher values of the fluid quality after the expansion.

Hypothesis and Assumptions
The main assumptions, that are common to all the considered cases, are collected in Table 1. In detail, the temperature of NG coming from the grid (T 1 ) is taken equal to 20 • C, regardless of value of the supply pressure. The temperature values at the outlet of inter and after-cooler on both the compression lines are set at 30 • C. Each compressor is characterized by polytropic efficiency of 0.63 and electromechanical efficiency equal to 0.90. These values were imposed following previous evaluations on the performance of volumetric compressors for this kind of applications discussed in a previous Authors' study [24].
Regarding the cooling section, a temperature equal to −50 • C is selected at the outlet of the chiller unit (T 8 ), whose Energy Efficiency Ratio (EER) is imposed equal to 1.1. The influence of the chiller outlet temperature on the cycle performance was evaluated by the Authors in [24], highlighting that this unequivocally affects the overall consumption of the circuit, which decreases as the outlet temperature decreases. The effectiveness of both the regenerative heat exchangers is assumed equal to 0.70, while the pressure drops through the various heat exchangers (including inter and after-coolers) are assumed equal to 2% [24]. Finally, the turbo-expander, when included in the layout, works with an isentropic efficiency equal to 0.70. The natural gas composition of the stream coming from the grid is assumed as 100% CH 4 . This approximation has negligible effects on the thermodynamic performance, as already underlined in previous Authors' study [24].

Optimization Results
The results of the design parametric analysis are shown in this section. For each supply pressure scenario, the best corresponding modified layout is presented, along with the set of parameters resulting from the optimization. The main difference in the cycle layout regards the secondary compression line C2, which can vary its configuration depending on the pressure ratio that needs to realize. Both the cases with the JT valve and with the turbo-expander are considered, in order to quantify the advantage in adopting an expansion device to produce mechanical power.
The resulting layouts are reported in Figure 2a-d, relative to the supply pressure scenarios HPN, MHPN, MLPN and LPN, respectively.
The high supply pressure (Figure 2a) concerns the main pipelines of the NG distribution network. It presents the advantage of reducing the pressure ratio operated by the primary compression line, thus decreasing the corresponding compression work. At the same time, the pressure difference between the storage pressure and the supply pressure results very high, increasing pressure ratio and power required by the secondary compressor line. The latter is also more complex, involving two stages of cooling (inter and after-cooler).
The optimized layout for a supply pressure equal to 30 bar ( Figure 2b) hold the main difference with the previous case in the absence of the inter-cooler in the secondary compressor line. This is due to the lower difference between the storage pressure (p 10 ) and the supply pressure (p 1 ), that requires a lower compression work. The compressor C2 is now modelled with one unit refrigerated by an after-cooler.
The optimized layout of the medium-low supply pressure case (MLPN) provides for no compressor for the recirculation of NG vapor. Nevertheless, in Figure 2c the compression line C2 is still represented, since in the parametric analysis it was still necessary for storage pressures below 15 bar. Thus, the compressor C2 can be removed from the plant layout, determining a drop of power consumption and resulting the optimum configuration within the MLPN scenario. In addition, the costs of purchase and installation of the compression line are avoided, reducing the total investment for the liquefaction plant. On the other hand, clearly increasing the complexity and costs for the high-pressure storage facilities must be considered too.
In a low supply pressure case (LPN, Figure 2d), the secondary compression line C2 results unnecessarily within the whole parametric analysis because the pressure in the flash tank is always higher than the supply pressure. Rather, a lamination valve (V2) could be necessary to reduce the pressure in case of high storage pressure. The resulting layout is the simplest among the analyzed scenarios and this results in less complexity and number of components to be installed, less maintenance and lower investment cost.
The thermodynamic results of the sensitivity analysis are described by Figure 3a-d, that report the performance indicator e and the quality of the fluid at the flash tank (x ft ) as function of the maximum pressure (p 6 ) and of the storage pressure (p 10 ), for the four supply pressure scenarios. In the 3D plots, the JT valve and the turbo-expander cases are represented, respectively, by the upper surface and by the lower one. In Table 2 the specific values of the parameters that allow the lowest consumption for the liquefaction process are collected, i.e., the values representing the minimum points of each surface in Figure 3a-d. Finally, in Figure 4a-d, the thermodynamic cycles related to each process are traced in pressure-enthalpy diagram. The following comments can be considered valid for all the four scenarios.   The impact of the adoption of the turbo-expander is significant, since it leads to a reduction of the total specific energy consumption between 27% and 18% (the percentage of energy savings decreases with the supply pressure reduction, see Table 2). The plots highlight a further advantage of replacing the lamination valve with an expander: the first solution, indeed, involves an isenthalpic process, while the second one causes a reduction of the enthalpy value, which implies also a lower quality of the natural gas at the end of the expansion, with consequent saving in energy consumption for the liquefaction process. However, the major impact on the specific energy consumption e is given by the storage pressure p 10 : in fact, regardless of the maximum pressure p 6 , the vapor quality in the flash tank decreases considerably with the increase of the storage pressure (p 10 ), reducing the fraction of recirculated NG as vapor and so the consumption of the compressors. This is due to the shape of the saturation curve and to the consistent reduction of the heat of vaporization, as the pressure rises from 3 bar to 15 bar. Obviously, also the power required by the compression chiller rises with the increment of the NG quality. It follows that, in all the scenarios, the optimized solution involves the maximum value of storage pressure within the considered range (3-15 bar).
The influence of maximum pressure of the circuit, p 6 , appears less intuitive: being the total specific consumption a result of the sum of three different components having opposite characteristic curves, the representative equation is not a monotonic function and has a global minimum. Looking at the pressure-enthalpy diagrams in Figure 4a-d, at given temperature at the inlet of the expansion device, the consumption of the compressors is affected by the pressure ratio, which increases with p 6 , and by the mass flow rate of NG over the entire cycle, that is decreased because of the quality reduction. This occurs for values of pressure up to 220 bar, above which the advantage is no more relevant. The same trend is obtained both with the expansion valve and with the turbine.

Discussion on the Proposed LNG Production Process
As already mentioned in the introduction of this work, the proposed LNG production process has been developed by Authors and presented in previous studies [23,24,26] with the main aim to realize a plug & play solution to be placed directly at vehicles and/or ships refueling stations. Indeed, small-scale LNG production facilities are currently rare and mainly realized in order to reduce the workload on gas reservoirs sites [26] or for offshore compact LNG production [27]. In this scenario, the idea of generating LNG directly at filling stations has been developed, in order to give a compact solution and to avoid the economic and environmental costs related to the LNG transportation by ship, contemporarily promoting LNG as fuel in transport sector. In addition, with respect to the processes usually employed in small-scale applications (i.e., N 2 expander cycle and single-mixed refrigerant process), the proposed process is aimed at the realization of simpler-and, as a consequence-less expensive solutions, which can be competitive in the future global energy scenario.
The proposed layout is general and can be easily adapted to different requirements in terms of boundary conditions. In particular, the secondary line (see Figure 1) can vary on the basis of the pressures of the natural gas network and of the storage tank, allowing to install the LNG production plant in a very wide range of conditions. In detail, the secondary line may include: (i) a compression train (with one or more compression stages, eventually inter-cooled) when the LNG storage pressure is lower than the NG network pressure, (ii) an isenthalpic valve when the LNG storage pressure is higher than the NG network pressure or (iii) no additional components when the LNG storage pressure and the NG network pressure differ only for the pressure drops through the line.
As it concerns the maturity level of technology, at present, a prototype of the LNG production process developed by the Authors has been realized and experimentally tested by the Italian Company Graf S.p.a., with the final aim to commercialize the technology [17]. Finally, the layout has been specifically developed for small-scale applications, namely for liquefaction facilities with a production capacity up to 1 MTPA (Mega Tons Per Annum). Below this value of capacity, the process is scalable and the simplicity of the solution can be maintained with economic convenience. On the other hand, for medium/large scale LNG production plants, the proposed process is not yet competitive with the commonly used technologies from an economic viewpoint.
In the following section, the calculation model, specifically developed for the proposed LNG production process, will be presented. In this respect, starting from this point, the expansion device considered for the continuation of the study is the lamination valve (JT layouts). To the Authors' opinion, indeed, this choice better reflects the current industrial state-of-the-art, since no commercial models of cryogenic expander are spread in the market yet. In more detail, at present, some prototypes of cryogenic expander are under development or at testing phase, as the advantage in saving consumption of replacing the JT valve with a turbine is indisputable. However, in order to give a larger scale applicability of the proposed approach, the developed off-design calculation model refers to JT expansion layouts for LNG production.

Engineering Design and off Design Performance Evaluation
The coupling of a small scale LNG production plant with non-programmable renewable generators involves the study of the off-design behavior of the system. This behavior is clearly influenced by the design point. For this reason, the off-design calculation model represents the third level of the developed design procedure after the definition of input, parameters and boundary conditions by the thermodynamic design and optimization phase (level 1, discussed in Section 2) and by the engineering design phase (level 2).
With reference to Figure 5, the thermodynamic design and optimization phase allows the definition of the plant layout and of the cycle thermodynamic parameters (pressures and temperatures) which guarantee the minimum value of total specific electric energy consumption as a function of the available NG supply pressure, assuming a range of variation of the maximum pressure and of the storage pressure. The output of this first evaluation level become the input of the calculation model for the engineering design of the LNG plant. In addition, the design values of the produced LNG mass flow rate and of the compressor rotating speed must be provided as calculation model input. It is important to point out that the engineering design model is based on the hypothesis to only use reciprocating compressor according to the evidence presented in [17,23].
The main outputs of the engineering design calculation (see Figure 5) are: • displacement of each compressor; • design volumetric efficiency of each compressor; • design isentropic efficiency of each compressor; • design shaft power of each compressor; • design NTU of each heat exchanger; • total absorbed electrical power; • design natural gas consumption.
These last results will be used as input for the last level of calculation (off-design calculation in Figure 5) that is able to analyze the behavior of the plant as a function of the boundary conditions. In more detail, the realized model was used in this study to calculate the LNG produced mass flow rate of the plant as a function of the available electrical power from a given non-programmable renewable generator.

Engineering Design Mathematical Model
In this study, a validated semi-empirical method for volumetric reciprocating compressors coming from the literature [28,29] is employed to calculate mass flow rate, power input and temperature at the end of each compression section. This method provides an integration between physical-based models for compressors at design speed and a polynomial approach for compressor efficiency at off-design speed.
First of all, as a starting point for the off-design calculation, the applied compressor model defines the involved physical quantities at the design point (Equations (3)-(11)), which can be calculated on the basis of the optimal results obtained by the thermodynamic analysis presented in Section 2 and of the equations presented in [28]. In more detail, based on the physical conditions (i.e., temperature and pressure) and on the mass flow rate at the inlet of each compressor ( where ρ = f (p, T) is the density of the fluid at the inlet of the compressor, expressed in [kg/m 3 ]. Furthermore, based on [28], the volumetric efficiency (η V,des ) can be calculated by the following equation: where p dis [kPa] is the discharge pressure (evaluated from the imposition of the compressor pressure ratio), p suc [kPa] is the suction pressure, k is the isentropic coefficient and b 1 , b 2 and dp are model parameters [28]. As a consequence, for known compressor speed (N des , [Hz]) and starting from the definition of the volumetric efficiency (Equation (5)), the compressor displacement (V disp , [m 3 ]) can be calculated by means of the Equation (6).
Thus, the compressor shaft power (W des , [kW]) can be calculated as: where a 1 , a 2 and a 3 [kPa] are three model parameters [29]. Furthermore, in Equation (7), a power loss is accounted (by means of the coefficient η o ): as evaluated in the semi-empirical model presented in [28], indeed, a power loss equal to the 2% is registered. Once evaluated, the compressor shaft power, the isentropic efficiency (η is,des ) can be calculated by means of the following equation: where h is is the discharge enthalpy in an isentropic process [kJ/kg], h suc is the suction enthalpy [kJ/kg], v suc = 1 ρ is the suction specific volume [m 3 /kg]. All these quantities are known from the previously presented thermodynamic analysis.
The evaluation of the discharge temperature (i.e., the temperature at the end of the compression section), finally, starts from the shell temperature (T shell , [K]), defined as shown in Equation (9): where c 1 [K], c 2 [K] and c 3 [K·kPa] are three model parameters [28].
Then, the heat transfer to the ambient ( . Q amb , [kW]) is evaluated by means of the following expression: . Q amb = US shell (T shell − T amb ) (10) where US shell [kW/K] is a model parameter and T amb is the ambient temperature [K]. Finally, the discharge temperature can be calculated from the discharge enthalpy (h dis , [kJ/kg]): The compressor model parameters employed in this study come from a dataset [30] of a reciprocating compressor and are listed in Appendix A.
Furthermore, as it regards heat exchangers, Equations (12)-(15) represent the design point. As well-known, the heat exchangers behavior at the design conditions is governed by the heat exchange Equation (12) Finally, by means of Equations (14) and (15), the heat exchangers in design conditions can be completely evaluated.
where ∆T ml,des is the mean logarithmic temperature difference [K], . C min,des is the minimum (between hot and cold sides) heat capacity [kW/K] and ∆T max is the maximum (between hot and cold sides) temperature difference [K].
where . C max,des is the maximum (between hot and cold sides) heat capacity [kW/K].

Off-Design Calculation Model
Once the physical quantities in design conditions have been calculated, the compressors off-design analysis is carried out by means of Equations (16) and (19). According to [28], the normalized volumetric and isentropic efficiencies are influenced only by the normalized compressor speed and are independent from the pressure ratio. Then, the two efficiencies can be represented as a function of the normalized compressor speed through a second order polynomial function (Equations (16) and (17)): where d 1 , d 2 and d 3 , e 1 , e 2 and e 3 are model parameters [29]. Consequently, the volumetric flow rate in off-design conditions can be calculated from Equations (5) and (16) as it follows: Furthermore, from Equations (8) and (17), the compressor shaft power in off-design conditions is evaluated as: As for the design section, the compressor model parameters come from the dataset reported in [30] and are listed in Appendix A.
Finally, it should be highlighted that the calculation model of the temperature at the end of the compression section in design conditions (Equations (12)-(15)) can be directly applied also to off-design conditions without any modification, since the model relies on system energy balance taking into account the heat transfer loss to the ambient.
On the other hand, in off-design conditions, the heat exchanger behavior cannot be described by ε-NTU method, but it is still governed by fluids energy balance and heat exchange equation, adjusted to the new operating conditions (Equations (20)- (22)). Evidently, considering the off-design of a specific heat exchanger, the exchange surface remains unchanged. As a consequence, the relationship between the thermal flows can be expressed only on the basis of the ratio between the products of global heat transfer coefficients and mean logarithmic temperature differences (Equation (20)): Being: where the subscript w is referred to the fluid with the worst convective heat transfer coefficient. Note that the exponent can slightly vary based on the considered fluids, the type of flow and the system geometry.
After the calculation of the global heat exchange coefficient, in off-design conditions (U), the NTU parameter can be determined as it follows: Finally, the heat exchanger effectiveness at off-design conditions (ε ) is still defined by Equation (15).

Case Study: A LNG Refueling Station for Maritime Applications
In this section, the calculation procedure presented in the previous paragraphs is applied to a real case study, consisting in the project-funded by the Emilia Romagna region and developed by the University of Bologna-to realize a LNG refueling station at the port dock of Ravenna (Italy) for maritime applications [31]. At present, a ferry is used to carry cars from one side to the other of the main canal of the Ravenna port (400 m wide) at a maximum rate of 10 runs per hour, for a total of 19 h per day in the whole year. The ferry is equipped with two diesel engines of 162 kW each.
The developed idea entails the modification of the two engines into dual fuel engines working with LNG and diesel fuels. From the design of the new propulsion system [31], the single-engine power is supposed to decrease down to 125 kW at 1500 rpm. The ferry will be further equipped with four storage tanks of 3300 L each, operating at a pressure equal to 3 bar, refueled approximatively every week. Therefore, the LNG system must produce less than 4000 kg per week.
Starting from this idea, and accounting for these constraints, the first objective was to properly design the required LNG production plant. Furthermore, it is supposed that all the electrical energy required to move the compressors and the compression chiller is provided by renewable energy sources (in particular, wind). In this respect, the evaluation of optimal size of the wind plant, which allows to weekly produce the required amount of LNG, is also object of the presented analysis. In more detail, in the next paragraphs, the definition of the LNG production plant configuration and operating parameters-based on the optimization analysis presented in the Section 2 of this paper-and the application of the developed calculation model in order to define the optimal coupling with a wind plant will be presented and discussed. To this respect, it should be highlighted that the integration between the LNG production process and a wind plant is particularly interesting in order to promote the renewable sources exploitation without posing issues to the national electrical grid (such as frequency instability, capacity challenges, etc.). However, this has to be intended just as a possibility, while the proposed LNG production process and the design and performance methodology are valid independently from the source of the electricity which feeds the process.
Future works, after the realization of the proposed integrated plant, will analyze the real plant operation and the relative experimental results.

LNG Production Plant Configuration and Assumptions
Based on the grid pressure levels available on the installation site, the value of 3 bar was chosen as supply pressure, thus the corresponding optimized layout resulting from the application of the thermodynamic design and optimization model is the one shown in Figure 2d with the Joule-Thompson valve for the LNG production. As already explained, indeed, the choice of considering the JT valve instead of the turbo-expander has been made considering the current technology development.
For the sake of completeness, the optimal values of the thermodynamic parameters assumed to design the LNG plant and kept constant for the whole off-design analysis are reported in Table 3.
As it regards the LNG production plant, it has been sized in order to be able to produce the required weekly amount of LNG with a security factor around 2 (corresponding to about 80 equivalent hours/week of operation). Based on these considerations, the LNG design mass flow rate was set at about 47 kg/h (0.013 kg/s). Table 3. Constant thermodynamic parameters for the design and off-design analysis (plant sections with reference to Figure 2d).

Parameter Value
NG composition 100% CH 4 T 1 20 • C p 1 3 bar T 2 −5 • C p 2 3 bar p 3 25 bar p 5 175 bar EER chiller 1.1 T 8 −50 • C p 11 15 bar The assumed compressor velocity, instead, was chosen equal to 1500 RPM according with the actual state of the art for what regards the compression of natural gas in the selected range of mass flow rate.
By the application of the engineering design calculation model, in the basis of the results of the thermodynamic analysis and of the values of LNG mass flow rate and compressor rotational speed, it follows a power plant in design condition with an electrical power consumption slightly higher than 43 kW. Some of the results of the developed calculation are presented in Table 4, while an example of engineering design procedure can be found in the Appendix A of this paper. On the basis of the previous results, the application of the off-design model allows the drawn the curves in Figures 6 and 7.
In more detail, Figure 6 shows the electric energy consumption of the compressors and of the chiller separately, as a function of the total electrical load of the system (i.e., the electric power required by the whole liquefaction process). As it can be noted, starting from the design point, the calculation model shows upper and lower operating limits for off-design conditions. In particular, the operating range can vary from 40% to 200% of the design electric load (this range derives from the operating range of the compressors). As it can be noted, the upper limit seems to be a quite high value and probably in real operation it cannot be reached, but it is presented here in order to show the validity range of the model. Furthermore, in Figure 7 the LNG produced mass flow rate are shown as a function of the plant load: an increase in this quantity with the electric load increase can be observed (with a variation relating to the design point from a minimum of −53% to a maximum of about +62%), since the electric energy available to the compressors is higher.

Wind and LNG Plants Integration
Once the configuration and the size of the LNG plant and the off-design behavior are known, the objective was the determination of the optimal size of a wind plant to be integrated into the system in order to provide the required electricity for the liquefaction process. Since, for the specific application, a certain amount of LNG must be guaranteed every week, the optimal wind size has been determined accounting for the plant operation during a whole year.
As a consequence, with this purpose, the production data of existing wind farms and the ones available from TERNA [32], related to the whole Italian territory, have been considered. Thus, starting from these data, the normalized wind production profile (on an hourly basis) has been obtained ( Figure 8). Then, a dedicated code, developed in VBA, has been developed in order to consider the coupling between wind and LNG plants: the objective was the determination of the wind size, which allows contemporarily to respect the constraint of LNG production (around 4000 kg/week every week per year) and to minimize the electricity introduction into the national electric grid. The developed code, indeed, allows to evaluate the yearly operating time of the LNG production process, the surplus of electric energy not used by the plant-and then introduced into the electric grid-and the total specific electric energy consumption as a function of the wind plant size. In more detail, the code-for each considered wind plant size-analyzes the coupling of the wind and LNG plants for a whole year, hourly based: -if, in the considered hour, the electrical production of the wind plant is lower than the lower operating limit of the LNG production process, then all the electric energy is introduced into the electric grid; -if, in the considered hour, the electrical production of the wind plant is within the operating range of the LNG plant, then all the electric energy is employed in the process; -if, in the considered hour, the electrical production of the wind plant is higher than the upper operating limit of the LNG production process, then the difference between the wind electric energy and the maximum electric load of the LNG plant is introduced into the electric grid.
The results of the coupling LNG and wind plants analysis are shown in Figure 9, in terms of the yearly LNG production (red curve) and electricity surplus to be sold to the grid (green curve), as a function of the wind size. As it can be seen, the wind size which allows to produce the LNG amount yearly required (about 4000 kg every week for a total of about 208,000 kg/year) is equal to about 78 kW. In this case, the electric energy surplus to be yearly sold to the grid results equal to a value slightly lower than 33 MWh/year. For wind sizes lower than 78 kW, the mass of LNG produced does not fit with the ferry needs constraint and, consequently, should not be considered for this case study. On the other hand, for higher wind sizes, an excess of LNG production is seen with respect to the case study constraint. However, it should be highlighted that, for wind sizes higher than the one selected as the optimal for the presented case study, a further decrease in electricity excess to be sold to the network is registered. This behavior can be explained considering the upper limit of the operational range of the compressors. Finally, in order to evaluate the LNG plant hours of operation, the wind production during a year is shown in Figure 10. As it can be seen, the LNG plant-during a whole year of operation-works for 187 h/year in design conditions and for 4720 h/year in off-design, split between 1455 h/year in which the electrical input is higher than the rated electrical power and 3265 h/year in which the electrical input is lower. On the other hand, for 3853 h/year the LNG plant is shut-down because the wind production is below the lower operating limit of the LNG plant. On the other hand, as a confirmation of the results shown in Figure 9 for the electricity yearly sold to the network, the upper operating limit is not reached with the selected wind size.

Concluding Remarks
The liquefied natural gas (LNG) can be adopted as an alternative to traditional fuels in the transport sector, to reduce the related costs and pollutant emissions. The adoption of small-scale liquefaction systems directly installed at the filling station is a solution for avoiding transportation costs, which have large impact on the whole LNG value chain. In this study, a numerical investigation on small-scale LNG production plants is performed, with the aim of assessing the best conditions that minimize the total consumption of the process, both in design and off-design behavior. A three-levels calculation code has been developed, in order to (i) optimize the layout and thermodynamic parameters, (ii) develop the engineering design and dimensioning of each component of the plant and (iii) analyze the plant operation in off-design conditions (represented in this study by the variation of the electric power input, due to the coupling with renewable non-programmable sources). In more detail, at the first level of the proposed procedure, a thermodynamic analysis is conducted for selecting the optimal layout and design parameters (optimal pressure and temperatures in each section of the plant), for various scenarios of the supply pressure available at the plant site (between 3 bar and 70 bar) and varying the maximum and the storage pressures. It was obtained that the configuration that assures the lowest total specific consumption is the one with 70 bars of supply pressure, but it must be considered that the related layout involves more components and complexity. The best storage pressure has proven to be 15 bar for all the supply pressure levels. Then, the second level of the procedure allows to determine the engineering design of the plant, including the specifics of the various components (compressors, inter/after-cooler and other heat exchangers), starting from the layout and parameters optimized in the previous procedure level. Finally, at the third level, the off-design analysis is carried out, allowing to evaluate the produced LNG and the plant operation conditions depending on the electric power input. The mathematical model for the engineering design and for the off-design evaluation of compressors is based on a semi-empirical approach, integrating physical models and polynomial correlations. As it regards the heat exchangers, instead, the corresponding performance are evaluated by the ε-NTU method in design conditions and by its modification for off-design operation.
The approach is therefore applied to a specific case study, regarding the project in Ravenna (Italy) of installing a plug & play natural gas liquefaction plant at the port dock, for producing the demand for fuel of a ferry that works continuously in the port canal. In this application, a wind plant is chosen to provide the required electricity for the LNG process. The aim of this last section was to prove the effectiveness of the off-design model for assessing the optimal liquefaction process layout and working conditions and the wind generator size, starting from the production of LNG required by the ferry operation. In any case, the developed procedure is general and can be extended to similar case studies, in which an integration between a non-programmable renewable generator and a LNG production plant occurs.  In the following of this Appendix A, as an example of the application of the proposed calculation model to the presented case study of Section 4, the design evaluation of the first compressor (C1-1 in Figure 2d) is presented in detail with reference to Equations (3)- (8).
Considering a design mass flow rate and a density of the fluid at the inlet of each compressor respectively equal to 0.031 kg/s (necessary to produce 0.013 kg/s of LNG) and to 2.175 kg/m 3 (fluid density at 3 bar and −5 • C), the suction volumetric flow rate is: being p dis = 24.5 bar and p suc = 3 bar as a result of the application of the thermodynamic optimization analysis presented in Section 2. Furthermore, a value of the isentropic coefficient equal to k = 1.4 has been set.
As it relates the compressor speed, the compressor of the analyzed case study is characterized by a value equal to 1500 rpm, corresponding to N des =25 Hz. As a consequence the displacement results equal to: where h is = 1213 kJ/kg and h suc = 842 kJ/kg are known from the previously presented thermodynamic analysis.