Esempio n. 1
0
    def _run_openmoc(self):
        """Plot the flat source regions in the geometry."""

        # Run an eigenvalue calculation to setup FSR centroids
        super(PlotFSRsTestHarness, self)._run_openmoc()

        # Create a series of Matplotlib Figures / PIL Images for different
        # plotting parameters and append to figures list
        self.figures.append(
            plot_flat_source_regions(self.input_set.geometry,
                                     gridsize=100,
                                     get_figure=True))
        self.figures.append(
            plot_flat_source_regions(self.input_set.geometry,
                                     gridsize=100,
                                     get_figure=True,
                                     xlim=(0., 2.),
                                     ylim=(0., 2.)))
        self.figures.append(
            plot_flat_source_regions(self.input_set.geometry,
                                     gridsize=100,
                                     get_figure=True,
                                     centroids=True,
                                     marker_size=3))
        self.figures.append(
            plot_flat_source_regions(self.input_set.geometry,
                                     gridsize=100,
                                     get_figure=True,
                                     centroids=True,
                                     library='pil'))
Esempio n. 2
0
    def _run_openmoc(self):
        """Plot the flat source regions in the geometry."""

        # Run an eigenvalue calculation to setup FSR centroids
        super(PlotFSRsTestHarness, self)._run_openmoc()

        # Create a series of Matplotlib Figures / PIL Images for different
        # plotting parameters and append to figures list
        self.figures.append(
            plot_flat_source_regions(self.input_set.geometry, gridsize=100, 
                       get_figure=True))
        self.figures.append(
            plot_flat_source_regions(self.input_set.geometry, gridsize=100, 
                       get_figure=True, xlim=(0., 2.), ylim=(0., 2.)))
        self.figures.append(
            plot_flat_source_regions(self.input_set.geometry, gridsize=100, 
                       get_figure=True, centroids=True, marker_size=3))
        self.figures.append(
            plot_flat_source_regions(self.input_set.geometry, gridsize=100, 
                       get_figure=True, centroids=True, library='pil'))
Esempio n. 3
0
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')
Esempio n. 4
0
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')
Esempio n. 5
0
###############################################################################
########################   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])
Esempio n. 6
0
solver.setConvergenceThreshold(opts.tolerance)
solver.setNumThreads(opts.num_omp_threads)
solver.computeEigenvalue(opts.max_iters)
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,
Esempio n. 7
0
# plot geometry materials and cells
plotter.plot_cells(geometry, gridsize=1000)
plotter.plot_materials(geometry, gridsize=1000)

###############################################################################
########################   Creating the TrackGenerator   ######################
###############################################################################

log.py_printf('NORMAL', 'Initializing the track generator...')

track_generator = ModularTrackGenerator(geometry, num_azim, track_spacing)
track_generator.setLatticeStructure(25,25)
track_generator.generateTracks()

plotter.plot_flat_source_regions(geometry, gridsize=1000)

###############################################################################
###########################   Running a Simulation   ##########################
###############################################################################

solver = ModularCPUSolver(geometry, track_generator)
solver.setSourceConvergenceThreshold(tolerance)
solver.setNumThreads(num_threads)
solver.convergeSource(max_iters)
solver.printTimerReport()

###############################################################################
############################   Generating Plots   #############################
###############################################################################