Lattice Boltzmann method for simulation of time-dependent neutral particle transport

In this paper, a novel model is proposed to investigate the neutron transport in scattering and absorbing medium. This solution to the linear Boltzmann equation is expanded from the idea of lattice Boltzmann method (LBM) with the collision and streaming process. The theoretical derivation of lattice Boltzmann model for transient neutron transport problem is proposed for the first time. The fully implicit backward difference scheme is used to ensure the numerical stability, and relaxation time and equilibrium particle distribution function are obtained. To validate the new lattice Boltzmann model, the LBM formulation is tested for a homogenous media with different sources, and both transient and steady-state LBM results get a good agreement with the benchmark solutions.

[1]  William A. Wieselquist,et al.  A cell-local finite difference discretization of the low-order quasidiffusion equations for neutral particle transport on unstructured quadrilateral meshes , 2014, J. Comput. Phys..

[2]  A. Trabelsi,et al.  Method for converting in-situ gamma ray spectra of a portable Ge detector to an incident photon flux energy distribution based on Monte Carlo simulation , 2008 .

[3]  R. Furfaro,et al.  A hybrid method for the solution of linear Boltzmann equation , 2014 .

[4]  S. Van Criekingen,et al.  parafish: A parallel FE–PN neutron transport solver based on domain decomposition , 2011 .

[5]  Gilberto Orengo,et al.  A model of the 14 MeV neutrons source term, for numerical solution of the transport equation to be used in BNCT simulation , 2012 .

[6]  R. P. Smedley-Stevenson,et al.  Linear and quadratic octahedral wavelets on the sphere for angular discretisations of the Boltzmann transport equation , 2005 .

[7]  Benoit Forget,et al.  The OpenMOC method of characteristics neutral particle transport code , 2014 .

[8]  Wenyun Luo,et al.  Energy loss caused by shielding effect of steel cage outside source tube , 2007 .

[9]  Shiyi Chen,et al.  LATTICE BOLTZMANN METHOD FOR FLUID FLOWS , 2001 .

[10]  Gábor Házi,et al.  Large eddy simulation of subchannels using the lattice Boltzmann method , 2007 .

[11]  Benoit Forget,et al.  2D full-core Monte Carlo pin-by-pin burnup calculations with the RMC code , 2014 .

[12]  H. Tan,et al.  The lattice Boltzmann method for one-dimensional transient radiative transfer in graded index gray medium , 2014 .

[13]  Suizheng Qiu,et al.  The development and verification of thermal-hydraulic code on passive residual heat removal system of Chinese advanced PWR , 2006 .

[14]  Wim Sweldens,et al.  Lifting scheme: a new philosophy in biorthogonal wavelet constructions , 1995, Optics + Photonics.

[15]  I. Karlin,et al.  Lattices for the lattice Boltzmann method. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[16]  R. T. Ackroyd.  The why and how of finite elements , 1981 .

[17]  C. Pain,et al.  Chebyshev spectral hexahedral wavelets on the sphere for angular discretisations of the Boltzmann transport equation , 2008 .

[18]  S. Iqbal,et al.  A spatially adaptive grid-refinement approach for the finite element solution of the even-parity Boltzmann transport equation , 2007 .

[19]  John Mandrekas,et al.  GTNEUT: A code for the calculation of neutral particle transport in plasmas based on the Transmission and Escape Probability method , 2004, Comput. Phys. Commun..

[20]  Q. Ma,et al.  Numerical simulation of high-energy neutron radiation effect of scintillation fiber , 2008 .

[21]  Kamal Hadad,et al.  Cellular neural network to the spherical harmonics approximation of neutron transport equation in x–y geometry. Part I: Modeling and verification for time-independent solution , 2011 .

[22]  Gumersindo Verdú,et al.  Solution of the Lambda modes problem of a nuclear power reactor using an h–p finite element method , 2014 .

[23]  J. Boon.  The Lattice Boltzmann Equation for Fluid Dynamics and Beyond , 2003 .

[24]  J. Negele,et al.  Mean free path in a nucleus , 1981 .

[25]  Li-Shi Luo,et al.  Some Progress in Lattice Boltzmann Method. Part I. Nonuniform Mesh Grids , 1996 .

[26]  M. M. Selim,et al.  Solution of the stochastic transport equation of neutral particles with anisotropic scattering using RVT technique , 2009, Appl. Math. Comput..

[27]  Hongchun Wu,et al.  An improved three-dimensional wavelet-based method for solving the first-order Boltzmann transport equation , 2009 .

[28]  D. Henderson,et al.  Numerical benchmark solutions for time-dependent neutral particle transport in one-dimensional homogeneous media using integral transport , 2004 .

[29]  R. Furfaro,et al.  Hybrid-transport point kinetics for initially-critical multiplying systems , 2014 .

[30]  Subhash C. Mishra,et al.  Application of the lattice Boltzmann method for solving the energy equation of a 2-D transient conduction–radiation problem , 2005 .

[31]  B. Ganapol,et al.  A hybrid transport point-kinetic method for simulating source transients in subcritical systems , 2011 .

[32]  R. T. Ackroyd,et al.  Finite Element Methods for Particle Transport: Applications to Reactor and Radiation Physics , 1997 .

[33]  James Tickner,et al.  Arbitrary perturbations in Monte Carlo neutral-particle transport , 2014, Comput. Phys. Commun..

[34]  Jingwei Hu,et al.  An asymptotic-preserving scheme for the semiconductor Boltzmann equation toward the energy-transport limit , 2013, J. Comput. Phys..

[35]  Gábor Házi,et al.  Lattice Boltzmann methods for two-phase flow modeling , 2002 .

[36]  R. P. Smedley-Stevenson,et al.  Goal-based angular adaptivity applied to the spherical harmonics discretisation of the neutral particle transport equation , 2014 .

[37]  Liangzhi Cao,et al.  A new three-dimensional method of characteristics for the neutron transport calculation , 2011 .

[38]  Christopher C. Pain,et al.  Goal-based angular adaptivity applied to a wavelet-based discretisation of the neutral particle transport equation , 2015, J. Comput. Phys..

[39]  Kun Yue,et al.  M-C simulation of shielding effects of PE, LiH and graphite fibers under 1 MeV electrons and 20 MeV protons , 2008 .

[40]  Subhash C. Mishra,et al.  Solving transient conduction and radiation heat transfer problems using the lattice Boltzmann method and the finite volume method , 2007, J. Comput. Phys..

[41]  Graeme A. Bird,et al.  Definition of mean free path for real gases , 1983 .

[42]  Dmitri Priimak,et al.  Finite difference numerical method for the superlattice Boltzmann transport equation and case comparison of CPU(C) and GPU(CUDA) implementations , 2014, J. Comput. Phys..