Excore Modeling with VERAShift

Omnibus User Manual

“…The model for a Virtual Environment for Reactor Applications (VERA) [9] XML input can be run through the Omnibus front end. The material compositions will be built from the model.…”

“…The Omnibus interface enables other codes such as ADVANTG [7] and SWORD [8] to create, execute, and post-process Denovo problems, but their use is outside the scope of this document. Other radiation transport codes such as VERAShift [9] and SCALE [2] use the Shift and Denovo codes through a lower-level interface via an internal Omnibus-based API. Although these codes are also outside the scope of this manual, the features and some interfaces presented here may inform the use of downstream codes.…”

“…from omnibus.data import plot matids = problem['matids'] zmats = matids.xs_by_index(z=slice(1,)) zmats.xs_by_index(y=slice(5, 15)).data[:] = 1 zmats.xs_by_index(x=slice(0,15), y=slice(10,15)).data[:] = 2 zmats.xs_by_index(y=slice(3,5)).data[:] = 2 plot(matids.xs(x=0.5), edgecolor='gray', linewidth=0.25); src = problem['volsrc'] spectra = src['spectra'] num_groups = 3 spectra.resize(2, num_groups) spectra.get_spectrum(0)[:] = [0.6, 0.3, 0.1] spectra.get_spectrum(1)[:] = [0.4, 0.2, 0.4] srcids = src['ids'] (continues on next page) A-85 (continued from previous page)srcids.xs_by_index(y=slice(3,5), z=slice(1,9)). data[:] = 1 plot(srcids.xs(z=1.0)); strength = src['strength'] q0 = np.linspace(0.5, 1.5, 15) q1 = np.linspace(0.5, 1.5, 20) strength.xs_by_index(x=slice(0,15), y=slice(10,15), z=slice(1,10)).data[:] = q0 strength.xs_by_index(y=slice(3,5), z=slice(1,10)).data[:] = q1 plot(strength.xs(z=1.0));A-86Now that the problem is constructed, it can be written to disk for later use in Omnibus: import h5py from omnibus.data import dump_root with h5py.File("simple-3d.h5", 'w') as f:dump_root(problem, f) INFO: Omitting empty point source A.2.6.2 Problem B This 2D problem is constructed similarly to the previous one.…”