scholarly journals Analysis of Galerkin and streamline-diffusion FEMs on piecewise equidistant meshes for turning point problems exhibiting an interior layer

2018 ◽  
Vol 123 ◽  
pp. 121-136
Author(s):  
Simon Becher
2012 ◽  
Vol 12 (2) ◽  
pp. 206-220 ◽  
Author(s):  
Eugene O'Riordan ◽  
Jason Quinn

Abstract A finite difference scheme on special piecewise-uniform grids condensing in tA linear singularly perturbed interior turning point problem with a continuous convection coefficient is examined in this paper. Parameter uniform numerical methods composed of monotone finite difference operators and piecewise-uniform Shishkin meshes, are constructed and analysed for this class of problems. A refined Shishkin mesh is placed around the location of the interior layer and we consider disrupting the centre point of this fine mesh away from the point where the convection coefficient is zero. Numerical results are presented to illustrate the theoretical parameter-uniform error bounds established.


Author(s):  
Simon N. Gaulter ◽  
Nicholas R. T. Biggs

In this paper, we develop an asymptotic scheme to approximate the trapped mode solutions to the time harmonic wave equation in a three-dimensional waveguide with a smooth but otherwise arbitrarily shaped cross section and a single, slowly varying ‘bulge’, symmetric in the longitudinal direction. Extending previous research carried out in the two-dimensional case, we first use a WKBJ-type ansatz to identify the possible quasi-mode solutions that propagate only in the thicker region, and hence find a finite cut-on region of oscillatory behaviour and asymptotic decay elsewhere. The WKBJ expansions are used to identify a turning point between the cut-on and cut-off regions. We note that the expansions are non-uniform in an interior layer centred on this point, and we use the method of matched asymptotic expansions to connect the cut-on and cut-off regions within this layer. The behaviour of the expansions within the interior layer then motivates the construction of a uniformly valid asymptotic expansion. Finally, we use this expansion and the symmetry of the waveguide around the longitudinal centre, x =0, to extract trapped mode wavenumbers, which are compared with those found using a numerical scheme and seen to be extremely accurate, even to relatively large values of the small parameter.


Sign in / Sign up

Export Citation Format

Share Document