Exemplo n.º 1
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])
Exemplo n.º 2
0
###############################################################################
########################   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   #############################
###############################################################################

log.py_printf('NORMAL', 'Plotting data...')

plotter.plot_fluxes(geometry, solver, energy_groups=[1,2,3,4,5,6,7,8], gridsize=1000)
try:
    import openmoc.cuda as cuda

    # Check whether or not the machine contains a GPU
    if cuda.machineContainsGPU():
        solver = cuda.GPUSolver(geometry, track_generator)
    else:
        solver = ThreadPrivateSolver(geometry, track_generator)
        solver.setNumThreads(num_threads)

# If the openmoc.cuda module has not been built, use the standard solver
except:
    solver = ThreadPrivateSolver(geometry, track_generator, cmfd)
    solver.setNumThreads(num_threads)

solver.setSourceConvergenceThreshold(tolerance)
solver.convergeSource(max_iters)
solver.printTimerReport()


###############################################################################
##############################   Plotting Data   ##############################
###############################################################################


log.py_printf('NORMAL', 'Plotting data...')

plotter.plot_fluxes(geometry, solver, energy_groups=[1,7], gridsize=400)

log.py_printf('TITLE', 'Finished')
Exemplo n.º 4
0
log.py_printf('NORMAL', 'Initializing the track generator...')

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.setSourceConvergenceThreshold(tolerance)
solver.setNumThreads(num_threads)
solver.convergeSource(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_fluxes(geometry, solver, energy_groups=[1,2,3,4,5,6,7])

log.py_printf('TITLE', 'Finished')