Sweep over a 2D mesh. More...

Inheritance diagram for Sweep2D_mod:

Collaboration diagram for Sweep2D_mod:

Public Member Functions
function	Sweep2D_mod (in input, in mesh, in mat, in quadrature, in boundary, in equation)
	Class constructor.
function	sweep (in obj, in s, in g)
	Sweep the mesh for all angles.
Public Attributes
Property	d_input
Property	d_mesh
Property	d_mat
Property	d_quadrature
Property	d_equation
Property	d_alpha
Property	d_psi_horizontal
Property	d_boundary
Property	d_nx
Property	d_ny
Private Member Functions
function	x_bounds (in obj, in o)
function	y_bounds (in obj, in o)

Detailed Description

Sweep over a 2D mesh.

The within-group transport equation is

$\mathbf{T}\psi = Q \, ,$

where $\mathbf{T}$ is the streaming and collision operator and $ Q $ is a discrete representation of all source contributions.

To invert the operator $\mathbf{T}$ , we "sweep" over the mesh for all angles, which gives us updated angular fluxes in each cell. However, we don't store the angular flux, but rather add its contribution to the scalar flux directly since storing the angular flux is too expensive.

Constructor & Destructor Documentation

function Sweep2D_mod	(	in	input,
		in	mesh,
		in	mat,
		in	quadrature,
		in	boundary,
		in	equation
	)

Class constructor.

Parameters:

input	Number of fine mesh per x coarse mesh.
mesh	Number of fine mesh per y coarse mesh.
mat	Coarse mesh edges along x axis.
quadrature	Coarse mesh edges along y axis.
boundary	Coarse mesh material map.

Returns:: Instance of the Mesh class.

Member Function Documentation

function sweep	(	in	obj,
		in	s,
		in	g
	)

Sweep the mesh for all angles.

This performs the action of $\mathbf{T}^{-1}$ on a given discrete right hand side vector. Currently, a single vector applicable to all angles is provided, since we work only with isotropic sources and scattering.