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     

Numerical solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

OPTIMAL INVESTMENT IN THE PROTECTION OF A VULNERABLE BIOLOGICAL RESOURCE

It is assumed that the probability of destruction of a biological asset by natural hazards can be reduced through investment in protection. Specifically a model, in which the hazard rate depends on both the age of the asset and the accumulated invested protection capital, is assumed. The protection capital depreciates through time and its effectiveness in reducing the hazard rate is subject to diminishing returns. It is shown how the investment schedule to maximize the expected net present value of the asset can be determined using the methods of deterministic optimal control, with the survival probability regarded as a state variable. The optimal investment pattern involves “bang-bang-singular” control. A numerical scheme for determining jointly the optimal investment policy and the optimal harvest (or replacement) age is outlined and a numerical example involving forest fire protection is given.     

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 novel tridentate nitrogen donor as ligand in copper catalyzed ATRP of methyl methacrylate

A tridentate ligand, BPIEP: 2,6‐bis[1‐(2,6‐diisopropyl phenylimino) ethyl] pyridine, having central pyridine unit and two peripheral imine coordination sites was effectively employed in controlled/“living” radical polymerization of MMA at 90°C in toluene as solvent, Cu^IBr as catalyst, and ethyl‐2‐bromoisobutyrate (EBiB) as initiator resulting in well‐defined polymers with polydispersities M_w/M_n ≤ 1.23. The rate of polymerization follows first‐order kinetics, k_app = 3.4 × 10^?5 s^?1, indicating the presence of low radical concentration ([P*] ≤ 10^?8) throughout the reaction. The polymerization rate attains a maximum at a ligand‐to‐metal ratio of 2:1 in toluene at 90°C. The solvent concentration (v/v, with respect to monomer) has a significant effect on the polymerization kinetics. The polymerization is faster in polar solvents like, diphenylether, and anisole, as compared to toluene. Increasing the monomer concentration in toluene resulted in a better control of polymerization. The molecular weights (M_n,SEC) increased linearly with conversion and were found to be higher than predicted molecular (M_n,Cal). However, the polydispersity remained narrow, i.e., ≤1.23. The initiator efficiency at lower monomer concentration approaches a value of 0.7 in 110 min as compared to 0.5 in 330 min at higher monomer concentration. The aging of the copper salt complexed with BPIEP had a beneficial effect and resulted in polymers with narrow polydispersitities and higher conversion. PMMA obtained at room temperature in toluene (33%, v/v) gave PDI of 1.22 (M_n = 8500) in 48 h whereas, at 50°C the PDI is 1.18 (M_n = 10,300), which is achieved in 23 h. The plot of lnk_app versus 1/T gave an apparent activation energy of polymerization as (ΔE^≠_app) 58.29 KJ/mol and enthalpy of equilibrium (ΔH⁰_eq) to 28.8 KJ/mol. Reverse ATRP of MMA was successfully performed using AIBN in bulk as well as solution. The controlled nature of the polymerization reaction was established through kinetic studies and chain extension experiments. © 2005 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 43: 4996–5008, 2005     

Synthesis and Characterization of Silica@Copper Core–Shell Nanoparticles: Application for Conjugate Addition Reactions

Silica@copper (SiO₂@Cu) core–shell nanoparticles were synthesized and well characterized by XRD, TEM, AFM, XPS, UV/Vis, TGA–MS, and ICP–AES techniques. The synthesized SiO₂@Cu core–shell nanoparticles were employed as catalysts for the conjugate addition of amines to α,β‐unsaturated compounds in water to obtain β‐amino carbonyl compounds in excellent yields in shorter reaction times. Furthermore, the catalyst works well for hetero‐Michael addition reactions of heteroatom nucleophiles such as thiols to α,β‐unsaturated compounds. As the reaction is performed in water, it allows for easy recycling of the catalyst with consistent activity.