Bandwidth-based mesh adaptation in multiple dimensions

Spectral methods are becoming increasingly prevalent in solving time-varying partial differential equations due to their fast convergence properties. However, they typically use regular computational meshes that do not account for spatially varying resolution requirements. This can significantly inc...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics 2018-10, Vol.371, p.651-662
Main Authors: Wise, Elliott S., Cox, Ben T., Treeby, Bradley E.
Format: Article
Language:English
Subjects:
Adaptation
Adaptive moving mesh method
Bandwidth
Bandwidths
Boundary layer
Burgers equation
Computational physics
Computer simulation
Convergence
Curvature
Decay rate
Finite element method
Mathematical models
Monge-Ampere equation
Nodes
Optical wireless
Partial differential equations
Pseudospectral method
Spectral methods
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Spectral methods are becoming increasingly prevalent in solving time-varying partial differential equations due to their fast convergence properties. However, they typically use regular computational meshes that do not account for spatially varying resolution requirements. This can significantly increase the overall grid density when resolution requirements vary sharply over the modelled domain. Moving mesh methods offer a remedy for this, by allowing the position of mesh nodes to adapt to the simulated model solution. In this paper, a mesh specification is presented that is based on a local measure of the spatial bandwidth of the model solution. This addresses the rate of decay of the model solution's frequency components by producing high-sampling rates when this decay is slow. The spatial bandwidth is computed using a combination of the original solution and its Riesz transformed counterparts. It is then integrated into a Fourier spectral moving mesh method, using the parabolic Monge–Ampère equation for mesh control. This method is used to solve a multidimensional version of the viscous Burgers equation, and a heterogeneous advection equation. The performance of bandwidth-based mesh adaptation is compared with arclength- and curvature-based adaptation, and against a static mesh. These numerical experiments show that the bandwidth-based approach produces superior convergence rates, and hence requires fewer mesh nodes for a given level of solution accuracy.
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2018.06.009