Simulation of Heat Conduction in Complex Domains of Multi-material Composites using a Meshless Method

Naman Bartwal,Shantanu Shahane,Somnath Roy,Surya Pratap Vanka
DOI: https://doi.org/10.48550/arXiv.2209.02082
2022-09-06
Abstract:Several engineering applications involve complex materials with significant and discontinuous variations in thermophysical properties. These include materials for thermal storage, biological tissues with blood capillaries, etc. For such applications, numerical simulations must exercise care in not smearing the interfaces by interpolating variables across the interfaces of the subdomains. In this paper, we describe a high accuracy meshless method that uses domain decomposition and cloud-based interpolation of scattered data to solve the heat conduction equation in such situations. The polyharmonic spline (PHS) function with appended polynomial of prescribed degree is used for discretization. A flux balance condition is satisfied at the interface points and the clouds of interpolation points are restricted to be within respective domains. Compared with previously proposed meshless algorithms with domain decomposition, the cloud-based interpolations are numerically better conditioned, and achieve high accuracy through the appended polynomial. The accuracy of the algorithm is demonstrated in several two and three dimensional problems using manufactured solutions to the heat conduction equation with sharp discontinuity in thermal conductivity. Subsequently, we demonstrate the applicability of the algorithm to solve heat conduction in complex domains with practical boundary conditions and internal heat generation. Systematic computations with varying conductivity ratios, interpoint spacing and degree of appended polynomial are performed to investigate the accuracy of the algorithm.
Numerical Analysis
What problem does this paper attempt to address?
### What problems does this paper attempt to solve? This paper aims to solve the key problems encountered in heat conduction simulation in complex multi - material composite structures. Specifically, the paper focuses on complex materials with significant and discontinuous changes in thermophysical properties, such as composite materials for heat storage, biological tissues containing blood vessels, etc. These materials are characterized by large differences in properties such as thermal conductivity between different regions, sometimes even more than 100 times. #### Main challenges: 1. **Interface treatment**: In numerical simulation, the interfaces between different sub - domains must be carefully treated to avoid interface smearing caused by interpolating variables, thus affecting the accuracy of the results. 2. **High - precision requirement**: For problems with large - range and discontinuous property changes, traditional methods (such as the finite volume method) are difficult to achieve high precision, especially in cases where the property differences are extremely large. #### Solutions: The paper proposes an efficient algorithm based on the meshless method, using domain decomposition and cloud - based scattered point interpolation to solve the heat conduction equation. The main features of this method include: - **Polyharmonic spline (PHS) function**: Combined with polynomials of a specified order for discretization. - **Flux balance condition**: Satisfy the flux balance condition at the interface points to ensure the continuity of heat flow on both sides of the interface. - **Restricted interpolation point cloud**: The interpolation point cloud in each sub - domain is limited to that sub - domain, avoiding interpolation across interfaces. #### Research purposes: - **Verify algorithm accuracy**: Test the algorithm accuracy in two - dimensional and three - dimensional problems through manufactured solutions. - **Apply to practical problems**: Apply the algorithm to heat conduction problems in complex geometries with actual boundary conditions and internal heat generation. - **Systematic research**: Systematically study the accuracy of the algorithm by changing the thermal conductivity ratio, point spacing, and the order of the additional polynomial. In conclusion, this paper is committed to developing an efficient meshless method to accurately simulate the heat conduction phenomenon in multi - material composite structures with complex geometries and significant property differences.