Computer Programs

NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROGRAM 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 AUTHORS, 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 |
---|---|---|---|

MELODIE | NEA-1057/02 | Tested | 16-APR-1993 |

MELODIE | NEA-1057/03 | Report | 29-APR-1997 |

Machines used:

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

NEA-1057/02 | CRAY X-MP | CRAY X-MP |

NEA-1057/03 | UNIX W.S. |

[ top ]

3. DESCRIPTION OF PROGRAM OR FUNCTION

MELODIE assesses radiological consequences associated with a nuclear waste repository in a deep geological formation. MELODIE consists of three parts: MELO for deterministic evaluation; PREP and SPOP for sensitivity and uncertainty analysis:

MELO is a modular deterministic code composed mainly of a geosphere model adopted from METIS, developed at the Ecole Nationale Superieure des Mines de Paris; and of the ABRICOT biosphere model, developed at CEA/IPSN.

- The geosphere model is a 2D finite element code which first solves the groundwater flow equation (considered as a steady state), and then in a second step the transport (advection-diffusion-dispersion) equation for each radionuclide. The geosphere is taken into account as an equivalent porous medium with explicit description of the main fractures. The retardation of the radionuclides is modeled through the Kd concept, and nuclear decay chains are processed inside this model.

The boundary conditions represent:

. for the groundwater flow calculation, the hydraulic head or the groundwater flow rate;

. for the transport calculation, the activity released by the repository itself to the geosphere (using, for low soluble radionuclides, the solubility limit concept). Three cases may be taken into account for computing the value of this boundary condition: consequence of the degradation of the matrix (vitrified waste, ...); uniform release during a given period of time; waste already degradated at the initial time.

- The biosphere model converts the activity released at the water outlets to the corresponding individual dose through intake, inhalation, or external irradiation from the use of contaminated water in the main food chains. This model is compartimental and takes into account farming and dietetic considerations.

- PREP provides samples using the Latin Hypercube techniques for the distribution functions of the input variables, such as permeability, solubility limit value for radionuclides, retardation factors, etc.

- SPOP performs uncertainty and sensitivity analyses on the output of the MELO simulation runs. PREP and SPOP have been developed at the Joint Research Centre of the European Communities.

MELODIE assesses radiological consequences associated with a nuclear waste repository in a deep geological formation. MELODIE consists of three parts: MELO for deterministic evaluation; PREP and SPOP for sensitivity and uncertainty analysis:

MELO is a modular deterministic code composed mainly of a geosphere model adopted from METIS, developed at the Ecole Nationale Superieure des Mines de Paris; and of the ABRICOT biosphere model, developed at CEA/IPSN.

- The geosphere model is a 2D finite element code which first solves the groundwater flow equation (considered as a steady state), and then in a second step the transport (advection-diffusion-dispersion) equation for each radionuclide. The geosphere is taken into account as an equivalent porous medium with explicit description of the main fractures. The retardation of the radionuclides is modeled through the Kd concept, and nuclear decay chains are processed inside this model.

The boundary conditions represent:

. for the groundwater flow calculation, the hydraulic head or the groundwater flow rate;

. for the transport calculation, the activity released by the repository itself to the geosphere (using, for low soluble radionuclides, the solubility limit concept). Three cases may be taken into account for computing the value of this boundary condition: consequence of the degradation of the matrix (vitrified waste, ...); uniform release during a given period of time; waste already degradated at the initial time.

- The biosphere model converts the activity released at the water outlets to the corresponding individual dose through intake, inhalation, or external irradiation from the use of contaminated water in the main food chains. This model is compartimental and takes into account farming and dietetic considerations.

- PREP provides samples using the Latin Hypercube techniques for the distribution functions of the input variables, such as permeability, solubility limit value for radionuclides, retardation factors, etc.

- SPOP performs uncertainty and sensitivity analyses on the output of the MELO simulation runs. PREP and SPOP have been developed at the Joint Research Centre of the European Communities.

[ top ]

[ top ]

[ top ]

6. TYPICAL RUNNING TIME

The computation time of the two MELODIE test cases are the following:

- Test1 (MELODIE transport calculation): about 32 seconds;

- Test2 (MELODIE steady-state flow sensitivity calculation): about 22 seconds.

The PREP test problem setup for 100 runs and 14 variables requires 0.44 seconds.

The SPOP test problem setup for 100 runs and 14 variables requires 10.8 seconds.

NEA 1057/02: NEA-DB executed the sample cases of the programs included in this package on a CRAY-XMS. The followingt CPU times were required: PREP: 2sec.; MELODIE case 1: 4min; MELODIE case 2: 1min; SPOP: 27sec.

The computation time of the two MELODIE test cases are the following:

- Test1 (MELODIE transport calculation): about 32 seconds;

- Test2 (MELODIE steady-state flow sensitivity calculation): about 22 seconds.

The PREP test problem setup for 100 runs and 14 variables requires 0.44 seconds.

The SPOP test problem setup for 100 runs and 14 variables requires 10.8 seconds.

NEA 1057/02: NEA-DB executed the sample cases of the programs included in this package on a CRAY-XMS. The followingt CPU times were required: PREP: 2sec.; MELODIE case 1: 4min; MELODIE case 2: 1min; SPOP: 27sec.

[ top ]

[ top ]

[ top ]

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

NEA-1057/02 | 16-APR-1993 | Tested at NEADB |

NEA-1057/03 | 29-APR-1997 | Report Only |

[ top ]

NEA-1057/02, included references:

- Marc Bourgeois:MODELE MELODIE - Evaluation de l'Impact Radiologioque Post

Fermeture d'un Stockage de Dechets Radioactifs en Formation

Geologique Profonde

Volume I - Notice Theorique

Not/The No.0/48 Rev. A (Fevrier 1, 1991).

Volume II - Techniques de Programmation

Tec/Pro No.0/45 Rev. B (Fevrier 1, 1991).

Volume III - Notice d'Utilisation

Not/Uti No.0/99 Rev. A (Fevrier 12, 1991).

NEA-1057/03, included references:

- M. Bourgeois and C. Serres:Modelisation du Terme Source dans le Code de Calcul MELODIE

DPRE/SERGD Rapport technique 97/06

[ top ]

[ top ]

NEA-1057/02

UNICOS 6.1 with compiler CFT77 5.0.1.0 (CRAY XMS).[ top ]

[ top ]

[ top ]

NEA-1057/03

report DPRE/SERGD 97/06 REPPTNEA-1057/02

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

NEA1057_02.001 | Information file | 76 |

NEA1057_02.002 | JCL to run PREP, MELODIE and SPOP programs | 14 |

NEA1057_02.003 | PREP FORTRAN 77 source program | 1299 |

NEA1057_02.004 | PREP sample input test case | 85 |

NEA1057_02.005 | MELODIE FORTRAN 77 source program | 217 |

NEA1057_02.006 | Library FORTRAN 77 source subroutines | 15462 |

NEA1057_02.007 | CRAY emulation of IMSL routines | 1558 |

NEA1057_02.008 | MELODIE sample input test case 1 | 2055 |

NEA1057_02.009 | MELODIE sample input test case 2 | 2190 |

NEA1057_02.010 | SPOP FORTRAN 77 source program | 2102 |

NEA1057_02.011 | SPOP sample input test case | 312 |

NEA1057_02.012 | MELODIE sample output test case 1 | 635 |

NEA1057_02.013 | MELODIE test case 1 unit fort.16 | 4256 |

NEA1057_02.014 | MELODIE test case 1 unit fort.31 | 955 |

NEA1057_02.015 | MELODIE test case 1 unit fort.32 | 212 |

NEA1057_02.016 | MELODIE test case 1 unit fort.33 | 841 |

NEA1057_02.017 | MELODIE test case 1 unit fort.41 | 188 |

NEA1057_02.018 | MELODIE test case 1 unit fort.42 | 188 |

NEA1057_02.019 | MELODIE test case 1 unit fort.53 | 188 |

NEA1057_02.020 | MELODIE test case 1 unit fort.61 | 188 |

NEA1057_02.021 | MELODIE test case 1 unit fort.62 | 188 |

NEA1057_02.022 | MELODIE test case 1 unit fort.63 | 96 |

NEA1057_02.023 | MELODIE sample output test case 2 | 5036 |

NEA1057_02.024 | MELODIE output test case 2 unit FORT.13 | 1750 |

NEA1057_02.025 | MELODIE output test case 2 unit FORT.14 | 12054 |

NEA1057_02.026 | MELODIE output test case 2 unit FORT.15 | 13300 |

NEA1057_02.027 | MELODIE output test case 2 unit FORT.40 | 304 |

NEA1057_02.028 | PREP sample output test case | 506 |

NEA1057_02.029 | PREP sample output unit FORT.50 | 300 |

NEA1057_02.030 | SPOP sample output test case | 1987 |

NEA1057_02.031 | SPOP sample output test case unit FORT.33 | 10 |

NEA1057_02.032 | SPOP sample output test case unit FORT.41 | 105 |

NEA1057_02.033 | SPOP sample output test case unit FORT.42 | 110 |

NEA1057_02.034 | SPOP sample output test case unit FORT.52 | 29 |

NEA1057_02.035 | SPOP sample output test case unit FORT.53 | 74 |

NEA1057_02.036 | SPOP sample output test case unit FORT.54 | 74 |

Keywords: convection, diffusion, dispersions, finite elements, ground water, radioactive waste storage, radionuclide migration, sensitivity analysis, solubility, statistics, systems analysis, two-dimensional.