Behavior of the Diamond Difference and LowOrder Nodal Numerical Transport Methods in the Thick Diffusion Limit for Slab Geometry
Abstract
The objective of this work is to investigate the thick diffusion limit of various spatial discretizations of the onedimensional, steadystate, monoenergetic, discrete ordinates neutron transport equation. This work specifically addresses the two lowest order nodal methods, AHOTN0 and AHOTN1, as well as reconsiders the asymptotic limit of the Diamond Difference method. The asymptotic analyses of the AHOTN0 and AHOTN1 nodal methods show that AHOTN0 does not possess the thick diffusion limit for cell edge or cell average fluxes except under very limiting conditions, which is to be expected considering the AHOTN0 method limits to the Step method in the thick diffusion limit. The AHOTN1 method, which uses a linear incell representation of the flux, was shown to possess the thick diffusion limit for both cell average and cell edge fluxes. The thick diffusion limit of the DD method, including the boundary conditions, was derived entirely in terms of cell average scalar fluxes. It was shown that, for vacuum boundaries, only when σ _{t}, h, and Q are constant and σ _{a} = 0 is the asymptotic limit of the DD method close to the finitedifferenced diffusion equation in the system interior, and that the boundary conditions between the systems willmore »
 Authors:
 Pennsylvania State Univ., University Park, PA (United States)
 Publication Date:
 Research Org.:
 Knolls Atomic Power Laboratory (KAPL), Niskayuna, NY (United States)
 Sponsoring Org.:
 USDOE
 OSTI Identifier:
 903208
 Report Number(s):
 LM07K035
TRN: US200719%%642
 DOE Contract Number:
 AC1200SN39357
 Resource Type:
 Thesis/Dissertation
 Country of Publication:
 United States
 Language:
 English
 Subject:
 42 ENGINEERING; BOUNDARY CONDITIONS; CROSS SECTIONS; DIAMONDS; DIFFUSION; DIFFUSION EQUATIONS; DISCRETE ORDINATE METHOD; EXACT SOLUTIONS; GEOMETRY; NEUTRON TRANSPORT; SCALARS; SCATTERING; SERIES EXPANSION; TRANSPORT
Citation Formats
Gill, Daniel Fury. Behavior of the Diamond Difference and LowOrder Nodal Numerical Transport Methods in the Thick Diffusion Limit for Slab Geometry. United States: N. p., 2007.
Web. doi:10.2172/903208.
Gill, Daniel Fury. Behavior of the Diamond Difference and LowOrder Nodal Numerical Transport Methods in the Thick Diffusion Limit for Slab Geometry. United States. doi:10.2172/903208.
Gill, Daniel Fury. Tue .
"Behavior of the Diamond Difference and LowOrder Nodal Numerical Transport Methods in the Thick Diffusion Limit for Slab Geometry". United States.
doi:10.2172/903208. https://www.osti.gov/servlets/purl/903208.
@article{osti_903208,
title = {Behavior of the Diamond Difference and LowOrder Nodal Numerical Transport Methods in the Thick Diffusion Limit for Slab Geometry},
author = {Gill, Daniel Fury},
abstractNote = {The objective of this work is to investigate the thick diffusion limit of various spatial discretizations of the onedimensional, steadystate, monoenergetic, discrete ordinates neutron transport equation. This work specifically addresses the two lowest order nodal methods, AHOTN0 and AHOTN1, as well as reconsiders the asymptotic limit of the Diamond Difference method. The asymptotic analyses of the AHOTN0 and AHOTN1 nodal methods show that AHOTN0 does not possess the thick diffusion limit for cell edge or cell average fluxes except under very limiting conditions, which is to be expected considering the AHOTN0 method limits to the Step method in the thick diffusion limit. The AHOTN1 method, which uses a linear incell representation of the flux, was shown to possess the thick diffusion limit for both cell average and cell edge fluxes. The thick diffusion limit of the DD method, including the boundary conditions, was derived entirely in terms of cell average scalar fluxes. It was shown that, for vacuum boundaries, only when σt, h, and Q are constant and σa = 0 is the asymptotic limit of the DD method close to the finitedifferenced diffusion equation in the system interior, and that the boundary conditions between the systems will only agree in the absence of an external source. For a homogeneous medium an effective diffusion coefficient was shown to be present, which was responsible for causing numeric diffusion in certain cases. A technique was presented to correct the numeric diffusion in the interior by altering certain problem parameters. Numerical errors introduced by the boundary conditions and material interfaces were also explored for a tworegion problem using the Diamond Difference method. A discrete diffusion solution which exactly solves the onedimensional diffusion equation in a homogeneous region with constant cross sections and a uniform external source was also developed and shown to be equal to the finitedifferenced diffusion discretization for c = 1 in the system interior, where again the boundary conditions again only agree in the absence of an external source. It was also shown that for c < 1 the exact discrete diffusion solution is written in terms of hyperbolic functions, with expressions which limit to the exact solution for the c = 1 case as c → 1. Finally, a transport discretization is developed which reproduces the exact S2 solution for the case of a purely scattering homogeneous region with vacuum boundary conditions, and an extension to the discretization for the case of c < 1 is found by considering a Taylor series expansion of the exact answer centered at c = 0.},
doi = {10.2172/903208},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Tue May 01 00:00:00 EDT 2007},
month = {Tue May 01 00:00:00 EDT 2007}
}

Nodal approximations of varying order by energy group for solving the diffusion equation
The neutron flux across the nuclear reactor core is of interest to reactor designers and others. The diffusion equation, an integrodifferential equation in space and energy, is commonly used to determine the flux level. However, the solution of a simplified version of this equation when automated is very time consuming. Since the flux level changes with time, in general, this calculation must be made repeatedly. Therefore solution techniques that speed the calculation while maintaining accuracy are desirable. One factor that contributes to the solution time is the spatial flux shape approximation used. It is common practice to use the samemore » 
Development of advanced nodal diffusion methods for modern computer architectures
A family of highly efficient multidimensional multigroup advanced neutrondiffusion nodal methods, ILLICO, were implemented on sequential, vector, and vectorconcurrent computers. Threedimensional realistic benchmark problems can be solved in vectorized mode in less than 0.73 s (33.86 Mflops) on a Cray XMP/48. Vectorconcurrent implementations yield speedups as high as 9.19 on an Alliant FX/8. These results show that the ILLICO method preserves essentially all of its speed advantage over finitedifference methods. A selfconsistent higherorder nodal diffusion method was developed and implemented. Nodal methods for global nuclear reactor multigroup diffusion calculations which account explicitly for heterogeneities in the assembly nuclear properties weremore » 
Nodal methods for problems in fluid mechanics and neutron transport
A new highaccuracy, coarsemesh, nodal integral approach is developed for the efficient numerical solution of linear partial differential equations. It is shown that various special cases of this general nodal integral approach correspond to several high efficiency nodal methods developed recently for the numerical solution of neutron diffusion and neutron transport problems. The new approach is extended to the nonlinear NavierStokes equations of fluid mechanics; its extension to these equations leads to a new computational method, the nodal integral method which is implemented for the numerical solution of these equations. Application to several test problems demonstrates the superior computational efficiencymore »