OPEN ACCESS
The paper aims at the optimization of the design of a biomass boiler under the inspiration of the Constructal Law by Bejan. The boiler is of the smoke tubes typology, fuel being biomass pellets. The smoke tubes are 16 and are placed in a staggered configuration. A model is built in MATLAB environment, based on empirical correlations and the mean log temperature methodology. The analysis is based on the development of a wide parametric analysis that involves variations of diameters, numbers and positioning of the tubes. Results are based on the concept of the overall performance coefficient methodology and investigate both the pressure drops variation and the thermal power generated in the different configurations’ boiler. However, since the aim of the boiler is to guarantee a fixed thermal power output to the users, results are also presented by fixing the thermal power output and adding as variable the volume of tubes adopted as an indicator to quantify materials required to achieve the same output. While from the point of view of pressure losses there is a clear best solution, this is not the case if the material uses minimization and pressure drops minimization are jointly analyzed. This leads to the need to a tradeoff in performance to identify a best case.
boiler, constructal law, modelling, parametric analysis, overall performance coefficient
The "Constructal theory", proposed in 1997 by Bejan is based on a simple principle: for fluids, the best geometry is the one that exerts the lesser resistance to motion; instead for the transfer of thermal energy, the best geometry is the one that generates the least thermal resistance[1]. Constructal theory is the view that the generation of flow configuration is a universal phenomenon of all physics, which is covered by a law of physics (the Constructal law):
“For a finitesize flow system to persist in time (to live) it must evolve such that it provides greater and greater access to the currents that flow through it”[2].
Also, Bejan explains the connection between “Constructal law” and the thermodynamics of flow systems by developing analytical and graphical formulations to maximize flow access in systems with heat and fluid flow irreversibility and freedom to change configuration [3].
Constructal theory has been applied to various fields, from heat and mass transfer applications in engineered systems to biology, geophysics, social dynamics and economics fields.In particular, there are several research works of Constructal Law on heat exchangers [410] and on the optimal shaping of fins with application to heat exchangers [11].
For example, the objective of [12] is the design of a heat exchanger that takes maximum advantage of high heat transfer density achieved by the use of smallscale channels with laminar flow. It is shown that the spacing of the elemental and firstconstruct channels can be optimized such that the overall pumping power required is minimal.
In [13], a Constructal Law analysis was conducted on a blast furnace ironmaking process. The process is optimized and the results show that the hot metal yield and useful energy are, respectively, increased by 3.13% and 2.66% after generalized Constructal optimization.
In [14] the problem of determining the configuration of a stream that must be heated in a fixed space while using minimal fuel is proposed, applied to the heating of a stream of solid metal in a furnace. The results show that the distribution of heaters per unit area depends on the shape of the floor and the density of heaters per unit area varies as x−0.8/p(x).
In other cases, it has been shown how to increase the power generation over the area occupied by the solar chimney configuration[15]. In this case it was found that the power generated per unit of land area is proportional to the length scale of the power plant.
Although energy applications of the Constructal Law are widespread in literature, there is only a limited availability of studies on the specific topic of optimizing the thermal performances of a boiler.
The investigation performed is based on an optimization of the design of a biomass boiler through the application of the Constructal Law.
The work approaches also the perspective of ecodesign, including the use of materials required to achieve a fixed heating power as parameter of the optimization.
In particular the aim of the work is to apply the Constructal Law and the adaptive overall performance coefficient developed by Lorenzini et al. [16] to perform an optimization of the geometrical design of a pellet boiler through a performance index dependent on the overall pressure drop and the quantity of building materials.
The following paragraphs discuss the technical characteristics of the actual boiler, the mathematical model implemented and the methodology chosen for the study.
2.1The boiler
The analysis is based on the study of an existing boiler. Its most relevant technical and geometrical data are reported in Table 1. Figure 1 shows a picture of the boiler.
Table 1. Geometrical and technical data of the boiler
Characteristics 
Unit 
Value 
Diameter 
mm 
50 
Length 
m 
0.6 
Width 
m 
0.6 
Height 
m 
0.6 
Number of tubes 
 
16 
Flow rate smoke tubes 
kg/s 
0.829 
Flow rate water 
kg/s 
0.245 
Thermal power 
kW 
46 
Efficiency 
 
89% 
The boiler is of the smoke tubes typology, fuel being biomass pellets. The smoke tubes are 16 and are placed in a staggered configuration.
Figure 1. Photo of the boiler analyzed
Biomass pellets are burnt in the combustion chamber, hot smokes enter the tubes visible in Figure 1 and leave from a chimney placed on the back of the boiler. Water exchanges heat with the smoke tubes entering from the left in cross flow configuration.
2.2 Mathematical modeling
To analyze the existing system, a stationary mathematical model was developed in MATLAB environment,implementing empirical correlations and the method of the logmean temperature difference to determinate the heat transfer between flows.
The model can undergo change of the boundary conditions,geometrical configuration of the heat exchanger, variations of the physical  chemical characteristics of the fluids, while respecting the laws of the thermofluid dynamics. It can be applied to any kind of similar boilers of the same kind using the same working principles.
The model implements heat transfer in the flow across a bank of tubes through the energy balance equation for a cross flow heat exchanger (Eq. 1) [17].
$\dot{Q}=\left(\dot{m} c_{p}\right)_{t}\left(T_{t, i}T_{t, o}\right)=\left(\dot{m} c_{p}\right)_{s}\left(T_{s, o}T_{s, i}\right)$ (1)
where Q is the heat duty in kW, is the mass flow rate [kg/s], cp is the specific heat capacity [kJ(kgK)1], T is the temperature [K], the subscripts t and s stand for the tubeside and shellside, respectively; and the subscripts i and o stand for the inlet and outlet of tube or shell side, respectively.
Moreover, the heat transferred in the heat exchanger can be calculated from the following equation 2.
$\dot{Q}=F U A \mathcal{A} T_{L M}$ (2)
where F, U,$\Delta T_{L M}$ and A are respectively a correction factor that it is a function of the geometrical arrangement and of the temperatures of both flow, heat transfer coefficient [W (m2K)1], the logarithmic mean temperature difference between the two flows [K] [18,19] and the heat transfer area [m2 ].
Eq. 3 shows the calculation of logarithmic mean temperature difference.
$\Delta \mathrm{T}_{\mathrm{LM}}=\frac{\left(\mathrm{T}_{\mathrm{s}, \mathrm{i}}\mathrm{T}_{\mathrm{t}, \mathrm{o}}\right)\left(\mathrm{T}_{\mathrm{s}, \mathrm{o}}\mathrm{T}_{\mathrm{t}, \mathrm{i}}\right)}{\ln \left(\left(\mathrm{T}_{\mathrm{s}, \mathrm{i}}\mathrm{T}_{\mathrm{t}, \mathrm{o}}\right) /\left(\mathrm{T}_{\mathrm{s}, \mathrm{o}}\mathrm{T}_{\mathrm{t}, \mathrm{i}}\right)\right)}$ (3)
The value of heat transfer area (A) is calculated in Eq.4.
$A=\pi \cdot d_{o} \cdot L \cdot N$ (4)
where L is length of tubes [m], do [m] is outside diameter of the tube and N is the total tube number
The shell volumetric hydraulic diameter (Dh)is calculated as follow [19]:
$D_{h}=\frac{4 A_{0} W}{A}$ (5)
where A0 is the passage area [m2 ] of the crossflow to bank of tubes and it is calculated in Eq. 6 and W is the width of boiler [m].
$A_{o}=\left(P_{T}d_{o}\right) \cdot L \cdot N_{r}$ (6)
where N is the tube number, do [m] is outside diameter of the tube, ST[m] and SL [m] are, respectively, the transversal and longitudinal pitch and PT and PL are the dimensionless pitch,calculated as ratio between pitch and outside diameter and the subscript r stand for the rows of tubes in longitudinal direction.
The thickness of the pipes is equal to 5 mm, as in the real case.
The thermal transmittance (U) is calculated using Eq. 7 and it is expressed in W/(m2 K).
$\frac{1}{U}=\frac{1}{h_{o}}+\mathrm{R}_{\mathrm{f}, \mathrm{o}}^{\prime \prime}+\frac{\mathrm{d}_{\mathrm{o}} \cdot \ln \left(\mathrm{d}_{\mathrm{o}} / \mathrm{d}_{\mathrm{i}}\right)}{2 \mathrm{k}}$$+\mathrm{R}_{\mathrm{f}, \mathrm{i}}^{\prime \prime} \frac{\mathrm{d}_{o}}{\mathrm{d}_{\mathrm{i}}}$$+\frac{1}{h_{i}} \frac{\mathrm{d}_{\mathrm{o}}}{\mathrm{d}_{\mathrm{i}}}$ (7)
where R” is the fouling factor [(m2 K)/ W], and h [W (m2 K)1 ] is the convection heat transfer coefficient and it is calculated as Eq. 8.
$h=\frac{N u \cdot k}{D_{h}}$ (8)
where Nu is the Nusselt number, k [W/(m K)] is the conduction heat transfer coefficient and Dh [m] is the hydraulic diameter.
This model calculates the Nusselt number and the pressure drop for bank of tubes as in the model of Zhukauskas as in [19] (Eqs. 9 and 10 respectively):
$N u_{c o l d}=C_{1} C_{2} \operatorname{Re}_{D, \max }^{m} \operatorname{Pr}^{n}\left(\frac{\operatorname{Pr}}{\operatorname{Pr}_{w}}\right)^{1 / 4}$ (9)
$\Delta P_{c o l d}=\chi N_{r}\left(\rho \frac{u_{\max , \text {cold}}^{2}}{2}\right) f_{\text {cold}}$ (10)
where C1, C2, m, n, are dimensionless parameters that change in function of physical and geometrical characteristic of cross flow heat exchanger arrangement. Correction factor C1 varies as function of Reynolds (Re) and arrangement of tubes.Instead, f and c are the friction factor and its correction factor [18,19]. If the number of rows of tubes is lower than 16, it is necessary to use the correction factor C2.
Pressure drops internal to the tubes are calculated as Eq. 11.
$\Delta P_{h o t}=f_{h o t} \frac{u_{h o t}^{2}}{2} \frac{L}{d_{i}}$ (11)
where f is the friction factor function of Re; u is the speed [ms1]; L [m] is length of tubes and di is the inner diameter of tubes.
The hot smokes stream inlet temperature (Thot,in) and the cold stream inlet temperature (Tcold,in) are fixed, since no modification is suggested to the combustion chamber and no modification in the water supply is expected during the use phase.
Mass flow rate is specified at the inlets of the hot and the cold stream. The fluid velocity is calculated as in Eqs. 12 and 13:
$u_{h o t}=\frac{\dot{m}_{h o t}}{\rho_{h o t} A_{t}}$ (12)
$u_{c o l d}=\frac{\dot{m}_{c o l d}}{\rho_{c o l d} A_{s h e l l}}$ (13)
where u is the velocity of the flow [ms1] and r is the density of the flow [kg m3].
The Reynolds number ReD,max for the foregoing correlation is based on the maximum fluid velocity, calculated in Eq. 14,occurring within the tube bank in Eq. 15.
$\mathrm{Re}_{h o t}=\frac{(\rho u)_{h o t} d_{i}}{\mu}$ (14)
$\operatorname{Re}_{D, \max , \operatorname{cold}}=\frac{\rho_{\operatorname{cold}} u_{\max } d_{o}}{\mu}$ (15)
For the fluid in the shell side, the maximum reachable speed depends on its geometrical shape configuration. The maximum speed is calculated through Eq. 16.
$u_{\max }=\left\{\begin{array}{cc}{\frac{S_{T}}{S_{T}d_{o}} u_{c o l d}} & {\text { for } S_{D} \geq \frac{S_{T}+d_{o}}{2}} \\ {\frac{S_{T}}{2\left(S_{D}d_{o}\right)} u_{c o l d}} & {\text { for } S_{D}<\frac{S_{T}+d_{o}}{2}}\end{array}\right.$ (16)
Eq. 17 is the diagonal distance between the centers of adjacent tubes, also called diagonal pitch (SD).
$S_{D} \equiv \sqrt{S_{L}^{2}+\left(\frac{S_{T}}{2}\right)^{2}}$ (17)
where ST[m] and SL [m] are, respectively, the transversal and longitudinal pitch.
2.3 Parametric analysis
The objective for the case studies is the minimization of the global pressure drop and the minimization of volume of materials in the boiler through a Constructal Law approach.
To determine the optimal Constructal Law design of the boiler, a parametric analysis was performed on the main geometrical parameters of the boiler.
In particular:
1. The tubes layout was varied into: staggered arrangement (as in the existing design, in Figure 2) or aligned arrangement (as in Figure 3);
2. The number of tubes (N) is varied between five discrete values: 8, 12, 16, 20 and 24;
3. The diameter of the tubes (50 mm) is varied between 30 and 150 mm.
The boiler has been designed with uniform element tube spacing and straight tubes as shown.
The geometrical models for the heat exchanger used can be referred to two different arrangements of the bank tubes. A tube bank usually consists of many rows of tubes both parallel and perpendicular to the flow direction. The tubes may be arranged in staggered or aligned configurations as shown in (Figure 2) and (Figure 3), respectively.
The boiler control volume is shown in Figure 2 and Figure 3 including also geometrical nomenclature (the length L, the height H and the width W).
The aligned configuration was modeled by integrating Eq.18 in substitution of Eq. 16 for the staggered configuration.
$u_{\max }=\frac{S_{T}}{S_{T}d_{o}} u_{c o l d}$ (18)
Similar formulations of the dimensionless parameters (C1,C2, m, n, f and c) are implemented for the aligned configuration in [19].
The thickness of the pipes is considered constant for all simulation, indicated with twall, equal to 2,5 mm, as in the real case.
Figure 2. Basic mathematical model of boiler (real case)
Figure 3. Basic mathematical model of boiler (Case 7)
Initially ten cases were defined, characterized by increasing values of number of tubes and different geometric arrangements. Table 2 resumes the parametric features of 10 cases, including the existing case
These cases have been chosen to not alter the boundary characteristics of the original boiler, in fact, the height, the depth and the width do not vary throughout the parametric analysis.
Applying a variation to the number of tubes without changing the size of the boiler, a maximum variation in the diameter of the tubes was fixed so as to allow a section of sufficient passage of water equal to 10% of the distance between the centers of the tubes, as described in [18].
Table 2. Design parameter
Case study 
Geometrical configuration 
Number of tubes 
Real case 
Staggered 
16 
#1 
Staggered 
8 
#2 
Staggered 
12 
#3 
Staggered 
20 
#4 
Staggered 
24 
#5 
Inline 
8 
#6 
Inline 
12 
#7 
Inline 
16 
#8 
Inline 
20 
#9 
Inline 
24 
The analysis of the results is divided in two contributions at first: pressure losses and heat transfer analyses; secondly, a comparison of the results with optimal diameter analysis and overall performance coefficient is discussed.
3.1 Pressure drops
The first stage concerns the study of overall pressure drops in both the tubes and the shell. Figure 4 and 5 show the trends of pressure drops generated, respectively in the shellside and in the tubeside, with increasing inner diameter.
In order to clarify how the diameter growth affects the two arrangements (aligned and staggered), in Figure 4 and Figure 5 pressure drops are represented with a logarithmic scale. The graphs show, respectively, the pressure drop that occur internally both in the cold and hot side for the 10 cases.
Figure 4. Pressure drops – shell side
Thanks to the smaller number of deviations to which the fluid is subject, the crossflow configuration with the aligned tubes generates lower pressure drops in the shell side Overall pressure drops (Eq. 19) is the total of the pressure drops generated inside and outside the bank of tubes.compared to those generated by the heat exchanger with staggered tubes. For example, with equal numbers of tubes,DPoverall = DPcold + DPhot(19)the real case generates pressure drops greater than 2x103 MPa if compared to Case 7; then, when the diameter is approximately 100mm, the curves are very close. The highest values are traced for Case 3 and 4 that are only to a limited extent higher than the real case. The lowest results are reported for case 5 (aligned with 8 tubes).
Figure 5. Pressure drop – flow inside tube
Figure 5 describes instead the pressure drops that are generated in the tubes while increasing their number, from 8 in Case 1 and 5 to Case 4 and 9 with 24 tubes. Generally, pressure drops decrease with the increase of the number of pipes. The results depend strongly on the number of tubes; the trends of pressure drops of the aligned cases are comparable to those of the staggered cases.
Figure 6. Overall pressure drops Figure 6 shows the trends of overall pressure drops.
Overall pressure drops (Eq. 19) is the total of the pressure drops generated inside and outside the bank of tubes.
$\Delta P_{\text {overall }}=\Delta P_{\text {cold }}+\Delta P_{\text {hot }}$ (19)
The real case presents the global minimum of overall pressure drops close to the actual diameter of 50 mm. In cases 3 and 4, however, the lower value is for diameter values close to 3040 mm. Other cases have the same trend; but they have the minimum in the vicinity of larger diameters.
3.2 Thermal analysis
The second analysis concerns the study of the heat exchanged by the boiler as function of the tubes’ diameter.
From the point of view of the sizing of the boiler, it is necessary that the machine is able to meet the required thermal load. For this reason, Figure 7 shows heat power as function of the diameter of the tubes.
Typically, the convection coefficient of a row increases with increasing row number, after which there is little change in flow conditions and hence in the convection coefficient. For large SL, the influence of upstream rows decreases, and heat transfer in the downstream rows is not enhanced. For this reason, operation of aligned tube banks with ST/SL < 0.7 is undesirable [18].
Figure 7. Thermal Power transfer
The dotted line called "Ref" shows the thermal power required by the user (46 kW). Its representation is needed to understand if and where the trend of thermal power developed in each case reaches the required thermal power As can be observed from Figure 7, the trend of the thermal power for the various cases, intersects the reference line for different diameters. The case, requiring the least diameter to cover the thermal load, is the case 4. Case 9, also, required a diameter equal to 78mm, but the optimum diameter exceeds the threshold diameter by failing to ensure the normal course of the cold fluid.
Table 3. Summary of the optimal data for the analyzed configurations
Real Case 
Case #1 
Case #2 
Case #3 
Case #4 
Case #5 
Case #6 
Case #7 
Case #8 
Case #9 

di 
[mm] 
50 
96 
73 
38 
37 
104 
98 
91 
84 
78 
do 
[mm] 
55 
101 
78 
43 
42 
109 
103 
96 
89 
83 
Number of tubes 
 
16 
8 
12 
20 
24 
8 
12 
16 
20 
24 
Pressure Drop shell side 
105[MPa] 
119.70 
302.48 
143.14 
119.75 
90.26 
1256.49 
352.53 
76.11 
21.09 
6.71 
Pressure Drop tubeside 
105[MPa] 
60.33 
4.59 
12.35 
188.40 
185.06 
3.12 
3.00 
3.41 
4.18 
5.16 
Overall pressure drop 
105[MPa] 
180.02 
307.07 
155.49 
308.15 
275.33 
1259.61 
355.53 
79.52 
25.27 
11.87 
Thermal Power 
[kW] 
46.00 
46.16 
46.13 
46.06 
46.08 
46.55 
46.24 
46.17 
46.01 
46.08 
Epsilon 
% 
88.53 
88.84 
88.77 
88.65 
88.68 
89.59 
89.00 
88.85 
88.55 
88.68 
Exchange surface 
[m2] 
1.66 
1.52 
1.76 
1.62 
1.90 
1.64 
2.33 
2.90 
3.36 
3.75 
3.3 Optimal inner diameter
The optimal diameter was calculated for each configuration on the basis of the required heat load (46 kW).
These graphs show that although case 5 and case 6 are characterized by small pressure drops both in the shell side and in the tube side, they are able to cover the thermal user requirements only with diameters larger than 90 mm. This would cause very large overall pressure drops as in Figure 6.
If the diameter of the tubes are considered constant, the resulting pressure drops and thermal power are reported in Table 3.For the 10 cases analyzed it was determined that the optimal diameter allows to obtain the pressure drops and the thermal power required as shown in Table 4.
The geometry, the number and arrangements of longitudinal pipes affect both the pressure drops and the amount of heat exchanged in the boiler.
Table 4. Results for constant inner diameter (50 mm)
Case 
di 
Overall pressure drop 
Thermal Power 
[m] 
105[MPa] 
[kW] 

Real Case 
0.05 
180.02 
46.00 
Case #1 
0.05 
106.17 
29.11 
Case #2 
0.05 
93.12 
41.61 
Case #3 
0.05 
381.24 
47.50 
Case #4 
0.05 
372.55 
47.81 
Case #5 
0.05 
105.04 
9.87 
Case #6 
0.05 
75.95 
16.24 
Case #7 
0.05 
60.36 
22.78 
Case #8 
0.05 
50.55 
28.37 
Case #9 
0.05 
43.78 
33.75 
The optimal diameter for cases 3 and 4, is equal to 38 and 37 mm respectively, but the pressure drops which are generated are higher than those that are generated in the real case. Increasing the internal diameter would allow for further nominal heat output of the machine and consequently lower drops of tube side pressure.
In any case, the actual boiler has pressure drops in line with the average of the cases.The best cases are the 2, 7, 8 and 9 which have pressure drops respectively equal 14%, 54%, 86% and 93% less than real case.
3.4 Overall performance coefficient application
In [16] a concept of performance analysis determined by the contribution of the pressure loss and heat removed, applied to fins in heat exchanger, was developed .
In this work, the Overall Performance Coefficient was applied to quantify the relative relevance of pressure drops and heat exchanged in the parametric analysis proposed.
The overall performance coefficient is calculated as in Eq.19:
$P_{1}^{*}=\alpha \tilde{Q}_{i}+(1\alpha) \frac{1}{\Delta P_{i}^{*}}$ (20)
where Qi is the dimensionless heat transfer in the ith case, calculated as ratio between heat transfer in the ith case and heat transfer in the real case [16];while DP* is a mean dimensionless relative loss of pressure in the ith case respect real case and a is called relevance, varying between 0 and 1, representing the weight of heat transfer maximization with respect to loss of pressure minimization, depending on the context and on the use of the heat exchanger. The results are shown in Figure 8.
Because of the pressure drop detected in the cases 7, 8 and 9 smaller than case of reference, it has opted for a graphical representation that does not show these trends for entire in Figure 8, thus obtaining a clearer understanding of the results. The results report that, for optimal diameter, the trend of cases is not constant. The Table 4 shows how the arrangement of the pipes inside the heat exchanger results in the variation of the pressure drops and the efficiency of the boiler. If a is 0, the overall performance coefficient, for the case 7, results 2.26, for the case 8 and case 9 are, respectively,
7.12 and 15.17; others are approximately equal to the reference case. If in the first part of the graph, the curves differ from each other, instead if a is 1, the trends converge towards the real case, in fact, the first prerequisite of the paper is the respect of the thermal power by varying the geometry. All curves converge on the right side for a=1 on the same point, that corresponds to the required thermal power of the machine, while the best performing solution would seem Case 9 (Aligned with 24 tubes), or in other words the scenario with the lowest pressure drops.
In fact, the best case is the minor overall pressure drops to cover thermal power required.
Figure 8. Overall performance coefficient in function of heat transfer
Figure 9. Overall performance (material) coefficient in function of volume of tubes
However, since the base assumption of the analysis is that the prerequisite of the pellet burner is to be able to cover and guarantee 46 kW of heat power, it makes sense to also re interpret the Overall performance coefficient and to better adapt it to the context. In fact, if the heat power achievable by the boiler is fixed and the idea of reducing the amount of materials – and indirectly the amount of energy used to produce them  needed to achieve the same results is introduced, it is possible to state that Overall performance coefficient in the second scenario can be calculated as in the following Eq.20:
$P_{2}^{*}=\beta \frac{1}{\tilde{V}_{i}}+(1\beta) \frac{1}{\Delta P_{i}}$ (21)
where Vi, is the dimensionless volume of tubes in the ith case, calculated as ratio between volume of tubes in the ith case and volume of tubes in the real case [16] and b can be defined relevance of materials, varying between 0 and 1,representing the weight of volume of materials minimization with respect to loss of pressure minimization.
In particular, if b is set equal to 1 this means that, for the technical use requested, P * is equal to the inverse ratio of dimensionless volume, if b is set equal to 0 this means that P *is equal to the dimensionless inverse ratio of pressure drops. Each value of relevance between 0 and 1 will describe intermediate needs.
It is desirable under these assumptions, that P *is as high as possible.
Results mark in this case, a heterogeneous scenario with diverging results. There appears not to be a best option from all points of view but the optimal solution needs to be identified through tradeoffs. In particular there isn’t a better global case respect to the real case, in fact minor pressure drops are connected to major volume of material of tubes. From this point of view it appears that the real case could be considered a good tradeoff between material use and pressure drops. In fact, Case 2 has values close to those of the real case in terms of both volume and pressure drops, presenting lower pressure drop and greater quantities of material used, instead Cases 1 and 3 have reversed trends. Cases 4, 5 and 6, have higher pressure drops and volume of materials employed.
The paper reports a parametric analysis under multiple points of view considering the Costructal Law as main inspiration for the design variation.
The parametric analysis considers the variation of most geometrical parameters of the boiler, including the diameter size, the positioning and the numbering of the tubes, while keeping the boiler volume constrained to the existing machine.
Results are presented in terms of heat power variation and overall pressure drops in both the shell and tubes sides. The application of the Overall performance coefficient leads to a comprehensive perspective on the results, identifying that the best performing solution is the aligned configuration with 24 tubes; but if volume of materials of tubes is introduced as variable, the results are more heterogeneous, because there isn’t an overall best case if compared to the real case and the best solution can be identified only through tradeoffs.
The analysis proposed in the paper marks one of the first applications to boilers of Constructal law and introduces the material uses perspective as a parameter to achieve a design optimization of the system under examination.
It is particularly needed to take care in the selection of relevance used in the overall performance coefficient method since the results are heavily influenced by it.
Moreover, the paper marks the introduction of an eco design inspired optimization in the Overall performance coefficient calculation as one of the two optimization pillars in the methodology.
A 
Heat transfer surface (m2) 
C 
Correction factor for Nusselt number 
cp 
Specific heat capacity at constant pressure (J kg−1 k−1) 
d 
Diameter (m) 
F 
LMTD correction factor 
f 
Friction coefficient () 
h 
Convective heat transfer coefficient for one fluid (W m−2 K−1) 
H 
Height 
k 
Thermal conductivity (W m−1 K−1) 
L 
Length of the heat exchanger channel (m) 
M 
Mass flow rate (kg s−1) 
N 
Number of the tubes 
Nu 
Nusselt number 
P* 1 
Overall performance coefficient function of heat transfer 
P* 2 
Overall performance coefficient function of materials 
PL 
Dimensionless longitudinal pintch 
Pr 
Prandtl number 
PT 
Dimensionless trasfersal pintch 
Q 
Heat flow (kW) 
R" f 
Fouling factor (m2 K/W) 
Re 
Reynolds number 
SD 
Diagonal pintch (m) 
SL 
Longitudinal pintch (m) 
SL 
Trasfersal pintch (m) 
T 
Temperature (K) 
U 
Overall heat transfer coefficient (W m−2 K−1) 
u 
Velocity of flow (m s−1) 
W 
Width (m) 
ΔP 
Pressure drop (MPa) 
ΔPoverall 
Overall pressure drop (MPa) 
ΔTlm 
Mean logtemperature difference (K) 
Greek symbols 

a 
Relative 
b 
Relative of materials 
m 
Viscosity (kg m−1 s−1) 
mw 
Viscosity in wall temperature (kg m−1 s−1) 
r 
Density of the fluid (kg m−3) 
Subscripts 

cold Cold fluid 
[1] Manjunath K., Kaushik S.C. (2014). Second law thermodynamic study of heat exchangers: a review,Renewable and Sustainable Energy Reviews, Vol. 40,pp. 348374.
[2] Bejan A. (1997). Constructaltheory network of conducting paths for cooling a heat generating volume,International Journal of Heat and Mass Transfer, Vol.40, pp. 799816.
[3] Bejan A., Lorente S. (2004). The constructal law and the thermodynamics of flow systems with configuration, International Journal of Heat and Mass Transfer, Vol. 47, pp. 32033214.
[4] Chen L., Feng H., Xie Z., Sun F. (2015). Thermal efficiency maximization for H and Xshaped heat exchangers based on constructal theory, Applied Thermal Engineering, Vol. 91, pp. 456462.
[5] Yang J., Fan A., Liu W., Jacobi A.M. (2014).Optimization of shellandtube heat exchangers conforming to TEMA standards with designs motivated by constructal theory, Energy Conversion and Management, Vol. 78, pp. 468476.
[6] Manjunath K., Kaushik S.C. (2014). The second law analysis of an unbalanced constructal heat exchanger,International Journal of Green Energy, Vol. 11, pp.173192.
[7] Manjunath K., Kaushik S.C. (2014). Entropy generation and thermoeconomic analysis of constructal heat exchanger, Heat Transfer  Asian Research, Vol. 43, pp. 3960.
[8] Zimparov V.D., da Silva A.K., Bejan A. (2006).Constructal treeshaped parallel flow heat exchangers,International Journal of Heat and Mass Transfer, Vol.49, pp. 45584566.
[9] Da Silva A.K., Bejan A. (2006). Dendritic counterflow heat exchanger experiments, International Journal of Thermal Sciences, Vol. 45, pp. 860869.
[10] Da Silva A.K., Lorente S., Bejan A. (2004).Constructal multiscale treeshaped heat exchangers,Journal of Applied Physics, Vol. 96, pp. 17091718.
[11] Lorenzini G., Rocha L.A.O. (2009). Geometric optimization of TYshaped cavity according to Constructal design, International Journal of Heat and Mass Transfer, Vol. 52, pp. 46834688.
[12] A. Bejan (2002). Dendritic constructal heat exchanger with smallscale crossflows and largerscales counterflows, International Journal of Heat and Mass Transfer, Vol. 45, pp. 46074620.
[13] Liu X., Chen L., Feng H., Qin X., Sun F. (2016).Constructal design of a blast furnace ironmaking process based on multiobjective optimization, Energy,Vol. 109, pp. 137151.
[14] Kang D.H., Lorente S., Bejan A. (2010). Constructal architecture for heating a stream by convection,International Journal of Heat and Mass Transfer, Vol.53, pp. 22482255.
[15] Koonsrisuk A., Lorente S., Bejan A. (2010).Constructal solar chimney configuration, International Journal of Heat and Mass Transfer, Vol. 53, pp. 327333.
[16] Lorenzini G., Moretti S. (2009). A. Bejan’s constructal theory approach to the overall optimization of heat exchanging finned modules with air in forced convection and laminar flow condition, Journal of Heat Transfer, Vol. 131, p. 081801.
[17] Khan W.A., Culham J.R., Yovanovich M.M. (2006).Convection heat transfer from tube banks in crossflow:Analytical approach, International Journal of Heat and Mass Transfer, Vol. 49, pp. 48314838.
[18] Incropera F.P., DeWitt D., Bergman T. (2011).Fundamentals of heat and mass transfer, J.W. Sons,Ed., seventh ed., pp. 468476,517593.
[19] Žukauskas A., Ulinskas R. (1985). Efficiency parameters for heat transfer in tube banks, Heat Transfer Engineering, Vol. 6, pp. 1925.