Moltres init: cross-code or experimental comparison about publications HOT 20 CLOSED

MPACT - MSR features?

from publications.

OK, here's how I'd like to address this. I'd run a 2D axisymmetric case in both coupled and decoupled mode, both steady state and transient, and a coarse and fine mesh, for 8 cases total. This would give info about the effects of coupling. In particular, the transient cases should illustrate the importance of fully coupling, especially with higher reactivity insertions. I'm just not sure what transient case would be best... probably the control rod movement.

from publications.

@lindsayad it seems like this issue sort of coincides with issue #43. Maybe I could set up the transient cases, and then you modify them for both a decoupled and a fully coupled solution? If that's not too much burden for you, that is.

from publications.

I would contend that the title of this issue doesn't really match the associated reviewer quote. Instead as you say I think this quote belongs with issue #43. More closely associated with the title of this issue is in the reviewer's summary:

The paper would represent an interesting contribution to the scientific community but it requires revision. In particular the authors should deepen their bibliographical research, and include a comparison either with experimental results or with other state-of-the-art tools.

There is certainly overlap. I guess in a perfect world with unlimited time we would try to produce the "best" results with Moltres any way we can. When there are trustworthy experimental results, there is no ambiguity about what "best" means; the calculations should reproduce the experimental results. However, if there aren't trustworthy experimental results then who is to say what is best? How can we really say what model is high fidelity and what isn't? We use a more a priori model, a finer mesh, higher degree polynomials, full physics coupling? Yes things like that are more probably more likely to produce a more physical result but we just don't know without real results.

I also don't really know how the reviewer got "slow performance" in his/her head. Is he assuming that because we use a lot of cores?

Normally problems of few million degrees of freedom can be handled with few tens of cores.

My response to that would be: if you have the resources, why not use them? Why solve something in an hour if you have the parallelizability to solve it in minutes?!

So returning to practical matters @gridley, I think it would be great if you (or anybody) could try and do some comparison to experimental or other codes' results using our default fully coupled treatment. That addresses the title of this issue in my mind. In terms of "high fidelity", there is no reason that full coupling shouldn't be just as or more accurate than a segregated approach. Then to address issue #43 we can compare segregated to monolithic calculations (which we can then also of course use in comparison to the experimental/cross-code results). So yes definitely some overlap, but I think we can keep both open. Waddya think? (after my long diatribe)

from publications.

I should say that if you iterate between physics and use the same tolerances in a segregated as in a monolithic solve, then the accuracy of the two solution methods should be the same. It's just that in very tight physics couplings, the iterative segregated approach may not converge. Also the resources used will probably be different. Fully coupled may use a lot more memory for example, but may also be faster. We won't know that until we try. I think resource comparison falls under the umbrella of #43

from publications.

So far only found measurements for control rod worth in zero power state: http://moltensalt.org/references/static/downloads/pdf/ORNL-TM-1626.pdf
Keep digging

@lindsayad @gridley Why you don't want to cross-verify moltres with Cammi's Comsol model?

from publications.

I found some neutron flux profile, which slightly lower than reference used in moltres paper (fast flux ~11 instead of 14). Looks like it some simulation (not an experiment) but, probably, useful:

Source: http://www.thmfgrcs.com/ORNL-3872.pdf

from publications.

from publications.

Yea, I still try to find experimental data, 22 semiannual reports, 300-40 page each, take time :)
Don't believe in conspiracy theory but looks like someone intentionally didn't include interesting flux and temperature measurements in those reports...

from publications.

from publications.

So, @lindsayad @andrewryh, what do you al lsay that we reproduce Cammi results for the single channel? The mesh would be simple, and IIRC we already have all of the physics coded in that they use. I think the biggest difference would be computing cross sections from parsed functions rather than interpolation tables. Also, they used a k-\epsilon turbulence model, which @lindsayad got working, right?

There are a few things we could reproduce, and not all would be necessarily needed. It seems that they calculate feedback coefficients in an infinite reflecting channel in eigenvalue calculation mode, which would be relatively easy after getting the mesh made since we could use Cammi's group constants.

In terms of the fluids part, we could use Alex's Petrov-Galerkin scheme that allows P1 lagrange velocity variables (which I have no clue how to use, is this online?).

Well, this seems like a lot, but we have a lot of time and could set up a small timeline if needed. 😄

from publications.

@gridley (you didn't ask me, but ... ) This is a good idea and worthwile endeavor, but I have 2 concerns.

1. I don't think reproducing the Cammi single channel result completely satisfies the desire expressed by the reviewer (who desired comparison to high fidelity simulations or experiment). It would be a helpful example, certainly, but experimental validation or high-fidelity-code-to-code verification. While the single channel model represents a modern simulation with coupled physics, our feature claims in this paper are related to more complex geometries (2-D axisymmetric many channel and 3-D many channel), so they should be verified against complex geometries, IMO.

2. Alex tried to replicate Cammi and failed to get agreement early on in Fall 2016, likely based on miscommunication in their paper which was never clarified despite many attempts to clarify it with the authors (is that this same single channel example, @lindsayad ? I forgot.)

from publications.

from publications.

from publications.

Well, I have contact both Manuele and Dr. Collins over the past two weeks, and they seem hesitant to let go of any data at the moment. So, on this front, we may be stuck unless some experimental results are found. We could possibly use these cross sections I found from Cammi in a presentation here and see if our cross sections look anything like that. @lindsayad, where are the group constants you generated, out of curiosity?

from publications.

from publications.

Hm, right, I see how your results would be supercritical because, judging from property_file_dir/newt_fuel_NSF.txt, your NSF cross sections for fast and thermal groups are respectively 22% and 30% higher than Cammi's presentation. I'm going to plug this single channel into Serpent and see what happens. I'd imagine some discrepancy could arise due to collapsing XSs in the supercritical spectrum rather than critical, and Serpent definitely does leakage corrections before averaging over the spectrum. As for NEWT, no clue.

from publications.

Well, OK, even if your case was quite supercritical @lindsayad, Cammi et. al. did add artrificial neutron poison, and from what I've found they don't provide how much Sigma_a was perturbed by. It may have been large.

from publications.

from publications.

Not finished for this paper per se, but we left comments on this point.

Closed by #66

from publications.