Numerical Simulation of 3D-Thermal Field with Phase Transformation

以激光熔凝表面强韧化处理为背景,应用空间的弹塑性有限单元和高精度的数值算法、同时考虑材料组织性能的变化来模拟材料的温度场。主要研究激光熔凝加工中瞬时温度场数值模拟,同时考虑相变潜热的影响,为第二步热应力场及残余应力的数值模拟做准备。最后用算例验证了模型的正确性,并给出了不同时刻温度场的分布。

Teracluster LSSC-II - Its Designing Principles and Applications in Large Scale Numerical Simulations

The teracluster LSSC-II installed at the State Key Laboratory of Scientific and Engineering Computing, Chinese Academy of Sciences is one of the most powerful PC clusters in China. It has a peek performance of 2Tflops. With a Linpack performance of 1.04Tflops, it is ranked at the 43rd place in the 20th TOP500 List (November 2002), 51st place in the 21st TOP500 List (June 2003), and the 82nd place in the 22nd TOP500 List (November 2003) with a new Linpack performance of 1.3Tflops. In this paper, we present some design principles of this cluster, as well as its applications in some largescale numerical simulations.

Numerical Investigation of Toroidal Shock Wave Focusing in a Cylindrical Chamber

In this paper, focusing of a toroidal shock wave propagating from an annular shock tube into a cylindrical chamber was investigated numerically with the dispersion controlled dissipation (DCD) scheme. The first case for an incident Mach number of 1.5 was conducted and compared with experiments for validation. Then, several cases were calculated for higher incident Mach numbers varying from 2.0 to 5.0, and complicated flow structures were observed. The numerical study was mainly focused on two aspects: focusing process and flow structures. The process, including diffraction, focusing, and reflection, is displayed to reveal the focusing mechanism, and the flow structures at different incident. Mach numbers are used to demonstrate shock reflection styles and focusing characteristics.

Numerical Simulation of Transient Thermal Field in Laser Melting Process

Numerical simulation of thermal field was studied in laser processing. The 3-D finite element model of transient thermal calculation is given by thermal conductive equation. The effects of phase transformation latent are considered. Numerical example is given to verify the model. Finally the real example of transient thermal field is given.

Numerical Simulation of Pulsed Liquid Fluidized Bed and Its Experimental Validation

The pulsed liquid fluidized bed was studied using numerical simulation and experimental methods, The area-averaged two-fluid model (TFM) was used to simulate the pulsed fluidization. The bed expansion and collapse processes were simulated first and the phenomena obtained from the calculation were consistent with our previous experiments and observations. In the pulsed fluidization, the variation of bed height, the variations of particle velocity and concentration distribution were obtained and analyzed. Experiments were carried out to validate the simulation results. The pressure variation with time at different locations was measured using pressure transducers and compared with the simulated results. The variations of bed height and particle concentration distribution were recorded using a digital video camera recorder. The results were consistent with the simulation results as a whole.

Numerical Simulation on the Micro-Scale Bending Induced by Pulsed Laser Beam

利用自行研制的含热传导、冲击动力学大、变形有限元程序,模拟了小尺寸梁在脉冲激光加热条件下的变形过程。在此基础上,利用商用程序模拟了冷却及残余应力的产生,研究了激光参数（强度及分布）等对于微弯曲的影响。数值模拟结果与文献中的实验观察相吻合。

Numerical simulation of temperature field evolution of metallic alloy surface due to pulsed laser treatment

This paper performed a numerical simulation on temperature field evolution for the surface layer of a metallic alloy subjected to pulsed Nd:YAG laser treatment. The enthalpy method was adopted to solve the moving boundary problem, I.e. Stefan problem. Computational results were obtained to show the temperature field evolution. Effects of latent heat and mushy zone width on the temperature field were investigated. The results also show very high values of temperature gradient and cooling rate, which are typical characteristics during the solidification process.

Damage Localization, Sensitivity of Energy Release and the Catastrophe Transition

Large earthquakes can be viewed as catastrophic ruptures in the earth’s crust. There are two common features prior to the catastrophe transition in heterogeneous media. One is damage localization and the other is critical sensitivity; both of which are related to a cascade of damage coalescence. In this paper, in an attempt to reveal the physics underlying the catastrophe transition, analytic analysis based on mean-field approximation of a heterogeneous medium as well as numerical simulations using a network model are presented. Both the emergence of damage localization and the sensitivity of energy release are examined to explore the inherent statistical precursors prior to the eventual catastrophic rupture. Emergence of damage localization, as predicted by the mean-field analysis, is consistent with observations of the evolution of damage patterns. It is confirmed that precursors can be extracted from the time-series of energy release according to its sensitivity to increasing crustal stress. As a major result, present research indicates that the catastrophe transition and the critical point hypothesis (CPH) of earthquakes are interrelated. The results suggest there may be two cross-checking precursors of large earthquakes: damage localization and critical sensitivity.

Behaviour analysis of numerical solutions and simulation of coherent structures in compressible mixing layers

For understanding the correctness of simulations the behaviour of numerical solutions is analysed, Tn order to improve the accuracy of solutions three methods are presented. The method with GVC (group velocity control) is used to simulate coherent structures in compressible mixing layers. The effect of initial conditions for the mixing layer with convective Mach number 0.8 on coherent structures is discussed. For the given initial conditions two types of coherent structures in the mixing layer are obtained.

Influence of Coriolis force on bulk flows during horizontal Bridgman growth on a centrifuge: a numerical study

Three-dimensional and time-dependent numerical simulations are performed For melt convection in horizontal Bridgman crystal growth tinder high gravity conditions by means of a centrifuge. The numerical results show that Coriolis Force can cause a stabilizing effect on the fluctuations of the melt flow under a specific relation direction and relation rates of the centrifuge as reported in previous experiments (Ma et al., Materials Processing in High Gravity, Plenum Press, New York, 1994, p. 61). The present simulation provides details of the now features associated with the effect of the Coriolis force. There are also some differences between the present three-dimensional and former two-dimensional numerical solutions particularly in the prediction of the critical conditions and flow patterns.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

Numerical Investigation On Detonation Cell Evolution In A Channel With Area-Changing Cross Section

The two-dimensional cellular detonation propagating in a channel with area-changing cross section was numerically simulated with the dispersion-controlled dissipative scheme and a detailed chemical reaction model. Effects of the flow expansion and compression on the cellular detonation cell were investigated to illustrate the mechanism of the transverse wave development and the cellular detonation cell evolution. By examining gas composition variations behind the leading shock, the chemical reaction rate, the reaction zone length, and thermodynamic parameters, two kinds of the abnormal detonation waves were identified. To explore their development mechanism, chemical reactions, reflected shocks and rarefaction waves were discussed, which interact with each other and affect the cellular detonation in different ways.

Numerical Study Of Heat Transfer Mechanism In Turbulent Supercritical Co2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Numerical Simulation Of Horizontal Migration Of Proppant

The horizontal migration of proppant was numerically investigated with a two-fluid model, in which the interaction between fracturing fluid and proppant, along with that among proppants was taken into account through interphase forces. The migration process and the volumetric concentration of the proppant were examined under various conditions, and the. averaged volumetric concentration of the proppant was obtained. The present research might be useful in the process design of the hydraulic fracturing in the oilfields.