A local high-order doubly asymptotic open boundary for diffusion in a semi-infinite layer

A high-order open boundary for transient diffusion in a semi-infinite homogeneous layer is developed. The method of separation of variables is used to derive a relationship between the modal function and the flux at the near field/far field boundary in the Fourier domain. The resulting equation in terms of the modal impedance coefficient is solved by expanding the latter into a doubly asymptotic series of continued fractions. As a result, the open boundary condition in the Fourier domain is represented by a system of algebraic equations in terms of i@w. This corresponds to a system of fractional differential equations of degree @a=0.5 in the time-domain. This temporally global formulation is transformed into a local description by introducing internal variables. The resulting local high-order open boundary condition is highly accurate, as is demonstrated by a number of heat transfer examples. A significant gain in accuracy is obtained in comparison with existing singly-asymptotic formulations at no additional computational cost.

[1]  Thomas L. Geers,et al.  Doubly asymptotic approximations for transient elastodynamics , 1997 .

[2]  Kai Diethelm,et al.  An investigation of some nonclassical methods for the numerical approximation of Caputo-type fractional derivatives , 2008, Numerical Algorithms.

[3]  Alex H. Barnett,et al.  A fast numerical method for time-resolved photon diffusion in general stratified turbid media , 2004 .

[4]  Laurence Halpern,et al.  Absorbing boundary conditions for diffusion equations , 1995 .

[5]  I. Podlubny.  Fractional differential equations , 1998 .

[6]  Xiaonan Wu,et al.  Adaptive absorbing boundary conditions for Schrödinger-type equations: Application to nonlinear and multi-dimensional problems , 2007, J. Comput. Phys..

[7]  Christophe Besse,et al.  A Review of Transparent and Artificial Boundary Conditions Techniques for Linear and Nonlinear Schrödinger Equations , 2008 .

[8]  Philip Underwood,et al.  Doubly asymptotic, boundary-element analysis of dynamic soil-structure interaction. Topical report 1 June 1975--31 March 1978 , 1978 .

[9]  Frank J. Rizzo,et al.  A method of solution for certain problems of transient heat conduction , 1970 .

[10]  Zhi-Zhong Sun,et al.  The stability and convergence of a difference scheme for the Schrödinger equation on an infinite domain by using artificial boundary conditions , 2006, J. Comput. Phys..

[11]  Houde Han,et al.  A class of artificial boundary conditions for heat equation in unbounded domains , 2002 .

[12]  T. Hagstrom.  Boundary conditions at outflow for a problem with transport and diffusion , 1987 .

[13]  Luis E. Suarez,et al.  A Comparison of Numerical Methods Applied to a Fractional Model of Damping Materials , 1999 .

[14]  Eric Dubach.  Artificial boundary conditions for diffusion equations: numerical study , 1996 .

[15]  R. A. Adey,et al.  Numerical methods in thermal problems: Edited by R. W. Lewis and K. Morgan Proceedings of the First International Conference held at University College, Swansea, 2–6 July 1979 , 1980 .

[16]  Houde Han,et al.  Exact and approximating boundary conditions for the parabolic problems on unbounded domains , 2002 .

[17]  Leslie Greengard,et al.  On the numerical solution of the heat equation I: Fast solvers in free space , 2007, J. Comput. Phys..

[18]  John P. Wolf,et al.  Consistent lumped‐parameter models for unbounded soil: Physical representation , 1991 .

[19]  Thomas L. Geers,et al.  Doubly asymptotic approximations for transient motions of submerged structures , 1978 .

[20]  Song-Ping Zhu,et al.  Solving linear diffusion equations with the dual reciprocity method in Laplace space , 1994 .

[21]  Chunxiong Zheng,et al.  Exact artificial boundary conditions for problems with periodic structures , 2008, J. Comput. Phys..

[22]  S. Tsynkov.  Numerical solution of problems on unbounded domains. a review , 1998 .

[23]  Johannes Tausch,et al.  A fast method for solving the heat equation by layer potentials , 2007, J. Comput. Phys..

[24]  Nobuyuki Shimizu,et al.  Numerical Algorithm for Dynamic Problems Involving Fractional Operators , 1998 .

[25]  Christophe Besse,et al.  Unconditionally stable discretization schemes of non-reflecting boundary conditions for the one-dimensional Schrödinger equation , 2003 .

[26]  Leslie Greengard,et al.  Fast evaluation of nonreflecting boundary conditions for the Schrödinger equation in one dimension , 2004 .

[27]  O. Agrawal,et al.  A Numerical Scheme for Dynamic Systems Containing Fractional Derivatives , 2002 .

[28]  James M. Kelly,et al.  Application of fractional derivatives to seismic analysis of base‐isolated models , 1990 .

[29]  António Tadeu,et al.  Boundary element method analyses of transient heat conduction in an unbounded solid layer containing inclusions , 2004 .

[30]  Gary F. Dargush,et al.  Application of the boundary element method to transient heat conduction , 1991 .

[31]  L. Debnath,et al.  Integral Transforms and Their Applications, Second Edition , 2006 .

[32]  John P. Wolf,et al.  Consistent lumped‐parameter models for unbounded soil: Frequency‐independent stiffness, damping and mass matrices , 1991 .

[33]  Peter Ruge,et al.  Consistent Modelling of Infinite Beams by Fractional Dynamics , 2004 .

[34]  Chongmin Song,et al.  A continued‐fraction‐based high‐order transmitting boundary for wave propagation in unbounded domains of arbitrary geometry , 2008, International Journal for Numerical Methods in Engineering.

[35]  Xavier Antoine,et al.  Numerical accuracy of a Padé‐type non‐reflecting boundary condition for the finite element solution of acoustic scattering problems at high‐frequency , 2005 .

[36]  Carlos Alberto Brebbia,et al.  The dual reciprocity boundary element formulation for diffusion problems , 1987 .

[37]  Xiaonan Wu,et al.  Convergence of difference scheme for heat equation in unbounded domains using artificial boundary conditions , 2004 .

[38]  Eduardo Kausel,et al.  Heat diffusion in layered media via the thin‐layer method , 2009 .

[39]  Jian-Fei Lu,et al.  Wave field simulation for heterogeneous porous media with singular memory drag force , 2005 .

[40]  Carolin Birk,et al.  A continued-fraction approach for transient diffusion in unbounded medium , 2009 .

[41]  Thomas L. Geers.  Singly and Doubly Asymptotic Computational Boundaries , 1998 .

[42]  Gary H. Brooke,et al.  Rational square-root approximations for parabolic equation algorithms , 1997 .

[43]  Lothar Gaul,et al.  On a critique of a numerical scheme for the calculation of fractionally damped dynamical systems , 2006 .

[44]  Sondipon Adhikari,et al.  Direct time-domain integration method for exponentially damped linear systems , 2004 .

[45]  W. Ames.  Mathematics in Science and Engineering , 1999 .

[46]  D. Givoli.  Numerical Methods for Problems in Infinite Domains , 1992 .

[47]  Thomas L. Geers,et al.  THIRD-ORDER DOUBLY ASYMPTOTIC APPROXIMATIONS FOR COMPUTATIONAL ACOUSTICS , 2000 .

[48]  Gao Lin,et al.  High‐order doubly asymptotic open boundaries for scalar wave equation , 2009 .

[49]  P. Ruge,et al.  Treatment of dynamic systems with fractional derivatives without evaluating memory-integrals , 2002 .

[50]  Carlos Alberto Brebbia,et al.  A formulation of the boundary element method for axisymmetric transient heat conduction , 1981 .

[51]  Richard Paul Shaw,et al.  An integral equation approach to diffusion , 1974 .

[52]  Thomas Hagstrom,et al.  Asymptotic expansions and boundary conditions for time-dependent problems , 1986 .