Semi‐analytical integration of the 8‐node plane element stiffness matrix using symbolic computation

The semi‐analytical integration of an 8‐node plane strain finite element stiffness matrix is presented in this work. The element is assumed to be super‐parametric, having straight sides. Before carrying out the integration, the integral expressions are classified into several groups, thus avoiding duplication of calculations. Symbolic manipulation and integration is used to obtain the basic formulae to evaluate the stiffness matrix. Then, the resulting expressions are postprocessed, optimized, and simplified in order to reduce the computation time. Maple symbolic‐manipulation software was used to generate the closed expressions and to develop the corresponding Fortran code. Comparisons between semi‐analytical integration and numerical integration were made. It was demonstrated that semi‐analytical integration required less CPU time than conventional numerical integration (using Gaussian‐Legendre quadrature) to obtain the stiffness matrix. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Cochran's statistical theorem for outer inverses of matrices and matrix quadratic forms

We extend the matrix version of Cochran's statistical theorem to outer inverses of a matrix. As applications, we investigate the Wishartness and independence of matrix quadratic forms for Kronecker product covariance structures.     

Matrix effects in the determination of butyltin compounds in environmental samples by GC AA after hydride generation

Experiments for the determination of mono-, di and tri-butyltin (MBT, DBT and TBT) by hydride generation/gas chromatography/atomic absorption spectrometry in various matrices (sediment, suspended matter, mussel, algae and water) have revealed that poor butyltin recoveries are obtained in sediments displaying high sulphur and hydrocarbon contents; very poor recoveries were also observed for TBT in sediments with high chlorophyll pigment contents as well as in algal samples. It was however not clear whether the hydride generatin was inhibited by these infering compounds, as was previously assumed in the case of hydrocarbons, or whether interferences affected the atomization rate. Further studies were performed to solve this problem in order to validate this method in the case of analyses of, for example, oil-contaminated sediment and algae. This paper presents the results obtained. It is concluded here that the poor recoveries were due to an inhibition of hydride generation rather than to interference at the atomization stage.     

Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

A “Multiplicative Coboundary” theorem for some sequences of random matrices

Schmidt(¹⁰) showed that if the family of distributions of the partial sums of a strictly stationary random sequence is tight, then that random sequence is a coboundary. Here an analog of that result is proved for some sequences of random matrices, with partial sums replaced by matrix products.     

Quantitative separation of zinc traces from cadmium matrices by solid-phase extraction with polyurethane foam

A system for separation of zinc traces from large amounts of cadmium is proposed in this paper. It is based on the solid-phase extraction of the zinc in the form of thiocyanate complexes by the polyurethane foam. The following parameters were studied: effect of pH and of the thiocyanate concentration on the zinc extraction, shaking time required for quantitative extraction, amount of PU foam necessary for complete extraction, conditions for the separation of zinc from cadmium, influence of other cations and anions on the zinc sorption by PU foam, and required conditions for back extraction of zinc from the PU foam. The results show that zinc traces can be separated from large amounts of cadmium at pH 3.0±0.50, with the range of thiocyanate concentration from 0.15 to 0.20 mol l⁻¹, and the shaking time of 5 min. The back extraction of zinc can be done by shaking it with water for 10 min. Calcium, barium, strontium, magnesium, aluminum, nickel and iron(II) are efficiently separated. Iron(III), copper(II) and cobalt(II) are extracted simultaneously with zinc, but the iron reduction with ascorbic acid and the use of citrate to mask copper(II) and cobalt(II) increase the selectivity of the zinc extraction. The anions nitrate, chloride, sulfate, acetate, thiosulphate, tartarate, oxalate, fluoride, citrate, and carbonate do not affect the zinc extraction. Phosphate and EDTA must be absent. The method proposed was applied to determine zinc in cadmium salts using 4-(2-pyridylazo)-resorcinol (PAR) as a spectrophotometric reagent. The result achieved did not show significant difference in the accuracy and precision (95% confidence level) with those obtained by ICP–AES analysis.     

The computation of the representation matrices of the generators of the unitary group

A method is presented for the efficient computation of the representation matrices of the unitary group, U(n) in the Gelfand—Tsetlin basis (corresponding to the usual spin-symmetry adapted basis for an N electron CI). The present scheme is conceptually and computationally attractive in that it is formulated directly in terms of Weyl tableaux and also that it permits simultaneous basis vector generation and matrix element evaluation. In addition the basis vectors are ordered so that subsequent restriction to the three dimensional rotation group is facilitated. An illustrative example is also presented.Taken in part from a thesis submitted to the University of London in partial fulfilment of the requirements for the degree of PhD.     

Catalysis in Organic Solvents with Lipase Immobilized by Sol-Gel Technique

A series of immobilized lipases were obtained by sol-gel process, using silica prepolymers prepared from tetramethoxysilane, methyltrimethoxysilane, propyltrimethoxysilane and (3-aminopropyl)triethoxysilane. The activities of these biocatalysts were compared with the lipase adsorbed on poly(methylhydroxysiloxane) and encapsulated into a silicone rubber, lipase entrapped in nanoporous silica matrix and commercial sol-gel lipase. Model reactions were the esterification of stearic acid and Corey lactone bisalcohol (an intermediate of prostaglandin synthesis). The positive effect of hydrophobic-hydrophilic interface, created by the addition of organosilanes, on the activity of biocatalysts was partially reduced by decreasing specific surface of mesopores. Hydrophobic solvents increased the activity of the lipase entrapped in tetramethoxysilane–methyltrimethoxysilane prepolymer in the sequence acetone < toluene < benzene < decane < hexane. The activity of silicone rubber-encapsulated biocatalysts was proportional to polymer swelling in organic solvents (hexane > toluene > acetone).     

STRUCTURES OF CIRCULANT INVERSE M-MATRICES

In this paper, we present a useful result on the structures of circulant inverse Mmatrices. It is shown that if the n × n nonnegative circulant matrix A = Circ[c0, c1,… , c（n- 1）] is not a positive matrix and not equal to c0I, then A is an inverse M-matrix if and only if there exists a positive integer k, which is a proper factor of n, such that cjk 〉 0 for j=0,1…, [n-k/k], the other ci are zero and Circ[co, ck,… , c（n-k）] is an inverse M-matrix. The result is then extended to the so-called generalized circulant inverse M-matrices.     

一种线性模型的精确转换方法

基于广义范德蒙矩阵的逆矩阵的计算方法,导出线性系统的连续时间模型和离散时间模型在含有重极点时的严格的直接换算关系,其结果可以直接表示为两者的分子多项式系数之间的关系,并且,对不同的调制器具有统一的转换形式.