Adaptive Mesh Refinement FEM for Damage Evolution of Heterogeneous Brittle Media

Damage evolution of heterogeneous brittle media involves a wide range of length scales. The coupling between these length scales underlies the mechanism of damage evolution and rupture. However, few of previous numerical algorithms consider the effects of the trans-scale coupling effectively. In this paper, an adaptive mesh refinement FEM algorithm is developed to simulate this trans-scale coupling. The adaptive serendipity element is implemented in this algorithm, and several special discontinuous base functions are created to avoid the incompatible displacement between the elements. Both the benchmark and a typical numerical example under quasi-static loading are given to justify the effectiveness of this model. The numerical results reproduce a series of characteristics of damage and rupture in heterogeneous brittle media.

Difference schemes on non-uniform mesh and their application

High order accurate schemes are needed to simulate the multi-scale complex flow fields to get fine structures in simulation of the complex flows with large gradient of fluid parameters near the wall, and schemes on non-uniform mesh are desirable for many CFD (computational fluid dynamics) workers. The construction methods of difference approximations and several difference approximations on non-uniform mesh are presented. The accuracy of the methods and the influence of stretch ratio of the neighbor mesh increment on accuracy are discussed. Some comments on these methods are given, and comparison of the accuracy of the results obtained by schemes based on both non-uniform mesh and coordinate transformation is made, and some numerical examples with non-uniform mesh are presented.

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.

钠原子激光增强电离光谱研究—以石墨杯为原子化器 原子吸收光谱法对发样中Zn、Cu、Mn、Ap的测定

、I钠原子激光增强电离光谱(LEIS)方法的研究-以石墨杯为原子化器在LEIS方法中，最常见的原子化器是火焰。但由于火焰背景噪声严重且难以克服，在火焰原子化过程中，雾化和热离解不充分，仅有10~(-2)%的分析溶液参与吸收以及火焰气体使测定元素受到高度稀释等不利因素的影响，火焰原子化限制了LEIS方法灵敏度的进一步提高。考虑到石墨炉原子化器较火焰具有取样量少，绝对灵敏度高；样品（包括固态、液态）可直接引入石墨炉内；不会发生如同火焰中所存在的干扰效应；蒸发效率和原子化效率较高，几乎全部样品都能参与吸收等优点，本工作在已建立火焰LEIS方法的实验基础上，将原子化器改换为石墨杯进行了钠原子LEIS方法的研究。到目前为止，国内外仅有的几篇有关石黑炉LEIS的研究报告中，都报导了该方法对钠原子的检出限估计可达到10~(-14)-10~(-15)克，由于此项研究尚处于探索性研究阶段，故有关方法性的系统研究几乎还未见报导。本工作在未使用任何放大器的情况下（实验条件限制）对影响钠原子LEIS信号强度的诸因互进行了实验观察。主要包括：钠原子化条件；激光束位置、阳极电压、激光输出能量、电极位置以及激光脉冲重复率对LEIS信号强度的影响等。并绘制了校准曲线，统计方法的相对标准偏差分别为11%（高浓度）18.2%（低浓度），在现有仪器条件下，还不能测出检出限，测定下限为3*10~(-9)克。对固体粉末直接进行了尝试，检测下限为5*10~(-8)克，进样是为5毫克。在进一步的研究工作中，如有条件使用低噪声的放大器及Bxear积分器，选择门检时间窗，或采用分步激发等手段，估计本方法定会达到预想的高灵敏度，检敏度至少提高了个数量级。对石墨炉原子化LEIS法来说，似比较详细的研究报告，截至实验停止时还未见报导。II原子吸收光谱法对发样中Zn、Cu、Mn、Al的测定发中微量元素ZN、Cu、Mn均属人体必需元素，与人体的生长发育和多种生理功能，临床医学等方面有着极为密切的关系，而Al则被认为是异致某种疾病的元素之一。本工作报告了用火焰法测定Zn、Cu；石墨炉法测定Al、Mn的结果，其中，对Al的测定，为摆脱基体干扰，加入改进剂Mg(NO_3)_2，并采用平台石墨炉进行试验，得到了线性较好的工作曲线，但在实际测定时，由于实验条件的限制，只能采用一般石墨管加基体改进剂对少娄样品中Al含量进行测定。Zn、Cu、Al三种元素由标准曲线法测定；而Mn由于Fe的干扰无法消除而采用标准加入法测定，并因此限制了测定样品数。Cu、ZN、Mn三种元素的回收率分别为102.8%, 99.7%, 102.5%，变异系数为9.6%, 11.3%, 9.7%，对本地居民发中（30个发样）Zn、Cu含量进行测定，Zn、Cu的含量范围为148-318ppm，7.2-15ppm，并计算了Zn/Cu比。本方法对发样中四种元素的测定结果与ICP法进行对照。两种方法测定结果吻合得较好。

OPTICAL AND STRUCTURAL-PROPERTIES OF AP-MOVPE GAINP/GAAS HETEROSTRUCTURES

Lattice matched GaInP/GaAs heterostructures were grown by atmospheric pressure-metal organic vapor phase epitaxy (AP-MOVPE). Compositional intermixing of As/P and Ga/In near the heterointerfaces was studied by photoluminescence (PL) spectroscopy. Indium segregation, memory effect of In into GaAs and the carry-over of As in the GaInP layer during the growth process were considered as three major factors giving rise to the anomalous emissions in the PL spectra. Both thermal annealing and zinc doping strongly enhanced the compositional interdiffusion near the heterointerfaces.

A multi-moment finite volume formulation for shallow water equations on unstructured mesh

A novel and accurate finite volume method has been presented to solve the shallow water equations on unstructured grid in plane geometry. In addition to the volume integrated average (VIA moment) for each mesh cell, the point values (PV moment) defined on cell boundary are also treated as the model variables. The volume integrated average is updated via a finite volume formulation, and thus is numerically conserved, while the point value is computed by a point-wise Riemann solver. The cell-wise local interpolation reconstruction is built based on both the VIA and the PV moments, which results in a scheme of almost third order accuracy. Efforts have also been made to formulate the source term of the bottom topography in a way to balance the numerical flux function to satisfy the so-called C-property. The proposed numerical model is validated by numerical tests in comparison with other methods reported in the literature. (C) 2010 Elsevier Inc. All rights reserved.