Numerical solutions of multi-dimensional partial differential equations using an adaptive wavelet method

by Damrongsak Wirasaet In this work, we describe an adaptive wavelet method for the solution of timeindependent and time-dependent partial differential equations in d-dimensions. The method is based on d-dimensional interpolating wavelets constructed from tensor products of 1-D interpolating wavelets. The connection between interpolating wavelets and dyadic grid points, and the fact that wavelet amplitudes indicate the local regularity of solutions are used in the construction of a computational grid of irregular points. Operations, such as the wavelet transform, its inverse, and interpolation are performed efficiently. In the spatial discretization, the derivative approximation on the irregular grid is obtained by means of consistent finite differences. An extension of the adaptive method to problems defined on more complicated domains is achieved by a domain transformation technique. For time-independent problems, the method is tested on 2and 3-D Poisson and Helmholtz problems with exact manufactured solutions in order to numerically study the connection between the order of the wavelet, the order of finite difference, the threshold values, and the accuracy of numerical solutions. The combustion of a 2-D flame ball-vortex interaction is used as a test problem for the time-independent algorithm. Application of adaptive method to incompressible Navier-Stokes equations is accomplished through the use of the Chorin Damrongsak Wirasaet projection method for time discretization. We apply the algorithm to simulate the flow in the 2-D lid-driven cavity at moderate Reynolds numbers and in the 2-D differentially-heated cavity at high Rayleigh numbers, and in the 3-D differentially heated cavity for various values of Rayleigh numbers. It is found that numerical results, while requiring a relatively small number of degrees of freedom, are in good agreement with the most accurate results available in the literature.

[1]  Sigal Gottlieb,et al.  Spectral Methods , 2019, Numerical Methods for Diffusion Phenomena in Building Physics.

[2]  Wolfgang Dahmen,et al.  Fast computation of adaptive wavelet expansions , 2007, Numerische Mathematik.

[3]  Michael Griebel,et al.  Adaptive sparse grid multilevel methods for elliptic PDEs based on finite differences , 1998, Computing.

[4]  Samuel Paolucci,et al.  Wavelet‐based adaptive multiresolution computation of viscous reactive flows , 2006 .

[5]  Samuel Paolucci,et al.  The Application of an Adaptive Wavelet Method to the 3-D Natural-Convection Flow in a Differentially Heated Cavity , 2006 .

[6]  Wolfgang Dahmen,et al.  Convergence of Adaptive Wavelet Methods for Goal-Oriented Error Estimation , 2006 .

[7]  Samuel Paolucci,et al.  Adaptive Wavelet Method for Incompressible Flows in Complex Domains , 2005 .

[8]  Oleg V. Vasilyev,et al.  An Adaptive Wavelet Collocation Method for Fluid-Structure Interaction at High Reynolds Numbers , 2005, SIAM J. Sci. Comput..

[9]  Application of an Adaptive Wavelet Method to Natural-Convection Flow in a Differentially Heated Cavity , 2005 .

[10]  Takeo S. Saitoh,et al.  Benchmark solutions for natural convection in a cubic cavity using the high-order time–space method , 2004 .

[11]  Olivier Roussel,et al.  A conservative fully adaptive multiresolution algorithm for parabolic PDEs , 2003 .

[12]  Oleg V. Vasilyev,et al.  Solving Multi-dimensional Evolution Problems with Localized Structures using Second Generation Wavelets , 2003 .

[13]  Wolfgang Dahmen,et al.  Sparse Evaluation of Compositions of Functions Using Multiscale Expansions , 2003, SIAM J. Math. Anal..

[14]  Albert Cohen,et al.  Fully adaptive multiresolution finite volume schemes for conservation laws , 2003, Math. Comput..

[15]  Albert Cohen.  Adaptive Methods for PDE's Wavelets or Mesh Refinement? , 2002 .

[16]  Mark A. Christon,et al.  Computational predictability of time‐dependent natural convection flows in enclosures (including a benchmark solution) , 2002 .

[17]  Stefano Berrone,et al.  A REALIZATION OF A WAVELET GALERKIN METHOD ON NONTRIVIAL DOMAINS , 2002 .

[18]  Wolfgang Dahmen,et al.  Wavelet Least Squares Methods for Boundary Value Problems , 2001, SIAM J. Numer. Anal..

[19]  Y. Rastigejev.  Multiscale computations with a wavelet-adaptive algorithm , 2002 .

[20]  C. Pozrikidis.  A note on the regularization of the discrete Poisson-Neumann problem , 2001 .

[21]  I. Daubechies,et al.  Tree Approximation and Optimal Encoding , 2001 .

[22]  Jaroslav Mackerle,et al.  Error estimates and adaptive finite element methods , 2001 .

[23]  W. Dahmen.  Wavelet methods for PDEs — some recent developments , 2001 .

[24]  Wolfgang Dahmen,et al.  Adaptive Wavelet Schemes for Elliptic Problems - Implementation and Numerical Experiments , 2001, SIAM J. Sci. Comput..

[25]  Rosa Donat,et al.  Point Value Multiscale Algorithms for 2D Compressible Flows , 2001, SIAM J. Sci. Comput..

[26]  Wolfgang Dahmen,et al.  Appending boundary conditions by Lagrange multipliers: Analysis of the LBB condition , 2001, Numerische Mathematik.

[27]  Wolfgang Dahmen,et al.  Adaptive wavelet methods for elliptic operator equations: Convergence rates , 2001, Math. Comput..

[28]  P. Oswald.  On Preconditioners for Sparse Grid Discretizations , 2001 .

[29]  G. Labrosse,et al.  A FIRST INCURSION INTO THE 3D STRUCTURE OF NATURAL CONVECTION OF AIR IN A DIFFERENTIALLY HEATED CUBIC CAVITY, FROM ACCURATE NUMERICAL SOLUTIONS , 2000 .

[30]  Frank Koster.  A Proof of the Consistency of the Finite Difference Technique on Sparse Grids , 2000, Computing.

[31]  Albert Cohen,et al.  Wavelet adaptive method for second order elliptic problems: boundary conditions and domain decomposition , 2000, Numerische Mathematik.

[32]  Wolfgang Dahmen,et al.  Nonlinear functionals of wavelet expansions – adaptive reconstruction and fast evaluation , 2000, Numerische Mathematik.

[33]  Silvia Bertoluzza,et al.  On the adaptive computation of integrals of wavelets , 2000 .

[34]  Kai Schneider,et al.  Numerical simulation of a mixing layer in an adaptive wavelet basis , 2000 .

[35]  Michael Griebel,et al.  Adaptive Wavelet Solvers for the Unsteady Incompressible Navier-Stokes Equations , 2000 .

[36]  Wolfgang Dahmen,et al.  Composite wavelet bases for operator equations , 1999, Math. Comput..

[37]  Mats Holmström,et al.  Solving Hyperbolic PDEs Using Interpolating Wavelets , 1999, SIAM J. Sci. Comput..

[38]  Johan Walden.  Filter Bank Methods for Hyperbolic PDEs , 1999 .

[39]  Reinhard H Stephan Dahlke.  Adaptive Wavelet Methods for Saddle Point Problems , 1999 .

[40]  Philippe Angot,et al.  A penalization method to take into account obstacles in incompressible viscous flows , 1999, Numerische Mathematik.

[41]  Karsten Urban,et al.  Adaptive Wavelet Galerkin Methods on Distorted Domains: Setup of the Algebraic System , 1999 .

[42]  Rémi Abgrall,et al.  Multiresolution Representation in Unstructured Meshes , 1998 .

[43]  Anthony Randriamampianina,et al.  An improved projection scheme applied to pseudospectral methods for the incompressible Navier-Stokes equations , 1998 .

[44]  Pascal Monasse,et al.  Orthonormal wavelet bases adapted for partial differential equations with boundary conditions , 1998 .

[45]  O. Botella,et al.  BENCHMARK SPECTRAL RESULTS ON THE LID-DRIVEN CAVITY FLOW , 1998 .

[46]  Masud Behnia,et al.  From onset of unsteadiness to chaos in a differentially heated square cavity , 1998, Journal of Fluid Mechanics.

[47]  Wei Cai,et al.  An Adaptive Spline Wavelet ADI (SW-ADI) Method for Two-Dimensional Reaction-Diffusion Equations , 1998 .

[48]  Rob Stevenson.  Piecewise linear (pre-)wavelets on non-uniform meshes , 1998 .

[49]  Michael Griebel,et al.  Numerical Simulation in Fluid Dynamics: A Practical Introduction , 1997 .

[50]  Silvia Bertoluzza.  An adaptive collocation method based on interpolating wavelets , 1997 .

[51]  O. Vasilyev,et al.  A Fast Adaptive Wavelet Collocation Algorithm for Multidimensional PDEs , 1997 .

[52]  Claudio Canuto,et al.  The wavelet element method. Part II: Realization and additional features in 2D and 3D , 1997 .

[53]  Claudio Canuto,et al.  The wavelet element method. Part I: Construction and analysis. , 1997 .

[54]  Wolfgang Dahmen,et al.  Stable multiscale bases and local error estimation for elliptic problems , 1997 .

[55]  Jochen Fröhlich,et al.  An Adaptive Wavelet-Vaguelette Algorithm for the Solution of PDEs , 1997 .

[56]  W. Dahmen.  Wavelet and multiscale methods for operator equations , 1997, Acta Numerica.

[57]  G. Chiavassa,et al.  On the Effective Construction of Compactly Supported Wavelets Satisfying Homogenous Boundary Conditions on the Interval. , 1997 .

[58]  Jianzhong Wang,et al.  Adaptive multiresolution collocation methods for initial boundary value problems of nonlinear PDEs , 1996 .

[59]  O. Vasilyev,et al.  A Dynamically Adaptive Multilevel Wavelet Collocation Method for Solving Partial Differential Equations in a Finite Domain , 1996 .

[60]  Silvia Bertoluzza.  Adaptive wavelet collocation method for the solution of Burgers equation , 1996 .

[61]  Philippe H. Tchamitchian.  Wavelets, functions, and operators , 1996 .

[62]  Wolfgang Dahmen,et al.  Multiscale methods for pseudodifferential equations , 1996 .

[63]  S. Paolucci,et al.  A Petrov-Galerkin method for the direct simulation of fully enclosed flows , 1996 .

[64]  Barna L. Bihari,et al.  Multiresolution Schemes for Conservation Laws with Viscosity , 1996 .

[65]  S. Bertoluzza,et al.  A Wavelet Collocation Method for the Numerical Solution of Partial Differential Equations , 1996 .

[66]  Wim Sweldens,et al.  The lifting scheme: a construction of second generation wavelets , 1998 .

[67]  Silvia Bertoluzza.  Adaptive wavelet collocation for the solution of steady-state equations , 1995, Defense, Security, and Sensing.

[68]  Björn Sjögreen,et al.  Numerical experiments with the multiresolution scheme for the compressible Euler equations , 1995 .

[69]  Leland Jameson,et al.  On the spline-based wavelet differentiation matrix , 1995 .

[70]  Wim Sweldens,et al.  The Construction and Application of Wavelets in Numerical Analysis , 1995 .

[71]  Wim Sweldens,et al.  An Overview of Wavelet Based Multiresolution Analyses , 1994, SIAM Rev..

[72]  William D. Henshaw,et al.  A Fourth-Order Accurate Method for the Incompressible Navier-Stokes Equations on Overlapping Grids , 1994 .

[73]  R.J.A. Janssen.  Instabilities in natural-convection flows in cavities , 1994 .

[74]  Leland Jameson,et al.  ON THE WAVELET OPTIMIZED FINITE DIFFERENCE METHOD , 1994 .

[75]  Silvia Bertoluzza,et al.  Wavelet Methods for the Numerical Solution of Boundary Value Problems on the Interval , 1994 .

[76]  Leland Jameson,et al.  ON THE DAUBECHIES-BASED WAVELET DIFFERENTIATION MATRIX , 1993 .

[77]  I. Daubechies,et al.  Wavelets on the Interval and Fast Wavelet Transforms , 1993 .

[78]  Naoki Saito,et al.  Multiresolution representations using the autocorrelation functions of compactly supported wavelets , 1993, IEEE Trans. Signal Process..

[79]  Wei Cai,et al.  ADAPTIVE WAVELET COLLOCATION METHODS FOR INITIAL VALUE BOUNDARY PROBLEMS OF NONLINEAR PDE''S , 1993 .

[80]  I. Weinreich,et al.  Wavelet-Galerkin methods: An adapted biorthogonal wavelet basis , 1993 .

[81]  Sam Qian,et al.  Wavelets and the Numerical Solution of Partial Differential Equations , 1993 .

[82]  Using wavelets to solve the Burgers equation: A comparative study. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[83]  Jinchao Xu,et al.  Galerkin-wavelet methods for two-point boundary value problems , 1992 .

[84]  S. Jaffard.  Wavelet methods for fast resolution of elliptic problems , 1992 .

[85]  Ingrid Daubechies,et al.  Ten Lectures on Wavelets , 1992 .

[86]  GermanyNumerische Mathematik,et al.  Multilevel Preconditioning , 1992 .

[87]  C.-C. Jay Kuo,et al.  A Domain Decomposition Preconditioner Based on a Change to a Multilevel Nodal Basis , 1991, SIAM J. Sci. Comput..

[88]  P. LeQuéré,et al.  Accurate solutions to the square thermally driven cavity at high Rayleigh number , 1991 .

[89]  Bakhtier Farouk,et al.  A numerical study of three-dimensional natural convection in a differentially heated cubical enclosure , 1991 .

[90]  Jacques Liandrat,et al.  Resolution of the 1D regularized Burgers equation using a spatial wavelet approximation , 1990 .

[91]  Neil A. Kelson.  The laminar boundary layer regime for natural convection of air in a square cavity , 1990 .

[92]  Youcef Saad,et al.  A Basic Tool Kit for Sparse Matrix Computations , 1990 .

[93]  A. Latto,et al.  Compactly supported wavelets and the numerical solution of Burgers' equation , 1990 .

[94]  Gilles Deslauriers,et al.  Symmetric iterative interpolation processes , 1989 .

[95]  William F. Mitchell,et al.  A comparison of adaptive refinement techniques for elliptic problems , 1989, TOMS.

[96]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[97]  S. Paolucci,et al.  Transition to chaos in a differentially heated vertical cavity , 1989, Journal of Fluid Mechanics.

[98]  M. Cross,et al.  Numerical grid generation—foundations and applications , 1986 .

[99]  S. Paolucci,et al.  Natural convection in an enclosed vertical air layer with large horizontal temperature differences , 1986, Journal of Fluid Mechanics.

[100]  S. Vanka.  Block-implicit multigrid solution of Navier-Stokes equations in primitive variables , 1986 .

[101]  Harry Yserentant,et al.  On the multi-level splitting of finite element spaces , 1986 .

[102]  G. Lauriat,et al.  A new formulation of the SADI method for the prediction of natural convection flows in cavities , 1985 .

[103]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[104]  P. Gresho,et al.  Finite element simulations of thermally induced convection in an enclosed cavity , 1980 .

[105]  W. Roger Briley,et al.  Numerical method for predicting three-dimensional steady viscous flow in ducts , 1974 .

[106]  G. Beylkin,et al.  On the representation of operators in bases of compactly supported wavelets , 1992 .