This paper investigates the thermal properties of a new type of hollow sphere structures. For this new type, the sphere shell is perforated by several holes in order to open the inner sphere volume and surface. The effective thermal conductivity of perforated sphere structures in several kinds of arrangements is numerically evaluated for different hole diameters. The results are compared to classical configurations without perforation. In the scope of this study, three-dimensional finite element analysis is used in order to investigate simple cubic, body-centered cubic, face-centered cubic and hexagonal unit cell models. A linear behavior was found for the heat conductivity of different hole diameters for several kinds of arrangements when the results are plotted over the average density

Hossein Hosseini, Seyed Mohammad

Merkel, Markus

ÄOchsner, Andreas

English

Öchsner, Andreas

Hosseini, Seyed Mohammad Hossein

University of Newcastle's Digital Repository

Numerical simulation of sintered perforated hollow sphere structures (PHSS) to investigate thermal conductivity

Universiti Teknologi Malaysia Institutional Repository

N um erical Sim ulation of Sintered Perforated  Hollow Sphere S tructu res (PHSS) to  Investigate T herm al C onductivityAndreas Ochsner *, Seyed Mohammad Hossein Hosseini ] Markus Merkel *Abstract— This paper investigates the thermal properties of a new type of hollow sphere structures. For this new type, the sphere shell is perforated by several holes in order to open the inner sphere vol­ume and surface. The effective thermal conductiv­ity of perforated sphere structures in several kinds of arrangements is numerically evaluated for differ­ent hole diameters. The results are compared to classical configurations without perforation. In the scope of this study, three-dimensional finite element analysis is used in order to investigate simple cubic, body-centered cubic, face-centered cubic and hexago­nal unit cell models. A  linear behavior was found for the heat conductivity of different hole diameters for several kinds of arrangements when the results are plotted over the average density.Keywords: Cellular m aterials, Thermal conductivity, Finite element method.1 IntroductionHollow sphere structures (HSS) are novel lightweight ma­terials within the group of cellular metals (such as metal foams) which are characterised by high specific stiffness, the ability to absorb high amounts of energy at a rel­atively low stress levels, potential for noise control, vi­bration damping and thermal insulation, cf. Fig. 1. A combination of these different properties opens a wide field of potential multifunctional applications e.g. in au­tom otive or aerospace industry. T ypical functional ap­plications of cellular metals in the scope of heat transfer are heat exchangers [1, 2, 3], fire retardancc systems [4] or thermal insulation. The high thermal insulation capa­bility of hollow sphere structure has been addressed in a US Patent [5] by Schneider and co-workers. Baumeister*Department of Applied Mechanics, Technical University of Malaysia, 81310 Skudai, Malaysia & Centre for Mass and Ther­mal Transport in Engineering Materials, The University of New­castle, Callaghan, New South Wales 2308, Australia, Email: An­dreas.Oechsner@gmail.comt Department of Applied Mechanics, Technical University of Malaysia, 81310 Skudai, Malaysia, Email: rsg.931@gmail.com* Department of Mechanical Engineering, University of Applied Sciences Aalen, 73430 Aalen, Germany, Email: Markus. Merkel® ht w-aalen. dcand colleagues [6] investigated the linear thermal expan­sion coefficient of corundum based hollow sphere compos­ites (HSC) using thermomechanical analysis. They found that the thermal behaviour of HSC is mainly governed by the epoxy resin. Lu and Kou [7] conducted a comprehen­sive numerical and experimental study based on unit cells of homogeneous hollow sphere structures. However, only automatically generated finite element meshes were used and the applied approach does not allow the considera­tion of different material combinations. In recent paper[8], the thermal conductivity of adhesively bonded and sintered HSS was numerically investigated while in paper[9], the influence of material parameters and geometrical properties of syntactic (hollow spheres completely em­bedded in a matrix) and partially bonded HSS (spheres joined in localised contact points) was analysed. A com­parison between analytical, numerical and experimental approaches is given in [10]. Fiedler et al. reported in [11] recent advances in the prediction of the thermal proper­ties of syntactic metallic hollow sphere structures. They described the application of the Finite Element and Lat­tice Monte Carlo Method in the case of syntactic peri­odic and random hollow sphere structures. Manufactur­ing techniques for perforated hollow sphere structures are actually under development [12] and we use an idealised model structure to clarify basic effects of the perforation on physical properties.a) b)Figure 1: Single hollow spheres: a) closed surface (com­mon configuration); b) with perforated surface (new de­velopment); (©by Glatt GmbH, Dresden, Germany).In the present study, perforated hollow spheres arranged in several types of periodic pattern are numerically in­vestigated based on the finite element method (unit cell approach) and the effective thermal conductivity is com-pared to structures without holes. The influence of the hole diameter and the geometry on the effective con­ductivity is analysed within a parametric computational study.2 ManufacturingA powder metallurgy based manufacturing process en­ables the production of metallic hollow spheres of defined geometry [13]. This technology brings a significant re­duction in costs in comparison to earlier applied galvanic methods and all materials suitable for sintering can be ap­plied. EPS (expanded polystyrol) spheres are coated with a metal powdcr-binder suspension by turbulence coating. The green spheres produced can either be sintered sep­arately to manufacture single hollow spheres or be pre­compacted and sintered in bulk (cf. Fig. 2) thus creat- ingsintering necks between adjacent spheres [14].Figure 2: Manufacturing processes of single hollowspheres and hollow sphere structures [12].Depending on the parameters of the sintering process the micro-porosity of the sintered cell wall can be adjusted. In a subsequent debindering process, the EPS spheres are pyrolised. The increase of the carbon content of the sintered metal by the diffusion of the incinerated binder and polymer causes degradation of mechanical properties and corrosion resistance. Special reducing processes are required to reduce this effect [15]. Various joining tech­nologies such as sintering, soldering and adhering can be used to assemble the single hollow spheres to interdepen­dent structures [16, 17].3 Perforated Hollow Sphere StructureThe major idea of introducing a perforation, i.e. holes of circular cross section in the sphere shells, is to make the inner sphere surface and volume usable. In this paper, the holes are defined in such a way that the largest pos­sible hole can be located between the sintering areas in a primitive cubic arrangement. In subsequent steps, the size of this initial hole was gradually reduced in the other models with different types of arrangements. Our mod­eling approach is restricted to the simple case where theholes do not intersect with the sintering area. The addi­tional inner surface may be used for chemical reactions in the case of catalysts or the additional inner volume may have a positive influence on the dissipation of e.g. acoustic waves. A simplified arrangement of perforated hollow sphere structures in a primitive cubic pattern is shown in Fig. 3.hollow spheresintering areaFigure 3: Schematic representation of primitive cubic sphere arrangements: a) sintered classical hollow sphere structure; b) sintered perforated hollow spheres struc­ture.Looking at a specific example where perforated spheres (outer diameter 2.7 mm and shell thickness of 0.1 mm, outer hole radius 0.68 mm) are considered, the available free surface increases by ~57% and the free volume in­crease by '■'■'259% while the bulk volume or the weight is reduced by ~33%. In addition, we can state that the porosity is only increased by ~3.9%.3.1 M odeling o f Perforated  H ollow Sphere StructuresSince the aim of this research is to highlight the differ­ence between perforated and entire spheres (cf. Fig. 3) to investigate the influence of perforation on the heat con­ductivity, different CAD models were used. By means of an image processing software, a series of micrographs was analysed and geometrical values for different model structures are derived [18], cf. Fig. 4,8 and Tab. 1.Figure 4: a) Photo of a sintered hollow sphere structure; b) Schematic sketch of a sintered hollow sphere structure.In order to reduce the size of the system of equations, the so-called unit-cell approach is chosen where the com­putational structure is reduced from a larger or infinite amount of randomly arranged spheres to a single unit-cell which is commonly based on typical spacc groups as in crystallography [19]. This is a first estimation for a finite element analysis that can give reasonable results. Com­mon unit cell approaches are based on primitive cubic, face-centered, body-centered or hexagonal arrangements, cf. Fig. 5.4 Finite Element Approach: Geometry and MeshFinite element analysis demands the discretisation of the geometries by subdividing them into geometrically simple finite elements. In the scope of these analyses, primar­ily three-dimensional geometries are considered. For the three dimensional meshing, tetrahedral or hexahedral el­ements can be employed. Earlier investigations [20, 21] have shown that hexahedral elements yield superior per­formance. Therefore, the geometry of the structures is discretised based on regular hexahedral elements. This approach is much more time-consuming, but it is impor­tant in order to achieve more accurate simulations. The meshing of the micro-structure of MHSS requires the de­composition of these complex geometries into simple sub­geometries. Based on the obtained fragments, meshing algorithms can be applied in order to obtain regular hexa­hedral meshes [22]. Therefore, these sub-geometries have to exhibit a principally cubical shape. Figure 6 demon­strates the decomposition of a syntactic face centered cu­bic geometry.Figure 5: Unit cell models for sintered spheres: a) primi­tive cubic (entire UC); b) hexagonal (entire UC); c) face- centered cubic (one-eight of UC); d) body-centered (one- eight of UC).Figure 6: Generation of the computational FCC model (one-twenty four of UC): a) CAD geometry; b) subdivi­sion of the solid geometry in simple parts; c) entire mesh.In order to minimise the number of fragments, symme­tries are exploited. In the example shown, only ^  of the final model sub-geometries need to be meshed and, by mirroring rest of fragments are obtained which assemble to the target geometry (cf. Fig. 7).5 Finite Element M ethod and Basics of Heat TransferA typical finite element mesh including boundary condi­tions is shown in Fig. 8. Constant temperature bound­ary conditions (rl \  = const.) are prescribed at the facesDimension Value [mm]Figure 7: Generation of the computational FCC model: a) one-twenty four of UC; b) one-twelve of UC; c) one- eight of UC.of the connective elements on the left and right side of the model. In order to generate a heat flux through the structure, only T\ /  T2 must be fulfilled.Dsbsdst2.660.61.470.1Table 1: Sphere dimensions.k = QA  A T (1)The area A  of the unit cell and the spatial distance Ay are given by the geometry (cf. Fig. 8), respectively the temperature gradient AT = T2 — T\ by the boundary conditions, only the heat flux Q remains to be deter­mined. This is done by summing up all nodal values of the reaction heat flux with a user subroutine at the left or right face where a temperature boundary condition is prescribed. Within the relevant temperature range, the contribution of thermal radiation to the heat transfer is low [4]. Furthermore, contributions from gaseous con­duction and convection are neglected. The whole model consists of 51056 elements in the case of the largest hole with a radius of 0.68 mm. Subsequent models where the hole diameter was reduced to 75, 50 and 25% where gen­erated. Figure 9 summarises the influence of mesh density on outcomes results for heat conductivity, calculated by the finite element code. One can observe that increasing the mesh density results in a stable value of conductiv­ity. Thus, as indicated in the figure, a choice of 27292 elements is reasonable to calculate the heat conductivity.T, = const. 0.06— QFigure 8: Finite element mesh and boundary conditions of a sintered perforated primitive cubic unit cell.According to the assumed symmetry and simplifications (no radiative or convective effects), the ratio of the heat flux perpendicular to all remaining surfaces is zero. The thermal properties of the considered base material is for steel k Rt.e e \  = 0.05 W/mm-K (AISI 8000).The evaluation of the effective thermal conductivities within the finite element approach is based on Fourier’s law where the area-related conductivity is defined by0.05 -£  0.040.00-------- 0u ----------------- ---------chosen model— ©---------PC wilhvut IjuIc k, =  0.05 W/(.iim K)10000 20000 30000Total number of nodes40000Figure 9: Results of mesh refinement analysis.All the simulations were done with the commercial finite element package MSC.Marc ©(MSC Software Corpora­tion, Santa Ana, CA, USA).6 Average Values for Cellular MetalsThe average density (Eq. 2) is an common value to make a statement on the density, mass, weight and furtherlightweight parameters. Physical properties are com­monly described as a function of their average density[23]Q»ve — s^p ( )^where g3p is the density of the solid material from which the cells are made. It should be noted that the volume of the unit cell (UC) can be expressed as the sum of solid and free volume as:V\JC “  VL0 +  Vfree ( 3 )A further important field is the comparison between the experimental tests with a rather random arrangement of spheres and the simulation results with the idealised ar­rangement of spheres. Real hollow sphere structures show a rather random arrangement of the spheres. On the other hand, many modelling (i.e. this paper) approaches are based on periodic geometries and structures which possess different packing densities such as primitive cu­bic (pc), body-centered cubic (bcc), face-centered cubic (fee) and hexagonal closest (he) (cf. Fig. 5). Knowing the packing density of the real structure, one may assign the proper topology for modelling the structure or use the result to interpolate model calculations based on the simplified periodic structures [24]. The results are plot­ted over the average density to be comparable with the experimental results which are currently under develop­ment.7 ResultsFigure 10 shows the normalised conductivity over the av­erage density for sintered hollow sphere structures. It can be seen that the heat conductivity is increasing when the average density is increasing (the hole diameter is de­creasing). In addition, one can observe that the results for different types of arrangements including primitive cubic (PC), face centred cubic(FCC), body centred cu­bic (BCC) and hexagonal (HEX) have a linear behaviour depending on the average density. That means that the thermal conductivity is much more influenced by the av­erage density than by different types of arrangements.Also the influence of the hollow spheres5 wall thickness on the thermal conductivity has been studied for sintered PHSS. Decreasing the wall thickness reduces the thermal conductivity, e.g. when the wall thickness decreases from Q. 1 mm to 0.07 mm, the thermal conductivity is also de­creased. Table 2 summarises this reduction for different configuration.Average density (kg/dm3)Figure 10: Normalised effective thermal conductivity as a function of relative density for different hole arrange­ments in the case of sintered spheres.Hole diameter FCC PC BCC HEX0ds 31.96% 31.67% 31.8% 28.74%0.254 31.93% 31.76% 35.13% 37.50%0.54 31.86% 32.12% 33.16% 38.10%0.754 40.66% 32.38% 40.64% 32.49%1 4 29.55% 33.10% 53.16% 40.36%Table 2: Reduction of thermal conductivity when the wall thickness of the hollow sphere is reduced from 0.1 mm to 0.07 mm for different hole diameters (Id is the case where the largest possible hole is located between the sintering areas with rfs =  1.47 mm, cf. Fig. 8)8 OutlookIn the scope of this research project, the effective ther­mal conductivity of perforated hollow sphere structures has been numerically estimated and compared to configu­rations without holes for the case that hollow spheres are sintered together. Introducing holes in the sphere shells clearly reduces the thermal conductivity of the HSS. Also reducing of the wall thickness in hollow spheres reduce the thermal conductivity. In addition, it is true to say that the thermal conductivity is much more influenced by the average density than by different types of arrangements. The numerical approach has been based on unit cells study. Future investigations must clarify the influence of the hole arrangement and geometrical modifications.References[1] Lu, W., Zhao, C.Y., Tassou, S.A., “Thermal Anal­ysis on Metal-Foam Filled Heat Exchangers. Part I: Metal-Foam Filled Pipes, “ Int J Heat Mass Trun, V49, N15-16, pp. 2751-2761, 7/06.[2] Lu, T.J., Stone, H.A., Ashby, M.F., “Heat Transfer in Open-Cell Metal Foams,w, Acta Mater, V46, N10, pp. 3619-3635, 6/98.[3] Boomsma, K., Poulilcakos, D., Zwick, F., “Metal Foams as Compact High Performance Heat Ex­changers,w, Mech Mater, V35, N12, pp. 1161-1176, 12/03.[4] Lu, T J , r Chen, C.r “Thermal Transport and Fire Retardance Properties of Cellular Aluminium Al­loys,M Acta Mater, V47, N5, pp. 1469-1485, 3/99.[5] Schneider, L., Boehm, A., Korhammer, C., Scholl, R., Voigtsberger, B., Stephani, G., US Patent 6501784. (2002).[6] Baumeister, E., Klaeger, S., Kaldos, A., “Lightweight, Hollow-Sphere-Composite (HSC) Materials for Mechanical Engineering Applica­tions, “ J  Mater Process Tech, V155, Special Issue: Part 2 Sp. Iss. SI, pp. 1839-1846, 11/04.[7] Lu, K.T., Kou, H J,, “Combined Boundray avd In­ertia Effects for Fully-Developed Mixed Convection in a Vertical Chanel Enbeded in Porous-Media, “ Int Commun Heat Mass, V20, N3, pp. 333-345, 6/93.♦[8] Fiedler, T., Ochsner, A., “On the Thermal Con­ductivity of Adhesively Bonded and Sintered Hollow Sphere Structures (HSS),“ Mat Sci Forum V553, pp. 39-44, 8/07.[9] Fiedler, T., Ochsner, A .,“Influence of the Morphol­ogy of Joining on the Heat Transfer Properties of Pe­riodic Metal Hollow Sphere, “ Mat Sci Forum V553, pp. 39-44, 8/07.[10] Fiedler, T., SolorzanO, E., Ochsner, A., “Numerical and Experimental Analysis of the Thermal Conduc­tivity of Metallic Hollow Sphere Structures,“ Mater Lett, V62, N8-9, pp. 1204-1207, 3/08.[11] Fiedler, T., Ochsner, A., Belova, I.V., Murch, G.E., “Recent Advances in the Prediction of the Ther­mal Properties of Syntactic Metallic Hollow Sphere Structures/* Adv Eng Mater, V 10, N4, pp. 361-365, 4/08.[12] Glatt GmbH, Binzen, Germany, private communicar tion (20Q7).[13] Jackel, M., German Patent, 1987,724 (3),156.[14] Jackel, M., German Patent, 1983, 210(3),770.[15] Studnitzky, T., Andersen, O., Cellular Metals and Polymers, T^ans Tech Publications, 2005.[16] Rousset A., Bonino, J.P., Blottiere, Y., Rossignol, G., FVench Patent, 1987, 707(8),440;[17] Degischer, H.P., Kriszt, B., editors, Handbook of Cel­lular Metals, WILEY-VCH,2002.[18] Veyhl, C., Winkler, R., Merkel, M., ochsner, A., “Structural Characterisation of Diffusion-Bonded Hollow Sphere Structure, “ Defect and Diffusion Fo­rum, V280-281, pp. 85-96, 11/08.[19] de Graef, M., McHenry, M., Structure of Materials, An Introduction to Crystallography, Diffraction and Symmetry, Cambridge: Cambridge University Press; 1997.[20] Benzley, S.E., Perry, E„ Merkle, K., Clark, B., Sjaardema, G.F., “A comparison of allhexagonal and alltetrahedral finite element meshes for elastic and elastic-plastic analysis, “ Fourth International Mesh- ing Roundtable, , Albuquerque, New Mexico, pp. 179b 191, 10/95.[21] Fiedler, T., Sturm, B., Ochsner, A., Gracio, J., Kuhn, G.. “Modelling the mechanical behaviour of adhesively bonded and sintered hollow sphere struc­tures, “ Mech Compos Mater, V42, N6, pp. 559-570. 11/06.[22] Fiedler, T., PhD Thesis, University of Aveiro, Por­tugal, 2007.[23] Gibson, L.J., Ashby, M.F., Cellular solids: struc­tures & properties, Cambridge University Press, Cambridge, 1997.[24] Veyhl. C., Master Thesis, University of Applied Sci­ence Aallen, Germany, 2008.

http://eprints.utm.my/id/eprint/12247/1/MarkusMerkel2009_NumericalSimulationofSinteredPerforatedHollow.pdf

Numerical simulation of sintered perforated hollow sphere structures (PHSS) to investigate thermal conductivity

Abstract

Similar works

Full text

Available Versions

University of Newcastle's Digital Repository

Universiti Teknologi Malaysia Institutional Repository