A Numerical Model of Rapid Solidification Processing of Ni-Al Alloy in Planar Flow Casting

A numerical model has been developed for simulating the rapid solidification processing (RSP) of Ni-Al alloy in order to predict the resultant phase composition semi-quantitatively during RSP. The present model couples the initial nucleation temperature evaluating method based on the time dependent nucleation theory, and solidified volume fraction calculation model based on the kinetics model of dendrite growth in undercooled melt. This model has been applied to predict the cooling curve and the volume fraction of solidified phases of Ni-Al alloy in planar flow casting. The numerical results agree with the experimental results semi-quantitatively.

Numerical Simulation of Transient Thermal Field and Residual Stress in Laser Melting Process

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

An analysis of collective damage for short fatigue cracks based on equilibrium of crack numerical density

Collective damage of short fatigue cracks was analyzed in the light of equilibrium of crack numerical density. With the estimation of crack growth rate and crack nucleation rate, the solution of the equilibrium equation was studied to reveal the distinct feature of saturation distribution for crack numerical density. The critical time that characterized the transition of short and long-crack regimes was estimated, in which the influences of grain size and grain-boundary obstacle effect were investigated. Furthermore, the total number of cracks and the first order of damage moment were discussed.

Numerical Simulation of Gaseous Detonation of H_2-O_2 Mixture with Detailed Chemical Reaction Model

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

Direct Numerical Simulation of Passive Scalar in Decaying Compressible Turbulence

The passive scalars in the decaying compressible turbulence with the initial Reynolds number (defined by Taylor scale and RMS velocity) Re=72, the initial turbulent Mach numbers (defined by RMS velocity and mean sound speed) Mt=0.2-0.9, and the Schmidt numbers of passive scalar Sc=2-10 are numerically simulated by using a 7th order upwind difference scheme and 8th order group velocity control scheme. The computed results are validated with different numerical methods and different mesh sizes. The Batchelor scaling with k(-1) range is found in scalar spectra. The passive scalar spectra decay faster with the increasing turbulent Mach number. The extended self-similarity (ESS) is found in the passive scalar of compressible turbulence.

Numerical simulation of standing solitons and their interaction

Standing soliton was studied by numerical simulation of ifs governing equation, a cubic Schrodiger equation with a complex conjugate term, which was derived by Miles and was accepted. The value of linear damping in Miles equation was studied. Calculations showed that linear damping effects strongly on the formation of a standing soliton and Laedke and Spatschek stable condition is only a necessary condition, but not a sufficient one. The interaction of two standing solitons was simulated. Simulations showed that the interaction pattern depends on system parameters. Calculations for the different initial condition and its development indicated that a stable standing soliton can be fanned only for proper initial disturbance, otherwise the disturbance will disappear or develop into several solitons.

Numerical and Analytical Study on the Pull-in Instability of Micro-Structure under Electrostatic Loading

The one-mode analysis method on the pull-in instability of micro-structure under electrostatic loading is presented. Taylor series are used to expand the electrostatic loading term in the one-mode analysis method, which makes analytical solution available. The one-mode analysis is the combination of Galerkin method and Cardan solution of cubic equation. The one-mode analysis offers a direct computation method on the pull-in voltage and displacement. In low axial loading range, it shows little difference with the established multi-mode analysis on predicting the pull-in voltages for three different structures (cantilever, clamped-clamped beams and the plate with four edges simply-supported) studied here. For numerical multi-mode analysis, we also show that using the structural symmetry to select the symmetric mode can greatly reduce both the computation effort and the numerical fluctuation.

Numerical approximation improves well survey calculation

It is now possible to improve the precision of well survey calculations by order of magnitude with numerical approximation.

Although the most precise method of simulating and calculating a wellbore trajectory generally requires more calculation than other, less-accurate methods, the wider use of computers in oil fields now eliminates this as an obstacle.

The results of various calculations show that there is a deviation of more than 10 m among the different methods of calculation for a directional well of 3,000 m.¹ Consequently, it is important to improve the precision and reliability of survey calculation-the fundamental, necessary work of quantitatively monitoring and controlling wellbore trajectories.

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

Numerical simulation of magnetic field design for damping thermocapillary convection in a floating half-zone

The magnetic fields produced by electrical coils are designed for damping the the thermocapillary convection in a floating half-zone in microgravity. The fields are designed specially to reduce the flow near the free surface and then in the melt zone by adjusting the longitudinal coil positions close to the melt zone. The effects of the designed magnetic fields on reducing the flow velocity and temperature distribution non-uniformity in the melt zone are stronger than those of the case of an uniform longitudinal magnetic field obtained by numerical simulation, particularly at the melt-rod interface. It brings fundamental insights into the heat and mass transfer control at the solidification interface by the magnetic field design for crystal growth by the floating full-zone method.

An analysis of subgrid-resolved scale interactions with use of results from direct numerical simulations

Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range; By decomposing the subgrid energy transfer and nonlinear interaction into 'forward' and 'backward' groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.

Numerical Simulation of Transient Thermal Features in the Remelting Processing of Stainless Steel

应用有限元方法对层流等离子体射流不锈钢表面重熔工艺中的瞬态热物理现象进行了数值模拟研究.针对不同加热距离,确定了材料熔化和凝固过程中的瞬态温度场、温度梯度和凝固率的时间和空间分布特征.通过引入等效温度面积密度概念,研究了不锈钢重熔热处理的适合条件.结果表明,9～13mm的范围是较为适宜的加热距离,该结果与试验观察基本符合.

Scheme Construction with Numerical Flux Residual Correction (NFRC) and Group Velocity Control (GVC)

For simulating multi-scale complex flow fields like turbulent flows, the high order accurate schemes are preferred. In this paper, a scheme construction with numerical flux residual correction (NFRC) is presented. Any order accurate difference approximation can be obtained with the NFRC. To improve the resolution of the shock, the constructed schemes are modified with group velocity control (GVC) and weighted group velocity control (WGVC). The method of scheme construction is simple, and it is used to solve practical problems.

Numerical simulation of rock failure and earthquake process on mesoscopic scale

On the basis of the lattice model of MORA and PLACE, Discrete Element Method, and Molecular Dynamics approach, another kind of numerical model is developed. The model consists of a 2-D set of particles linked by three kinds of interactions and arranged into triangular lattice. After the fracture criterion and rules of changes between linking states are given, the particle positions, velocities and accelerations at every time step are calculated using a finite-difference scheme, and the configuration of particles can be gained step by step. Using this model, realistic fracture simulations of brittle solid (especially under pressure) and simulation of earthquake dynamics are made.

Numerical Solutions for the Transient Flow in the Homogenous Closed Circle Reservoirs

There are many fault block fields in China. A fault block field consists of fault pools. The small fault pools can be viewed as the closed circle reservoirs in some case. In order to know the pressure change of the developed formation and provide the formation data for developing the fault block fields reasonably, the transient flow should be researched. In this paper, we use the automatic mesh generation technology and the finite element method to solve the transient flow problem for the well located in the closed circle reservoir, especially for the well located in an arbitrary position in the closed circle reservoir. The pressure diffusion process is visualized and the well-location factor concept is first proposed in this paper. The typical curves of pressure vs time for the well with different well-location factors are presented. By comparing numerical results with the analytical solutions of the well located in the center of the closed circle reservoir, the numerical method is verified.