Research Papers

On Solving One-Dimensional Partial Differential Equations With Spatially Dependent Variables Using the Wavelet-Galerkin Method

[+] Author and Article Information
Simon Jones

Assistant Professor
Department of Mechanical Engineering,
Rose-Hulman Institute of Technology,
Terre Haute, IN 47803
e-mail: jones5@rose-hulman.edu

Mathias Legrand

Assistant Professor
Department of Mechanical Engineering,
McGill University,
Montreal, QC H3C 1PA, Canada
e-mail: mathias.legrand@mcgill.ca

1Corresponding author.

Manuscript received August 29, 2012; final manuscript received January 19, 2013; accepted manuscript posted February 12, 2013; published online August 21, 2013. Editor: Yonggang Huang.

J. Appl. Mech 80(6), 061012 (Aug 21, 2013) (7 pages) Paper No: JAM-12-1428; doi: 10.1115/1.4023637 History: Received August 29, 2012; Revised January 19, 2013; Accepted February 12, 2013

The discrete orthogonal wavelet-Galerkin method is illustrated as an effective method for solving partial differential equations (PDE's) with spatially varying parameters on a bounded interval. Daubechies scaling functions provide a concise but adaptable set of basis functions and allow for implementation of varied loading and boundary conditions. These basis functions can also effectively describe C0 continuous parameter spatial dependence on bounded domains. Doing so allows the PDE to be discretized as a set of linear equations composed of known inner products which can be stored for efficient parametric analyses. Solution schemes for both free and forced PDE's are developed; natural frequencies, mode shapes, and frequency response functions for an Euler–Bernoulli beam with piecewise varying thickness are calculated. The wavelet-Galerkin approach is shown to converge to the first four natural frequencies at a rate greater than that of the linear finite element approach; mode shapes and frequency response functions converge similarly.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Gockenbach, M., 2011, Partial Differential Equations: Analytical and Numerical Methods, SIAM USA, Philadelphia.
Fletcher, C., 1984, Computational Galerkin Methods, Springer-Verlag, New York.
Lapidus, L., and Pinder, G., 1982, Numerical Solution of Partial Differential Equations in Science and Engineering, Wiley, New York.
Mallat, S., 2009, A Wavelet Tour of Signal Processing, Academic, New York.
Strang, G., and Nguyen, T., 1996, Wavelets and Filter Banks, Wellesley-Cambridge Press, Wellesley, MA.
Cohen, A., Daubechies, I., and Vial, P., 1993. “Wavelets on the Interval and Fast Wavelet Transforms,” Appl. Comput. Harmon. Anal., 1, pp. 54–81. [CrossRef]
Daubechies, I., 1992, Ten Lecturers on Wavelets, SIAM, Philadelphia.
Williams, J., and Amaratunga, K., 1994, “Introduction to Wavelets in Engineering,” Int. J. Numer. Methods Eng., 37, pp. 2365–2388. [CrossRef]
Amaratunga, K., and Williams, J., 1997, “Wavelet-Galerkin Solution of Boundary Value Problems,” Arch. Comput. Methods Eng., 4(3), pp. 243–285. [CrossRef]
Ala, G., Silvestre, M. L. D., Francomano, E., and Tortorici, A., 2003, “An Advanced Numerical Model in Solving Thin-Wire Integral Equations by Using Semi-Orthogonal Compactly Supported Spline Wavelets,” IEEE Trans. Electron., 45, pp. 218–228.
Joly, P., Maday, Y., and Perrier, V., 1994, “Towards a Method for Solving Partial Differential Equations by Using Wavelet Packet Bases,” Comput. Methods Appl. Mech. Eng., 116, pp. 301–307. [CrossRef]
Beylkin, G., and Keiser, J., 1997, “On the Adaptive Numerical Solution of Nonlinear Partial Differential Equations in Wavelet Bases,” J. Comput. Phys., 132(2), pp. 233–259. [CrossRef]
Chen, M., Hwang, C., and Shih, Y., 1996, “The Computation of Wavelet-Galerkin Approximation on a Bounded Interval,” Int. J. Numer. Methods Eng., 39, pp. 2921–2944. [CrossRef]
Lakestani, M., Rassaghi, M., and Dehghan, M., 2006, “Numerical Solution of the Controlled Duffing Oscillator by Semi-Orthogonal Spline Wavelets,” Phys. Scr., 74, pp. 362–366. [CrossRef]
Lakestani, M., and Dehghan, M., 2006, “The Solution of a Second-Order Nonlinear Differential Equation With Neumann Boundary Conditions Using Semi-Orthogonal B-Spline Wavelets,” Int. J. Comput. Math., 83, pp. 685–694. [CrossRef]
Restrepo, J. M., and Leaf, G. K., 1993, “Inner Product Computations Using Periodized Daubechies Wavelets,” Int. J. Numer. Methods Eng., 40(19), pp. 3557–3578. [CrossRef]
Pernot, S., and Lamarque, C.-H., 2001, “A Wavelet-Galerkin Procedure to Investigate Time-Periodic Systems: Transient Vibration and Stability Analysis,” J. Sound Vib., 245, pp. 845–875. [CrossRef]
Beylkin, G., 1992, “On the Representation of Operators in Bases of Compactly Supported Wavelets,” SIAM (Soc. Ind. Appl. Math.) J. Numer. Anal., 6(6), pp. 1716-1740. [CrossRef]
Romine, C. H., and Peyton, B. W., 1997, “Computing Connection Coefficients of Compactly Supported Wavelets on Bounded Intervals,” U.S. Department of Energy, Report No. ORNL/TM-13413.
Plonka, G., Selig, K., and Tasche, M., 1995, “On the Construction of Wavelets on a Bounded Interval,” Adv. Comput. Math., 4, pp. 357–388. [CrossRef]
Hein, H., and Feklistova, L., 2011, “Free Vibrations of Non-Uniform and Axially Functionally Graded Beams Using Haar Wavelets,” Eng. Sci., 33(12), pp. 3696–3701. [CrossRef]
Boyd, J. P., 2001, Chebyshev and Fourier Spectral Methods, Dover, New York.
Welstead, S., 1999, Fractal and Wavelet Image Compression Techniques, SPIE, Bellingham, WA.
Latto, A., Resnikoff, H., and Tenenbaum, E., 1992, “The Evaluation of Connection Coefficients of Compactly Supported Wavelets,” Proceedings of the French–USA Workshop on Wavelets and Turbulence, Princeton University, Princeton, NJ, June, Springer-Verlag, New York.
Doyle, J., 1997, Wave Propagation in Structures, Springer-Verlag, New York.
Zhang, T., Tian, Y.-C., Tadé, M., and Utomo, J., 2007, “Comments on ‘The Computation of Wavelet-Galerkin Approximation on a Bounded Interval,” Int. J. Numer. Methods Eng., 72, pp. 244–251. [CrossRef]
Humar, J., 1990, Dynamics of Structures, Prentice-Hall, Upper Saddle River, NJ.


Grahic Jump Location
Fig. 1

Thickness variation h(x) of the beam

Grahic Jump Location
Fig. 2

Convergence of the first four natural frequencies using the finite element method (FEM) and wavelet-Galerkin method (W-G)

Grahic Jump Location
Fig. 3

First four mode shapes predicted by the finite element method (FEM) and wavelet-Galerkin method (W-G) using 128 degrees of freedom (J = 7)

Grahic Jump Location
Fig. 4

Frequency response function at midspan (x = 1 m) under a 10Nm harmonic moment acting at x = 2 m



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In