Efficient high order method for differential equations in unbounded domains using generalized coordinate transformation

Faisal Mumtaz, Hamed Saidaoui, Fahhad Alharbi

Research output: Contribution to journalArticle

Abstract

Many physical phenomena are localized in space, where physical quantities of interest extend to large distances (effectively to infinity) where they vanish smoothly. Thus, its modeling necessitates the imposition of vanishing boundary conditions (BCs) to the formulated differential equations (DEs). Numerical solutions for such DEs can be complicated as they involve operations on physical quantities expanding infinitely in space. Currently, many approaches have been implemented in this context, where they mostly use techniques like domain truncation, using functions intrinsic to unbounded domains, and domain transformation. However, since the actual solutions have a decaying nature, solving such problems with high accuracy and efficiency while taking into consideration a wide range of decaying rates and oscillation frequencies is a major challenge. In this work, we present an efficient high-order method, based on Galerkin method, to solve DEs in unbounded domains using properly designed domain transformations where the physical unbounded domain (−∞,∞) is mapped into a computational bounded domain (0,1). The mapped sine series is then used to solve the differential equation in the computational domain (0,1) where the vanishing quantities at infinities are mapped to 0 (corresponding to physical −∞) and 1 (corresponding to physical ∞). The designed transformations maintain the orthogonality in both physical and computational spaces. Moreover, using sine series basis allows – in most cases – analytical integrations. For verification, the proposed method is applied to solve 1) Helmholtz equation in one-dimension, where the exact solutions for various cases (decaying exponential or algebraic function) are known, and 2) three-dimensional Schrödinger wave equation with harmonic oscillator. For the Helmholtz problem, the error in the approximate solution was less than 10−40 using 100 basis functions, while for the second problem, the error is less than 10−25 using the same number of bases. All results and the analyses demonstrate how the properly designed mapped basis sets can be used to develop efficient high-order methods to solve DEs with vanishing boundary conditions and with wide range of decaying rates and oscillation frequencies simultaneously.

Original languageEnglish
Pages (from-to)275-289
Number of pages15
JournalJournal of Computational Physics
Volume381
DOIs
Publication statusPublished - 15 Mar 2019

Fingerprint

coordinate transformations
Differential equations
differential equations
sine series
Boundary conditions
Helmholtz equation
infinity
Galerkin methods
Wave equations
boundary conditions
oscillations
Helmholtz equations
problem solving
Galerkin method
orthogonality
exponential functions
harmonic oscillators
wave equations
approximation

Keywords

  • Coordinate transformation
  • Helmholtz equation
  • High order method
  • Simple harmonic oscillator
  • Spectral method
  • Unbounded domains

ASJC Scopus subject areas

  • Physics and Astronomy (miscellaneous)
  • Computer Science Applications

Cite this

Efficient high order method for differential equations in unbounded domains using generalized coordinate transformation. / Mumtaz, Faisal; Saidaoui, Hamed; Alharbi, Fahhad.

In: Journal of Computational Physics, Vol. 381, 15.03.2019, p. 275-289.

Research output: Contribution to journalArticle

@article{287906bf5dd445689b6cd16f14bd9483,
title = "Efficient high order method for differential equations in unbounded domains using generalized coordinate transformation",
abstract = "Many physical phenomena are localized in space, where physical quantities of interest extend to large distances (effectively to infinity) where they vanish smoothly. Thus, its modeling necessitates the imposition of vanishing boundary conditions (BCs) to the formulated differential equations (DEs). Numerical solutions for such DEs can be complicated as they involve operations on physical quantities expanding infinitely in space. Currently, many approaches have been implemented in this context, where they mostly use techniques like domain truncation, using functions intrinsic to unbounded domains, and domain transformation. However, since the actual solutions have a decaying nature, solving such problems with high accuracy and efficiency while taking into consideration a wide range of decaying rates and oscillation frequencies is a major challenge. In this work, we present an efficient high-order method, based on Galerkin method, to solve DEs in unbounded domains using properly designed domain transformations where the physical unbounded domain (−∞,∞) is mapped into a computational bounded domain (0,1). The mapped sine series is then used to solve the differential equation in the computational domain (0,1) where the vanishing quantities at infinities are mapped to 0 (corresponding to physical −∞) and 1 (corresponding to physical ∞). The designed transformations maintain the orthogonality in both physical and computational spaces. Moreover, using sine series basis allows – in most cases – analytical integrations. For verification, the proposed method is applied to solve 1) Helmholtz equation in one-dimension, where the exact solutions for various cases (decaying exponential or algebraic function) are known, and 2) three-dimensional Schr{\"o}dinger wave equation with harmonic oscillator. For the Helmholtz problem, the error in the approximate solution was less than 10−40 using 100 basis functions, while for the second problem, the error is less than 10−25 using the same number of bases. All results and the analyses demonstrate how the properly designed mapped basis sets can be used to develop efficient high-order methods to solve DEs with vanishing boundary conditions and with wide range of decaying rates and oscillation frequencies simultaneously.",
keywords = "Coordinate transformation, Helmholtz equation, High order method, Simple harmonic oscillator, Spectral method, Unbounded domains",
author = "Faisal Mumtaz and Hamed Saidaoui and Fahhad Alharbi",
year = "2019",
month = "3",
day = "15",
doi = "10.1016/j.jcp.2018.12.030",
language = "English",
volume = "381",
pages = "275--289",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - Efficient high order method for differential equations in unbounded domains using generalized coordinate transformation

AU - Mumtaz, Faisal

AU - Saidaoui, Hamed

AU - Alharbi, Fahhad

PY - 2019/3/15

Y1 - 2019/3/15

N2 - Many physical phenomena are localized in space, where physical quantities of interest extend to large distances (effectively to infinity) where they vanish smoothly. Thus, its modeling necessitates the imposition of vanishing boundary conditions (BCs) to the formulated differential equations (DEs). Numerical solutions for such DEs can be complicated as they involve operations on physical quantities expanding infinitely in space. Currently, many approaches have been implemented in this context, where they mostly use techniques like domain truncation, using functions intrinsic to unbounded domains, and domain transformation. However, since the actual solutions have a decaying nature, solving such problems with high accuracy and efficiency while taking into consideration a wide range of decaying rates and oscillation frequencies is a major challenge. In this work, we present an efficient high-order method, based on Galerkin method, to solve DEs in unbounded domains using properly designed domain transformations where the physical unbounded domain (−∞,∞) is mapped into a computational bounded domain (0,1). The mapped sine series is then used to solve the differential equation in the computational domain (0,1) where the vanishing quantities at infinities are mapped to 0 (corresponding to physical −∞) and 1 (corresponding to physical ∞). The designed transformations maintain the orthogonality in both physical and computational spaces. Moreover, using sine series basis allows – in most cases – analytical integrations. For verification, the proposed method is applied to solve 1) Helmholtz equation in one-dimension, where the exact solutions for various cases (decaying exponential or algebraic function) are known, and 2) three-dimensional Schrödinger wave equation with harmonic oscillator. For the Helmholtz problem, the error in the approximate solution was less than 10−40 using 100 basis functions, while for the second problem, the error is less than 10−25 using the same number of bases. All results and the analyses demonstrate how the properly designed mapped basis sets can be used to develop efficient high-order methods to solve DEs with vanishing boundary conditions and with wide range of decaying rates and oscillation frequencies simultaneously.

AB - Many physical phenomena are localized in space, where physical quantities of interest extend to large distances (effectively to infinity) where they vanish smoothly. Thus, its modeling necessitates the imposition of vanishing boundary conditions (BCs) to the formulated differential equations (DEs). Numerical solutions for such DEs can be complicated as they involve operations on physical quantities expanding infinitely in space. Currently, many approaches have been implemented in this context, where they mostly use techniques like domain truncation, using functions intrinsic to unbounded domains, and domain transformation. However, since the actual solutions have a decaying nature, solving such problems with high accuracy and efficiency while taking into consideration a wide range of decaying rates and oscillation frequencies is a major challenge. In this work, we present an efficient high-order method, based on Galerkin method, to solve DEs in unbounded domains using properly designed domain transformations where the physical unbounded domain (−∞,∞) is mapped into a computational bounded domain (0,1). The mapped sine series is then used to solve the differential equation in the computational domain (0,1) where the vanishing quantities at infinities are mapped to 0 (corresponding to physical −∞) and 1 (corresponding to physical ∞). The designed transformations maintain the orthogonality in both physical and computational spaces. Moreover, using sine series basis allows – in most cases – analytical integrations. For verification, the proposed method is applied to solve 1) Helmholtz equation in one-dimension, where the exact solutions for various cases (decaying exponential or algebraic function) are known, and 2) three-dimensional Schrödinger wave equation with harmonic oscillator. For the Helmholtz problem, the error in the approximate solution was less than 10−40 using 100 basis functions, while for the second problem, the error is less than 10−25 using the same number of bases. All results and the analyses demonstrate how the properly designed mapped basis sets can be used to develop efficient high-order methods to solve DEs with vanishing boundary conditions and with wide range of decaying rates and oscillation frequencies simultaneously.

KW - Coordinate transformation

KW - Helmholtz equation

KW - High order method

KW - Simple harmonic oscillator

KW - Spectral method

KW - Unbounded domains

UR - http://www.scopus.com/inward/record.url?scp=85060269540&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85060269540&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2018.12.030

DO - 10.1016/j.jcp.2018.12.030

M3 - Article

VL - 381

SP - 275

EP - 289

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -