Hello,
For one of my use cases I’m trying to create a 17x17 assembly (reflective) with steady-state conditions EXCEPT for a changing fuel temperature distribution across the assembly. I can create the assembly just fine, and I have figured out how to set the temperature distributions for each fuel pin within each cell correctly, if the get_all_cells() function is to be believed. However, when running a full calculation, the tallied pin powers still show a completely flat distribution, like there’s absolutely no temperature variation. Intuitively this shouldn’t be the case, as my distribution is quite large (delta of ~500 K give or take); I think I should be seeing some variation. Is there a flag I’m not setting correctly? I can provide the input xml files if necessary.
@fshriver Apologies for the late response here. Temperature settings are controlled by:
settings = openmc.Settings()
settings.temperature = {...}
If you want cross section data to be interpolated in temperature, you should set
settings.temperature = {'method': 'interpolation'}
Give that a try and see if it does lead to a non-uniform power distribution, If you’re still having issues, do feel free to provide your Python code / input XML files as that may help to diagnose further.