HIGH PERFORMANCE SPARSE SOLVER FOR UNSYMMETRICAL LINEAR EQUATIONS WITH OUT-OF-CORE STRATEGIES AND ITS APPLICATION ON MESHLESS METHODS

A new direct method for solving unsymmetrical sparse linear systems(USLS) arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin (MLPG) method need to solve large USLS. The proposed solution method for unsymmetrical case performs factorization processes symmetrically on the upper and lower triangular portion of matrix, which differs from previous work based on general unsymmetrical process, and attains higher performance. It is shown that the solution algorithm for USLS can be simply derived from the existing approaches for the symmetrical case. The new matrix factorization algorithm in our method can be implemented easily by modifying a standard JKI symmetrical matrix factorization code. Multi-blocked out-of-core strategies were also developed to expand the solution scale. The approach convincingly increases the speed of the solution process, which is demonstrated with the numerical tests.     

横观各向同性体应力波传播问题的特征分析

本文根据广义特征理论,对横观各向同性体轴对称弹性波传播问题进行了理论分析。提出了一个简化特征分析的方法,这种方法将成为应用特征线法求解各类各向异性应力波问题的一个突破口。     

时一空守恒元解元（CE／SE）方法综述

时一空守恒元解元方法是近年来兴起的一种全新的高分辨率守恒型方程计算方法．它具有物理概念清晰、计算精度高和格式构造简单等优点,是一种具有广阔发展前景的计算方法．本文详细地介绍了cE／ss方法的基本原理、发展历史、应用情况和最新进展、并指出了当前研究的不足和发展方向．     

Gradient-augmented hybrid interface capturing method for incompressible two-phase flow

Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method(GAHM) is presented for incompressible two-phase flow.A front tracking method(FTM) is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives(normals and curvatures) and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.     

横观各向同性含液饱和多孔介质中应力波传播的特征分析

根据广义特征理论,对横观各向同性含液饱和多孔介质中应力波传播特性进行了特征分析．给出了特征曲面的微分方程以及沿次特征线的相容条件,得到了波阵面的解析表达式．详细地讨论了应力波在横观各向同性含液饱和多孔介质中传播时,其速度曲面和波阵面的形状及性质．分析结果亦表明,纯固体中应力波传播的特征方程,是含液饱和多孔介质中应力波特征方程的特例．     

Solutions for a Mode Ⅲ Growing Crack in a Piezoelectric Plane Under Two Kinds of Electric Boundary Conditions ＊

The solutions for a mode Ⅲ crack growing along an arbitrary propagation path in a piezoelectric plane are studied under the impermeable surface condition and the electrical contact surface condition respectively. According to the two kinds of electric boundary conditions, the Hilbert and Riemann boundary value problems in a half-plane including opening smooth arc are obtained from the theoretical analysis. Moreover, the equipollence of the solution formed under these two electric boundaries is proved, and unified solutions for the stress and electric displacement distribution in the crack-tipfield of the piezoelectric plane are achieved.     

气相爆轰波数值模拟中化学反应模型研究

 基于改进的时 空守恒元解元算法对气相爆轰波数值模拟中3种常用化学反应模型(二步模型,基元反应模型和Sichel的二步模型)进行了考察。对平面爆轰波和具有胞格结构的爆轰波进行了数值模拟,并对数值结果进行了比较和讨论。结果表明3种化学反应模型得到的爆轰参数准确性有所差异,但得到的胞格结构均能和实验结果较好吻合。3种化学反应模型在爆轰波数值模拟中各有优缺点,应视具体问题决定使用哪种化学反应模型。     

Atomic diffusion across Ni_(50)Ti_(50)Cu explosive welding interface:Diffusion layer thickness and atomic concentration distribution

Molecular dynamics simulations are carried out to study atomic diffusion in the explosive welding process of Ni50Ti50–Cu(at.%). By using a hybrid method which combines molecular dynamics simulation and classical diffusion theory, the thickness of the diffusion layer and the atomic concentration distribution across the welding interface are obtained.The results indicate that the concentration distribution curves at different times have a geometric similarity. According to the geometric similarity, the atomic concentration distribution at any time in explosive welding can be calculated. Ni50Ti50–Cu explosive welding and scanning electron microscope experiments are done to verify the results. The simulation results and the experimental results are in good agreement.     

Tuning of band-gap of phononic crystals with initial confining pressure

<正>The mechanism of the shift of the band-gap in phononic crystal(PC) with different initial confining pressures is studied experimentally and numerically.The experimental results and numerical analysis simultaneously indicate that the confining pressure can efficiently tune the location in and the width of the band-gap.The present work provides a basis for tuning the band-gap of phononic crystal in engineering applications.     

含孔正交各向异性复合材料板的动应力集中

运用一种改进的非结构化四边形格子法,对含孔正交各向异性板条受面内冲击拉伸时弹性应力波的传播过程和孔边的动应力集中进行了研究．非结构化格子法采用与有限元类似的网格剖分方法,并基于围绕每个节点的积分平衡方程,并自然满足复杂边界的自由边界条件．计算中不需存储刚度矩阵,因而计算速度快、效率高、节省内存,在解决应力波传播问题中具有显著的优越性．通过对多种工况进行数值模拟,分析了材料的各向异性性质、纤维方向、孔径比、加载脉冲周期等参数对孔边动应力的影响,得到了一些规律性的结果．并与现有实验结果进行对比,验证了该方法的有效性．