A temporally local method for the numerical solution of transient diffusion problems in unbounded domains was proposed by combining the scaled boundary finite element method and a novel solution procedure for fractional differential equations. The scaled boundary finite element method was employed to model the unbounded domain. In the Fourier domain (ω), an equation of the stiffness matrix for diffusion representing the flux–temperature relationship at the discretized near field/far field interface was established. A continued-fraction solution in terms of √iω was obtained. By using the continued-fraction solution and introducing auxiliary variables, the flux–temperature relationship was formulated as a system of linear equations in √iω. In the time-domain, it was interpreted as a system of fractional differential equations of degree α = ½. To eliminate the computationally expensive convolution integral, the fractional differential equation was transformed to a system of first-order differential equations. Numerical examples of two- and three-dimensional heat conductions demonstrated the accuracy of the proposed method. The computational cost of both the temporally global and local approach for transient analysis was examined.

A Continued-Fraction Approach for Transient Diffusion in Unbounded Medium. C.Birk, C.Song: Computer Methods in Applied Mechanics and Engineering, 2009, 198[33-36], 2576-90