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.     

Laguerre expansion on the Heisenberg group and Fourier-Bessel transform on C~n Dedicated to Professor Sheng GONG on the occasion of his 75th birthday

Given a principal value convolution on the Heisenberg group Hn = Cn×R, we study the relation between its Laguerre expansion and the Fourier-Bessel expansion of its limit on Cn. We also calculate the Dirichlet kernel for the Laguerre expansion on the group Hn.     

High energy asymptotics of the scattering amplitude for the Schrödinger equation

We find an explicit function approximating at high energies the kernel of the scattering matrix with arbitrary accuracy. Moreover, the same function gives all diagonal singularities of the kernel of the scattering matrix in the angular variables. This paper is dedicated to Jean-Michel Combes on the occasion of his sixtieth birthday.     

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.     

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

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

Stress concentration and effective stiffness of aligned fiber reinforced composite with anisotropic constituents

The paper addresses the problem of calculation of the local stress field and effective elastic properties of a unidirectional fiber reinforced composite with anisotropic constituents. For this aim, the representative unit cell approach has been utilized. The micro geometry of the composite is modeled by a periodic structure with a unit cell containing multiple circular fibers. The number of fibers is sufficient to account for the micro structure statistics of composite. A new method based on the multipole expansion technique is developed to obtain the exact series solution for the micro stress field. The method combines the principle of superposition, technique of complex potentials and some new results in the theory of special functions. A proper choice of potentials and new results for their series expansions allow one to reduce the boundary-value problem for the multiple-connected domain to an ordinary, well-posed set of linear algebraic equations. This reduction provides high numerical efficiency of the developed method. Exact expressions for the components of the effective stiffness tensor have been obtained by analytical averaging of the strain and stress fields.     

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.