Abstract
In this paper we discuss spectral approximations of the Poisson equation in deformed quadrilateral domains. High order polynomial approximations are used for both the solution and the representation of the geometry. Following an isoparametric approach, the four edges of the computational domain are first parametrized using high order polynomial interpolation. Transfinite interpolation is then used to construct the mapping from the square reference domain to the physical domain. Through a series of numerical examples we show the importance of representing the boundary of the domain in a careful way; the choice of interpolation points along the edges of the physical domain may significantly effect the overall discretization error. One way to ensure good interpolation points along an edge is based on the following criteria: (i) the points should be on the exact curve; (ii) the derivative of the exact curve and the interpolant should coincide at the internal points along the edge. Following this approach, we demonstrate that the discretization error for the Poisson problem may decay exponentially fast even when the boundary has low regularity.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References
Bernardi, C., Maday, Y.: Spectral methods. In: Ciarlet, P.G., Lions, J.L. (eds.) Handbook of Numerical Analysis. Volume V: Techniques of Scientific Computing (Part 2), pp. 209–485. Elsevier, Amsterdam (1997)
Berrut, J.-P., Trefethen, L.N.: Barycentric Lagrange interpolation. SIAM Rev. 46(3), 501–517 (2004)
Canuto, C., Hussaini, M.Y., Quarteroni, A., Zang, T.A.: Spectral Methods, Fundamentals in Single Domains. Springer, Berlin (2006)
Canuto, C., Hussaini, M.Y., Quarteroni, A., Zang, T.A.: Spectral Methods, Evolution to Complex Geometries and Applications to Fluid Dynamics. Springer, Berlin (2007)
Carey, G.F., Dinh, H.T.: Grading functions and mesh redistribution. SIAM J. Numer. Anal. 22(5), 1028–1040 (1985)
Degen, W.L.F.: Geometric Hermite interpolation. Comput. Aided Geom. Des. 22(7), 573–592 (2005). In memoriam Josef Hoschek
Deville, M.O., Fischer, P.F., Mund, E.H.: High-Order Methods for Incompressible Fluid Flow. Cambridge University Press, Cambridge (2002)
Fornberg, B.: A Practical Guide to Pseudospectral Methods. Cambridge (1998)
Gordon, W., Hall, C.: Construction of curvilinear co-ordinate systems and applications to mesh generation. Int. J. Numer. Methods Eng. 7(4), 461–477 (1973)
Hesthaven, J., Warburton, T.: Nodal Discontinuous Galerkin Methods. Springer, Berlin (2008)
Karniadakis, G., Sherwin, S.: Spectral/hp Element Methods for CFD. Oxford (2005)
Maday, Y., Patera, A.T.: Spectral element methods for the Navier–Stokes equations. In: Noor, A.K., Oden, J.T. (eds.) State of the Art Surveys in Computational Mechanics, pp. 71–143. ASME, New York (1989)
Mørken, K., Scherer, K.: A general framework for high-accuracy parametric interpolation. Math. Comput. 66(217), 237–260 (1997)
O’Neill, B.: Elementary Differential Geometry. Academic Press, New York (2006)
Rababah, A.: High order approximation method for curves. Comput. Aided Geom. Des. 12(1), 89–102 (1995)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Bjøntegaard, T., Rønquist, E.M. & Tråsdahl, Ø. Spectral Approximation of Partial Differential Equations in Highly Distorted Domains. J Sci Comput 52, 603–618 (2012). https://doi.org/10.1007/s10915-011-9561-8
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-011-9561-8