Numerical solution of time-dependent two-dimensional P1 equation of neutron transport.

概要

論文の詳細を見る
The time-dependent P1 equation for two-dimensional neutron transport is numerically solved by a finite difference approximation of the explicit form along the bicharacteristics of the P1 equation. Applying von Neumann's stability condition to this numerical procedure in an infinite space, we can derive the condition necessary for the solution to be stable. This condition is that the mesh widths satisfy the inequality 0<λ ?? √3/2 with λ=time mesh Δt/space mesh Δx. or Δz, where the time t is measured in units of inverse neutron speed 1/υ. The sufficient stability condition on the ratio λ is to be determined by numerical experiments. It has been found that the upper bound of λ becomes larger for smaller values of space mesh width. In respect of the stability of numerical solution, the P1 approximation is more advantageous than the diffusion approximation. Transient behavior of neutron flux distribution due to a stationary neutron source is numerically determined assuming zero initial values. After the transient state terminates, the steady state distribution is obtained.
一般社団法人日本原子力学会の論文