Evaluation of the effects of the metal foams geometrical features on thermal and fluid-dynamical behavior in forced convection

Metal foams are a material, featuring interesting characteristics for the aeronautical and automotive fields because of their low specific weight, high thermal properties, and mechanical performances. In particular, this paper deals with thermal and fluid dynamic study of 24 open-cell aluminum EN43500 (AlSi10MnMg) metal foams produced by indirect additive manufacturing (I-AM), combining 3D printing and metal casting to obtain a controllable morphology. A study of foam behavior function of the morphological features (pores per inch (PPI), branch thickness (r), and edges morphology (smooth-regular)) was performed. The samples produced were heated by radiation and tested in an open wind circuit gallery to measure the fluid dynamic properties such as pressure drop (Δp), inertial coefficient (f), and permeability (k), in an air forced convection flow. The thermal characterization was performed evaluating both the theoretical (kth) and effective (keff) thermal conductivity of the foams. Also, the global heat transfer coefficient (HTCglobal) was evaluated with different airflow rates. Analysis of variance (ANoVA) was performed to figure out which geometrical parameters are significant during both thermal and fluid dynamic processes. The results obtained show how the controllable foam morphology can affect the involved parameters, leading to an ad hoc design for industrial applications that require high thermo-fluid-dynamical performances.


Introduction
Metal foams are a class of cellular materials with cells randomly oriented and homogeneous in size and shape. Due to their high surface area per unit volume, many authors discussed the properties for a large number of engineering applications, such as condenser for geothermal power plants, groundwater pollution and filters, chemical reactors, vibration control, and energy management, combustion chambers, strain isolation, catalytic beds, air conditioning, and refrigeration systems [1][2][3][4][5][6][7].
Because of the high surface-volume ratio, the open-cell cellular materials offer excellent performance in heat exchange, due to their unusual structure that enhances the flow mixing ability and the production of turbulence in the fluid flow when passing through the foam structure. In particular, Buonomo et al. [8,9] evaluated numerically the improvement of the energy performance ratio (EPR) for different Reynolds numbers for an aluminum foam included in a heat exchanger. Also, Gankapatnam et al. [10] and De Carvalho et al. [11] investigated numerically respectively the influence of the pore scale on pressure drops and thermal behavior. Instead, Baiocco et al. [12] developed a prototype of a heat exchanger with open-cell metal foams used as an exchange surface, focusing on the evaluation of the heat transfer coefficient. Also, metal foams represent a suitable substrate for the electroplating processes despite their complex structure. This allows a localized surface modification and an enhancement of their features [13][14][15][16].
Considering all the aforementioned structural and thermal properties, it is important to classify and evaluate the physical properties of the foams.
Different production processes can be used to obtain structures with similar morphologies, densities, and cell size. In melt gas injection (MGI), bubbles are created through air embedding into molten metal. These bubbles rise to the surface and create a foam [17]. Powder metallurgy (PM) consists of mixing metal powders, a blowing agent, and ceramic particles aimed at modifying the molten metal viscosity and stabilizing the foam. The mix is compacted to form a dense precursor, which is heated above the melting temperature of the metal. The blowing agent is chosen to decompose at a lower temperature, to release gas before the powder melting. The porous structure is preserved by fast cooling to room temperature [18,19].
A different way to produce metal foams is investment casting (IC). The IC process consists of liquid metal pouring into a ceramic cavity designed to reproduce the desired part through wax patterns [20][21][22]. Jiang et al. have proposed another casting process for the ceramic cavity, which is the expendable pattern shell casting process with vacuum and low pressure (EPSC-VL) to produce thin-walled aluminum, magnesium alloy, and A356 aluminum alloy castings [23][24][25][26]. Thanks to the gas pressure and vacuum pressure during casting, both the filling ability and feeding capacity of the molten metal are improved with this kind of process. Recently, a new method for the production of metal foams called friction stir processing (FSP) has been developed [27][28][29].
In this process, a mixture of blowing agent powder (TiH 2 ) and stabilization agent powder (Al 2 O 3 ) is used to form a composite through FSP. This composite, when heated up, converts into an aluminum metal foam, due to the evolution of hydrogen gas from the blowing agent. In this production process, it is not possible to exert a thorough control of the foam morphology. However, the regulation of the parameters can lead to the overall control of the samples in terms of pores per inch (PPI) [30].
To overtake this issue, an additive manufacturing technology (AM) can be used. AM allows the realization of components with complex geometries, starting from a 3-D computeraided design (CAD) system [31,32]. AM represents an effective and viable alternative in the production of items characterized by complex geometries when the standard subtractive processes are too expensive. The samples realized with this kind of process are valuable for structural applications [33][34][35][36][37]. In particular, Calleja et al. [38][39][40] worked on the improvement of laser cladding technology for the realization of structural customized parts with structural performances and reduced weight. The papers are focused on the optimization of the steps to be followed before accomplishing five-axis laser cladding operations by taking into account the machine kinematics.
Anyway, limitations such as costly machinery, difficult parameter settings, and the reduced size of the components produced restrain the application of this technology. Also, for this kind of reticular structures, the remotion of the internal supports for structural sustain and heat dissipation could be too hard if not impossible to achieve. To avoid these critical issues, an indirect additive manufacturing process can be exploited. Almonti et al. produced a foam structure made of a foundry resin using a 3D stereolithography laser printer. Following, this model was reproduced with a casting process to achieve a metal foam with the designed structure [41,42].
Since the cellular structure model is produced with CAD and printed by overlapping micrometric layers, it was possible to control different foam parameters such as PPI, branch thickness (r), and branch connection nodes typology. The connection between nodes was made with sharp edges (regular) or smoothed (smooth). This was not possible with the previous aforementioned production processes, where the cavity displacement was random. Therefore, the work aims to relate these geometrical parameters (PPI, r, smooth-regular) to the thermal and fluid-dynamical behavior in forced convection heat exchange process. The characterization of the foams was performed by measuring the pressure drops (Δp), global heat transfer coefficient (HTC global ), permeability coefficient (k), and inertia factor ( f ).
For the evaluation of thermal conductivity (k th ), the heat conduction phenomenon in the branches was considered significant. Indeed, between the heat convection in the gaseous and solid phase, there is a difference of about 5 orders of magnitude. In this work, the theoretical value of thermal conductivity k th was calculated in accordance with the literature [43].
Finally, an analysis of variance (ANoVA) was carried out. There is an extremely vast literature where the ANoVA was used to understand the significant factors on experimental results [44][45][46]. In the present work, it was applied to evaluate if the geometrical parameters and their linear combination affect the airflow.

Materials and methods
Different foams made of EN43500 aluminum alloy with dimensions 20 × 20 × 60 mm were realized and characterized. The samples produced were characterized by different PPI values and branches thickness (r) as showed in Table 1; also, two different kinds of edges (smooth/regular) were considered. The structure of the samples was designed with the ad hoc method described in references [41,42]. A Spatial Voronoi Tessellation with the Rhinoceros plug-in GrassHopper was used to design the structure. A cell-by-cell construction iterative algorithm was used to build the external Tessellation. After defining the volume, it was populated with a cloud of points which define the cell cores, based on the number of PPI desired. At the end of this process, different thickness of the structural branches and smoothing was assigned. The CAD file obtained after a cleanup of the repeated entities and open surfaces is shown in Fig. 1. Once the CAD is obtained, it is guaranteed the reproducibility of the foams since the same geometry can be used for replicas or future developments.
The sample models were realized with a photosensitive resin FUSIA 444 using a stereolithography 3D printer XFAB 2000 DWS. The printing area is a cylinder with a diameter of 180 mm. The printing process is realized with a Solid State BluEdge Be-1300X laser, in a slicing range between 10 and 100 μm The choice of the resin was made by considering the suitability for casting thin and solid models. The reticular structure of the foams allows not to use supports in the inner part of the models during the printing process [47]. Following the 3D printing stage, the samples were polymerized for 1 h in a UV oven to complete the foundry model of the foam, as showed in Fig. 2.
When the reticular structure was achieved, it was exploited to realize the plaster mold for the aluminum foam casting. The plaster suspension was produced by mixing plaster powder and water at 25°C for 4 h, to obtain a homogenous compound. Thereafter, a vibrating platform and a vacuum machine were exploited to avoid air inclusions in the suspension, to prevent defects during the casting. After, the suspension was inserted in the foundry cylinder to form the plaster mold. Finally, the plaster mold was dried at room temperature for 2 h since residual humidity in the plaster mold may lead to a shift of the model during the heating phase. A thermal cycle up to 730°C was imposed on the plaster mold to solidify the plaster die and remove the resin. Also, the plaster must reach the casting temperature to reduce thermal stress during the foundry process. The aluminum EN43500 (AlSi10MnMg) used for the production of the samples was heated up to a casting temperature of 730°C and mixed with low-frequency pulses. The casting process was realized with the foundry machine ASEG Galloni G5 which enables the aluminum melting in the inert atmosphere (Argon). Subsequently, the molten metal fills the mold through a vacuum. Following the casting process, the plaster was removed from the cylinder after a 6-h cooling. Then, the pouring channels were removed and the foam samples obtained were removed from their sprues and cut to the dimension of 20 × 20 × 60 mm, as showed in Fig. 3.
A schematic diagram of the manufacturing process for the realization of the tested foams is shown in (Fig. 4). The experimental setup was an open-circuit wind gallery, presented in Fig. 5. Two steel pipes of 1.5 m with a square section of 20  × 20 mm were linked to a sample holder of ABS, the flanges were connected by screws. The foam was placed in the sample holder in a cavity with the same dimension of the pipes, to guarantee the continuity of the airflow. A hermetic window of transparent PVC was realized on the side of the holder, for monitoring the temperature of the foam with a thermal camera. To measure the inner and the outer temperature of the airflow, two thermocouple type K were placed before and after the sample holder. The temperature values were functional for the evaluation of the heat absorbed by the air while flowing across the foam.
The pressure drop (Δp) was measured close to the sample holder with a differential pressure gauge PCE-P01. Due to avoid measurement errors on Δp and achieve good repeatability, a test without any foam was performed. The porosity was obtained by calculating the relative density (ρ*ρ −1 ): where ρ* is the density obtained experimentally from the ratio between mass and volume and ρ is the material density.  [48] calculated the heat transfer surface by modeling the structure by a twodimensional array of hexagonal cells where the fibers were from the sides of the hexagons. In this paper, since the foam structures were designed with CAD models, the exchange surface, number, and mean diameter of cells were evaluated through the CAD software (Table 2). A more focused analysis was performed by calculating the permeability coefficient (k) and the inertia factor ( f ). Different authors discussed the calculation of k and f since they are strongly related to the structure of the foam. Beavers and Sparrow [49] were among the first authors who quantified pressure drops and f for nickel foams. The subsequent study of Vafai and Tien [50] measured k and f for "Foametal" and found out reference values for the studied foams. Experiments on high porosity metal foams (0.89 < ε < 0.96) were conducted also by Pack et al. [51], who analyzed pore densities corresponding to 10, 20, and 40 PPI.
For the calculation of the aforementioned coefficients, the momentum equation for fully developed flow with small viscous effects is [47]: where L is the foam length, v is the air flow velocity, and a and b are coefficients obtained by a linear regression with R 2 > 99.5%. Through a and b, the k and f values were evaluated by means of the following equations: A compressor produced the airflow while the inlet speed was regulated with two control valves and measured with a PCE-433 hot wire anemometer. An analysis with 3 different velocities (5 m s −1 , 7.5 m s −1 , 10 m s −1 ) was performed.
In our work, the HTC global was evaluated according to Mancin et al. [52]. Particularly, they performed a thermal balance between the supplied electrical power provided P EL and the heat flow rate was verified.
In this paper, a 50-W led lamp was placed under the sample holder to heath through irradiation. The test lasts 120 s, in particular in the first 90 s the sample was heated up; in the last 30 s, the compressor is turned on and generated the airflow. The steady-state condition of the flow was reached after 2 s and guaranteed by the hot wire anemometer. The temperature was monitored with a thermal imager; in particular, 9 core  points were considered to have a higher accuracy of the measurements as shown in Fig. 6. Thermal imaging is strongly dependent on the emissivity term (ε) in the Stephan-Boltzmann equation: where P is the power radiated by the sample, A is the surface of the sample, ε the emissivity of the sample, and T the measured temperature. The (ε) of the samples was changed and approximated as a black body (i.e., ε = 1) by coating the tested foams with graphite powder before temperature measurements with the thermal imager.
From the experimental measurements, it was possible to define the power absorbed by the air (P EL ) as: Fig. 8 Δp=L v at various PPI branch thickness (r) and structure (smooth/regular) vs velocity inlet (v)  whereṁ air is the air mass flow, c p;air is the specific heat of the air at constant pressure, T air, out is the air outlet temperature, and T air, in is the air inlet temperature. The HTC global is defined as: where A is the heat transfer surface, ΔT ml is the logarithm mean temperature defined as: Fig. 9 Results of HTC global at various air mass flow rates for all foams at various PPI, diameter ligaments (r), and geometry (smooth/regular) where T w,in, and T w,out are respectively the inlet and outlet temperatures on the foam walls. T w,in was calculated as the medium value of the cursors 1-2-3, while T w,out was calculated as the medium value of the cursors 7-8-9 as shown in Fig. 6.
The theoretical thermal conductivity (k th ) was obtained as follow: the fluid-phase contribution is considered a product between porosity (ε) and thermal conductivity of the fluid (k f ) and solid-phase contribution is considered the product between the remaining volume part of the porous material (1 − ε) and thermal conductivity of metal (k s ).  The previous formulation does not take into account the morphology of the structure but only the porosity of the sample under consideration.
Instead, the effective thermal conductivity (k eff ) was evaluated by applying the following equation: whereQ is the heat power exchanged for radiation by the UV lamp, A is the heat transfer surface, T s is the starting environment temperature, T e the temperature at the ending of the heating, and s is the distance between fronts defined as ( Fig.  6) follows: top front (cursor 1-4-7), middle front (cursor 2-5-8), and the lower front (cursor 3-6-9).
The measurements produced for the analysis of variance (ANoVA) was structured in a preliminary plan as showed in Table 3.

Results
The foams realized through I-AM are reported in Table 4.
With the 10 PPI foams, the reduced cell size did not allow the removal of the plaster from the core of the structure, displayed in Fig. 7. For this reason, characterization tests on 10 PPI foams were not taken into consideration since the flow would be affected by the plaster entrapped into the samples. Also, the HTC global measurements would be influenced due to Fig. 11 Main effects and interaction plots for Δp the different thermal exchange coefficient of the plaster compared to pure aluminum. Figure 8 shows the normalized pressure drop ( Δp=L v ) plotted concerning the inlet velocity, PPI value, and branch thickness.
The trends show a ( Δp=L v ) increment with PPI value, branch thickness, and velocity inlet. This result can be explained because as the PPI and r increase, the foam internal structure increase in complexity, which means that the turbulence develops easier. Also, r and the PPI increment reduce the flow section and, therefore, an increment of the Δp. The measures were repeated 3 times to guarantee a higher accuracy and are in good agreement with the literature [52].
Through the linear regression shown in Fig. 8, the a and b coefficients were calculated. Successively, k and f were calculated considering μ air and ρ air at environment condition, as shown in Table 5. The results achieved are shown in Table 6 and agree with the reported literature [48,53]. The linear evolution of pressure drop values can be explained since the foam geometrical characteristics analyzed changes, also, linearly. Figure 9 reports the HTC global versus air mass flow rate curves as a function of PPI, nodes geometry, and r. These measurements were performed to state the influence of the foam morphology on the heat transfer. All the trends obtained matches the literature results [52,53].
The trends of HTC global increase with branch thickness. The influence of r is significant for lower PPI foams in which the higher thickness of the branches influences the fluid-dynamical behavior by increasing the turbulence of the flow. For the 8 PPI ones, since the geometry is already complex, the thermal Fig. 12 Main effects and interaction plots for HTC global performance is not influenced by this parameter. Plots for 8 PPI foams show a different slope and behavior at r 0.3 compared to r 0.5, r 0.7. Also, the thermal performance of all foams increases at a higher air mass flow rate. This result can be explained since the increasing of the turbulence in the flow easies the thermal dissipation, which translates into a higher HTC global . The graphs in Fig. 9 show that the influence of the branch thickness on the HTC global is higher between r 0.3 and 0.5, while lower between r 0.5 and 0.7 mm. This indicates that the thermal inertia of the specimens increases in function of their mass. The influence of the geometry edges (regular/smooth) was analyzed; in particular, smoothed edges showed higher HTC global increase than regular edges. This could be justified by the particular conformation of the smoothed branches since they have a similar shape of conical fins and favor a uniform temperature condition, and a more effective thermal exchange.
Regarding the thermal conductivity, theoretical values k th were calculated through the knowledge of ε. The working fluid is air so k f is 0.026 W m −1 K −1 and samples were constituted by aluminum so k s is 290 W m −1 K −1 . The results are shown in Table 7.  To achieve a conductivity value related to the experimental data, the temperatures recorded by the 9 cursors were divided into 3 temperature fronts (top-middle-bottom). The distance between fronts of 7 mm was identified, and A was assumed as the area crossed by thermal power, which is equal to the plate placed as a base of the samples (20 × 60 mm 2 ).
According to previous trials, the power transmitted for the conductive effect is about 80% of emitted power, to take into account the heat losses of the system. Figure 10 shows the difference between trends of k th and k eff .
The results obtained are in agreement with the theoretical values calculated and have percentage errors in a range of 10 to 20%. This can be explained considering the presence of fluid inside the porosities even when the compressor is turned off. The air tends to favor the establishment of convective thermal exchange on the surface of the foams and this inhibits the visibility of the conductive thermal exchange alone.
The results from the ANoVA analysis show which effects are significant in term of Δp, HTC, k, and f.
In particular, in Table 8, the P value for the first two aforementioned parameters is shown. On linear interaction, every effect is significant (P value < 0.05) except the geometry one in terms of pressure drop. It was observed a slight decrease of the drops with smoother geometry ascribable to lower turbulence. This result could be explained because the geometry change is focused only on the edges, causing little geometrical alterations that may not influence the pressure drop. The geometry of the edges does bit play a crucial role; in the 2-and 3-way interaction where it is involved, significant effects were not observed.
The main effects and interaction plots are also shown in Fig. 11 and Fig. 12.
From the HTC global , the main effect plot can be noticed that there is an exponential increase with higher PPI, which means that for thermal applications, higher PPI foams are better. The ANoVA analysis for permeability and inertia coefficient (k, f) did not take into account the inlet velocity, since both coefficients derive from a linear regression on this parameter. For this reason, the analysis had only 2 ·3 2 levels and only linear and 2-way interactions are present, as shown in Table 9. Figures 13 and 14 show the main effect and the interaction plots for k and f. The significant effect value in terms of f is the PPI value. The results show that the f value is highly influenced by 8 PPI foams. For the 8 PPI samples, there is the influence of r as showed from the 2-way interaction plots. r results significantly only in terms of k in the linear interaction. The other interactions are not significant, since the P value is higher than 0.05.

Conclusions
In this paper, the influence of the metal foam morphology on its thermo-fluid dynamic behavior was investigated. In particular, 24 open-cell aluminum EN43500 (AlSi10MnMg) foams were produced through I-AM. With this process, it was possible to realize the reticular structures which characterize the foams without internal supports, unlike the other metallic AM standard processes (EBM, SLM). Moreover, it was possible to reproduce the exact geometries obtained through CAD, to evaluate experimentally the influence of the shape parameters and guarantee the repeatability of the experimental results.
The tested foams differ in PPI, branch thickness, and edge morphology. Tests were carried out in an experimental wind gallery and measurements were performed according to the experimental plan. Due to technological limits, 10 PPI foams could not be tested since it was not possible to remove all plaster from the internal structure.
The results obtained in term of HTC global and Δp show that the trends for different inlet velocity and air mass flow rate are linear and increase with the PPI, branch thickness, and velocity inlet. The max thermal performance is reached by 8 PPI smoothed foam with r 0.7 with 650 W m −2 K −1 . This can be explained as the foam internal structure increases, the turbulence in the flow develops easier, which translates in higher drops and a better thermal exchange. It is notable that with the increase of the PPI value, the branch thickness (r) and the edge morphology (smooth-regular) do not influence the HTC global (2% variation on HTC global at higher velocity). Permeability k and inertia coefficient f of the foams were also calculated for all samples.
In the conductive analysis, the thermal conductivity was calculated with 2 different methods. The first takes only into account the porosity, giving us the theoretical values of k th . The second takes into account the real geometry of the foam and measured parameters during the experimental tests (k eff ). Trends for both coefficients were plotted, and have percentage errors in a range of 10-20%.
An ANoVA analysis was performed to check the significance of the analyzed parameters. Results showed that only the geometry of edges is not a significant factor for Δp and HTC global , although it still has a role in the permeability coefficient k. It was also noted that the branch thickness does not influence f.
Authors' contributions All the authors contributed equally to the various aspects of this work: Daniele Almonti: methodology, formal analysis, software, writing (original draft), investigation Gabriele Baiocco: validation, conceptualization, writing (original draft), formal analysis Emanuele Mingione: formal analysis, investigation, writing (original draft), methodology, software Nadia Ucciardello: conceptualization, methodology, validation, formal analysis Funding Open access funding provided by Università degli Studi di Roma Tor Vergata within the CRUI-CARE Agreement. The presented work was partially supported by the University of Campania Luigi Vanvitelli through the VALERE program (VAnviteLli pEr la RicErca, 2019).
Data availability The processed data generated and analyzed during the current study are not available due to using in another publication, but are available from the corresponding author on reasonable request.

Compliance with ethical standards
Conflict of interest The authors declare that they have no competing interests.

Code availability Not applicable
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.