An Iteratively Adaptive Multiscale Finite Element Method with Application to Interface Problems

以作者查詢圖書館館藏

、以作者查詢臺灣博碩士

、以作者查詢全國書目

、勘誤回報

、線上人數：41

、訪客IP：3.21.248.47

姓名

姚建州(Chien-Chou Yao) 查詢紙本館藏

畢業系所

數學系

論文名稱

(An Iteratively Adaptive Multiscale Finite Element Method with Application to Interface Problems)

相關論文

★ 非線性塊狀高斯消去牛頓演算法在噴嘴流體的應用	★ 以平行 Newton-Krylov-Schwarz 演算法解 Poisson-Boltzmann 方程式的有限元素解在膠體科學上的應用
★ 最小平方有限元素法求解對流擴散方程以及使用Bubble函數的改良	★ Bifurcation Analysis of Incompressible Sudden Expansion Flows Using Parallel Computing
★ Parallel Jacobi-Davidson Algorithms and Software Developments for Polynomial Eigenvalue Problems in Quantum Dot Simulation	★ An Inexact Newton Method for Drift-DiffusionModel in Semiconductor Device Simulations
★ Numerical Simulation of Three-dimensional Blood Flows in Arteries Using Domain Decomposition Based Scientific Software Packages in Parallel Computers	★ A Parallel Fully Coupled Implicit Domain Decomposition Method for the Stabilized Finite Element Solution of Three-dimensional Unsteady Incompressible Navier-Stokes Equations
★ A Study for Linear Stability Analysis of Incompressible Flows on Parallel Computers	★ Parallel Computation of Acoustic Eigenvalue Problems Using a Polynomial Jacobi-Davidson Method
★ Numerical Study of Algebraic Multigrid Methods for Solving Linear/Nonlinear Elliptic Problems on Sequential and Parallel Computers	★ A Parallel Multilevel Semi-implicit Scheme of Fluid Modeling for Numerical Low-Temperature Plasma Simulation
★ Performance Comparison of Two PETSc-based Eigensolvers for Quadratic PDE Problems	★ A Parallel Two-level Polynomial Jacobi-Davidson Algorithm for Large Sparse Dissipative Acoustic Eigenvalue Problems
★ A Full Space Lagrange-Newton-Krylov Algorithm for Minimum Time Trajectory Optimization	★ Parallel Two-level Patient-specific Numerical Simulation of Three-dimensional Rheological Blood Flows in Branching Arteries

檔案

[Endnote RIS 格式]

[Bibtex 格式]

[相關文章]

[文章引用]

[完整記錄]

[館藏目錄]

[檢視]

[下載]

本電子論文使用權限為同意立即開放。
已達開放權限電子全文僅授權使用者為學術研究之目的，進行個人非營利性質之檢索、閱讀、列印。
請遵守中華民國著作權法之相關規定，切勿任意重製、散佈、改作、轉貼、播送，以免觸法。

摘要(中)

本論文是解橢圓介面問題(elliptic interface problem). 為了有更準確的結果，我們用多尺度有限元素法(multiscale finite element method）來求解。多尺度有限元素法主要由兩個部分組成，分別是局部(local) 多尺度基底函數和這些基底函數所組成的總體粗網格(global coarse formulation). 由於多尺度基底函數的準確與否會直接影響數值結果，因此設定多尺度基底函數邊界會是重要環節。本研究旨是為多尺度有限元素法提出一個好的離散方式，求解橢圓介面問題的數值PDE. 此方式透過迭代過程對局部粗網格邊界值更新得以讓粗網格(coarse-grid) 與細網隔(fine-grid)間資訊做交換，得到原始PDE的解做為設定下一次迭代局部粗網格所須邊界條件的資訊，避免花費昂貴的計算成本直接解細網格的資訊。另外，迭代法的操作對於數值誤差有消除的效果，因此每次計算得到的細網格數值不須要求高精確度的計算結果，透過迭代的過程得以修正誤差。
然而，為了得到更準確的細網格數值解做為局部粗網格邊界條件的設定，在每次迭代中所得到細網格的數值解會透過數值迭代法做幾步smoothing的迭代運算以削去局部粗網格邊界點上震盪的高頻誤差(high frequency error)，如此也有利加速得到數值結果的收斂。針對被介面曲線通過的局部細網格而言，為了有更吻合的數值估算結果，將採取immersed finite element基底的方式做計算。
在我們數值結果中，i-ApMsFEM 這方式能夠有效消除介面上和粗網格邊界上的誤差。此外L2-norm可達二次收斂和H1-norm為一次收斂的結果; 誤差值也能和係數比例(contrast ratios)呈獨立關係，這些結果都符合理論敘述。

摘要(英)

This thesis is to develop an accurate and effective numerical scheme for solving elliptic interface problem based on a framework of the multiscale finite element method. The two major ingredients of multiscale finite element method are the construction of multiscale basis
functions on coarse elements and the global coarse formulation using these basis functions. The multiscale basis functions are defined to satisfy the original governing partial differential equation restricted in a coarse element with some proper boundary condition. The choice of the boundary condition for multiscale basis functions is important since it affects the accuracy of the global solution significantly. In this work, we propose an adaptive multiscale finite element method through an iterative process to improve the boundary settings for the multiscale basis functions, namely i-ApMsFEM. To demonstrate the capability of the proposed i-ApMsFEM to a variety of high contrast elliptic interface problems, we consider a series of testing problems. An immersed finite element method is used for constructing numerically multiscale basis and for formulating the coarse-grid problem. Our numerical
results show that i-ApMsFEM effectively eliminate the errors along the interface curve and the boundaries of coarse elements to produce highly accurate solutions. In addition, the second-order rate convergence in the L2-norm (or first-order in the H1-norm) is achieved which is independent of the contrast ratios.

關鍵字(中)

★ 多尺度有限元素法
★ 橢圓介面問題
★ 自適應邊界條件

關鍵字(英)

★ Multiscale finite elements
★ elliptic interface problem
★ adaptive boundary condition

論文目次

1 Introduction . . . . . . . . . . . . . . . . . . . . . 1
2 1D elliptic interface problem . . . . . . . . . . . . . 3
2.1 1D model problem . . . . . . . . . . . . . . . . . . 3
2.2 Weak formulation for the elliptic interface problem 3
2.3 Galerkin finite element formulation . . . . . . . . 9
2.4 Immersed finite element formulation . . . . . . . . 11
2.5 Multiscale finite element formulation . . . . . . . 14
2.6 Numerical integration . . . . . . . . . . . . . . . 19
2.7 Numerical results for 1D elliptic interface problem 22
3 2D elliptic interface problem 34
3.1 2D model problem . . . . . . . . . . . . . . . . . 34
3.2 Weak formulation for the elliptic interface problem 35
3.3 Galerkin finite element formulation . . . . . . . . 37
3.4 Nonconforming immersed finite element formulation . 38
3.5 Numerical integration . . . . . . . . . . . . . . . 40
4 Iteratively adaptive multiscale finite element method ..45
5 Numerical experiments and discussion. . . . . . . . . . 52
6 Conclusion. . . . . . . . . . . . . . . . . . . . . . . 82
Bibliography. . . . . . . . . . . . . . . . . . . . . . . 83
A Multiscale finite element formulation for Explicit formulation . . . . . . . . . . . . . . . . . . . . . . . 84
B The coefficients of immersed finite element basis for one-dimensional model . . . . . . . . . . . . . . . . . . . . 87
C The error of the solution evaluated for one-dimensional model . . . . . . . . . . . . . . . . . . . . . . . . . . 90
D The error of single inclusion experiment model . . . . .92
E Algorithm . . . . . . . . . . . . . . . . . . . . . . . 94

參考文獻

[1] Eric B Becker, Graham F Carey, and John Tinsley Oden. Finite elements, an introduction: Volume i. ., 258, page 1981, 1981.
[2] Wenqiang Feng. Immersed finite element method for interface problems with algebraic multigrid solver. 2013.
[3] Xiaoming He, Tao Lin, and Yanping Lin. Immersed finite element methods for elliptic interface problems with non-homogeneous jump conditions. Int.J. Numer. Anal. Model, 8(2):284–301, 2011.
[4] Thomas JR Hughes. The finite element method: linear static and dynamic finite element analysis. Courier Dover Publications, 2012.
[5] Zhilin Li and Kazufumi Ito. The immersed interface method: numerical solutions of PDEs involving interfaces and irregular domains, volume 33.
Siam, 2006.
[6] Zhilin Li, Tao Lin, and Xiaohui Wu. New cartesian grid methods for interface
problems using the finite element formulation. Numerische Mathematik, 96(1):61–98, 2003.
[7] Yousef Saad. Iterative methods for sparse linear systems. Siam, 2003.
[8] Ulrich Trottenberg, Cornelius W Oosterlee, and Anton Schuller. Multigrid. Academic press, 2000.
[9] Xu Zhang. Nonconforming Immersed Finite Element Methods for Interface Problems. PhD thesis, Virginia Polytechnic Institute and State University, 2013.

指導教授

黃楓南(Feng-Nan Hwang)

審核日期

2014-8-27

推文