|Program name||Package id||Status||Status date|
|Package ID||Orig. computer||Test computer|
|CCC-0759/05||Linux-based PC,PC Windows|
TITAN is a time-independent deterministic radiation transport simulation code in 3-D Cartesian geometry. The hybrid approach in the TITAN code allows different transport solvers (Sn or Ray-tracing) to be applied in different regions. TITAN solves both the k-effective and fixed source forward/adjoint problems. It has been benchmarked on a number of OECD/NEA benchmark problems:
3D Radiation Transport Benchmarks for Simple Geometries with Void Regions (Kobayashi problem).
Benchmark on Deterministic Transport Calculations Without Spatial Homogenisation MOX Fuel Assembly 3-D Extension Case (C5G7 MOX).
Benchmarking the Accuracy of Solution of 3-Dimensional Transport Codes and Methods over a Range in Parameter Space.
(A summary report can be found at http://18.104.22.168/report/).
The extra sweep with fictitious quadrature technique enables TITAN to simulate the SPECT (Single Photon Emission Computed Tomography) projection images. SPECT simulations generally are performed by Monte Carlo codes. A benchmark using the TITAN code on a torso phantom has been established. In the TITAN SPECT simulation, collimators are not explicitly simulated, instead, a technique, called Circular Ordinate Splitting, is used to simulate the collimator blurring effects.
TITAN numerically solves the time-independent first order transport equation (Linear Boltzmann Equation) using a hybrid Discrete Ordinate (Sn) and Ray-tracing method. Two transport solvers, an Sn Solver and a ray-tracing solver, are integrated in the TITAN code. Both solvers work on the coarse mesh level in Cartesian geometry. Generally, a TITAN problem model contains more than one coarse mesh. This allows users to apply different solvers to different coarse mesh. This feature can be useful for problems containing a large region of low scattering medium. In such regions, the Sn method requires finer angular and spatial meshing and becomes less efficient. TITAN's ray-tracing solver is more efficient to solve the transport equation in such regions. The ray-tracing solver is essentially a 3-D Method of Characteristics solver, only it applies to an individual coarse mesh, instead of the whole spatial domain. Currently the ray-solver applies only on coarse mesh with one material region, and the total cross-section of the material should be close to zero to qualify as 'low scattering' medium. For a multi-region regular coarse mesh, the Sn solver should be used.
TITAN uses the object oriented programming paradigm with a layered and modular structure. Some features of the code include:
Integrated SN and ray-tracing solvers
Shared scattering source kernel allowing arbitrary order anisotropic scattering
Block-oriented data structure allowing localized quadrature sets and solvers
Layered code structure
Level-symmetric and Pn-Tn quadrature sets
Incorporation of three ordinate splitting techniques (rectangular, local PN-TN, and circular)
Fast and memory-efficient spatial and angular projections on the interfaces of coarse meshes by using sparse projection matrix
wave-front interface flux handling
A binary I/O library to visualize and post-process data with TECPLOT
Extra Sweep technique with the fictitious quadrature technique for calculations of angular fluxes along arbitrary directions
What's new in Version 1.29:
Added the Lobatto-Chebyshev quadrature which includes a direction along an axis
Improved angular source handling
Improved parallel performance
(Refer to change.log in the package for more details)
|Package ID||Status date||Status|
Ce Yi and A. Haghighat:
"A Three-Dimensional Block-Oriented Hybrid Discrete Ordinates and Characteristics Method," Nucl. Sci. Eng., 164:3, 221-247 (March 2010).
|Package ID||Computer language|
|CCC-0759/05||FORTRAN 2003, FORTRAN-90, FORTRAN-95|
Both TITAN and PENMSHXP are written in FORTRAN 90/95/2003 language. The FORTRAN 2003 features in the codes are limited to some user-defined data objects and some I/O statements. Most up-to-date FORTRAN compliers should be able to compile both codes. Portland FORTRAN compiler (ver. 8.0 or higher), Intel FORTRAN compilers (ver. 10.0 or higher), and GCC/GFORTRAN (ver. 4.3) are tested in WINDOWS/LINUX 32/64 bit systems. The codes also can be compiled on MAC OS X with Inter FORTRAN for MAC or GFORTRAN/GCC installed. For the parallel TITAN version, an MPI2 implementation is required, such as MPICH2, Microsoft MPI, Intel MPI, or OpenMPI on WINDOWS/LINUX. Further instructions on how to compile the codes are included in the package.
For PENMSHXP v2.64b, it is the user's responsibility to ensure that access to DISLIN is in compliance with the makers of the DISLIN library (http://www.mps.mpg.de/dislin/). PENMSHXP also utilizes TECPLOT software if available for 3-D rendering with macro support http://www.tecplot.com/. A DISLIN-free version of PENMSHXP (ver. 2.61a) is also included. Ver 2.61a does not generate graphic output files. Instead, it generates a binary data file, which can be loaded with TECPLOT to render various plots.
Keywords: SN method, adjoint, complex geometry, deterministic radiation transport, discrete ordinates, ray-tracing, three-dimensional.