A virtual boundary method with improved computational efficiency using a multi‐grid method

The flow around spherical, solid objects is considered. The boundary conditions on the solid boundaries have been applied by replacing the boundary with a surface force distribution on the surface, such that the required boundary conditions are satisfied. The velocity on the boundary is determined by extrapolation from the flow field. The source terms are determined iteratively, as part of the solution. They are then averaged and are smoothed out to nearby computational grid points. A multi‐grid scheme has been used to enhance the computational efficiency of the solution of the force equations. The method has been evaluated for flow around both moving and stationary spherical objects at very low and intermediate Reynolds numbers. The results shows a second order accuracy of the method both at creeping flow and at Re=100. The multi‐grid scheme is shown to enhance the convergence rate up to a factor 10 as compared to single grid approach. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling of a single drop impact onto liquid film using particle method

The process of single liquid drop impact on thin liquid surface is numerically simulated with moving particle semi‐implicit method. The mathematical model involves gravity, viscosity and surface tension. The model is validated by the simulation of the experimental cases. It is found that the dynamic processes after impact are sensitive to the liquid pool depth and the initial drop velocity. In the cases that the initial drop velocity is low, the drop will be merged with the liquid pool and no big splash is seen. If the initial drop velocity is high enough, the dynamic process depends on the liquid depth. If the liquid film is very thin, a bowl‐shaped thin crown is formed immediately after the impact. The total crown subsequently expands outward and breaks into many tiny droplets. When the thickness of the liquid film increases, the direction of the liquid crown becomes normal to the surface and the crown propagates outward. It is also found that the radius of the crown is described by a square function of time: r_C = [c(t ? t₀)]^0.5. When the liquid film is thick enough, a crown and a deep cavity inside it are formed shortly after the impact. The bottom of the cavity is initially oblate and then the base grows downward to form a sharp corner and subsequently the corner moves downward. Copyright © 2004 John Wiley & Sons, Ltd.     

Least‐squares meshfree method for incompressible Navier–Stokes problems

A least‐squares meshfree method based on the first‐order velocity–pressure–vorticity formulation for two‐dimensional incompressible Navier–Stokes problem is presented. The convective term is linearized by successive substitution or Newton's method. The discretization of all governing equations is implemented by the least‐squares method. Equal‐order moving least‐squares approximation is employed with Gauss quadrature in the background cells. The boundary conditions are enforced by the penalty method. The matrix‐free element‐by‐element Jacobi preconditioned conjugate method is applied to solve the discretized linear systems. Cavity flow for steady Navier–Stokes problem and the flow over a square obstacle for time‐dependent Navier–Stokes problem are investigated for the presented least‐squares meshfree method. The effects of inaccurate integration on the accuracy of the solution are investigated. 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.     

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     

Effects of poly[styrene‐b‐(ethylene‐co‐butylene)‐b‐styrene] on the charge distributions of low‐density polyethylene/polystyrene blends

The effect of the triblock copolymer poly[styrene‐b‐(ethylene‐co‐butylene)‐b‐styrene] (SEBS) on the formation of the space charge of immiscible low‐density polyethylene (LDPE)/polystyrene (PS) blends was investigated. Blends of 70/30 (wt %) LDPE/PS were prepared through melt blending in an internal mixer at a blend temperature of 220 °C. The amount of charge that accumulated in the 70% LDPE/30% PS blends decreased when the SEBS content increased up to 10 wt %. For compatibilized and uncompatibilized blends, no significant change in the degree of crystallinity of LDPE in the blends was observed, and so the effect of crystallization on the space charge distribution could be excluded. Morphological observations showed that the addition of SEBS resulted in a domain size reduction of the dispersed PS phase and better interfacial adhesion between the LDPE and PS phases. The location of SEBS at a domain interface enabled charges to migrate from one phase to the other via the domain interface and, therefore, resulted in a significant decrease in the amount of space charge for the LDPE/PS blends with SEBS. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 2813–2820, 2004     

NMR investigation of concentrated alumina and silica slurries

The longitudinal relaxation times (T₁) of water in concentrated silica and alumina slurries were measured as a function of solids content. It was shown that the results could be fit very well with a two-phase fast-exchange model between free and surface-bound water. As expected, values of T₁ for bound water were in the order of 20–2000 times lower than that for free water, indicating a higher effective viscosity of the surface-bound water. The strength of the interaction depended on the particular surface, and all of the aluminas examined interacted more strongly with water than the two silicas studied, which themselves differed considerably. The chemical mechanical polishing (CMP) removal rate of tantalum by silica slurries was shown to be directly correlated with the interaction parameters, derived from the NMR relation times rather than with total surface hydroxyl group concentration.     

X波段六腔渡越管振荡器的高频特性研究

 从麦克斯韦方程出发，采用时域有限差分法（FDTD）和离散傅里叶变换（DFFT）相结合的方法，通过数值计算得出了六腔开放式谐振腔中前四个谐振频率和场分布，计算出的谐振频率与实验测量结果基本相同。比较了开放腔和封闭腔谐振频率，验证了TEM波吸收边界条件，并在实际编程计算中得以应用。计算结果为六腔渡越管振荡器的机理研究提供了依据。