Abstract

Fully developed laminar flow and heat transfer performance of V-pattern folded core sandwich structures subjected to forced convection are numerically investigated. Based on the periodic feature of folded core, a methodology to simulate the fully developed flow and heat transfer state with a few unit cells is developed. The influences of geometric parameters on the coolant flow resistance and heat transfer behaviour are evaluated. Afterwards, approximate models are established based on design to optimize the heat dissipation efficiency. Design variables considered are geometric parameters of the V-pattern folded core, while the corresponding relative density is taken as a constraint. The simulation results show that the presented methodology makes it possible to analyse the fully developed regime with less computing resources. It is found that both pressure drop and heat transfer rate increase with the tortuosity of the flow channel and the inclination of the folded core. The optimization generates a design with a mildly tortuous feature in the streamwise direction to gain a maximum of heat dissipation efficiency under the current conditions.

1. Introduction

Sandwich structures are widely used in aerospace industry due to their lightweight, high specific strength, and high specific stiffness characteristics [1]. A folded core sandwich structure features various geometric parameters together with its open cell geometry, which gives it the potential for simultaneous load bearing and heat dissipation [2, 3]. It consists of two face sheets and a three-dimensional core which can be folded from a planar base material along certain geometric lines [4, 5]. The folded core can be manufactured from a variety of sheet materials and can be tailored into different geometries meeting different functional requirements. In addition, the open cell design of folded core allows internal fluid transport, further improving the heat dissipation capacity [3]. The mechanical and impact behaviours of V-pattern folded core sandwich structure have been paid much more attention and studied both experimentally and numerically in recent years, while little work exists on its active cooling capabilities [3, 57].

The fluid flow and heat transfer characteristics of sandwich structures with honeycomb, prismatic, and lattice cores have been extensively addressed in recent years. Wen et al. [8, 9] investigated the effects of topology on convective heat dissipation within metallic cellular structures both experimentally and numerically; also, an optimization study based on analytical results was performed. The results showed that cellular sandwich heat sinks having triangular cells were superior to square cells, both thermally and mechanically. Kim et al. [10] carried out experimental and numerical studies for lattice-frame materials to identify the local fluid flow features and their influence on the local and overall pressure and heat transfer characteristics. It was found that the horseshoe vortices were formed in front of the vertices of the LFM and the arch-shaped vortices were formed behind the vertices. The delay of the flow separation might lead to an increase in the overall heat transfer. Gao and Sun [11, 12] studied the thermal control of lattice truss core sandwich structures subjected to an internal thermal load under forced water convection. The results suggested that with an increase of inclination angle of truss, the heat transfer efficiency was improved, leading to a reduction of structural maximum temperature. Rakow and Waas [13] conducted a series of numerical and experimental studies to investigate the heat transfer and structural response of actively cooled metal foam sandwich panels. Feng et al. [14] numerically investigated the heat transfer performance of finned metal foam heat sinks subjected to uniform impinging flow and quantified the effects of all relevant configurational parameters of the heat sink upon the thermal performance. It was demonstrated that the finned metal foam heat sinks outperformed the plate-fin heat sinks on the basis of given weight or given pumping power. Shenoy et al. [15] also introduced forced convective flow and heat transfer characteristics in a wavy horizontal channel. Rashidi et al. [16, 17] studied the condensation flow inside rough and smooth nanochannels by attaching the roughness on the lower walls and investigated the effect of rough surface on the interface properties of the condensation flow for different roughness shapes. Results indicated that the condensation process of flow can be empowered by enhancement of the roughness height. Khoshvaght-Aliabadi et al. [1821] performed experimental studies for forced convective flows of different nanofluids through a corrugated wavy channel and discussed the effects of different geometrical parameters and composition of fluid mixtures. Subsequently, laminar convection through the straight minichannel and wavy minichannel with various cross-section geometries was studied numerically, and the results depicted higher values of the heat transfer rate and pumping power for the wavy minichannel compared to the straight minichannel. However, as a promising candidate for multifunctional design, studies on the heat transfer performance of the folded core are still limited, which makes it necessary and indispensable to investigate their heat dissipation capacity.

The present study establishes an approach to simulate the fully developed heat transfer problem of V-pattern folded core sandwich structures. Both friction factor and Nusselt number are evaluated using the method of computational fluid dynamics (CFD). The effects of the folded core geometric parameters on the pressure and heat transfer characteristics are investigated to gain an initial insight into the problem. Finally, the heat dissipation efficiency is optimized against relative density followed by a short discussion on the optimal results.

2. Numerical Model

2.1. V-Pattern Folded Core Sandwich Structures

A typical sandwich structure with V-pattern folded cores is shown in Figure 1. The top-face sheet is subjected to a uniform heat flux while the bottom-face sheet is assumed to be adiabatic. Forced convection is performed in channels (along the direction) formed by the V-pattern folded cores and face sheets, which means that most heat is dissipated and taken away. Specifically, two different types of channels are identified in the structure, i.e., the channel with the triangular cross-section and the channel with the inverted triangular cross-section.

Typically, a unit cell of V-pattern folded core can be completely described by a set of four independent geometric parameters as shown in Figure 2(a). Using terms in [22], the geometric model is defined by the folded core height , half of the zigzag lines pitch , the saw-tooth lines pitch 2 L, and the zigzag lines amplitude . In reality, when core’s web thickness is involved, the folding edge consists more of a plateau than a sharp edge, which reduces the curvature and the extent of potential imperfections. It is assumed that the folded core has uniform web thickness, as shown in Figure 2(b).

2.2. Boundary Conditions and Material Properties

A mass flow rate of 0.0016 kg/s with a constant static temperature of 300 K entering through the streamwise periodic boundary is prescribed for all cases. The direction of the inlet flow is perpendicular to the inlet cross-section. The top-face sheet is subjected to a constant uniform heat flux of 20000 W/m2 while adiabatic wall is set for the bottom-face sheet.

Aluminium 2024-T851 is used for both folded core and face sheets. Kerosene [23] is chosen as the coolant. Since influences of geometric parameters are mainly investigated in this paper, constant material properties are adopted. Material properties are shown in Table 1.

2.3. Data Reduction Parameters

For fully developed flow and heat transfer in a periodic channel, information obtained from a representative unit cell can be utilized to describe the whole structure. To evaluate the flow and heat transfer characteristics, three nondimensional quantities (the Reynolds number , the Darcy friction factor , and the Nusselt number ) are introduced. and based on the folded core height are defined as where is the mean inlet velocity, and are the density and viscosity of the coolant, respectively, and refers to the pressure drop per unit length (i.e., 2S).

At a streamwise location , the perimeter averaged heat transfer coefficient is defined as where and are peripherally averaged wall heat flux and temperature, respectively, and is the bulk fluid temperature. The overall mean heat transfer coefficient in one unit cell can then be obtained from the spatial average of as

Therefore, the corresponding based on folded core height can be defined as

2.4. Computational Grid

Hexahedral grid was used for the whole computational domain, and meshes are refined around the fluid-solid interface to better resolve the boundary layer. The fluid domain and the solid domain were divided by common nodes, which could improve the accuracy of data transmission at the fluid-solid interface. Convergence studies were conducted to establish an estimated accuracy of the results and to establish mesh independence of the numerical solutions. The influence was examined on the basis of the and dependence on the element number. The values of and were taken from the central unit for a typical three-unit cell length V-pattern folded core panel with a uniform velocity and temperature profile at the inlet boundary. The results of the convergence study are presented in Figure 3, and it can be seen that a monotonic convergence toward the final solution is displayed in the figure. In addition, convergence studies were conducted for several other Reynolds numbers and yielded comparable results. The mesh size with 694650 volume elements was eventually settled upon for the model as the resulting solution differed only marginally from a finer mesh. Illustration of a typical computational mesh is shown in Figure 4.

3. Computational Method

3.1. Governing Equations

In the present study, some assumptions are considered in the numerical simulation: (1) the flow and heat transfer are steady; (2) the fluid is incompressible and Newtonian; (3) natural convection and radiation heat transfer are neglected; (4) viscous dissipation heat is not considered; (5) thermophysical properties of the material do not change with temperature. Based on above assumptions, the governing equations can be expressed as follows.

Continuity equation where is the fluid velocity vector.

Momentum equation where is the fluid density, stands for the pressure, and refers to the viscous stress under incompressible condition and can be expressed as

Energy equation where is the fluid temperature and and stand for the specific heat and thermal conductivity of the fluid, respectively.

For the solid part, the steady heat conduction differential equation can be expressed as where is the solid temperature and is the thermal conductivity of the solid.

The temperature and heat flux are continuous in the fluid-solid interface and can be expressed as

3.2. Numerical Procedure

A three-dimensional model is considered to model fully developed, steady-state laminar fluid flow and heat transfer within V-pattern folded core sandwich structures. Due to the periodic characteristic of the V-pattern folded core, flow across the computational unit cell of the folded core sandwich panel is assumed to be fully developed both hydrodynamically and thermally. The assumption is reasonable for a unit cell placed apart from the inlet and outlet region. Thus, both flow and heat transfer calculations in geometries are periodically repetitive in the streamwise direction and spanwise direction . Taking a single representative unit cell from the structure and imposing periodic boundary conditions in two directions is desired to minimize the computational expense. A pressure drop occurs across translationally periodic boundaries in the streamwise direction while no pressure drop occurs across the periodic planes in the spanwise direction. This can be summed up to a periodic heat transfer problem for a wavy structure with two different types of periodic boundary conditions.

Therefore, a two-step simulation method is developed to achieve a fully developed periodic solution, and details are shown in Figure 5. ANSYS CFX is adopted to perform numerical simulations, where a vertex-based finite volume method is used to solve the Navier-Stokes and energy equations. All simulations are conducted using the high resolution advection scheme. RMS residual below 10-6 and global imbalance below 0.01 are taken as the convergence criteria.

First, fully developed periodic flow is solved using ANSYS CFX without heat transfer. Two different types of periodic boundary conditions are defined in two directions with a mass flow rate specified in the streamwise direction. Both velocity and pressure field are computed in a single simulation.

Sequentially, the heat transfer is solved leaving the flow field unaltered. To determine the heat transfer characteristics of the structure under thermal fully developed conditions, an iterative procedure is required. The energy transport equation is first solved with a uniform inlet temperature and the velocity field obtained from the previous step. Afterwards, the outlet temperature field of the first calculation is scaled by the difference between the mass flow averaged outlet temperature and the initial inlet temperature and then reapplied to the inlet for the second iteration. This process is repeated several times until a thermally developed state is achieved. Here, a stopping criteria based on the average Nusselt number is given by

However, the temperature gradient information was not included during iterations [24]. Heat fluxes across the solid periodic interface were also neglected as adiabatic boundary conditions were used. To overcome this, a number of repeating units of geometry need to be included in order to ensure that the heat transfer is fully developed in the unit of interest. Generally, three repeating units are sufficient when all data are extracted from the central unit, as shown in Figure 6. All above procedures are integrated using Isight [25], which makes a good preparation for later analyses and optimization.

For cases of constant wall heat flux, which is mainly considered in this paper, temperature profiles are similar in shape at streamwise locations separated by the period length . And the temperature difference between two locations is proportional to the wall heat flux [26]. To gain an initial understanding of the issue, a typical three-unit folded core structure (see Figure 6) is analysed using the previously mentioned two-step method. Temperature profiles of the inlet and outlet for the first heat transfer analysis and iteration five are shown in Figures 7(a) and 7(b), respectively. It is seen that the shapes of the temperature profiles at the entrance and exit are similar after five iterations. It can be considered that the thermally fully developed state has been reached at the 14th cell.

3.3. Method Validation

In order to confirm the accuracy of the method, the Fanning friction factor and Nusselt number for laminar fully developed flow and heat transfer in a straight channel with an equilateral triangular cross-section were calculated. Table 2 shows a comparison of product and with values taken from reference [27]. The agreement is greater than 99% for both parameters.

4. Results and Discussion

Folded core sandwich structure has an excellent designability due to its various geometric parameters. The influences of geometric parameters on the flow and heat transfer characteristics are investigated. The folded core sandwich structure affects the coolant flow mainly by the flow channel’s tortuosity and cross-section, which can be characterized by the geometric parameters , , , and , respectively. Two independent parameters and , as shown in Figure 2(a), with a constant folded core unit cell length () and height () are investigated.

Figure 8 shows the flow and heat transfer characteristic variations with the folded core tortuosity along the flow direction. In particular, the wavy channel turns into a straight one when equals to zero. It is seen that both pressure drop and heat transfer performance increase with the angle , which is quite obvious since flows in tortuous channels normally give rise to pressure drop penalty and heat transfer enhancement relative to flow in straight channels. It is noted that the heat transfer enhancement is greater than the pressure drop penalty at small turning angles. However, with the continuing increase of the tortuosity, the pressure drop grows rapidly while the heat transfer enhancement slows down.

Figure 9 shows the influence of the inclined angle on the flow and heat transfer characteristics. It is seen that both indexes increase with the cross-section changing from an obtuse triangle to an acute triangle. The increase rate of the pressure drop is greater than that of the heat transfer performance before the inclined angle reaches about 50°. From Figure 8, structures with a small are more likely to offer efficient heat dissipation. However, it should be noted that constant mass flow rate and folded core height are used in this investigation. Thus, the coolant inlet velocity, as well as the Reynolds number, varies with the cross-sectional area. There might be a different relation among these parameters investigated at a constant Reynolds number, which is not included in the current study.

5. Optimum Design

Normally, a structure with good heat transfer behaviour is inclined to suffer a relatively large pressure loss, which makes it necessary to conduct the optimum design for efficient heat dissipation. The heat dissipation performance of the structure can be easily gauged by the ratio of the average Nusselt number to the friction factor [28]. Consequently, the corresponding objective function of the optimization problem is simply given by

All four geometric parameters are considered as design variables while the folded core and face sheet thicknesses are taken as a fixed value of 1 mm. Based on the previous study of the influence of geometric parameters on the flow and heat transfer performance, a specific restriction which would avoid severely tortuous flow channel or extremely obtuse-triangular and acute-triangular cross-section was imposed when setting the ranges of design variables. The corresponding parameter ranges are shown in Table 3. The relative density of the folded core is taken as a constraint to prevent the structure from being too heavy; thus, an upper bound of core’s relative density is specified.

The optimization problem can be succinctly described as:

First, 200 samples were selected in the design space by using Latin Hypercube Sampling (LHS) technique. These samples were simulated using the previously proposed CFD model and then used to construct approximate models based on the radial basis functions (RBF). 20 random points were used to evaluate the accuracy of the approximate models. Details are shown in Table 4, where RMSE represents root mean square error and refers to the coefficient of determination. It is seen that the approximate models are of high precision and can be used to predict the structural performances.

Afterwards, RBF approximate models were invoked to replace the time-consuming CFD models during optimization where sequential quadratic programming was employed for seeking optimal solutions. The optimal geometric dimensions are shown in Table 5. It is seen that is relatively small to gain a higher , as the pressure loss increases dramatically for a severely tortuous channel. An additional CFD analysis for the optimal parameters was also conducted. The comparison of the RBF predictions and CFD results is shown in Table 6. It is evident that the agreement is excellent, with less than 1% difference for all interested parameters, which indicates that the approximate models have fairly good accuracy.

Specifically, the CFD analysis using the previously mentioned two-step simulation approach shows that eight iterations are needed for the optimal parameters to satisfy the stopping criteria. Figure 10 shows how the perimeter averaged Nusselt number varies through the central unit along the flow direction with increasing iterations. It is seen that the curves were almost coincide with each other after five iterations, which further supports the proposed simulation methodology.

6. Conclusions

In this paper, the heat transfer performance and flow state of coolant in a folded core sandwich structure subjected to forced convection conditions have been numerically investigated. A two-step approach to model fully developed laminar flow and heat transfer process is presented. The key conclusions can be drawn as follows: (1)The proposed two-step simulation method can simulate periodic heat transfer with only a few unit cells, which greatly reduces the amount of calculation in every single run. Simulation of a straight channel with an equilateral triangular cross-section was used to validate the accuracy and feasibility of the methodology(2)The influences of geometric parameters on the flow and heat transfer characteristics are investigated. Two independent parameters and were introduced, which characterize the tortuosity of the channel and the angle of the cross-sectional triangle, respectively. The results show that with the continuing increase of the tortuosity, the pressure drop grows rapidly while the heat transfer enhancement slows down(3)Finally, the geometric parameters of folded core, , , , and , were optimized based on an RBF approximation model with the objective of higher heat transfer efficiency and lower relative density. The optimization results show that a structure with a gently tortuous flow channel dissipates heat more efficient for conditions used in the current study

Nomenclature

:Fanning friction factor
:Specific heat capacity, J·kg-1·K-1
:Darcy friction factor
:Heat transfer coefficient, W·m-2·K-1
:Folded core height, m
:Thermal conductivity, W·m-1·K-1
:Pressure, Pa
:Heat flux, W·m-2
:Reynolds number
:Temperature, K
:Velocity, m·s-1
:Velocity vector, m·s-1
:Dynamic viscosity, N·s·m-2
:Density, kg·m-3
:Viscous stress, Pa
:Angle, degree.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the Funding of Jiangsu Innovation Program for Graduate Education (grant no. CXLX13_163), the Fundamental Research Funds for the Central Universities (grant no. NZ2016101), and a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).