Implementation Strategy for the Numerical Efficiency Improvement of the Multiscale Interpolation Wavelet-Galerkin Method

  • Seo Jeong Hun (School of Mechanical and Aerospace Engineering and National Creative Research Initiatives Center for Multiscale Design, Seoul National University) ;
  • Earmme Taemin (School of Mechanical and Aerospace Engineering and National Creative Research Initiatives Center for Multiscale Design, Seoul National University) ;
  • Jang Gang-Won (School of Mechanical Engineering, Kunsan National University) ;
  • Kim Yoon Young (School of Mechanical and Aerospace Engineering and National Creative Research Initiatives Center for Multiscale Design, Seoul National University)
  • 발행 : 2006.01.01

초록

The multi scale wavelet-Galerkin method implemented in an adaptive manner has an advantage of obtaining accurate solutions with a substantially reduced number of interpolation points. The method is becoming popular, but its numerical efficiency still needs improvement. The objectives of this investigation are to present a new numerical scheme to improve the performance of the multi scale adaptive wavelet-Galerkin method and to give detailed implementation procedure. Specifically, the subdomain technique suitable for multiscale methods is developed and implemented. When the standard wavelet-Galerkin method is implemented without domain subdivision, the interaction between very long scale wavelets and very short scale wavelets leads to a poorly-sparse system matrix, which considerably worsens numerical efficiency for large-sized problems. The performance of the developed strategy is checked in terms of numerical costs such as the CPU time and memory size. Since the detailed implementation procedure including preprocessing and stiffness matrix construction is given, researchers having experiences in standard finite element implementation may be able to extend the multi scale method further or utilize some features of the multiscale method in their own applications.

키워드

참고문헌

  1. Basic Linear Algebra Subprograms Technical (BLAST) Forum, 1997, DRAFT Document for the Basic Linear Algebra Subprograms (BLAS) Standard
  2. Bertoluzza, S., 1997, 'An Adaptive Collocation Method Based on Interpolating Wavelets,' Multiscale Wavelet Methods for PDEs, Dahmen W, Kurdila AJ, Oswald P (eds): Academic Press, San Diego, pp. 109-135
  3. Carney, S., Heroux, M. A., Li, G., Pozo, R., Remington, K. A. and Wu, K., 1996, A Revised Proposal for Sparse BLAS Toolkit, SPARKER Working Note #3
  4. Christon, M. A. and Roach, D. W., 2000, 'The Numerical Performance of Wavelets for PDEs: the Multi-scale Finite Element,' Computational Mechanics, Vol. 25, pp. 230-244 https://doi.org/10.1007/s004660050472
  5. Cohen, A. and Masson. R., 1999, 'Wavelet Methods for Second Order Elliptic Problems-Preconditioning and Adaptivity,' SIAM J. Sci. Comp., Vol. 21, pp. 1006-1026 https://doi.org/10.1137/S1064827597330613
  6. Cohen, A., Dahmen, W. and DeVore, R., 1998, Adaptive Wavelet Methods for Elliptic Operator Equations-Convergence Rates, RWTH Aachen, IGPM Preprint No. 165
  7. Dahmen, W., 2001, 'Wavelet Methods for PDEs-Some Recent Developments,' Journal of Computational and Applied Mathematics, Vol. 128, pp. 133-185 https://doi.org/10.1016/S0377-0427(00)00511-2
  8. Diaz, A.R., 1999, 'A Wavelet-Galerkin Scheme for Analysis of Large-Scale Problems on Simple Domains,' Int. J. Numer. Meth. Engng., Vol. 44, pp.1599-1616 https://doi.org/10.1002/(SICI)1097-0207(19990420)44:11<1599::AID-NME556>3.0.CO;2-P
  9. Glowinski, R., Pan, T. W., Wells, R. O. Jr. and Zhou X., 1994, Wavelet and Finite Element Solutions for the Neumann Problem Using Fictitious Domains, Technical Report, Computational Mathematics Lab, Rice University
  10. Jang, G. W., 2004, Shape and Topology Optimization Using the Multiscale Method, Ph.D. Thesis, Seoul National University
  11. Jang, G. W., Kim, J. E. and Kim, Y. Y., 2004a, 'Multiscale Galerkin Method Using Interpolation Wavelets for Two-Dimensional Elliptic Problems in General Domains,' Int. J. Numer. Meth. Engng., Vol. 59, pp. 225-253 https://doi.org/10.1002/nme.872
  12. Jang, G. W., Kim, Y. Y. and Choi, K. K., 2004b, 'Remesh-Free Shape Optimization Using the Wavelet-Galerkin Method,' International Journal of Solids and Structures, Vol. 41, pp. 6465-6483 https://doi.org/10.1016/j.ijsolstr.2004.05.010
  13. Latto, A., Resnikoff, H. K. and Tenenbaum, E., 1991, 'The Evaluation of Connection Coefficients of Compactly Supported Wavelets,' In Proceedings of the USA-French Workshop on Wavelets and Turbulance, Princeton University
  14. Mallat, S., 1998, A Wavelet Tour of Signal Processing. Academic Press: London
  15. Saad, Y., 1994, SPARSKIT: A Basic Tool-Kit for Sparse Matrix Computations. VERSION2