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     

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     

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.     

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.     

束晕-混沌的复杂性理论与控制方法及其应用前景

本文系统论述涉及强流加速器等强流离子束装置中产生的束晕-混沌的复杂性理论与控制方法及其应用前景。强流离子束在核材料生产与增殖、洁净核能、放射性废物嬗变、放射性药物生产、重离子聚变、高能物理、核科学与工程、国防与民用工业和医疗等许多方面都有极其重要的应用潜力和诱人的发展前景。尤其是，近年来强流加速器驱动的放射性洁净核能系统是国内外关注的热门课题，因为它比常规核电更安全、更干净、更便宜。但是，强流离子束形成的束晕-混沌的复杂性现象已引起了国内外广泛关注，需要加以抑制、控制和消除这类现象，解决这一难题已经成为强流离子束应用中的关键问题之一。目前不仅必须深入研究这类束晕-混沌的复杂特性及其产生的物理机制，而且需要研究如何实现对束晕-混沌的有效控制，并寻求和发展其新理论、新方法和新技术。这就向强流离子束物理和非线性-复杂性科学及其技术提出了一系列极富挑战性的新课题。本文结合国内外的研究概况，根据我们多年来的研究成果，特别是我们首创性地提出了一些束晕-混沌的有效控制方法，它们包括：非线性反馈控制法，小波反馈控制法，变结构控制法，延迟反馈控制法，参数自适应控制法等，进行重点的介绍。对上述课题当前的主要进展及相关问题进行系统的总结和比较全面综述的评论。最后，指出该领域今后的研究方向，以推动这个崭新领域的深入研究和应用发展。     

Construction algorithms for polynomial lattice rules for multivariate integration

We introduce a new construction algorithm for digital nets for integration in certain weighted tensor product Hilbert spaces. The first weighted Hilbert space we consider is based on Walsh functions. Dick and Pillichshammer calculated the worst-case error for integration using digital nets for this space. Here we extend this result to a special construction method for digital nets based on polynomials over finite fields. This result allows us to find polynomials which yield a small worst-case error by computer search. We prove an upper bound on the worst-case error for digital nets obtained by such a search algorithm which shows that the convergence rate is best possible and that strong tractability holds under some condition on the weights.

We extend the results for the weighted Hilbert space based on Walsh functions to weighted Sobolev spaces. In this case we use randomly digitally shifted digital nets. The construction principle is the same as before, only the worst-case error is slightly different. Again digital nets obtained from our search algorithm yield a worst-case error achieving the optimal rate of convergence and as before strong tractability holds under some condition on the weights. These results show that such a construction of digital nets yields the until now best known results of this kind and that our construction methods are comparable to the construction methods known for lattice rules.

We conclude the article with numerical results comparing the expected worst-case error for randomly digitally shifted digital nets with those for randomly shifted lattice rules.

     

Zero Dissipative, Explicit Numerov-Type Methods for Second Order IVPs with Oscillating Solutions

New explicit, zero dissipative, hybrid Numerov type methods are presented in this paper. We derive these methods using an alternative which avoids the use of costly high accuracy interpolatory nodes. We only need the Taylor expansion at some internal points then. The method is of sixth algebraic order at a cost of seven stages per step while their phase lag order is fourteen. The zero dissipation condition is satisfied, so the methods possess an non empty interval of periodicity. Numerical results over some well known problems in physics and mechanics indicate the superiority of the new method.     

Local Feasible QP-Free Algorithms for the Constrained Minimization of SC1 Functions

We consider the problem of minimizing an SC¹ function subject to inequality constraints. We propose a local algorithm whose distinguishing features are that: (a) a fast convergence rate is achieved under reasonable assumptions that do not include strict complementarity at the solution; (b) the solution of only linear systems is required at each iteration; (c) all the points generated are feasible. After analyzing a basic Newton algorithm, we propose some variants aimed at reducing the computational costs and, in particular, we consider a quasi-Newton version of the algorithm.