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.     

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 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     

乌桕类可可脂结晶过程中表现体积变化

为了解乌桕类可可脂(CTCBE)缓冷结晶横向胀罐爆裂原因，研究CTCBE结晶过程中的膨胀特性，采用流体静力法测试了不同等温结晶条件下CTCBE结晶形成的表观密度和表观体积及其横纵向膨胀状况。结果表明，CTCBE缓冷结晶后形成外部和中部两个晶区。在5-25℃的各等温结晶范围内，当结晶温度增加时，外部结晶区域减小、表观密度变化不大；中部结晶区域增大、表观密度明显减小，CTCBE的表观体积、横向、纵向膨胀率均增大，膨胀横向大于纵向。进一步表明CTCBE在自然缓冷固化的表观体积根本上由其结晶温度决定。为有效防止横向胀罐，自然缓冷固化温度应低于10℃。     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

Statistical thermodynamics evaluation of polymer–polymer miscibility

The Simha and Somcynsky (S–S) statistical thermodynamics theory was used to compute the solubility parameters as a function of temperature and pressure [δ = δ(T, P)], for a series of polymer melts. The characteristic scaling parameters required for this task, P*, T*, and V*, were extracted from the pressure–temperature–volume (PVT) data. To determine the potential polymer–polymer miscibility, the dependence of δ versus T (at ambient pressure) was computed for 17 polymers. Close proximity of the δ versus T curves for four miscible polymer pairs: PPE/PS, PS/PVME, and PC/PMMA signaled the usefulness of this approach. It is noteworthy, that the tabulated solubility parameters (derived from the solution data under ambient conditions) propounded the immiscibility of the PVC/PVAc pair. The computed values of δ also suggested miscibility for polymer pairs of unknown miscibility, namely PPE/PVC, PPE/PVAc, and PET/PSF. In recognizing the limitations of the solubility parameter approach (the omission of several thermodynamic contributions), these preliminary results are auspicious because they indicate a new route for estimating the miscibility of any polymeric material at a given temperature and pressure. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 2909–2915, 2004     

The Modeling of Velocity Enhancement in Polymer Flooding

In single-phase polymer flooding experiments it has repeatedly been observed that the average velocity of the polymer molecules is higher than the average velocity of the water molecules. This effect is incorporated in many conventional Enhanced Oil Recovery (EOR) simulators by the introduction of a constant velocity enhancement factor. In this paper we show that, in absence of dispersion, a constant enhancement factor in the mathematical model for two-phase polymer flow (Buckley--Leverett displacement) leads to ill-posedness of the model equations. We propose a saturation dependent enhancement factor, derived from a model based on percolation concepts, for which this problem does not occur.     

The effect of gamma irradiation and shelf aging in air on the oxidation of ultra-high molecular weight polyethylene

This study has investigated the effect of shelf aging, for up to one year in air, on the properties of gamma-irradiated ultra-high molecular weight polyethylene (UHMWPE). A variety of techniques were used to characterize the properties of treated samples. Differential scanning calorimetery (DSC) was used to characterize the morphology. The extent of cross-linking in a polymer network was detected by swelling measurements. The durometer hardness test was used to measure the relative hardness of this material, and changes in density were also measured. Results from all these measurements were combined to explain the changes in the microstructure of the aged, irradiated UHMWPE. This study shows that crystallinity is increased with radiation dose and with aging due to chain scission, which leads to a reduction in the molecular weight of the material. This allows the chains to rearrange to form crystalline regions. Positron annihilation lifetime spectroscopy confirms these conclusions. Fractional free volumes have been deduced from lifetime parameters, which correlate with the data obtained by the other techniques.