def _run_openmoc(self): """Run OpenMOC and plot the spatial fluxes in the geometry.""" # Run an eigenvalue calculation super(PlotSpatialFluxesTestHarness, self)._run_openmoc() # Specify energy groups for which to plot the spatial flux energy_groups = [1, 3, 5, 7] # Create a series of Matplotlib Figures / PIL Images for different # plotting parameters and append to figures list self.figures.extend( plot_spatial_fluxes(self.solver, gridsize=100, get_figure=True, energy_groups=energy_groups)) self.figures.extend( plot_spatial_fluxes(self.solver, gridsize=100, get_figure=True, xlim=(0., 2.), ylim=(0., 2.), energy_groups=energy_groups)) self.figures.extend( plot_spatial_fluxes(self.solver, gridsize=100, get_figure=True, energy_groups=energy_groups, 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')
solver.printTimerReport() ############################################################################### ############################ Generating Plots ############################# ############################################################################### process.compute_fission_rates(solver, use_hdf5=False) plotter.plot_materials(geometry, gridsize=500, plane='xy') plotter.plot_materials(geometry, gridsize=500, plane='xz', offset=-10.0) plotter.plot_materials(geometry, gridsize=500, plane='yz') plotter.plot_cells(geometry, gridsize=500) plotter.plot_flat_source_regions(geometry, gridsize=500, plane='xy') plotter.plot_flat_source_regions(geometry, gridsize=500, plane='xz') plotter.plot_flat_source_regions(geometry, gridsize=500, plane='yz') plotter.plot_spatial_fluxes(solver, energy_groups=[1, 2], gridsize=500, plane='xy', offset=0.) plotter.plot_spatial_fluxes(solver, energy_groups=[1, 2], gridsize=500, plane='xz', offset=0.) plotter.plot_spatial_fluxes(solver, energy_groups=[1, 2], gridsize=500, plane='yz', offset=0.)