Finite Element Method and Natural Boundary Reduction
Finite Element Method and Natural Boundary Reduction
复制标题
DOI:
--
复制
发表时间:
2010
期刊:
影响因子:
--
通讯作者:
Feng KASOc;Feng Eang
中科院分区:
文献类型:
--
作者:
Feng KASOc;Feng Eang
One of the major advances in numerical methods for partial differential equations made in the recent twenty years is the finite element method (FEM). The method is based on the variational formulation of elliptic equations and on the triangulated approximations. The first component, the variational principle, is an old one ajid leads to the classical BayleighBitz method, which, though successful in the past, suffers from numerical instability and geometric inflexibility, originating from the analytic approximations adopted, but unnoticed in the pre-computer times due to the limited size and complexity of the problems then attacked. The second component, the triangulated local approximations, used but not exploited in full in the finite difference methods, is more elementary and much older. Dating back to ancient times, it was for a long time overshadowed by the later achievements in analytic approximations, but revived eventually due to its innate stability and flexibility, which becomes important in the computer era. A judicious combination of the two old components, conventionally in juxtaposition, gives rise to the FEM, an innovation of general applicability, especially suited for problems of great complexity as well as for computer usage. In FEM, all the essential properties of elliptic operators, e.g., symmetry, coerciveness and locality are well preserved after discretization. This leads, on the one hand, to an efficient computational scheme and, on the other hand, to a sound theoretical foundation, on which the Sobolev space theory of elliptic equations is invoked in a natural way, ensuring the reliability of the method in practice. Moreover, the logic of FEM is simple, intuitive and easy to be implemented on the computer,