crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-29 06:35:48 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Set cmfd tally type and estimator through C API

Browse source
This commit is contained in:
Shikhar Kumar 2018-08-16 20:12:04 -04:00
parent d80501ea17
commit 2f43fa3960
2 changed files with 11 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

3
examples/cmfd_testing/basic/capi/run_openmc_cmfd.py
View file

@ -7,15 +7,14 @@ cmfd_mesh.lower_left = [-10, -1, -1]
cmfd_mesh.upper_right = [10, 1, 1]
cmfd_mesh.dimension = [10, 1, 1]
cmfd_mesh.albedo = [0., 0., 1., 1., 1., 1.]
cmfd_mesh.energy = [0.001, 0.01, 0.1]
# Initialize CMFDRun object
cmfd_run = openmc.CMFDRun()
# Set all runtime parameters (cmfd_mesh, tolerances, tally_resets, etc)
# All error checking done under the hood when setter function called
cmfd_run.cmfd_mesh = cmfd_mesh
cmfd_run.cmfd_reset = [5,10]
# Run CMFD
cmfd_run.run()

12
openmc/cmfd.py
View file

@ -637,7 +637,7 @@ class CMFDRun(object):
self._cmfd_shift = 1.e-6
self._cmfd_spectral = 0.
self._cmfd_stol = 1.e-8
self._cmfd_reset = [5,10]
self._cmfd_reset = []
self._cmfd_write_matrices = False
# External variables used during runtime but users don't have control over
@ -899,7 +899,7 @@ class CMFDRun(object):
#call time_cmfdsolve % reset()
# Initialize all numpy arrays used for cmfd solver
self._allocate_cmfd(n_batches)
self._allocate_cmfd()
def _read_cmfd_input(self):
print(' Configuring CMFD parameters for simulation')
@ -1085,6 +1085,8 @@ class CMFDRun(object):
tally.filters = [mesh_filter]
# Set scores for tally
tally.scores = ['flux', 'total']
tally.type = 'volume'
tally.estimator = 'analog'
# Set attributes of CMFD neutron production tally
if i == 1:
@ -1095,6 +1097,8 @@ class CMFDRun(object):
tally.filters = [mesh_filter]
# Set scores for tally
tally.scores = ['nu-scatter', 'nu-fission']
tally.type = 'volume'
tally.estimator = 'analog'
# Set attributes of CMFD surface current tally
if i == 2:
@ -1105,6 +1109,8 @@ class CMFDRun(object):
tally.filters = [meshsurface_filter]
# Set scores for tally
tally.scores = ['current']
tally.type = 'mesh-surface'
tally.estimator = 'analog'
# Set attributes of CMFD P1 scatter tally
if i == 3:
@ -1115,6 +1121,8 @@ class CMFDRun(object):
tally.filters = [mesh_filter, legendre_filter]
# Set scores for tally
tally.scores = ['scatter']
tally.type = 'volume'
tally.estimator = 'analog'
# Set all tallies to be active from beginning
tally.active = True