from openmoc import *
import openmoc.log as log
import openmoc.plotter as plotter
import openmoc.materialize as materialize
from openmoc.options import Options
###############################################################################
####################### Main Simulation Parameters ########################
###############################################################################
options = Options()
num_threads = options.getNumThreads()
track_spacing = options.getTrackSpacing()
num_azim = options.getNumAzimAngles()
tolerance = options.getTolerance()
max_iters = options.getMaxIterations()
log.set_log_level("NORMAL")
log.py_printf("TITLE", "Simulating the OECD's C5G7 Benchmark Problem...")
###############################################################################
########################### Creating Materials ############################
###############################################################################
log.py_printf("NORMAL", "Importing materials data from HDF5...")
materials = materialize.materialize("../../c5g7-materials.h5")
import numpy
from openmoc import *
import openmoc.log as log
from openmoc.options import Options
###############################################################################
####################### Main Simulation Parameters ########################
###############################################################################
options = Options()
num_threads = options.getNumThreads()
track_spacing = options.getTrackSpacing()
num_azim = options.getNumAzimAngles()
tolerance = options.getTolerance()
max_iters = options.getMaxIterations()
log.set_log_level('NORMAL')
log.py_printf('TITLE', 'Simulating HW3 from Fall 2010 22.212...')
###############################################################################
########################### Creating Materials ############################
###############################################################################
log.py_printf('NORMAL', 'Creating materials...')
fuel = Material(1)
moderator = Material(2)
import openmoc
import openmoc.log as log
import openmoc.plotter as plotter
from openmoc.options import Options
from geometry import geometry
###############################################################################
####################### Main Simulation Parameters ########################
###############################################################################
options = Options()
num_threads = options.num_omp_threads
azim_spacing = options.azim_spacing
num_azim = options.num_azim
polar_spacing = options.polar_spacing
num_polar = options.num_polar
tolerance = options.tolerance
max_iters = options.max_iters
###############################################################################
######################## Creating the TrackGenerator ######################
###############################################################################
log.py_printf('NORMAL', 'Initializing the track generator...')
track_generator = openmoc.TrackGenerator3D(geometry, num_azim, num_polar,
azim_spacing, polar_spacing)
track_generator.setNumThreads(num_threads)
import openmoc
import openmoc.plotter as plotter
import openmoc.process as process
from openmoc.options import Options
from lattices import lattices, universes, cells
import numpy as np
refines = 1
###############################################################################
####################### Main Simulation Parameters ########################
###############################################################################
opts = Options()
###############################################################################
########################### Creating Lattices #############################
###############################################################################
c = universes['Core']
v = universes['Void']
a = universes['Control Rod']
r = universes['Reflector']
lattices['Root'].setWidth(width_x=5.0 / refines,
width_y=5.0 / refines,
width_z=5.0 / refines)
lattices['Root'].setUniverses(
[[np.repeat([r, r, r, r, r], refines).tolist()] * 4 * refines +
[np.repeat([r, r, r, a, r], refines).tolist()] * refines] * 2 * refines +
[[np.repeat([r, r, r, r, r], refines).tolist()] * 2 * refines +
from openmoc import *
import openmoc.log as log
import openmoc.plotter as plotter
import openmoc.materialize as materialize
from openmoc.options import Options
###############################################################################
####################### Main Simulation Parameters ########################
###############################################################################
options = Options()
num_threads = options.getNumThreads()
track_spacing = options.getTrackSpacing()
num_azim = options.getNumAzimAngles()
tolerance = options.getTolerance()
max_iters = options.getMaxIterations()
acceleration = options.getCmfdAcceleration()
relax_factor = options.getCmfdRelaxationFactor()
mesh_level = options.getCmfdMeshLevel()
log.set_log_level('NORMAL')
log.py_printf('TITLE', 'Simulating the OECD\'s C5G7 Benchmark Problem...')
###############################################################################
########################### Creating Materials ############################
###############################################################################
log.py_printf('NORMAL', 'Importing materials data from HDF5...')
#import numpy
#import matplotlib.pyplot as plt
from openmoc import *
import openmoc.log as log
import openmoc.plotter as plotter
import openmoc.materialize as materialize
import openmoc.process as process
from openmoc.options import Options
###############################################################################
####################### Main Simulation Parameters ########################
###############################################################################
options = Options()
num_threads = options.getNumThreads()
track_spacing = options.getTrackSpacing()
num_azim = options.getNumAzimAngles()
tolerance = options.getTolerance()
max_iters = options.getMaxIterations()
acceleration = options.getCmfdAcceleration()
relax_factor = options.getCmfdRelaxationFactor()
mesh_level = options.getCmfdMeshLevel()
log.set_log_level('NORMAL')
log.py_printf('TITLE', 'Simulating a Design-a-Critical Reactor Problem...')