A numerical experimental study of inverse preconditioning for the parallel iterative solution to 3D finite element flow equations

Integration of the subsurface flow equation by finite elements (FE) in space and finite differences (FD) in time requires the repeated solution to sparse symmetric positive definite systems of linear equations. Iterative techniques based on preconditioned conjugate gradients (PCG) are one of the most attractive tool to solve the problem on sequential computers. A present challenge is to make PCG attractive in a parallel computing environment as well. To this aim a key factor is the development of an efficient parallel preconditioner. FSAI (factorized sparse approximate inverse) and enlarged FSAI relying on the approximate inverse of the coefficient matrix appears to be a most promising parallel preconditioner. In the present paper PCG using FSAI, diagonal and pARMS (parallel algebraic recursive multilevel solvers) preconditioners is implemented on the IBM SP4/512 and CLX/768 supercomputers with up to 32 processors to solve underground flow problems of a large size. The results show that FSAI may allow for a parallel relative efficiency larger than 50% on the largest problems with p=32 processors. Moreover, FSAI turns out to be significantly less expensive and more robust than pARMS. Finally, it is shown that for p in the upper range may be much improved if PCG–FSAI is implemented on CLX.     

Quantification of trace elements in human serum fractions by liquid chromatography and inductively coupled plasma mass spectrometry

On‐line coupling of LC and ICP‐MS has been used for fractionation and detection of species of Cu, Fe, I, Se and Zn in human serum. It has been shown that anion exchange chromatography provided better separation capability (both intra‐ and inter‐element) than size‐exclusion chromatography. The mobile phases for ion exchange chromatography consisted of Tris–HNO₃ buffer and ammonium salt (nitrate, acetate or formate). Formate was found to be the best mobile phase counter ion, enabling good chromatographic separation, and is acceptable for mass spectrometry too. The quantitative evaluation of element concentrations adhering to individual fractions was performed by the peak area normalization method. The repeatability of results ranged from 3 to 15% (depending on the element concentration level) and represented the main part of the result uncertainty. The accuracy of Cu and Zn fraction determinations was confirmed by comparison with the isotope dilution technique. Copyright © 2006 John Wiley & Sons, Ltd.     

新疆野生黑加仑与马林果实的营养成分分析

通过对新疆野生黑加仑和野生马林的脂肪油，8种微量元素，17种氨基酸以及维生素C的含量进行测定得出：黑加仑脂肪油含量为17．9％，马林的脂肪油含量为22．3％，其中，黑加仑脂肪油中，r－亚麻酸的含量为19．02％，而马林中，其含量为9．40％，黑加仑中K，Na,Ca,Cu,Zn,Fe的含量比马林中的含量高，而Mn,Mg元素的含量却比马林含量低；马林的维生素C含量为16．22mg/g，黑加仑的含量为108mg/g，黑加仑的氨基酸的总量为1．63％，马林为1．14％。     

A mixed‐FEM formulation for nonlinear incompressible elasticity in the plane

This article deals with an expanded mixed finite element formulation, based on the Hu‐Washizu principle, for a nonlinear incompressible material in the plane. We follow our related previous works and introduce both the stress and the strain tensors as further unknowns, which yields a two‐fold saddle point operator equation as the corresponding variational formulation. A slight generalization of the classical Babu?ka‐Brezzi's theory is applied to prove unique solvability of the continuous and discrete formulations, and to derive the corresponding a priori error analysis. An extension of the well‐known PEERS space is used to define an stable associated Galerkin scheme. Finally, we provide an a posteriori error analysis based on the classical Bank‐Weiser approach. © 2002 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 18: 105–128, 2002     

Solid phase extraction for analysis of biogenic carbonates by electrothermal vaporization inductively coupled plasma mass spectrometry (ETV-ICP-MS): an investigation of rare earth element signatures in otolith microchemistry

Uptake of trace elements into fish otoliths is governed by several factors such as life histories and environment in addition to stock and species differences. In an attempt to elucidate the elemental signatures of rare earth elements (REEs) in otoliths, a solid phase extraction (SPE) protocol was used in combination with electrothermal vaporization (ETV) as a sample introduction procedure for the determinations by inductively coupled plasma quadrupole mass spectrometry (ICP-MS). Effects of various parameters, such as carrier gas flow rate, atomization temperature and chemical modification, were examined for optimization of the conditions by ETV-ICP-MS. Atomization was achieved at 2800 °C. Lower temperatures (i.e. 2600 °C) resulted in severe memory problems due to incomplete atomization. Palladium was used as a chemical modifier. It was found that an increase in Pd concentration up to 0.5 μg in the injection volume (70 μl) led up to four-fold enhancement in the integrated signals. This phenomenon is attributed to the carrier effect of Pd rather than the stabilization since no significant losses were observed for high temperature drying around 700 °C even in the absence of Pd. Preconcentration was performed on-line at pH 5 by using a mini-column of Toyopearl AF-Chelate 650M chelating resin, which also eliminated the calcium matrix of otolith solutions. After preconcentration of 6.4 ml of solution, the concentrate was collected in 0.65 ml of 0.5% (v/v) HNO₃ in autosampler cups, and then analyzed by ETV-ICP-MS. The method was validated with the analysis of a fish otolith certified reference material (CRM) of emperor snapper, and then applied to samples. Results obtained from otoliths of fish captured in the same habitat indicated that otolith rare earth element concentrations are more dependent on environmental conditions of the habitat than on species differences.     

Mass‐lumping or not mass‐lumping for eigenvalue problems

In this article we analyze the effect of mass‐lumping in the linear triangular finite element approximation of second‐order elliptic eigenvalue problems. We prove that the eigenvalue obtained by using mass‐lumping is always below the one obtained with exact integration. For singular eigenfunctions, as those arising in non convex polygons, we prove that the eigenvalue obtained with mass‐lumping is above the exact eigenvalue when the mesh size is small enough. So, we conclude that the use of mass‐lumping is convenient in the singular case. When the eigenfunction is smooth several numerical experiments suggest that the eigenvalue computed with mass‐lumping is below the exact one if the mesh is not too coarse. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 653–664, 2003     

Volcanic monitoring for radon and chemical species in the soil and in spring water samples

Soil radon has been monitored at two fixed stations in the northern flank of Popocatepetl Volcano, a high risk volcano located 60 km SE from Mexico City. Water samples from three springs were also studied for radon as well as major and trace elements. Radon in the soil was recorded using track detectors. Radon in the water samples was evaluated using the liquid scintillation method and an Alphaguard. The major elements were determined through conventional chemical methods and trace elements using an ICP-MS equipment. Soil radon levels were low, indicating a moderate diffuse degassing through the flanks of the volcano. Groundwater radon had almost no relation with the eruptive stages. Water chemistry was stable in the reported time (2000–2002).     

A finite element-spherical harmonics radiation transport model for photon migration in turbid media

In this paper, we solve the steady-state form of the Boltzmann transport equation in homogeneous and heterogeneous tissue-like media with a finite element-spherical harmonics (FE-P_N) radiation transport method. We compare FE-transport and diffusion solutions in terms of the ratio of absorption to reduced scattering coefficient, (μ_a/μ_s′) and the anisotropy factor g. Two different scattering phase function formulas are employed to model anisotropic scattering in the slab media with high g-value. Influence of void-like heterogeneities, and of their boundaries with the surrounding medium on the transport of photons are also examined.     

Mathematical and numerical analysis of a simplified time-dependent viscoelastic flow

A time-dependent model corresponding to an Oldroyd-B viscoelastic fluid is considered, the convective terms being disregarded. Global existence in time is proved in Banach spaces provided the data are small enough, using the implicit function theorem and a maximum regularity property for a three fields Stokes problem. A finite element discretization in space is then proposed. Existence of the numerical solution is proved for small data, so as a priori error estimates, using again an implicit function theorem. Supported by the Swiss National Science Foundation. Fellowship PBEL2–114311.     

数码相机原理与系统设计研究

介绍了数码相机的工作原理 ,分析了数码相机的系统结构 ,提出了设计数码相机的处理流程。典型的数码相机系统有镜头、闪光灯、光学取景器、LCD显示屏、图像数据存储扩展设备接口、图像数据传输接口、供电系统以及核心处理器等八个主要模块。数码相机的数据流向从图像传感器开始 ,止于图像数据的存储和传输。数据流的处理主要包括模数转换、光学黑电平钳位、针对镜头的边缘畸变的运算修正、坏像素处理、白平衡处理、伽马校正、色彩合成处理、边缘检测 (锐度检测 )和伪彩色检测 (伪彩色抑制 )、JPEG压缩和图像存储器等模块