I think your boundary conditions has some problems.
Transmission boundary between assemblies
Vacuum boundary outside of the geometry &
Periodic boundary condition should be used in the two planes which is covering your 1/6 th portion of the geometry.
Also, if the geometry is infinite in z-axis you can use to planes in z-axis to bind it with reflective boundary conditions. About the vacuum boundary condition, it depends on your geometry. If it is infinite without any leakage then you should use reflective boundary condition instead of vacuum boundary as mentioned in point 2 above.
I found the reasion is that I used complicated surface for the bottom of the geometry, geometry.xml input is given at the end.
the openmc runs well when I change the bottom surface to a y-plane, but this boundary is not what i need for my benchmark problem. by the way, my benchmark problem is simple 5 rings hexagonal problem with pin pitch sqrt(3.0), the boundary condition for all outer surface is reflective.
then my question is how to model the benchmark problem?
thank you very much.
I changed the geometry to this, then OpenMC can run well.
Follow the example, create two universes one for fuel & another with moderator (outer universe). start from ring 5 instead of ring 7 & this will give you an assembly with 5 rings. Place the universes according to your model description. Give lattice pitch sqrt (3) instead of 1.275. Put the value for edge length of the assembly according to your assembly pitch.
Sorry to say that my problem is a multi-pin problem (61 pins in total) WITHOUT assembly edge, so that my boundary is the outest pins’ out side. this problem is defined not similar to real hexagonal assembly in order to test hexagonal deterministic method (e.g. Hexagonal Sn method or SP3 method).
the forum system just allow me to upload xml file, so the XS python file and geometry.xml and materials.xml file are attached. the Reference keff of this problem (without Gd) is about 1.2183. cuurently I can’t accurately incorperate reflective boundary condition as I want, so the result from OpenMC is 1.1797, this value is equal to the result for all the pins being fuel only; I guess this is because the boundary conditon setted in my geometry.xml bypass the outest ring of water pins.
could developer kindly reply if OpenMC can model this simple problem?
It is just a 7-pin hexagonal problem with reflective outer boundary conditon, the difficulty for me is how to model it WITHOUT “assembly edge”.
Thanks you very much.
I just want to use it as a benchmark result for my deterministic code developed.