다차원 2계 중성자 수송방정식의 방향근사를 위한 영역상수법

Piecewise-Constant Method for Angular Approximation for the Second-Order Multidimensional Neutron Transport Equations

  • 노태완 (홍익대학교 기초과학부 응용수학과)
  • 발행 : 2007.02.28

초록

특정한 방향에 대해 방향중성자속(angular neutron flux)을 정의하는 방향차분 방정식(discrete-ordinates or $S_{N}$ equation)과 달리 방향변수를 구분된 방향영역에 대하여 적분한 값을 사용하고, 해당 방향영역 내에서 방향중성자속이 일정하다고 가정하는 영역상수법(piecewise-constant method)을 개발하였다. 기존 방향차분법과 본 연구에서 개발된 영역상수법을 1계 수송방정식(1'st-order Boltzmann transport equation)과 2계 우성 방정식(even-parity equation)에 적용하여 방향차분 방정식인 $S_{N}$ 방정식과 유사 방향차분방정식($S_{N}$-like equation)인 $PC_{N}$ 방정식을 유도하였다. 우성 방정식에 영역상수법을 적용한 경우 기존 방향차분법의 단점인 광첨두 현상(ray effect)이 현저히 감소함을 확인하였는데 이는 우성 방정식의 혼합 미분항의 기여도가 작아지기 때문인 것으로 판단된다. 이러한 이론은 우성 방정식에서 혼합 미분항이 제거된 단순우성 방정식(simplified even-parity equation)을 사용하는 경우 광첨두 현상이 완전 제거 또는 극단적으로 감소되었던 이전의 결과를 이론적으로 설명한다.

The piecewise constant angular approximation is developed to replace the conventional angular quadrature sets in the solution of the second-order, multi-dimensional $S_{N}$ neutron transport equations. The newly generated quadrature sets by this method substantially mitigate ray effects and can be used in the same manner as the conventional quadrature sets are used. The discrete-ordinates and the piecewise-constant approximations are applied to both the first-order Boltzmann and the second-order form of neutron transport equations in treating angular variables. The result is that the mitigation of ray effects is only achieved by the piecewise-constant method, in which new angular quadratures are generated by integrating angle variables over the specified region. In other sense, the newly generated angular quadratures turn out to decrease the contribution of mixed-derivative terms in the even-parity equation that is one of the second-order neutron transport equation. This result can be interpreted as the entire elimination or substantial mitigation of ray effect are possible in the simplified even-parity equation which has no mixed-derivative terms.

키워드

참고문헌

  1. Lewis, E.E.; Miller, Jr. W.F. Computational Methods of Neutron Transport, John Wiley & Sons, 1993
  2. Azmy, Y.Y.; 'Multidimensional Nodal Transport Methods for Multiple-Instruction Multiple-Data, Distributed Memory Machines', Trans. Am. Nucl. Soc., 1988, 56, 292
  3. Lawrence, R.D. 'Progress in Nodal Methods for the Solution of the Neutron Diffusion and Transport Equations', Prog. Nucl. Energy, 1986, 17, 271 https://doi.org/10.1016/0149-1970(86)90034-X
  4. Haghighat, A.; Azmy, Y.Y. 'Parallelization of a Spherical SN Algorithm Based on the Spatial Domain Decomposition', Proceedings of the Conference on Advances in Mathematics, Computation, and Reactor Physics, Pittsburgh, PA, 1991, April 28-May 1
  5. Miller, Jr. W.F.; Noh, T.W. 'Finite-Differences versus Finite elements in Slab Geometry, Even-Parity Transport Theory', Transport Theory and Statistical Physics, 1993, 22
  6. Noh, T.W.; Miller, Jr. W.F.; Morel, J.E. 'Improved Approximations Applied to the $S_N$ Even-Parity Trasport Equations', Trans. Am. Nucl. Soc., 1993, 69, 214
  7. Noh, T.W., Miller, Jr. W.F.; Morel, J.E. 'The Even- Parity and Simplified Even-Parity Transport Equation in Two-Dimensional X-Y Geometry', Nucl. Sci. Eng., 1996, 123, 38 https://doi.org/10.13182/NSE96-A24211
  8. Morel, J.E.; McGhee, J.M. 'A Self-Adjoint Angular Flux Equation', Nucl. Sci. Eng., 1999, 132, 312-325 https://doi.org/10.13182/NSE132-312
  9. Lathrop, K.D. 'Ray Effects in Discrete Ordinates Equations', Nucl. Sci. Eng., 1968, 32, 357 https://doi.org/10.13182/NSE68-4
  10. Lathrop, K.D. 'Remedies for Ray effects', Nucl. Sci. Eng., 1971, 45, 255 https://doi.org/10.13182/NSE45-03-255
  11. Larsen, E.W.; McGhee, J.M.; Morel, J.E. 'The Simplified PN Equations as an Asymptotic Limit of the Transport Equation', Trans. Am. Nucl. Soc., 1992, 66, 231