nodal high order methods on unstructured grids we present a convergent high order accurate scheme for the solution of linear conservation laws in geometrically complex domains as our main example we include a detailed development and analysis of a scheme for the time domain solution of maxwell s equations in a three dimensional domain the fully unstructured spatial discretization is made possible by the use of a high order nodal basis employing multivariate lagrange polynomials defined on the triangles and tetrahedra while the equations themselves are satisfied in a discontinuous galerkin form with the boundary conditions being enforced weakly through a penalty term accuracy stability and convergence of the semidiscrete approximation to maxwell s equations is established rigorously and bounds on the growth of the global divergence error are provided concerns related to efficient implementations are discussed in detail this sets the stage for the presentation of examples verifying the theoretical results and illustrating the versatility flexibility and robustness when solving two and three dimensional benchmark problems in computational electromagnetics pure scattering as well as penetration is discussed and high parallel performance of the scheme is demonstrated

discrete dispersion relation for hp version finite element approximation at high wave number the dispersive properties of high order finite element schemes are analyzed in the setting of the helmholtz equation and an explicit form of the discrete dispersion relation is obtained for elements of arbitrary order it is shown that the numerical dispersion displays three different types of behavior depending on the size of the order of the method relative to the mesh size and the wave number quantitative estimates are obtained for the behavior and rates of decay of the dispersion error in the differing regimes all estimates are fully explicit and are shown to be sharp limits are obtained where transitions between the different regimes occur and used to provide guidelines for the selection of the mesh size and the polynomial order in terms of the wave number so that the dispersion error is controlled

finite element dispersion of cylindrical and spherical acoustic waves a reference solution for cylindrical waves is obtained by wkb approximation finite element representations of both cylindrical and spherical waves satisfy the resolution dependent plane wave dispersion relation and approximate the desired decay rates independent of mesh resolution

the p and hp versions of the finite element method for problems with boundary layers we study the uniform approximation of boundary layer functions exp x d for x 0 1 d 0 1 by the p and hp versions of the finite element method for the p version with fixed mesh we prove super exponential convergence in the range p 1 2 e 2d we also establish for this version an overall convergence rate of o p 1 ln p in the energy norm error which is uniform in d and show that this rate is sharp up to the ln p term when robust estimates uniform in d 0 1 are considered for the p version with variable mesh i e the hp version we show that exponential convergence uniform in d 0 1 is achieved by taking the first element at the boundary layer to be of size o pd numerical experiments for a model elliptic singular perturbation problem show good agreement with our convergence estimates even when few degrees of freedom are used and when d is as small as e g 10 8 they also illustrate the superiority of the hp approach over other methods including a low order h version with optimal exponential mesh refinement the estimates established in this paper are also applicable in the context of corresponding spectral element methods

modeling of electromagnetic absorption scattering problems using hp adaptive finite elements a model problem for the steady state form of maxwell s equations is considered the problem is formulated in a weak form using a lagrange multiplier laying down a foundation for a general class of novel hp adaptive fe approximations a convergence proof for affine elements is presented the proposed method is illustrated and verified with a series of 2d experiments including elements with curved boundaries and nonhomogeneous media

a comparison of high order polynomial and wave based methods for helmholtz problems the application of computational modelling to wave propagation problems is hindered by the dispersion error introduced by the discretisation two common strategies to address this issue are to use high order polynomial shape functions e g hp fem or to use physics based or trefftz methods where the shape functions are local solutions of the problem typically plane waves both strategies have been actively developed over the past decades and both have demonstrated their benefits compared to conventional finite element methods but they have yet to be compared in this paper a high order polynomial method p fem with lobatto polynomials and the wave based discontinuous galerkin method are compared for two dimensional helmholtz problems a number of different benchmark problems are used to perform a detailed and systematic assessment of the relative merits of these two methods in terms of interpolation properties performance and conditioning it is generally assumed that a wave based method naturally provides better accuracy compared to polynomial methods since the plane waves or bessel functions used in these methods are exact solutions of the helmholtz equation results indicate that this expectation does not necessarily translate into a clear benefit and that the differences in performance accuracy and conditioning are more nuanced than generally assumed the high order polynomial method can in fact deliver comparable and in some cases superior performance compared to the wave based dgm in addition to benchmarking the intrinsic computational performance of these methods a number of practical issues associated with realistic applications are also discussed

dispersion minimized mass for isogeometric analysis abstract we introduce the dispersion minimized mass for isogeometric analysis to approximate the structural vibration which we model as a second order differential eigenvalue problem the dispersion minimized mass reduces the eigenvalue error significantly from the optimum order of 2 p to the superconvergence order of 2 p 2 for the p th order isogeometric elements with maximum continuity which in return leads to a more accurate method we first establish the dispersion error where the leading error term is explicitly written in terms of the stiffness and mass entries for arbitrary polynomial order isogeometric elements we derive the dispersion minimized mass in one dimension by solving a p dimensional local matrix problem for the p th order approximation and then extend it to multiple dimensions on tensor product grids we show that the dispersion minimized mass can also be obtained by approximating the mass matrix using optimally blended quadratures we generalize the lower order quadrature blending results to arbitrary polynomial order isogeometric approximations as well as to arbitrary quadrature rules various numerical examples validate the eigenvalue and eigenfunction error estimates

high order discontinuous galerkin methods for the elastodynamics equation on polygonal and polyhedral meshes abstract we propose and analyze a high order discontinuous galerkin finite element method for the approximate solution of wave propagation problems modeled by the elastodynamics equations on computational meshes made by polygonal and polyhedral elements we analyze the well posedness of the resulting formulation prove h p version error a priori estimates and present a dispersion analysis showing that polygonal meshes behave as classical simplicial quadrilateral grids in terms of dispersion properties the theoretical estimates are confirmed through various two dimensional numerical verifications

finite element solution of the helmholtz equation with high wave number part ii the h p version of the fem in this paper which is part ii in a series of two the investigation of the galerkin finite element solution to the helmholtz equation is continued while part i contained results on the h version with piecewise linear approximation the present part deals with approximation spaces of order p ge 1 as in part i the results are presented on a one dimensional model problem with dirichlet robin boundary conditions in particular there are proven stability estimates both with respect to data of higher regularity and data that is bounded in lower norms the estimates are shown both for the continuous and the discrete spaces under consideration further there is proven a result on the phase difference between the exact and the galerkin finite element solutions for arbitrary p that had been previously conjectured from numerical experiments these results and further preparatory statements are then employed to show error estimates for the galerkin finite element method fem it becomes evident that the error estimate for higher approximation can with certain assumptions on the data be written in the same form as the piecewise linear case namely as the sum of the error of best approximation plus a pollution term that is of the order of the phase difference the paper is concluded with a numerical evaluation

an improved subspace selection algorithm for meshless collocation methods choosing data points is a common problem for researchers who employ various meshless methods for solving partial differential equations on the one hand high accuracy is always desired on the other ill conditioning problems of the resultant matrices which may lead to unstable algorithms prevent some researchers from using meshless methods for example the optimal placements of source points in the method of fundamental solutions or of the centers in the radial basis functions method are always unclear intuitively such optimal locations will depend on many factors the partial differential equations the domain the trial basis used i e the employed method itself 1pt the computational precisions some user defined parameters and so on such complexity makes the hope of having an optimal centers placement unpromising in this paper we provide a data dependent algorithm that adaptively selects centers based on all the other variables copyright 2009 john wiley sons ltd

estimation of the dispersion error in the numerical wave number of standard and stabilized finite element approximations of the helmholtz equation an estimator for the error in the wave number is presented in the context of finite element approximations of the helmholtz equation the proposed estimate is an extension of the ideas introduced in steffens and d iez comput methods appl mech engng 2009 198 1389 1400 in the previous work the error assessment technique was developed for standard galerkin approximations here the methodology is extended to deal also with stabilized approximations of the helmholtz equation thus the accuracy of the stabilized solutions is analyzed including also their sensitivity to the stabilization parameters depending on the mesh topology the procedure builds up an inexpensive approximation of the exact solution using post processing techniques standard in error estimation analysis from which the estimate of the error in the wave number is computed using a simple closed expression the recovery technique used in steffens and diez comput methods appl mech engng 2009 198 1389 1400 is based in a polynomial least squares fitting here a new recovery strategy is introduced using exponential in a complex setup trigonometric local approximations respecting the nature of the solution of the wave problem copyright 2011 john wiley sons ltd

error analysis for the numerical evaluation of the diagonal forms of the scalar spherical addition theorem the numerical solution of wave scattering from large objects or from a large cluster of scatterers requires excessive computational resources and it becomes necessary to use approximate but fast methods such as the fast multipole method however since these methods are only approximate it is important to have an estimate for the error introduced in such calculations an analysis of the error for the fast multipole method is presented and estimates for truncation and numerical integration errors are obtained the error caused by polynomial interpolation in a multilevel fast multipole algorithm is also analyzed the total error introduced in a multilevel implementation is also investigated numerically

a discrete commutator theory for the consistency and phase error analysis of semi discrete c 0 finite element approximations to the linear transport equation a novel method based on a discrete commutator for the analysis of consistency error and phase relations for semi discrete continuous finite element approximation of the one way wave equation is presented the technique generalizes to arbitrary dimension accommodates the use of compatible quadratures does not require the use of complex calculations is applicable on non uniform mesh geometries and is especially useful when conventional taylor series or fourier approaches are intractable following the theory the analysis method is demonstrated for several test cases

dispersion and pollution of the fem solution for the helmholtz equation in one two and three dimensions for high wave numbers the helmholtz equation suffers the so called pollution effect this effect is directly related to the dispersion a method to measure the dispersion on any numerical method related to the classical galerkin fem is presented this method does not require to compute the numerical solution of the problem and is extremely fast numerical results on the classical galerkin fem p method is compared to modified methods presented in the literature a study of the influence of the topology triangles is also carried out the efficiency of the different methods is compared the numerical results in two of the mesh and for square elements show that the high order elements control the dispersion well the most effective modified method is the qsfem 1 2 but it is also very complicated in the general setting the residual free bubble 3 4 is effective in one dimension but not in higher dimensions the least square method 1 5 approach lowers the dispersion but relatively little the results for triangular meshes show that the best topology is the criss cross pattern copyright 1999 john wiley sons ltd

integration of hp adaptivity and a two grid solver for electromagnetic problems we present implementation details and analyze convergence of a two grid solver forming the core of a fully automatic hp adaptive strategy for electromagnetic problems the solver delivers a solution for a fine grid obtained from an arbitrary coarse hp grid by a global hp refinement the classical v cycle algorithm combines an overlapping block jacobi smoother with optimal relaxation and a direct solve on the coarse grid a theoretical analysis of the two grid solver is illustrated with numerical experiments several electromagnetic applications show the efficiency of combining the fully automatic hp adaptive strategy with the two grid solver

the influence of the mass matrix on the dispersive nature of the semi discrete second order wave equation abstract the application of discrete solution methods to the second order wave equation can yield a dispersive representation of the non dispersive wave propagation problem resulting in a phase speed that depends not only upon the wavelength of the signal being propagated but also upon the direction of propagation in this work the dependence of the dispersive errors on the wave propagation direction mesh aspect ratio and wave number is investigated with the goal of understanding and hopefully reducing the phase and group errors associated with the two dimensional bilinear finite element an analysis of the dispersive effects associated with the consistent row sum lumped and higher order mass matrices has led to a reduced coupling penta diagonal mass matrix that yields improved phase and group errors with respect to wavelength and propagation direction the influence of row sum lumping the finite element mass matrix is demonstrated to always introduce lagging phase and group error while a linear combination of the lumped and consistent mass matrices i e a higher order mass matrix is shown to improve the dispersion characteristics of both the reduced and full integration element using a 2 5 phase error metric the higher order mass matrix can reduce the mesh resolution requirements relative to the lumped mass matrix by nearly a factor of two in one dimension a factor of 4 in two dimensions similarly the reduced coupling mass matrix can reduce the grid requirements by a factor of approximately four relative to the lumped mass matrix a factor of 16 in two dimensions although the reduced coupling mass matrix requires additional storage relative to the lumped mass the factor of 16 reduction in grid requirements promises tremendous overall savings for computing a fixed wavelength spectrum or alternatively for resolving shorter wavelengths with a fixed grid resolution

goal oriented adaptivity using unconventional error representations for the multidimensional helmholtz equation in goal oriented adaptivity goa the error in the quantity of interest qoi is represented using the error functions of the direct and adjoint problems this error representation is subsequently bounded above by element wise error indicators that are used to drive optimal refinements in this work we propose to replace in the error representation the adjoint problem by an alternative operator the main advantage of the proposed approach is that when judiciously selecting such alternative operator the corresponding upper bound of the error representation becomes sharper leading to a more efficient goa while the method can be applied to a variety of problems we focus here on two and threedimensional 2d and 3d helmholtz problems we show via extensive numerical experimentation that the upper bounds provided by the alternative error representations are sharper than the classical ones and lead to a more robust p adaptive process we also provide guidelines for finding operators delivering sharp error representation upper bounds we further extend the results to a convectiondominated diffusion problem as well as to problems with discontinuous material coefficients finally we consider a sonic logging while drilling lwd problem to illustrate the applicability of the proposed method

robustness and dispersion analysis of the partition of unity finite element method applied to the helmholtz equation abstract the partition of unity finite element method pufem has been widely used for the numerical simulation of the helmholtz equation in different physical settings in fact it is a numerical pollution free alternative method to the classical piecewise polynomial based finite element methods taking into account a plane wave enrichment of the piecewise linear finite element method the main goal of this work is focused on the derivation of the numerical dispersion relation and the robustness analysis of the pufem discretization when a spurious perturbation is presented in the wave number value used in the enrichment definition from the one dimensional helmholtz equation the discrete wave number is estimated based on a bloch s wave analysis and a priori error estimates are computed explicitly in terms of the mesh size the wave number and the perturbation value