New ligands for the Fe(III)‐mediated reverse atom transfer radical polymerization of methyl methacrylate

A series of (di)picolinic acids and their derivates are investigated as novel complexing tridentate or bidentate ligands in the iron‐mediated reverse atom transfer radical polymerization of methyl methacrylate in N,N‐dimethylformamide at 100 °C with 2,2′‐azobisisobutyrontrile as an initiator. The polymerization rates and polydispersity indices (1.32–1.8) of the resulting polymers are dependent on the structures of the ligands employed. Different iron complexes may be involved in iron‐mediated reverse atom transfer radical polymerization, depending on the type of acid used. ¹H NMR spectroscopy has been used to study the structure of the resulting polymers. Chain‐extension reactions have been performed to further confirm the living nature of this catalytic system. © 2006 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 44: 2912–2921, 2006     

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.     

Multilevel preconditioners for a quadrature Galerkin solution of a biharmonic problem

Efficient multilevel preconditioners are developed and analyzed for the quadrature finite element Galerkin approximation of the biharmonic Dirichlet problem. The quadrature scheme is formulated using the Bogner–Fox–Schmit rectangular element and the product two‐point Gaussian quadrature. The proposed additive and multiplicative preconditioners are uniformly spectrally equivalent to the operator of the quadrature scheme. The preconditioners are implemented by optimal algorithms, and they are used to accelerate convergence of the preconditioned conjugate gradient method. Numerical results are presented demonstrating efficiency of the preconditioners. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Efficient computation of natural convection in a concentric annulus between an outer square cylinder and an inner circular cylinder

In this work, the natural convection in a concentric annulus between a cold outer square cylinder and a heated inner circular cylinder is simulated using the differential quadrature (DQ) method. The vorticity‐stream function formulation is used as the governing equation, and the coordinate transformation technique is introduced in the DQ computation. It is shown in this paper that the outer square boundary can be approximated by a super elliptic function. As a result, the coordinate transformation from the physical domain to the computational domain is set up by an analytical expression, and all the geometrical parameters can be computed exactly. Numerical results for Rayleigh numbers range from 104 to 106 and aspect ratios between 1.67 and 5.0 are presented, which are in a good agreement with available data in the literature. It is found that both the aspect ratio and the Rayleigh number are critical to the patterns of flow and thermal fields. The present study suggests that a critical aspect ratio may exist at high Rayleigh number to distinguish the flow and thermal patterns. Copyright © 2002 John Wiley & Sons, Ltd.     

Acoustic reflection from a baffled elastic/piezoelectric plate

A finite difference/boundary integral procedure to determine the acoustic reflected pressure from a fluid-loaded bi-laminated plate is described. The bi-laminate is composed of a piezoelectric layer and an elastic layer in contact with the fluid, and is held by an acoustically hard baffle. In the numerical model, the fluid pressure at fluid/solid interface is replaced by a continuum of point sources weighted by the normal acceleration of the elastic plate, and the governing equation system is solved in the solid domain. With the normal acceleration found, the reflected pressure in the fluid is determined by an integral expression involving the Green's function. It is demonstrated that an appropriate applied voltage potential across the piezoelectric layer has the effect of cancelling either the reflected or scattered pressure of the plate at any chosen field points in the fluid. Project supported by the National Natural Science Foundation of China (No. 10172039).     

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.     

对求解非线性方程方法的探索

BaiCL在2001年提出了双曲函数法，并把解设成具体的双曲函数，即sinhw和coshw的线性组合．本文把解设成是由2个函数f(ξ)和g(ξ)组成的线性组合，它们可以有多种取法，从而使方程更直接、更有效．借助数学软件Maple，用改进后的方法和吴文俊消元法，求解Bai C L文中的一个例子，获得了包含Bai文结果在内的更为丰富、精确的行波解．     

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