EP2215617A1

EP2215617A1 - Reactor dosimetry applications using a parallel 3-d radiation transport code

Info

Publication number: EP2215617A1
Application number: EP08833567A
Authority: EP
Inventors: Gianluca Longoni
Original assignee: Westinghouse Electric Co LLC; Westinghouse Electric Corp
Current assignee: Westinghouse Electric Co LLC; CBS Corp
Priority date: 2007-09-27
Filing date: 2008-09-25
Publication date: 2010-08-11
Also published as: CN101861609B; JP2010540936A; EP2215617A4; KR101545035B1; CN101861609A; WO2009042747A1; KR20100066572A; JP5460601B2

Abstract

The invention relates generally to a method for the calculation of radiation field distributions employing a new parallel 3-D radiation transport code and, a multi-processor computer architecture. The code solves algorithms using a domain decomposition approach. For example, angular and spatial domains can be partitioned into subsets and, the subsets can be independently allocated and processed.

Description

REACTOR DOSIMETRY APPLICATIONS USING A PARALLEL 3-D RADIATION TRANSPORT CODE

CROSS-REFERENCE TO RELATED APPLICATIONS

This provisional patent application claims priority to U.S. Provisional Patent Application Number 60/975,525, which was filed on September 27, 2007.

BACKGROUND OF THE INVENTION

1. Field of the Invention

[0001] The present invention relates generally to the calculation of radiation field distributions, and is particularly useful in predicting neutron-dosimetry responses for nuclear reactor cavities and internal components.

2. Background Information

[0002] Various methodologies can be used to obtain numerical solutions of the

Linearized Boltzmann Equation (LBE) for neutron and gamma radiation transport applications. The discrete ordinates method (SN) is one such methodology used, in particular, in the nuclear engineering field. The numerical solution of the SN equations is achieved through the concurrent discretization of the phase space, i.e., angular, spatial and energy domains. The concurrent discretization of the phase space leads to a large number of unknowns in the SN equations and therefore, extensive computational resources are required to solve this problem.

[0003] For large 3-D neutron and gamma transport applications, the main memory required to generate a numerical solution of the LBE using SN equations may exceed current computational capabilities of a typical single-processor workstation. For example, the solution of a full 3-D neutron transport problem for a typical 2-loop Pressurized Water Reactor (PWR), characterized by approximately 1.5 million spatial meshes, an S₈ quadrature set, a P₃ expansion of the scattering kernel, and 47-neutron energy groups, can lead to a main memory requirement of approximately 45 GByte. The significant computational resources required may preclude the use of single-processor workstations to solve such problems. [0004] It would be desirable to overcome these difficulties, by developing a new solution algorithm(s) for the SN equations to take advantage of multi-processor computing architectures, i.e., distributed memory architectures. For example, it would be desirable to configure a number of physically independent workstations linked together via a network backbone, to establish what is generally referred to as a cluster computing environment. This type of computing platform has found widespread applications in recent years especially in the fields of scientific computing and large scale numerical simulation. However, it is necessary to devise specialized algorithms in order to exploit the capabilities of a cluster environment.

[0005] Thus, there is room for improvement in a set of solution algorithms for the

SN equations to take advantage of multi-processor computing architectures. There is also room for improvement in a methodology to obtain numerical solutions of the LBE for the calculation of radiation field distributions such as neutron and gamma radiation field distributions. Moreover, there is room for improvement of a methodology to predict dosimetry responses in an accurate and efficient manner for application in nuclear reactors.

SUMMARY OF THE INVENTION

[0006] A method for calculating a radiation field distribution, comprising applying a 3-D radiation transport computer code, said code comprising domain decomposition algorithms including domains selected from the group consisting of angular and spatial domains, wherein said domains are allocated and processed independently on a multi-processor computer architecture.

[0007] A computer program for calculating a radiation field distribution. The program including a code segment that when executed divides the angular and spatial domains into subsets, independently allocates the subsets, and processes the subsets on a multi-processor architecture. BRIEF DESCRIPTION OF THE DRAWINGS

[0008] Figure Ia shows the geometry and material distribution of a 3-D transport model for a 2-loop PWR. Figure Ib shows a 2-D section of the model on the x-y plane at z=0.0 cm for a 2-loop PWR.

[0009] Figure 2 shows the measured versus calculated (M/C) ratios of the dosimetry data calculated using the Directional Theta Weighted adaptive differencing scheme.

[0010] Figures 3a, 3b and 3c show the M/C ratios for a corrected reactor pressure vessel having a corrected thickness as compared to a non-corrected thickness for the core top location (3a), the core mid-plane location (3b) and the core bottom location (3c).

[0011] Figure 4 shows the speed-up obtained for a range of processors (e.g., up to

20) using different domain decomposition strategies.

[0012] Figure 5 shows a flow chart of an embodiment of the present invention, wherein the 3-D radiation transport computer program is implemented to generate a radiation field distribution.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0013] The present invention relates to a method of calculating radiation field distribution for a system. The radiation fields can include neutron and gamma radiation. In an aspect of the invention, ex-vessel neutron dosimetry responses in the cavity of a nuclear reactor can be calculated. The type of nuclear reactor is not limiting and can include a variety of commercial designs known in the art. Suitable reactors can include but are not limited to Pressurized Water Reactors (PWRs) and Boiling Water Reactors (BWRs). For simplicity of disclosure, this aspect of the invention will be described with reference to a 2-loop commercial PWR. Computer modeling of the PWR to generate dosimetry responses are used in the design and operation of the PWR. [0014] The method of the present invention includes application of a 3-D parallel radiation transport code referred to herein as RAPTOR-M3G (RApid Parallel Transport

Of Radiation - Multiple 3D Geometries). The transport code provides a set of parallel algorithms for solving the SN equations. The method is based on domain decomposition algorithms, where the spatial, angular and/or energy domains are partitioned into subsets which can be independently allocated and processed on a multi-processor architecture. Examples of suitable 3-D parallel deterministic transport codes which are known include PENTRAN™ and PARTISN (Sjoden G. E. and Haghighat A., "PENTRAN - Parallel Environment Neutralparticle TRANport in 3-D Cartesian Geometry," Proceedings of the Joint International Conference on Mathematical Methods and Supercomputing for Nuclear Applications, Vol. 1, pp. 232-234, Saratoga Springs, NY (1997)). As compared to traditional single-processor applications, the method of the present invention reduces the computational load as well as the memory requirement per processor, yielding an efficient solution methodology for large 3-D problems.

[0015] The RAPTOR-M3G computer code is developed in Fortran 90 using the

Message Passing Interface (MPI) parallel libraries. (Gropp W., Lusk E., and Skjellum A., Using MPI Portable Parallel Programming with the Message Passing Interface, The MIT Press, Cambridge, Massachusetts (1999)). Several features of the RAPTOR-M3G include the following:

• Solution of the multi-group S_N equations on 3-D Cartesian (RAPTOR- XYZ) and cylindrical geometries (RAPTOR-RTZ) (M. A. Hunter, G. Longoni, and S. L. Anderson, "Extension of RAPTOR-M3G to r-θ-z geometry for use in reactor dosimetry applications," Proceedings of the 13th International Symposium on Reactor Dosimetry, The Netherlands (2008)) on non-uniform orthogonal structured meshes;

• Spatial, angular, and coupled spatial/angular domain decomposition algorithms;

• Positive definite weighted differencing schemes: Zero/Theta Weighted, and Directional Theta Weighted;

• Automatic generation of level-symmetric quadrature sets up to order 20 (Longoni G. et al., "Investigation of New Quadrature Sets for Discrete Ordinates Method with Application to Non-Conventional Problems," Transactions of the American Nuclear Society, Vol. 84, pp. 224-226 (2001)); • Parallel memory: allows local allocation of the spatial and angular sub- domains, therefore reducing the memory requirements per processor;

• Parallel tasking: concurrent solution of the S_N equations on multiple processors to reduce computational time as compared to single processor technology;

• Parallel I/O: each processor locally accesses its storage devices to reduce I/O time; and

• Compatibility and integration with BOT3P (R. Orsi, "Potential Enhanced Performances in Radiation Transport Analysis on Structured Mesh Grids Made Available by BOT3P," Nuclear Science and Engineering, Vol. 157, pp. 110-116 (2007)), an automated mesh generator, and GIP, a multi-group cross-sections pre-processor.

[0016] Figure Ia shows the geometry of a 3-D transport model for an embodiment of the present invention, e.g., a 2-loop PWR. The PWR can include a 12- foot nuclear core, thermal shield design, and a 3-inch reactor cavity air gap. The model geometry includes a core-water mixture, core shroud, core barrel, thermal shield, Reactor Pressure Vessel (RPV) including stainless-steel liner, and reflective insulation. The RPV of a PWR typically has a generally cylindrical shape and is closed at both ends, e.g., by a bottom head and a removable top head. The upper and lower internals regions above and below the reactor core are modeled using a steel-water mixture. The lower internals of the RPV include the core barrel (i.e., a core support structure). The core barrel is enclosed with the thermal shield between the core barrel and an inner wall of the RPV. In some instances, neutron pads are used in lieu of the thermal shield. The core shroud sets inside the core barrel. An annular downcomer surrounds the reactor core barrel. Cooling fluid, typically water, is circulated into the downcomer.

[0017] Figure Ib shows a 2-D section of the model on the x-y plane at z=0.0 cm for a 2-loop PWR. Also shown is a material distribution in the PWR. The model geometry and the mesh discretization are generated using the BOT3P code, version 5.2. The model extends from 0.0 cm to 245.0 cm along the x-, and y-axis, and from -200.0 cm to 200.0 cm along the z-axis. A uniform mesh is applied throughout the model; a mesh size of 2.0 x 2.0 x 4.0 cm is specified along the x-, y- and z-axes, respectively, yielding a total of 1,464,100 meshes.

[0018] The cross sections for the material mixtures in the transport model are processed using the BUGLE-96 cross sections library (RSICC Data Library Collection BUGLE-96, "Coupled 47 Neutron, 20 Gamma-Ray Group Cross Section Library Derived from ENDF/B-VI for LWR Shielding and Pressure Vessel Dosimetry Applications," Oak Ridge National Laboratory, Oak Ridge, TN (1999)) and the GIP computer code, part of the DOORS package (RSICC Computer Code Collection DOORS 3.2a, "One-, Two- and Three-Dimensional Discrete Ordinates Neutron/Photon Transport Code System," Oak Ridge National Laboratory, Oak Ridge, TN (2003)). An Sg level symmetric quadrature set and a P₃ spherical harmonics expansion of the scattering kernel is used for the transport calculations. A system of passive neutron detectors can be installed in the reactor cavity air gap between the reflective insulation and the pressure vessel. The dosimetry system can provide accurate information relative to the fast neutron exposure over the beltline region of the reactor vessel. Pure metal foils can be installed in the reactor cavity, encased in an aluminum shell, which minimizes distortions of the fast neutron spectrum, effectively yielding a free-field measurement. The neutron dosimeters installed in the reactor cavity air gap are not explicitly defined in the transport model.

[0019] An aspect of the present invention includes a domain decomposition algorithm(s) for the discretization of the SN equations that uses an approach wherein the spatial and/or angular energy domains are partitioned into subsets which can be independently allocated and processed on multi-processor architectures.

[0020] The spatial and angular discretization of the S_N equations, as well as the angular domain decomposition algorithm described herein is specific for the 3D Cartesian XYZ version of the code. The formulation of the SN equations developed for RAPTOR-RTZ is different from RAPTOR-XYZ due to the presence of the scattering redistribution term.

[0021] The phase space of the SN equations is discretized, i.e., angle, space and energy. Thus, the resulting set of linear algebraic equations is suitable for solution on a digital computer. The energy domain is discretized using the multigroup approach into a number of discrete intervals, i.e., g=l ...G, starting with the highest energy particles (g=l), and ending with the lowest (g=G). The transport equation (i.e., LBE) in the multigroup approximation is formulated in Eq. (1).

(1) [0022] The angular domain is discretized by considering a finite set of directions and by applying an appropriate quadrature integration scheme. Each discrete direction can be visualized as a point on the surface of a unit sphere with an associated surface area which mathematically corresponds to the weight of the quadrature scheme. The combination of the discrete directions and the corresponding weights is referred to as quadrature set. In general, quadrature sets must satisfy a number of conditions in order to be accurate and mathematically determined; several approaches can be employed, e.g., level-symmetric quadrature set (LQn) and Legendre polynomial based quadrature sets (Longoni G. and Haghighat A., "Development of New Quadrature Sets with the Ordinate Splitting Technique," Proceedings of the ANS International Meeting on Mathematical Methods for Nuclear Applications (M&C 2001), Salt Lake City, UT, September 9-13, 2001, American Nuclear Society, Inc., La Grange Park, IL (2001)). The quadrature sets developed and implemented in RAPTOR-M3G are based on the LQn method. [0023] The spatial variable can be discretized with several techniques, e.g., finite difference and finite element methods. The formulation developed in RAPTOR-M3G is based on the finite difference approach which includes partitioning the spatial domain into computational cells, e.g., fine meshes, where the cross sections are assumed constant within each cell. In 3D Carstesian geometry, the angular flux at the cell-center location is evaluated using Eq. (2).

[0024] In Eq. (2), the angle and energy dependence are denoted by the indices m and g, respectively. The term q^_k represents the sum of the scattering, fission and external sources at cell-center. The indices i, j, k represent the cell-center values, and the weights ajj,k,m,g, bij,_k,_m,g, ^an^ ^cij,k,m,g are restricted to the range between 0.5 and 1.0; RAPTOR-M3G utilizes the Theta- Weighted (TW), Zero- Weighted (ZW), or the adaptive Directional Theta Weighted (DTW) differencing schemes to calculate the weights during the transport sweep (B. Petrovic and A. Haghighat, "New Directional Theta- Weighted S_N Differencing Scheme and Its Application to Pressure Vessel Fluence Calculations," Proceedings of the 1996 Radiation Protection and Shielding Topical Meeting, Falmouth, MA, Vol. 1, pp. 3-10 (1996)).

[0025] The S_N equations are solved by going through each direction starting from the boundary of the problem domain; this solution process is also referred to as transport sweep. The angular flux defined at center-cell locations is evaluated starting from boundary conditions or from the boundary angular flux previously calculated in adjacent cells. The cell-center angular flux is calculated using Eq. (2). The angular flux exiting the computational cell is calculated using additional relationships referred to as the "differencing schemes".

[0026] The transport sweep is performed within an iterative process which is termed source iteration, also known as fixed point iteration, or Richardson iteration. This process is continued until an appropriate convergence criterion is satisfied, i.e., the relative error on the scalar flux in any norm between two iterations is below a certain cutoff value (Adams M. L. and Larsen E. W., "Fast Iterative Methods for Discrete- Ordinates Particle Transport Calculations," Progress in Nuclear Energy, Vol. 40, n. 1

(2002)). For radiation shielding calculations this cutoff value is generally set to 1.Oe^'3 or 1.Oe^"4. [0027] The parallel algorithms developed in RAPTOR-M3G are based on the decomposition of the angular and/or spatial domains on a network of processors. RAPTOR-M3G creates a virtual topology based a number of processors allocated to the angular and spatial domains, specified as P_a and P_s respectively. The total number of processors required for any decomposition is P_n = P_a * P_s- Based on this information, the network of processors is mapped on the spatial and angular domains creating a virtual topology which associates each processor to its local sub-domain.

[0028] The angular domain is partitioned on an octant basis, where the processors specified on the angular domain, are sequentially assigned to the local octants. The local number of octants allocated per processor is given by Eq. (3).

^F° O)

[0029] The transport sweep is locally performed on Ni_oct octants on P_a processors; an MPI communicator for the angular domain is used to synchronize the angular flux among processors and to account for reflective boundary conditions.

[0030] The spatial domain is partitioned along the z-axis by sequentially assigning the P_s processors to a number of x-y planes. The total number of fine meshes along the z-axis, i.e., km, is partitioned on P_s processors; a mapping array, i.e., kmloc, is used to assign the x-y planes to the P_s processors. The number of x-y planes assigned to the P_s processors is arbitrary; however, the condition in Eq. (4) needs to be satisfied in order to define a spatial decomposition that is topologically consistent with the problem geometry.

[0031] The flexibility to map the processors on the spatial domain to an arbitrary number of x-y planes can depend on the fact that the number of z-planes may not be exactly divisible by the number of processors on the spatial domain. An uneven partitioning of the x-y planes on the P_s processors can lead to processor load imbalance with consequent loss in performance. In the present invention, a hybrid angular/spatial decomposition strategy may be applied to overcome this difficulty. The hybrid decomposition includes the combination of the angular and spatial domains to include concurrent partitioning of these domains. Hybrid domain decomposition is further described in the Example 1 below.

[0032] Figure 5 shows a flow chart of an embodiment of the present invention, wherein the 3-D radiation transport computer program is implemented to generate a radiation field distribution. This embodiment includes obtaining geometrical and material information on the system to be modeled. The information can be obtained from various sources such as but not limited to nuclear reactor drawings. An adequate S_N order required for the calculation can then be selected. Further, the P_n expansion order to the angular domain (i.e., scattering kernel and angular flux) can be selected. An appropriate differencing scheme (i.e., TW, ZW or DTW) can be selected as well. A 3-D model of the system can be generated and discretized using an appropriate mesh generator for the Cartesian XYZ or RTZ geometries. A cross section table is generated for each material by mixing them with an appropriate data set such as BUGLE-96. The number of processors and the corresponding decomposition strategy is selected for the problem to be solved. The calculation is then performed. Following calculation, postprocessing and analysis of the results generated may occur.

[0033] The RAPTOR-M3G computer code used in the present invention provides an accurate and efficient (e.g., reduced computational time) solution of radiation transport problems. In an aspect of the present invention, ex-vessel neutron dosimetry responses are calculated for a cavity in a nuclear reactor vessel. In comparison with actual measurements, the fast neutron reactions in the reactor cavity air gap of a 2-loop PWR calculated by RAPTOR-M3G were 96% accurate on average. Further, solution of the transport problem was obtained in approximately 106 minutes of clock time on a 20- processor computer cluster using a hybrid angular/spatial domain decomposition strategy.

[0034] While specific embodiments of the invention have been described in detail, it will be appreciated by those skilled in the art that various modifications and alternatives to those details could be developed in light of the overall teachings of the disclosure. For instance, an aspect of the present invention described herein relates to the nuclear industry and in particular, nuclear reactors. However, the present invention can also be used in a wide range of other applications such as the medical field. For example, the present invention can be used to determine the dose of radiation delivered to a patient for the treatment and/or cure of cancer. Accordingly, the particular arrangements disclosed are meant to be illustrative only and not limiting as to the scope of the invention which is to be given the full breadth of the appended claims and any and all equivalents thereof.

EXAMPLES

[0035] Example 1 - RAPTOR-M3G Parallel Performance Analysis

[0036] The transport calculations discussed in the Examples were performed with

RAPTOR-M3G running on a 20 processors computer cluster, i.e., EAGLE-I. The cluster was composed of 5 nodes with 2 dual-core dual processor AMD Opteron 64-bit architecture. The cluster total memory, i.e., RAM, available was 40 GByte; the network interconnection was characterized by 1 GBhVs bandwidth. With this hardware configuration, RAPTOR-M3G completed a full 3D transport calculation for a 2-loop PWR in approximately 106 minutes on 20 processors. No significant differences in performance were observed using the DTW, TW, or ZW differencing schemes.

[0037] Moreover, a simple test problem was set up to analyze the parallel performance of the code. The test problem consisted of 50 x 50 x 50 cm box with a uniform distributed fixed source, discretized with a 1 cm uniform mesh. An S 8 quadrature set and PO isotropic scattering was used, along with a one energy group cross section set. The wall-clock time, speed-up, and parallel efficiency were used to evaluate the parallel performance of RAPTOR-M3G. Speed-up and parallel efficiency were defined as in Eqs. 5 and 6, respectively. S_p = T_s / T_p (5) η_p = Sp / N_p (6) wherein, T_s and T_p are the wall-clock times required by the single-processor and multi-processor calculations, respectively. N_p is the number of processors utilized to achieve the wall-clock time T_p. Figure 4 shows a comparison of the speed-up obtained up to 20 processors using different decomposition strategies.

[0038] The speed-up obtained with the space decomposition was progressively reduced as the number of processors was increased. It was believed that this behavior was due to a finer computational granularity per processor; as the space domain was decomposed into smaller sub-domains, the number of operation per processor was reduced while the inter-processor communication time was increased; therefore leading to a reduced performance. The network data transfer among the nodes was generally the limiting factor on distributed memory architectures. The larger number of iterations required to converge the problem further contributed to reduce the performance of the space decomposition strategy. However, the hybrid decomposition where the angular and spatial domains were concurrently partitioned yielded better results. It was believed that this behavior was due to the coarser computational granularity induced by this decomposition; also for the hybrid decomposition, the number of iterations required to converge the problem did not increase as much as the space decomposition.

[0039] Example 2 - Comparison of Measured Dosimetry Responses With

RAPTOR-M3G Calculated Responses

[0040] A comparison was conducted between measured dosimetry responses and corresponding predictions obtained with RAPTOR-M3G. The IRDF-2002 Dosimetry Library (I. Kodeli and A. Trkov, "Validation of the IRDF-2002 Dosimetry Library," Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, Vol. 57, Issue 3, pp. 664-681 (2007)) was used to generate the calculated dosimetry responses for the neutron reactions listed in Table 1.

[0041] Measured dosimetry responses for the reactions listed in Table 1 were compared with the responses calculated by RAPT0R-M3G.

[0042] Table 1. Neutron Reactions Measured Bv The Dosimetry System

Material Reaction

Copper ⁶³Cu(n,α)⁶⁰Co

Iron ⁵⁴Fe(n,p)⁵⁴Mn

Nickel ⁵⁸Ni(n,p)⁵⁸Co

Uranium ²³⁸U(n,f)¹³⁷Cs

Neptunium ²³⁷Np(n,f)¹³⁷Cs

[0043] The reactions listed in Table 1 were measured using Cadmium shielded metal foils; therefore, the thermal component of the neutron spectrum was suppressed.

[0044] The measured responses were obtained at four azimuthal locations, i.e., 0°,

15°, 30°, and 45° at core mid-plane in the reactor cavity air gap. Since 2-loop PWR reactors are generally characterized by peak fast neutron fluence at the 0° position, due to close proximity of the nuclear fuel with the RPV, additional measurement were obtained at this location. Specifically at 0°, measurements were obtained axially at the top and bottom of the active core. Initially the calculated dosimetry responses were found to consistently over-predict the measured data.

[0045] Further investigation revealed that the RPV thickness used in the transport model was smaller than what was being measured during the In-Service Inspection (ISI) of the Reactor Pressure Vessel (RPV), confirming the initial findings. The new RPV thickness was introduced into the transport model, and the accuracy of the calculated dosimetry data improved on average by -8%. The measured vs. calculated (M/C) ratios of the dosimetry data calculated using the DTW adaptive differencing scheme are shown in Figure 2. As shown in Figure 2, the M/C ratios were consistently located within a 10% range at each location and for every dosimetry material. The over-prediction at the 30° and 45° azimuthal positions could be reduced by using a non-uniform mesh refinement at these locations, where the curvature of the system becomes more relevant. The average M/C ratio over all dosimetry locations was 0.96.

[0046] Figures 3a through 3c present the M/C values obtained with the ISI- corrected RPV thickness as compared to the value obtained without the thickness correction. The comparison was conducted and presented for all the dosimetry specimens located at the 0° azimuthal location. The corrected RPV thickness using the ISI measurements improved the accuracy of the calculated responses at every dosimetry location. Similar results were obtained also at 15°, 30°, and 45° azimuthal locations.

[0047] The measured and calculated reaction rates at each dosimetry location, and the M/C ratio using the RPV thickness corrected with ISI measurements are shown in Table 2. The average M/C ratio for the reactions listed in Table 2 across all the dosimetry locations was 0.96.

[0048] Table 2. Measured and Calculated Reaction Rates Obtained with the

DTW Differencing Scheme

Claims

What is claimed is:

1. A method for calculating a radiation field distribution, comprising: applying a 3-D radiation transport computer code, said code comprising domain decomposition algorithms including domains selected from the group consisting of angular and spatial domains, wherein said domains are allocated and processed on a multi-processor computer architecture.

2. The method of claim 1 wherein said radiation field distribution includes radiation fields selected from the group consisting of neutron radiation and gamma radiation.

3. The method of claim 1 wherein said radiation field distribution is for a cavity of a nuclear reactor.

4. The method of claim 3 wherein said nuclear reactor is a 2-loop pressurized water reactor.

5. The method of claim 1 wherein said algorithm is structured to numerically solve a Linearized Boltzmann Equation.

6. The method of claim 1 wherein said multi-processor computer architecture comprises a number of physically independent workstations linked together by a network connection.

7. The method of claim 1 further comprising partitioning said domains into subsets which are independently allocated and processed on said multi-processor architecture.

8. The method of claim 5 wherein said partitioning of said domains includes using parallel memory on multiple processors.

9. The method of claim 1 wherein said radiation field distribution includes ex- vessel neutron dosimetry responses in a nuclear reactor cavity.

10. The method of claim 7 wherein said dosimetry responses are obtained on a 3-D Cartesian geometry.

11. The method of claim 1 wherein said algorithms are structured to solve concurrently SN equations using parallel algorithms.

12. The method of claim 1 wherein each processor of said multi-processor computer architecture locally accesses its own devices selected from the group consisting of storage device and local memory.

13. The method of claim 1 further comprising discretizing said energy domain into a number of discrete intervals.

14. The method of claim 1 further comprising partitioning said spatial domain into computational cells.

15. The method of claim 14 wherein the assumed cross sections within each cell are constant.

16. The method of claim 1 wherein the domain decomposition is a hybrid decomposition wherein the angular and spatial domains are concurrently partitioned.

17. The method of claim 16 wherein the number of processors is up to and including 20.

18 The method of claim 17 wherein the problem is solved in a time period of less than 2 hours.

19. The method of claim 17 wherein a solution derived by using said method has an accuracy of 90% or greater when compared with calculated values.

20. A computer program for calculating a radiation field distribution, said program comprising a code segment that when executed divides the angular and spatial domains into subsets, independently allocates said subsets; and processes said subsets on a multiprocessor architecture.