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 analysis of the spiral Couette flow of a rarefied gas between coaxial cylinders

An implicit quasi-monotone second-order accurate method is proposed for analyzing the spiral Couette flow of a rarefied gas between coaxial cylinders. The basic advantages of the method over the conventional method of stationry iterations are that the former is conservative with respect to the collision integral, has a simple software implementation for any types of boundary conditions, and applies to a wide range of Knudsen numbers.     

LocalH-theorem for a kinetic variational theory

A localH-theorem is derived for a recently proposed extension of Enskog kinetic theory to a dense model fluid composed of particles with interactions extending beyond a hard core.On leave from: Katedra Fizyki, Uniwersytetu Szczecinskiego, 70-451 Szczecin, Poland.     

铁铝催化体系催化丙烯酸丁酯聚合

研究了Fe(acac)3-Al(i-Bu)3-8-羟基喹啉(acac=乙酰丙酮)催化体系催化丙烯酸丁酯(BA)聚合，考察了聚合规律，用凝胶渗透色谱研究了聚合物分子量和分子量分布．动力学研究表明聚合反应对单体浓度呈一级关系，表观活化能为13．9kJ／mol．     

可压稠密气固两相流动过程的模拟计算

基于稠密气体分子运动论和颗粒动理学,建立可压稠密气固两相流动模型。采用梯度模拟来考虑气相可压缩性对气相湍流的影响。模拟计算表明气固两相射流速度沿轴向和径向减小,颗粒浓度下降。气固两相射流具有高的颗粒温度,呈现强烈的气固两相湍流流动特性。     

Direct numerical simulation of settling of a large solid particle during bioconvection

Settling of a large solid particle in bioconvection flow caused by gyrotactic microorganisms is investigated. The particle is released from the top of the bioconvection chamber; its settling pattern depends on whether it is released in the centre of the bioconvection plume or at its periphery. The Chimera method is utilized; a subgrid is generated around a moving particle. The method suggested by Liu and Wang (Comput. Fluid 2004; 33 :223–255) is further developed to account for the presence of a moving boundary in the streamfunction‐vorticity formulation using the finite‐difference method. A number of cases for different release positions of the particle are computed. It is demonstrated that bioconvection can either accelerate or decelerate settling of the particle depending on the initial position of the particle relative to the plume centre. It is also shown that the particle impacts bioconvection plume by changing its shape and location in the chamber. Copyright © 2005 John Wiley & Sons, Ltd.     

A numerical and experimental investigation of the interactions between a non‐uniform planar array of incompressible free jets

The interaction between multiple incompressible air jets has been studied numerically and experimentally. The numerical predictions have been first validated using experimental data for a single jet configuration. The spreading features of five unequal jets in the configuration of one larger central jet surrounded by four smaller equi‐distant jets, have been studied, for different lateral spacing ratios of 1.5, 2.0 and 2.5 and a central jet Reynolds number of 1.24×105 (corresponding to a Mach number of 0.16). Flow of five equal jets has also been simulated, for the sake of comparison. The jet interactions commence at an axial distance of about 3–4 diameters and complete by an axial distance of about 10 diameters for the lowest spacing ratio of 1.5. For larger spacing ratios, the length required for the start and completion of jet interaction increase. Peripheral jets bend more towards the central jet and merge at a smaller distance, when their sizes are smaller than that of the central jet. The entrainment ratio for multiple jets is higher than that for a single jet. Excellent agreement is observed between the experimental data and theoretical predictions for both mean flow field and turbulent quantities, at regions away from the jet inlet. The potential core length and initial jet development, however, are not predicted very accurately due to differences in the assumed and actual velocity profiles at the jet inlet. Copyright © 2004 John Wiley & Sons, Ltd.     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

Dependence of the kinetics of the anionic polymerization of methyl methacrylate on the concentration in active centers

Studies on the anionic polymerization of methyl methacrylate in tetrahydrofuran and in the presence of sparteine have revealed a beneficial effect due to this additive, resulting in a decrease in the extent of termination. Better control of the definition of the polymers formed can thus be achieved in the presence of this additive. On the other hand, macromolecular engineering requires a range of active species concentrations lower than 10^?3 mol L^?1 and particularly the synthesis of polymers of high molar masses. For a better understanding of the mechanism of chain growth under such concentration conditions, the kinetics of polymerization have been investigated with a technique based on adiabatic calorimetry. Sparteine has been found to lack sufficient cation‐binding power to prevent the propagating enolate ion pairs from aggregating. The rate constant of propagation of nonaggregated species has been estimated, as well as the aggregation constant of equilibrium. For very low initiator concentrations, termination reactions have been shown to profoundly alter the control of the polymerization and to prevent a quantitative monomer conversion. Theoretical maximal conversions have been calculated from kinetic data and compare well with the experimental values. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4964–4975, 2004     

Heterochain model for simultaneous long‐chain branching and crosslinking. III. Multicomponent polymerization

The matrix formula developed in the context of heterochain theory, M?_w = M?_wp + WF ( I ? M )^?1 S , was applied to describe the molecular weight development during free‐radical multicomponent polymerization. All of the required probabilistic parameters are expressed in terms of the kinetic‐rate constants and the various concentrations associated with them. In free‐radical polymerization, the number of heterochain types, N, needs to be extrapolated to infinity, and such extrapolation is conducted with only three different N values. This matrix formula can be used as a benchmark test if other approximate approaches can give reasonable estimates of the weight‐average molecular weights. The moment equations with the average pseudo‐kinetic‐rate constants for branching and crosslinking reactions may provide poor estimates when the copolymer composition drift during polymerization is very significant. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 2801–2812, 2004