Block ID not functioning for PKAConstant & PKAFunction

[ChrocheMisawa]

Dear all,

    I recently found this application and found it particularly interesting, because my master's thesis was to couple MOOSE with MCNP and TRIM(in a more robust and simple way), and calculate the damage(dpa) caused in a multi-phase grain geometry. Therefore, I was exciting to see if I can reproduce my work done in the past with this new app.
    So I started working with simple test cases, first with just a PKAGun shooting ions to a homogenous material, which worked perfectly fine and agrees with the normal TRIM results. However, when I was trying to create an isotropic ion source(with PKAConstant or PKAFunction) in just a local region with specific block ID, the calculation gives results that indicate PKAs were still being distributed across the entire mesh.
    Is it due the application itself, or am I writing the input file wrong? I have attached my input file below, it would be great if I can get any help.

Thank you, Zeyu Chen

[Mesh] [./generator] type = GeneratedMeshGenerator dim = 3 nx = 50 ny = 1 nz = 1 xmin = -50
xmax = 50 ymin = -1 ymax = 1 zmin = -1 zmax = 1 elem_type = HEX8 [] [block1] type = SubdomainBoundingBoxGenerator input = generator bottom_left = '-50 -1 -1' top_right = '-1 1 1' block_id = 0 [] [block2] type = SubdomainBoundingBoxGenerator input = block1 bottom_left = '-1 -1 -1' top_right = '1 1 1' block_id = 1 [] [block3] type = SubdomainBoundingBoxGenerator input = block2 bottom_left = '1 -1 -1' top_right = '50 1 1' block_id = 2 []

[]

[Mesh]

type = MyTRIMMesh

dim = 3

nx = 50

ny = 1

nz = 1

xmin = -50

xmax = 50

ymin = -1

ymax = 1

zmin = -1

zmax = 1

block = 1

[]

[Problem] kernel_coverage_check = false []

[Variables] [./dummys] [../] []

[AuxVariables] [./Be] order = CONSTANT family = MONOMIAL initial_condition = 1 [../]

[./vac] order = CONSTANT family = MONOMIAL [../] []

[AuxKernels] [./vac] variable = vac type = MyTRIMElementResultAux runner = runner ivar = 0 defect = VAC [../] []

[BCs] [./Periodic] [./all] auto_direction = 'x y z' [../] [../]

[]

[UserObjects] [./constant] type = PKAConstant E = 2.73e6 Z = 1 m = 3 pka_rate = 1 block = 1 use_displaced_mesh = true

force_preic =true

# seed = 1

[../]

[./constant]

type = PKAGun

E = 2.73e6

Z = 1

m = 3

num_pkas = 10000

point = '-50 0 0'

direction = '-1 0 0'

[../]

[./rasterizer] type = MyTRIMRasterizer var = 'Be' M = '9' Z = '4'

Edisp = '16.3'

length_unit = MICROMETER
periodic_var = dummys
site_volume = 0.00797
pka_generator = constant

[../] [./runner] type = MyTRIMElementRun rasterizer = rasterizer [../] []

[Postprocessors] [./total_vacancies] type = ElementIntegralVariablePostprocessor variable = vac [../] []

[Executioner] type = Transient num_steps = 2 nl_abs_tol = 1e-10 []

[Outputs] exodus = true []

[dschwen]

Block restriction is not implemented for the BCMC. It should not be hard to add if it is required.

[ChrocheMisawa]

That explains, thank you for your reply!

I think it would be nice to be implemented it so that users can study the impact just from some ion-emitting material. Such as the plot in the link below: https://drive.google.com/file/d/1m6nVBEBNSAZ0oUktCaaUrwlBDMSTKOqF/view?usp=sharing

Also, I'm not sure how MyTRIM handles the BCMC, but it might be nice if the application can also gather information on particles' end location (instead of just PKAs or Energy Deposition), such that the user can study some features such as Helium accumulation on certain geometries.

Thank you and best wishes, Zeyu Chen