A Meshless Finite Difference Method for Conjugate Heat Conduction Problems

Author(s):  
Chandrashekhar Varanasi ◽  
Jayathi Y. Murthy ◽  
Sanjay Mathur

In recent years, there has been a great deal of interest in developing meshless methods for computational fluid dynamics (CFD) applications. In this paper, a meshless finite difference method is developed for solving conjugate heat transfer problems in complex geometries. Traditional finite difference methods (FDMs) have been restricted to an orthogonal or a body-fitted distribution of points. However, the Taylor series upon which the FDM is based is valid at any location in the neighborhood of the point about which the expansion is carried out. Exploiting this fact, and starting with an unstructured distribution of mesh points, derivatives are evaluated using a weighted least squares procedure. The system of equations that results from this discretization can be represented by a sparse matrix. This system is solved with an algebraic multigrid (AMG) solver. The implementation of Neumann, Dirichlet and mixed boundary conditions within this framework is described, as well as the handling of conjugate heat transfer. The method is verified through application to classical heat conduction problems with known analytical solutions. It is then applied to the solution of conjugate heat transfer problems in complex geometries, and the solutions so obtained are compared with more conventional unstructured finite volume methods. Metrics for accuracy are provided and future extensions are discussed.

2010 ◽  
Vol 132 (8) ◽  
Author(s):  
Chandrashekhar Varanasi ◽  
Jayathi Y. Murthy ◽  
Sanjay Mathur

A meshless finite difference method is developed for solving conjugate heat transfer problems. Starting with an arbitrary distribution of mesh points, derivatives are evaluated using a weighted least-squares procedure. The resulting system of algebraic equations is sparse and is solved using an algebraic multigrid method. The implementation of the Neumann, Dirichlet, and mixed boundary conditions within this framework is described. For conjugate heat transfer problems, continuity of the heat flux and temperature are imposed on mesh points at multimaterial interfaces. The method is verified through application to classical heat conduction problems with known analytical solutions. It is then applied to the solution of conjugate heat transfer problems in complex geometries, and the solutions so obtained are compared with more conventional unstructured finite volume methods. The method improves on existing meshless methods for conjugate heat conduction by eliminating spurious oscillations previously observed. Metrics for accuracy are provided and future extensions are discussed.


Author(s):  
Imam Basuki ◽  
C Cari ◽  
A Suparmi

<p class="Normal1"><strong><em>Abstract: </em></strong><em>Partial Differential Equations (PDP) Laplace equation can be applied to the heat conduction. Heat conduction is a process that if two materials or two-part temperature material is contacted with another it will pass heat transfer. Conduction of heat in a triangle shaped object has a mathematical model in Cartesian coordinates. However, to facilitate the calculation, the mathematical model of heat conduction is transformed into the coordinates of the triangle. PDP numerical solution of Laplace solved using the finite difference method. Simulations performed on a triangle with some angle values α and β</em></p><p class="Normal1"><strong><em> </em></strong></p><p class="Normal1"><strong><em>Keywords:</em></strong><em>  heat transfer, triangle coordinates system.</em></p><p class="Normal1"><em> </em></p><p class="Normal1"><strong>Abstrak</strong> Persamaan Diferensial Parsial (PDP) Laplace  dapat diaplikasikan pada persamaan konduksi panas. Konduksi panas adalah suatu proses yang jika dua materi atau dua bagian materi temperaturnya disentuhkan dengan yang lainnya maka akan terjadilah perpindahan panas. Konduksi panas pada benda berbentuk segitiga mempunyai model matematika dalam koordinat cartesius. Namun untuk memudahkan perhitungan, model matematika konduksi panas tersebut ditransformasikan ke dalam koordinat segitiga. Penyelesaian numerik dari PDP Laplace diselesaikan menggunakan metode beda hingga. Simulasi dilakukan pada segitiga dengan beberapa nilai sudut  dan  </p><p class="Normal1"><strong> </strong></p><p class="Normal1"><strong>Kata kunci :</strong> perpindahan panas, sistem koordinat segitiga.</p>


2015 ◽  
Vol 20 (4) ◽  
pp. 733-756
Author(s):  
S.K. Deb Nath ◽  
N.K. Peyada

Abstract In the present study, we have developed a code using Matlab software for solving a rectangular aluminum plate having void, notch, at different boundary conditions discretizing a two dimensional (2D) heat conduction equation by the finite difference technique. We have solved a 2D mixed boundary heat conduction problem analytically using Fourier integrals (Deb Nath et al., 2006; 2007; 2007; Deb Nath and Ahmed, 2008; Deb Nath, 2008; Deb Nath and Afsar, 2009; Deb Nath and Ahmed, 2009; 2009; Deb Nath et al., 2010; Deb Nath, 2013) and the same problem is also solved using the present code developed by the finite difference technique (Ahmed et al., 2005; Deb Nath, 2002; Deb Nath et al., 2008; Ahmed and Deb Nath, 2009; Deb Nath et al., 2011; Mohiuddin et al., 2012). To verify the soundness of the present heat conduction code results using the finite difference method, the distribution of temperature at some sections of a 2D heated plate obtained by the analytical method is compared with those of the plate obtained by the present finite difference method. Interpolation technique is used as an example when the boundary of the plate does not pass through the discretized grid points of the plate. Sometimes hot and cold fluids are passed through rectangular channels in industries and many types of technical equipment. The distribution of temperature of plates including notches, slots with different temperature boundary conditions are studied. Transient heat transfer in several pure metallic plates is also studied to find out the required time to reach equilibrium temperature. So, this study will help find design parameters of such structures.


Author(s):  
Lucas Peixoto ◽  
Ane Lis Marocki ◽  
Celso Vieira Junior ◽  
Viviana Mariani

Sign in / Sign up

Export Citation Format

Share Document