C. Scalone,L. Einkemmer,J. Kusch,R. J. McClarren
Abstract:Computing the dominant eigenvalue is important in nuclear systems as it determines the stability of the system (i.e. whether the system is sub or supercritical). Recently, the work of Kusch, Whewell, McClarren and Frank \cite{KWMF} showed that performing a low-rank approximation can be very effective in reducing the high memory requirement and computational cost of such problems. In this work, we propose a rank adaptive approach that changes the rank during the inverse power iteration. This allows us to progressively increase the rank (i.e. changing the fidelity of the model) as we get closer to convergence, thereby further reducing computational cost. We then exploit this multi-fidelity approach to optimize a simplified nuclear reactor. In this case the system is parameterized and the values of the parameters that give criticality are sought.
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
This paper mainly focuses on **the development of efficient numerical algorithms in nuclear reactor criticality calculations**. Specifically, the authors propose a **multi - fidelity adaptive dynamic low - rank optimization algorithm** to solve the fission criticality problem. The following are the main objectives and problems of this paper:
1. **Calculate the dominant eigenvalue (effective multiplication factor)**:
- In a nuclear system, calculating the dominant eigenvalue (i.e., the effective multiplication factor \( k_{\text{eff}} \)) is very important because it determines the stability of the system (i.e., whether the system is sub - critical, critical, or super - critical). Mathematically,
\[
k_{\text{eff}} = 1
\]
indicates that the system is in a critical state; \( k_{\text{eff}}>1 \) indicates a super - critical state; \( k_{\text{eff}} < 1 \) indicates a sub - critical state.
2. **Reduce computational cost and memory requirements**:
- Traditional numerical methods require a large amount of memory and computational resources when calculating the dominant eigenvalue, especially in the case of high - dimensional or complex geometric structures. The authors propose a **low - rank approximation** method to reduce these requirements. The core idea of low - rank approximation is to simplify the calculation by reducing the rank of the matrix, thereby significantly reducing memory usage and computation time.
3. **Adaptively adjust the rank to improve accuracy**:
- To further improve computational efficiency, the authors introduce an **adaptive rank adjustment technique**. During the iteration process, as the convergence is approached, the rank (i.e., the fidelity of the model) is gradually increased, thereby further reducing the computational cost while ensuring accuracy.
4. **Optimize parameters to reach the critical state**:
- The authors also propose an optimization algorithm based on low - rank approximation to optimize the design parameters of the nuclear reactor (such as material thickness, etc.), so that the system reaches the critical state. This involves solving the following optimization problem:
\[
\alpha_1^*, \dots, \alpha_n^*=\arg\min_{\alpha_1, \dots, \alpha_n\geq0} f(\alpha_1, \dots, \alpha_n)
\]
where \( f(\alpha_1, \dots, \alpha_n) = |k_{\text{eff}}(\alpha_1, \dots, \alpha_n)-k_{\text{eff}}^*| \), and \( k_{\text{eff}}^* = 1 \) represents the critical state.
### Summary
This paper aims to develop an efficient numerical algorithm. Through the multi - fidelity adaptive dynamic low - rank method, it solves the problems of high computational cost and memory requirements in nuclear reactor criticality calculations and optimizes the design parameters to achieve the critical state. This method not only improves computational efficiency but also ensures the accuracy of the results.