Approximate reflective MOC boundary condition class. More...

Inheritance diagram for AppxReflectiveTrack:

Collaboration diagram for AppxReflectiveTrack:

Public Member Functions
function	AppxReflectiveTrack (in boundary, in mesh, in quadrature, in side)
	Class constructor.
function	update (in this)
	Update the boundary flux.
function	initialize (in this)
	Initialize the boundary condition.
Public Attributes
Property	d_index
Property	d_number_tracks
Property	d_group_a
	Index of all tracks as a function of local angle index (increasing x)
Property	d_group_x
	Index of all tracks as a function of spatial point (increasing phi)
Property	d_n_space
Property	d_n_angle
Property	d_order_angle
	Angular expansion order. Maximum is number angles - 1.
Property	d_order_space
	Spatial expansion order. Maximum is number of points - 1.
Property	d_basis_space
	Discrete Legendre polynomial spatial basis.
Property	d_basis_angle
	Discrete Legendre polynomial angular basis.
Protected Member Functions
function	build_index (in this)
	Build simple table of what I grab and what I set.
function	build_basis (in this)
	Build my basis.
function	build_group (in this)
	Build space and angle index groups.
function	expand (in this, in psi_exact)
	Expand the reflected flux and approximate.
Protected Attributes
Property	d_input
Property	d_boundary
	Boundary flux thisect.
Property	d_mesh
	Mesh.
Property	d_quadrature
	Quadrature.
Property	d_side

Detailed Description

Approximate reflective MOC boundary condition class.

The reflective condition is approximated via an (n,m)th order expansion in space(n) and azimuthal angle(m).

Constructor & Destructor Documentation

function AppxReflectiveTrack	(	in	boundary,
		in	mesh,
		in	quadrature,
		in	side
	)

Class constructor.

More detailed description of what the constructor does.

Parameters:

boundary	Boundary flux class.
side	Surface identifier.

Returns:: Instance of the AppxReflectiveTrack class.

Member Function Documentation

function build_basis ( in this ) [protected]

Build my basis.

For m azimuthal angles per quadrant, we need 1+log_3 m spatial basis sets. For n=9 and n=27, we need 3 and 4 sets. The sets have 1, 3, 9, and 27 spatial points. The lone point isn't really a basis, though...

We need 3 and 4 angle sets, too. For n=27, there are 7 angles per quadrant, or 14 over the relevant half-space. Then 10, 6, 4, and 2.

function build_group ( in this ) [protected]

Build space and angle index groups.

This builds two sets of indices. The first contains a vector of all tracks at a particular spatial point, ordered by increasing angle. Hence, an angular basis can be used directly on the values at that point.

The second is a list of all tracks for a given angle, ordered by increasing (local) spatial location. Hence, the spatial basis can be used directly for this set of tracks.

function build_index ( in this ) [protected]

Build simple table of what I grab and what I set.

function expand	(	in	this,
		in	psi_exact
	)		`[protected]`

Expand the reflected flux and approximate.

This takes the exact reflected flux and replaces it with an (n,m)th order approximation in space and angle.

Todo:: Study the effect of approximating space or angle first. It doesn't matter in a full expansion, but it does for low orders as they are not independent variables.

function initialize ( in this ) [inherited]

Initialize the boundary condition.

For fixed boundaries, this is the only relevant call. Other need not implement this.

Reimplemented in IsotropicBC.

function update ( in this )

Update the boundary flux.