Numerically stable Monte Carlo-burnup-thermal hydraulic coupling schemes

This paper presents stochastic implicit coupling method intended for use in Monte-Carlo (MC) based reactor analysis systems that include burnup and thermal hydraulic (TH) feedbacks. Both feedbacks are essential for accurate modeling of advanced reactor designs and analyses of associated fuel cycles. In particular, we investigate the effect of different burnup-TH coupling schemes on the numerical stability and accuracy of coupled MC calculations. First, we present the beginning of time step method which is the most commonly used. The accuracy of this method depends on the time step length and it is only conditionally stable. This work demonstrates that even for relatively short time steps, this method can be numerically unstable. Namely, the spatial distribution of neutronic and thermal hydraulic parameters, such as nuclide densities and temperatures, exhibit oscillatory behavior. To address the numerical stability issue, new implicit stochastic methods are proposed. The methods solve the depletion and TH problems simultaneously and use under-relaxation to speed up convergence. These methods are numerically stable and accurate even for relatively large time steps and require less computation time than the existing methods. © 2013 Elsevier Ltd. All rights reserved.

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.

薄层水流水力特性试验研究

对粘砂床面的薄层水流进行定床试验,研究了不同流量、坡度下薄层水流的水力学参数(流态、水深、流速及阻力系数)随坡度和流量的变化规律。结果表明:水流雷诺数Re随着坡度和流量的增大而增大,但Re主要是由流量决定的,坡度对Re的影响不大,水流流态基本上在过渡流区和紊流区;平均水深-单宽流量,流速-单宽流量,水流阻力系数-单宽流量均成很好的幂函数形式,平均水深和流速均随单宽流量的增大而增大,阻力系数随着单宽流量的增大而减小。

Studies of a scintillator-bar detector for neutron wall at externaltarget facility

To achieve a better time resolution of a scintillator-bar detector for a neutron wall at the external target facility of HIRFL-CSR,we have carried out a detailed study of the photomultiplier,the wrapping material and the coupling media. The timing properties of a scintillator-bar detector have been studied in detail with cosmic rays using a high and low level signal coincidence. A time resolution of 80 ps has been achieved in the center of the scintillator-bar detector.