任意多连通域网格自动生成算法及其应用

给出了一种任意平面域内三角形网格自动生成的算法，本算法是根据下面两个事实而设计出来的：１、任意平面多连通域总可以剖分成若干个三角形单元。２、平面域内的活节点数目随着单元剖分总会趋近于零。文中还介绍了几种网格性能改善的方法，整个程序结构紧凑，需要人工输入的数据量少，生成的网格质量较高，而且处稳定性能好。     

三维自适应FE-SPH耦合算法在多层间隔金属靶侵彻问题中的应用

鉴于有限元算法不能有效地模拟侵彻过程所产生的金属碎片, 本文中基于三维自适应FE-SPH耦合算法的基本理论, 自主开发了模拟多层间隔金属靶侵彻问题的三维FE-SPH耦合计算程序。该程序采用四面体单元对多层间隔金属靶侵彻模型进行初始离散, 计算过程中, 当四面体单元等效塑性应变超过某一设定值时, 单元自动转化为SPH粒子, 并引入有限单元-粒子接触算法和耦合算法, 实现大变形和破碎区域采用SPH方法计算, 克服有限元法单元畸变存在的问题。多层间隔靶侵彻算例分析表明, 三维FE-SPH耦合计算程序采用等效塑性应变作为转化判据计算结果较稳定, 并且能够有效地再现侵彻过程中所产生的碎片, 能够模拟侵彻碎片对后层靶的毁伤效应。     

二阶矩两相湍流模型中的改进算法

建立了融合SIMPLEC算法在内的考虑了浓度修正影响的颗粒压力修正方程。提出了二阶矩湍流模型下考虑浓度修正值影响的两相湍流流动的算法,并将它和没有考虑浓度修正值影响的二阶矩湍流模型进行了对照。结果表明在二阶矩模型中是否考虑浓度修正影响会不同程度地影响流场的速度以及浓度等参数分布,考虑浓度修正影响的二阶矩湍流模型更能有效地预测稠密两相湍流流动。     

MA对偶-信赖域算法在非线性不等式约束优化问题中的应用研究

针对含有非线性不等式约束条件的优化问题,提出了MA对偶-信赖域算法。在每次迭代过程中,基于信赖域方法和问题的逼近属性,构造了原优化问题中目标函数和约束函数的移动渐进线函数,由此建立简单的子优化问题。运用对偶方法求解子问题得到原优化问题的下降方向,再用线搜索方法取得搜索步长,最后得到下一步的迭代点。应用数学推理证明了该算法的全局收敛性。以悬臂梁最小柔度问题为例,应用MA对偶-信赖域算法对优化问题进行了求解,数值算例的结果表明,MA对偶-信赖域算法在求解非线性约束优化问题时比MMA和GCMMA算法的迭代次数少,收敛速度快。     

无网格算法在多段翼型流动计算中的应用

研究了一种求解欧拉方程的无网格算法，发展出了一套布点及点云自动生成的方法；在点云离散的基础上，采用最小二乘法求解矛盾方程的方法来求取空间导数，进而获得数值通量；采用四步龙格-库塔方法进行时间推进，并引入当地时间步长和残值光顺等加速收敛措施。通过对NA-CA0012翼型的跨音速流动和多段翼型复杂绕流的数值模拟，验证了上述无网格算法的正确性和实用性。     

CAD软件在工程地质三维建模中的应用

如何快速、准确地建立地质体的三维模型一直是众多岩土工程数值模拟工作者所面临的难题。虽然三维地学模拟软件具有很好的三维地质建模能力,但是由于数据结构的差异,采用他们现行三维地学模拟软件建立的地质模型难以导入数值模拟分析软件中,以为相应工程问题的数值模拟服务。目前,随着各种CAD、CAM软件行业的的飞速发展,涌现出了许多优秀的三维建模软件,而且这些软件大都与现行数值分析软件有着良好的数据接口功能。据此,本文提出了采用现行CAD软件来建立工程地质体的三维模型,使得建立的模型达到既"可视"又"可算"的目的。将其应用于云南某高速公路边坡的三维建模中,证明了该法具有方便、快捷和合理等优点。     

RNGк—ε模式在工程湍流数值计算中的应用

自Yakhot和Orszag提出RNGк-ε湍流模式以来，很多学者对其进行了修正和发展。并将其应用于某些实际湍流问题的数值模拟。取得了一些与实验近似一致的结果。本文主要从工程湍流计算的角度出发。结合作者的部分研究工作。对近年来RNGк-ε模式在湍流流动数值研究中的应用现状及进展情况进行了总结，指出了该模式存在的问题，并对其应用前景进行了展望。     

RNGκ-ε模式在工程湍流数值计算中的应用

自Yakhot和Orszag提出RNG k-ε湍流模式以来,很多学者对其进行了修正和发展,并将其应用于某些实际湍流问题的数值模拟,取得了一些与实验近似一致的结果。本文主要从工程湍流计算的角度出发,结合作者的部分研究工作,对近年来RNG k-ε模式在湍流流动数值研究中的应用现状及进展情况进行了总结,指出了该模式存在的问题,并对其应用前景进行了展望。     

拉格朗日松弛算法在平台罗经多故障诊断技术中的应用研究

最优多故障诊断问题是一个NP-hard问题。针对平台罗经这一复杂系统，采用有向图描述元件与测试点间的因果依赖关系，并建立系统的多信号模型。在考虑元件发生故障的先验概率的前提下，提出一种基于拉格朗日算法（LRA）和子梯度优化算法（SOA）近最优多故障诊断算法，并在某型平台罗经的方位稳定系统多故障诊断中得到应用。结果表明：该方法相对于传统的单故障诊断方法诊断速度快，适用于平台罗经这类大型复杂系统的多故障诊断。     

A parallel implicit domain decomposition algorithm for the large eddy simulation of incompressible turbulent flows on 3D unstructured meshes

We present a parallel fully implicit algorithm for the large eddy simulation (LES) of incompressible turbulent flows on unstructured meshes in three dimensions. The LES governing equations are discretized by a stabilized Galerkin finite element method in space and an implicit second-order backward differentiation scheme in time. To efficiently solve the resulting large nonlinear systems, we present a highly parallel Newton-Krylov-Schwarz algorithm based on domain decomposition techniques. Analytic Jacobian is applied in order to obtain the best achievable performance. Two benchmark problems of lid-driven cavity and flow passing a square cylinder are employed to validate the proposed algorithm. We then apply the algorithm to the LES of turbulent flows passing a full-size high-speed train with realistic geometry and operating conditions. The numerical results show that the algorithm is both accurate and efficient and exhibits a good scalability and parallel efficiency with tens of millions of degrees of freedom on a computer with up to 4096 processors. To understand the numerical behavior of the proposed fully implicit scheme, we study several important issues, including the choices of linear solvers, the overlapping size of the subdomains, and, especially, the accuracy of the Jacobian matrix. The results show that an exact Jacobian is necessary for the efficiency and the robustness of the proposed LES solver.     

Simulation of 2D external viscous flows by means of a domain decomposition method using an influence matrix technique

Two-dimensional external viscous flows are numerically approximated by means of a domain decomposition method which combines a vortex method and a finite difference method. The vortex method is used in the flow region which is dominated by convective effects, whereas the finite difference method is used in the flow region where viscous diffusion effects are dominant. An influence matrix technique combined with the uniformity condition of the pressure is used to enforce the tangential velocity boundary condition. Comparisons between numerical and experimental data show that the method is well adapted for simulating two-dimensional flows.     

Investigation of temperature statistic characteristics in turbulent water flow with unsteady heat transfer

The results of an experimental study of a temperature field and its statistical characteristics in turbulent water flow upon a sudden change of heat generation in the channel wall are reported. Measurements were performed in 5 mm × 40 mm, 10 mm × 40 mm, and 20 mm × 40 mm channels in the regions of thermal stabilization and stabilized heat transfer at Reynolds numbers of (0.8–6.8) × 10⁴. The measurement results are generalized using a dimensionless time scale. The results of the calculation of heat transfer coefficients at unsteady heat transfer are presented.     

Numerical modeling of the turbulent viscous-gas flow past a cylinder with a circular vortex cell in the presence of suction from the central body surface

The compressibility effect on the cylinder drag reduction due to air suction through the surface of a central body in a circular vortex cell is estimated on the basis of the solution of the steady Reynolds equations closed by the shear stress transfer model, together with the continuity, energy, and state equations.     

无网格法在三维复杂超音速流场模拟中的应用研究

对最小二乘无网格方法在含复杂外形三维超音速流场中的应用进行了研究.选用分解法求解采用最小二乘法得到的对称方程组,针对最小二乘无网格方法的计算特点生成近似正交均匀分布的离散点,对B1AC2R常规导弹超音速流场采用最小二乘无网格方法进行了无粘数值模拟,计算了B1AC2R常规导弹在不同攻角下的轴向力、法向力及俯仰力矩系数,并将数值结果与实验结果进行了比较.结果表明,最小二乘无网格方法在求解含复杂外形超音速流场时具有较高的准确度,将其应用于三维含复杂外形超音速流场的模拟是完全可行的.     

应变花埋丝技术在齿科修复学上的应用研究

本文通过采用三维应变花埋丝技术，对固定义齿修复的牙周应力分布进行了定量分析研究，绘制了修复前后牙周的应力分布随基牙倾斜角度的变化曲线，测得下颌第二磨牙用作固定桥基牙时近中倾斜的最大角度为７０度，从生物力学的角度证明，固定义齿修复有利于倾斜基牙牙周组织健康、为齿科临床上广泛开展倾斜基牙固定修复提供了依据。     

On the interfacial roughness scale in turbulent stratified two-phase flow: 3D lattice Boltzmann numerical simulations with forced turbulence and surfactant

Numerical 3D simulations of turbulent, stratified two-phase shear flow with a surfactant laden interface were used to test and develop a phenomenological interfacial roughness scale model where the energy required to deform the interface (buoyancy, interfacial tension, and viscous work) is proportional to the turbulent kinetic energy adjacent to the interface.The turbulence was forced in the upper and lower liquids in the simulations, to emulate the interfacial dynamics without requiring (prohibitively) large simulation domains and Reynolds numbers. The addition of surfactant lead to an increased roughness scale (for the same turbulent kinetic energy) due to the introduction of interfacial dilatational elasticity that suppressed horizontal motion parallel to the interface, and enhanced the vertical motion.The phenomenological roughness scale model was not fully developed for dilatational elasticity in this work, but we proposed a source term that represents surfactant induced pressure fluctuations near the interface. This source term should be developed further to account for the relation between surfactant density fluctuations and turbulence adjacent to the interface. We foresee that the roughness scale model can be used as a basis for more general interfacial closure relations in Reynolds averaged turbulence models, where also mobile surfactant is accounted for.     

自由面势流问题的域外奇点边界元法及其数值误差分析

讨论了域外奇点边界元法在自由面势流问题计算中的作用，并以连续及离散Fourier分析对该方法(就m阶面元的一般情况)进行数值误差分析，导出了计及面元阶数、奇点至自由面垂向距离、配置点移动、差分格式等因素影响的数值误差一般表达式。从理论上证明了自由面势流问题计算中采用域外奇点法可改善离散产生的数值色散误差并能结合配置点前移(向上游)等方法以数值满足辐射条件。     

用刚体元方法求空间极限平衡问题的上限解

刚体元方法是研究块体系的变形和运动的一种数值方法,可以用它来研究弹塑性物体的极限平衡问题,文中用通用的二维刚体元程序计算了正方形柱体顶面中心区域作用方形基础载荷情况的上限解,结果表明,这种方法较好的精度和广泛的适用性。