Solving Hyperbolic PDEs Using Interpolating Wavelets
algorithmBurgers' equationfinite differencesfirst-order linear advection-diffusion equationhyperbolic initial boundary value problemshyperbolic PDEsinterpolating wavelet transforminterpolatory subdivision schemeinverse sparse wavelet transformRunge-Kutta methodsparse point representationwavelet representation
Initial value problems for first-order hyperbolic systems (35L45) KdV equations (Korteweg-de Vries equations) (35Q53) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Method of lines for initial value and initial-boundary value problems involving PDEs (65M20) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Numerical methods for wavelets (65T60)
The wavelet representation of a function comprises a finite linear combination of scaling functions and wavelets. The representation may be compressed by deleting wavelet coefficients less than some threshold \(\epsilon\), to leave \(N_s\) significant coefficients, with \(N_s\) generally much smaller than \(N\), the original number of coefficients. Such a procedure has the advantage of leading to significant savings in storage and computing time. In this paper, the author reports on a procedure for exploiting the sparsity of wavelet representations in solving time-dependent hyperbolic partial differential equations. The method uses a sparse point representation, which is based on an interpolating wavelet transform using polynomial interpolation on dyadic grids. Centred finite differences are used to approximate spatial derivatives, and a Runge-Kutta method is used for time integration. A detailed description is provided of an interpolatory subdivision scheme, which is an algorithm through which function values \(f_{j,k}\), defined on the set of dyadic grids \(V_j = \{ x_{j,k} \in R: x_{j,k} = 2^{-j}k\), \(k \in Z \}\), may be extended to all points in \(V_{j+1}\). Essentially, values at the even-numbered points \(x_{j+1,2k}\) already exist, while values at the odd-numbered points are obtained by polynomial interpolation from the values at even-numbered points. Used recursively, the algorithm generates function values on a fine grid, given values on a coarse grid. The reverse process, of computing function values on a fine grid, given those on a coarse grid, entails computing the difference between the known function value and the value predicted by interpolation on the coarser grid, to obtain the wavelet coefficients \(d_{j,k}\). This process provides an algorithm for computing the full wavelet representation from function values on a fine grid, or in other words, the interpolating wavelet transform. The process for obtaining the inverse transform is straightforward. Due to the one-to-one correspondence between wavelet coefficients and point values, the process of transforming \(N_s\) wavelet coefficients to \(N_s\) point values is well-defined. The process is referred to as the inverse sparse wavelet transform leading to a sparse point representation (SPR). The SPR is used to solve approximately hyperbolic initial boundary value problems on domains in \(R^n\) that are \(n\)-cubes. The algorithm entails transforming and thresholding the initial value function, extending the SPR in space and frequency, and discretising to obtain a system of ordinary differential equations, which are solved by a Runge-Kutta method. A number of examples is treated. First, the author examines the errors in approximating various test functions in one dimension by SPRs. Next, he considers the quality of approximations arising after multiplication and differentiation of functions. He then solves the first-order linear advection-diffusion equation with initial conditions in the form of a superposition of a sine function and a Gaussian, and Burgers' equation with initial condition a sine function. Corresponding two-dimensional versions of the one-dimensional examples are also presented. It is found that the proposed method outperforms a finite difference method on a uniform grid, in the sense that it is accurate, and faster.
- Multiscale cell-based coarsening for discontinuous problems
- Space-time adaptive multiresolution methods for hyperbolic conservation laws: Applications to compressible Euler equations
- ENO adaptive method for solving one-dimensional conservation laws
- High-order compact finite-difference scheme for singularly-perturbed reaction-diffusion problems on a new mesh of Shishkin type
- Fast visualization of volume emissions using conservative subdivision.
- Adaptive multiresolution approach for two-dimensional PDEs.
- The existence and uniqueness of solution to wavelet collocation
- An adaptive wavelet collocation method for solving optimal control of elliptic variational inequalities of the obstacle type
- Multiresolution based adaptive schemes for second order hyperbolic PDEs in elastodynamic problems
- Adaptive multiresolution approach for solution of hyperbolic PDEs
- On a wavelet-based method for the numerical simulation of wave propagation
- Interpolating wavelet collocation method of time dependent Maxwell's equations: Characterization of electrically large optical waveguide discontinuities.
- Adaptive wavelet collocation methods for image segmentation using TV-Allen-Cahn type models
- Review of wavelet methods for the solution of reaction-diffusion problems in science and engineering
- Wave propagation in generalized thermo-poro-elastic media via wavelet-based cell-adaptive central high resolution schemes using UNO limiters
- Modeling the long-range wave propagation by a split-step wavelet method
- Multiscale polynomial-based high-order central high resolution schemes
- Multiresolution-based adaptive central high resolution schemes for modeling of nonlinear propagating fronts
- Multiresolution symplectic scheme for wave propagation in complex media
- Fully adaptive multiresolution schemes for strongly degenerate parabolic equations with discontinuous flux
- Simultaneous space-time adaptive wavelet solution of nonlinear parabolic differential equations
- Simulation of two-dimensional elastodynamic problems using a new adaptive physics-based method
- A modified wavelet approximation for multi-resolution AWCM in simulating nonlinear vibration of MDOF systems
- An adaptive wavelet method for nonlinear partial differential equations with applications to dynamic damage modeling
- Comparison of adaptive multiresolution and adaptive mesh refinement applied to simulations of the compressible Euler equations
- Approximate solution of modified Camassa-Holm and Degasperis-Procesi equations using wavelet optimized finite difference method
- scientific article; zbMATH DE number 2134532 (Why is no real title available?)
- Anisotropy in wavelet-based phase field models
- A three-step wavelet Galerkin method for parabolic and hyperbolic partial differential equations
- WAMR: an adaptive wavelet method for the simulation of compressible reacting flow. Part I: Accuracy and efficiency of algorithm
- A wavelet Taylor-Galerkin method for parabolic and hyperbolic partial differential equations
- scientific article; zbMATH DE number 5538074 (Why is no real title available?)
- Improvement in computational efficiency of Euler equations via a modified sparse point representation method
- Filter Bank Methods for Hyperbolic PDEs
- Solving Multi-dimensional Evolution Problems with Localized Structures using Second Generation Wavelets
- Massively parallel simulations of binary black hole intermediate-mass-ratio inspirals
- Six-Dimensional Adaptive Simulation of the Vlasov Equations Using a Hierarchical Basis
- Highly accurate Coiflet wavelet-homotopy solution of Jeffery-Hamel problem at extreme parameters
- Solving shallow water equations by the sparse point representation method
- A Wavelet-Adaptive Method for Multiscale Simulation of Turbulent Flows in Flying Insects
- A wavelet method for hyperbolic partial differential equations
- TIME ACCURATE FAST THREE-STEP WAVELET-GALERKIN METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS
- A local pseudo arc-length method for hyperbolic conservation laws
- WAVELET ADAPTIVE MULTIRESOLUTION REPRESENTATION: APPLICATIONS TO VISCOUS MULTISCALE FLOW SIMULATIONS
- An adaptive multilevel wavelet collocation method for elliptic problems
- Second-generation wavelet collocation method for the solution of partial differential equations
- An efficient adaptive wavelet method for pricing time-fractional American option variational inequality
- Review of solution methodologies for structural analysis of composites
- Simulating 2D waves propagation in elastic solid media using wavelet based adaptive method
- Multiscale central high-resolution schemes with different types of extended slope limiters on wavelet-based adapted cells
- Automatic adaptive and targeted localized refinement in transient structural dynamic simulations using the multiresolution finite wavelet domain method
- Interpolatory wavelets for manifold-valued data
- Wavelet optimized finite difference method with non-static re-gridding
- Adaptive multiresolution WENO schemes for multi-species kinematic flow models
- Grid structure impact in sparse point representation of derivatives
This page was built for publication: Solving Hyperbolic PDEs Using Interpolating Wavelets
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4702388)