Computer Programs

NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROBLEM OR FUNCTION, METHOD OF SOLUTION, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, UNUSUAL FEATURES OF THE PROGRAM, RELATED AND AUXILIARY PROGRAMS, STATUS, REFERENCES, MACHINE REQUIREMENTS, LANGUAGE, OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED, OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS, NAME AND ESTABLISHMENT OF AUTHOR, MATERIAL, CATEGORIES

[ top ]

[ top ]

To submit a request, click below on the link of the version you wish to order. Rules for end-users are
available here.

Program name | Package id | Status | Status date |
---|---|---|---|

TRIDENT | CCC-0293/01 | Tested | 08-APR-1982 |

TRIDENT | CCC-0293/02 | Tested | 01-JUL-1978 |

Machines used:

Package ID | Orig. computer | Test computer |
---|---|---|

CCC-0293/01 | CDC 7600 | CDC 7600 |

CCC-0293/02 | IBM 370 series | IBM 370 series |

[ top ]

3. DESCRIPTION OF PROBLEM OR FUNCTION

TRIDENT solves the two- dimensional-multigroup-transport equations in rectangular (x-y) and cylindrical (r-z) geometries using a regular triangular mesh.

Regular and adjoint, inhomogeneous and homogeneous (K-eff and eigenvalue searches) problems subject to vacuum, reflective, white, or source boundary conditions are solved. General anisotropic scattering is allowed and anisotropic-distributed sources are permitted.

TRIDENT solves the two- dimensional-multigroup-transport equations in rectangular (x-y) and cylindrical (r-z) geometries using a regular triangular mesh.

Regular and adjoint, inhomogeneous and homogeneous (K-eff and eigenvalue searches) problems subject to vacuum, reflective, white, or source boundary conditions are solved. General anisotropic scattering is allowed and anisotropic-distributed sources are permitted.

[ top ]

4. METHOD OF SOLUTION

The discrete-ordinates approximation is used for the neutron directional variables. An option is included to append a fictitious source to the discrete-ordinates equations that is defined such that spherical-harmonics solutions (in x-y geometry) or spherical-harmonics-like solutions (in r-z geometry) are obtained. This option is useful for cases in which ray-effect distortions are severe. A spatial-finite-element method is used in which the angular flux is expressed as a linear polynomial in each triangle that is discontinuous at triangle boundaries. Both inner (within-group) and outer iteration cycles are accelerated by either whole-system or fine-mesh rebalance.

Provision is made for creation of standard interface output files for Sn constants, angle-integrated (scalar) fluxes, and angular fluxes. Standard interface input files for Sn constants, inhomogeneous sources, cross sections, and the scalar flux may be read. Subroutines DRED and DRIT perform information transfers between LCM and random disk. Data transfer between large and small core are performed by CRED and CRIT. Sequential binary operations are performed by SEEK, REED, and RITE. Flexible edit options as well as a dump and restart capability are provided.

The discrete-ordinates approximation is used for the neutron directional variables. An option is included to append a fictitious source to the discrete-ordinates equations that is defined such that spherical-harmonics solutions (in x-y geometry) or spherical-harmonics-like solutions (in r-z geometry) are obtained. This option is useful for cases in which ray-effect distortions are severe. A spatial-finite-element method is used in which the angular flux is expressed as a linear polynomial in each triangle that is discontinuous at triangle boundaries. Both inner (within-group) and outer iteration cycles are accelerated by either whole-system or fine-mesh rebalance.

Provision is made for creation of standard interface output files for Sn constants, angle-integrated (scalar) fluxes, and angular fluxes. Standard interface input files for Sn constants, inhomogeneous sources, cross sections, and the scalar flux may be read. Subroutines DRED and DRIT perform information transfers between LCM and random disk. Data transfer between large and small core are performed by CRED and CRIT. Sequential binary operations are performed by SEEK, REED, and RITE. Flexible edit options as well as a dump and restart capability are provided.

[ top ]

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

Variable

dimensioning is used so that any combination of problem parameters leading to a container array less than MAXLEN can be accommodated. On CDC Machines MAXLEN can be about 40 000 words of Small Core Memory (SCM) and Large Core Memory (LCM) is used for most group- dependent data. On IBM machines, TRIDENT executes in double precision (8 bytes per Word) so that MAXLEN can be on the order of 100 000 or more. Although most problems can be fast core contained on the IBM machine and SCM and LCM contained on the CDC machine, an option exists to automatically overflow the storage of selected large arrays to random disk if necessary.

Variable

dimensioning is used so that any combination of problem parameters leading to a container array less than MAXLEN can be accommodated. On CDC Machines MAXLEN can be about 40 000 words of Small Core Memory (SCM) and Large Core Memory (LCM) is used for most group- dependent data. On IBM machines, TRIDENT executes in double precision (8 bytes per Word) so that MAXLEN can be on the order of 100 000 or more. Although most problems can be fast core contained on the IBM machine and SCM and LCM contained on the CDC machine, an option exists to automatically overflow the storage of selected large arrays to random disk if necessary.

[ top ]

6. TYPICAL RUNNING TIME

A six group, S2, 1700 triangle, K-eff calculation of an EBR-II core modeled in x-y geometry requires approximately 4.4 min of CDC-7600 time. This is approximately the same as the TRIPLET running time. Running times vary almost linearly with the total number of unknowns. When the ray-effect correction option is used, the running time is increased by a factor of approximately two, depending upon the problem.

A six group, S2, 1700 triangle, K-eff calculation of an EBR-II core modeled in x-y geometry requires approximately 4.4 min of CDC-7600 time. This is approximately the same as the TRIPLET running time. Running times vary almost linearly with the total number of unknowns. When the ray-effect correction option is used, the running time is increased by a factor of approximately two, depending upon the problem.

CCC-0293/01

NEA-DB executed the test case included in the package on CDC 7600 in 192 seconds of CPU time.[ top ]

[ top ]

[ top ]

Package ID | Status date | Status |
---|---|---|

CCC-0293/01 | 08-APR-1982 | Tested at NEADB |

CCC-0293/02 | 01-JUL-1978 | Tested at NEADB |

[ top ]

10. REFERENCES

- R. Douglas O'Dell and Raymond E. Alcouffe:

Transport Calculations for Nuclear Analyses: Theory and

Guidelines for Effective Use of Transport Codes

LA-10983-MS and UC-32 (September 1987)

- R. Douglas O'Dell and Raymond E. Alcouffe:

Transport Calculations for Nuclear Analyses: Theory and

Guidelines for Effective Use of Transport Codes

LA-10983-MS and UC-32 (September 1987)

CCC-0293/01, included references:

- T.J. Seed, W.F. Miller, F.W. Brinkley:"TRIDENT; A Two-Dimensional, Multigroup, Triangular Mesh

Discrete Ordinates, Explicit Neutron Transport Code"

LA-6735-M, (March 1977)

- Daniel Carson:

"Core to Core Conversion by Daniel Carson"

ANL-K250S-1 (Revision, February 1970)

CCC-0293/02, included references:

- T.J. Seed, W.F. Miller, F.W. Brinkley:"TRIDENT; A Two-Dimensional, Multigroup, Triangular Mesh

Discrete Ordinates, Explicit Neutron Transport Code"

LA-6735-M, (March 1977)

- Daniel Carson:

"Core to Core Conversion by Daniel Carson"

ANL-K250S-1 (Revision, February 1970)

[ top ]

11. MACHINE REQUIREMENTS

A maximum of three interface units, open at once, and three system input/output units are required. One random- disk unit is required if the automatic-data-overflow option is used. A large bulk memory is assumed by the code such as the CDC-7600 LCM. This function may, however, be replaced by a portion of fast-core memory with IBM 360 machines or by disk, drum, or tape storage.

A maximum of three interface units, open at once, and three system input/output units are required. One random- disk unit is required if the automatic-data-overflow option is used. A large bulk memory is assumed by the code such as the CDC-7600 LCM. This function may, however, be replaced by a portion of fast-core memory with IBM 360 machines or by disk, drum, or tape storage.

CCC-0293/01

To run the test case on CDC 7600, main storage require- ments are: 62000 octal words of SCM; 241000 octal words of LCM.[ top ]

Package ID | Computer language |
---|---|

CCC-0293/01 | FORTRAN-IV |

CCC-0293/02 | FORTRAN+ASSEMBLER |

[ top ]

CCC-0293/01

7000 SCOPE 2.1 LVL 533-A0.[ top ]

[ top ]

[ top ]

CCC-0293/01

File name | File description | Records |
---|---|---|

CCC0293_01.003 | TRIDENT INFORMATION FILE | 49 |

CCC0293_01.004 | TRIDENT SOURCE PROGRAM (FORTRAN-4) | 11561 |

CCC0293_01.005 | SYSTEM SUBROUTINES (FORTRAN-4) | 315 |

CCC0293_01.006 | ECRDWR SUBROUTINE | 12 |

CCC0293_01.007 | TRIDENT TEST CASE INPUT DATA | 571 |

CCC0293_01.008 | TRIDENT JCL | 20 |

CCC0293_01.009 | TRIDENT TEST CASE PRINTED OUTPUT | 13815 |

CCC-0293/02

File name | File description | Records |
---|---|---|

CCC0293_02.001 | SOURCE PROGRAM (F4,EBCDIC) | 11972 |

CCC0293_02.002 | SOURCE PROGRAM (ASSEMBLER,EBCDIC) | 694 |

CCC0293_02.003 | DD CARDS FOR THE EXECUTION | 28 |

CCC0293_02.004 | OVERLAY CARDS | 30 |

CCC0293_02.005 | SAMPLE PROBLEM INPUT DATA | 571 |

CCC0293_02.006 | SAMPLE PROBLEM PRINTED OUTPUT | 17445 |

Keywords: anisotropic scattering, discrete ordinate method, multigroup, r-z, transport theory, x-y.