diff --git a/tests/postprocessors/gold/neutron_leakage_out.csv b/tests/postprocessors/gold/neutron_leakage_out.csv
new file mode 100644
index 0000000000..b4aa4f334f
--- /dev/null
+++ b/tests/postprocessors/gold/neutron_leakage_out.csv
@@ -0,0 +1,3 @@
+time,group1_leakage,group2_leakage,total_leakage
+0,0,0,0
+1,1,10,11
diff --git a/tests/postprocessors/neutron_leakage.i b/tests/postprocessors/neutron_leakage.i
new file mode 100644
index 0000000000..bc2feebd77
--- /dev/null
+++ b/tests/postprocessors/neutron_leakage.i
@@ -0,0 +1,103 @@
+[GlobalParams]
+  group_fluxes = 'group1 group2'
+  num_groups = 2
+  num_precursor_groups = 0
+  use_exp_form = false
+  temperature = 900
+  sss2_input = true
+[]
+
+[Mesh]
+  [gmg]
+    type = GeneratedMeshGenerator
+    dim = 2
+    nx = 2
+    ny = 2
+    xmax = 1
+    ymax = 1
+  []
+[]
+
+[Variables]
+  [group1]
+  []
+  [group2]
+  []
+[]
+
+[Kernels]
+  [group1_diffusion]
+    type = GroupDiffusion
+    variable = group1
+    group_number = 1
+  []
+  [group2_diffusion]
+    type = GroupDiffusion
+    variable = group2
+    group_number = 2
+  []
+[]
+
+[BCs]
+  [group1_left]
+    type = DirichletBC
+    variable = group1
+    boundary = left
+    value = 3
+  []
+  [group1_right]
+    type = DirichletBC
+    variable = group1
+    boundary = right
+    value = 2
+  []
+  [group2_left]
+    type = DirichletBC
+    variable = group2
+    boundary = left
+    value = 30
+  []
+  [group2_right]
+    type = DirichletBC
+    variable = group2
+    boundary = right
+    value = 20
+  []
+[]
+
+[Materials]
+  [mat]
+    type = MoltresJsonMaterial
+    base_file = mat.json
+    material_key = 'mat1'
+    interp_type = 'none'
+    group_constants = 'DIFFCOEF'
+  []
+[]
+
+[Executioner]
+  type = Steady
+[]
+
+[Postprocessors]
+  [group1_leakage]
+    type = NeutronLeakage
+    variable = group1
+    boundary = 'right'
+    group_number = 1
+  []
+  [group2_leakage]
+    type = NeutronLeakage
+    variable = group2
+    boundary = 'right'
+    group_number = 2
+  []
+  [total_leakage]
+    type = TotalNeutronLeakage
+    boundary = 'right'
+  []
+[]
+
+[Outputs]
+  csv = true
+[]
diff --git a/tests/postprocessors/tests b/tests/postprocessors/tests
index 019113fee1..3f7e389682 100644
--- a/tests/postprocessors/tests
+++ b/tests/postprocessors/tests
@@ -1,16 +1,20 @@
 [Tests]
-  [./side_weighted_integral]
-    type = 'Exodiff'
+  [side_weighted_integral]
+    type = Exodiff
     input = 'side_weighted_integral.i'
     exodiff = 'side_weighted_integral_out.e'
-
     requirement = 'The system shall compute the weighted integral sum of a variable over a side.' 
-  [../]
-  [./side_weighted_integral_RZ]
-    type = 'Exodiff'
+  []
+  [side_weighted_integral_RZ]
+    type = Exodiff
     input = 'side_weighted_integral_RZ.i'
     exodiff = 'side_weighted_integral_RZ_out.e'
-
     requirement = 'The system shall compute the weighted integral sum of a variable over a side in RZ coordinates.' 
-  [../]
+  []
+  [neutron_leakage]
+    type = CSVDiff
+    input = 'neutron_leakage.i'
+    csvdiff = 'neutron_leakage_out.csv'
+    requirement = 'The system shall compute group-wise and total neutron leakage'
+  []
 []