Lie group invariant finite difference schemes for the neutron diffusion equation
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Finite difference techniques are used to solve a variety of differential equations. For the neutron diffusion equation, the typical local truncation error for standard finite difference approximation is on the order of the mesh spacing squared. To improve the accuracy of the finite difference approximation of the diffusion equation, the invariance properties of the original differential equation have been incorporated into the finite difference equations. Using the concept of an invariant difference operator, the invariant difference approximations of the multi-group neutron diffusion equation were determined in one-dimensional slab and two-dimensional Cartesian coordinates, for multiple region problems. These invariant difference equations were defined to lie upon a cell edged mesh as opposed to the standard difference equations, which lie upon a cell centered mesh. Results for a variety of source approximations showed that the invariant difference equations were able to determine the eigenvalue with greater accuracy, for a given mesh spacing, than the standard difference approximation. The local truncation errors for these invariant difference schemes were found to be highly dependent upon the source approximation used, and the type of source distribution played a greater role in determining the accuracy of the invariant difference scheme than the local truncation error.
- Research Organization:
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- W-7405-ENG-36
- OSTI ID:
- 10165908
- Report Number(s):
- LA-12791-T; ON: DE94014980
- Resource Relation:
- Other Information: DN: Thesis submitted to University of Illinois, Urbana-Champaign, IL; TH: Thesis (Ph.D); PBD: Jun 1994
- Country of Publication:
- United States
- Language:
- English
Similar Records
11-th order of accuracy for numerical solution of 3-D Poisson equation with irregular interfaces on unfitted Cartesian meshes
Numerical solution of second-order boundary value problems on nonuniform meshes