With dynamic coupling model of fractured anisotropic full-tensor permeability, a complete three-dimensional three-phase fluid-solid coupling model for anisotropic ultra-low permeability reservoir is established. Using modular explicit coupling iteration, we establish a fluid & solid coupling seepage numerical simulation method considering fracture deformation and corresponding coupling program. It shows that for fractured reservoir, with change of reservoir pressure in process of mining, rock deformation and anisotropy of reservoir medium strain can make various deformation degree in multiple sets of cracks, cause nonlinear change of fracture permeability, which makes principal value of reservoir fracture permeability tensor rotating, changing flow direction, affecting reservoir development effect.

With gas-kinetic unified algorithm (GKUA) based on Boltzmann equation for various flows regimes, rotational non-equilibrium effect is investigated in Rykov model, in which spin movement of gas molecules is described with moment of inertia. Numerical method of Boltzmann model equation involving rotational non-equilibrium effect is developed. Nitrogen shock wave structure, flows around a 2D blunt-head body and flows around a 3D tine bicone are simulated to validate the method in various flows regimes.

To investigate effect of rotational degree of freedom on gas flows covering various flow regimes,Boltzmann-Rykov model is studied and reduced velocity distribution functions are introduced with quadrature to rotational degree of freedom of molecular velocity distribution function.A single gas-kinetic model equation is translated into simultaneous equations of three reduced velocity distribution functions at discrete velocity ordinate points with discrete velocity ordinate method and numerical integration technique.One-and two-dimensional Boltzmann-Rykov model equations for diatomic gases are computed with finite-difference method of computational fluid dynamics.One-dimensional shock-tube and the two-dimensional flow past erect plate are analyzed in the whole range of Knudsen numbers.Reliability of gas-kinetic unified algorithm(GKUA) is validated in solving one-and two-dimensional flows from free molecular flow to continuum regimes.It is indicated that the gas rarefaction degree and molecular inner degree of freedom affect flow field greatly.And rarefied gas flows with high Knudsen numbers take serious non-equilibrium effect.

An adaptive time integration algorithm is proposed for ADE(alternating direction explicit based on dimension splitting) scheme.Based on the algorithm,a parallel adaptive program for multi-material hydrodynamics is developed on JASMIN.An implosion experiment in ICF(inertial confinement fusion) is simulated on 512 processors.Both simulation results and performance analysis demonstrate correctness and efficiency of the algorithm and the code.

The eigenvalue and the two-particle scattering matrix are obtained in a spin-ladder model. They are solved with the coordinate Bethe Ansatz method in integrable. The eigenvalue equation is formed with the Schrödinger equation. Then the eigenvalue is solved as the wave function is given. Finally, the two-particle scattering matrix is calculated by the eigenvalue equation and the continuous condition of the wave function. The eigenvalues of one-particle, two-particle and N₀-particle are solved, and the two-particle scattering matrix is obtained.

Neutron transport numerical simulation of different energy of neutron sources is performed in hypothetical nuclear warhead models.Three methods of actively warhead detection are analyzed and the feasibility of actively detecting warheads with three types of neutron sources is studied.With 14.1 MeV pulse neutron source,distribution of delayed neutron flow with time is studied.Results show that the warheads can be detected without exposing the sensitive information of nuclear weapons by detecting the delayed neutrons.

The γ ray transport process in some pits of warhead models has been studied.Since the penetrating γ ray spectra are different for different models,the feasibility of using the γ ray from fissile materials as the fingerprint of the warhead has been analyzed.By comparing the characteristics of γ ray's intensities,the thickness of bare ball of Wg.Pu can be estimated if its thickness is less than 1.5 cm.

A proposal of 200MW Inertial Confinement Fusion Fission Hybrid reactor which does not produce plutonium is presented on the basis of one dimensional discrete ordinate transport calculation and burmup calculation by using code BISON 1.5^[4].The time dependent results of energy deposit,²³⁸U(n,γ),239Pu(n,α),⁶Li(n,α) reac-tion rates and K_eff are given.Such a novel reactor proposed here not only could maintain high output power but also could keep critical safety.

Open boundary condition (OBC) is one of the most challenge subjects in the field of computational method for infinite domain. There are already many different kinds of OBCs in wave problems. But, they still lack practical comparisons of effectiveness. This article devotes to such comparisons of several important OBCs at the artificial boundaries for the nonlinear wave equation Boussinesq equation. A test example is given on the diffraction of solitary wave by a hump at the sea bottom. 3-D scattering wave diagrams are also given. The effectiveness is analyzed with instructive conclusion.