Class CBCSweepChunk

Defined in File cbc_sweep_chunk.h

Inheritance Relationships

Base Type

public opensn::SweepChunk (Class SweepChunk)

Class Documentation

class CBCSweepChunk : public opensn::SweepChunk 

Implements the core sweep operation for a single cell within the cell-by-cell (CBC) sweep algorithm.

This class is responsible for performing the discrete ordinates transport calculation on a given cell for all angles and groups managed by its current AngleSet It interacts with a CBC_FLUDS object to obtain upwind angular flux data (from local neighbors, MPI remote buffers, or boundaries) and to store outgoing angular flux data (to local neighbors or MPI send buffers)

Public Functions

CBCSweepChunk(DiscreteOrdinatesProblem &problem, LBSGroupset &groupset)

virtual void SetAngleSet(AngleSet &angle_set) override: Set the current AngleSet.

virtual void SetCell(Cell const *cell_ptr, AngleSet &angle_set) override: Set the current cell to be swept.

virtual void Sweep(AngleSet &angle_set) override

Performs the discrete ordinates sweep calculation for the currently set cell, for all angles and groups within the provided AngleSet.

It:

Assembles the local transport equation system for each angle and group
Retrieves upwind angular fluxes from local neighbors, remote locations (via MPI data managed by CBC_FLUDS), or boundaries
Solves the local system for the outgoing angular fluxes at the cell nodes
Updates the global scalar flux moments
If save_angular_flux is true, stores the computed angular fluxes into the global angular flux vector
Propagates outgoing angular fluxes to local downwind neighbors or stages them for MPI transmission to remote downwind neighbors

Protected Attributes

CBC_FLUDS *fluds_

size_t gs_size_

unsigned int gs_gi_

size_t num_angles_in_as_

unsigned int group_stride_

size_t group_angle_stride_

bool surface_source_active_

const Cell *cell_

std::uint32_t cell_local_id_

const CellMapping *cell_mapping_

CellLBSView *cell_transport_view_

size_t cell_num_faces_

size_t cell_num_nodes_

DenseMatrix<Vector3> G_

DenseMatrix<double> M_

std::vector<DenseMatrix<double>> M_surf_

std::vector<Vector<double>> IntS_shapeI_