ML23093A160
ML23093A160 | |
Person / Time | |
---|---|
Issue date: | 04/14/2023 |
From: | Hathaway A NRC/RES/DSA |
To: | |
References | |
Download: ML23093A160 (1) | |
Text
MelMACCS - Converting MELCOR Output to MACCS Input Trey Hathaway, PhD U.S. Nuclear Regulatory Commission Office of Nuclear Regulatory Research Accident Analysis Branch Alfred.Hathaway@nrc.gov 14th European MELCOR User Group Meeting Ljubljana, Slovenia 12-14 April, 2023
MelMACCS
- Post-processing MELCOR accident progression results for MACCS consequence analysis calculations
- MACCS considers plume segments
- e.g., release fractions for chemical groups per plume segments
- MACCS typically considers 9 chemical classes
- Equivalent to first 9 MELCOR chemical classes
- MelMACCS combines CsI (Class 16) and CsMoO4 (Class 17) into Cs (Class 2), I (Class 4), and Mo (Class 7) using mass fractions
- This comparison is based on MelMACCS 2.0.1 NOTE: This presentation will attempt to focus on information passed between MELCOR analysis to the MACCS analysis 2
Radionuclide Inventory
- SCALE-ORIGEN used to generate radionuclide mass, activity, and decay heat
- MELCOR needs radionuclide mass and decay heat
- Radionuclide results are combined into chemical classes
- MACCS needs radionuclide mass and activity results
- Radionuclides are mapped to the chemical classes
- MelMACCS Inventory file
- Used to generate activities for MACCS
- Inventory files have three groups of radionuclides
- FISSION, ACTIVATION, ACTINIDE
- Provide radionuclide masses and activities for each group
- Sum of radionuclide masses in the inventory file within a chemical class must equal the MELCOR class mass
- MelMACCS scales results to match MELCOR class mass 3
Release Flowpath Characteristics
- User defines characteristics of each flowpath to the environment
- Building Height
- (SigmaY-Initial)
- (SigmaZ-Initial)
- Flowpaths are defined by the MELCOR analyst
- Must confirm building characteristics at release point 4
Ground Height
- MACCS references elevation relative to ground level MELCOR Flowpath elevation MACCS Flowpath elevation
- MELCOR has a relative 0.0 m elevation defined by the analyst (e.g.
bottom of the inner surface of RPV set to 0.0 m)
Ground
- Define a correction Level factor to adjust MELCOR flowpath elevations to ground level MELCOR 0.0 m elevation MACCS Reference Height 5
and
- Defines the initial plume dimensions in the crosswind ( ) and vertical dimensions ( ) for each plume released
- Plume dimensions at a downwind distance of zero
- Estimates dispersion from building wake effects at the initial release
- Serve as boundary conditions for calculating virtual source location
- =
4.3
- is the building width at the release point
- =
2.15
- is the building height at the release point 6
FL_MACCS FL_MACCS MACCSN NFL MACCSname MACCSnum MCCSFP DIRFL DF
- MELCOR outputs the cumulative mass released through a flow path for each particle size bin in kg
- MACCS-R-M-RE-C.N
- R - flow path number
- C - Class number
- N - particle size bin number
- Particle size bin 0 reflects vapor phase
- NOTE: Can be used to force MELCOR to write data to the plot file
- MACCS expects flow to the environment
- DF parameter on FL_MACCS card controls if output is MACCS flow path (i.e. to environment)
- Default is MACCS (i.e. flow to the environment)
- If not defined appropriately, MelMACCS can double count the radionuclide release leading to an error in the MACCS input
- Additional output: SCRAM time, Number of particle bins and size of the particles, Chemical groups and initial mass of chemical groups, Flow path release height, fluid temperature, sensible heat, fluid flow rate, molecular weight 7
Plume Segments
- MACCS performs consequence calculations using plume segments of chemical groups
- Plume Segments have a start time (RDPDELAY) and a duration (RDPLUDUR)
- A reference time corrects for differences between SCRAM time and accident initiation
- + and + + is equivalent to MELCOR time
- Plume segments define information related to the source term and the location of the release
- Location of release: Building height (WEBUILDH), Building width (WEBUILDW), Building length (WEBUILDL), Incident wind angle (WEBUILDA), Trapped plume release height (RDPHTRAP), Initial Sigma-Y width (SIGYINIT), Initial Sigma-Z width (SIGZINIT), Height of release (RDPLHITE)
- Source Term: Particle size distribution of chemical classes, Start of plume segment release (RDPDELAY), Plume segment duration (RDPLUDUR), Release of sensible heat (RDPLHEAT),
Mass flow rate (RDPLMFLA), Gas density (RDPLMDEN), Release fraction of chemical classes (RDRELFRC) 8
Release Fraction (RDRELFRC)
- RDRELFRC - mass of a chemical class released during a plume segment through a flowpath divided by the initial chemical class mass
- MELCOR Outputs used: MACCS-R-M-RE-C.N MACCS-INITIAL-MASS(C)
- R is the flow path, C is the chemical class, and N is the particle size bin (0-10)
- Particle size bin 0 represents vapor
- Sum of MACCS-R-M-RE-C.N for a flow path and chemical class reflects the total chemical class mass released through a flow path
- Dividing the mass released during the plume segment (i.e. released for PLUDUR at PDELAY) yields the release fraction for that plume segment 9
Release Fraction Comparison 10
Particle Size Distribution (RDPSDIST)
- RDPSDIST - Fraction of a chemical group associated with a particle size bin
- MELCOR Outputs used:
MACCS-R-M-RE-C.N
- The total mass released in a particle size bin of a chemical class for all flow paths divided by the total mass released by a class
- N is for particle bins 1 through 10
- Neglects vapor
- , =
11
Deposition Velocity
- Two options available
- Gravitational Settling
- Dependent on aerosol size
- Independent of wind conditions or surface roughness
- Expert Elicitation Data
- Wind speed and surface roughness considered as variables
- Two wind speeds (2 and 5 m/s) and three surface types with representative surface roughness
- Choose representative wind speed, surface roughness, and quantile of distribution
- If aerodynamic particle diameter
( ) is greater than cutoff diameter the calculation is switched to gravitational settling
- Output value is larger value of the value calculated using gravitational settling versus the value calculated using expert elicitation 12
Sensible Heat (RDPLHEAT)
- RDPLHEAT - Rate of sensible heat released
- Used to determine buoyancy flux using the Power Model
- Plume rise occurs only when wind speed at release of plume segment is less than critical wind speed calculated using sensible heat
- MELCOR Outputs Used:
MACCS-R-PLHEAT.0
,2 ,1
- =
2 1
- , 2 = . 0 2
- 2 = + +
- , 1 = . 0 1
- 1 = +
13
Average Plume Flow Rate (RDPLMFLA)
- RDPLMFLA - Average flow rate for plume segment in units of kg/s
- Used to calculate buoyancy flux using Density and Flow model
- MELCOR Outputs used:
MACCS-R-PLMFLO.0 MACCS-R-PLMWT.0
- Converts MELCOR fluid mass flow rate in moles to kg by multiplying by average molecular weight
- = 1
+1 2
- = . 0
- 1 = . 0 1
- = . 0
- 1 = . 0 1 1
- =
2 1 =1
- 2 = + +
- 1 = +
14
Plume Gas Density (RDPLMDEN)
- RDPLMDEN - Average gas density over all time steps in a plume segment
- Used to calculate buoyancy flux using Density and Flow model
- MELCOR Outputs used: MACCS-R-PLMFLO.0 MACCS-R-PLMWT.0 MACCS-R-PLTEMP.0
.0
- =
.0 1
- =
2 1 1
+ 1 1 2
- PLMDEN=
- = 1.012 3
- = 8.314105
- Ideal Gas Equation: pV = n
=
15
Activity Calculation
- MACCS code converts plume segment release fractions of chemical classes to released isotopic activity
- MelMACCS inventory defines radionuclide mass (kg) and activity (Ci)
- Radionuclides are mapped to chemical classes
- Sum of the radionuclide masses in a chemical class is equal to ELMMAS of the MELCOR DCH_EL card
=
16
Activity Calculation
- MACCS prints the radionuclide activity per plume segment as a data block in the output file
- Data block heading:
RELEASED INVENTORY (Ci)
OF ALL PLUMES heading
- The end of the plume segment is equal to + +
- Sum plume segments with same ending time to determine the activity released at that time from that flowpath 17
Activity Calculation 18