about the use of standard integration schemes for x fem in solid mechanics plasticity abstract the introduction of plasticity in an extended finite element x fem framework is very recent and arises several technical and scientific difficulties due to the moving position of integration points as the crack proceeds if a fixed grid is not used currently the existing analyses in plasticity with x fem are limited to monotonic plasticity approximated by non linear elasticity using the hrr field decomposition however such an analysis cannot be applied in fatigue where cyclic plasticity plays a major role in crack propagation the aim of this study is to tackle one of the difficulties induced by cyclic plastic computation the projection of the internal variables during crack propagation the discontinuous heaviside function in the elements cut by the crack is replaced by a new continuous enrichment function this new function is polynomial and its degrees and coefficients depend on the type of element and on the position of the interface in the element the automatic determination of the new enrichment and its introduction in a linear and non linear framework have been implemented for linear elements in two dimensions and three dimensions and tested with convergence tests in elasticity and plasticity

almost invariant subspaces an approach to high gain feedback design part ii almost conditionally invariant subspaces in part i of this paper almost controlled invariant subspaces were studied in this part we will consider their duals the almost conditionally invariant subspaces these concepts give immediately by dualization of the almost disturbance decoupling control by state feedback the solution of the almost disturbance decoupled estimation problem finally we consider the problem of approximate disturbance decoupling by measurement feedback and it is shown that this problem is solvable to any arbitrary degree of accuracy if and only if 1 almost disturbance decoupling by state feedback and 2 almost disturbance decoupled estimation of the to be controlled output are both possible

almost disturbance decoupling by a proportional derivative state feedback law abstract in this paper the almost disturbance decoupling problem addp for a multivariable linear system is studied by means of a proportional derivative pd state feedback law u fx r dot x a condition for solvability of the addp by state feedback has been given by willems ieee trans aut control ac 26 235 1981 who also indicated that the solution of the addp requires high gain state feedback f it is shown here that under the condition given by willems the addp also admits a solution by a pd law involving finite maps it is also shown that the type of pd law considered here preserves some structural properties of linear systems

fast and accurate computation of hypersingular integrals in galerkin surface integral equation formulations via the direct evaluation method hypersingular 4 d integrals arising in the galerkin discretization of surface integral equation formulations are computed by means of the direct evaluation method the proposed scheme extends the basic idea of the singularity cancellation methods usually employed for the regularization of the singular integral kernel by utilizing a series of coordinate transformations combined with a reordering of the integrations the overall algebraic manipulation results in smooth 2 d integrals that can be easily evaluated via standard quadrature rules finally the reduction of the dimensionality of the original integrals together with the smooth behavior of the associated integrands lead up to unmatched accuracy and efficiency

feedforward control pid control laws and almost invariant subspaces we define feedforward control as a control policy in which the exogenous disturbances are known for all time at the moment when the control is applied it is shown that disturbance decoupling by feedforward control is possible iff it is possible by pid control or iff approximate disturbance decoupling by state feedback is possible

finite element solution of the helmholtz equation with high wave number part i the h version of the fem abstract the paper addresses the properties of finite element solutions for the helmholtz equation the h version of the finite element method with piecewise linear approximation is applied to a one dimensional model problem new results are shown on stability and error estimation of the discrete model in all propositions assumptions are made on the magnitude of hk only where k is the wavelength and h is the stepwidth of the fe mesh previous analytical results had been shown with the assumption that k 2 h is small for medium and high wavenumber these results do not cover the meshsizes that are applied in practical applications the main estimate reveals that the error in h 1 norm of discrete solutions for the helmholtz equation is polluted when k 2 h is not small the error is then not quasioptimal i e the relation of the fe error to the error of best approximation generally depends on the wavenumber k it is noted that the pollution term in the relative error is of the same order as the phase lead of the numerical solution in the result of this analysis thorough and rigorous understanding of error behavior throughout the range of convergence is gained numerical results are presented that show sharpness of the error estimates and highlight some phenomena of the discrete solution behavior the h p version of the fem is studied in part ii

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

economical symmetrical quadrature rules for complete polynomials over a square domain it is of interest in numerical analysis to develop symmetrical quadrature rules for integration of complete polynomial functions over a square domain with minimum computational effort gaussian product quadrature rules integrate such functions with maximum effort symmetrical quadrature rules are developed and presented for integration of complete polynomial functions up to 21st order with minimum computational effort

extra dof free and linearly independent enrichments in gfem abstract a gfem without extra dof is developed without the extra dof the resulting linear system becomes independent in size of the order of local approximation the resulting pu approximation becomes free of linear dependence and the resulting stiffness matrix becomes good conditioned numerical studies show the new gfem s excellent convergence properties and excellent stability

integration of singular enrichment functions in the generalized extended finite element method for three dimensional problems a mapping method is developed to integrate weak singularities which result from enrichment functions in the generalized extended finite element method the integration scheme is applicable to 2d and 3d problems including arbitrarily shaped triangles and tetrahedra implementation of the proposed scheme in existing codes is straightforward numerical examples for 2d and 3d problems demonstrate the accuracy and convergence properties of the technique copyright 2008 john wiley sons ltd

implicit solution of the unsteady euler equations for high order accurate discontinuous galerkin discretizations efficient solution techniques for high order accurate time dependent problems are investigated for solving the two dimensional non linear euler equations in this work the spatial discretization consists of a high order accurate discontinuous galerkin dg approach implicit time integration techniques are considered exclusively in order to avoid the stability restrictions of explicit methods standard backwards differencing methods bdf1 and bdf2 as well as a second order crank nicholson cn2 and a fourth order implicit runge kutta irk4 scheme are considered in an attempt to balance the spatial and temporal accuracy of the overall approach the implicit system arising at each time step is solved using a p multigrid approach which is shown to produce h independent convergence rates while remaining relatively insensitive to the time step size the crank nicholson methodology although not l stable demonstrates superior performance compared to the bdf2 scheme for the problems chosen in this work however the fourth order accurate implicit runge kutta scheme is found to be the most efficient in terms of computational cost for a given accuracy level as compared to the lower order schemes in spite of the added cost per time step and the benefits of this scheme increase for tighter error tolerances

application of group representation theory to derive hermite interpolation polynomials on a triangle current methods used to devise sets of hermite interpolation polynomials of minimal order that ensure c n continuity across triangular element boundaries in two dimensions are not readily extensible to higher dimensions the extension of such methods is especially difficult when the number of degrees of freedom afforded by data at points is different from the number of degrees of freedom determined by the coefficients of a complete polynomial basis to a particular order this work introduces a formalism based on group representation theory that can accomplish this task in general the method is introduced through the derivation of c 1 continuous hermite polynomials that interpolate data at the three vertices of an equilateral triangular element these interpolation polynomials are reported here for the first time the polynomials derived here are compared to the standard polynomials defined in a right triangle by using the two sets of polynomials to solve the laplace equation over finite elements the methodology presented here is of use in higher dimensional elements when the complete polynomial degrees of freedom exceed the total c n degrees of freedom at the nodes

spectral element method for elastic and acoustic waves in frequency domain numerical techniques in time domain are widespread in seismic and acoustic modeling in some applications however frequency domain techniques can be advantageous over the time domain approach when narrow band results are desired especially if multiple sources can be handled more conveniently in the frequency domain moreover the medium attenuation effects can be more accurately and conveniently modeled in the frequency domain in this paper we present a spectral element method sem in frequency domain to simulate elastic and acoustic waves in anisotropic heterogeneous and lossy media the sem is based upon the finite element framework and has exponential convergence because of the use of gll basis functions the anisotropic perfectly matched layer is employed to truncate the boundary for unbounded problems compared with the conventional finite element method the number of unknowns in the sem is significantly reduced and higher order accuracy is obtained due to its spectral accuracy to account for the acoustic solid interaction the domain decomposition method ddm based upon the discontinuous galerkin spectral element method is proposed numerical experiments show the proposed method can be an efficient alternative for accurate calculation of elastic and acoustic waves in frequency domain

stability and accuracy of finite element methods for flow acoustics i general theory and application to one dimensional propagation the dispersion properties of finite element models for aeroacoustic propagation based on the convected scalar helmholtz equation and on the galbrun equation are examined the current study focusses on the effect of the mean flow on the dispersion and amplitude errors present in the discrete numerical solutions a general two dimensional dispersion analysis is presented for the discrete problem on a regular unbounded mesh and results are presented for the particular case of one dimensional acoustic propagation in which the wave direction is aligned with the mean flow the magnitude and sign of the mean flow is shown to have a significant effect on the accuracy of the numerical schemes quadratic helmholtz elements in particular are shown to be much less effective for downstream as opposed to upstream propagation even when the effect of wave shortening or elongation due to the mean flow is taken into account these trends are also observed in solutions obtained for simple test problems on finite meshes a similar analysis of two dimensional propagation is presented in an accompanying article

from h to p efficiently implementing finite and spectral hp element methods to achieve optimal performance for low and high order discretisations the spectral hp element method can be considered as bridging the gap between the traditionally low order finite element method on one side and spectral methods on the other side consequently a major challenge which arises in implementing the spectral hp element methods is to design algorithms that perform efficiently for both low and high order spectral hp discretisations as well as discretisations in the intermediate regime in this paper we explain how the judicious use of different implementation strategies can be employed to achieve high efficiency across a wide range of polynomial orders furthermore based upon this efficient implementation we analyse which spectral hp discretisation which specific combination of mesh size h and polynomial order p minimises the computational cost to solve an elliptic problem up to a predefined level of accuracy we investigate this question for a set of both smooth and non smooth problems

discontinuous galerkin methods with plane waves for time harmonic problems a general framework for discontinuous galerkin methods in the frequency domain with numerical flux is presented the main feature of the method is the use of plane waves instead of polynomials to approximate the solution in each element the method is formulated for a general system of linear hyperbolic equations and is applied to problems of aeroacoustic propagation by solving the two dimensional linearized euler equations it is found that the method requires only a small number of elements per wavelength to obtain accurate solutions and that it is more efficient than high order drp schemes in addition the conditioning of the method is found to be high but not critical in practice it is shown that the ultra weak variational formulation is in fact a subset of the present discontinuous galerkin method a special extension of the method is devised in order to deal with singular solutions generated by point sources like monopoles or dipoles aeroacoustic problems with non uniform flows are also considered and results are presented for the sound radiated from a two dimensional jet

entropy stable shock capturing space time discontinuous galerkin schemes for systems of conservation laws we present a streamline diffusion shock capturing spacetime discontinuous galerkin dg method to approximate nonlinear systems of conservation laws in several space dimensions the degrees of freedom are in terms of the entropy variables and the numerical flux functions are the entropy stable finite volume fluxes we show entropy stability of the formally arbitrarily high order accurate method for a general system of conservation laws furthermore we prove that the approximate solutions converge to the entropy measure valued solutions for nonlinear systems of conservation laws convergence to entropy solutions for scalar conservation laws and for linear symmetrizable systems is also shown numerical experiments are presented to illustrate the robustness of the proposed schemes

a paradigm for higher order polygonal elements in finite elasticity using a gradient correction scheme abstract recent studies have demonstrated that polygonal elements possess great potential in the study of nonlinear elastic materials under finite deformations on the one hand these elements are well suited to model complex microstructures e g particulate microstructures and microstructures involving different length scales and incorporating periodic boundary conditions on the other hand polygonal elements are found to be more tolerant to large localized deformations than the standard finite elements and to produce more accurate results in bending and shear with mixed formulations lower order mixed polygonal elements are also shown to be numerically stable on voronoi type meshes without any additional stabilization treatment however polygonal elements generally suffer from persistent consistency errors under mesh refinement with the commonly used numerical integration schemes as a result non convergent finite element results typically occur which severely limit their applications in this work a general gradient correction scheme is adopted that restores the polynomial consistency by adding a minimal perturbation to the gradient of the displacement field with the correction scheme the recovery of optimal convergence for solutions of displacement based and mixed formulations with both linear and quadratic displacement interpolants is confirmed by numerical studies of several boundary value problems in finite elasticity in addition for mixed polygonal elements the various choices of the pressure field approximations are discussed and their performance on stability and accuracy are numerically investigated we present applications of those elements in physically based examples including a study of filled elastomers with interphasial effect and a qualitative comparison with cavitation experiments for fiber reinforced elastomers

the generation of triangular meshes for nurbs enhanced fem summary this paper presents the first method that enables the fully automatic generation of triangular meshes suitable for the so called non uniform rational b spline nurbs enhanced finite element method nefem the meshes generated with the proposed approach account for the computer aided design boundary representation of the domain given by nurbs curves the characteristic element size is completely independent of the geometric complexity and of the presence of very small geometric features the proposed strategy allows to circumvent the time consuming process of de featuring complex geometric models before a finite element mesh suitable for the analysis can be produced a generalisation of the original definition of a nefem element is also proposed enabling to treat more complicated elements with an edge defined by several nurbs curves or more than one edge defined by different nurbs three examples of increasing difficulty demonstrate the applicability of the proposed approach and illustrate the advantages compared with those of traditional finite element mesh generators finally a simulation of an electromagnetic scattering problem is considered to show the applicability of the generated meshes for finite element analysis 2016 the authors international journal for numerical methods in engineering published by john wiley sons ltd

efficient symmetrical cubature rules for complete polynomials of high degree over the unit cube the fully symmetrical cubature rules developed and presented are at present the most economically available for numerical integration over the unit cube of complete polynomial functions up to and including nineteenth order conventional gaussian triple product integration formulae by comparison require more than twice as many integration points and weights

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

symmetric triangle quadrature rules for arbitrary functions abstract despite extensive research on symmetric polynomial quadrature rules for triangles as well as approaches to their calculation few studies have focused on non polynomial functions particularly on their integration using symmetric triangle rules in this paper we present two approaches to computing symmetric triangle rules for singular integrands by developing rules that can integrate arbitrary functions the first approach is well suited for a moderate amount of points and retains much of the efficiency of polynomial quadrature rules the second approach better addresses large amounts of points though it is less efficient than the first approach we demonstrate the effectiveness of both approaches on singular integrands which can often yield relative errors two orders of magnitude less than those from polynomial quadrature rules

matrix pencil characterization of almost a z invariant subspaces a classification of geometric concepts the equivalence between the algebraic matrix pencil characterization of the sub spaces of the extended geometric theory and their dynamic characterization 13 established as a result a complete classification of almost a b invariant a b invariant almost controllability and controllability subspaces is derived in terms of matrix pencil invariants the frequency propagation aspects of infinite spectrum a b invariant subspaces are investigated and it is shown that they are limits of closed loop eigenspaces with arbitrarily largo eigenvalues finally the importance of the infinite frequency subspaces in the study of the asymptotic behaviour of the closed loop eigenspaces and eigenvalues under scalar gain output feedback is discussed

volumetric testing parallel to the boundary surface for a nonconforming discretization of the electric field integral equation the volumetric monopolar rwg discretization of the electric field integral equation efie imposes no continuity constraint across edges in the surface discretization around a closed conductor the current is expanded with the monopolar rwg set and the electric field is tested over a set of tetrahedral elements attached to the boundary surface this scheme is facet oriented and therefore well suited for the scattering analysis of nonconformal meshes or composite objects the observed accuracy though is only competitive with respect to the rwg discretization for a restricted range of heights of the tetrahedral elements in this communication we introduce a novel implementation of the volumetric monopolar rwg discretization of the efie with testing over a set of wedges we show with rcs and near field results that this scheme offers improved accuracy for a wider range of heights than the approach with tetrahedral testing the application of the wedge testing to the even surface odd volumetric monopolar rwg discretization of the efie edge oriented and therefore less versatile shows similar accuracy as with tetrahedral testing which is a sign of robustness

symmetric quadrature rules for tetrahedra based on a cubic close packed lattice arrangement a family of quadrature rules for integration over tetrahedral volumes is developed the underlying structure of the rules is based on the cubic close packed ccp lattice arrangement using 1 4 10 20 35 and 56 quadrature points the rules are characterized by rapid convergence positive weights and symmetry each rule is an optimal approximation in the sense that lower order terms have zero contribution to the truncation error and the leading order error term is minimized quadrature formulas up to order 9 are presented with relevant numerical examples

the partition of unity finite element method for short wave acoustic propagation on non uniform potential flows a novel numerical method is proposed for modelling time harmonic acoustic propagation of short wavelength disturbances on non uniform potential flows the method is based on the partition of unity finite element method in which a local basis of discrete plane waves is used to enrich the conventional finite element approximation space the basis functions are local solutions of the governing equations they are able to represent accurately the highly oscillatory behaviour of the solution within each element while taking into account the convective effect of the flow and the spatial variation in local sound speed when the flow is non uniform many wavelengths can be included within a single element leading to ultra sparse meshes results presented in this article will demonstrate that accurate solutions can be obtained in this way for a greatly reduced number of degrees of freedom when compared to conventional element or grid based schemes numerical results for lined uniform two dimensional ducts and for non uniform axisymmetric ducts are presented to indicate the accuracy and performance which can be achieved numerical studies indicate that the pollution effect associated with cumulative dispersion error in conventional finite element schemes is largely eliminated

special short wave elements for flow acoustics a numerical formulation based on the partition of unity method pum is proposed for modelling the propagation of short acoustic waves on irrotational mean flows the method seeks to reduce the pollution error which exists in conventional fe schemes at high frequencies by using a local basis which is enriched by plane wave solutions of the convected helmholtz equation initially the method is demonstrated with reference to a one dimensional model consisting of a variable area converging diverging duct with mean flow next a simple two dimensional model of a straight duct with uniform flow is considered in both cases the accuracy and the conditioning of the numerical solution is investigated for ranges of frequency and mach number characteristic of aero engine bypass ducts

a finite element semi lagrangian method with l2 interpolation summary high order accurate methods for convection dominated problems have the potential to reduce the computational effort required for a given order of solution accuracy the state of the art in this field is more advanced for eulerian methods than for semi lagrangian slag methods in this paper we introduce a new slag method that is based on combining the modified method of characteristics with a high order interpolating procedure the method employs the finite element method on triangular meshes for the spatial discretization an l2 interpolation procedure is developed by tracking the feet of the characteristic lines from the integration nodes numerical results are illustrated for a linear advection diffusion equation with known analytical solution and for the viscous burgers equation the computed results support our expectations for a robust and highly accurate finite element slag method copyright 2012 john wiley sons ltd

on the crossing rule abstract the paper gives conditions on a family of matrices which guarantee that some matrix in the family will have a multiple eigenvalue in particular the main theorem states exactly which dimensions admit k dimensional subspaces of matrices for which all nonzero elements have distinct eigenvalues this question arises naturally in the theory of first order hyperbolic systems of partial differential equations the main theorem in this context tells exactly for which integers n an n x n system in k space variables may be strictly hyperbolic author

dispersion analysis and error estimation of galerkin finite element methods for the helmholtz equation when applying numerical methods for the computation of stationary waves from the helmholtz equation one obtains numerical waves that are dispersive also in non dispersive media the numerical wave displays a phase velocity that depends on the parameter k of the helmholtz equation in dispersion analysis the phase difference between the exact and the numerical solutions is investigated in this paper the authors recent result on the phase difference for one dimensional problems is numerically evaluated and discussed in the context of other work directed to this topic it is then shown that previous error estimates in h1 norm are of nondispersive character but hold for medium or high wavenumber on extremely refined mesh only on the other hand recently proven error estimates for constant resolution contain a pollution term with certain assumptions on the exact solution this term is of the order of the phase difference thus a link is established between the results of dispersion analysis and the results of numerical analysis throughout the paper the presentation and discussion of theoretical results is accompanied by numerical evaluation of several model problems special attention is given to the performance of the galerkin method with a higher order of polynomial approximation p h p version

an accurate approach for the calculation of mom matrix elements the method of moments solution on integral equations for electromagnetic scattering and radiation problems with rao wilton glisson rwg and n spl times rwg basis functions requires the calculation of singular double surface integrals over triangular subdomains in this paper an approach is presented to deal with the singularity and almost singularity in which the triangle subdomain is divided into several subtriangles and then the singular or almost singular integral is calculated over each of the subtriangles this approach is mainly based on the projection process and the local polar frame it can be used to deal with the singular double surface integrals of order 1 r and 1 r sup 2 in which vector functions may be of zeroth order and higher order numerical examples are provided to verify the validity of the method presented in this paper for dealing with the singularity and almost singularity

high degree efficient symmetrical gaussian quadrature rules for the triangle gaussian quadrature is required for the computation of matrices based on the isoparametric formulztion of the finite element method a brief review of existing quadrature rules for the triangle is given and the method for the determination of high degree efficient symmetrical rules for the triangle is discussed new quadrature rules of degree 12 20 are presented and a short fortran program is included

generic non selfadjoint zakharov shabat operators in this paper we develop tools to study within a family of non selfadjoint operators l varphi depending on a parameter varphi in a real hilbert space those with partially simple spectrum as a case study we consider the zakharov shabat operators l varphi appearing in the lax pair of the focusing nls on the circle in particular the main result implies that the set of potentials varphi of sobolev class h n n ge 0 so that all non real eigenvalues of l varphi are simple is path connected and dense

on the application of two symmetric gauss legendre quadrature rules for composite numerical integration over a triangular surface this paper first presents a gauss legendre quadrature rule for the evaluation of i t f x y d x d y where f x y is an analytic function in x y and t is the standard triangular surface x y 0 x y 1 x y 1 in the two space x y we transform this integral into an equivalent integral s f x y frac x y d d where s is the 2 square in space 1 1 we then apply the one dimensional gauss legendre quadrature rules in and variables to arrive at an efficient quadrature rules with new weight coefficients and new sampling points then a second gauss legendre quadrature rule of composite type is obtained this rule is derived by discretising t into three new triangles tic i 1 2 3 of equal size which are obtained by joining centroid of t c 1 3 1 3 to the three vertices of t by use of affine transformations defined over each tic and the linearity property of integrals leads to the result i underover i 1 3 tic f x y d x d y frac 1 3 t g x y d x d y where g x y i 1n n f xic x y yic x y and x xic x y and y yic x y refer to affine transformations which map each tic into t the standard triangular surface we then write t g x y d x d y s g x y frac x y d d and a composite rule of integration is thus obtained we next propose the discretisation of the standard triangular surface t into n2 right isosceles triangular surfaces ti i 1 1 n2 each of which has an area equal to 1 2 n2 units we have again shown that the use of affine transformation over each ti and the use of linearity property of integrals lead to the result t f x y d x d y underover i 1 n n ti f x y d x d y frac 1 n2 t h x y d x d y where h x y i 1n n f xi x y yi x y and x xi x y y yi x y refer to affine transformations which map each ti in x y space into t a standard triangular surface t in the x y space we can now apply the two rules earlier derived to the integral t h x y d x d y this amounts to application of composite numerical integration of t into n2 and 3n2 triangles of equal sizes respectively we can now apply the rules which are derived earlier to the evaluation of the integral t f x y d x d y and each of these procedures converges to the exact value of the integral t f x y d x d y for sufficiently large value of n and the convergence is much faster for higher order rules we have demonstrated this aspect by applying the above composite integration method to some typical integrals 2007 elsevier inc all rights reserved

indirect finite element evaluation of two dimensional finite part integral using fourier transformation a broad class of engineering problems in fracture mechanics thermal fluid transport and electromagnetic theory involve the evaluation of two dimensional finite part integrals of the form a method for evaluation of such integrals is developed by deriving an equivalent integral using fourier transformation this equivalent integral does not involve a kernel with singular behaviour consequently standard numerical integration methodologies with conventional analytical evaluation techniques can be used in the finite element computations the accuracy and convergence of the developed numerical procedure are successfully demonstrated by numerical examples for planar fracture geometries

implicit runge kutta methods for the discretization of time domain integral equations implicit runge kutta based schemes are proposed for solving the time domain integral equations of electromagnetic theory the proposed technique maps a laplace domain equation to a domain equation using the butcher tableau of the implicit runge kutta scheme a discrete time domain system is recovered by computing the inverse transform numerically the resulting technique is capable of third or fifth order accuracy in time and is stable independent of time step numerical results illustrate the accuracy and stability of the technique

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

adjoint sensitivity formulation for discontinuous galerkin discretizations in unsteady inviscid flow problems this paper presents an unsteady discrete adjoint algorithm for high order implicit discontinuous galerkin discretizations in time dependent inviscid flow problems the major function of the adjoint approach is to obtain the sensitivity information in a time dependent functional output which in turn is used to drive an unsteady shape optimization process to deliver a minimum of the objective functional a gradient based optimization strategy is investigated in which the sensitivity derivatives of the objective functional with respect to input variables are formulated in the context of high order discontinuous galerkin discretizations while special emphasis is given to the variations and linearizations of curvilinear boundary elements implicit temporal discretizations consisting of a second order backward euler scheme and a fourth order implicit runge kutta scheme are considered exclusively in this work where the corresponding adjoint problem is required to be solved in a backward time integration manner due to the associated transpose operation two numerical examples for the unsteady shape design techniques are presented to verify the derived sensitivity formulations and to demonstrate the performance of the adjoint approach the first involves an inverse shape optimization case by matching a time dependent target pressure profile for a two dimensional periodic vortical gust impinging on a rae 2822 airfoil and the second considers minimization of the acoustic noise produced by subsonic flow over a naca0012 airfoil with a 0 03c thick blunt trailing edge

generalized gaussian quadrature rules for discontinuities and crack singularities in the extended finite element method new gaussian integration schemes are presented for the efficient and accurate evaluation of weak form integrals in the extended finite element method for discontinuous functions we construct gauss like quadrature rules over arbitrarily shaped elements in two dimensions without the need for partitioning the finite element a point elimination algorithm is used in the construction of the quadratures which ensures that the final quadratures have minimal number of gauss points for weakly singular integrands we apply a polar transformation that eliminates the singularity so that the integration can be performed efficiently and accurately numerical examples in elastic fracture using the extended finite element method are presented to illustrate the performance of the new integration techniques

application of an ultra weak variational formulation of elliptic pdes to the two dimensional helmholtz problem a new technique to solve elliptic linear pdes called ultra weak variational formulation uwvf in this paper is introduced in b despres c r acad sci paris 318 1994 pp 939 944 this paper is devoted to an evaluation of the potentialities of this technique it is applied to a model wave problem the two dimensional helmholtz problem the new method is presented in three parts following the same style of presentation as the classical one of the finite elements method even though they are definitely conceptually different methods the first part is committed to the variational formulation and to the continuous problem the second part defines the discretization process using a galerkin procedure the third part actually studies the efficiency of the technique from the order of convergence point of view this is achieved using theoretical proofs and a series of numerical experiments in particular it is proven and shown the order of convergence is lower bounded by a linear function of the number of degrees of freedom an application to scattering problems is presented in a fourth part

generalized gaussian quadrature rules on arbitrary polygons in this paper we present a numerical algorithm based on group theory and numerical optimization to compute efficient quadrature rules for integration of bivariate polynomials over arbitrary polygons these quadratures have desirable properties such as positivity of weights and interiority of nodes and can readily be used as software libraries where numerical integration within planar polygons is required we have used this algorithm for the construction of symmetric and non symmetric quadrature rules over convex and concave polygons while in the case of symmetric quadratures our results are comparable to available rules the proposed algorithm has the advantage of being flexible enough so that it can be applied to arbitrary planar regions for the integration of generalized classes of functions to demonstrate the efficiency of the new quadrature rules we have tested them for the integration of rational polygonal shape functions over a regular hexagon for a relative error of 10 8 in the computation of stiffness matrix entries one needs at least 198 evaluation points when the region is partitioned whereas 85 points suffice with our quadrature rule copyright 2009 john wiley sons ltd

the mini mixed finite element for the stokes problem an experimental investigation abstract o h 3 2 superconvergence in pressure and velocity has been experimentally investigated for the two dimensional stokes problem discretized with the mini mixed finite element even though the classic mixed finite element theory for the mini element guarantees linear convergence for the total error recent theoretical results indicate that superconvergence of order o h 3 2 in pressure and of the linear part of the computed velocity to the piecewise linear nodal interpolation of the exact velocity is in fact possible with structured three directional triangular meshes the numerical experiments presented here suggest a more general validity of o h 3 2 superconvergence possibly to automatically generated and unstructured triangulations in addition the approximating properties of the complete computed velocity have been compared with the approximating properties of the piecewise linear part of the computed velocity finding that the former is generally closer to the exact velocity whereas the latter conserves mass better

a new triangular and tetrahedral basis for high order hp finite element methods in this paper we describe the foundations of a new hierarchical modal basis suitable for high order hp finite element discretizations on unstructured meshes it is based on a generalized tensor product of mixed weight jacobi polynomials the generalized tensor product property leads to a low operation count with the use of sum factorization techniques variable p order expansions in each element are readily implemented which is a crucial property for efficient adaptive discretizations numerical examples demonstrate the exponential convergence for smooth solutions and the ability of this formulation to handle easily very complex two and three dmensional computational domains employing standard meshes

distribution of the discretization and algebraic error in numerical solution of partial differential equations abstract in the adaptive numerical solution of partial differential equations local mesh refinement is used together with a posteriori error analysis in order to equilibrate the discretization error distribution over the domain since the discretized algebraic problems are not solved exactly a natural question is whether the spatial distribution of the algebraic error is analogous to the spatial distribution of the discretization error the main goal of this paper is to illustrate using standard boundary value model problems that this may not hold on the contrary the algebraic error can have large local components which can significantly dominate the total error in some parts of the domain the illustrated phenomenon is of general significance and it is not restricted to some particular problems or dimensions to our knowledge the discrepancy between the spatial distribution of the discretization and algebraic errors has not been studied in detail elsewhere

a new c1 finite element full heptic a new c1 finite element complete in 7th degree polynomial basis is described and compared with other c1 elements in an eigen elements analysis

on the application of two gauss legendre quadrature rules for composite numerical integration over a tetrahedral region in this paper we first present a gauss legendre quadrature rule for the evaluation of i f f t ff x y z dxdydz where f x y z is an analytic function in x y z and tis the standard tetrahedral region x y z 0 x y z i x y z 1 in three space x y z we then use a transformation x x y y and z z to change the integral into an equivalent integral i 1 1 1 1 1 1 f x y x y z d d d over the standard 2 cube in space 1 1 we then apply the one dimensional gauss legendre quadrature rules in and variables to arrive at an efficient quadrature rule with new weight coefficients and new sampling points then a second gauss legendre quadrature rule of composite type is obtained this rule is derived by discretising the tetrahedral region t into four new tetrahedra t c i i 1 2 3 4 of equal size which are obtained by joining the centroid of t c 1 4 1 4 1 4 to the four vertices of t by use of the affine transformations defined over each t c i and the linearity property of integrals leads to the result i 4 i 1 tci f x y z dxdydz 1 4 t g x y z dxdydz where g x y z 1 p 3 f x k x y z y k x y z z k x y z x k x k x y z y k y k x y z and z k z k x y z refer to an affine transformations which map each t c i into the standard tetrahedral region t we then write i t g x y z dxdydz 1 0 1 0 1 0 g x y z x y z d d d and a composite rule of integration is thus obtained we next propose the discretisation of the standard tetrahedral region t into p 3 tetrahedra t i i 1 1 p 3 each of which has volume equal to 1 6p 3 units we have again shown that the use of affine transformations over each t and the use of linearity property of integrals leads to the result t f x y z dxdydz p3 i 1 t c i f x y z dxdydz p3 1 t p f x p y p z p dx p dy p dz p 1 p3 t h x y z dxdydz where p3 h x y z 1 f x p x y z y p x y z z p x y z x p x p x y z y p y p x y z and z p z p x y z refer to the affine transformations which map each t i in x p y p z p space into a standard tetrahedron t in the x y z space we can now apply the two rules earlier derived to the integral t h x y z dxdydz this amounts to the application of composite numerical integration of t into p3 and 4p3 tetrahedra of equal sizes we have demonstrated this aspect by applying the above composite integration method to some typical triple integrals

singularity subtraction technique for high order polynomial vector basis functions on planar triangles in this paper a singularity subtraction technique is developed for computing the impedance matrix elements of various electromagnetic surface integral equation formulations with the galerkin method and high order basis functions analytical closed form formulas for computing surface integrals with r r sup n n spl ges 3 singularities times polynomial nodal shape functions of arbitrary order on a planar triangle are presented

a high order locally corrected nystrom scheme for charge current integral equations a high order nystrom scheme is developed for electromagnetic em scattering problems of smooth objects the implementation is based on different well known charge current formulations this study aims to develop a novel quadrature scheme to compute different singular integrals that appear on integral equation methods in ems the proposed quadrature scheme is used on a locally corrected nystrom method and applied to charge current integral equations to obtain a fast accurate and a low frequency stable method to solve em scattering problems the proposed method allows obtaining accuracy in the range 4 10 digits in the near electric and magnetic scattered field for certain simple geometries and a wide range of frequencies we also do a preliminary research on nonsmooth geometries obtaining a certain degree of accuracy depending on the particular charge current integral equation used

on the assignability of infinite root loci in almost disturbance decoupling in this paper we investigate in detail the almost disturbance decoupling problem with state feedback an alternative proof of the sufficiency part of the original result by willoms 1981 on this problem requiring the construction of a sequence of feedback matrices is given it is shown that one can achieve almost disturbance decoupling by letting some of the closed loop eigenvalues run to infinity according to a spatial arrangement determined by the infinite zero structure the proofs of our results use the fact that almost controllability subspaces can be approximated by controlled invariant subspaces our results on the closed loop eigenstructure are stated in terms of such approximating sequences we also prove a generalization of a theorem by wonham 1979 used in the context of minimal order observer design

a convergence analysis of an h version finite element method with high order elements for two dimensional elasto plasticity problems in this paper we will give an h version finite element method for a two dimensional nonlinear elasto plasticity problem a family of admissible constitutive laws based on the so called gauge function method is introduced first and then a high order h version semidiscretization scheme is presented the existence and uniqueness of the solution for the semidiscrete problem are guaranteed by using some special properties of the constitutive law and finally we will show that as the maximum element size h longrightarrow 0 the solution of the semidiscrete problem will converge to the solution of the continuous problem the high order h version discretization scheme introduced here is unusual if the partition of the spatial space only has rectangles or parallelograms involved then there would not be any limit on the element degree however if the partition of the spatial space has some triangular elements then only certain combinations of finite element spaces for displacement and stress functions can be used the discretization scheme also provides a useful idea for applications of hp version or high order h version finite element methods for two dimensional problems where the elasto plastic body is not a polygon such as a disk or an annulus

assessment of finite and spectral element shape functions for efficient iterative simulations of interior acoustics several shape function families currently used in the spectral element and the finite element method fem are reviewed widely known shape functions such as the conventional lagrange functions various p fem shapes and recent developments in the area of spectral methods are compared especially shape functions derived from bernstein polynomials are discussed due to their ease of formulation for almost any geometric entity their applicability and efficiency for acoustic simulations are assessed in particular the efficiency is measured in terms of performance of current krylov solvers

hyperbolic conservation laws with umbilic degeneracy i in this paper a compactness framework for approximate solutions to nonlinear hyperbolic systems with umbilic degeneracy is established by combining techniques of compensated compactness with some classical methods and by a detailed analysis of a highly singular equation of euler poisson darboux type then this framework is successfully applied to prove the convergence of the viscosity method and to prove the existence of global entropy solutions for the cauchy problem with large initial data for a canonical class of the systems with quadratic flux form

mixed inner outer iteration technique based surface integral equations for fast solving em scattering from penetrable objects in this paper a mixed inner outer iteration technique based on surface integral equation is proposed for solving scattering problems from penetrable and homogeneous dielectric structures the most important novelty of the proposed technique is combining the tangential poggio miller chang harrington wu tsai as outer iteration and the electric magnetic current combined field integral equation as inner iteration respectively in this case the new method offers the guarantee of high accuracy with faster convergence especially for electrically dielectric structures in large scale furthermore it presents an excellent compatibility with sparse approximate inverse preconditioning numerical results are given to demonstrate the accuracy efficiency and compatibility of the proposed method

second order accurate derivatives and integration schemes for meshfree methods summary the consistency condition for the nodal derivatives in traditional meshfree galerkin methods is only the differentiation of the approximation consistency dac one missing part is the consistency between a nodal shape function and its derivatives in terms of the divergence theorem in numerical forms in this paper a consistency framework for the meshfree nodal derivatives including the dac and the discrete divergence consistency ddc is proposed the summation of the linear ddc over the whole computational domain leads to the so called integration constraint in the literature a three point integration scheme using background triangle elements is developed in which the corrected derivatives are computed by the satisfaction of the quadratic ddc we prove that such smoothed derivatives also meet the quadratic dac and therefore the proposed scheme possesses the quadratic consistency that leads to its name qc3 numerical results show that qc3 is the only method that can pass both the linear and the quadratic patch tests and achieves the best performances for all the four examples in terms of stability convergence accuracy and efficiency among all the tested methods particularly it shows a huge improvement for the existing linearly consistent one point integration method in some examples copyright 2012 john wiley sons ltd

proper generalized decomposition for parameterized helmholtz problems in heterogeneous and unbounded domains application to harbor agitation abstract solving the helmholtz equation for a large number of input data in an heterogeneous media and unbounded domain still represents a challenge this is due to the particular nature of the helmholtz operator and the sensibility of the solution to small variations of the data here a reduced order model is used to determine the scattered solution everywhere in the domain for any incoming wave direction and frequency moreover this is applied to a real engineering problem water agitation inside real harbors for low to mid high frequencies the proper generalized decomposition pgd model reduction approach is used to obtain a separable representation of the solution at any point and for any incoming wave direction and frequency here its applicability to such a problem is discussed and demonstrated more precisely the contributions of the paper include the pgd implementation into a perfectly matched layer framework to model the unbounded domain and the separability of the operator which is addressed here using an efficient higher order projection scheme then the performance of the pgd in this framework is discussed and improved using the higher order projection and a petrov galerkin approach to construct the separated basis moreover the efficiency of the higher order projection scheme is demonstrated and compared with the higher order singular value decomposition

monomial cubature rules since stroud a compilation a bibliography of references to cubature rules which have appeared since the publication of stroud s book 1971 is presented the standard regions that are treated in this paper are the n cube the n simplex the n sphere and the entire space

on the adaptive coupling of fem and bem in 2 d elasticity this paper concerns the combination of the finite element method fem and the boundary element method bem using the symmetric coupling as a model problem in two dimensions we consider the hencky material a certain nonlinear elastic material in a bounded domain with navier lame differential equation in the unbounded complementary domain using some boundary integral operators the problem is rewritten such that the galerkin procedure leads to a fem bem coupling and quasi optimally convergent discrete solutions beside this a priori information we derive an a posteriori error estimate which allows up to a constant factor the error control in the energy norm since information about the singularities of the solution is not available a priori in many situation and having in mind the goal of an automatic mesh refinement we state adaptive algorithms for the h version of the fem bem coupling illustrating numerical results are included

nonstrictly hyperbolic nonlinear systems on etudie le caractere bien pose du probleme aux valeurs initiales localement en temps pour une classe de systemes hyperboliques non lineaires

information transfer between incompatible finite element meshes application to coupled thermo viscoelasticity this article is concerned with information transfer between non matching finite element meshes such a feature is not new in the literature but we focus herein on a geometric approach to transfer solution fields in order to be as application independent as possible moreover the dedicated case concerns the transfer of finite element fields defined at integration points of the meshes and allows iterative exchange of fields in both directions to do so we propose an extension to the mortar technique that fulfills these goals and that does not suffer from a high computational cost the application is the simulation of a strongly coupled thermo viscoelastic problem with phase transition solved with a partitioning technique

numerical approximation of internal discontinuity interface problems this work focuses on the finite element discretization of boundary value problems whose solution features either a discontinuity or a discontinuous conormal derivative across an interface inside the computational domain the interface is characterized via a level set function the discontinuities are accounted for by using suitable extension operators whose numerical implementation requires a very low computational effort after carrying out the error analysis numerical results to validate our approach are presented in one two and three dimensions

general behavior and shear center location of prismatic anisotropic beams via power series abstract the behavior of a tip loaded anisotropic cantilever beam with an arbitrary cross section is studied using saint venant s semi inverse method along with a power series solution for the local in plane and out of plane deformation warping functions the power series coefficients are determined by solving a set of variationally derived linear algebraic equations using the resulting three dimensional displacement solutions the shear deformation associated with applied tip loads is investigated as well as the shear center location two different definitions of the shear center are presented for anisotropic beams by extending existing approaches developed for isotropic sections both of the extended definitions reveal the linear dependency of the shear center location on beam length numerical results are presented for three different cross sections ellipse triangle naca 0012 airfoil and two different materials al 6061 t6 off angle high strength graphite epoxy fibers

on structural shape optimization using an embedding domain discretization technique summary this contribution presents a novel approach to structural shape optimization that relies on an embedding domain discretization technique the evolving shape design is embedded within a uniform finite element background mesh which is then used for the solution of the physical state problem throughout the course of the optimization we consider a boundary tracking procedure based on adaptive mesh refinement to separate between interior elements exterior elements and elements intersected by the physical domain boundary a selective domain integration procedure is employed to account for the geometric mismatch between the uniform embedding domain discretization and the evolving structural component thereby we avoid the need to provide a finite element mesh that conforms to the structural component for every design iteration as it is the case for a standard lagrangian approach to structural shape optimization still we adopt an explicit shape parametrization that allows for a direct manipulation of boundary vertices for the design evolution process in order to avoid irregular and impracticable design updates we consider a geometric regularization technique to render feasible descent directions for the course of the optimization this article is protected by copyright all rights reserved

2 adic valuations of certain ratios of products of factorials and applications abstract we prove the conjecture of falikman friedland loewy on the parity of the degrees of projective varieties of n n complex symmetric matrices of rank at most k we also characterize the parity of the degrees of projective varieties of n n complex skew symmetric matrices of rank at most 2 p we give recursive relations which determine the parity of the degrees of projective varieties of m n complex matrices of rank at most k in the case the degrees of these varieties are odd we characterize the minimal dimensions of subspaces of n n skew symmetric real matrices and of m n real matrices containing a nonzero matrix of rank at most k the parity questions studied here are also of combinatorial interest since they concern the parity of the number of plane partitions contained in a given box on the one hand and the parity of the number of symplectic tableaux of rectangular shape on the other hand

fast matlab assembly of fem matrices in 2d and 3d we propose an effective and flexible way to assemble finite element stiffness and mass matrices in matlab we apply this for problems discretized by edge finite elements typical edge finite elements are raviart thomas elements used in discretizations of h div spaces and nedelec elements in discretizations of h curl spaces we explain vectorization ideas and comment on a freely available matlab code which is fast and scalable with respect to time

the discontinuous enrichment method we propose a finite element based discretization method in which the standard polynomial field is enriched within each element by a nonconforming field that is added to it the enrichment contains free space solutions of the homogeneous differential equation that are not represented by the underlying polynomial field continuity of the enrichment across element interfaces is enforced weakly by lagrange multipliers in this manner we expect to attain high coarse mesh accuracy without significant degradation of conditioning assuring good performance of the computation at any mesh resolution examples of application to acoustics and advection diffusion are presented

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

edged based smoothed point interpolation method for acoustic radiation with perfectly matched layer abstract it is known that for exterior acoustics the domain based finite element method always encounters the accuracy loss for high wavenumbers and the difficulty in dealing with the sommerfeld radiation condition in this work by using the polynomial and radial bases with specific node selection schemes the point interpolation method with the generalized gradient smoothing technique is employed to improve the accuracy for acoustic computations the appropriately softened acoustical stiffness of the discrete model is constructed based on the edge based smoothing domains via the gradient smoothing technique in addition the perfectly matched layer technique enables the smoothed point interpolation method to investigate exterior acoustics the complex coordinate transformation is utilized to determine the proper parameters for the perfectly matched layer through theoretical and numerical analyses several typical acoustic problems are analyzed revealing that better accuracy is obtained by the present method compared with the finite element method especially for high wavenumbers

computation of moderate degree fully symmetric cubature rules on the triangle using symmetric polynomials and algebraic solving a novel method is presented for expressing the moment equations involved in computing fully symmetric cubature rules on the triangle by using symmetric polynomials to represent the two kinds of invariance inherent in these rules this method results in a system of polynomial equations that is amenable to solution using algebraic solving techniques using grobner bases rules of degree up to 15 are computed and presented some of them new and with all their points inside the triangle since all solutions to the polynomial system are computed it is for the first time possible to prove whether a given rule type results in specific rules of a given quality it is thus proved that for degrees up to 14 there are no non fortuitous rules that can improve on the presented results for degree 10 an example is also provided showing how the proposed method can be used to exclude the existence of better fortuitous rules as well

nurbs enhanced finite element method nefem the development of nurbs enhanced finite element method nefem is revisited this technique allows a seamless integration of the cad boundary representation of the domain and the finite element method fem the importance of the geometrical model in finite element simulations is addressed and the benefits and potential of nefem are discussed and compared with respect to other curved finite element techniques

on numerical integration in isogeometric subdivision methods for pdes on surfaces abstract subdivision surfaces offer great flexibility in capturing irregular topologies combined with higher order smoothness for instance loop and catmull clark subdivision schemes provide c 2 smoothness everywhere except at extraordinary vertices where the generated surfaces are c 1 smooth the combination of flexibility and smoothness leads to their frequent use in geometric modeling and makes them an ideal candidate for performing isogeometric analysis of higher order problems with irregular topologies e g thin shell problems in the neighborhood of extraordinary vertices however the resulting surfaces and the functions defined on them are non polynomial thus numerical integration based on quadrature has to be performed carefully to preserve the expected higher order consistency this paper presents a detailed case study of different quadrature schemes for isogeometric discretizations of partial differential equations on closed surfaces with loop s subdivision scheme as model problems we consider elliptic equations with the laplace beltrami and the surface bi laplacian operator as well as the associated eigenvalue problem for the laplace beltrami operator a particular emphasis is on the robustness of the approach in the vicinity of extraordinary vertices based on a series of numerical experiments different quadrature schemes are compared mid edge quadrature which can easily be implemented via a lookup table turns out to be a preferable choice due to its robustness and efficiency

analysis of high order finite elements for convected wave propagation in this paper we examine the performance of high order finite element methods fem for aeroacoustic propagation based on the convected helmholtz equation a methodology is presented to measure the dispersion and amplitude errors of the p fem including non interpolating shape functions such as bubble shape functions a series of simple test cases are also presented to support the results of the dispersion analysis the main conclusion is that the properties of p fem that make its strength for standard acoustics e g exponential p convergence low dispersion error remain present for flow acoustics as well however the flow has a noticeable effect on the accuracy of the numerical solution even when the change in wavelength due to the mean flow is accounted for and an approximation of the dispersion error is proposed to describe the influence of the mean flow also discussed is the so called aliasing effect which can reduce the accuracy of the solution in the case of downstream propagation this can be avoided by an appropriate choice of mesh resolution

machine precision evaluation of singular and nearly singular potential integrals by use of gauss quadrature formulas for rational functions a new technique for machine precision evaluation of singular and nearly singular potential integrals with 1 r singularities is presented the numerical quadrature scheme is based on a new rational expression for the integrands obtained by a cancellation procedure in particular by using library routines for gauss quadrature of rational functions readily available in the literature this new expression permits the exact numerical integration of singular static potentials associated with polynomial source distributions the rules to achieve the desired numerical accuracy for singular and nearly singular static and dynamic potential integrals are presented and discussed and several numerical examples are provided

adjoint based h p adaptive discontinuous galerkin methods for the 2d compressible euler equations in this paper we investigate and present an adaptive discontinuous galerkin algorithm driven by an adjoint based error estimation technique for the inviscid compressible euler equations this approach requires the numerical approximations for the flow i e primal problem and the adjoint i e dual problem which corresponds to a particular simulation objective output of interest the convergence of these two problems is accelerated by an hp multigrid solver which makes use of an element gauss seidel smoother on each level of the multigrid sequence the error estimation of the output functional results in a spatial error distribution which is used to drive an adaptive refinement strategy which may include local mesh subdivision h refinement local modification of discretization orders p enrichment and the combination of both approaches known as hp refinement the selection between h and p refinement in the hp adaptation approach is made based on a smoothness indicator applied to the most recently available flow solution values numerical results for the inviscid compressible flow over an idealized four element airfoil geometry demonstrate that both pure h refinement and pure p enrichment algorithms achieve equivalent error reductions at each adaptation cycle compared to a uniform refinement approach but requiring fewer degrees of freedom the proposed hp adaptive refinement strategy is capable of obtaining exponential error convergence in terms of degrees of freedom and results in significant savings in computational cost a high speed flow test case is used to demonstrate the ability of the hp refinement approach for capturing strong shocks or discontinuities while improving functional accuracy

a spectral finite element approach to modeling soft solids excited with high frequency harmonic loads abstract an approach for efficient and accurate finite element analysis of harmonically excited soft solids using high order spectral finite elements is presented and evaluated the helmholtz type equations used to model such systems suffer from additional numerical error known as pollution when excitation frequency becomes high relative to stiffness i e high wave number which is the case for example for soft tissues subject to ultrasound excitations the use of high order polynomial elements allows for a reduction in this pollution error but requires additional consideration to counteract runge s phenomenon and or poor linear system conditioning which has led to the use of spectral element approaches this work examines in detail the computational benefits and practical applicability of high order spectral elements for such problems the spectral elements examined are tensor product elements i e quad or brick elements of high order lagrangian polynomials with non uniformly distributed gauss lobatto legendre nodal points a shear plane wave example is presented to show the dependence of the accuracy and computational expense of high order elements on wave number then a convergence study for a viscoelastic acoustic structure interaction finite element model of an actual ultrasound driven vibroacoustic experiment is shown the number of degrees of freedom required for a given accuracy level was found to consistently decrease with increasing element order however the computationally optimal element order was found to strongly depend on the wave number