We present DG method for solving two-way coupling compressible gas-solid two-phase flow. Equations of both phases are discreted simultaneously, including convection term and source term. Splitting technique to discretize governing equations separately is avoided. Numerical flux of both phases is based on approximate Riemann solver. Dusty-gas shock tube problem with particles in low pressure section is simulated. Comparisons of equilibrium flow and frozen flow are made. Influence of particles in gas and interaction rules between two phases in relaxation zone behind shock are studied. It found that mass ratio of particles determines last equilibrium state and particles diameter determines transition process of two-phase flow from nonequilibrium to equilibrium flow. Namely, different diameter particles correspond to different relaxation time and distance. It shows that the numerical method proposed is reliable. It lays a foundation for more complicated gas-solid two-phase flow problems.

A study on oblique shock wave reflection in condensed matter is carried out by means of numerical simulation and theoretical analysis.Runge-Kutta control volume discontinuous finite element method is used to solve Euler equations.Equations of state for condensed matter adopt "stiffen gas" formulas.Patterns of oblique shock wave reflection in condensed matter are discussed.A shock polar theory is employed in analyzing critical agles of transition from regular reflection to irregular reflection.It gives states of reflected shock wave.Numerical results and shock polar solutions are compared and typical oblique shock wave reflections are obtained.

We simulate radiation effects of solar cosmic protons on spacecraft and satellites with a Monte Carlo software GEANT4.Damage induced by proton irradiation and shielding of Al layer are calculated.Linear energy transfer and range in materials agree well with referenced data. Energy deposited in incident direction shows a Bragg curve and nonelastic interactions show important impact.In Si recoil atoms are mainly distributed along incident proton trajectory and isolated recoils distribute along lateral direction with lower concentration. Al shielding layer shows a shielding on high energetic protons.Shielding is hardly changed as layer thickness is larger than 10mm.On the contrary,secondary production increases obviously.

Runge-Kutta control volume discontinuous finite element method (RKCVDFEM) is proposed to solve numerically hyperbolic conservation laws,in which space discretization is based on control volume finite element method (CVFEM) while time discretization is based on a second order accurate TVB Runge-Kutta technique.Piecewise linear function space is chosen as finite element space.The scheme is total variation bounded (TVB) and is formally second order accurate in space and time.Numerical examples show that numerical solution converges to the entropy solution,and order of convergence is optimal for smooth solution.Compared with numerical results of Runge-Kutta discontinuous Galerkin method (RKDGM) those of RKCVDFEM are better.

We describe an unsteady flow model of blood in arterioles of normal and tumor tissues,and show effect of external magnetic field on flow field of blood considering coupling of magnetic field and flow. Traditional single-phase model is modified by a discrete method,and dynamical simulation with a finite element route is used to investigate transmission of magnetic nano-drug in arteriole.The result indicates that the capture of nano-drug by magnetic field is entirely different from that of microdrug by magnetic field.The trajectoriy of nano-drug is primarily coincident with the streamline of blood flow.However,the effect of magnetic field on the blood flow is obvious due to the presence of magnetic nano-drug.Several whirlpools are formed near the magnet.In tumor arteriole these whirlpools stay for a long time and the magnetic nano-drug is efficiently captured by external magnetic field.

An in-core fuel management code package HEX-ICFM is developed for Xi'an Pulsed Reactor.The orthogonal design model and the loading pattern optimization program HEX-ORTH are studied.The cell calculation is performed using WIMS-D/4.A two dimensional hexagonal geometry multigroup nodal theory code SIXTUS-2 is used for core diffusion calculation.The core physical parameters are calculated using HEX-ICFM for the first cycle of Xi'an Pulsed Reactor.With the fuel of end of cycle 3 and 30 fresh fuel rods the optimum loading patterns are calculated for the objective function max(K^BOC_eff)using HEX-ORTH.

The sampling method by coupling criticality source direction biasing and the exponential transform is established for duct shielding calculation.In sample problem studying,the optimal source direction biasing parameter p₁、exponential transform parameter p₂ are supposed and the variance of duct calculating results by using this method is shown to reduce efficiently.Correction calculation for the duct of Nuclear Power Institute of China (NPIC) pulsed reactor shows the results are more close to the measured ones.

It gives the pulsing parameter formulae based on Fuchs-Hansen's adiabatic model, calculates the kinetic parameters of Xi'an Pulsed Reactor when various reactivities are inserted, and analyses the influences of wat er channel and Mo-Tc target on the core safety parameters.

By using lattice calculation code WIMS-D/4,hexagonal nodal diffusion code SIXTUS and the WIMS-N1/N2 library which include the data of the nuclide H in ZrH,Er-166 and Er-167,nuclear safety parameters are analyzed of uranium zirconium hydride pulsed reactor core,power peaking factor and fuel temperature coefficient.

The block-up coefficient technique and the finite volume approach are used to solve numerically the complete N-S equations in a square enclosure with spacers. A pressure-velocity correction algorithm (SIMPLE-C, Semi-Implicit Method for Pressure Linked Equation-Consistent) is adopted in the direct calculation of all the physical variables. The applicability, feasibility and efficiency of the block-up coefficient technique applied to simulate three-dimensional viscous flows are investigated.The primary results have shown that the block-up coefficient technique is indeed a simple and convenient numerical approach, with encouraging prospects, for computing the velocity, pressure and other physical quantities for the flows in a channel with complex boundary shape.