Highly accurate finite element method for one-dimensional elliptic interface problems (Q960292)

The authors develop a fourth order accurate finite element based method for the solution of the one-dimensional elliptic interface problem \[ -(\beta u_x)_x+ru= f_c+\nu\delta_\alpha,\quad x\in\Omega=(0,1),\quad u(0)=u(1)=0 \] where \(\beta\) and \(r\) are smooth functions on both \(\Omega^-=(0,\alpha)\) and \(\Omega^+=(\alpha,1)\). The interface jump conditions associated with the discontinuities in the coefficients and singularities of the forces are derived and used to appropriately modify the basis functions such that the jump conditions are satisfied either exactly or approximately. Some numerical experiments are presented.