Approximate Inverse Preconditioners for Some Large Dense Random Electrostatic Interaction Matrices

A sparse mesh-neighbour based approximate inverse preconditioner is proposed for a type of dense matrices whose entries come from the evaluation of a slowly decaying free space Green’s function at randomly placed points in a unit cell. By approximating distant potential fields originating at closely spaced sources in a certain way, the preconditioner is given properties similar to, or better than, those of a standard least squares approximate inverse preconditioner while its setup cost is only that of a diagonal block approximate inverse preconditioner. Numerical experiments on iterative solutions of linear systems with up to four million unknowns illustrate how the new preconditioner drastically outperforms standard approximate inverse preconditioners of otherwise similar construction, and especially so when the preconditioners are very sparse. AMS subject classification (2000) 65F10, 65R20, 65F35, 78A30     

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     

Acoustic reflection from a baffled elastic/piezoelectric plate

A finite difference/boundary integral procedure to determine the acoustic reflected pressure from a fluid-loaded bi-laminated plate is described. The bi-laminate is composed of a piezoelectric layer and an elastic layer in contact with the fluid, and is held by an acoustically hard baffle. In the numerical model, the fluid pressure at fluid/solid interface is replaced by a continuum of point sources weighted by the normal acceleration of the elastic plate, and the governing equation system is solved in the solid domain. With the normal acceleration found, the reflected pressure in the fluid is determined by an integral expression involving the Green's function. It is demonstrated that an appropriate applied voltage potential across the piezoelectric layer has the effect of cancelling either the reflected or scattered pressure of the plate at any chosen field points in the fluid. Project supported by the National Natural Science Foundation of China (No. 10172039).     

Numerical solution of the three‐dimensional parabolic equation with an integral condition

Developement of numerical methods for obtaining approximate solutions to the three dimensional diffusion equation with an integral condition will be carried out. The numerical techniques discussed are based on the fully explicit (1,7) finite difference technique and the fully implicit (7,1) finite difference method and the (7,7) Crank‐Nicolson type finite difference formula. The new developed methods are tested on a problem. Truncation error analysis and numerical examples are used to illustrate the accuracy of the new algorithms. The results of numerical testing show that the numerical methods based on the finite difference techniques discussed in the present article produce good results. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 193–202, 2002; DOI 10.1002/num.1040     

Some degenerate and quasilinear parabolic systems not in divergence form

This paper deals with positive solutions of degenerate and quasilinear parabolic systems not in divergence form: u_t=u^p(Δu+av), v_t=v^q(Δv+bu), with null Dirichlet boundary conditions and positive initial conditions, where p, q, a and b are all positive constants. The local existence and uniqueness of classical solution are proved. Moreover, it will be proved that all solutions exist globally if and only if ab?λ₁², where λ₁ is the first eigenvalue of −Δ in Ω with homogeneous Dirichlet boundary condition.     

Correlation between the Mechanical Properties and the Amount of Desorbed Water for Composites Based on Low-Density Polyethylene and Linen Yarn Production Waste

The effect of the amount of desorbed water on the mechanical properties of composites based on low-density polyethylene and linen yarn production waste (LW) is analyzed by statistical methods. It is shown that the amount of absorbed water decreases during the desorption process at room temperature both for specimens modified and unmodified with diphenylmethane diisocyanate (DIC.) The most sensitive to the action of water is the elastic modulus, which decreases considerably under the effect of water and is fully restored in the desorption process. The tensile strain also increases with the amount of absorbed water. It is found that the elastic modulus of the unmodified composite correlates linearly with the amount of desorbed water. Between the amount of desorbed water and the tensile strain, as well the specific work of deformation, a negative linear correlation is revealed. After water desorption, all strength and deformation characteristics of both the modified and unmodified composites are fully recovered.__________Translated from Mekhanika Kompozitnykh Materialov, Vol. 41, No. 4, pp. 515–524, July–August, 2005.     

Convergence analysis of a modified inexact implicit method for general mixed monotone variational inequalities

We consider a useful modification of the inexact implicit method with a variable parameter in Wang et al. J Optim Theory 111: 431–443 (2001) for generalized mixed monotone variational inequalities. One of the contributions of the proposed method in this paper is that the restrictions imposed on the variable parameter are weaker than the ones in Wang et al. J Optim Theory 111: 431–443 (2001). Another contribution is that we establish a sufficient and necessary condition for the convergence of the proposed method to a solution of the general mixed monotone variational inequality.     

Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

An AFM, XPS and electrochemical study of molecular electroactive monolayers formed by wet chemistry functionalization of H-terminated Si(1 0 0) with vinylferrocene

Molecular electroactive monolayers have been produced from vinylferrocene (VFC) via light-assisted surface anchoring to H-terminated n- and p-Si(1 0 0) wafers prepared via wet chemistry, in a controlled atmosphere. The resulting Si-C bound hybrids have been characterized by means of XPS and AFM. Their performance as semiconductor functionalized electrodes and their surface composition have been followed by combining electrochemical and XPS measurements on the same samples, before and after use in an electrochemical cell. White-light photoactivated anchoring at short (1 h) exposure times has resulted in a mild route, with a very limited impact on the initial quality of the silicon substrate. In fact, the functionalized Si surface results negligibly oxidized, and the C/Fe atomic ratio is close to the value expected for the pure molecular species. The VFC/Si hybrids can be described as (η⁵-C₅H₅)Fe²⁺(η⁵-C₅H₄)-CH₂-CH₂-Si species, on the basis of XPS results. Electrochemical methods have been applied in order to investigate the role played by a robust, covalent Si-C anchoring mode towards substrate-molecule electronic communication, a crucial issue for a perspective development of molecular electronics devices. The response found from cyclic voltammograms for p-Si(1 0 0) functionalized electrodes, run in the dark and under illumination, has shown that the electron transfer is not limited by the number of charge carriers, confirming the occurrence of electron transfer via the Si valence band. The hybrids have shown a noticeable electrochemical stability and reversibility under cyclic voltammetry (cv), and the trend in peak current intensity vs. the scan rate was linear. The molecule-Si bond is preserved even after thousands of voltammetric cycles, although the surface coverage, evaluated from cv and XPS, decreases in the same sequence. An increasingly larger surface concentration of Fe³⁺ at the expenses of Fe²⁺ redox centers has been found at increasing number of cv’s, experimentally associated with the growth of silicon oxide. Surface SiO⁻ groups from deprotonated silanol termination, induced by the electrochemical treatments, are proposed as the associated counterions for the Fe³⁺ species. They could be responsible for the observed decrease in the electron transfer rate constant with electrode ageing.     

Substructuring preconditioners for saddle-point problems arising from Maxwell's equations in three dimensions

This paper is concerned with the saddle-point problems arising from edge element discretizations of Maxwell's equations in a general three dimensional nonconvex polyhedral domain. A new augmented technique is first introduced to transform the problems into equivalent augmented saddle-point systems so that they can be solved by some existing preconditioned iterative methods. Then some substructuring preconditioners are proposed, with very simple coarse solvers, for the augmented saddle-point systems. With the preconditioners, the condition numbers of the preconditioned systems are nearly optimal; namely, they grow only as the logarithm of the ratio between the subdomain diameter and the finite element mesh size.