×

High-order accurate conjugate heat transfer solutions with a finite volume method in anisotropic meshes with application in polymer processing. (English) Zbl 07767262

Summary: The computational modeling has become an indispensable tool to support the engineering design and control of many industrial processes. In polymer processing applications, the simulation of conjugate heat transfer phenomena is critical as rigorous temperature control is necessary to ensure that the produced parts meet the required specifications. In this article, a high-order accurate finite volume method is proposed to improve the numerical accuracy and the computational efficiency of conjugate heat transfer simulations with application in polymer processing. Anisotropic meshes are investigated to significantly reduce the number of unknowns in convection-dominated problems where the higher temperature variations occur perpendicularly to curved boundaries and interfaces. The reconstruction for off-site data method based on polynomial reconstructions is employed to fulfill the prescribed boundary and interface conditions solely using polygonal meshes to avoid the limitations of curved mesh approaches. A code verification benchmark based on manufactured solutions proves that the proposed method provides a fourth-order of convergence and is computationally more cost-effective than the classical second-order of convergence. Moreover, meshes with higher aspect ratios improve the calculation efficiency but suffer from a small accuracy penalty due to conditioning deterioration. Comparison with the classical finite volume discretization techniques, widely implemented in commercial and open-source software packages, is also provided. The applicability and performance of the proposed method are further supported with a practical case study for the sheet extrusion cooling stage.
{© 2021 John Wiley & Sons Ltd.}

MSC:

65Nxx Numerical methods for partial differential equations, boundary value problems
76Mxx Basic methods in fluid mechanics
76Dxx Incompressible viscous fluids

Software:

OpenFOAM
Full Text: DOI

References:

[1] RauwendaalC, Gonzalez‐NunezR, RodrigueD. Polymer processing: extrusion. Encyclopedia of Polymer Science and Technology. John Wiley & Sons; 2017:1‐67.
[2] RosatoDV. Extruding Plastics: A Practical Processing Handbook. Springer Science & Business Media; 1998.
[3] BairdDG, ColliasDI. Polymer Processing: Principles and Design. John Wiley & Sons; 2014.
[4] Abdel‐BaryEM (ed.), ed. Handbook of Plastic Films I. Smithers Rapra Publishing; 2003.
[5] VlachopoulosJ, PolychronopoulosN, TanifujiS, MüllerJP. Chapter 4: Flat film and sheet dies. In: CarneiroOS (ed.), NóbregaJM (ed.), eds. Design of Extrusion Forming Tools. Smithers Rapra; 2012:113‐140.
[6] MacauleyNJ, Harkin‐JonesEMA, MurphyWR. The influence of extrusion parameters on the mechanical properties of polypropylene sheet. Polym Eng Sci. 1998;38(4):662‐673.
[7] SollogoubC, FelderE, DemayY, AgassantJF, DeparisP, MiklerN. Thermomechanical analysis and modeling of the extrusion coating process. Polym Eng Sci. 2008;48(8):1634‐1648.
[8] CarneiroOS, NóbregaJM. Chapter 1: Main issues in the design of extrusion tools. In: CarneiroOS (ed.), NóbregaJM (ed.), eds. Design of Extrusion Forming Tools. Smithers Rapra; 2012:1‐36.
[9] WhitefieldML. Sheet extrusion for today’s products. J Past Film Sheet. 1990;6:10‐16.
[10] NizamiJ. Stability analysis and controller design for polymer sheet extrusion. J Vib Control. 2000;6:1083‐1105.
[11] SmithDE. Design sensitivity analysis and optimization for polymer sheet extrusion and mold filling processes. Int J Numer Meth Eng. 2003;57:1381‐1411. · Zbl 1062.76509
[12] YangH. Conjugate thermal simulation for sheet extrusion die. Polym Eng Sci. 2014;54(3):682‐694.
[13] LinP, JaluriaY. Heat transfer and solidification of polymer melt flow in a channel. Polym Eng Sci. 1997;37(7):1247‐1258.
[14] LinP, JaluriaY. Conjugate transport in polymer melt flow through extrusion dies. Polym Eng Sci. 1997;37(9):1582‐1595.
[15] LinP, JaluriaY. Conjugate thermal transport in the channel of an extruder for non‐Newtonian fluids. Int J Heat Mass Transf. 1998;41:3239‐3253.
[16] WapperomP, HassagerO. Numerical simulation of wire‐coating: the influence of temperature boundary conditions. Polym Eng Sci. 1999;39(10):2007‐2018.
[17] NóbregaJM, CarneiroOS, CovasJA, PinhoFT, OliveiraPJ. Design of calibrators for extruded profiles Part I: modeling the thermal interchanges. Polym Eng Sci. 2004;44(12):2216‐2228.
[18] MousseauP, DelaunayD, LefèvreN. Analysis of the heat transfer in PVC profiles during the extrusion calibration/cooling step. Int Polym Process. 2009;24:122‐132.
[19] CarneiroOS, NóbregaJM, MotaAR, SilvaC. Prototype and methodology for the characterization of the polymer‐calibrator interface heat transfer coefficient. Polym Test. 2013;32(6):1154‐1161.
[20] MurashovMV, PaninSD. Numerical modelling of contact heat transfer problem with work hardened rough surfaces. Int J Heat Mass Transf. 2015;90:72‐80.
[21] CocoA, RussoG. Second order finite‐difference ghost‐point multigrid methods for elliptic problems with discontinuous coefficients on an arbitrary interface. J Comput Phys. 2018;361:299‐330. · Zbl 1422.65306
[22] XingY, SongL, HeX, QiuC. A generalized finite difference method for solving elliptic interface problems. Math Comput Simul. 2020;178:109‐124. · Zbl 1515.65272
[23] VarsakelisS, MarichalY. Numerical approximation of elliptic interface problems via isoparametric finite element methods. Comput Math Appl. 2014;68:1945‐1962. · Zbl 1369.65150
[24] BadiaS, VerdugoF, MarínAF. The aggregated unfitted finite element method for elliptic problems. Comput Methods Appl Mech Eng. 2018;336:533‐553. · Zbl 1440.65175
[25] GuoR, LinT. An immersed finite element method for elliptic interface problems in three dimensions. J Comput Phys. 2020;414:109478. · Zbl 1440.65207
[26] MuL, ZhangX. An immersed weak Galerkin method for elliptic interface problems. J Comput Appl Math. 2019;362:471‐483. · Zbl 1422.65404
[27] LatigeM, CollinT, GalliceG. A second order Cartesian finite volume method for elliptic interface and embedded Dirichlet problems. Comput Fluids. 2013;83:70‐76. · Zbl 1290.65100
[28] CaoY, WangB, XiaK, WeiG. Finite volume formulation of the MIB method for elliptic interface problems. J Comput Appl Math. 2017;321:60‐77. · Zbl 1366.65096
[29] MarquesF, ClainS, MachadoGJ, MartinsB, CarneiroOS, NóbregaJM. A new energy conservation scheme for the numeric study of the heat transfer in profile extrusion calibration. Heat Mass Transf. 2017;53:2901‐2913.
[30] ZhouH, ShengZ, YuanG. A finite volume method preserving maximum principle for the diffusion equations with imperfect interface. Appl Numer Math. 2020;158:314‐335. · Zbl 1451.65124
[31] GuoK, LiL, XiaoG, AuYeungN, MeiR. Lattice Boltzmann method for conjugate heat and mass transfer with interfacial jump conditions. Int J Heat Mass Transf. 2015;88:306‐322.
[32] LuJH, LeiHY, DaiCS. A unified thermal lattice Boltzmann equation for conjugate heat transfer problem. Int J Heat Mass Transf. 2018;126:1275‐1286.
[33] PatankarSV. Numerical Heat Transfer and Fluid Flow. Series in Computational Methods in Mechanics and Thermal Sciences. 1st ed.Hemisphere Publishing; 1980. · Zbl 0521.76003
[34] WangZJ, FidkowskiK, AbgrallR, et al. High‐order CFD methods: current status and perspective. Int J Numer Methods Fluids. 2013;72:811‐845. · Zbl 1455.76007
[35] WangZJ. High‐order computational fluid dynamics tools for aircraft design. Phil Trans R Soc A. 2014;372:20130318.
[36] CaiZ, ThornberB. An internal penalty discontinuous Galerkin method for simulating conjugate heat transfer in a closed cavity. Int J Numer Methods Fluids. 2018;87(3):134‐159.
[37] CangianiA, GeorgoulisE, SabawiY. Adaptive discontinuous Galerkin methods for elliptic interface problems. Math Comput. 2018;87(314):2675‐2707. · Zbl 1397.65252
[38] HuynhLNT, NguyenNC, PeraireJ, KhooBC. A high‐order hybridizable discontinuous Galerkin method for elliptic interface problems. Int J Numer Meth Eng. 2013;93(2):183‐200. · Zbl 1352.65513
[39] WangB, KhooBC. Hybridizable discontinuous Galerkin method (HDG) for Stokes interface flow. J Comput Phys. 2013;247:262‐278. · Zbl 1349.76075
[40] QiuW, SolanoM, VegaP. A high order HDG method for curved‐interface problems via approximations from straight triangulations. J Sci Comput. 2016;69(3):1384‐1407. · Zbl 1371.65121
[41] PaipuriM, TiagoC, Fernández‐MéndezS. Coupling of continuous and hybridizable discontinuous Galerkin methods: application to conjugate heat transfer problem. J Sci Comput. 2019;78:321‐350. · Zbl 1422.65407
[42] SunH, DarmofalDL. An adaptive simplex cut‐cell method for high‐order discontinuous Galerkin discretizations of elliptic interface problems and conjugate heat transfer problems. J Comput Phys. 2014;278:445‐468. · Zbl 1349.80034
[43] OjedaSM, SunH, AllmarasSR, DarmofalDL. An adaptive simplex cut‐cell method for high‐order discontinuous Galerkin discretizations of conjugate heat transfer problems. Int J Numer Meth Eng. 2017;110:350‐378. · Zbl 1375.80014
[44] MuL, WangJ, WeiG, YeX, ZhaoS. Weak Galerkin methods for second order elliptic interface problems. J Comput Phys. 2013;250:106‐125. · Zbl 1349.65472
[45] MuL, WangJ, YeX, ZhaoS. A new weak Galerkin finite element method for elliptic interface problems. J Comput Phys. 2016;325:157‐173. · Zbl 1380.65383
[46] GuzmánJ, SánchezMA, SarkisM. Higher‐order finite element methods for elliptic problems with interfaces. ESAIM Math Model Numer Anal. 2016;50(5):1561‐1583. · Zbl 1353.65120
[47] HuangP, WuH, XiaoY. An unfitted interface penalty finite element method for elliptic interface problems. Comput Meth Appl Mech Eng. 2017;323:439‐460. · Zbl 1439.74422
[48] LehrenfeldC, ReuskenA. Analysis of a high‐order unfitted finite element method for elliptic interface problems. IMA J Numer Anal. 2017;38(3):1351‐1387. · Zbl 1462.65193
[49] JiH, ChenJ, LiZ. A high‐order source removal finite element method for a class of elliptic interface problems. Appl Numer Math. 2018;130:112‐130. · Zbl 1397.65270
[50] WuH, XiaoY. An unfitted hp‐interface penalty finite element method for elliptic interface problems. J Comput Math. 2019;37(3):316‐339. · Zbl 1449.65326
[51] CheungJ, PeregoM, BochevP, GunzburgerM. Optimally accurate high‐order finite element methods for polytopial approximations of domains with smooth boundaries. Math Comput. 2019;88:2187‐2219. · Zbl 1417.65186
[52] CheungJ, GunzburgerM, BochevP, PeregoM. An optimally convergent higher‐order finite element coupling method for interface and domain decomposition problems. Results Appl Math. 2020;6:100094. · Zbl 1453.65427
[53] ZhouYC, ZhaoS, FeigM, WeiGW. High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources. J Comput Phys. 2006;213(1):1‐30. · Zbl 1089.65117
[54] ZhongX. A new high‐order immersed interface method for solving elliptic equations with imbedded interface of discontinuity. J Comput Phys. 2007;225(1):1066‐1099. · Zbl 1343.65130
[55] AngelovaIT, VulkovLG. High‐order finite difference schemes for elliptic problems with intersecting interfaces. Appl Math Comput. 2007;187(2):824‐843. · Zbl 1120.65331
[56] MarquesAN, NaveJ‐C, RosalesRR. High order solution of Poisson problems with piecewise constant coefficients and interface jumps. J Comput Phys. 2017;335:497‐515. · Zbl 1380.65224
[57] LiuJ‐K, ZhengZ‐S. Efficient high‐order immersed interface methods for heat equations with interfaces. Appl Math Mech. 2014;35(9):1189‐1202. · Zbl 1298.65124
[58] ElhaddadM, ZanderN, BogN, et al. Multi‐level hp‐finite cell method for embedded interface problems with application in biomechanics. Int J Numer Meth Biomed Eng. 2018;34(4):e2951.
[59] LehrenfeldC. High order unfitted finite element method on level set domains using isoparametric mappings. Comput Meth Appl Mech Eng. 2016;300:716‐733. · Zbl 1425.65168
[60] BoscheriW, DumbserM. High order accurate direct arbitrary‐Lagrangian‐Eulerian ADER‐WENO finite volume schemes on moving curvilinear unstructured meshes. Comput Fluids. 2016;136:48‐66. · Zbl 1390.76398
[61] GaoX, HuangS, CuiM, et al. Element differential method for solving general heat conduction problems. Int J Heat Mass Transf. 2017;115:882‐894.
[62] DuvigneauR. CAD‐consistent adaptive refinement using a NURBS‐based discontinuous Galerkin method. Int J Numer Methods Fluids. 2020;92(9):1096‐1117.
[63] PezzanoS, DuvigneauR. A NURBS‐based discontinuous Galerkin method for conservation laws with high‐order moving meshes. J Comput Phys. 2021;434:110093. · Zbl 07508514
[64] LiP, LiuJ, LinG, ZhangP, YangG. A NURBS‐based scaled boundary finite element method for the analysis of heat conduction problems with heat fluxes and temperatures on side‐faces. Int J Heat Mass Transf. 2017;113:764‐779.
[65] Ollivier‐GoochC, AltenaMV. A high‐order accurate unstructured mesh finite‐volume scheme for the advection‐diffusion equation. J Comput Phys. 2002;181(2):729‐752. · Zbl 1178.76251
[66] JalaliA, Ollivier‐GoochCF. Higher‐order finite volume solution reconstruction on highly anisotropic meshes. Proceedings of the 21st AIAA Computational Fluid Dynamics Conference; 2013.
[67] BoularasA, ClainS, BaudoinF. A sixth‐order finite volume method for diffusion problem with curved boundaries. Appl Math Model. 2017;42:401‐422. · Zbl 1443.65291
[68] Ollivier‐GoochC, NejatA, MichalakC. On obtaining high‐order finite‐volume solutions to the Euler equations on unstructured meshes. AIAA Paper. 2007;2007‐4464.
[69] Ollivier‐GoochC, NejatA, MichalakC. Obtaining and verifying high‐order unstructured finite volume solutions to the Euler equations. AIAA J. 2009;47(9):2105‐2120.
[70] MichalakC, Ollivier‐GoochC. Unstructured high‐order accurate finite volume solutions of the Navier‐Stokes equations. AIAA Paper. 2009;2009‐2954.
[71] CostaR, ClainS, LoubèreR, MachadoGJ. Very high‐order accurate finite volume scheme on curved boundaries for the two‐dimensional steady‐state convection‐diffusion equation with Dirichlet condition. Appl Math Model. 2018;54:752‐767. · Zbl 1480.65318
[72] CostaR, NóbregaJM, ClainS, MachadoGJ, LoubèreR. Very high‐order accurate finite volume scheme for the convection‐diffusion equation with general boundary conditions on arbitrary curved boundaries. Int J Numer Meth Eng. 2019;117(2):188‐220. · Zbl 07865099
[73] CostaR, NóbregaJM, ClainS, MachadoGJ. Very high‐order accurate polygonal mesh finite volume scheme for conjugate heat transfer problems with curved interfaces and imperfect contacts. Comput Meth Appl Mech Engrg. 2019;357:112560. · Zbl 1442.74007
[74] HablaF, FernandesC, MaierM, et al. Development and validation of a model for the temperature distribution in the extrusion calibration stage. Appl Therm Eng. 2016;100:538‐552.
[75] ClainS, MachadoGJ, NóbregaJM, PereiraRMS. A sixth‐order finite volume method for the convection‐diffusion problem with discontinuous coefficients. Comput Meth Appl Mech Eng. 2013;267:43‐64. · Zbl 1286.80003
[76] MarićT, HöpkenJ, MooneyKG. The OpenFOAM technology primer (OpenFOAM‐v2012.1). Zenodo; 2021.
[77] VälikangasT. Conjugate heat transfer in OpenFOAM. In: NilssonH (ed.), ed. Proceedings of CFD with OpenSource Software; Chalmers University of Technology; 2016.
[78] AbbassiME, LahayeD, VuikK. Modelling turbulent combustion coupled with conjugate heat transfer in OpenFOAM. In: VermolenFJ (ed.), VuikC (ed.), eds. Numerical Mathematics and Advanced Applications ENUMATH 2019. Lecture Notes in Computational Science and Engineering; Springer, Cham; 2019;139:1137‐1145. · Zbl 1478.80003
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.