A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

On Series of Ordinals and Combinatorics

This paper deals mainly with generalizations of results in finitary combinatorics to infinite ordinals. It is well-known that for finite ordinals ∑_bT<αβ is the number of 2-element subsets of an α-element set. It is shown here that for any well-ordered set of arbitrary infinite order type α, ∑_bT<αβ is the ordinal of the set M of 2-element subsets, where M is ordered in some natural way. The result is then extended to evaluating the ordinal of the set of all n-element subsets for each natural number n ≥ 2. Moreover, series ∑_β<αf(β) are investigated and evaluated, where α is a limit ordinal and the function f belongs to a certain class of functions containing polynomials with natural number coefficients. The tools developed for this result can be extended to cover all infinite α, but the case of finite α appears to be quite problematic.     

Method of lines with boundary elements for 1‐D transient diffusion‐reaction problems

Time‐dependent differential equations can be solved using the concept of method of lines (MOL) together with the boundary element (BE) representation for the spatial linear part of the equation. The BE method alleviates the need for spatial discretization and casts the problem in an integral format. Hence errors associated with the numerical approximation of the spatial derivatives are totally eliminated. An element level local cubic approximation is used for the variable at each time step to facilitate the time marching and the nonlinear terms are represented in a semi‐implicit manner by a local linearization at each time step. The accuracy of the method has been illustrated on a number of test problems of engineering significance. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Cochran's statistical theorem for outer inverses of matrices and matrix quadratic forms

We extend the matrix version of Cochran's statistical theorem to outer inverses of a matrix. As applications, we investigate the Wishartness and independence of matrix quadratic forms for Kronecker product covariance structures.     

Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

A WKB analysis of the buckling condition for a cylindrical shell of arbitrary thickness subjected to an external pressure

We consider the plane-strain buckling of a cylindrical shellof arbitrary thickness which is made of a Varga material andis subjected to an external hydrostatic pressure on its outersurface. The WKB method is used to solve the eigenvalue problemthat results from the linear bifurcation analysis. We show thatthe circular cross-section buckles into a non-circular shapeat a value of µ₁ which depends on A₁/A₂ and a mode number,where A₁ and A₂ are the undeformed inner and outer radii, andµ₁ is the ratio of the deformed inner radius to A₁. Inthe large mode number limit, we find that the dependence ofµ₁ on A₁/A₂ has a boundary layer structure: it is constantover almost the entire region of 0 < A₁/A₂ < 1 and decreasessharply from this constant value to unity as A₁/A₂ tends tounity. Our asymptotic results for A₁ – 1 = O(1) and A₁– 1 = O(1/n) are shown to agree with the numerical resultsobtained by using the compound matrix method.     

Absorbing boundary conditions for reaction-diffusion equations

We treat here of the question of absorbing boundary conditionsfor nonlinear diffusion equations. We use the conditions designedfor the linear equation, we prove them to be well posed forthe nonlinear problem, and through numerical experiments thatthey are well suited for reaction–diffusion equations.     

A posteriori error estimates for a time discrete scheme for a phase-field system of Penrose-Fife type

A time discrete scheme is used to approximate the solution toa phase field system of Penrose–Fife type with a non-conservedorder parameter. An a posteriori error estimate is presentedthat allows the estimation of the difference between continuousand semidiscrete solutions by quantities that can be calculatedfrom the approximation and given data.     

Sulfonated poly(aryl ether ketone)s containing naphthalene moieties obtained by direct copolymerization as novel polymers for proton exchange membranes

A series of sulfonated poly(aryl ether ketone)s (SPAEKs) were prepared by aromatic nucleophilic polycondensation of 2,6‐dihydroxynaphthalene with 5,5′‐carbonyl‐bis(2‐fluorobenzenesulfonate) and 4,4′‐difluorobenzophenone. The structure and degree of sulfonation (DS) of the SPAEKs were characterized using ¹H NMR spectroscopy. The experimentally observed DS values were close to the expected values derived from the starting material ratios. The thermal stabilities of the SPAEKs were characterized by thermogravimetric analysis, which showed that in acid and sodium salt forms they were thermally stable in air up to about 240 and 380 °C, respectively. Transparent membranes cast from the directly polymerized SPAEKs exhibited good mechanical properties in both dry and hydrated states. The dependence of water uptake and of membrane swelling on the DS at different temperatures was studied. SPAEK membranes with a DS from 0.72 to 1.60 maintained adequate mechanical properties after immersion in water at 80 °C for 24 h. The proton conductivity of SPAEK membranes with different degrees of sulfonation was measured as a function of temperature. The proton conductivity of the SPAEK films increased with increased DS, and the highest room temperature conductivity (4.2 × 10^?2 S/cm) was recorded for a SPAEK membrane with a DS of 1.60, which further increased to 1.1 × 10^?1 S/cm at 80 °C. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 2866–2876, 2004     

UV‐ignited frontal polymerization of an epoxy resin

By combining frontal polymerization and radical‐induced cationic polymerization, it was possible to cure thick samples of an epoxy monomer bleached by UV light. The effect of the relative amounts of cationic photoinitiator and radical initiator was thoroughly investigated and was related to the front's velocity and its maximum temperature. The materials obtained were characterized by quantitative conversion also in the deeper layers, not reached by UV light. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 2066–2072, 2004