A parallel and scalable algorithm for ADI method with pre-propagation and message vectorization

It is hard to implement the ADI method in an efficient way on distributed-memory parallel computers. We propose “P-scheme” which parallelizes a tridiagonal linear system of equations for the ADI method, but its effectiveness is limited to the cases where the problem size is large enough mainly because of the communication cost of the propagation phase of the scheme.

In order to overcome this difficulty, we propose an improved version of the P-scheme with “message vectorization” which aggregates several communication messages into one and alleviates the communication cost. Also we evaluate the effectiveness of message vectorization for the ADI method and show that the improved version of the P-scheme works well even for smaller problems and linear and super-linear speedups can be achieved for 8194 × 8194 and 16,386 × 16,386 problems, respectively.     

基于光敏自组装多层膜的液晶光控取向膜

采用layer-by-layer(LBL)的方法，将基板依次在具有光敏双键的季铵盐水溶液和聚乙烯基苯磺酸钠水溶液中浸泡，石英基板表面形成了自组装多层沉积膜。经线性偏振紫外光辐照后，LBL多层膜表现出明显的紫外吸收各向异性。将两片带有各向异性薄膜的基板做成平行液晶器件。在正交的偏光显微镜下观察，发现取得均一、稳定的取向效果。该种光控取向膜有望用作LCoS微显示的取向层。     

Application of piece‐wise linear weight functions for 2D 8‐node quadrilateral element in contact problems

The present study is a continuation of our previous work with the aim to reduce problems caused by standard higher order elements in contact problems. The difficulties can be attributed to the inherent property of the Galerkin method which gives uneven distributions of nodal forces resulting in oscillating contact pressures. The proposed remedy is use of piece‐wise linear weight functions. The methods to establish stiffness and/or mass matrix for 8‐node quadrilateral element in 2D are presented, i.e. the condensing and direct procedures. The energy and nodal displacement error norms are also checked to establish the convergence ratio. Interpretation of calculated contact pressures is discussed. Two new 2D 8‐node quadrilateral elements, QUAD8C and QUAD8D, are derived and tested in many examples, which show their good performance in contact problems. Copyright © 2004 John Wiley & Sons, Ltd.     

电子表格EXCEL5.0在选煤线性规划问题中的应用

采用电子表格软件ＥＸＣＥＬ５．０可以方便地解决许多选煤生产中的线性规划问题。本文通过选煤过程线性规划实例介绍ＥＸＣＥＬ５．０在此领域的应用方法。     

红外遥感用低温制冷器的进展

低温制冷器与红外探测器配套，是军事现代化、高新技术等急需发展的关键产品。国内研制技术也已接近国际先进水平，基本掌握了小型化、低振动、长寿命、细方形漆包线制作、间隙密封、消除电磁干扰等相关关键技术。本文列出了已研制过的直线电机谐振压缩机的主要技术性能。希望主管、用户、研制、生产各方通力合作，应用－改进－再应用－再改进，实现我国低温制冷器的产业化。     

An eight‐node hybrid‐stress solid‐shell element for geometric non‐linear analysis of elastic shells

This paper presents eight‐node solid‐shell elements for geometric non‐linear analysis of elastic shells. To subdue shear, trapezoidal and thickness locking, the assumed natural strain method and an ad hoc modified generalized laminate stiffness matrix are employed. A selectively reduced integrated element is formulated with its membrane and bending shear strain components taken to be constant and equal to the ones evaluated at the element centroid. With the generalized stresses arising from the modified generalized laminate stiffness matrix assumed to be independent from the ones obtained from the displacement, an extended Hellinger–Reissner functional can be derived. By choosing the assumed generalized stresses similar to the assumed stresses of a previous solid element, a hybrid‐stress solid‐shell element is formulated. Commonly employed geometric non‐linear homogeneous and laminated shell problems are attempted and our results are close to those of other state‐of‐the‐art elements. Moreover, the hybrid‐stress element converges more readily than the selectively reduced integrated element in all benchmark problems. Copyright © 2002 John Wiley & Sons, Ltd.     

The effect of particle distribution on damage formation in particulate reinforced metal matrix composites deformed in compression

Image analysis results are reported on the generation of damage in particulate reinforced metal matrix composites during compressive deformation. The technique allows the automated collection of data on the incidence of particle fracture and void formation in the matrix as a function of important microstructural parameters such as local particle volume fraction and particle size. There is a strong relationship between damage and the local volume fraction of the reinforcement proving that damage formation is accentuated in regions of particle clustering. With the SiC reinforced materials examined, there was observed to be a change in dominance of damage mechanism from particle fracture at low local volume fractions to void formation in the matrix within strongly clustered regions. The results are compared with finite element (FE) modelling of the compressive deformation of clustered particles using a simple cluster of equi-spaced particles. The FE results suggest that plastic flow is generally inhibited in clustered regions. In certain highly clustered configurations shielding is such that flow does not occur in the heart of the cluster even at high levels of average plastic strain. The modelling suggests that the change in dominance of damage mechanism is related to the dramatic increase in tensile hydrostatic stresses in the matrix with higher levels of particle clustering.     

The Logarithmic Image Processing Model: Connections with Human Brightness Perception and Contrast Estimators

The logarithmic image processing (LIP) model is amathematical framework based on abstract linear mathematicswhich provides a set of specific algebraic and functionaloperations that can be applied to the processing of intensityimages valued in a bounded range. The LIP model has been provedto be physically justified in the setting of transmitted lightand to be consistent with several laws and characteristics ofthe human visual system. Successful application examples havealso been reported in several image processing areas, e.g.,image enhancement, image restoration, three-dimensional imagereconstruction, edge detection and image segmentation.The aim of this article is to show that the LIP model is atractable mathematical framework for image processing which isconsistent with several laws and characteristics of humanbrightness perception. This is a survey article in the sensethat it presents (almost) previously published results in arevised, refined and self-contained form. First, an introductionto the LIP model is exposed. Emphasis will be especially placedon the initial motivation and goal, and on the scope of themodel. Then, an introductory summary of mathematicalfundamentals of the LIP model is detailed. Next, the articleaims at surveying the connections of the LIP model with severallaws and characteristics of human brightness perception, namelythe brightness scale inversion, saturation characteristic, Weber'sand Fechner's laws, and the psychophysical contrast notion. Finally,it is shown that the LIP model is a powerful and tractable framework for handling the contrast notion. This is done througha survey of several LIP-model-based contrast estimators associated with special subparts (point, pair of points,boundary, region) of intensity images, that are justified bothfrom a physical and mathematical point of view.     

线性时变系统的区间稳定性与鲁棒稳定性

本文应用向量比较定理研究线性时变系统的区间稳定性和具非线性时变摄动的线性时变系统的鲁棒稳定性，所得的新结果包含文献的一些主要结果作为特例，本文的研究方法说明向量比较方法是分析区间稳定性和鲁棒稳定性的一种自然而有力的工具。     

Best L1 Piecewise Monotonic Data Modelling

Let us consider n data measurements of a univariate process that have been altered by random errors. We assume that an underlying model function has a substantially smaller number of turning points than the observed ones. We propose algorithms that make least the sum of the moduli of the errors by requiring k monotonic sections, alternately increasing and decreasing, in the sequence of the smoothed values. The main difficulty in this calculation is that the optimal positions of the joins of the monotonic sections have to be found automatically among so many combinations that it is impossible to test each one separately. Moreover, the calculation seems to be very intractable to general optimization techniques because O(n^k) local minima can occur. It is shown that dynamic programming can be used for separating the data into optimal disjoint sections of adjacent data, where each section requires a single L₁ monotonic calculation. This procedure is highly efficient, requiring at most O(kn²) computer operations and O(n) best L₁ monotonic calculations to subranges of data for a global minimum.