Three-dimensional adjoint transport calculation module was integrated into S_N transport code ARES in which automated variance reduction parameters are generated based on consistent adjoint driven importance sampling method to accelerate calculation of MCNP5. It shows that automated variance reduction parameters are effective to improve MC calculational efficiency and to produce unbiased statistical results. Automated variance reduction technique with S_N function estimates particle importance more economically and accurately. It avoids obstacles of manual estimation and could be applied for MC simulation of large-scale, complicated shielding problems.

Hybrid boundary node method (hybrid BNM) is applied for thermal analysis of 3D composites. A new formulation is derived for inclusion-based composites, which reduces degrees of freedom (DOFs) compared with conventional multi-domain solver. A version of fast multipole method (FMM) is coupled with new formulation for large scale analysis. Numerical examples are presented to analyze thermal behavior of composites with many inclusions.

Smoothed particle hydrodynamics (SPH) is tested for simulating flow driven by high speed underwater body (FDUB).Developing history of supercavity induced by a body is predicated by SPH.Shape of developed supercavity derived by SPH is compared with result given by principle of independence of cavity sections expansion of Logvinovich.Ghost particle method and Dummy particle method,that remedy boundary deficiency of SPH,are combined to implement wall conditions.

We propose a method for data-visualization of smoothed particle hydrodynamics(SPH).Data calculated with SPH on convex or non convex continuum domain are processed.Particle set is discretized into a set of triangular elements by Delaunay triangulation,and only those elements on which nodes interact with each other are kept.With kept elements and data on their nodes a finite element interpolation is employed to derive values on interior points of the elements.A free surface is easily extracted with the method.Feasibility of the method is tested.Methods dealing with walls,severe particle splashing or multi-materials are provided,respectively.They supply feasible ways for realizing a general post-processing software package for SPH and even for gridless methods.

A multi-dimensional particle-in-cell code KLAP is introduced. Field ionization, impact ionization and two body collision are considered in a one-dimensional code. Moving window technology is used in a three-dimensional code to study laser wakefield acceleration in GeV region. Terahertz radiation, laser propagation in neutral gas medium, surface electron acceleration as well as GeV electron generation in laser wakefield acceleration are studied.