An analysis of over-relaxation in a kinetic approximation of systems of conservation laws

The over-relaxation approach is an alternative to the Jin–Xin relaxation method in order to apply the equilibrium source term in a more precise way. This is also a key ingredient of the lattice Boltzmann method for achieving second-order accuracy. In this work, we provide an analysis of the over-relaxation kinetic scheme. We compute its equivalent equation, which is particularly useful for devising stable boundary conditions for the hidden kinetic variables.     

Fused vs. spiro: Kinetic,not thermodynamic preference may direct the reaction of α-carbonyl oxonium ylides

Rh(II)-catalyzed decomposition of certain cyclic α-diazocarbonyl compounds in the presence of cyclic ethers has been shown to give bicyclic ring expansion products. These are thought to arise from a [1,4]-alkyl shift toward the carbonyl oxygen atom and are in contrast with the recently observed spirocyclic products of a Stevens-type [1,2]-alkyl shift within the postulated oxonium ylide intermediate. Quantum chemical calculations performed at the B3LYP/6-31G* level of theory showed that the former reaction pathway (toward fused bicycles) is kinetically preferred.     

The mechanism of ethylene polymerization by nickel salicylaldiminato catalysts - Agostic interactions and their kinetic isotope effects

The ethylene polymerization reaction of a neutral nickel catalyst was studied by DFT calculations at the Becke3LYP/6-31G(d) level of theory. As in related cases a β-agostic bond stabilizes the nickel alkyl ground states. Transition states for the insertion of the olefin show a distinct α-agostic interaction, which has not been observed for late metal polymerization catalysts before. An ethylene-alkyl complex was identified as the resting state of the reaction. The overall barrier height of the reaction amounts to 17.54 kcal/mol, which slightly increases to 17.60 kcal/mol for the polymerization of deuterated ethylene. Therefore, a small positive kinetic isotope effect (k_H/k_D = 1.09) can be calculated, which is caused by the α-agostic interaction in the transition state. A comparison to other late metal based polymerization systems reveals that the ethylene coordination step of highly active catalysts is significantly lower in energy compared to catalysts which are only moderately active.     

催化动力学光度法测定水中痕量铁

本文研究了ｐｈ５．０的ＨＡｃ－ＮａＡｃ介质中，痕量铁催化过氧化氢氧化紫脲酸的褪色反应，建立了催化动力学光度法测定痕量铁的新方法，方法的线性范围０－０．４０，０．４０－２．０μｇ／２５ｍＬ，检出限为８×１０＾－１１ｇ／ｍＬ，用于水样中痕量铁的测定，结果满意。     

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     

Finite Element Solution of Conical Diffraction Problems

This paper is devoted to the numerical study of diffraction by periodic structures of plane waves under oblique incidence. For this situation Maxwell's equations can be reduced to a system of two Helmholtz equations in R ² coupled via quasiperiodic transmission conditions on the piecewise smooth interfaces between different materials. The numerical analysis is based on a strongly elliptic variational formulation of the differential problem in a bounded periodic cell involving nonlocal boundary operators. We obtain existence and uniqueness results for discrete solutions and provide the corresponding error analysis.     

The Linear Fokker-Planck Equation for the Ornstein-Uhlenbeck Process as an (Almost) Nonlinear Kinetic Equation for an Isolated N-Particle System

It is long known that the Fokker-Planck equation with prescribed constant coefficients of diffusion and linear friction describes the ensemble average of the stochastic evolutions in velocity space of a Brownian test particle immersed in a heat bath of fixed temperature. Apparently, it is not so well known that the same partial differential equation, but now with constant coefficients which are functionals of the solution itself rather than being prescribed, describes the kinetic evolution (in the N→∞ limit) of an isolated N-particle system with certain stochastic interactions. Here we discuss in detail this recently discovered interpretation. An erratum to this article can be found at     

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 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.     

A complete boundary integral formulation for compressible Navier–Stokes equations

A complete boundary integral formulation for compressible Navier–Stokes equations with time discretization by operator splitting is developed using the fundamental solutions of the Helmholtz operator equation with different order. The numerical results for wall pressure and wall skin friction of two‐dimensional compressible laminar viscous flow around airfoils are in good agreement with field numerical methods. Copyright © 2004 John Wiley & Sons, Ltd.