The Logarithmic finite element ("LogFE") method is a novel finite element approach for solving boundary-value problems proposed in [1]. In contrast to the standard Ritz-Galerkin formulation, the shape functions are given on the logarithmic space of the deformation function, which is obtained by the exponentiation of a linear combination of the shape functions given by the degrees of freedom.Unlike many existing multigrid formulations, the LogFE method allows for a very smooth interpolation between nodal values on the coarse grid. It can thus avoid problems with regard to locking and convergence that appear in multigrid applications using only linear interpolation, especially for larger coarsening factors.We illustrate the use of the LogFE method as a coarse grid algorithm, in conjunction with an atomistic finite element method on the fine grid, for calculating the mechanical response of super carbon nanotubes.
MotivationMany multigrid models do not take account of the actual geometry of a structure. For algebraic multigrid methods, the explicit aim is to formulate an algorithm that can be applied independently of the geometry of the structure under consideration. While this significantly simplifies the generation of the coarse meshes, one of its drawbacks is that rather simple prolongation functions, such as linear interpolation, are generally being used to identify the coarse grid correction on the fine grid. This procedure often introduces spurious high frequency deformations (such as simple kinks in piecewise linear functions) that lead to a deterioration of the convergence characteristics on the fine grid. The errors on the fine grid then may lead to deficient estimates of the derivatives with regard to the objective function, which in turn impede on the convergence of the coarse grid or even lead to a breakdown of the optimization algorithm.In this respect, it is important to recall that, generally, a smoothing algorithm, i.e. not a global solver, is being employed on the fine grid. As a result, if the coarse grid experiences locking, that effects cannot be mitigated by the fine grid algorithm, which generally exhibits very slow convergence for low-frequency spatial residuals. As a result, the use of linear interpolation on the coarse grid may lead to locking effects in the global calculation, even for small coarsening factors. While higher order interpolation functions can avoid kinks, they can also introduce spurious high frequencies, e.g. oscillations, on the fine grid.The Logarithmic finite element ("LogFE") method [1-3] provides a smooth interpolation of the field variables on the coarse grid. It thus minimizes spurious high-frequency components on the fine grid and can avoid locking phenomena associated with the use of linear interpolation. As a result, it allows for an increase in the coarsening factor, leading to a significant reduction of the relative number of degrees of freedom on the coarse grid, compared to the fine grid. This reduction of the number of degrees of freedom on the coarse grid is...