Implementation of particle-in-cell stellar dynamics codes on the connection machine-2

The development of massively parallel supercomputers provides a unique opportunity to advance the state of the art inN-body simulations. TheseN-body codes are of great importance for simulations in stellar dynamics and plasma physics. For systems with long-range forces, such as gravity or electromagnetic forces, it is important to increase the number of particles toN 10⁷ particles. Significantly improved modeling ofN body systems can be expected by increasingN, arising from a more realistic representation of physical transport processes involving particle diffusion and energy and momentum transport. In addition, it will be possible to guarantee that physically significant portions of complex physical systems, such as Lindblad resonances of galaxies or current sheets in magnetospheres, will have an adequate population of particles for a realistic simulation. Particle-mesh (PM) and particle-particle particle-mesh (P³M) algorithms present the best prospects for the simulation of large-scaleN-body systems. As an example we present a two-dimensional PM simulation of a disk galaxy that we have developed on the Connection Machine-2, a massively parallel boolean hypercube supercomputer. The code is scalable to any CM-2 configuration available and, on the largest configuration, simulations withN = 128 M = 2²⁷ particles are possible in reasonable run times.     

Computing Green's function of the initial-boundary value problem for the wave equation in a layered cylinder

A new analytical method for the approximate computation of the time-dependent Green's function for the initial-boundary value problem of the three-dimensional wave equation on multi-layered bounded cylinder is suggested in this paper. The method is based on the derivation of eigenvalues and eigenfunctions for an ordinary differential equation with piecewise constant coefficients, and an approximate computation of Green's function in the form of the Fourier series with a finite number of terms relative to the orthogonal set of the derived eigenfunctions. The computational experiment confirms the robustness of the method.     

受控振动柱体对波浪影响的数值模拟研究

以不可压缩黏性流体的Navier-Stokes方程为控制方程,采用RNG k-ε湍流模型,利用VOF方法追踪运动的自由表面的方法,建立了三维数值波浪与振动柱体相互作用的模型.基于经典绕射理论验证了波浪绕固定柱体传播的数值模型的可靠性,再采用动边界技术赋予柱体振动规律,研究受控振动柱体对波浪传播的影响,并通过控制柱体初始振动相位变化,系统地给出了在柱体振动的不同初始相位下柱体前后波浪变化的规律.     

载流导体的二维温度场的分析

为分析机电设备导体中电流引起的热分布问题,建立了机电设备的简化有限元模型,对导体进行了自适应网格划分,采用有限元方法求解热导微分方程,得到导体中热分布情况,可以估算导体的形状和尺寸,为机电设备选择合适的导体提供了一种有效的设计方法。     

湍流流动及燃烧模拟研究进展

湍流流动往往伴随着燃烧过程。湍流燃烧现象广泛存在于自然界及诸多实际工程问题中，成为近几年的研究热点。现有的相关理论及模型纷繁复杂。本文对该领域研究进行了分类总结、综合评述。可对全面了解研究现状和进一步的发展前景提供方向性的指导和帮助。     

离子交换过程数学模型及仿真研究

根据离子交换原理建立了离子交换系统模型。并以钨离子交换为例。编写仿真软件求解该模型，获得淋洗曲线及各交换参数，通过与实验结果验证，计算机仿真数据与实验结果大致一致。     

HYSYS过程模拟软件在空分设备操作中的应用

描述了HYSYS过程模拟软件在空分设备操作中的应用,并介绍了软件应用过程中的体会。     

AOD冶炼脱碳反应模型与计算机模拟

AOD冶炼中,终点碳和终点温度的命中率是提高生产率和产品质量的关键。由初始温度,初始成分及补加合金量等原始数据,通过本模型可预报要达到目标终点碳所需要的脱碳时间以及相应的终点温度。本模型计算得到的终点碳与实际终点碳,在同等条件下,两者极为接近,可作为一般AOD冶炼中预报终点碳与脱碳时间关系之用。通过分级调整O2、Ar、N2气比例得到最佳脱碳效果,终点碳命中率可提高到98％以上,最高温度可以控制不超过1750℃,延长炉衬寿命。