Next Article in Journal
Research on Estimation Method of Fuel Cell Health State Based on Lumped Parameter Model
Previous Article in Journal
Exhaust Emissions and Energy Consumption Analysis of Conventional, Hybrid, and Electric Vehicles in Real Driving Cycles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Nonlinear Inverse Design Problem for a Pipe Type Heat Exchanger Equipped with Internal Z-Shape Lateral Fins and Ribs

Department of Systems and Naval Mechatronic Engineering, National Cheng Kung University, Tainan 701, Taiwan
*
Author to whom correspondence should be addressed.
Energies 2020, 13(23), 6424; https://doi.org/10.3390/en13236424
Submission received: 12 November 2020 / Revised: 30 November 2020 / Accepted: 2 December 2020 / Published: 4 December 2020

Abstract

:
A non-linear three-dimensional inverse shape design problem was investigated for a pipe type heat exchanger to estimate the design variables of continuous lateral ribs on internal Z-shape lateral fins for maximum thermal performance factor η. The design variables were considered as the positions, heights, and number of ribs while the physical properties of air were considered as a polynomial function of temperature; this makes the problem non-linear. The direct problem was solved using software package CFD-ACE+, and the Levenberg–Marquardt method (LMM) was utilized as the optimization tool because it has been proven to be a powerful algorithm for solving inverse problems. Z-shape lateral fins were found to be the best thermal performance among Z-shape, S-shape, and V-shape lateral fins. The objective of this study was to include continuous lateral ribs to Z-shape lateral fins to further improve η. Firstly, the numerical solutions of direct problem were solved using both polynomial and constant air properties and then compared with the corrected solutions to verify the necessity for using polynomial air properties. Then, four design cases, A, B, C and D, based on various design variables were conducted numerically, and the resultant η values were computed and compared. The results revealed that considering continuous lateral ribs on the surface of Z-shape lateral fins can indeed improve η value at the design working condition Re = 5000. η values of designs A, B and C were approximately 13% higher than that for Z-shape lateral fins, however, when the rib numbers were increased, i.e., design D, the value of η became only 11.5 % higher. This implies that more ribs will not guarantee higher η value.

1. Introduction

It is well known that if a furnace produced industrial waste heat that can be recycled by utilizing economizers or regenerators before it is discharged into the environment, the thermal efficiency of the system can thus be improved significantly. In addition, if the exhaust gas can be cooled down to a reasonably low temperature before its discharge then this will also be environmentally friendly.
Economizers and regenerators are basically the so-called pipe type heat exchangers, and internally longitudinal finned tubes are usually utilized to effectively enhance heat transfer for the tube side. They can be used in many practical industrial applications such as petroleum industries, boiler systems, power plants and chemical engineering industries.
Many researchers have investigated the research topic regarding internally longitudinal finned tubes problems. For instance, Webb [1] has discussed internally finned tubes in detail, and Fabbri [2,3] examined an optimal design problem in determining the lateral profile of longitudinal fins located in cylindrical tubes.
The pressure drop behaviors inside circular finned tubes under turbulent flow conditions were examined experimentally by El-Sayed et al. [4]. They have reported that, in the periodic fully developed region, the tube pressure drop with continuous fins is higher than that of inline arrangement fins and lower than that of staggered arrangement fins. The heat transfer performance of vertical internally finned tube under forced convection condition were investigated numerically by Al-Sarkh and Abu-Nada [5]. They have found the existence of an optimum combination of fin numbers and height to yield the best heat transfer performance.
Heat transfer performance of a three-dimensional heat exchanger with internally longitudinally finned tube with blocked core tube and stream wise wavy fin was examined numerically by Wang et al. [6]. The numerical solutions were verified by comparing them with the experimental data and the effects of both wave height and wave distance on heat transfer performance were investigated. The numerical results indicated that the Nusselt number and friction factor increase with the increase in the wave height, while they decrease with the increase in the wave distance. Yu and Tao [7] investigated the pressure drop and heat transfer characteristics of turbulent flow in annular tubes with internal wave-like longitudinal fins. Results revealed that heat transfer can be enhanced using wave-like finned tubes and wave number 20 can result in best performance.
Three different lateral fin profiles, i.e., S-shape, Z-shape and V-shape, of internally finned tubes with blocked core-tube was investigated and the heat transfer performances of those lateral fin profiles were examined by Wang et al. [8]. The hydro-thermal performances of the studied lateral fin profiles were computed and compared under three working conditions. It was found that Z-shape fins have the best performance.
Based on the research findings in Wang et al. [8], Duan et al. [9] extended the concept of Z-shape fins to construct the blossom shape internal fins and utilized in a double-tube structure internally finned tube. The hydro-thermal behaviors of the heat exchanger were analyzed, and the conclusion was that the blossom shape fin is more suitable for the exhaust gas heat recovery system than the wave-like fin.
Rib-roughened cooling passages are commonly used in heat exchangers; the utilization of ribs on the fin surface of heat exchanger introduces two heat transfer enhancement characteristics, i.e., an increase in heat transfer area and a significant increase in heat transfer coefficients.
Much research regarding this topic can be found in the open literature. For example, Chai et al. [10] numerically investigated the thermal performance of a microchannel heat sink with rectangular inserted ribs. They have considered the designs with different rib lengths and widths and arrangements, and concluded that when Re < 600, the interrupted microchannel with ribs was a good passive method for heat transfer enhancement. Xia et al. [11] studied the micro heat sink equipped with fan-shaped reentrant cavities and internal ribs. Result indicated that the effect of rib height was stronger than the individual effect of the size or arrangement of the reentrant cavity for Re > 300.
Ahmed [12] investigated the optimal design of a ribbed flat-plate fin heat sink. In this work, the sizes, positions, numbers, and orientations of ribs were considered as the design variables and the objective was to obtain the optimal thermal performance heat sink. Results indicated that a ribbed plate fin heat sink provides thermal performance of 1.55 times greater than a plate fin heat sink under the studied conditions.
Based on the above preliminary reviews, it is expected that the thermal performance factor of a pipe type heat exchanger can be improved by including continuous lateral ribs to internal Z-shape lateral fins. The objective of this shape design problem is to determine the optimal sizes and number of continuous lateral ribs to yield the optimal thermal performance factor under fixed fin volume constraints with an efficient optimization tool.
The Levenberg–Marquardt method (LMM) [13] has been considered as an efficient optimization tool for shape design problems. For instance, Chen and Huang [14] used LMM and commercial code SHIPFLOW to estimate the optimal shape of the after hull based on the desired wake distribution. Huang and Li [15] determined the optimal boundary geometry of a three-dimensional conductive body for producing boundary isotherms with LMM. Huang and Chen [16] used LMM to estimate the optimal shape and perforation diameters of a perforated pin-fin array module. A three-dimensional wavy-shaped fin array design problem was considered by Huang and Tung [17], and the objective was to estimate the optimal wavy-shaped fin array with deformed sinusoidal function to minimize the average temperature of the base plate of fin array under a fixed fin volume constraint. Therefore, LMM was chosen in this study as the optimization algorithm.

2. Problem Formulation

A 3-D Pipe Type Heat Exchanger with blocked core-tube and internal Z-shape lateral fins with ribs (PTHE-ZR) was studied in the present inverse design problem to illustrate the methodology for estimating the optimal design variables of lateral ribs using LMM and CFD-ACE+ [18].
Figure 1 indicates the physical model of PTHE with internal Z-type lateral fins (PTHE-Z) listed with the parameter index. PTHE-Z is composed of an annulus tube and Z-shape lateral fins, the tube length is 400 mm, and there are 20 waves located in the cross-section of internally finned tube. The diameters are considered as do = 16 mm, Di = 32 mm, Do = 35 mm and δ= 0.2 mm, t = 0.1 mm, L = 7.6 mm, Lt = 400 mm, and N = 20. The inlet air temperature and wall temperature are taken as Tin = 373 K and Tw = 303 K, respectively [9], inlet pressure is considered as Pin = 0.2 MPa, and thermal conductivity of copper is kcopper = 387 W/m-K.
Here, Ω represents the computational domain and {Ω} = {Ω1 Ω2}, where Ω1 indicates the solid regions (tube and Z-shape lateral fin materials) and Ω2 is the working air flow region. The outer boundary surfaces of PTHE are subjected to constant wall temperature boundary conditions.
Due to the axial symmetry of PTHE-Z, for simplicity of computations, the computational domain is chosen as 1/40 of the original PTHE-Z and its 3-D computational domain and grid system are given in Figure 2a,b, respectively.
The conduction equation for the 3-D solid domain Ω1 is given by:
[ 2 T s ( Ω 1 ) x 2 + 2 T s ( Ω 1 ) y 2 + 2 T s ( Ω 1 ) z 2 ] = 0 ;   in   Ω 1  
where Ts represents the solid, i.e., tube and fin, material temperature distributions in Ω1.
It is assumed that air is a 3-D steady-state incompressible flow in region Ω2 and the thermal properties of air are assumed to be a function of temperature; this makes the problem non-linear. In addition, due to force convection and moderate air temperature conditions, buoyancy and radiative heat transfer effects are both neglected.
The 3-D equations of continuity, momentum, and energy in the steady-state turbulent flow with the standard k–ε model, turbulent kinetic energy and turbulent energy dissipation rate are illustrated in Equations (2)–(6), respectively [19]:
ρ u i x i   =   0
ρ u j u i x j   =   p x i + x j [ μ t ( u i x j + u j x i ) ]
ρ u j T a x j   =   x j [ ( μ l σ l + μ t σ t ) T a x j ]
ρ u j k x j   =   x j ( μ t σ k k x j ) + μ t ( u i x j + u j x i ) u i x j ρ ε
ρ u j ε x j   =   x j ( μ t σ ε ε x j ) + C 1 μ t ε k ( u i x j + u j x i ) u i x j C 2 ρ ε 2 k
where μl and μt are the kinematic viscosity and turbulent viscosity, respectively, and σl and σt stand for the kinematic Schmidt number and turbulent Schmidt number, respectively, with μt = Cμ(k2/ε), C1 = 1.44, C2 = 1.92, Cμ = 0.09, σk = 1.0 and σε = 1.3. Ta represents air temperature distribution in Ω2. It is assumed a non-slip boundary condition, i.e., ui = 0, applied on all pipe surfaces. The 3-D hydro-thermal conjugate problem in domain Ω of PTHE-Z is solved using CFD-ACE+, and the physical properties of air are assumed as a function of air temperature Ta.
The functions of density ρ(Ta), specific heat Cp(Ta), viscosity μl(Ta) and conductivity ka(Ta) of air are given below [9]:
ρ ( T a ) = 6 × 10 6 × ( T a 273 ) 2 0.0058 × ( T a 273 ) + 2.4028   ( k g / m 3 )
C p ( T a ) = 3 × 10 9 × ( T a 273 ) 3 + 2 × 10 6 × ( T a 273 ) 2   0.0001 × ( T a 273 ) + 1.004   ( k J / k g K )
μ l ( T a ) = 1 × 10 6 × ( T a 273 ) 3 0.0007 × ( T a 273 ) 2   + 0.1679 × ( T a 273 ) + 10.282   k g / m s
k a ( T a ) = 3 × 10 6 × ( T a 273 ) 2 + 0.0084 × ( T a 273 )   + 2.3802   W / m K
Equations (7)–(10) indicate that the values of thermal properties depend on the value of Ta, and this makes the problem nonlinear. Equations (1)–(10) are utilized to calculate the pressure, velocity and temperature distributions for PTHE-Z when all the boundary conditions and the size of the internal Z-shape lateral fins are known. Thereafter, the average outlet air temperature Tout can be calculated.
The Reynolds number Re, Nusselt number Nu, and friction factor f are defined as follows:
Re = ρ u i n D h μ l
N u = h D h k a
f = Δ P ( ρ u 2 / 2 ) ( L t / D h )
where u i n and Δp are the inlet velocity and total pressure drop, respectively, Dh denotes hydraulic diameter and is defined clearly in [8], and Lt is the tube length with internal fins. In addition, the heat transfer coefficient h is defined as:
h = m ˙ C p ( T in T out ) A Δ T LMTD
In Equation (14), A denotes the total heat transfer area, and ΔTLMTD represents the logarithmic mean temperature difference and is defined as:
Δ T LMTD = ( T in T w ) ( T out T w ) ln ( T in T w T out T w )
Here, Tin and Tout represent the inlet and average outlet air temperatures of heat exchanger, respectively, and Tw is the tube wall temperature. Based on the constant pumping power constraint, the thermal performance factor η of PTHE-ZR is defined as the ratio between the convective heat transfer coefficients of PTHE-ZR, hZR, and PTHE-Z, hZ. Finally, η of PTHE-ZR can be obtained below [20,21]:
η = h ZR h Z p p = N u ZR N u Z p p = N u ZR N u Z f ZR f Z 1 3
Under identical power conditions, η 1 indicates that the increase in the enhancement ratio of heat transfer rate (NuZR/NuZ) is larger than the increase in friction losses (fZR/fZ)−1/3. A higher η will result in better heat dissipation performance under the same energy consumption conditions.

3. The Inverse Design Problem of Ribs

Wang et al. [8] utilized a k–ε turbulence model with a wall function method to investigate the hydro-thermal performance of PTHE with internal lateral fins and blocked core-tube. PTHE with S-shape, Z-shape, and V-shape internally lateral fin profiles, i.e., PTHE-S, PTHE-Z and PTHE-V, were investigated and compared. They have concluded that PTHE-Z achieved the best hydro-thermal performance.
It is well known that roughness elements such as ribs can be used on heat exchanger surfaces to promote turbulence and enhance convective heat transfer. These various types of ribs are mounted on fin surfaces and have been widely used in engineering applications to enhance heat transfer of heat exchangers, mixing chambers, turbine blade cooling and electronic equipment cooling.
The design concept of this study was to include some pairs of rectangular lateral ribs on the surface of Z-shape fin to improve the thermal performance factor η of PTHE-ZR.
The proper design variables need be identified prior to the design process, and only effective design variables need to be included to avoid the long computational time of the design process. The following four design cases will be considered and investigated.

3.1. Design A

In design A, it is assumed that there is one pair of rib and the position W and height H of rib are considered as the design variables. Figure 3a–c indicates the geometry of the isometric view, cross-sectional view, and computational domain of PTHE-ZR, respectively. The design variables W and H can be rearranged in the following expression:
Θ =   { Θ 1 , Θ 2 } = { W , H }
Figure 3 implies that the fin position and height of upper and lower ribs are identical. In addition, the volume of each lateral fin with ribs, Vol, is obtainable using the following equation:
Vol   =   ( L × δ + 2 × H × t ) × L t
Once H is determined, the thickness of lateral fin δ can be obtained using volume constraint condition Equation (18) and the geometry of Figure 3c can be constructed.

3.2. Design B

In design B, it is also assumed that there is one pair of ribs, however, the position W1 and the height H1 of the upper rib, and the position W2 and the height H2 of the lower rib are considered as the design variables. Figure 4a–c indicates the geometry of the isometric view, cross-sectional view, and computational domain of PTHE-ZR, respectively. The design variables W1, W2, H1 and H2 are arranged as follows:
Θ =   { Θ 1 , Θ 2 , Θ 3 , Θ 4 } = { W 1 , H 1 , W 2 , H 2 }
Figure 4c indicates that there are four design variables in design B. The volume of each lateral fin with ribs, Vol, can be computed with the equation below:
Vol   =   [ L × δ + H 1 + H 2 × t ] × L t
Once H1 and H2 are obtained, the thickness of the lateral fin δ can be computed from Equation (20) and the geometry of PTHE-ZR can then be plotted.

3.3. Design C

In design C, two pairs of ribs are considered and Figure 5a–c plots the geometry of the isometric view, cross-sectional view, and computational domain of the PTHE-ZR in design C, respectively. It can be learned from Figure 5c that the positions W1 and W2 and the heights H1 and H2 of ribs are considered as the design variables, i.e., there are also four design variables in design C.
The design variables W1, W2, H1 and H2 are arranged as follows:
Θ =   { Θ 1 , Θ 2 , Θ 3 , Θ 4 } = { W 1 , H 1 , W 2 , H 2 }
The volume of each lateral fin with ribs, Vol, can be calculated using the following expression:
Vol   =   [ L × δ + 2 × H 1 + H 2 × t ] × L t
Once H1 and H2 are obtained, the thickness of the lateral fin δ can be obtained using Equation (22) and the geometry of PTHE-ZR in design C can be plotted.

3.4. Design D

Three pairs of ribs are considered in design D, and Figure 6a–c illustrates the geometry of the isometric view, cross-sectional view, and computational domain of PTHE-ZR in design D, respectively. It is clear that the positions W1, W2, and W3 and the heights H1, H2, and H3 of ribs are considered as the design variables, i.e., there are six design variables in design D.
The design variables W1, W2, W3, H1, H2 and H3 can be arranged as follows:
Θ =   { Θ 1 , Θ 2 , Θ 3 , Θ 4 , Θ 5 , Θ 6 } = { W 1 , H 1 , W 2 , H 2 , W 3 , H 3 }
The volume of each lateral fin with ribs, Vol, can be calculated using the following expression:
Vol   =   [ L × δ + 2 × H 1 + H 2 + H 3 × t ] × L t
Once H1, H2 and H3 are estimated, the thickness of lateral fin δ can be obtained using Equation (24) and the geometry of PTHE-ZR in design D can be plotted.

3.5. The Cost Function

The purpose of this work was to obtain the optimal value of η(Θ) by determining the optimal shape of PTHE-ZR. Letting the desired η(Θ) of PTHE-ZR be denoted as ϕ, the fin shape design problem can be stated as follows: utilize ϕ to design PTHE-ZR for yielding optimal thermal performance factor η(Θ).
The cost function of the present shape design problem is as follows:
J [ Ω ( Θ ) ]   = [ η ( Θ ) ϕ ] 2 = A T A
Here, Ω denotes the computational domain of the PTHE-ZR and is a function of the design variable Θ. η denotes the estimated thermal performance factor of PTHE-ZR, and its value is calculated from the solution of the direct problem by using an updated geometry of PTHE-ZR. In addition, the volume constraints, Equations (18), (20), (22) and (24) for designs A, B, C and D, must be satisfied.

4. The Optimization Algorithm

The cost function is minimized with respect to the design variable Θ to obtain the following equation:
J [ Ω ( Θ ) ] Θ = η ( Θ ) Θ η ϕ = 0
Equation (26) is linearized by expanding η(Θ) in a Taylor series and retaining only the first-order terms. Thereafter, a weighting coefficient μn is included in the resultant expression to speed up the rate of convergence. Finally, the Levenberg–Marquardt method [13] can be obtained as:
( F   + μ n I ) Δ Θ = D
F = ΨTΨ
D = ΨTA
Δ Θ   = Θ n + 1 Θ n
Here, I denotes the identity matrix, the superscripts T and n represent the transpose matrix and iteration index, respectively, and Ψ is the Jacobian matrix and is given as:
Ψ   = η Θ T
The Jacobian matrix given in Equation (31) can be obtained by perturbing the design variables Θ one at a time and calculating the resultant change in the thermal performance factor η from the solution of direct problems.
To yield an expression for suitable iterative calculation, Equation (27) is rearranged and expressed as:
Θ n + 1   = Θ n + ( Ψ T Ψ + μ n I ) 1 Ψ T ( η ϕ )
In the beginning, the steepest-descent method (ηn ) is applied, then the value of μn decreases, and eventually Newton’s method (μn = 0) is utilized to compute the optimal solutions [13]. A sequence of direct problems is solved by CFD-ACE+ to update the design variables for PTHE-ZR by minimizing the difference between η and ϕ in this optimization algorithm.

5. Computational Steps

The iteration steps for obtaining the optimal design variables of this work are summarized below:
Step 1:
Use the original design variables Θ0 as the initial guesses to begin the calculations.
Step 2:
Solve the direct problems and then the estimated η can be obtained from Equation (16).
Step 3:
The Jacobian matrix is computed with Equation (31).
Step 4:
Update Θ n+1 from Equation (32).
Step 5:
If ( J n + 1 J n / J n + 1 ) < ε is not satisfied, return to Step 2 and repeat.

6. Results and Discussion

The optimal design algorithm with LMM was considered in the present study to estimate the design variables of the ribs for the internal Z-shape lateral fins in PTHE-ZR. To illustrate the algorithm of LMM in estimating the design variables of the lateral ribs by minimizing the cost function of PTHE-ZR, the optimal shape design problem is described below.
Copper with the thermal conductivity k = 387 W/(m − K) was chosen for PTHE-ZR materials and the working conditions were the same as those given by Wang et al. [8]. In addition, the temperature-dependent physical properties of air were identical to those reported by Duan et al. [9].
It is important to verify the accuracy of the computed numerical data because it plays a critical role in the shape design problem. If the solutions obtained by CFD-ACE+ are not accurate enough, the optimal rib dimensions of PTHE-ZR cannot be designed properly and correctly. The following benchmark model is utilized for the above-mentioned purpose, and the numerical results will be compared with the corrected results given by Wang et al. [8].
The geometry for PTHE-Z is illustrated in Figure 1 and is considered as the initial design of this study. Based on Equation (16), the thermal performance factor of PTHE-Z is taken as η = 1.0 and the following dimensions are adopted:
  • do = 16 mm, Di = 32 mm, Do = 35 mm, Tw = 303 K, Tin = 373 K,
  • Pin = 0.2 MPa, kcopper = 387 W/m-K, δ = 0.2 mm, t = 0.1 mm,
  • L = 7.6 mm, Lt = 400 mm, N = 20
The volume of each internal Z-type lateral fin can be obtained as Vol = L × δ × Lt = 608 mm3. Additionally, temperature-dependent physical properties of air, Equations (7)–(10), were utilized.
The average outlet air temperature of PTHE-Z, Tout, was calculated utilizing CFD-ACE+ package with four different grid numbers 158,400, 323,400, 653,400 and 983,400 at Re = 5000, and the computed results were obtained as Tout = 304.1, 305.6, 306.8 and 307.5 K, respectively. The relative error of Tout between grid numbers 653,400 and 983,400 was calculated as 0.23%, which is very small; thus, grid number 653,400 was considered for all computations in this work.
In addition, in order to verify the necessity for using polynomial temperature-dependent physical properties of air, Tout and f using constant physical properties measured at T = 373 K (inlet temperature), 306 K (outlet temperature) and 340 K (averaged temperature), together with utilizing the ideal gas condition were computed with different Re values, and the results were then compared with those using polynomial physical properties and the corrected equations reported by Wang et al. [8].
Figure 7a,b illustrates the comparisons of Tout and f among six cases. Figure 7a indicates that at higher Re, the values of Tout obtained by using polynomial temperature-dependent physical properties of air matched best with corrected results and Figure 7b shows that the values of f obtained by using polynomial properties of air always best fit with corrected results. The validity of using temperature-dependent physical properties of air is thus verified. The computed temperature and pressure distributions of PTHE-Z at z = 100 mm are plotted in Figure 8a and Figure 9a, respectively.

6.1. Optimal Design of PTHE-ZR

The objective of the present study was to examine whether the use of lateral ribs can further improve η of the PTHE-Z. A 20% increase in η from PTHE-Z was required, i.e., φ = 1.2. After executing the design algorithm, it was difficult to reach this desired φ value; however, the optimal design variables could still be obtained for maximum η by using the LMM. The following four design cases were investigated in sequence to identify the most effective design pattern. The stopping criterion ε = 10−5 and volume of fin and ribs Vol = 608 mm3 were adopted in this study.

6.2. Design A

In design A, it is assumed that there is one pair of ribs and the upper and lower ribs have same position and height, i.e., the design variables are considered as W and H. The initial guesses were taken as W = 4.50 mm and H = 0.3 mm, by considering Re = 5000, φ = 1.2 and ε = 10−5, and after five iterations the estimated design variables were obtained as W = 5.16 mm and H = 0.516 mm, and the fin thickness reduced to η = 0.18641 mm. In addition, Nu, f and η were calculated as 148.4, 0.046 and 1.13, respectively, and are listed in Table 1. The computed temperature and pressure distributions at z = 100 mm are plotted in Figure 8b and Figure 9b, respectively.
It is clear from Figure 8a,b that by utilizing one pair of lateral ribs in PTHE-ZR, heat dissipation ability in PTHE-ZR is greatly improved than in PTHE-Z, as a result, temperature distribution in Figure 8b is observed lower than that in Figure 8a and therefore Nu of PTHE-ZR is obtained higher than that of PTHE-Z. It is also observed that the optimal position of the ribs of PTHE-ZR is located approximately at the position with the highest air temperature in the upper and lower zones of PTHE-Z. Utilizing intruded cool ribs into the hottest region of air will definitely reduce the average temperature of air effectively.
Due to the added surface area of the ribs, the pressure distribution of PTHE-ZR given in Figure 9b is found to be higher than that of PTHE-Z given in Figure 9a. Therefore, the computed friction f of PTHE-ZR is larger than that of PTHE-Z. However, when considering the objective function of this study, i.e., to maximize the thermal performance factor η, it is calculated as η = 1.13 with the present optimal design shape of ribs, i.e., η of PTHE-ZR is 13% higher than that of PTHE-Z under fixed fin volume constraint. It verifies that the design concept of this work is correct.

6.3. Design B

In design A, it was found that there is a significant increase in η by considering one pair of rib with identical position and height. In design B the upper and lower ribs can have different positions, W1 and W2, and heights, H1 and H2, and hope that η can be increased further.
The design variables are now considered as Θ = {W1, H1, W2, H2}. The optimal design variables of design A is used as the initial guesses of design B, i.e., W1o = W2o = 5.16 mm and H1 = H2 = 0.516 mm. By performing the LMM, the optimal solutions can be obtained as W1 = 5.01 mm, W2, = 5.44 mm, H1 = 0.50 mm, and H2 = 0.55 mm after four iterations, and the fin thickness is obtained as δ = 0.18623 mm. The resultant characteristics are obtained as Nu = 150, f = 0.0472 and η = 1.1321. The calculated temperature and pressure distributions at z = 100 mm are plotted in Figure 8c and Figure 9c, respectively, and the data are summarized in Table 1. It is clear from Figure 8b,c that temperature distributions of design B PTHE-ZR are lower than those of design A PTHE-ZR, due to the flexibility of the designing positions and heights of upper and lower ribs.
It is found that (1) W1 and W2 move to the left and right sides from the initial position, respectively, and (2) H1 and H2 become shorter and longer than the initial height, respectively, to yield better thermal performance factor η.
The reason for the first observation is that due to structure asymmetry, the hottest region of the lower passage is located left-hand-side of the upper passage. In design A, the positions of upper and lower ribs are identical, therefore the optimal position W is located approximately between the hottest regions of upper and lower passage. In design B, W1 and W2 are estimated by LMM separately to result in the best heat dissipation performance, and therefore W1 and W2 move to the hottest regions of lower and upper passages, respectively. In addition, the vertical distance of the hottest point of the upper passage to the internally lateral fin is greater than that of the lower passage, and as a result, H2 should be estimated higher than H1, and that is the explanation of the second observation.
In design B, Nu and f are both increased when compared with design A and the improvement of thermal performance factor η is not significant because it is calculated as η = 1.1321 with the estimated optimal design shape of ribs, i.e., η of design B PTHE-ZR is 13.21% higher than that of PTHE-Z under fixed fin volume constraint, but only 0.21% higher than that in design A.

6.4. Design C

Two pairs of ribs were considered in design C and the position and height of ribs in each pair were assumed to be identical. Figure 5c plots the geometry of the design, i.e., the design variables are taken as Θ = {W1, H1, W2, H2}.
The initial guesses of design C are assumed as Θ0 = {1.52, 0.30, 4.65, 0.3}. By performing the LMM, after four iterations the optimal solutions can be obtained as W1 = 3.90 mm, W2, = 5.54 mm, H1 = 0.31 mm, H2 = 0.57 mm and η = 0.17682 mm. The resultant characteristics are obtained as Nu = 153, f = 0.0498 and η = 1.1343. The calculated temperature and pressure distributions at z = 100 mm are plotted in Figure 8d and Figure 9d, respectively, and the data are summarized in Table 1.
Figure 8d indicates that the second pair of ribs is located at the original hottest region and the first pair of ribs is located at the resultant second hottest region to enhance heat dissipation performance, and the height of ribs in the second pair is higher than that in the first pair for the same reason stated previously. In design C, Nu and f are both increased when compared with design B, again the improvement of η is not significant because it is 13.43% higher than that of PTHE-Z under fixed fin volume constraint, but only 0.43% higher than that in design B.
The comparisons among PTHE-Z and PTHE-ZRs indicate that η of the PTHE-ZRs can be improved significantly, however, by varying rib positions and heights and increasing number of ribs from case A to B and to C, η is improved gradually but insignificantly from designs A to B and to C.
Can η be further increased by increasing number of pair of ribs? To answer this, the following design D needs be examined.

6.5. Design D

Three pairs of ribs were utilized in design D and the position and height of ribs in each pair were assumed to be identical. Figure 6c plots the geometry of the design, i.e., the design variables are taken as Θ = {W1, H1, W2, H2, W3, H3}.
The initial guesses of design D are assumed as Θ0 = {3.90, 0.31, 5.54, 0.57, 6.50, 0.57}. By performing the LMM, after three iterations, the optimal solutions can be obtained as W1 = 3.85 mm, W2 = 5.13 mm, W3 = 6.14 mm, H1 = 0.31, H2 = 0.56 mm, H3 = 0.55 mm, and δ = 0.16274 mm. The resultant characteristics are obtained as Nu = 154.3, f = 0.0539 and η = 1.1147. The calculated temperature and pressure distributions at z = 100 mm are plotted in Figure 8e and Figure 9e, respectively, and the data are summarized in Table 1.
Figure 8e and Figure 9e illustrate that utilizing three pairs of ribs can enhance heat dissipation but also increase pressure drop. The calculated thermal performance factor η is even smaller than that in design A, which implies that in design D the increase in the enhancement ratio of (NuZR/NuZ)(fZR/fZ)-1/3 is smaller than that in design A, and it can be concluded that more pairs of ribs cannot guarantee a higher thermal performance factor η.
The complete cross-sectional view of temperature and pressure distributions for various designs at z = 100 mm are illustrated in Figure 10. It is obvious from Figure 10 that temperatures are decreasing (i.e., Nu values are increasing) and pressures are increasing (i.e., f values are increasing) monotonically from PTHE-Z to design D. However, with regards to the concern of η, it shows that the η of design D is the smallest among four designs of PTHE-ZR and there is no significant difference of η among designs A, B and C. Due to manufacturing costs consideration, simpler fin structures are preferred, therefore design A is recommended for practical engineering applications.

6.6. Hydro-Thermal Characteristics of PTHE-ZRs at Re = 3000 and 7000

Next, the computed values of Nu, f and η of the above PTHE-Z and designed PTHE-ZRs with Re = 3000 and 7000 are summarized in Table 1. The variations in Nu, f and η with Re for various designs are plotted in Figure 11a–c, respectively. As expected, Nu and f both increase as rib number increases, and designs C and D have the largest and smallest η among PTHE-ZRs, respectively. Again, η for optimal designs A to D are all greater than 1.0 for the tested range of Re. This implies that when PTHE-ZRs are not working at the design condition Re = 5000, the hydro-thermal characteristics are still similar to those working at the design condition.
Finally, it is concluded that the validity of the optimal design for PTHE-ZR has been demonstrated, and this technology can readily be utilized in heat exchanger designs. In addition, wavy-shape fins may have better heat dissipation performance than straight fine, and the determination of optimal internal wavy-shape lateral fins can thus be investigated in the future.

7. Conclusions

An inverse shape design problem for pipe type heat exchangers was investigated in this work to determine the optimal shape of continuous lateral ribs on internal Z-shape lateral fins for maximum η with LMM and CFD-ACE+. The validity of using temperature-dependent polynomial function of air properties were verified by comparing the computed Tout and f with the corrected solutions [8], then four design cases at working condition Re = 5000 were considered based on various positions, heights and numbers of continuous lateral ribs. Results indicated that by considering lateral ribs on the surface of Z-shape lateral fins, Nu and f both increase as the number of ribs increase; however, this is not the case for η. The improvements of η in designs A to C PTHE-ZRs (i.e., with two to four ribs) is about 13% higher than that of PTHE-Z, however, is only 11.5% higher in design D. This implies that more ribs will not guarantee a higher η value. In addition, the values of η for designs A, B and C are approximately the same, due to the fact that simpler fin structures lead to lower manufacturing costs, therefore, design A is the best design for industrial applications. Finally, the hydro-thermal characteristics for off-design working conditions at Re = 3000 and 7000, which are not the design conditions, were considered and Nu, f and η ware computed. Results revealed that the trends of Nu, f and η at Re = 3000 and 7000 are similar to those working at the design condition Re = 5000. This indicates that the optimal shapes of continuous lateral ribs can be utilized at off-design conditions.

Author Contributions

Conceptualization, C.-H.H.; methodology, C.-H.H.; software, C.-Y.K.; validation, C.-H.H. and C.-Y.K.; formal analysis, C.-H.H. and C.-Y.K.; investigation, C.-H.H. and C.-Y.K.; writing—original draft preparation, C.-H.H. and C.-Y.K.; writing—review and editing, C.-H.H.; visualization, C.-Y.K.; supervision, C.-H.H.; project administration, C.-H.H.; funding acquisition, C.-H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part through the Ministry of Science and Technology, Taiwan (TW), Grant number, MOST-109-2221-E-006-048-MY3.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Atotal heat transfer area (mm2)
Cp(Ta)specific heat capacity of air (kJ/kg-K)
Dhhydraulic diameter
ffriction factor
hheat transfer coefficient (W/m2-K)
Hheight of rib (mm)
Jthe cost function
ka(Ta)thermal conductivity of air (W/m-K)
Ltube length with fins (mm)
Ltfin length (mm)
Nnumber of waves
NuNusselt number
ReReynolds number
trib thickness (mm)
Tssolid material temperature (K)
Taair temperature (K)
Twwall temperature (K)
Tininlet air temperature (K)
Toutoutlet air temperature (K)
Volvolume of fins and ribs (mm3)
Wrib position (mm)
Greek symbols
Ώcomputational domain
Θdesign variable
ΨJacobian matrix
ϕdesired thermal performance factor
δfin thickness (mm)
ρ(Ta)air density (kg/m3)
μweighting coefficient
μl(Ta)air dynamic viscosity (kg/m-s)
ηthermal performance factor
εstopping criterion
ΔPpressure drop (Pa)
ΔTLMTDlogarithmic mean temperature difference (K)

References

  1. Webb, R.L. Principles of Enhanced Heat Transfer; John Wiley & Sons, Inc.: New York, NY, USA, 1994. [Google Scholar]
  2. Fabbri, G. Heat transfer optimization in internally finned tubes under laminar flow conditions. Int. J. Heat Mass Transf. 1998, 41, 1243–1253. [Google Scholar] [CrossRef]
  3. Fabbri, G. Optimum profiles for asymmetrical longitudinal fins in cylindrical ducts. Int. J. Heat Mass Transf. 1999, 42, 511–523. [Google Scholar] [CrossRef]
  4. EI-Sayed, S.A.; EL-Sayed, S.A.; Abdel-Hamid, M.E.; Sadoun, M.N. Experimental study of turbulent flow inside a circular tube with longitudinal interrupted fins in the stream wise direction. Exp. Therm. Fluid Sci. 1997, 15, 1–15. [Google Scholar] [CrossRef]
  5. Al-Sarkhi, A.; Abu-Nada, E. Characteristics of forced convection heat transfer in vertical internally finned tube. Int. Commun. Heat Mass Transf. 2005, 32, 557–564. [Google Scholar] [CrossRef]
  6. Wang, Q.; Lin, M.; Zeng, M.; Tian, L. Investigation of Turbulent Flow and Heat Transfer in Periodic Wavy Channel of Internally Finned Tube with Blocked Core Tube. J. Heat Transf. 2008, 130, 061801. [Google Scholar] [CrossRef]
  7. Yu, B.; Tao, W.-Q. Pressure drop and heat transfer characteristics of turbulent flow in annular tubes with internal wave-like longitudinal fins. Heat Mass Transf. 2004, 40, 643–651. [Google Scholar] [CrossRef]
  8. Wang, Q.; Lin, M.; Zeng, M. Effect of lateral fin profiles on turbulent flow and heat transfer performance of internally finned tubes. Appl. Therm. Eng. 2009, 29, 3006–3013. [Google Scholar] [CrossRef]
  9. Duan, L.; Ling, X.; Peng, H. Flow and heat transfer characteristics of a double-tube structure internal finned tube with blossom shape internal fins. Appl. Therm. Eng. 2018, 128, 1102–1115. [Google Scholar] [CrossRef]
  10. Chai, L.; Xia, G.; Zhou, M.; Li, J.; Qi, J. Optimum thermal design of interrupted microchannel heat sink with rectangular ribs in the transverse microchambers. Appl. Therm. Eng. 2013, 51, 880–889. [Google Scholar] [CrossRef]
  11. Xia, G.; Zhai, Y.; Cui, Z. Numerical investigation of thermal enhancement in a micro heat sink with fan-shaped reentrant cavities and internal ribs. Appl. Therm. Eng. 2013, 58, 52–60. [Google Scholar] [CrossRef]
  12. Ahmed, H. Optimization of thermal design of ribbed flat-plate fin heat sink. Appl. Therm. Eng. 2016, 102, 1422–1432. [Google Scholar] [CrossRef]
  13. Marquardt, D.W. An Algorithm for Least-Squares Estimation of Nonlinear Parameters. J. Soc. Ind. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
  14. Chen, P.F.; Huang, C.H. An Inverse Hull Design Problem in Optimizing the Desired Wake of Ships. J. Ship Res. 2002, 46, 138–147. [Google Scholar]
  15. Huang, C.-H.; Li, Y.-Y. A three-dimensional shape design problem in determining the boundary geometry to yield isotherms. Numer. Heat Transf. Part A Appl. 2019, 76, 517–532. [Google Scholar] [CrossRef]
  16. Huang, C.-H.; Chen, M.-H. An estimation of the optimum shape and perforation diameters for pin fin arrays. Int. J. Heat Mass Transf. 2019, 131, 72–84. [Google Scholar] [CrossRef]
  17. Huang, C.-H.; Tung, P.-W. Numerical and experimental studies on an optimum Fin design problem to determine the deformed wavy-shaped heat sinks. Int. J. Therm. Sci. 2020, 151, 106282. [Google Scholar] [CrossRef]
  18. ESI-CFD Inc. CFD-RC User’s Manual. Available online: https://myesi.esi-group.com/downloads/software-documentation/ace-suite-2020.5-user-guide-cfd-viscart-online (accessed on 1 November 2020).
  19. Yang, Y.-T.; Peng, H.-S. Numerical study of the heat sink with un-uniform fin width designs. Int. J. Heat Mass Transf. 2009, 52, 3473–3480. [Google Scholar] [CrossRef]
  20. Rahimi, M.; Shabanian, S.R.; Alsairafi, A.A. Experimental and CFD studies on heat transfer and friction factor characteristics of a tube equipped with modified twisted tape inserts. Chem. Eng. Process. Process Intensif. 2009, 48, 762–770. [Google Scholar] [CrossRef]
  21. Eiamsaard, P.; Piriyarungroj, N.; Thianpong, C. A case study on thermal performance assessment of a heat exchanger tube equipped with regularly-spaced twisted tapes as swirl generators. Case Stud. Therm. Eng. 2014, 3, 86–102. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The physical model of PTHE-Z. (a) three-dimensional view, (b) cross-sectional view and (c) computational region.
Figure 1. The physical model of PTHE-Z. (a) three-dimensional view, (b) cross-sectional view and (c) computational region.
Energies 13 06424 g001
Figure 2. Grid system: (a) cross-sectional view and (b) three-dimensional view.
Figure 2. Grid system: (a) cross-sectional view and (b) three-dimensional view.
Energies 13 06424 g002
Figure 3. The physical model of design A PTHE-ZR. (a) three-dimensional view, (b) cross-sectional view and (c) computational region.
Figure 3. The physical model of design A PTHE-ZR. (a) three-dimensional view, (b) cross-sectional view and (c) computational region.
Energies 13 06424 g003
Figure 4. The physical model of design B PTHE-ZR. (a) three-dimensional view, (b) cross-sectional view and (c) computational region.
Figure 4. The physical model of design B PTHE-ZR. (a) three-dimensional view, (b) cross-sectional view and (c) computational region.
Energies 13 06424 g004
Figure 5. The physical model of design C PTHE-ZR. (a) three-dimensional view, (b) cross-sectional view and (c) computational region.
Figure 5. The physical model of design C PTHE-ZR. (a) three-dimensional view, (b) cross-sectional view and (c) computational region.
Energies 13 06424 g005
Figure 6. The physical model of design D PTHE-ZR. (a) three-dimensional view, (b) cross-sectional view and (c) computational region.
Figure 6. The physical model of design D PTHE-ZR. (a) three-dimensional view, (b) cross-sectional view and (c) computational region.
Energies 13 06424 g006
Figure 7. The comparisons of (a) Tout and (b) f for PTHE-Z with the corrected values [8] using five different physical properties of air.
Figure 7. The comparisons of (a) Tout and (b) f for PTHE-Z with the corrected values [8] using five different physical properties of air.
Energies 13 06424 g007
Figure 8. The comparisons of temperature distributions for (a) PTHE-Z, (b) design A (c) design B (d) design C and (e) design D at z = 100 mm.
Figure 8. The comparisons of temperature distributions for (a) PTHE-Z, (b) design A (c) design B (d) design C and (e) design D at z = 100 mm.
Energies 13 06424 g008
Figure 9. The comparisons of pressure distributions for (a) PTHE-Z, (b) design A (c) design B (d) design C and (e) design D at z = 100 mm.
Figure 9. The comparisons of pressure distributions for (a) PTHE-Z, (b) design A (c) design B (d) design C and (e) design D at z = 100 mm.
Energies 13 06424 g009
Figure 10. The comparisons of cross-sectional view of (a) temperature and (b) pressure distributions for various designs at z = 100 mm.
Figure 10. The comparisons of cross-sectional view of (a) temperature and (b) pressure distributions for various designs at z = 100 mm.
Energies 13 06424 g010
Figure 11. The variations of (a) Nu, (b) f and (c) η with Re for various designs.
Figure 11. The variations of (a) Nu, (b) f and (c) η with Re for various designs.
Energies 13 06424 g011
Table 1. Design variables and optimal results of PTHE-ZRs with various Re.
Table 1. Design variables and optimal results of PTHE-ZRs with various Re.
Design
Cases
Design
Variables
(mm)
Re = 3000Re = 5000
(Design Condition)
Re = 7000
Nuf η Nuf η Nuf η
PTHE-ZN/A94.20.03201112.00.02851124.40.02401
Case A of
PTHE-ZR
W = 5.16, h = 0.516
δ = 0.18641
118.10.06051.0990148.40.04601.1300166.50.03821.1462
Case B of
PTHE-ZR
W1 = 5.01, h1 = 0.50
W1 = 5.44, h2 = 0.55
δ = 0.18623
118.80.06141.0998150.00.04721.1321169.10.03951.1516
Case C of
PTHE-ZR
W1 = 3.90, h1 = 0.31
W2 = 5.54, h2 = 0.57
δ = 0.17682
122.20.06511.1097153.00.04981.1343172.10.04161.1518
Case D of
PTHE-ZR
W1 = 3.85, h1 = 0.31
W2 = 5.13, h2 = 0.56
W3 = 6.14, h3 = 0.55
δ = 0.16274
123.90.07061.0952154.30.05391.1147173.70.04531.1294
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, C.-H.; Kuo, C.-Y. A Nonlinear Inverse Design Problem for a Pipe Type Heat Exchanger Equipped with Internal Z-Shape Lateral Fins and Ribs. Energies 2020, 13, 6424. https://doi.org/10.3390/en13236424

AMA Style

Huang C-H, Kuo C-Y. A Nonlinear Inverse Design Problem for a Pipe Type Heat Exchanger Equipped with Internal Z-Shape Lateral Fins and Ribs. Energies. 2020; 13(23):6424. https://doi.org/10.3390/en13236424

Chicago/Turabian Style

Huang, Cheng-Hung, and Chih-Yang Kuo. 2020. "A Nonlinear Inverse Design Problem for a Pipe Type Heat Exchanger Equipped with Internal Z-Shape Lateral Fins and Ribs" Energies 13, no. 23: 6424. https://doi.org/10.3390/en13236424

APA Style

Huang, C. -H., & Kuo, C. -Y. (2020). A Nonlinear Inverse Design Problem for a Pipe Type Heat Exchanger Equipped with Internal Z-Shape Lateral Fins and Ribs. Energies, 13(23), 6424. https://doi.org/10.3390/en13236424

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