test/modules/LinearBoltzmannSolvers/Transport_Keigen/KEigenvalueTransport1D_1G.lua

-- 1D 1G KEigenvalue::Solver test using power iteration
-- Test: Final k-eigenvalue: 0.9995433
num_procs = 4
 
-- NOTE: For command line inputs, specify as:
--       variable=[[argument]]
 
--############################################### Check num_procs
if (check_num_procs == nil and chi_number_of_processes ~= num_procs) then
    chiLog(LOG_0ERROR,"Incorrect amount of processors. " ..
                      "Expected "..tostring(num_procs)..
                      ". Pass check_num_procs=false to override if possible.")
    os.exit(false)
end
 
chiMPIBarrier()
 
-- ##################################################
-- ##### Parameters #####
-- ##################################################
 
-- Mesh variables
if (L == nil) then L = 100.0 end
if (n_cells == nil) then n_cells = 50 end
 
-- Transport angle information
if (n_angles == nil) then n_angles = 16 end
if (scat_order == nil) then scat_order = 0 end
 
-- k-eigenvalue iteration parameters
if (kes_max_iterations == nil) then kes_max_iterations = 5000 end
if (kes_tolerance == nil) then kes_tolerance = 1e-8 end
 
-- Source iteration parameters
if (si_max_iterations == nil) then si_max_iterations = 500 end
if (si_tolerance == nil) then si_tolerance = 1e-8 end
 
-- Delayed neutrons
if (use_precursors == nil) then use_precursors = true end
 
 
-- ##################################################
-- ##### Run problem #####
-- ##################################################
 
--############################################### Setup mesh
nodes = {}
dx = L/n_cells
for i=0,n_cells do
  nodes[i+1] = i*dx
end
 
meshgen1 = chi_mesh.OrthogonalMeshGenerator.Create({ node_sets = {nodes} })
chi_mesh.MeshGenerator.Execute(meshgen1)
 
--############################################### Set Material IDs
chiVolumeMesherSetMatIDToAll(0)
 
--############################################### Add materials
materials = {}
materials[1] = chiPhysicsAddMaterial("Fissile Material")
 
chiPhysicsMaterialAddProperty(materials[1], TRANSPORT_XSECTIONS)
 
xs_file = "simple_fissile.cxs"
chiPhysicsMaterialSetProperty(materials[1], TRANSPORT_XSECTIONS,
                              CHI_XSFILE, xs_file)
 
--############################################### Setup Physics
num_groups = 1
lbs_block =
{
  num_groups = num_groups,
  groupsets =
  {
    {
      groups_from_to = {0, num_groups-1},
      angular_quadrature_handle =
        chiCreateProductQuadrature(GAUSS_LEGENDRE,n_angles),
      inner_linear_method = "gmres",
      l_max_its = si_max_iterations,
      l_abs_tol = si_tolerance,
    }
  }
}
 
lbs_options =
{
  scattering_order = scat_order,
 
  use_precursors = use_precursors ,
 
  verbose_inner_iterations = false,
  verbose_outer_iterations = true,
}
 
phys = lbs.DiscreteOrdinatesSolver.Create(lbs_block)
lbs.SetOptions(phys, lbs_options)
 
k_solver0 = lbs.XXNonLinearKEigen.Create
({
  lbs_solver_handle = phys,
  nl_max_its = kes_max_iterations,
  nl_abs_tol = kes_tolerance
})
chiSolverInitialize(k_solver0)
chiSolverExecute(k_solver0)
 
--############################################### Get field functions
--############################################### Line plot
--############################################### Volume integrations
--############################################### Exports
--############################################### Plots