A 66 line heat transfer finite element code to highlight the dual approach. (English) Zbl 1443.65192

Summary: The paper presents a compact Matlab implementation of a finite element code performing dual analysis of heat conduction problems. After a short presentation of the concerned variational principles, the conductivity matrices of the dual models are developed explicitly. In the next step, a Matlab procedure is written. It shows in 66 statements how it is possible to perform the dual analysis using the unique formalism based on a global conductivity matrix. This fully general method allows using the same finite element code to run the dual analyses. After showing the convergence properties, a second application is developed in order to show how easily the code can be modified. The last part of the paper, addressed to skilled readers, provides useful comments on the heat problem formulation. The Matlab procedure is proposed as a free code for educational purposes and can be extracted directly from the document using the standard copy and paste function.


65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
65M22 Numerical solution of discretized equations for initial value and initial-boundary value problems involving PDEs
80M10 Finite element, Galerkin and related methods applied to problems in thermodynamics and heat transfer


top.m; Matlab
Full Text: DOI


[1] Beckers, B.; Beckers, P., (Beckers, B., Radiative Simulation Methods, in Solar Energy at Urban Scale (2012), John Wiley and Sons, Inc.), 205-236, (Chapter 10)
[2] Beckers, B.; Beckers, P., Sky vault partition for computing daylight availability and shortwave energy budget on an urban scale, Light. Res. Technol., 46, 716-728 (2014)
[3] Planas, P.; Flett, D. C., ESARAD: From R&D to industrial utilization, ESA Bull., 93, 93, 61-67 (1998)
[5] Bramble, J. H.; Payne, L. E., A priori bounds in the first boundary value problem in elasticity, J. Res. Natl. Bur. Stand. Sec. B Math. Math. Phys., 65, 4 (1961) · Zbl 0104.35102
[6] Sigmund, O., A 99 line topology optimization code written in Matlab, Struct. Multidiscip. Optim., 21, 120-127 (2001)
[7] Lewis, R. W.; Nithiarasu, P.; Seetharamu, K. N., Fundamentals of the Finite Element Method for Heat and Fluid Flow, 356 (2004), John Wiley & Sons Ltd.
[8] Fraeijs de Veubeke, B. M.; Hogge, M. A., Dual analysis for heat conduction problems by finite elements, Internat. J. Numer. Methods Engrg., 5, 1, 65-82 (1972) · Zbl 0251.65061
[9] Courant, R.; Hilbert, D., Methods of Mathematical Physics Vol. 1, 575 (1953), Wiley Interscience: Wiley Interscience New York
[10] Beckers, P., Les fonctions de tension dans la méthode des éléments finis (1972), University of Liège: University of Liège Belgium, (Doctoral thesis)
[11] Debongnie, J.-F.; Zhong, H. G.; Beckers, P., Dual analysis with general boundary conditions, Comput. Methods Appl. Mech. Engrg., 122, 1-2, 183-192 (1995) · Zbl 0851.73057
[12] Ladevèze, P.; Maunder, E. A.W., A general methodology for recovering equilibrating finite element tractions and stress fields for plate and solid elements, Comput. Assist. Mech. Eng. Sci., 4, 533-548 (1997) · Zbl 0974.74060
[13] Maunder, E. A.W.; Zhong, H. G.; Beckers, P., A posteriori error estimators related to equilibrium defaults of finite element solutions for elastostatic problems, Finite Elem. Anal. Des., 26, 171-192 (1997) · Zbl 0917.73074
[14] Kempeneers, M.; Debongnie, J. F.; Beckers, P., Pure equilibrium tetrahedral finite elements for global error estimation by dual analysis, Internat. J. Numer. Methods Engrg., 81, 4, 513-536 (2010) · Zbl 1183.74284
[15] Pavot, S.; Florentin, E.; Rougeot, P.; Champaney, L., Finite element verification in the case of missing data, Finite Elem. Anal. Des., 64, 90-96 (2013) · Zbl 1282.74094
[16] Coorevits, P.; Ladevèze, P.; Pelle, J.-P., An automatic procedure with a control of accuracy for finite element analysis in 2D elasticity, Comput. Methods Appl. Mech. Engrg., 121, 91-120 (1995) · Zbl 0851.73056
[17] Southwell, R. V., On the analogues relating flexure and extension of flat plates, Quart. J. Mech. Appl. Math., 3, 3, 257-270 (1950) · Zbl 0041.10401
[18] Needham, T., Visual Complex Analysis, 610 (1997), Oxford University Press
[19] Beckers, B.; Beckers, P., Reconciliation of Geometry and Perception in Radiation Physics, 192 (2014), John Wiley and Sons
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. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.