Coupled parabolic-hyperbolic system often appears in the studies of thermoelasticity, magnetoelasticity, biological problems and radiation hydrodynamics with high temperature. In this paper, we investigate a problem of thermal diffusion in elastic body with moving boundary. Three numerical methods, two uncoupled and one coupled, all with quadratic convergence order in time and space, are presented and compared in relation to the execution time. Tables and figures of the approximate solution are shown to verify the efficiency and feasibility of the proposed method. In addition, we show that the numerical results are consistent with the theoretical results. The approximate numerical solutions are calculated using the finite element method in spatial variable and finite difference method in time. To compare different numerical schemes, we show the convergence rates through numerical simulations with calculated error and the execution time.
KEYWORDScoupled and uncoupled method, finite elements method, moving boundary, numerical simulation, thermal diffusion
INTRODUCTIONThe thermal diffusion in elastic body is governed by a coupled system of hyperbolic-parabolic equations. In [1], the theory was developed using the coupled thermoelastic model. Uniqueness and reciprocity theorems for the generalized problem of thermal diffusion in elastic body in isotropic media, was proved in [2], under restrictive assumptions on the elastic coefficients. In these works, one