adaptive multigrid fe methods an optimal way to solve pdes l.
Download
Skip this Video
Loading SlideShow in 5 Seconds..
Adaptive Multigrid FE Methods -- An optimal way to solve PDEs PowerPoint Presentation
Download Presentation
Adaptive Multigrid FE Methods -- An optimal way to solve PDEs

Loading in 2 Seconds...

play fullscreen
1 / 34
hart

Adaptive Multigrid FE Methods -- An optimal way to solve PDEs - PowerPoint PPT Presentation

154 Views
Download Presentation
Adaptive Multigrid FE Methods -- An optimal way to solve PDEs
An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

  1. Adaptive Multigrid FE Methods -- An optimal way to solve PDEs Zhiming Chen Institute of Computational Mathematics Chinese Academy of Sciences Beijing 100080

  2. Adaptive Concept Refine: more nodes around singularities Coarsen: less nodes in smooth region

  3. The adaptive method finds the solution of given tolerance on a self-generated mesh according to the properties of the solution (singularities,oscillations). • An “optimal mesh” is the mesh on which the error is approximately the same on each element. This motivates the error equi-distribution strategy. • The adaptive FEM based on a posteriori error estimates provides a systematic way to refine or coarsen the mesh according to the local a posteriori error estimators on each element.

  4. A Linear Elliptic Problem Elliptic problem with piecewise constant coefficients: Variational problem:

  5. Discrete Problem: find such that A priori error estimate:

  6. A posteriori error estimate (Babuska & Miller, 1987) The error indicator where

  7. Theorem (Verfürth, 1992): We have where

  8. Adaptive Algorithm Solve → Estimate → Refine/Coarsen Error equi-distribution strategy where tolerance, constant , number of elements in

  9. Numerical Experiments where and Exact solution (Kellogg)

  10. The exact solution for

  11. FEM with uniform mesh 128x128 mesh: 512x512 mesh: 1024x1024 mesh: Convergence rate: A priori error analysis implies that one must introduce nodes in each space direction to bring the energy error under 0.1.

  12. The surface plot of the relative error The maximum of the relative error is 0.2368.

  13. The adaptive mesh of 2673 nodes. The energy error is 0.07451.

  14. The surface plot of the adaptive solution and the relative error . The maximum of the relative error is0.0188.

  15. Definition Let be the sequence of FE solutions generated by the adaptive algorithm. The meshes and the associated numerical complexity are called quasi-optimal if are valid asymptotically. DOFs(k) is the number of degree of freedoms of the mesh .

  16. Quasi-optimality of the estimators. The quasi-optimal decay is indicated by the dotted line of slope –1/2.

  17. Gauss-Seidel Iteration Method

  18. 5-point finite difference scheme

  19. Multigrid V-cycle Algorithm

  20. Adaptive Multigrid Method ◆ Local relaxation: Gauss-Seidel relaxation performed only on new nodes and their immediate neighboring nodes ◆ Each multigrid iteration requires only O(N) number of operations ◆ Theorem (Wu and Chen): We have ◆ Numerical Example

  21. Applications • Continuous casting problem Chen, Nochetto and Schmidt (2000) • Wave scattering by periodic structure Chen and Haijun Wu (2002) • Convection diffusion problem Chen and Guanghua Ji (2003)

  22. Continuous Casting problem

  23. 例子:振荡铸钢速度系数 变化铸钢速度: 速度 v(t), 单元个数, 时间步长

  24. t=0.05 和 t=0.07 时的网格和温度。

  25. 周期结构上电磁波的散射问题

  26. A Linear Convection Diffusion Problem Rotating Cylinder problem:

  27. Convergence Rate Epsilon=10e-3 Epsilon=10e-5

  28. Thank you !