Computation_Unit

Structure used to consolidate the cross-sections, geometry, solvers and sources, execute transport calculations and extract its results.

Mandatory field(s)

• cross_sections::Cross_Sections : cross-section library.
• geometry::Geometry : geometry.
• solvers::Solvers : solvers.
• sources::Sources : fixed sources.

Optional field(s) - with default values

• N/A

set_cross_sections(this::Computation_Unit,cross_sections::Cross_Sections)

Assigns the cross-sections library to the computation unit.

Input Argument(s)

• this::Computation_Unit : computation unit.
• cross_sections::Cross_Sections : cross-sections library.

Output Argument(s)

N/A

Examples

julia> cs = Cross_Sections()
julia> ... # Define cross-sections properties and generate multigroup cross-sections.
julia> cu = Computation_Unit()
julia> cu.set_cross_sections(cs)

set_geometry(this::Computation_Unit,geometry::Geometry)

Assigns the geometry to the computation unit.

Input Argument(s)

• this::Computation_Unit : computation unit.
• geometry::Geometry : geometry.

Output Argument(s)

N/A

Examples

julia> geo = Geometry()
julia> ... # Define geometry and its properties
julia> cu = Computation_Unit()
julia> cu.set_geometry(geo)

setsolvers(this::ComputationUnit,solvers::Solvers)

Assigns the solvers to the computation unit.

Input Argument(s)

• this::Computation_Unit : computation unit.
• solvers::Solvers : collection of solvers per particle.

Output Argument(s)

N/A

Examples

julia> ms = Solvers()
julia> ... # Define all the discretization solvers and their properties
julia> cu = Computation_Unit()
julia> cu.set_solvers(ms)

set_sources(this::Computation_Unit,sources::Fixed_Sources)

Assigns the fixed sources to the computation unit.

Input Argument(s)

• this::Computation_Unit : computation unit.
• sources::Fixed_Sources : collection of fixed sources.

Output Argument(s)

N/A

Examples

julia> fs = Fixed_Sources()
julia> ... # Define all the fixed sources and their properties
julia> cu = Computation_Unit()
julia> cu.set_sources(fs)

run(this::Computation_Unit)

Execute transport calculations and obtain the flux solution.

Input Argument(s)

• this::Computation_Unit : computation unit.

Output Argument(s)

N/A

Examples

julia> cu = Computation_Unit()
julia> ... # Define the cross-sections, geometry, fixed sources and discretization solvers
julia> cu.run()

get_voxels_position(this::Computation_Unit,axis::String)

Get the mid-point voxels position along the specified axis.

Input Argument(s)

• this::Computation_Unit : computation unit.
• axis::String : axis, which can takes the following values:
  • boundary = "x" : along x-axis
  • boundary = "y" : along y-axis
  • boundary = "z" : along z-axis

Output Argument(s)

• x::Vector{Float64} : mid-point voxels position along the specified axis.

Examples

julia> cu = Computation_Unit()
julia> ... # Define computation unit and run it.
julia> x = cu.get_voxels_position("x")

get_energies(this::Computation_Unit,particle::Particle)

Get the mid-point energy in each group for the specified particle.

Input Argument(s)

• this::Computation_Unit : computation unit.
• particle::Particle : particle.

Output Argument(s)

• E::Vector{Float64} : mid-point energy in each group for the specified particle.

Examples

julia> electron = Electron() # Particle to be transported
julia> cu = Computation_Unit()
julia> ... # Define computation unit and run it.
julia> E = cu.get_voxels_position(electron)

get_flux(this::Computation_Unit,particle::Particle)

Get the array containing the flux in each voxels and in each energy group for the specified particle.

Input Argument(s)

• this::Computation_Unit : computation unit.
• particle::Particle : particle.

Output Argument(s)

• flux::Array{Float64} : flux array.

Examples

julia> electron = Electron() # Particle to be transported
julia> cu = Computation_Unit()
julia> ... # Define computation unit and run it.
julia> flux = cu.get_flux(electron)

get_energy_deposition(this::Computation_Unit,particle::Particle)

Get the array containing the energy deposition in each voxels by a given particle.

Input Argument(s)

• this::Computation_Unit : computation unit.
• particle::Particle : particle.

Output Argument(s)

• energy_deposition::Array{Float64} : energy deposition array.

Examples

julia> electron = Electron() # Particle to be transported
julia> cu = Computation_Unit()
julia> ... # Define computation unit and run it.
julia> electron_energy_deposition = cu.get_energy_deposition(electron)

get_energy_deposition(this::Computation_Unit)

Get the array containing the total energy deposition in each voxels.

Input Argument(s)

• this::Computation_Unit : computation unit.

Output Argument(s)

• energy_deposition::Array{Float64} : energy deposition array.

Examples

julia> cu = Computation_Unit()
julia> ... # Define computation unit and run it.
julia> energy_deposition = cu.get_energy_deposition()

get_charge_deposition(this::Computation_Unit,particle::Particle)

Get the array containing the charge deposition in each voxels by a given particle.

Input Argument(s)

• this::Computation_Unit : computation unit.
• particle::Particle : particle.

Output Argument(s)

• charge_deposition::Array{Float64} : charge deposition array.

Examples

julia> electron = Electron() # Particle to be transported
julia> cu = Computation_Unit()
julia> ... # Define computation unit and run it.
julia> electron_charge_deposition = cu.get_charge_deposition(electron)

get_charge_deposition(this::Computation_Unit)

Get the array containing the total charge deposition in each voxels.

Input Argument(s)

• this::Computation_Unit : computation unit.

Output Argument(s)

• charge_deposition::Array{Float64} : charge deposition array.

Examples

julia> electron = Electron() # Particle to be transported
julia> cu = Computation_Unit()
julia> ... # Define computation unit and run it.
julia> charge_deposition = cu.get_charge_deposition()