Abstract:Ever since the pioneering work of Ihlenburg and Babuska for the one-dimensional Helmholtz equation (1995, Comput. Math. Appl., 30, 9-37), a lot of pre-asymptotic analyses, including dispersion and pollution analyses, have been carried out for finite element methods (FEMs) for the Helmholtz equation, but no rigorous pre-asymptotic error analysis of the FEM for the Helmholtz equation in two and three dimensions has been done yet. This paper addresses pre-asymptotic stability and error estimates of some continuous interior penalty finite element methods (CIP-FEM) and the FEM using piecewise linear polynomials for the Helmholtz equation in two and three dimensions. The CIP-FEM, which was first proposed by Douglas and Dupont for elliptic and parabolic problems in the 1970s and then successfully applied to convection-dominated problems as a stabilization technique, modifies the bilinear form of the FEM by adding a least squares term penalizing the jump of the gradient of the discrete solution at mesh interfaces, while the penalty parameters are chosen as complex numbers here instead of real numbers as usual. In this paper, it is proved that if the penalty parameter is a pure imaginary number i gamma with 0 < gamma <= C, then the CIP-FEM is stable (hence well posed) without any mesh constraint. Moreover, the error of the CIP-FEM in the H-1 norm is bounded by C(1)kh + C(2)k(3)h(2) when k(3)h(2) <= C-0 and by C(1)kh + C-2/gamma when k(3)h(2)> C-0 and <= C-3, where k is the wave number, h is the mesh size, and the C's are some positive constants independent of k, h and gamma. By taking gamma -> 0+ in the CIP-FEM it is shown that if k(3)h(2) <= C-0, then the FEM attains a unique solution and satisfies the error estimate C(1)kh + C(2)k(3)h(2) in the H-1 norm. Previous analytical results for the classical linear FEM in higher dimensions require the assumption that k(2)h is small, but it is too strict for a medium or high wave number. Optimal order L-2-error estimates are also derived. Numerical results are provided to verify the theoretical findings and illustrate the performance of the two methods. In particular, it is shown that the pollution effect may be greatly reduced by tuning the penalty parameters in the CIP-FEM. The theoretical and numerical results reveal that the CIP-FEM has advantages over the FEM and is worthy of being investigated further. The pre-asymptotic error analysis of the hp version of the CIP-FEM and the FEM is studied in Part II.

An optimized CIP-FEM to reduce the pollution errors for the Helmholtz equation on a general unstructured mesh

Pre-asymptotic Error Analysis of CIP-FEM and FEM for the Helmholtz Equation with High Wave Number. Part I: Linear Version

Dispersion Analysis of CIP-FEM for the Helmholtz Equation

Continuous Interior Penalty Finite Element Method for Helmholtz Equation with High Wave Number: One Dimensional Analysis

A CIP-based Numerical Simulation of Wave Interaction with a Fluid-Filled Membrane Submerged Breakwater

Numerical Simulation of Water Entry of Two-Dimensional Structures with Complex Geometry Using a CIP-based Model

Simulations of Interaction Between Flow and a Thin Elastic Beam Using Cip-Based Model with Adaptive Mesh

Numerical Simulation of Wave Interaction with Coastal Structures Using a CIP-based Method

A Meshless Method for Some Inverse Problems Associated with the Helmholtz Equation

Higher-order FEM and CIP-FEM for Helmholtz equation with high wave number and perfectly matched layer truncation

An Improved CIP-based Numerical Model for Simulating Free-Surface Flow with Adaptive Mesh

Pre-asymptotic Error Analysis of CIP-FEM and FEM for Helmholtz Equation with High Wave Number. Part II: $hp$ Version

Adaptive FEM for Helmholtz Equation with Large Wavenumber

Efficient Implementation of the CPML in 3-D Hybrid Implicit-Explicit FDTD Method

A Conservative Finite Volume-Particle Hybrid Method for Simulation of Incompressible Interfacial Flow

Adaptive Finite Element Method for a Nonlinear Helmholtz Equation with High Wave Number

Single Field HIE FDTD Solver with Nonuniform Grid Stratagem for More Fine Structure

An Interface Penalty Finite Element Method for Elliptic Interface Problems on Piecewise Meshes

A three-dimensional parallel CIP-based model for wave-structure interaction

A cell‐based smoothed radial point interpolation method with virtual nodes for three‐dimensional mid‐frequency acoustic problems

Continuous Interior Penalty Finite Element Methods for the Time-Harmonic Maxwell Equation with High Wave Number