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.     

Structure of solid monolayers and multilayers of n -hexane on graphite

We present all-atom molecular dynamics simulations ofn-hexane on the basal plane of graphite at monolayer and multilayer coverages. In keeping with experimental data, we find the presence of ordered adsorbed layers both at single monolayer coverage and when the adsorbed layer coexists with excess liquid adsorbate. Using a simulation method that does not impose any particular periodicity on the adsorbed layer, we quantitatively compare our results to the results of neutron diffraction experiments and find a structural transition from a uniaxially incommensurate lattice to a fully commensurate structure on increasing the coverage from a monolayer to a multilayer. The zig-zag backbone planes of all the alkane molecules lie parallel to the graphite surface at the multilayer coverage, while a few molecules are observed to attain the perpendicular orientation at monolayer coverage. Dedicated to Professor C N R Rao on his 70th birthday     

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

为了解乌桕类可可脂(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     

Karhunen–Loeve representations of turbulent channel flows using the method of snapshots

For three‐dimensional flows with one inhomogeneous spatial coordinate and two periodic directions, the Karhunen–Loeve procedure is typically formulated as a spatial eigenvalue problem. This is normally referred to as the direct method (DM). Here we derive an equivalent formulation in which the eigenvalue problem is formulated in the temporal coordinate. It is shown that this so‐called method of snapshots (MOS) has some numerical advantages when compared to the DM. In particular, the MOS can be formulated purely as a matrix composed of scalars, thus avoiding the need to construct a matrix of matrices as in the DM. In addition, the MOS avoids the need for so‐called weight functions, which emerge in the DM as a result of the non‐uniform grid typically employed in the inhomogeneous direction. The avoidance of such weight functions, which may exhibit singular behaviour, guarantees satisfaction of the boundary conditions. The MOS is applied to data sets recently obtained from the direct simulation of turbulence in a channel in which viscoelasticity is imparted to the fluid using a Giesekus model. The analysis reveals a steep drop in the dimensionality of the turbulence as viscoelasticity is increased. This is consistent with the results that have been obtained with other viscoelastic models, thus revealing an essential generic feature of polymer‐induced drag reduced turbulent flows. Published in 2006 by John Wiley & Sons, Ltd.     

镍-酸性铬兰K络合物吸附波

用微分脉冲溶出伏安法研究了镍(11)与酸性铬兰K(ACBK)所生成的络合物在悬汞电极上的电化学特性及电化学反应机理，发现电极反应为镍与酸性铬兰K所生成的络合物吸附在电极表面后发生的不可逆还原。在HCl-硼砂介质中(pH8．67)，络合物在-316mV处有一吸附还原峰，其峰电流与镍浓度在0．001～0．010μ/mL范围内呈现良好的线性关系，其最低检出浓度为0．0005μg／mL。     

Synthesis and adsorption properties of gold nanoparticles within pores of surface‐functional porous polymer microspheres

Mesoporous polymer microspheres with gold (Au) nanoparticles inside their pores were prepared considering their surface functionality and porosity. The Au/polymer composite microspheres prepared were characterized by transmission electron microscope (TEM), X‐ray diffraction (XRD), and Brunauer–Emmett–Teller (BET) techniques. The results showed that the adsorption of Au nanoparticles could be increased by imparting the pore structure and surface‐functional groups into the supporting polymer microspheres (in this study, poly (ethylene glycol dimethacrylate‐co‐acrylonitrile) and poly (EGDMA‐co‐AN) system). Above all, from this study, it was established that the porosity of the polymer microspheres is the most important factor that determines the distribution and adsorption amount of face‐centered cubic (fcc) Au nanoparticles in the final products. Our study showed that the continuous adsorption of Au nanoparticles with the aid of the large surface area and surface interaction sites formed more favorably the Au/polymer composite microspheres. The BET measurements of Au/poly(EGDMA‐co‐AN) composite microspheres reveals that the adsorption of Au nanoparticles into the pores kept the pore structure intact and made it more porous. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5627–5635, 2004