随机对流扩散方程的数值仿真

本文针对对流一扩散随机过程在随机输入（即随机输运和源项）,作用下进行数值仿真。我们先将对流扩散随机微分方程中的随机函数采用有限项截断的多项式浑沌展开（Polynomial Chaos Expansion）展开,再由Galerkin映射法得到求解浑沌展开系数的确定性方程组。这是一个在物理空间包含多尺度解的大方程组。为此我...     

三维不规则区域热传导问题无网格方法的数值模拟

配点型无网格法是纯无网格法,它不需要任何背景网格,效率高。本文用加权最小二乘配点方法(Weighted Least-Squares Collocation Method-WLSCM)计算不规则区域热传导问题,形函数采用径向基函数近似。通过二维具有分析解的实例表明WLSCM方法精度高,稳定性好且具有较高的计算效率。此外,将WLSCM方法应用于工程中常见的三维不规则区域热传导问题,结果表明:WLSCM方法的计算结果与FLUENT的计算结果符合很好。     

中能重离子碰撞中同位旋相关核子-核子碰撞截面和同位旋相关平均场的探针

在最近几年我们通过使用同位旋相关的量子分子动力学模型 (IQMD)系统的研究了同位旋相关的平均场和介质中核子 核子 (N N)碰撞截面对中能重离子碰撞 (HIC)中碎裂和耗散的同位旋效应。我们发现原子核阻止R和Qzz,中等质量碎片多重性Nimf和质子(中子 )发射数Np(Nn)敏感的依赖于介质中N N碰撞截面的同位旋效应 ,而弱的依赖于同位旋相关的平均场 (对称势 ) ,这些物理量作为提取相对高能范围缺中子系统的同位旋相关介质中N N碰撞截面的探针。我们也可以通过相对低能区到 1 5 0MeV/u的前平衡核子发射中质比来提取关于对称势的知识和讨论它的同位旋依赖性。     

空间网格结构损伤探测的研究进展

空间网格结构损伤探测技术是目前结构工程领域研究的热点和难点．根据国内外的最新文献信息，讨论了损伤探测技术在空间网格结构中的应用及最新的研究进展，对现有文献中的方法进行了分类和评述．重点介绍了目前比较适用的柔度法和模型修正方法，对该领域中尚存在的问题进行了讨论，对今后的发展方向及趋势进行了展望．     

一种全四边形网格生成方法——改进模板法

首先对全四边形单元网格自动剖分算法中的模板法进行了探讨，并提出了相应的改进方法。在此基础上提出了一种新的全四边形单元网格自动生成方法。该方法允许在两个方向上存在网格疏密过渡，并可以提高单元的密度要求自动计算亲单元每条边上的结点数，有效地对局部实施加密处理。     

计算钝感炸药爆轰过程的一种新反应率模型

根据钝感炸药爆轰过程中含有激化过程和慢反应,建立了一种新的反应率模型。与其他反应率模型相比,这种反应率模型可以应用于较粗的网格。在每厘米50个网格条件下,炸药驱动铝飞片和钽飞片的自由面计算结果与实验结果很接近。同时,应用此反应率模型计算了钝感炸药驱动LiF过程,在每厘米50个网格的条件下,炸药与LiF间速度的计算值同实验值接近,且误差随网格尺寸变小而变小。这些表明,此反应率模型能够在较粗的网格条件下,比较准确地描述钝感炸药驱动飞片过程,有利于在工程实际中应用。     

完全变换法在无网格伽辽金方法中的应用

由于移动最小二乘形函数一般不具有常规有限元或边界元形函数所具有的插值特征，本质边界条件的处理成为无网格伽辽金法实施中的一个难点。本文通过建立节点位移和广义位移之间的关系对移动最小二乘形函数进行修正，给出了修正的移动最小二乘形函数；以二维问题为例，对完全交换法在无网格伽辽金方法中的应用进行了研究，实现了本质边界条件在节点处的精确施加。数值计算结果表明该方法不仅简单合理，而且具有较高的精度、收敛性和稳定性。     

利用无网格方法分析斜拉破坏钢筋混凝土梁

从传统的混凝土分离裂缝模型和弥散裂缝模型出发，使用不同的模型模拟不同发展阶段的混凝土裂缝，并利用无网格方法可以灵活布置节点和边界的优点，使分离裂缝模型的使用得到很大的简化。通过在宏观裂缝表面布置基于试验结果的混凝土与混凝土界面单元，正确地模拟了裂缝表面抗剪能力的变化。算例表明，使用本文方法可以准确模拟斜拉破坏混凝土梁的破坏过程，且精度和数值稳定性要高于传统的有限元弥散裂缝方法。     

氢氧混合气体爆轰波的真实化学反应模型数值模拟

采用高精度的ENO格式和基于基元化学反应的真实化学反应模型求解氢氧混合气体一维爆轰波的精细结构。采用直接起爆方法得到稳定传播的爆轰波 ,计算的爆轰波阵面参数和实验相当符合。对爆轰波反应区化学反应的研究表明 ,参与反应的不同组分具有不同类型的变化特征。网格尺寸影响的研究表明 ,计算结果的精度随着网格尺寸的增加而增加 ,并能保持较好的收敛性。移动网格研究结果表明 ,网格运动速度和爆轰速度接近时 ,两者的相互作用对计算结果产生一定影响。     

THE COMPACTION OF TIME-DEPENDENT VISCOUS GRANULAR MATERIALS CONSIDERING INERTIAL FORCES

In the present paper, compactions of time-dependent viscous granular materials are simulated step by step using the automatic adaptive mesh generation schemes. Inertial forces of the viscous incompressible aggregates are taken into account. The corresponding conservation equations, the weighted-integral formulations, and penalty finite element model are investigated. The fully discrete finite element equations for the simulation are derived. Polygonal particles of aggregates are simplified as mixed three-node and four-node elements. The automatic adaptive mesh generation schemes include contact detection algorithms, and mesh upgrade schemes. Solutions of the numerical simulation are in good agreement with some results from literatures. With minor modification, the proposed numerical model can be applied in several industries, including the pharmaceutical, ceramic, food, and household product manufacturing.