From 6d5ba773fae665b7f59fd4e952f794288eb58127 Mon Sep 17 00:00:00 2001 From: Connor Moore Date: Fri, 10 Jul 2026 10:22:53 +0200 Subject: Upped particles per batch and added entropy tracking --- krusty.py | 14 +++++++++++--- 1 file changed, 11 insertions(+), 3 deletions(-) diff --git a/krusty.py b/krusty.py index 9e0ea06..a29dcb7 100755 --- a/krusty.py +++ b/krusty.py @@ -5,6 +5,7 @@ import argparse import hashlib import sys import time +import numpy as np # KRUSTY Model in OpenMC # C.Moore, @@ -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() -- cgit v1.2.3