Based on C-J (Chapman-Jouguet) theory different reference states of gaseous detonation products and unreacted explosives in detonation problem are considered. According to these reference states, specific Mie-Grüneisen EOS (equation of states) is selected. As chemical reaction process is neglected, a zero-thickness section of guided shock wave exists as an interface in front of the detonation wave. In numerical simulation, evolution of detonation wave includes two parts:Propagation of wave section, as well as interaction with unreacted medium. In the propagation process, speed is defined as the constant detonation speed, and detonation products forms instantly. In the interaction process, Mie-Grüneisen mixture model is employed to simulate continuous impact of detonation wave. With Mie-Grüneisen EOS, as well as the Mie-Grüneisen mixture model, motion of detonation wave is simulated well. Comparing with related theoretical data and numerical results good performance was found.

Numerical study based on a wake oscillator model is conducted to determine response performance of vortex-induced vibration (VIV) on a rigid cylinder. A coupled model of structural oscillator of a rigid cylinder and wake oscillator is established, and then the model is discretized and solved based on a standard central finite difference method of the second order. VIV response characteristics including dimensionless displacement, dimensionless lift, and frequency ratio and lock-in region varied with mass ratios and damping ratios are compared. It can be found that the numerical method simulates response performances of vortex-induced vibration on a rigid cylinder very well. As dimensionless mass ratio increases, onset of lock-in region occurs later, lock-out point occurs earlier, and so the lock-in region becomes smaller.

As a high-temperature phase of Pu, δ-Pu doped with little Ga can stay to room temperature. Crystal structure and electronic structure of systems with different contents of Ga are calculated with density functional theory (DFT) method. Calculations mainly include lattice constant, density, formation energy, density of states(DOS), electron density and Mulliken population. It shows that within studied doping scope, lattice constant of systems decrease and density of systems increase with increase of content of Ga, while stability of system with 6.25% content of Ga is superior to that with 3.125% and 12.5% content of Ga. Ga-doping enhanced locality and strengthens bonding ability of electrons, which, to some degree, reveals electronic mechanism of Ga stabilizing δ-Pu. Character of Pu-Ga bond is metallic state, and interaction is mainly contributed by Pu 7s, 6p, 6d and Ga 4s, 4p orbital electrons. While interaction is relatively weak, doped systems maintain fine mechanical properties and machining performance. Contribution of Ga to stability of δ-Pu lies mainly in its improvement on bonding performance of Pu, instead of its immediate bonding effect with Pu.

A multi-heat source ring-shaped heat-supply network is not a plane net and supply-return pipe net is dissymmetric under faulty condition. We present a method for topology structure of spatial heat-supply network based on graph theory. With object-oriented programming,it defines graph and includes correlation arithmetic according to basic structure information of supply-return pipe net and builds incidence matrix and fundamental circuit matrix automatically. It creates a spatial heat-supply network "tree"which includes information of heat source and heat users and adopts breadth-first search based on plane graph. According the tree,basic incidence matrix and fundamental circuit matrix of spatial heat-supply network are established. In simulating heating network of Harbin development zone,a program of multi-heat source ring-shaped heat-supply network in topology structure is build and computing performance is verified. It provides a method for complex heat-supply network hydrodynamic condition analysis and dissymmetric heatsupply network calculation based on arithmetic of graph and adapts to various model pipe nets.

Taking vorticity equation as a controlling equation, we simulate three vortexes developing progress with initial condition of Gaussian probability vorticity distribution. The wavelet expression of vorticity contains scaling coefficient term and wavelet coefficient tenn. The latter can be divided into a significant wavelet coefficient term and a trivial wavelet coefficient term reference to a given critical value. The critical value contains practical flow information. The vorticity is reconstructed only with scaling coefficients and significant wavelet coefficients. It saves computing time and capture most part of entrophy. Numerical results show that scaling coefficients together with significant wavelet coefficients is less than 10% of total wavelet coefficients, but they contains more than 99% of total enstrophy.