diff options
Diffstat (limited to 'krusty.py')
| -rwxr-xr-x | krusty.py | 14 |
1 files changed, 11 insertions, 3 deletions
@@ -5,6 +5,7 @@ import argparse import hashlib import sys import time +import numpy as np # KRUSTY Model in OpenMC # C.Moore, <connor.moore@psi.ch> @@ -264,7 +265,7 @@ geo_obj.export_to_xml() ## Settings ## set_obj = openmc.Settings() -set_obj.particles = 20000 +set_obj.particles = 100000 set_obj.batches = 4500 set_obj.inactive = 500 set_obj.temperature = {"method": "interpolation"} @@ -272,12 +273,19 @@ set_obj.temperature = {"method": "interpolation"} set_src = openmc.IndependentSource() set_src.space = openmc.stats.Point((0.0,0.0,30.0)) set_src.space = openmc.stats.CylindricalIndependent( - r = openmc.stats.Uniform(geo_air_inner_cyl.r,geo_fuel_cyl.r), - z = openmc.stats.Uniform(geo_z_core_bot.z0,geo_z_core_top.z0), + r = openmc.stats.Uniform(0.0, geo_rad_cyl[-1].r), + z = openmc.stats.Uniform(geo_z_abs_bot.z0,geo_z_abs_top.z0), phi = openmc.stats.Uniform(0.0,2*math.pi), origin = (0.0,0.0,0.0) ) +set_entropy_mesh = openmc.RegularMesh() +set_entropy_mesh.lower_left = (-geo_rad_cyl[-1].r, -geo_rad_cyl[-1].r, geo_z_abs_bot.z0) +set_entropy_mesh.upper_right = (geo_rad_cyl[-1].r, geo_rad_cyl[-1].r, geo_z_abs_top.z0) +set_entropy_mesh.dimension = (30,30,30) + +set_obj.entropy_mesh=set_entropy_mesh + #set_obj.verbosity=10 set_obj.export_to_xml() |