Need help. Failed to reproduce the tutorial "csg_to_cad" of Cardinal

User Support
1

I am following a multi-physics coupling tutorial “csg_to_cad” of Cardinal. I don’t make any modifications to the code provided in the tutorial, and the mesh file is also generated following the instruction from tutorial. An error occurred when I tried to run this coupled code, and part of the detailed information is as following:

Time Step 1, time = 1, dt = 1

Skinning geometry into 4 temperature bins, 1 density bins, and 1 block bins... 

Re-generating OpenMC model with 53 cells... 

Mapping between 97620 MOOSE elements and 56 OpenMC cells (on 2 coordinate levels):
------------------------------------------------------------------------------
|              | # T Elems | # rho Elems | # T+rho Elems | # Uncoupled Elems |
------------------------------------------------------------------------------
| MOOSE mesh   |     97620 |           0 |             0 |                 0 |
| OpenMC cells |     97620 |           0 |             0 |                 0 |
------------------------------------------------------------------------------


 ===================>  OPENMC SUBDOMAIN MATERIAL MAPPING  <====================

      Subdomain:  Subdomain name; if unnamed, we show the ID
       Material:  OpenMC material name(s) in this subdomain; if unnamed, we
                  show the ID. If N duplicate material names, we show the
                  number in ( ).

------------------------------------
| Subdomain |       Material       |
------------------------------------
| 1         |  U-7.65Mo Fuel, VOID |
------------------------------------

Sending temperature to OpenMC cells... 
 Sent cell-averaged min/max (K): 800, 800
 Running OpenMC with 10000 particles per batch...
Executing filter editors...
Executing tally editors...
 Maximum neutron transport energy: 20000000 eV for U235
 Initializing source particles...

 ====================>     K EIGENVALUE SIMULATION     <====================

  Bat./Gen.      k            Average k
  =========   ========   ====================
 WARNING: No intersection found with DAGMC cell 627, filled with material 1
        1/1    0.71049
        2/1    0.95294
        3/1    0.96429
        4/1    0.97124
        5/1    0.95848
        6/1    0.95099
        7/1    0.95922
        8/1    0.95660
        9/1    0.96724
       10/1    0.97117
 WARNING: No intersection found with DAGMC cell 627, filled with material 1
       11/1    0.96293
       12/1    0.95076
       13/1    0.95921
       14/1    0.96144
       15/1    0.97190
       16/1    0.96344
       17/1    0.96541
       18/1    0.96644
 WARNING: No intersection found with DAGMC cell 627, filled with material 1
 WARNING: No intersection found with DAGMC cell 627, filled with material 1
       19/1    0.95665
       20/1    0.96901
       21/1    0.96026
       22/1    0.96099
       23/1    0.96805
 WARNING: No intersection found with DAGMC cell 627, filled with material 1
       24/1    0.96340
       25/1    0.95383
       26/1    0.92693
 WARNING: No intersection found with DAGMC cell 627, filled with material 1
       27/1    0.95250
       28/1    0.96456
       29/1    0.93986
 WARNING: No intersection found with DAGMC cell 627, filled with material 1
       30/1    0.95844
       31/1    0.97379
       32/1    0.95184
       33/1    0.97142
       34/1    0.96659
       35/1    0.95084
       36/1    0.96859
       37/1    0.96617
       38/1    0.94527
       39/1    0.96657
       40/1    0.95819
       41/1    0.95503
       42/1    0.96281
       43/1    0.94553
       44/1    0.95178
       45/1    0.95947
       46/1    0.95743
       47/1    0.94619
 WARNING: No intersection found with DAGMC cell 627, filled with material 1
       48/1    0.97706
       49/1    0.94506
       50/1    0.96280
       51/1    0.94494
 WARNING: No intersection found with DAGMC cell 627, filled with material 1
       52/1    0.98690    0.96592 +/- 0.02098
       53/1    0.96575    0.96586 +/- 0.01211
       54/1    0.98315    0.97019 +/- 0.00959
       55/1    0.95110    0.96637 +/- 0.00835
 WARNING: No intersection found with DAGMC cell 627, filled with material 1
 ERROR: Maximum number of lost particles has been reached.
Abort(-1) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, -1) - process 0
段错误 (核心已转储)

What confused me is that when I tried a less coarser mesh(about 70,000 elements in total) the openmc run normally, but another error occurred:

*** ERROR ***
/home/mechaos/Jnotebook/sim_std/openmc.i:3.1:
The following occurred in the Problem 'MOOSE Problem' of type OpenMCCellAverageProblem.

flux tallies do not match the global flux tally:
 Global value: 8564.38
 Tally sum:    660.323
 Difference:   7904.06

This means that the tallies created by Cardinal are missing some hits over the domain.
You can turn off this check by setting 'check_tally_sum' to false.

Additionally, when I tried run the code with commandcardinal-opt -i openmc.i --n-threads=4, openmc immediately threw the particles lost error, it seems that the number of threads could affect the execution. I don’t know if anyone else has also tried this tutorial and encountered the same problem, and I honestly don’t know exactly where the problem lies.
I sincerely hope someone could take a look at this post and help me find out the reason. Thanks in advance.

1 Like