Using cold EGR method with variable venturi and turbocharger has very significant effect, simultaneously on the reduction of NO_{X} and grime simultaneously. EGR cooler is one of the most important parts in the cold EGR circuit. In this paper optimum design of cooler for working in different percentages of EGR, besides determination of optimum temperature of exhausted gases, efficiency growth, weight reduction, dimension and expenditures reduction, sediment reduction, and optimum performance by using gasoil which has significant amounts of brimstone, is investigated and optimized.
Cold EGR; NOX; Cooler
A_{to} total heat transfer area (m^{2})
A_{o} external heat transfer area of tube (m^{2})
A_{i} internal heat transfer area of tube (m^{2})
A_{w} logarithmic mean heat transfer area (m^{2})
C_{ps} specific heat of shell fluid (cooler fluid) (J/kg K)
C_{p} specific heat of gas (J/kg K)
D_{e} equivalent hydraulic diameter (m)
D_{i} tube internal diameter (m)
D_{o} tube external diameter (m)
D_{s} shell internal diameter (m)
F correction factor of LMTD at multi pass flows
h_{id} heat transfer coefficient of shell in ideal flow model (W/m^{2} K)
h_{o} heat transfer coefficient in shell side (W/m^{2} K)
h_{t} heat transfer coefficient in tube side (W/m^{2} K)
J_{l} correction factor for effects of wall seals
J_{b} correction factor for cross flows
J_{s} correction factor for walls distance at inlet and output sections
J_{r} correction factor for Reynolds number less than 100
J_{c} correction factor for fluid shear effect near the wall
k_{t} thermal conductivity of tube side fluid (W/m K)
k_{w} thermal conductivity of wall (W/m K)
N_{t} number of tubes
NTU number of heat transfer unit
 mass flow rate of shell side fluid
Q actual heat transfer rate (W)
U_{o} overall heat transfer coefficient (W/m^{2} K)
R_{i} tube internal fouling resistance (m^{2} K/W)
R_{o} tube external fouling resistance (m^{2} K/W)
t_{w} thickness of tube wall (m)
 cooling fluid viscosity (water) (N s/m^{2})
 hot fluid viscosity (gas) (N s/m^{2})
 exchanger effectiveness
 temperature dependent viscosity correction factor
The most important pollutants of diesel engines are Nitrogen Oxides (NO_{X}) and grime particles. One of most effective methods for NO_{X} pollutants reduction is restoring exhausted gases (EGR) into the combustion chamber; this action would be possible by three following methods: Thermal, Chemical and Dilution [1].
By cooling EGR gases (Cold EGR), more percent of exhausted NO_{X} pollutants will be decreased. In cold EGR, restored gases are cooled by a heat exchanger and amount of absorbed heat by EGR will be increased [2] and [3]. In other side, by decreasing the restored gases temperature, durability of EGR valves and other parts which are depended on EGR will be increased.
In diesel engines, most of produced Nitrogen Oxides are thermal, and by increasing combustion temperature, the rate of NO production will be higher. In other side, carbon suspended particles are the results of imperfect combustion and reduction of temperature. However, the most important effort of pollutants reduction and estimating new environmental standards in diesel engines is to reduce Nitrogen Oxides and grime simultaneously.
For diminishing Nitrogen Oxides and particles with no brunt to fuel economics, engine systems should be designed and controlled carefully to acquire optimum conditions of inlet temperature and pressure, time and amount of fuel injection, and quantity and schedule of input EGR. Prediction of all actual optimization variables is a very complicated and important obligation.
Numerous experiences on different EGR systems of heavy diesel engines, in scope of measuring performance affections, special fuel consumption and pollution, reveal that cold EGR (synthesis of high pressure exhausted gases and high pressure input air) or in other words short circuit EGR method with a variable venturi, with increased injection pressure (common rail) and turbocharger using, in comparison with the low pressure EGR system (synthesis of hot exhausted gases and input air before turbocharger), or in the other words long circuit EGR system, is more effective for NO_{X} reduction and obtaining minimum amount of increased specific fuel consumption. Also short circuit EGR obstructs sediments in compressor and intercooler (if there is); in other side retardation of injection, increases the specific fuel consumption [4] and [5]. Fig. 1a shows these two systems.

Figure 1. (A) High pressure EGR system and (B) low pressure EGR system.

Briefly about reduction of NO_{X} and grime in cold EGR, we can mention two important factors: firstly input load temperature reduction, lengthens combustion retardation and consequently induces more homogeneity between fuel steam and Oxygen synthesis; and secondly trapped air mass in the cylinder in IVC will increase which induces λ increment and oxygen concentration. Input Oxygen density increment will reduce dilution effects. Also in this situation cylinder pressure increment and consequently fortification of compress and combustion stroke cause effective pressure optimization which will optimize fuel consumption.
One of important sections of cold EGR design is designing and manufacturing of heat exchanger for gas cooling, which since now will be called EGR cooler. EGR coolers are mostly in shelltube type [6], [7] and [8]; in this case gas flows through tubes and coolant (engine jacket water) flows into the shell. This facilitates tubes cleanings and controlling of cooler exhausted gas temperature.
Inspected results reveal that although tubes numbers increment will increase the shell dimension (which will decrease heat transfer coefficient), in other side will increase heat exchanger and heat transfer area and decrease pressuredrop [7]. However, tubes numbers increment will increase weight and production expenditures.
By using passive heat transfer methods such as spiral tubes we can increase heat transfer efficiency and decrease the length of heat exchanger [9]. Fig. 2 shows four types of enhanced tubes and simple tubes, and compares them with regard to rate of heat transfer, pressuredrop and manufacturing cost [8]. Also a specimen of spiral tubes is shown in Fig. 3. About heat transfer increment, spiral tubes have more efficiency in comparison with other extended tubes; also they have appropriate production cost and pressuredrop.

Figure 2. Comparison of different tubes which are used in heat exchangers [8].


Figure 3. Specimen of spiral tube [10].

Using spiral tubes will be very effective in scope of equalizing fluid speed in the shell, local boiling and thermal stresses obstruction, and intensifying perpetuity of cooler [8].
Pressuredrop at the gas, and subsequently speed decrement of it in tubes, is a design restrictive factor. With regard to EGR flow different percentages, Reynolds number ranges between 2000 and 9000 at the cooler; it is between transient flow ranges which does not have precise figured results for calculating heat transfer coefficient. So estimating friction coefficient by using analytic methods will not get a trusted quantity.
Also we cannot calculate actual pressuredrop during entrance diffuser. Moreover a stagnation pressure point will be created at return section. This happens at the end of diffuser, flow striking section to the wall, and boundary layer contact at conical section; which has a very important effect on the amount of crossing flow, particularly at side tubes. Therefore in scope of investigation of diffuser entrance diameter effects on the uniform and effective distribution of returning gases into all tubes, and cooler performance estimation, some calculations are done by CFD method.
Significant portion of pressuredrop in diffuser is due to flow shear stress at returning section.
Despite rational pressuredrop, appropriateness with working conditions gets designed cooler to have higher efficiency in comparison with typical coolers [6], [7] and [8].
Governing equations encompass continuity, momentum, energy, correlated pressuredrop and heat transfer relations. Applying different EGR percentages and presence of various design parameters such as arrangements, sizes, number of passages, and type of tubes leads us to use trial and error method which brings complexity in design; therefore, a code is written in Fortran 90.
Elementary design for different EGR percentages firstly is done by LMTD method and tubular exchanger manufacturers’ association standard (TEMA) and then is probed and proved by the number of transfer units’ method. Also permissive amount of pressuredrop is investigated by related equations in the tubes.
In the LMTD method heat transfer area is calculated by the following equation [11]:

(1) 
is LMTD of Counter flow and is described as follows:

(2) 
Major calculations in EGR cooler design are related to U_{o} calculation; this coefficient is depended on fouling resistance and heat transfer coefficients of shell and tube sides.

(3) 

(4) 
In cooler tubes, flow is in transient conditions and there is no precise solution for heat transfer coefficient calculation. One of the best graphs to obtain this coefficient in tubes is Tate & Sieder which is shown in Fig. 4.

Figure 4. Heat transfer coefficient in the tube by Tate & Sieder [12].

We can arrive to the following relations by fitting this curve in this section for different proportions of L/D, where L is the length of tube and D is the internal diameter of tube.

(5) 
h_{i} is heat transfer coefficient in tube and is acquired by following equation:

(6) 
where
In the method effectiveness () is described as follows:

(7) 

(8) 
Basic equation for mean heat transfer coefficient calculation in shell side is explained on the basis of Bell–Delaware method [12] as follows:

(9) 
where h_{id} is heat transfer coefficient for completely cross flow with ideal tube set and ’s are correction factors for seals and others:

(10) 
Factors calculations and complete information for heat transfer coefficient and pressuredrop calculations in shell and tube sides are mentioned in [11] and [12] references.
For figuring out number of tubes, if shell diameter is less than 230 mm following relations should be used. For triangular arrangement Eq. (11) and for quadrangular and turned quadrangular arrangement Eq. (12) are used.

(11) 

(12) 
For one and two tube passes is 1 and 0.92 respectively. Calculations reveal that number of tubes for triangular arrangement and circular arrangement is equal.
With regard to problems such as transient flow in tubes, shortness of cooler length, and effects of entrance length, and what are mentioned before, usage of computational fluid dynamics method (CFD) for flow specimens estimation, temperature distribution and pressure in cooler, is very advantageous. Computational fluid dynamics software which is used in this paper is Fluent. In simulations pressure basis algorithms are used for continuous solving of converted differential equations in numerical form. Used algorithm here is SIMPLE. In every solving section, independence of solving trend from mesh form is assessed.
For numerical results validation, firstly a comparison between experimental data [7] and numerical results for a cooler tube simulation is done; then in scope of heat transfer and pressuredrop investigation, simulations for three entrance shapes for the whole cooler are probed. In Fig. 5 experimental and numerical results in tube length for two sections of center and 0.9 mm near the wall are exhibited. In this simulation, results are acquired for inlet gas at 250 °C and 0.25 g/s, in which temperature of cooler water is supposed to be permanent at 90 °C, for this condition acquired Reynolds number is 2000. Also exhausted gas specifications are supposed to be mean specifications of combustion productions.

Figure 5. Gas temperature diversion through a simple tube.

In next stage numerical analysis for cooler and three types of diffuser (cylindrical, nuzzle and conical) with flows of 100 g/s, 150 g/s and 200 g/s is done and compared with the experimental results. Results comparisons are shown in Fig. 6 in pressuredrop form, on flow basis, for a quarter prototype. Moreover, for more information, outlet temperature of cooler, related to flow rate, is shown in Fig. 7. Internal diameter of shell is 78 mm and tube number is 70. Chart shows a good coincidence between experimental and numerical results which maximum deviation is about 7%.

Figure 6. Pressure drop related to flow rate.


Figure 7. Outlet temperature of cooler, related to flow rate.

Spiral tubes usage prompts growth of heat transfer efficiency in tube side and reduction of cooler size and weight. However, it will induce pressuredrop increment in tube side. Some optimizations are done for reducing this negative effect to an acceptable amount. Effective parameters of Spiral tube design on flow specifications and heat transfer, comprise tube pitch, spiral angle, spiral height, and spiral numbers.
As shown in Fig. 8, spiral wall induces fluid spinning near the wall; at this region fluid velocity exceeds plain tube velocity; also boundary layer is thinner than plain tube. Subsequently heat transfer coefficient near the wall is larger and as a result heat transfer rate is greater than before. Results are acquired for a tube with internal diameter of 7.8 mm, pitch of 1.4 mm, length of 200 mm and gas flow of 2.21 kg/h.

Figure 8. A & Bvelocity contours at plain tube and one way spiral tube.

With regard to negative effects of pressuredrop on heat transfer rate in cooler tubes, effects of spiral numbers increment on pressuredrop and heat transfer are probed. As shown in Fig. 9, increase of spiral numbers from one to two (or four), with constant hydraulics diameter, induces central flow symmetry and decreases fluid velocity in comparison with one spiral tube.

Figure 9. A & Bvelocity contours at 2 and 4 way spiral tube.

Also by increasing heat transfer area because of increment of spiral numbers significant decrement is not seen in thermal efficiency in comparison with one way spiral tube.
Obtained results are shown in Figure 10 and Figure 11.

Figure 10. Spiral number effect on thermal effectiveness.


Figure 11. Spiral number effect on pressure drop.

Spiral pitch reduction and spiral height increment, cause heat transfer coefficient growth.
Sediments increment at the tube wall, decreases total heat transfer coefficient; particularly, because of brimstone sediments and condensation of water vapor, a very corrosive and deteriorative environment is created. By fluid velocity increment in tube, sediment creation will be decreased. Fluid swirl in tube, causes sediments reduction and prompts selfclearance.
Despite pressuredrop negative effects, by investigation in increment of numbers, and pitch and height of spirals for clear and fouled tubes, optimum design conditions are selected to attain best situation with regard to efficiency increment. This condition is shown in Fig. 12.

Figure 12. Influence of pitch and height of spiral on effectiveness and pressure drop for clean and fouled tubes.

Gas flow equal dispensation into tubes is one of important subjects in EGR cooler entrance head design. Fig. 13 shows circular arrangement for 19 tubes of 3/8 in. in the shell with 63 mm diameter. Tubes are numbered from 1 to 19 to facilitate tubes flow percentage indication in Fig. 14.

Figure 13. Circular arrangement with denoting tubes’ numbers in scope of flow calculation.


Figure 14. Influence of entrance head diameter to shell diameter on the tubes’ flow.

At low proportions of entrance diameter to shell diameter (D_{1}/D_{2}), flow percent decreases at side tubes and as a result, more flow is passed into the centrals. For example, at D_{1}/D_{2} = 0.4, flow velocity in the central tubes is 35% more than side tubes; this causes pressure reduction at central tubes and cooling water temperature increment near them, which can make local boiling at that section. Local boiling causes cooler vibration and thermal stresses in working time. Also flow rate decrement and subsequently turbulence decrement at side tubes, cause local heat transfer rate reduction. Fluid velocity reduction can prompt sediment increase too. All of these, cause decrement of effectiveness and lifetime of cooler. Fig. 15 shows influence of internal diffuser diameter to shell diameter ratios of 40% and 70%. Temperature increase in the central tubes, because of high flow rate of gas and consequently temperature raise of cooling fluid near them, is apparently visible.

Figure 15. Influence of diffuser entrance diameter to shell diameter ratio (40% & 70%).

In other side, if the proportion of D_{1}/D_{2} be less than before, we will need more length for diffuser to distribute flow appropriately and equally. This will oblige larger cooler length and more pressuredrop in gas side. Based on fittings’ type and EGR tube connectors, with regard to acquired results, suggested ratios of D_{1}/D_{2} are 60% to 75%. Of course amounts of D_{1} parameter are restricted regarding tubes’ diameter.
Cooler design for diesel engine of 4–248 M.F. Tractors, after performance and pollution tests and calculation of maximum quantity of EGR, is done and results are given in Table 1; tests reveal that cold EGR growth, more than 20%, will increase grime amount significantly. With regard to gasoil properties, high percent of brimstone in it, and based on numerous tests, optimum EGR quantity for this engine is 10%; but design is done in a case which works appropriately up to 15%.
Optimum cooler  Normal cooler  

Input gas flow (kg/s)  0.0153  0.015 
Input and output gas temperature (K)  473,723  473,723 
Tube internal diameter (mm)  7.2  7 
Tube external diameter (mm)  9.5  9.5 
Number of tubes  19  19 
Tube length (mm)  284  410 
Pitch proportion  1.25  1.25 
Tube pressure drop (Pa)  316  185 
Water flow (kg/s)  0.1809  0.1894 
Input and output water temperature (°C)  90–85  90–85 
Shell internal diameter  63  63 
High pressure EGR system, with regard to explained benefits of it, is installed on the engine. Because exhausted EGR gases are extracted for cooling operations before turbocharger, they have a temperature about 500 °C. With respect to tests, temperature of exhausted gases after turbocharger is 70–100 °C less than before. We suppose that thermodynamical properties of EGR gases are the same as exhausted engine gases mean specifications, which are combustion productions of sample gasoil with air. The used fuel (gasoil) chemical formulation is C_{12}H_{26.}
Cooling temperature range is very important. In scope of considering brimstone particles effects on the cooling EGR system, some hypothetical and experimental works are done. Researches exhibit that condensation of various brimstone compositions, in every temperature, is dependent on the chemical syntheses of exhausted gases; but dew point of exhausted gases sulfuric acid, by using a fuel with 370 ppm sulfur, is about 130 °C [13]. More temperature decrement causes sulfur and water vapor condensation at backflow gases; consequently sulfuric acid formation causes cooling tubes erosion and deterioration of valves, EGR system connectors, and combustion chamber. Also high temperatures of exhausted gases, impediments appropriate reduction of NO_{X} and grime. With regard to tests, cooler outlet proper temperature is about 170–210 °C.
Acquired results, for optimum cooler, exhibit weight and dimensions decrement in comparison with typical coolers; it is due to high surface density, effectiveness increment, and equal flow distribution into the tubes. Negative effects of pressuredrop are decreased significantly; moreover, cooler has a high perpetuity against thermal stresses.
One way spiral tube usage, increases cooler thermal effectiveness up to 10% in comparison with typical coolers and because of semiuniformity of flow in the tubes, using multi way spiral tubes causes pressuredrop reduction; moreover, it does not have high negative influence on the thermal efficiency, and it is the result of heat transfer area increment. Also increment of fluid spinning velocity helps to sediment reduction, selfcleaning, and heat transfer increment.
Pitch size reduction and spiral depth increment, are effective factors at heat transfer and pressuredrop growth.
Simulations, with regard to pitch size, height, number of spirals and sediment effects, are done to ideal and optimal tube design with balance between heat transfer effectiveness and pressuredrop at EGR cooler working range. Entrance diffuser and inlet diameter to shell diameter ratio, in scope of uniform flow distribution, thermal effectiveness increment, and pressuredrop reduction, are investigated. Standard diffuser shape with 60–70% diameter ratios gives best results.
Published on 12/04/17
Licence: Other
Are you one of the authors of this document?