Next Article in Journal
Optimal Sliding-Mode Control of Semi-Bridgeless Boost Converters Considering Power Factor Corrections
Previous Article in Journal
Progress in Research on Net-Zero-Carbon Cities: A Literature Review and Knowledge Framework
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Dual Fuel ICE-Based Micro-CHP System and Experimental Evaluation of Its Performance at Light Loads Using Natural Gas as Primary Fuel

by
Edisson S. Castaño Mesa
*,
Sebastián H. Quintana
and
Iván D. Bedoya
Grupo de Ciencia y Tecnología del Gas y Uso Racional de la Energía—GASURE, Departamento de Ingeniería Mecánica, Facultad de Ingeniería, Universidad de Antioquia, Calle 70 No. 53-108, Medellín 050010, Colombia
*
Author to whom correspondence should be addressed.
Energies 2023, 16(17), 6281; https://doi.org/10.3390/en16176281
Submission received: 20 April 2023 / Revised: 27 May 2023 / Accepted: 1 June 2023 / Published: 29 August 2023
(This article belongs to the Section I2: Energy and Combustion Science)

Abstract

:
This study presents the implementation of a micro-generation system and its operation procedure, based on a dual fuel diesel engine using natural gas as the primary fuel and conventional diesel as the pilot fuel. On the other hand, the evaluation and validation results by experimental testing of a model according to International Energy Agency—IEA—Annex 42, applied to dual fuel diesel micro-cogeneration system, are also presented. The control procedure for experimental operation depends of both inputs’ electric power generation demand and desired substitution level due a given natural gas availability. The heat recovery system of the micro-generation system uses a gas–liquid compact heat exchanger that was selected and implemented, where wasted heat from exhaust gases was transferred to liquid water as a cool fluid. Effective operation engine performance was determined by measurement of masses’ flow rate such as inlet air, diesel and natural gas, and also operation parameters such as electric power generation and recovered thermal power were measured. Electric power was generated by using an electric generator and then dissipated as heat by using an electric resistors bank with a dissipation capacity of 18 kW . Natural gas fuel was supplied and measured by using a sonic nozzle flowmeter; in addition, natural gas composition was close to 84.7 %   CH 4 , 0.74 %   CO 2 and 1.58 %   N 2 , with the rest of them as higher hydrocarbons. The highest overall efficiency (electric efficiency plus heat recovery efficiency) was 52.3 % at 14.4 kW of electric power generation and 0 % of substitution level. Several substitution levels were tested at each engine electric power generation, obtaining the maximum substitution level of 61 % at 17.7 kW of electric power generation. Finally, model prediction results were closed to experimental results, both stationary and transient. The maximum error presented was close to 20 % associated to thermal efficiency. However, errors for all other variables were lower than 10 % for most of micro-cogeneration system operation points.

1. Introduction

Compression ignition (CI) engines operating in dual-fuel mode are a reliable technological strategy for power generation, reduction of heavy fossil fuel consumption, and attenuation of carbon dioxide and soot emissions. The so-called “dual fuel engines” have been extensively tested around the world to quantify the mentioned benefits when compared with typical diesel fuel operation, mostly at constant load frameworks, leading to a well-established international knowledge on performance and emission characteristics of dual fuel engines for power generation. Dual fuel mode means that two fuels are used for heat release during combustion, and although many options could be considered for this purpose, usually a gaseous fuel named “main fuel” [1] provides around 60 % to 97 % of the heat released and diesel acts as “pilot fuel” or ignition source for the air-gaseous fuel mixture and provides the rest of the heating energy [2]. Diesel substitution is one of the most important parameters to evaluate the dual fuel performance, and a bigger share of gaseous fuel promotes premixed combustion in the air gaseous fuel mixture, leading to reductions in soot emissions [3,4]. Depending on how the main fuel is admitted, dual fuel engines are classified as “port-fuel injected” if this fuel is admitted along with air in the intake system, and “direct-fuel injected” if the main fuel is admitted into the cylinder when the piston is close to the top dead-center position [5]. In dual fuel engines, the intake system is usually designed to admit a premixed air gaseous fuel mixture; however, direct injection of the gaseous fuel during the compression stroke has been used in some research projects and commercial engines [6,7].
Natural gas as gaseous fuel is used in many applications such as furnaces, water heaters, cooking stoves, processing foods, air conditioning, generating steam in water boilers then used in steam turbines, power generation in gas turbines and internal combustion engines both stationary and vehicular [8]. The composition of natural gas depends generally on its source of extraction, thus its heating value usually varies from 26.000 kJ/kg (700 Btu/scf) to 59.000 kJ/kg ( 1.600 Btu/scf) [9]. On the other side, natural gas is mostly composed of >85% methane ( CH 4 ), 3–8% ethane ( C 2 H 6 ), <1% propane ( C 3 H 8 ), <2% heaver hydrocarbons as butane ( C 4 H 10 ) and pentane ( C 5 H 12 ), 1–2% carbon dioxide ( CO 2 ), 1–5% nitrogen ( N 2 ) and with minor constituents such as helium ( H 2 ) and hydrogen sulfide ( H 2 S ) [8,10]. Natural gas is considered a cleaner fuel compared with other heavier fossil fuels, thereby offering important environmental benefits because of a significantly lower production of particulate matter ( PM ), oxides of nitrogen ( NOx ) and hydrocarbons ( HC ) than conventional diesel and gasoline combustion engines [11,12,13].
Aiming to enhance the dual fuel diesel engine performance keeping high substitution levels, some strategies to improve the auxiliary components of the engine may be implemented such as a direct common-rail injection system, injection delay setting depending on load and speed, and even turbochargers [14]. However, dual fuel diesel engines could be seen as double power generation source systems, combining electric or mechanical power generation with thermal power generation, where this thermal power can be recovered from the water coolant and exhaust gases, allowing them to replace the installation of electric heaters and gas burner plants when combined application is needed [15,16]. Some processes can need heating, such as textile products drying, dehydration of vegetables or fruits, water and air heating or substances. The engine exhaust gases can achieve temperatures close to 150∼450 °C depending on the load and substitution level in dual fuel engines, having a high sensible heat recovery potential [17,18]. Thereby, the strategy to generate combined electric and thermal power from recovering the exhaust gases or water coolant heat can lead to increase the thermodynamic efficiency of the system, bringing a lower overall fuel consumption, where combined production of heat and electric power refers commonly to cogeneration [19].
Cogeneration usually refers to systems able to produce simultaneously both mechanical or electric energy and useful thermal energy from the same primary energy source, where the thermal energy is recovered from wasted heat from the combustion process of engines or turbines, even in some cases from high temperature fuel cells [20,21]. If thermal energy demand is taken into consideration, installed on-site electric power generation has an important thermodynamic efficiency advantage compared to centralized power generation, as wasted heat energy from the electric power generation process can be used [22]. Often, the centralized power generation plants do not use these wasted heat energies and these same energies must be ejected into the environment, representing economic loss coupled with greater specific emissions releasing. When an on-site power generation plant recovers wasted heat energy for using on the process itself, then the thermodynamic efficiency rises, becoming a fuel-economic advantage, and this thermodynamic efficiency is also named overall thermal efficiency, which relates electric power and thermal power to chemical fuel energy based on lower heating values [19,22].
In cogeneration, the overall thermal efficiency of energy conversion increases depending on the prime mover cycle, and internal combustion engines and Stirling engines have an electric efficiency close to 20–35% but obtain an overall thermal efficiency over 45% up to 80%, micro steam turbines have an electric efficiency close to 10–20% but obtaining an overall thermal efficiency over 65% up to 80%, micro gas turbines have an electric efficiency close to 15–30% but obtain an overall thermal efficiency over 60% up to 80% [23], and finally, fuel cells have an electric efficiency close to 30–70% or 25–40% but obtain an overall thermal efficiency over 60% up to 80% depending on the technology of the fuel cell, SOFC (Solid Oxide Fuel Cell) or PEMFC (Polymer Electrolyte Membrane Fuel Cell), respectively [24].
Models have been developed for applying and control of cogeneration systems depending on operation parameters such as electric efficiency, thermal efficiency associated to heat recovered, mass flow of heat recovery fluid, temperatures of heat recovery fluid, properties and mass flow rate of fuel and working fluid, and thermal and geometric parameters of both the heat recovery device and electric power generation system [25]. These models have to be as simple as possible, but also able to predict the real operation of the micro-cogeneration system. These models can be experimentally calibrated to try to obtain the operation maps and variables prediction as described above [26,27].
Therefore, this study arises as an experimental validation of a previous study, which develops and calibrates a dual fuel engine micro-cogeneration (DECOG) model having as the main input parameters the electric power generation and substitution level. The present study shows the implementation and experimental evaluation of a micro-cogeneration system and also its model validation for a dual fuel diesel-engine-based micro-cogeneration system (electric power close to 15 kW) using natural gas as the primary fuel and diesel as the pilot fuel, where the thermal power is recovered only from wasted heat from exhaust gases. The applied models for this study are dependent on electric power generation and substitution level as the main operating parameters in a dual fuel diesel engine. Experimental operation maps allow us to calibrate the models for subsequent evaluation and validation. Additionally, a transient validation and prediction capacity is presented based on the outlet heat recovery fluid temperature as the main control parameter for micro-cogeneration system implementation in real applications. Additionally, a description of the operation control procedure is explained, which depends on both inputs of electric power generation demand and substitution level for a given natural gas availability, then sensor and actuators are set electronically by a computer achieving as closed as possible the desired inputs and measuring output variables as recovered heat.

2. Models Theoretical Description

A model for an internal combustion engines-based cogeneration system has been developed initially within International Energy Agency (IEA) Annex 42 for usages applied to building simulation programs, but then modified and adapted by DECOG for a dual fuel diesel-engine-based micro-cogeneration system [28]. The model aims to accurately predict the electric and thermal power outputs of any internal-combustion-engines-based cogeneration system. The correct operation of the model requires a specified number of coefficients found experimentally, so that this model must be calibrated and validated by experimental results, which will have applicability only to the implemented engine type, implemented heat recovery devices and fuel for the testing operation [29]. The model defines the control volume for the cogeneration device as Figure 1 shows and it is explained below [25].
  • Energy conversion control volume (ECCV) involves the engine working fluid, exhaust gases and engine alternator, providing information from the combustion engine unit operation map to the model.
  • Combustion engine control volume (CECV) involves the aggregated thermal capacitance associated with the engine block and heat recovery device shells.
  • Heat recovery water control volume (HRWCV) involves the recovery water flowing through the devices and heat recovery devices in thermal contact.

2.1. Dual Fuel Diesel Engines

In dual fuel diesel engines operation, the gaseous fuel as the primary fuel enters together with intake air; however, it does not auto-ignite by itself. Ignition of that premixed gas is achieved when a small diesel fuel quantity is injected as a pilot fuel, timed at the end of the compression stroke [18]. The pilot fuel auto-ignites first, inducing mixture with the primary fuel, and combustion air ignites as well. Later on, gaseous fuel combustion occurs via a flame initiation at an unspecified location in the combustion chamber [30,31]. Thus, depending on the primary fuel quantity, a pilot fuel quantity must to be injected, which is why a substitution level Z is defined for dual fuel operation as the percentage between the quantity of diesel fuel reduction during dual mode and conventional diesel fuel quantity at a given power [5,32], as shown in Equation (1).
Z = m ˙ D , c m ˙ D m ˙ D , c · 100
where m ˙ D , c is the diesel mass flow rate in conventional operation, i.e., at 0 % substitution level, m ˙ D is the diesel mass flow rate in dual fuel operation mode. Due to dual fuel operation, equivalence ratio ϕ can be redefined as a function of both fuels, natural gas and diesel. Equation (2) shows the equivalence ratio for dual fuel operation mode [30].
ϕ = m ˙ D ϕ D , s t q + m ˙ N G ϕ N G , s t q m ˙ a
where ϕ D , s t q and ϕ N G , s t q are the stoichiometric fuel–air ratio for diesel and NG, respectively, and m ˙ N G and m ˙ a are the NG and combustion air mass flow rates, respectively.

2.2. Energy Conversion Control Volume

A steady state energy balance can be applied to the energy conversion control volume to obtain Equation (3). Such equation relates the air and fuel total enthalpy as equal to the sum between electric power generation, heat power production and exhaust gases total enthalpy.
H ˙ a + H ˙ f = P E + P T + H ˙ e x h
where H ˙ a , H ˙ f and H ˙ e x h are the air, fuel and exhaust gases total enthalpies, respectively, and P E and P T are the electric power generation and thermal power recovered, respectively. However, the model does not aim to complete characterization of the energy balance described by Equation (3). Instead, the engine performance is related to total chemical power P Q added by both gaseous and pilot fuels, so that natural gas and diesel are included as shown in Equations (4) and (5).
P E = η E P Q = η E m ˙ D L H V D + m ˙ N G L H V N G
P T = η T P Q = η T m ˙ D L H V D + m ˙ N G L H V N G
where η E and η T are the electric and thermal efficiencies, L H V i is the lower heating value of fuel species i, m ˙ i is the mass flow rate of fuel i, and subscripts D and N G represent the diesel and natural gas, respectively. In steady state operation, the electric and thermal efficiencies can be determined by using empirical correlations, which could relate variables as electric power, heat recovery water flow rate and temperature of the same heat recovery water [29]. Nevertheless, electric and thermal efficiency correlations can be redefined as a function of electric power and substitution level, since the current study applies a dual fuel diesel engine for the micro-cogeneration system, keeping as a constant the heat recovery water flow rate and also its inlet temperature as much as possible. Consequently, Equations (6) and (7) show the electric and thermal efficiency correlations, respectively. The heat-electric power ratio is calculated in Equation (8).
η E = a 1 + a 2 P E + a 3 P E 2 + a 4 Z + a 5 Z 2 + a 6 P E Z + a 7 P E Z 2 + a 8 P E 2 Z + a 9 P E 2 Z 2
η T = b 1 + b 2 P E + b 3 P E 2 + b 4 Z + b 5 Z 2 + b 6 P E Z + b 7 P E Z 2 + b 8 P E 2 Z + b 9 P E 2 Z 2
R H E = P T / P E = η T / η E
Additionally, for diesel, natural gas and air mass flow rate correlations are functions of the electric power and substitution level as shown in Equations (9)–(11), respectively.
m ˙ D = c 1 + c 2 P E + c 3 P E 2 + c 4 Z + c 5 Z 2 + c 6 P E Z + c 7 P E Z 2 + c 8 P E 2 Z + c 9 P E 2 Z 2
m ˙ G N = d 1 + d 2 P E + d 3 P E 2 + d 4 Z + d 5 Z 2 + d 6 P E Z + d 7 P E Z 2 + d 8 P E 2 Z + d 9 P E 2 Z 2
m ˙ a = e 1 + e 2 P E + e 3 P E 2 + e 4 Z + e 5 Z 2 + e 6 P E Z + e 7 P E Z 2 + e 8 P E 2 Z + e 9 P E 2 Z 2

2.3. Heat Transfer Considerations

Heat can be transferred from engine coolant and exhaust gases to the heat recovery water. However, in the present study, the heat is transferred only from exhaust gases to heat recovery water, so that the model assumes the heat transfer is proportional to the difference in temperature between the heat recovery water control volume and the combustion engine control volume, which in this case is replaced by exhaust gas temperature. Then, heat transferred is quantified by using an overall heat transfer coefficient as Equation (12) shows [25].
Q H X = U A H X T e x h T w , o
where U A H X is the overall heat transfer coefficient between HRWCV and CECV, T e x h is the exhaust gases temperature and T w , o is the heat recovery water outlet temperature. It assumes that the lost heat from the engine unit is proportional to the temperature between CECV and the surroundings, thus lost heat is obtained by using another overall heat transfer coefficient, as Equation (13) shows.
Q L = U A L T c e T s u r r
where T c e is the combustion engine control volume temperature and T s u r r is the surroundings temperature.

2.4. Combustion Engine Control Volume

The dynamic thermal behavior of a micro-cogeneration system is determined by the engine block thermal mass, encapsulated working fluid and internal heat transfer. A detailed model demands to characterize the individual thermal response of each component above listed, but this model allows us to represent those as a single control volume. Hence, the thermal energy stored within CECV is quantified by using an overall thermal capacitance M C c e based on a general energy balance via Equation (14) [33,34].
M C c e d T c e d t = U A H X T w , o T e x h + U A L T s u r r T c e + P T

2.5. Heat Recovery Water Control Volume

Additionally, the model assumes that elements of HRWCV can be represented by using a single overall thermal capacitance M C w that quantifies the thermal energy stored within based on an energy balance via Equation (15).
M C w d T w , o d t = m ˙ c p w T w , i T w , o + U A H X T e x h T w , o
where m ˙ c p w is the thermal capacitance for heat recovery water mass flow rate and T w , i is the heat recovery water inlet temperature.

3. Experimental Setup

A four cylinder, four stroke, water-cooled, naturally aspired, mechanical direct injection diesel engine was operated as a dual fuel engine with natural gas as primary fuel. The technical specifications of the experimental test engine are given in Table 1. The generation unit was linked to an electric resistor load bank with a capacity of 18 kW, an electrical network analyzer type EBCHQ 54115/PD76-24E-WFF was used to measure electric power with a three-phases power controller type NOVUS PCWE-3P-100.
Natural gas as primary fuel had the properties as shown in Table 2. Additionally, commercial diesel was used as a pilot fuel with a lower heating value ( L H V D ) of 43,200 kJ/kg and a stoichiometric fuel–air ratio ( ϕ D , s t q ) of 0.064935 kgDiesel / kgAir . The natural gas mass flow rate was controlled and measured by a sonic flow meter reaching choked flow conditions by flowing through convergent nozzles, ensuring supply conditions such as pressure and temperature. The supply pressure was measured by a gauge pressure electronic transducer type Novus NP-430D with range 0–10 bar and temperature by thermocouple type K. The diesel mass flow rate was measured by using a gravimetric method with a high precision electronic weighing scale type Shimadzu TX3202L. Air mass flow rate was estimated measuring the intake pressure by using an absolute pressure electronic transducer type Wika A-10 with range 0–6 bar, and intake temperature via thermocouple type K.
For the heat recovery water loop, the water mass flow rate was measured by a hall effect water flow meter sensor type YF-S201 with operation range between 1–30 L/min. The heat exchanger implemented to recover heat from exhaust gases was a compact-type build in aluminum, one side with elliptical tubes and the second side with thin fins. Additional temperatures were measured by thermocouples type K, such as heat recovery water inlet and outlet temperature, exhaust gas temperature, engine water coolant temperature and engine oil temperature. Engine block temperature was measured by using an infra-red thermometer type Fluke 568. A detailed schematic diagram of the micro-cogeneration device setup is shown in Figure 2.
The schematic diagram shows a secondary heat exchanger, which dissipates the heat recovered by the primary heat exchanger, allowing it to keep the heat recovery water inlet temperature as constant between 28–33 °C, thus eliminating such temperatures from model equations. At the same time, the heat recovery water flow rate was kept constant and close to 25 L/min ( 1.53 m 3 / h or 425 g / s in mass). It must be clear that heat is recovered only from wasted heat from exhaust gases without heat from the engine cooling system. Finally, the micro-cogeneration system was tested with different electric power generations, such as 2.5 , 7.5 , 10.6 , 14.4 and 17.7 kW, at several substitution levels. Electric power generations of 9 kW, 12.5 kW, 14.4 kW and 16 kW were then tested with random substitution levels to determine the prediction accuracy of the micro-cogeneration system.

4. Operation Methodology

4.1. Combustion Engine Requirements

According to electric power demand for a continuous operation (primer power), i.e., operation without intermittent under unlimited number of hours, an appropriated electric generator for medium and heavy duty has been selected. These medium and heavy duty electric generators operate normally at 1800 rpm (60 Hz) for a local electric network where dual fuel ICE-based micro-CHP systems would be installed.
Based on the previous, the engine was couple to the electric generator, being designed to operate at the same engine speed so that maximum electric power demand can be supplied. The maximum demanded electric power was set at 30 kW, then 10 % was set as over electric power generation, which results in a total of 33 kW. Additionally, knowing that a dual fuel ICE-based micro-CHP system would be installed in a city 1500 m above sea level having an atmospheric pressure of 84.7 kPa, the altitude power losses must be estimated to achieve the maximum demanded electric power, including electric generator inefficiency.
Rated power and partial load power were corrected according to standard SAE J1349 [35,36]. This method provides equations to calculate the deviations in inlet air and fuel supply conditions. In this study, only inlet air deviations were considered (altitude variation) since there are no density or viscosity values for the local diesel fuel. Corrections were made by Equations (16)–(18).
Q = 120,000 m ˙ d m V d N
f m = 0.2 , if   Q / R < 37.2 0.036 ( Q / R ) 1.14 , if   37.2 < Q / R < 65 1.2 , if   Q / R > 65
P b c = P b r F a = P b r P a t m 99 α 298 T a + 273 β f m
Based on previous equations and looking for the commercial stationary engine brands at the installation city, Cummins and John Deere were the best options according to factors such as quality, reliability, price and availability of spare parts. Several available options were reviewed, finding the best option to John Deere 4039D, to which the SAE J1349 was applied, obtaining the results as Table 3 shows.
Due to altitude and electric generator losses, maximum electric power is 34.8 kW as shown in Table 3, which is barely higher than the maximum demanded electric power of 33 kW (even with 10 % of over electric power). This engine has a built-in automatically controlled diesel fuel system by centrifugal weights that allows us to keep the engine speed at 1800 rpm.

4.2. Conventional Operation as Reference for Dual Operation

At first, it was needed to measure a conventional diesel operation map, which means Z = 0 , measuring variables such as fuel consumption, exhaust gases flow and temperatures. By knowing the conventional fuel consumption at several electric powers, then a natural gas flow rate can be initially estimated to operate the micro-cogeneration system at a desired substitution level and fixed electric power. Fuel consumption, when operated with conventional diesel, is shown in Figure 3. Inputs such as electric power and substitution level are entered into a computer interface that controls electronic actuators and read sensors for several variables.
Thereby, a equation for the tendency fitted curve can be obtained for diesel fuel mass flow rate at conventional diesel operation as a function of electric power generation, as Equation (19) shows. When natural gas is introduced to the intake system, a portion of the mass flow rate of diesel and its energy has to be substituted, but instead natural gas has to add approximately the same amount of energy what was substituted from diesel, if it is assumed that electric efficiency is kept as constant. The previous statement allows us to obtain Equation (20), which combines Equation (1) as well.
m ˙ D , c = 1.657 P E 2 + 13.210 P E + 759.68
m ˙ N G = m ˙ D , c Z L H V D L H V N G
Equation (20) tends to estimate a natural gas mass flow as a function of electric power generation and substitution level. As said above, it considers no electric power inefficiency changes due to variations of substitution level increasing, that is why, during micro-cogeneration system operation, both the desired substitution level and the actual substitution level can be seen, where the last one is calculated from the actual natural gas mass flow measured in real time, allowing it to be closer to desired value.

4.3. Heat Exchanger Selection Criteria

Considerations for heat exchanger selection have been taken such as a lower pumping work in the exhaust gases side to reduce back pressure on the engine, a large heat transfer area but with a short heat exchanger volume, a low cogeneration system installation area, existing design on the local market and low cost related to manufacturing and maintenance.
Based on previous considerations, a compact heat exchanger type cross flow tube-fin with both fluids unmixed seemed the best selection due to its large free-flow area for the exhaust gases passing through fins, high ratio between total heat transfer area and heat exchanger volume, generally given for these arrangements of heat exchangers with ratios higher than 700 m 2 / m 3 ), being highly commercial in the local market for automotive industry which means a low cost related to manufacturing, installation and maintenance as well. Therefore, selection criteria values for heat exchanger were taken from initial experimental data for conventional diesel operation, trying to obtain the maximum exhaust gases flow rate and its respective temperature, also by setting a desired heat recovery water flow rate and corresponding inlet and outlet temperatures as is shown in Table 4.
Then, using values in Table 4, initial fluid properties can be estimated for both heat recovery water and exhaust gases, which then allows us to apply an iterative method to obtain dimensionless characteristic numbers and convective heat transfer coefficients for both liquid and gas sides [37]. Subsequently, by having previous heat transfer coefficients and heat transfer conductance coefficients for material, a total heat transfer coefficient can be obtained for calculating a Number of Transfer Units symbolized N T U . The method proposed by Kays and Crawford, for an unmixed cross flow heat exchanger [38], is shown in Equation (21), also methods to determine values such as effectiveness ε and pressure losses Δ p can be found.
ε = 1 exp N T U 0.22 c exp c · N T U 0.78 1
Finally, based on the above-described methodology, the general characteristics of the selected compact heat exchanger are shown in Table 5. Additionally, estimated values for heat recovery system operation are shown in Table 5, such as thermal recovered power, heat exchanger effectiveness and even pressure losses for both sides.
As shown in Table 5, a high ratio between total heat transfer area and heat exchanger volume was achieved close to 1474 m2/m3, more than twice the common values. Additionally, an estimation of actual possible heat recovery of 22.2 kW, offers a high potential for combined power generation, both electric power as much as heat power, having low pressure drop for tubes sides and low back pressure is generated by heat exchanger to the engine (17.38 kPa), allowing us to conserve as high as possible electric power generation efficiency.

4.4. Gaseous and Pilot Fuel Supply Control

The natural gas fuel flow is controlled by a sonic flow meter, which operates as a function of supply pressure and temperature, fixing values for nozzle diameters and discharge coefficients. Knowing that supply fuel temperature was kept constant, the natural gas mass flow rate depends only on supply pressure [39]. The previous procedure can be inverted, so that supply pressure can be obtained by a desired natural gas flow rate, then by combining Equation (20) on actual choked flow rate equations, Equation (22) is obtained as a result.
p 0 = m ˙ D , c Z 1 A * C d L H V D L H V N G R N G T 0 γ γ + 1 2 γ + 1 2 γ 1
where A * is the nozzle critical area and C d its discharge coefficient. Thereby, the natural gas flow rate can be controlled depending only on electric power and substitution level. For this study, supply absolute pressure was fixed at a maximum of 7.3 bar, for reducing possibilities of NG leakages and safety work. Minimum supply absolute pressure was fixed at 1.6 bar, to warranty chocked (sonic) flow conditions. In parallel, four nozzles were installed and computationally controlled. Table 6 shows the substitution level operation range for each nozzle independently of electric power generation, then Figure 4 shows the operation map of the sonic nozzle flow meter for 0.8890 mm ( 0.035 in) throat diameter for a substitution level range from 25 % up to 55 % .
In both conventional and dual operation modes, the diesel fuel flow is controlled automatically by an automatic mechanical system of weights built inside the engine that allows us to keep the engine speed at 1800 rpm.

4.5. Control and Visualization Interface

Devices and equipment explained previously were linked by using a National Instruments data acquisition and control system. A real-time operation was allowed to modify operation parameters such as orifice selection for natural gas fuel flow, safety assurance of sonic flow meter, expected substitution level Z, visualization of working parameters of engine and heat recovery systems such as temperatures, diesel fuel flow, water flow, power and efficiency at operation, and also data recording. A screenshot of the Labview interface is shown at Figure 5.

5. Results and Discussion

The micro-cogeneration system operation maps were obtained in dual fuel mode at several substitution levels, keeping the water heat recovery mass flow rate and inlet temperature as constant, at steady state, thus, the models for diesel, natural gas and air mass flow rate, electric and thermal efficiency are valid at the same state, i.e., steady operation for a long time. Nonetheless, other results are presented when running at transient operation for a desired change of thermal or electric output power or even heat recovery water outlet temperature. Previous studies have calibrated the models for micro-cogeneration systems with second- or third-degree polynomials as a function of electric power depending on the desired fit with experimental data [34,40]. However, those models presented in previous studies do not include the substitution level as a factor, since they only operate with one fuel. The coefficients obtained for calibrated models regarding the electric and thermal efficiencies as well as diesel, natural gas and air mass flow rate are shown in Table 7.
As shown in Table 7, all correlations are functions of the electric power and substitution level. It is noted that no coefficients for a overall efficiency correlation are presented, but it must be quantified as the arithmetic sum of the two efficiencies. A similar procedure occurs for the electric and thermal power ratio as shown in Equation (8).

5.1. Steady State Model Results

The operation maps were obtained by testing the micro-cogeneration system to several electric powers, however, operation data for electric power of 9 kW, 12.5 kW, 14.4 kW and 16 kW were excluded out of model calibration, aiming to prove the prediction accuracy with them. Figure 6, Figure 7 and Figure 8 show the results obtained by models (continuous lines) and experimental data (markers) for diesel, natural gas and air mass flow rates in steady state for micro-cogeneration system operation, respectively.
Figure 6 shows that the air mass flow rate decreases when the substitution level increases, which is because the substitution level increases the natural gas flow rate displacing a volume portion of the engine intake air for combustion. At a constant substitution level, the air mass flow rate continues to decrease due to a higher electric power produces an intake temperature increasing [31,41]. Other important behavior occurs when the substitution level moves between 30 % to 50 % , holding the air mass flow rate as constant and close to 47.3 g/s. The result for natural gas mass flow rate is shown in Figure 7.
Figure 7 shows that the natural gas flow rate increases when the substitution level grows, because natural gas substitutes an important portion of the chemical energy supplied by diesel. However, an opposite effect occurs when the substitution level remains constant but electric power varies, having higher natural gas flow rates for lower electric powers, which can be associated with electric efficiency reductions when the engine operates at low loads. On the other hand, the diesel mass flow rate presents an expected tendency according to original engine operation as conventional diesel, having lower diesel mass flow rates when both electric power and substitution level decrease, as shown in Figure 8. From the above, engine electric power inefficiencies were compensated by the natural gas mass flow rate increasing. Knowing the incoming flows to engine intake system, the equivalence ratio can be obtained as shown in Figure 9. The equivalence ratio results agree with reported values on several studies of dual fuel diesel engines, increasing the values when higher electric power is demanded [17].
Now, the results for electric and thermal efficiencies are shown in Figure 10 and Figure 11, respectively. For higher electric powers, the higher electric efficiency is achieved. The highest value for electric efficiency was close to 27 % , but had a pronounced decreasing until 20 % when the substitution level increases from 20 % up to 60 % . For electric power from 2.5 kW up to 14.4 kW, maximum values of electric efficiency were achieved for conventional diesel engine operation (i.e., 0 % of substitution level) [42,43], then decreasing linearly when substitution level grows. Figure 11 shows the results for thermal efficiency, substitution level has a slight effect on thermal efficiency decreasing, having values from 24.5 % up to 30 % at electric power from 7.5 kW up to 17.7 kW. A thermal efficiency atypical behavior is obtained at 2.5 kW of electric power, with the highest value close to 32 % , but also with fast decreasing close to 23 % when the substitution level grows up to 25 % .
As a combination of electric and thermal efficiencies, the overall efficiency is shown in Figure 12, which shows a matching behavior to the electric efficiency shown in Figure 10, so that a higher overall efficiency is obtained for higher electric power, presenting marked drops for electric powers such as 2.5 kW, 7.5 kW and 10.6 kW, but less marked decreases for 14.4 kW and 17.7 kW of electric power. The maximum overall efficiency was close to 52 % . Then, the results for heat-electricity ratio are presented in Figure 13, showing that higher heat is recovered, relative to electric power, for lower electric power, presenting values between from 4 to 6 in case of 2.5 kW, meanwhile at higher electric power, the heat-electricity ratio tends to keep from 1 up to 2, i.e., the thermal power recovered is equal to electric power generation.

5.2. Steady State Model Random Validation

Trying to know the model accuracy at a random micro-cogeneration system operation, depending on electric power generation and natural gas availability that require changes of the substitution level, the micro-cogeneration system was operated at several random operation points. In that way, Table 8 lists the random operation points and their respective electric power and substitution level achieved. It must be known that random operation points were not used to predict model coefficients nor model calibration.
Thus, a comparison between experimental results and model predicted values can be shown with respective errors. Figure 14 shows efficiency prediction and measured values, and the maximum error was close to 20 % at 12.5 kW of electric power generation and 37 % of substitution level for thermal efficiency, which agrees with the results in Figure 11, where thermal efficiency presents the largest differences between predicted and measured value. Overall efficiency presents a similar trend because it is affected by thermal efficiency behavior but having a maximum error that is slightly lower, around 15 % . In the case of electric efficiency, the error remains lower than other efficiencies, with an average values close to 5 % .
Additionally, for mass flow rates such as diesel fuel and natural gas, and intake air, a comparison between experimental results and model predicted values are shown with their errors on Figure 15 and Figure 16, respectively. For the mass flow rate of fuels, maximum errors are around 9 % and 7 % for diesel and natural gas, respectively. However, along all random operations, the average error was close to 2 % for both fuels. Regarding to air mass flow rate, accurate results were obtained, having a maximum error around 1.5 % for the first operation point, but then going lower than 0.5 % .

5.3. Transient State Model Validation

Micro-CHP systems can offer wide applications, as required water outlet temperature T w , o in industry, heating processes, also water flow rate m ˙ w could change as well, so that transient validation aims to know the predictability for water outlet temperature. Additional temperatures were measured in micro-cogeneration system operation as well as steady state operation, as Figure 17 shows. Additionally, such operation points are shown in Table 9.
Energy balances are obtained for each control volume as it is shown by Equations (12)–(15) for any operation state, i.e., steady and transient of the micro-cogeneration system. It is necessary to keep in mind that, in the present study, only heat from exhaust gases was recovered, excluding the heat from engine water coolant, which allows us to obtain the above equations. Applying measured temperature values from Figure 17 on the above equations assuming steady state can be obtained results for both overall heat transfer coefficients U A H X and U A L , and even correlation as a function of electric power generation and substitution level. Thus, Equation (23) can be obtained, representing the calibrated model for overall heat transfer coefficient between HRWCV and CECV.
U A H X = 0.078468 0.003424 P E + 0.000133 P E 2 0.137199 Z + 0.377765 Z 2 + 0.035587 P E Z 0.080694 P E Z 2 0.001662 P E 2 Z + 0.003437 P E 2 Z 2
Transient state validation for micro-cogeneration system operation was carried out by comparison between experiments and modeling results, applying a step input of two random electric power levels for the combustion engine, corresponding from 9 kW to 12.5 kW. Data for each variable were collected per second during states as steady–transient–steady. The overall thermal capacitance M C w for HRWCV was equal to 180 kJ / K . Therefore, heat recovery water outlet temperature T w , o behavior can be obtained for both measured and modeling results, which are shown in Figure 18. Additionally, measured inlet water heat recovery temperature T w , i and exhaust gas temperature are shown in Figure 18.
Figure 18 shows the modeling result for the heat recovery water outlet temperature, which fits to measured values even when it is used as a constant overall heat transfer coefficient U A H X for a determined electric power and a substitution level by using Equation (23). Differences are obtained just at the beginning of the input steps, but such a difference decreases when the micro-cogeneration system is achieving the steady state operation.

6. Conclusions

For experimentation, it is necessary to properly choose the heat exchanger system to keep the lowest engine back pressure to hold the engine thermal efficiency as high as possible, low pumping energy, low pricing and maintenance ratio, high heat recovery properties, and no contamination of heat recovery water for possible uses in processes. Thus, compact heat exchanges contain enough features over other types of heat exchangers as shell and tubes, and even some ejectors considered in the first stage of study. The heat exchanger selected in this study offered pressure drops of 17.38 kPa ( 0.174 bar) and 11.37 kPa ( 0.114 bar) for the tube and fin sides, respectively. There was a high value for heat transfer area and volume ratio of 1474 m 2 / m 3 with no crossing flow to keep clean water.
In this study, thermal isolation of a engine exhaust system and exhaust gases diffuser was applied by using ceramic fiber blanket. The length of the engine exhaust gases system is close to 2.5 m, which induced high thermal energy losses prior to the heat exchanger.
A gravimetric method by using a weighing scale for liquid fuel allows an easier and precise measurement task, even though the engine has a large diesel flow return from the injection system to fuel tank, due mechanical pumps. However, precise ground level and vibration isolation must be assured for the weighing scale. On the other hand, for the measurement and supply of NG fuel flow, the sonic flow meter allows us to reduce the operation time at dual mode even without automatic close loop control. High accuracy and large mass flow range is featured adjusting only the gas supply pressure as probed in previous studies. The sonic flow meter can be improved by implementing electronic pressure control valves coupled to a PID strategy.
Now, regarding the thermal performance of dual fuel diesel engines, it lets to reduce the electric power generation costs related with heavy fossil fuel consumption as diesel being replaced by natural gas using or even biogas from local sources, more diesel replacement represents a higher substitution level for engine operating. The highest substitution level was around 60 % at an electric power generation of 17.7 kW. However, higher substitution levels (higher than 60 % ) were limited by the electric power generation control system, which does not allow us to achieve a combustion engine rated power or electric power generation greater than 17.7 kW. The maximum substitution level was 23 % at the lowest electric power generation of 2.5 kW, but presenting an unstable operation of the cogeneration system followed by the lowest electric efficiencies of 4 % .
A higher exhaust gas temperature is achieved for higher substitution levels; however, thermal efficiency decreases by small amounts. Thermal efficiency decreases by larger amounts when increasing the electric power generation compared with substitution level increases. The highest thermal efficiency was 33 % at the electric power generation of 10.6 kW and substitution level of 40 % . The highest electric efficiency was 28 % at the electric power generation of 17.7 kW and substitution level of 24 % . The overall efficiency presents similar behavior, reaching its maximum value of 53 % at 14.4 kW and a substitution level of 0 % . Recovering heat from engine water coolant could allow us to reach higher overall efficiency for cogeneration system, some studies report values close to 75 % .
Additionally, the maximum error presented was close to 20 % associated to thermal efficiency. However, errors for all other variables were lower than 10 % for most of the cogeneration system operation points. The modeling transient result for outlet water heat recovery temperature fits to the experimental measuring results. Therefore, calibrated empirical correlations make it possible to control a cogeneration system depending on both demand for electric power generation or thermal power requirements accurately.
Finally, for the next stages of this study, the authors want to implement not only heat from exhaust gases, but also from an engine cooling system, with an electronic close loop control for the water temperature and substitution level as set-points based on NG availability, likely using electronically controlled injection system as a common rail for the engine, allowing us to change the injection timing aiming to improve both electric and thermal efficiencies, and also allowing us to achieve higher substitution levels even at low electric powers.

Author Contributions

Conceptualization, E.S.C.M., S.H.Q. and I.D.B.; methodology, E.S.C.M. and S.H.Q.; interface, E.S.C.M.; validation, E.S.C.M. and S.H.Q.; formal analysis, E.S.C.M.; investigation, E.S.C.M.; writing—original draft preparation, E.S.C.M.; writing—review and editing, E.S.C.M., S.H.Q. and I.D.B.; supervision, I.D.B.; project administration, E.S.C.M.; funding acquisition, E.S.C.M. and I.D.B. All authors have read and agreed to the published version of the manuscript.

Funding

To Colombia Scientific Program within the framework of the call Ecosistema Científico (Contract No. FP44842- 218-2018) and to University of Antioquia thought the research project identified by the code PRG 2017-16268.

Data Availability Statement

The data and calculations presented in this study are not available due to restrictions established by the funders of the project, for which such data and calculations must be confidential for eight years to protect the possible registration of future developments.

Acknowledgments

The authors gratefully acknowledge the financial support provided by the Colombia Scientific Program within the framework of the call Ecosistema Científico (Contract No. FP44842- 218-2018). The authors also acknowledge the financial support provided by the University of Antioquia thought the research project identified by the code PRG 2017-16268.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations and Nomenclature

The following abbreviations, variables and symbols are used along this article:
CECVCombustion Engine Control Volume
CHPCombined Heat and Power
CICompression Ignition Engine
DECOGDual Fuel Engine Micro-Cogeneration
DIDirect Injection
ECCVEnergy Conversion Control Volume
HRWCVHeat Recovery Water Control Volume
ICEInternal Combustion Engine
IEAInternational Energy Agency
NGNatural Gas
PEMFCPolymer Electrolyte Membrane Fuel Cell
PIDProportional, Integral and Derivative Control
SOFCSolid Oxide Fuel Cell
α Pressure exponent
α f Heat transfer area-volume ratio m 2 / m 3
β Temperature exponent
γ Specific heat ratio for NG
ε Effectiveness
η e g Electric generator efficiency%
η E Electric efficiency%
η T Thermal efficiency%
η O Overall efficiency%
ϕ Equivalence ratio
ϕ D , s t q Diesel equivalence ratio kgDiesel / kgAir
ϕ N G , s t q Natural gas equivalence ratio kgNG / kgAir
Δ p t Tubes side pressure drop kPa
Δ p f Fins side pressure drop kPa
a i Coefficient i for η E polynomial
A * Nozzle critical area m 2
b i Coefficient i for η T polynomial
cFluids capacitance ratio
c i Coefficient i for m ˙ D polynomial
C d Discharge coefficient
d i Coefficient i for m ˙ N G polynomial
e i Coefficient i for m ˙ a polynomial
f m Engine factor
F a Altitude correction factor
H ˙ a Air total enthalpy kW
H ˙ f Fuel total enthalpy kW
H ˙ e x h Exhaust gases total enthalpy kW
L H V D Diesel lower heating value kJ / kg
L H V N G Diesel lower heating value kJ / kg
m ˙ a Air mass flow rate g / s
m ˙ d m Declared fuel consumption g / s
m ˙ D Diesel mass flow rate in dual mode mg / s
m ˙ D , c Diesel mass flow rate in conventional mode mg / s
m ˙ e x h Exhaust gases mass flow rate kg / s
m ˙ N G Flujo másico de combustible mg / s
m ˙ w Heat recovery water mass flow rate kg / s
m ˙ c p w Thermal capacitance for heat recovery water kW / K
M N G Natural gas molar mass kg / kmol
M C c e Overall thermal capacitance for CECV kJ / K
M C w Overall thermal capacitance for HRWCV kJ / K
N t Heat exchanger total tubes
n f Heat exchanger total fins
NEngine speed rpm
N T U Number of thermal units
p 0 Natural gas supply pressure bar
P a t m Inlet air supply pressure at test conditions kPa
P b c Corrected power kW
P b r Prime rated power kW
P E Electric power kW
P E c maximum electric corrected power kW
P T thermal power kW
P Q Chemical power kW
QFuel delivery Mg / L · cycle
Q a Actual possible heat recovery W
Q H X Heat rate to HRWCV kW
Q m a x Heat recovery potential W
Q L Heat lost from CECV to surroundings kW
RPressure ratio
R N G NG ideal gas constant kPa · m 3 / kg · K
R H E Heat-electric power ratio
T 0 Natural gas supply temperature K
T a Inlet air temperature at test conditions o C
T c e CECV temperature K
T e x h Exhaust gases temperature K
T s u r r Surroundings temperature K
T w , i Inlet water temperature K
T w , o Outlet water temperature K
U A H X Overall heat transfer coef. CECV to HRWCV kW / K
U A L Overall heat transfer coef. CECV to surr. kW / K
V d Cylinder capacity L
V t Heat exchanger volume m 3
V ˙ w Heat recovery water volumetric flow rate L / min
ZSubstitution level

References

  1. Shan, X.; Qian, Y.; Zhu, L.; Lu, X. Effects of EGR rate and hydrogen/carbon monoxide ratio on combustion and emission characteristics of biogas/diesel dual fuel combustion engine. Fuel 2016, 181, 1050–1057. [Google Scholar] [CrossRef]
  2. Wang, Z.; Zhao, Z.; Wang, D.; Tan, M.; Han, Y.; Liu, Z.; Dou, H. Impact of pilot diesel ignition mode on combustion and emissions characteristics of a diesel/natural gas dual fuel heavy-duty engine. Fuel 2016, 167, 248–256. [Google Scholar] [CrossRef]
  3. Dimitriou, P.; Tsujimura, T.; Suzuki, Y. Hydrogen-diesel dual-fuel engine optimization for CHP systems. Energy 2018, 160, 740–752. [Google Scholar] [CrossRef]
  4. Simio, L.D.; Gambino, M.; Iannaccone, S. Feasibility of a dual-fuel engine fuelled with waste vegetable oil and municipal organic fraction for power generation in urban areas. J. Combust. 2012, 2012, 176037. [Google Scholar] [CrossRef]
  5. Bedoya, I.D.; Arrieta, A.A.; Cadavid, F.J. Effects of mixing system and pilot fuel quality on diesel-biogas dual fuel engine performance. Bioresour. Technol. 2009, 100, 6624–6629. [Google Scholar] [CrossRef] [PubMed]
  6. Akkouche, N.; Loubar, K.; Nepveu, F.; Kadi, M.E.A.; Tazerout, M. Micro-combined heat and power using dual fuel engine and biogas from discontinuous anaerobic digestion. Energy Convers. Manag. 2020, 205, 112407. [Google Scholar] [CrossRef]
  7. Wierzbicki, S. Laboratory Control and Measurement System of a Dual-Fuel Compression Ignition Combustion Engine Operating in a Cogeneration System; Trans Tech Publications Ltd.: Stafa-Zurich, Switzerland, 2014; Volume 210, pp. 200–205. [Google Scholar] [CrossRef]
  8. Mokhatab, S.; Poe, W.A.; Speight, J.G. Handbook of Natural Gas-Transmission and Processing; Elsevier: Amsterdam, The Netherlands, 2006; p. 672. [Google Scholar]
  9. Gallagher, J.E. Natural Gas Measurement Handbook; Gulf Publishing Company: Houston, TX, USA, 2006; p. 496. [Google Scholar]
  10. Guo, B.; Ghalambor, A. Natural Gas Engineering Handbook; Gulf Publishing Company: Houston, TX, USA, 2005; p. 457. [Google Scholar]
  11. Walker, N.R.; Wissink, M.L.; DelVescovo, D.A.; Reitz, R.D. Natural Gas for High Load Dual-Fuel Reactivity Controlled Compression Ignition (RCCI) in Heavy-Duty Engines; ASME: New York, NY, USA, 2014; p. V001T03A016. [Google Scholar] [CrossRef]
  12. Hallquist, A.M.; Jerksjö, M.; Fallgren, H.; Westerlund, J.; Sjödin, A. Particle and gaseous emissions from individual diesel and CNG buses. Atmos. Chem. Phys. 2013, 13, 5337–5350. [Google Scholar] [CrossRef]
  13. Amirante, R.; Distaso, E.; Iorio, S.D.; Sementa, P.; Tamburrano, P.; Reitz, R.D.V.B.M. Effects of Natural Gas Composition on Performance and Regulated, Greenhouse Gas and Particulate Emissions in Spark-Ignition Engines. Energy Convers. Manag. 2017, 143, 338–347. [Google Scholar] [CrossRef]
  14. Zirngibl, S.; Wachtmeister, G. Extensive Investigation of a Common Rail Diesel Injector Regarding Injection Characteristics and the Resulting Influences on the Dual Fuel Pilot Injection Combustion Process; SAE International: Warrendale, PA, USA, 2016. [Google Scholar] [CrossRef]
  15. O’Connell, N.; Röll, A.; Lechner, R.; Luo, T.; Brautsch, M. PODE-blend as pilot fuel in a biomethane dual fuel engine: Experimental analysis of performance, combustion and emissions characteristics. Renew. Energy 2019, 143, 101–111. [Google Scholar] [CrossRef]
  16. Grosshans, G.; Litzler, M. Pielstick Experience with Dual-Fuel Engines and Cogeneration. J. Eng. Gas Turbines Power 1988, 110, 349–355. [Google Scholar] [CrossRef]
  17. Heywood, J.B. Internal Combustion Engine Fundamentals; McGraw-Hill Education: New York, NY, USA, 1988; p. 930. [Google Scholar]
  18. Lakshmanan, T.; Nagarajan, G. Experimental investigation on dual fuel operation of acetylene in a DI diesel engine. Fuel Process. Technol. 2010, 91, 496–503. [Google Scholar] [CrossRef]
  19. Alanne, K.; Söderholm, N.; Sirén, K. Implementation and validation of combustion engine micro-cogeneration routine for the simulation program ida-ice. In Proceedings of the Eleventh International IBPSA Conference, Glasgow, UK, 27–30 July 2009; pp. 33–40. [Google Scholar]
  20. Onovwiona, H.; Ugursal, V. Residential cogeneration systems: Review of the current technology. Renew. Sustain. Energy Rev. 2006, 10, 389–431. [Google Scholar] [CrossRef]
  21. Onovwiona, H.I.; Ugursal, V.I.; Fung, A.S. Modeling of internal combustion engine based cogeneration systems for residential applications. Appl. Therm. Eng. 2007, 27, 848–861. [Google Scholar] [CrossRef]
  22. Petchers, N. Combined Heating, Cooling and Power Handbook: Technologies and Applications; The Fairmont Press: Lilburn, GA, USA, 2003; p. 875. [Google Scholar] [CrossRef]
  23. Maghanki, M.M.; Ghobadian, B.; Najafi, G.; Galogah, R.J. Micro combined heat and power (MCHP) technologies and applications. Renew. Sustain. Energy Rev. 2013, 28, 510–524. [Google Scholar] [CrossRef]
  24. Murugan, S.; Horák, B. A review of micro combined heat and power systems for residential applications. Renew. Sustain. Energy Rev. 2016, 64, 144–162. [Google Scholar] [CrossRef]
  25. Beausoleil-Morrison, I.; Weber, A.; Maréchal, F.; Griffith, B.; Ferguson, A.; Kelly, N.J.; Marechal, F. Specifications for Modelling Fuel Cell and Combustion-Based Residential Cogeneration Devices within Whole-Building Simulation Programs. 2007. 111p. Available online: https://www.iea-ebc.org/Data/publications/EBC_Annex_42_Cogen_Model_Specifications.pdf (accessed on 31 May 2023).
  26. Beausoleil-Morrison, I.; Weber, A.; Maréchal, F.; Griffith, B.; Ferguson, A.; Kelly, N.J.; Marechal, F. Inter-Model Comparative Testing and Empirical Validation of Annex 42 Models for Residential Cogeneration Devices. 2007. 344p. Available online: https://central.bac-lac.gc.ca/.item?id=M154-14-5-2007E&op=pdf&app=Library (accessed on 31 May 2023).
  27. Arndt, U.; Beausoleil-Morrison, I.; Davis, M.; D’haeseleer, W.; Dorer, V.; Entchev, E.; Ferguson, A.; Gusdorf, J.; Kelly, N.; Manning, M.; et al. Experimental Investigation of Residential Cogeneration Devices and Calibration of Annex 42 Models; Natural Resources Canada: Ottawa, ON, Canada, 2007; p. 255. Available online: https://publications.gc.ca/collections/collection_2008/nrcan/M154-14-4-2007E.pdf (accessed on 31 May 2023).
  28. Quintana, S.H.; Mesa, E.S.C.; Bedoya, I.D. DECOG—A dual fuel engine micro-cogeneration model: Development and calibration. Appl. Therm. Eng. 2019, 151, 272–282. [Google Scholar] [CrossRef]
  29. Rosato, A.; Sibilio, S. Calibration and validation of a model for simulating thermal and electric performance of an internal combustion engine-based micro-cogeneration device. Appl. Therm. Eng. 2012, 45–46, 79–98. [Google Scholar] [CrossRef]
  30. Papagiannakis, R.G.; Rakopoulos, C.D.; Hountalas, D.T.; Rakopoulos, D.C. Emission characteristics of high speed, dual fuel, compression ignition engine operating in a wide range of natural gas/diesel fuel proportions. Fuel 2010, 89, 1397–1406. [Google Scholar] [CrossRef]
  31. Hegab, A.; La, A.; Shayler, P. Towards keeping diesel fuel supply and demand in balance: Dual-fuelling of diesel engines with natural gas. Renew. Sustain. Energy Rev. 2016, 70, 666–697. [Google Scholar] [CrossRef]
  32. Mustafi, N.N.; Raine, R.R.; Verhelst, S. Combustion and emissions characteristics of a dual fuel engine operated on alternative gaseous fuels. Fuel 2013, 109, 669–678. [Google Scholar] [CrossRef]
  33. Rosato, A.; Sibilio, S.; Ciampi, G. Energy, environmental and economic dynamic performance assessment of different micro-cogeneration systems in a residential application. Appl. Therm. Eng. 2013, 59, 599–617. [Google Scholar] [CrossRef]
  34. Rosato, A.; Sibilio, S. Energy performance of a micro-cogeneration device during transient and steady-state operation: Experiments and simulations. Appl. Therm. Eng. 2013, 52, 478–491. [Google Scholar] [CrossRef]
  35. J1349:1995; Engine Power Test Code—Spark Ignition and Compression Ignition—As Installed Net Power Rating. SAE International: Pittsburgh, PA, USA, 1995; Volume 4970.
  36. Sodré, J.R.; Soares, S.M. Comparison of engine power correction factors for varying atmospheric conditions. J. Braz. Soc. Mech. Sci. Eng. 2003, 25, 279–284. [Google Scholar] [CrossRef]
  37. Kays, W.M.W.M.; London, A.L.A.L. Compact Heat Exchangers; McGraw-Hill Book Company: New York, NY, USA, 1998; p. 335. [Google Scholar]
  38. Hesselgreaves, J. Compact Heat Exchangers: Selection, Design and Operation; Pergamon Press: Oxford, UK, 2001; p. 416. [Google Scholar]
  39. Quintana, S.; Castaño-Mesa, E.; Marín, S.; Bedoya, I.; Zapata, J. Measurement and control of natural gas mass flow in a dual-fuel engine operating at partial load through sonic nozzles. Flow Meas. Instrum. 2019, 65, 187–194. [Google Scholar] [CrossRef]
  40. Zheng, C.Y.; Wu, J.Y.; Zhai, X.Q.; Yang, G.; Wang, R.Z. Experimental and modeling investigation of an ICE (internal combustion engine) based micro-cogeneration device considering overheat protection controls. Energy 2016, 101, 447–461. [Google Scholar] [CrossRef]
  41. Gumus, M. Effects of volumetric efficiency on the performance and emissions characteristics of a dual fueled (gasoline and LPG) spark ignition engine. Fuel Process. Technol. 2011, 92, 1862–1867. [Google Scholar] [CrossRef]
  42. Atkins, R.D. An Introduction to Engine Testing and Development; SAE International: Pennsylvania, PA, USA, 2009; p. 308. [Google Scholar]
  43. Merker, G.P.; Schwarz, C.; Teichmann, R. Combustion Engines Development-Mixture Formation, Combustion, Emissions and Simulation; Springer: Cham, Switzerland, 2012; p. 660. [Google Scholar] [CrossRef]
Figure 1. Control volumes diagram of Annex 42 model by IEA [25,28].
Figure 1. Control volumes diagram of Annex 42 model by IEA [25,28].
Energies 16 06281 g001
Figure 2. Schematic diagram of the experimental micro-cogeneration system.
Figure 2. Schematic diagram of the experimental micro-cogeneration system.
Energies 16 06281 g002
Figure 3. Conventional diesel mass flow rate at several electric power generations.
Figure 3. Conventional diesel mass flow rate at several electric power generations.
Energies 16 06281 g003
Figure 4. Sonic flow meter operation for 25 Z < 55 ( 0.8890 mm throat diameter nozzle).
Figure 4. Sonic flow meter operation for 25 Z < 55 ( 0.8890 mm throat diameter nozzle).
Energies 16 06281 g004
Figure 5. Labview interface for data acquisition, monitoring and system control.
Figure 5. Labview interface for data acquisition, monitoring and system control.
Energies 16 06281 g005
Figure 6. Air mass flow rate at several electric power and substitution level.
Figure 6. Air mass flow rate at several electric power and substitution level.
Energies 16 06281 g006
Figure 7. Natural gas mass flow rate at several electric power and substitution level.
Figure 7. Natural gas mass flow rate at several electric power and substitution level.
Energies 16 06281 g007
Figure 8. Diesel mass flow rate at several electric power and substitution level.
Figure 8. Diesel mass flow rate at several electric power and substitution level.
Energies 16 06281 g008
Figure 9. Equivalence ratio at several electric power and substitution level.
Figure 9. Equivalence ratio at several electric power and substitution level.
Energies 16 06281 g009
Figure 10. Electric efficiency at several electric power generation and substitution level.
Figure 10. Electric efficiency at several electric power generation and substitution level.
Energies 16 06281 g010
Figure 11. Thermal efficiency at several electric power and substitution level.
Figure 11. Thermal efficiency at several electric power and substitution level.
Energies 16 06281 g011
Figure 12. Overall efficiency at several electric power and substitution level.
Figure 12. Overall efficiency at several electric power and substitution level.
Energies 16 06281 g012
Figure 13. Heat–electricity ratio at several electric power and substitution level.
Figure 13. Heat–electricity ratio at several electric power and substitution level.
Energies 16 06281 g013
Figure 14. Efficiencies comparison for measured and predicted values.
Figure 14. Efficiencies comparison for measured and predicted values.
Energies 16 06281 g014
Figure 15. Diesel and natural gas mass flow rates comparison for measured and predicted values.
Figure 15. Diesel and natural gas mass flow rates comparison for measured and predicted values.
Energies 16 06281 g015
Figure 16. Air mass flow rate comparison for measured and predicted values.
Figure 16. Air mass flow rate comparison for measured and predicted values.
Energies 16 06281 g016
Figure 17. Measured temperatures related to steady state operation.
Figure 17. Measured temperatures related to steady state operation.
Energies 16 06281 g017
Figure 18. Prediction of heat recovery water outlet temperature.
Figure 18. Prediction of heat recovery water outlet temperature.
Energies 16 06281 g018
Table 1. Technical engine specifications.
Table 1. Technical engine specifications.
CharacteristicValue
DesignationJohn Deere 4039D DI
Cylinders numberFour (4) cylinders
Cooling typeWater cooled
Charge aspirationNaturally aspired
Displacement3.9 L
Compression ratio17.8:1
Bore × Stroke106.5 mm × 110 mm
Rated power41 kW @ 1800 rpm
Table 2. Natural gas (NG) composition and properties.
Table 2. Natural gas (NG) composition and properties.
Species CH 4 C 2 H 6 C 3 H 8 C 4 H 10 CO 2 N 2
Vol. Comp. 84.70 8.820 3.110 1.050 0.740 1.580
Mass Comp. 71.54 13.97 7.224 3.215 1.719 2.335
DescriptionSymbolUnitsValue
NG gas constant R N G kPa · m 3 / kg · K 0.4389
Molar mass M N G kg / kmol 18.943
Stoich. fuel-air ratio ϕ N G , s t q kgNG / kgAir 0.0619
Specific heat ratio γ @ 25 C - 1.2760
Lower heating value L H V N G kJ / kg 47,249
Table 3. Values for electric power generation corrections.
Table 3. Values for electric power generation corrections.
DescriptionSymbolUnitsValue
Cylinder capacity V d L3.9
Engine speed for maximum powerNrpm1800
Declared fuel consumption m ˙ d m L/h (g/s)10.5 (2.48)
Pressure ratio *R1
Pressure exponent * α 1
Temperature exponent * β 0.7
Altitude correction factor F a 0.944
Electric generator efficiency η e g 0.88–0.92
Prime rated power (PRP) P b r kW41.0
Corrected power P b c kW38.7
Maximum electric corrected power P E c kW34.8
* Applied values for naturally aspired engines.
Table 4. Initial experimental data for conventional diesel and desired temperatures.
Table 4. Initial experimental data for conventional diesel and desired temperatures.
DescriptionSymbolUnitsValue
Inlet heat recovery water temperature T w , i ° C20
Outlet heat recovery water temperature T w , o ° C70
Exhaust gases temperature (average) T e x h ° C285
Exhaust gases mass flow rate m ˙ e x h kg/s0.239
Heat recovery water mass flow rate m ˙ w kg/s0.225
Heat recovery water volumetric flow rate V ˙ w L/min13.62
Table 5. Heat exchanger characteristics and estimations.
Table 5. Heat exchanger characteristics and estimations.
DescriptionSymbolUnitsValue
Heat exchanger volume V t m30.0018
Heat transfer area-volume ratio α f m2/m31474
Total tubes n t -29
Total fins n f -11,100
Tubes side pressure drop Δ p t kPa17.38
Fins side pressure drop Δ p f kPa11.37
Effectiveness ε -0.331
Heat recovery potential (ideal) Q m a x W67,242
Actual possible heat recovery Q a W22,195
Table 6. Sonic nozzle flow meter operation parameters.
Table 6. Sonic nozzle flow meter operation parameters.
Throat Diameter C d Operation Range
0.4064  mm ( 0.016  in) 0.8853 0 < Z < 15
0.6096  mm ( 0.024  in) 0.8982 15 Z < 25
0.8890  mm ( 0.035  in) 0.9219 25 Z < 55
1.3208  mm ( 0.052  in) 0.8999 55 Z < 70
Table 7. Coefficients for model correlations.
Table 7. Coefficients for model correlations.
Coefficients η E η T m ˙ D m ˙ NG m ˙ a
0.006248 0.331000 719.853 0.79160 50,139.1
P E 0.029212 0.004081 21.2416 0.0 110.647
P E 2 0.000832 0.000001 1.44269 0.0 0.0
Z 0.101548 1.011470 549.056 0.0 −14,331.6
Z 2 0.741322 1.766230 0.0 13,943.517,331.9
P E Z 0.071513 0.163051 55.3363 387.704 475.929
P E Z 2 0.148371 0.299428 0.0 2031.57 537.707
P E 2 Z 0.004091 0.006119 0.0 19.8360 0.0
P E 2 Z 2 0.007169 0.011386 0.0 78.1700 0.0
Table 8. Random operation points for model validation.
Table 8. Random operation points for model validation.
Operation Point12345678
P E [kW]9 12.5 12.5 14.4 14.4 14.4 1616
Z [-] 0.00 0.00 37.2 36.9 24.8 53.8 0.00 34.3
Table 9. Operation points for measured temperatures.
Table 9. Operation points for measured temperatures.
Operation Point12345678
P E [kW] 2.5 7.5 10.6 14.4 17.7 2.5 2.5 7.5
Z [-]00000 16.6 23.2 24.2
Operation Point 9 10 11 12 13 14 15 16
P E [kW] 7.5 10.6 14.4 14.4 14.4 17.7 17.7 17.7
Z [-] 28.6 30.6 24.8 36.9 53.8 24.6 45.9 60.7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Castaño Mesa, E.S.; Quintana, S.H.; Bedoya, I.D. Development of a Dual Fuel ICE-Based Micro-CHP System and Experimental Evaluation of Its Performance at Light Loads Using Natural Gas as Primary Fuel. Energies 2023, 16, 6281. https://doi.org/10.3390/en16176281

AMA Style

Castaño Mesa ES, Quintana SH, Bedoya ID. Development of a Dual Fuel ICE-Based Micro-CHP System and Experimental Evaluation of Its Performance at Light Loads Using Natural Gas as Primary Fuel. Energies. 2023; 16(17):6281. https://doi.org/10.3390/en16176281

Chicago/Turabian Style

Castaño Mesa, Edisson S., Sebastián H. Quintana, and Iván D. Bedoya. 2023. "Development of a Dual Fuel ICE-Based Micro-CHP System and Experimental Evaluation of Its Performance at Light Loads Using Natural Gas as Primary Fuel" Energies 16, no. 17: 6281. https://doi.org/10.3390/en16176281

APA Style

Castaño Mesa, E. S., Quintana, S. H., & Bedoya, I. D. (2023). Development of a Dual Fuel ICE-Based Micro-CHP System and Experimental Evaluation of Its Performance at Light Loads Using Natural Gas as Primary Fuel. Energies, 16(17), 6281. https://doi.org/10.3390/en16176281

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

Article Metrics

Back to TopTop