大量椭球颗粒随机分布三维区域的模拟及其四面体网格快速生成算法

针对颗粒为椭球形状并且随机分布的多相复合材料给出了三维情形下的计算机模拟方法,并对此复杂区域提供了一种四面体网格快速生成算法。此方法的基本思想:首先是根据材料的要求,给出随机颗粒材料的计算机模拟,然后对模拟后的颗粒边界初始化,得到一些初始剖分点,并给出剖分区域的初始网格(正四面体),接着把颗粒的初始边界点投入到剖分区域中,替换区域中靠此点近的点,并且采用加点寻找边界的方法和修正初始点方法来解决四面体穿刺边界问题,从而得到整个区域的四面体网格。最后用实例说明了算法的有效性。     

混合边界粒子法在分析平板弯曲问题中的应用

本文提出一种平板弯曲问题的数值方法—混合边界粒子法.该方法无需事先进行形状函数假设,直接从平板弯曲问题的基本微分方程出发,通过在局部区域内将微分方程转化为积分方程并利用数值积分,得到指定点的离散解.通过选择合理的指定点途径,最终使得任意指定点的数值可以用边界点上的未知量表达,从而大幅度减少了未知量的个数和系数的计算时间...     

基于光滑粒子流体动力学的滑油泵数值仿真

基于光滑粒子流体动力学方法,构建齿轮泵壳体及内部流体的粒子模型,对内啮合齿轮泵在不同工况下的流量特性进行数值模拟。首先设置了均匀分布的油泵出入口压强,计算得到的流量结果与试验结果吻合,且在中低转速下流量与转速呈线性关系;针对高转速下油泵流量降低的问题,通过适当减小周期模型内部的流体粒子数反映流量降低,获得了与试验值相吻合的连续转速流量的模拟结果;针对影响齿轮泵性能的间隙和空化现象,结合模型特点给出了相应的近似处理方法。通过以上研究,将SPH方法成功地应用于滑油泵问题的分析计算。     

区域重力异常辅助导航定位的统计分析

海洋区域的重力异常变化非常复杂,在利用重力异常数据进行水下载体的辅助导航的仿真中,有些航线我们可以得到较好的匹配结果,有些航线因为辅助区域的重力异常变化较小而不能提供较好辅助定位.利用区域内均匀分布的航线进行多航线的统计分析,可以给出比较客观的区域匹配结果.本文探讨了重力异常辅助惯性导航的相关技术,利用2'×2'的卫星测高重力异常数据,在纬度25°～130°,经度130°～135°的范围对SITAN算法的区域导航能力进行了统计分析.多条航线的统计结果表明,利用SITAN算法进行的重力异常辅助导航可以满足重力辅助导航要求.     

VISCOUS FLOWFIELD BASED ON DISCONTINUOUS BOUNDARY ELEMENT METHOD AND VORTEX METHOD

基于非协调边界元方法和涡方法的联合应用, 模拟了二维和三维黏性不可压缩流场. 计算中利用离散涡元对漩涡的产生、凝聚和输送过程进行模拟, 并将整体计算域分解为采用涡泡模拟的内部区域和用涡列模拟的数字边界层区域. 计算域中涡量场的拉伸和对流由Lagrangian涡方法模拟, 用随机走步模拟涡量场的扩散. 内部区域涡元涡量场速度由广义Biot-Savart公式计算, 势流场速度则采用非协调边界元方法计算. 非协调边界元将所有节点均取在光滑边界处, 从而避免了法向速度的不连续现象; 而对于系数矩阵不对称的大型边界元方程组,引入了非常高效的预处理循环型广义极小残余(the generalized minimum residual, GMRES)迭代算法, 使得边界元法的优势得到了充分发挥, 同时, 在内部涡元势流场计算中对近边界点采用了正则化算法, 该算法将奇异积分转化为沿单元围道上一系列线积分, 消除了势流计算中速度及速度梯度的奇异性. 二维、三维流场算例证明了所用方法的正确性, 也验证了该算法可以大幅度提高模拟精度和效率.     

SPH方法Delaunay三角刨分与自由液面重构

光滑粒子法(SPH)作为一种拉格朗日型无网格方法,兼具欧拉网格方法和拉格朗日网格方法的优势,已经成功应用于科学和工程的众多领域。SPH方法后处理一般基于无规则分布的粒子,不如网格类方法后处理简便、直接。另外,SPH方法模拟自由液面流动等问题时,通过粒子位置难以重构自由液面的准确位置。发展一种基于Delaunay三角刨分的SPH后处理方法,即先基于SPH粒子位置利用Delaunay三角刨分建立三角网格,然后将粒子信息转化成网格单元/节点信息,从而可以在三角网格上进行后处理,实现基于网格方法的后处理功能,并可以在三角网格上直接提取或重构自由液面。将本文的方法应用到液滴碰撞和溃坝流SPH模拟结果的后处理中,得到了非常好的结果,表明本文的方法有效可靠。     

平面图形几何性质的实用计算方法

借助格林公式和积分公式, 利用节点坐标给出了截面图形几何性质计算的一 种新方法. 利用给出的一般表达式, 对于任意直边边界区域都可以得到精确结果.     

太阳强烈辐射对大跨度球面网壳静力性能的影响

在太阳辐射下大跨度建筑物的外部温度场呈非均匀分布,阴面区域的温度低于阳面区域的温度.均匀温度场和非均匀温度场假定给出的分析结果不同.外部温度场受太阳辐射、风速、涂层材料及周围物体辐射反射率等因素的影响.采用简化方法实现大跨度空间钢结构在太阳辐射下外部温度场模拟,得到一个地区的球面网壳外部温度场的分布规律.使用有限元法分析外部非均匀温度场对网壳结构受力性能的影响.算例分析表明,均匀温度分布假定不适用于大跨度双层球面网壳结构的设计.     

2—D问题条形传递函数方法与有限元法的分区耦合

将条形传递函数法（ＳＤＴＦＭ）和有限元法（ＦＥＭ）结合起来,给出了一种求解弹性２－Ｄ问题的新方法。该方法通过把二维解区域解成多个子区域,利用ＳＤＴＦＭ建立矩形子区域（超级单元）基于边界点的刚度矩阵和结点力矢量,而对于其它几何形状的子区域则有限元法建立刚度矩阵和结些点力矢量,从而将ＳＤＴＦＭ推广到任意几何形状的平面区域,克服了ＳＤＴＦＭ只能用于规则几何平面区域的不足,与单纯用有限元法求解相比较,本文     

基于SPH方法模拟30mm线膛炮波纹罩压垮过程的力学参数研究

为研究高旋转对30 mm聚能装药破甲作用的影响机理及自旋补偿原理,采用LS-DYNA有限元软件中的SPH方法模拟30 mm线膛炮波纹罩的压垮过程,得到粒子的实际运动可分解成向心运动与绕中心圆周切线运动,提出压垮过程的4个阶段:压垮前期、缓冲期、波纹槽区域粒子速度增大期和中心粒子相互作用期。射流形成层沿逆时针方向旋转,而形成杵体的材料以相反方向旋转。结果表明:波纹罩特殊设计可以补偿旋转扰动对30 mm聚能装药侵彻作用的负面影响。     

A model of liquid phase sintering by the homogenization

We study the first stage of liquid phase sintering, when the particles rearrangement due to capillary forces is over. We give the boundary value problem satisfied by the displacement field of points of the medium in the phase of elastic compression of solid particles, for given capillary forces acting as a density of external forces, by using the homogenization method and we characterize the mechanical behavior of this constrained medium from the material properties of each elementary components.     

Investigation of intermittency measurement methods for transitional boundary layer flows

Three turbulent intermittency methods, namely the , TERA (turbulent energy recognition algorithm), and M-TERA (modified turbulent energy recognition algorithm) methods, for identifying the intermittent flow characteristics associated with boundary layer transition from laminar to turbulent were considered and compared. The data used were obtained from hot-wire measurements in transitional boundary layer flows on a concave surface with a 2-m radius of curvature and on a flat plate. Comparisons show that the and TERA methods are more sensitive to the choice of threshold constants than the M-TERA method. In terms of the intermittency distribution across the boundary layer, the values obtained by the and TERA methods are unrealistically high in the near-wall region, while those obtained by the M-TERA method are more realistic. In the outer boundary layer region and outside the boundary layer, the and M-TERA methods give reasonable intermittency values, whereas the TERA method produces unrealistically high values in the region outside the boundary layer. In addition, the M-TERA method provides a sharper definition of theend of transition.     

AGAIN DISCUSSING ABOUT SINGULAR PERTURBATION OF GENERAL BOUNDARY VALUE PROBLEM FOR HIGHER ORDER ELLIPTIC EQUATION CONTAINING TWO-PARAMETER

In this paper using the method of "The Two-Variable Expansion Procedure" we again discuss the construction of asymptotic expression of solution of general boundary value problem for higher order ellitptic equation containing two-parameter whose boundary condition is more general than [1]. We give asymptotic expression of solution as well as the estimation corresponding to the remainder term.     

二维边界层方程的迭代求解

针对二维边界层方程,提出了分析积分迭代法.首先将该方法应用于Blasius方程和Falkner-Skan方程的求解,数值计算结果稳定,计算精度高;然后对外部有势流不能达到自相似要求、必须二维求解的二维层流边界层问题,在分析积分迭代法中加上计算力学的松弛迭代法,形成了一套有效的算法.数值结果表明该方法用于层流边界层的计算是很有效的.     

An anti-plane shear problem

We consider an anti-plane shear of an elastic cylinder with a non-convex stored energy function. So, we look for solutions of a non-convex problem of the calculus of variations with Dirichlet boundary conditions. We give sufficient conditions on the boundary data to get existence or non-existence results for this non-convex problem. We also prove some uniqueness results for the relaxed problem associated with the initial problem.     

SOLVING THE FREE BOUNDARY PROBLEM IN CONTINUOUS CASTING BY USING BOUNDARY ELEMENT METHOD

ＳＯＬＶＩＮＧＴＨＥＦＲＥＥＢＯＵＮＤＡＲＹＰＲＯＢＬＥＭＩＮＣＯＮＴＩＮＵＯＵＳＣＡＳＴＩＮＧＢＹＵＳＩＮＧＢＯＵＮＤＡＲＹＥＬＥＭＥＮＴＭＥＴＨＯＤＬｉＹａｏｙｏｎｇ（李耀勇）；ＺｈａｎｇＺｈｉｌｉ（张自立）（ＲｅｃｅｉｖｅｄＪｕｎｅ，１８，１９...     

两相流显微PIV/PTV系统的开发

开发了一个能同时测量两相流中两相速度和细颗粒尺寸分布的显微PIV/PTV系统,其硬件系统包括大功率连续激光器、显微镜、高速摄像机;软件系统由改进的球形颗粒图像识别算法、各种图像处理算法和各种先进的PIV/PTV算法组成。其中改进的圆弧识别算法能够进行更准确地进行曲线分割而能对充满噪音并相互重叠的颗粒图像给出较好的识别结果。应用该PIV系统,可以在微秒和微米数量级上捕获细颗粒/气泡图像,并能较准确地同时得到两相速度、颗粒尺寸和浓度分布。对焚香可吸入颗粒物进行了速度和尺寸的同时测量,得到了较满意的结果。     

Integrating thermal-concentration smoothed profile with lattice Boltzmann methods for simulating sedimentation of nonisothermal circular particles

A thermal-concentration smoothed profile-lattice Boltzmann method is proposed to study the effect of the concentration field on the dynamic behavior of nonisothermal cylindrical particles during the sedimentation process. The velocity, temperature, and concentration equations are solved using the lattice Boltzmann method. Moreover, the smoothed profile method is employed to enforce the nonslip boundary condition as well as constant temperature and constant concentration boundary conditions at the particles surfaces. Moreover, the Boussinesq approximation is used to couple the velocities, temperatures, and concentrations fields. The proposed combined method is validated by comparing the present numerical results with those found in the literature, showing good consistency. Then, the effect of the concentration buoyancy on the behavior of nonisothermal particles is discussed. In addition, the effect of Prandtl, Schmidt, and thermal Grashof numbers on the settling process is investigated. The results show that, by adding the effect of concentration, the maximum settling velocity of hot particles is reduced more relative to the cold ones; accordingly, the cold particles are settled faster than the hot ones. Finally, the sedimentation of two particles in a container at high thermal Grashof is investigated. It is shown that, at high thermal Grashof, there is an intense competition between the buoyancy force and gravity for the hot particles. The buoyancy flow generated leads to the reversal of the drafting-kissing-tumbling motion of the hot particles, making the particles move upward.     

基于格子 Bhatnagar-Gross-Krook模型的地震压力波模拟

给出一种新的用于模拟地震压力波的格子Boltzmann模型. 通过使 用Chapman-Enskog展开和多重尺度技术，得到了一系列的格子Boltzmann方程和时间 尺度$t_0$上守恒律，给出了满足地震压力波方程所要求的高阶矩以及简单的平衡态分布 函数表达式. 数值结果表明这种方法可以用来模拟地震压力波.     

Identification of a moving boundary for a heat conduction problem in a multilayer medium

In this paper, we give a uniqueness theorem for the moving boundary of a heat problem in a composite medium. Through solving the Cauchy problem of heat equation in each subdomain, we finally find an approximation to the moving boundary for one-dimensional heat conduction problem in a multilayer medium. The numerical scheme is based on the use of the method of fundamental solutions and a discrete Tikhonov regularization technique with the generalized cross-validation choice rule for a regularization parameter. Numerical experiments for five examples show that our proposed method is effective and stable.