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 solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

毛细管气相色谱法测定中药中有机氯农药残留量

建立了江西省两种中药中有机氯农药残留量的气相色谱分析方法。样品以丙酮提取,在NaCl存在下,以正已烷进行液—液分配,提取液用H2SO4净化,采用DB—1701毛细管柱,GC—ECD检测有机氯农药残留量。最低检测限为1.0×10-3-5.0×10-3mg/kg,添加回收率为85.7%-115.2%,应用于实际样品中有机氯农药残留量的检测,获得了较为满意结果。     

Rhythmic growth of ring‐banded spherulites in blends of liquid crystalline methoxy‐poly(aryl ether ketone) and poly(aryl ether ether ketone)

Rhythmic growth of ring‐banded spherulites in blends of liquid crystalline methoxy‐poly(aryl ether ketone) (M‐PAEK) and poly(aryl ether ether ketone) (PEEK) has been investigated by means of differential scanning calorimetry (DSC), polarized light microscopy (PLM), and scanning electron microscopy (SEM) techniques. The measurements reveal that the formation of the rhythmically grown ring‐banded spherulites in the M‐PAEK/PEEK blends is strongly dependent on the blend composition. In the M‐PAEK‐rich blends, upon cooling, an unusual ring‐banded spherulite is formed, which is ascribed to structural discontinuity caused by a rhythmic radial growth. For the 50:50 M‐PAEK/PEEK blend, ring‐banded spherulites and individual PEEK spherulites coexist in the system. In the blends with PEEK as the predominant component, M‐PAEK is rejected into the boundary of PEEK spherulites. The cooling rate and crystallization temperature have great effect on the phase behavior, especially the ring‐banded spherulite formation in the blends. In addition, the effects of M‐PAEK phase transition rate and phase separation rate on banded spherulite formation is discussed. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 3011–3024, 2007     

Living/controlled radical copolymerization of N‐substituted maleimides with styrene in 1‐butyl‐3‐methylimidazolium hexafluorophosphate and anisole

The atom transfer radical copolymerization of N‐substituted maleimides such as N‐phenylmaleimide (PhMI), N‐cyclohexylmaleimide (ChMI), and N‐butylmaleimide (NBMI) with styrene initiated with dendritic polyarylether 2‐bromoisobutyrates in an ionic liquid, 1‐butyl‐3‐methylimidazolium hexafluorophosphate ([bmim][PF₆]), at room temperature and anisole at 110 °C was investigated. The dendritic‐linear block copolymers obtained in ionic liquid possessed well‐defined molecular weight and low polydispersity (1.05 < M_w/M_n < 1.32) and could be used as a macroinitiator for chain‐extension polymerization, suggesting the living nature of the reaction system. The ionic liquids containing catalyst could be recycled in the atom transfer radical polymerization systems without further treatment. Compared with polymerization conducted in anisole, the polymerization in ionic liquid had a stronger tendency for alternation. The tendency for alternation decreased in the order PhMI > NBMI > ChMI in [bmim][PF₆] and PhMI > ChMI > NBMI in anisole. © 2003 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 41: 2156–2165, 2003     

A test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.     

Computations of two passing‐by high‐speed trains by a relaxation overset‐grid algorithm

This paper presents a relaxation algorithm, which is based on the overset grid technology, an unsteady three‐dimensional Navier–Stokes flow solver, and an inner‐ and outer‐relaxation method, for simulation of the unsteady flows of moving high‐speed trains. The flow solutions on the overlapped grids can be accurately updated by introducing a grid tracking technique and the inner‐ and outer‐relaxation method. To evaluate the capability and solution accuracy of the present algorithm, the computational static pressure distribution of a single stationary TGV high‐speed train inside a long tunnel is investigated numerically, and is compared with the experimental data from low‐speed wind tunnel test. Further, the unsteady flows of two TGV high‐speed trains passing by each other inside a long tunnel and at the tunnel entrance are simulated. A series of time histories of pressure distributions and aerodynamic loads acting on the train and tunnel surfaces are depicted for detailed discussions. Copyright © 2004 John Wiley & Sons, Ltd.     

Assessment of two‐equation turbulence modelling for high Reynolds number hydrofoil flows

This paper presents an evaluation of the capability of turbulence models available in the commercial CFD code FLUENT 6.0 for their application to hydrofoil turbulent boundary layer separation flow at high Reynolds numbers. Four widely applied two‐equation RANS turbulence models were assessed through comparison with experimental data at Reynolds numbers of 8.284×106 and 1.657×107. They were the standard k–εmodel, the realizable k–εmodel, the standard k–ωmodel and the shear‐stress‐transport (SST) k–ωmodel. It has found that the realizable k–εturbulence model used with enhanced wall functions and near‐wall modelling techniques, consistently provides superior performance in predicting the flow characteristics around the hydrofoil. Copyright © 2004 John Wiley & Sons, Ltd.     

High‐order accurate numerical solutions of incompressible flows with the artificial compressibility method

A high‐order accurate, finite‐difference method for the numerical solution of incompressible flows is presented. This method is based on the artificial compressibility formulation of the incompressible Navier–Stokes equations. Fourth‐ or sixth‐order accurate discretizations of the metric terms and the convective fluxes are obtained using compact, centred schemes. The viscous terms are also discretized using fourth‐order accurate, centred finite differences. Implicit time marching is performed for both steady‐state and time‐accurate numerical solutions. High‐order, spectral‐type, low‐pass, compact filters are used to regularize the numerical solution and remove spurious modes arising from unresolved scales, non‐linearities, and inaccuracies in the application of boundary conditions. The accuracy and efficiency of the proposed method is demonstrated for test problems. 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.