def _run_openmoc(self):
"""Plot the cells in the geometry."""
# Create a series of Matplotlib Figures / PIL Images for different
# plotting parameters and append to figures list
self.figures.append(
plot_cells(self.input_set.geometry, gridsize=100,
get_figure=True))
self.figures.append(
plot_cells(self.input_set.geometry, gridsize=100,
zcoord=10., get_figure=True))
self.figures.append(
plot_cells(self.input_set.geometry, gridsize=100,
get_figure=True, xlim=(0., 2.), ylim=(0., 2.)))
self.figures.append(
plot_cells(self.input_set.geometry, gridsize=100,
get_figure=True, library='pil'))
def _run_openmoc(self):
"""Plot the cells in the geometry."""
# Create a series of Matplotlib Figures / PIL Images for different
# plotting parameters and append to figures list
self.figures.append(
plot_cells(self.input_set.geometry, gridsize=100, get_figure=True))
self.figures.append(
plot_cells(self.input_set.geometry,
gridsize=100,
zcoord=10.,
get_figure=True))
self.figures.append(
plot_cells(self.input_set.geometry,
gridsize=100,
get_figure=True,
xlim=(0., 2.),
ylim=(0., 2.)))
self.figures.append(
plot_cells(self.input_set.geometry,
gridsize=100,
get_figure=True,
library='pil'))
track_generator.setSegmentFormation(openmoc.OTF_STACKS)
track_generator.generateTracks()
###############################################################################
########################### Running a Simulation ##########################
###############################################################################
solver = openmoc.CPUSolver(track_generator)
solver.setNumThreads(num_threads)
solver.setConvergenceThreshold(tolerance)
solver.computeEigenvalue(max_iters)
solver.printTimerReport()
###############################################################################
# Generating Plots
###############################################################################
log.py_printf('NORMAL', 'Plotting data...')
plotter.plot_tracks(track_generator, plot_3D=True)
plotter.plot_materials(geometry, gridsize=500, plane='xy', offset=0.)
plotter.plot_cells(geometry, gridsize=500, plane='xy', offset=0.)
plotter.plot_flat_source_regions(geometry, gridsize=500, plane='xy', offset=0.)
plotter.plot_spatial_fluxes(solver, energy_groups=[1,2,3,4,5,6,7], \
plane='xy', offset=0.)
plotter.plot_energy_fluxes(solver, fsrs=range(geometry.getNumFSRs()))
log.py_printf('TITLE', 'Finished')
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)
track_generator.setSegmentFormation(openmoc.OTF_STACKS)
track_generator.generateTracks()
###############################################################################
# Running a Simulation
###############################################################################
solver = openmoc.CPUSolver(track_generator)
solver.setNumThreads(num_threads)
solver.setConvergenceThreshold(tolerance)
solver.computeEigenvalue(max_iters)
solver.printTimerReport()
###############################################################################
# Generating Plots
###############################################################################
log.py_printf('NORMAL', 'Plotting data...')
plotter.plot_materials(geometry, gridsize=500)
plotter.plot_cells(geometry, gridsize=500)
plotter.plot_flat_source_regions(geometry, gridsize=500)
plotter.plot_spatial_fluxes(solver, energy_groups=[1, 2, 3, 4, 5, 6, 7])
log.py_printf('TITLE', 'Finished')
###############################################################################
######################## Creating the TrackGenerator ######################
###############################################################################
log.py_printf('NORMAL', 'Initializing the track generator...')
track_generator = TrackGenerator(geometry, num_azim, track_spacing)
track_generator.generateTracks()
###############################################################################
########################### Running a Simulation ##########################
###############################################################################
solver = CPUSolver(geometry, track_generator)
solver.setSourceConvergenceThreshold(tolerance)
solver.setNumThreads(num_threads)
solver.convergeSource(max_iters)
solver.printTimerReport()
###############################################################################
############################ Generating Plots #############################
###############################################################################
log.py_printf('NORMAL', 'Plotting data...')
plotter.plot_cells(geometry)
plotter.plot_materials(geometry)
plotter.plot_flat_source_regions(geometry)
plotter.plot_fluxes(geometry, solver, energy_groups=[1,2,3,4,5,6,7,8])
track_generator = TrackGenerator(geometry, num_azim, track_spacing)
track_generator.setNumThreads(num_threads)
track_generator.generateTracks()
###############################################################################
########################### Running a Simulation ##########################
###############################################################################
solver = CPUSolver(geometry, track_generator)
solver.setNumThreads(num_threads)
solver.setSourceConvergenceThreshold(tolerance)
solver.convergeSource(max_iters)
solver.printTimerReport()
###############################################################################
############################ Generating Plots #############################
###############################################################################
log.py_printf('NORMAL', 'Plotting data...')
#plotter.plot_tracks(track_generator)
#plotter.plot_segments(track_generator)
#plotter.plot_materials(geometry, gridsize=500)
plotter.plot_cells(geometry, gridsize=500)
#plotter.plot_flat_source_regions(geometry, gridsize=500)
#plotter.plot_fluxes(geometry, solver, energy_groups=[1,2,3,4,5,6,7])
log.py_printf('TITLE', 'Finished')
