In this paper, mesh-free element free Galerkin (EFG) method is extended to solve two-dimensional potential flow problems. Two ideal fluid flow problems (i.e. flow over a rigid cylinder and flow over a sphere) have been formulated using variational approach. Penalty and Lagrange multiplier techniques have been utilized for the enforcement of essential boundary conditions. Four point Gauss quadra...

An efficient implementation of the high-order spectral volume (SV) method is presented for multi-dimensional conservation laws on unstructured grids. In the SV method, each simplex cell is called a spectral volume (SV), and the SV is further subdivided into polygonal (2D), or polyhedral (3D) control volumes (CVs) to support high-order data reconstructions. In the traditional implementation, Gau...

Computational methods are developed for generating Gauss-type quadrature formulae having nodes of arbitrary multiplicity at one or both end points of the interval of integration. Positivity properties of the boundary weights are investigated numerically, and related conjectures are formulated. Applications are made to moment-preserving spline approximation. AMS subject classification: 65D30.

An efficient implementation of the high-order spectral volume (SV) method is presented for multidimensional conservation laws on unstructured grids. In the SV method, each simplex grid cell is called a spectral volume (SV), and the SV is further subdivided into polygonal (2D), or polyhedral (3D) control volumes (CVs) to support high-order data reconstructions. In the traditional implementation,...

The RK5GL3 method is a numerical method for solving initial value problems in ordinary differential equations, and is based on a combination of a fifth-order Runge-Kutta method and 3-point Gauss-Legendre quadrature. In this paper we describe the propagation of local errors in this method, and show that the global order of RK5GL3 is expected to be six, one better than the underlying RungeKutta m...

This paper presents a Gauss Legendre quadrature method for numerical integration over the standard triangular surface: {(x, y) | 0 , 1, 1} x y x y ≤ ≤ + ≤ in the Cartesian two-dimensional (x, y) space. Mathematical transformation from (x, y) space to (ξ, η) space map the standard triangle in (x, y) space to a standard 2-square in (ξ, η) space: {(ξ, η)|–1 ≤ ξ, η ≤ 1}. This overcomes the difficul...

Let I[f ] = ∫ 1 −1 f(x) dx, where f ∈ C ∞(−1, 1), and let Gn[f ] = ∑n i=1 wnif(xni) be the n-point Gauss–Legendre quadrature approximation to I[f ]. In this paper, we derive an asymptotic expansion as n → ∞ for the error En[f ] = I[f ]−Gn[f ] when f(x) has general algebraic-logarithmic singularities at one or both endpoints. We assume that f(x) has asymptotic expansions of the forms f(x) ∼ ∞ ∑ ...

A HIERARCHICAL SPHERICAL RADIAL QUADRATURE ALGORITHM FOR MULTILEVEL GLMMs, GSMMs, and GENE PATHWAY ANALYSIS SEPTEMBER 2010 JACOB A. GAGNON, B.S., MASSACHUSETTS INSTITUTE OF TECHNOLOGY M.S., UNIVERSITY OF MASSACHUSETTS AT AMHERST Ph.D., UNIVERSITY OF MASSACHUSETTS AT AMHERST Directed by: Professor Anna Liu The first part of my thesis is concerned with estimation for longitudinal data using gener...