Abstract:Jacobian-free Newton Krylov (JFNK) is an attractive method to solve nonlinear equations in the nuclear engineering community, and has been successfully applied to steady-state neutron diffusion k-eigenvalue problems and multi-physics coupling problems. Preconditioning technique plays an important role in the JFNK algorithm, significantly affecting its computational efficiency. The key point is how to automatically construct a high-quality preconditioning matrix that can improve the convergence rate and perform the preconditioning matrix factorization efficiently and robustly. A reordering-based ILU(k) preconditioner is proposed to achieve the above objectives. In detail, the finite difference technique combined with the coloring algorithm is utilized to automatically construct a preconditioning matrix with low computational cost. Furthermore, the reordering algorithm is employed for the ILU(k) to reduce the additional non-zero elements and pursue robust computational performance. A 2D LRA neutron steady-state benchmark problem is used to evaluate the performance of the proposed preconditioning technique, and a steady-state neutron diffusion k-eigenvalue problem with thermal-hydraulic feedback is also utilized as a supplement. The results show that coloring algorithms can automatically and efficiently construct the preconditioning matrix. The computational efficiency of the FDP with coloring could be about 60 times higher than that of the preconditioner without the coloring algorithm. The reordering-based ILU(k) preconditioner shows excellent robustness, avoiding the effect of the fill-in level k choice in incomplete LU factorization. Moreover, its performances under different fill-in levels are comparable to the optimal computational cost with natural ordering.

New solution to ILU preconditioning parallelization and its applications in reservoir simulations

A Parallel Simulator for Massive Reservoir Models Utilizing Distributed-Memory Parallel Systems

Preconditioners for Parallel Reservoir Simulation.

Accelerating the iterative linear solver for reservoir simulation on multicore architectures

A Two-level GPU-Accelerated Incomplete LU Preconditioner for General Sparse Linear Systems

An Efficient and Robust ILU(k) Preconditioner for Steady-State Neutron Diffusion Problem Based on MOOSE

A staggered-grid multilevel incomplete LU for steady incompressible flows

A Parallelization Technique Based on Factor Combination and Graph Partitioning for General Incomplete LU Factorization

A parallel algorithm for unit commitment including AC power flow constraints

A Parallel Method With Hybrid Algorithms For Mixed Integer Nonlinear Programming

Parallel multilevel restricted Schwarz preconditioners for implicit simulation of subsurface flows with Peng-Robinson equation of state

23 European Symposium on Computer Aided Process Engineering: A Parallel Method with Hybrid Algorithms for Mixed Integer Nonlinear Programming

Regional-scale calculation of the LS factor using parallel processing

GPU Parallelization Nested Decomposition Method for Solving Large Linear Systems in Reservoir Numerical Simulation

Parallel Computing for Water Quality Equation Based on Jppf

Parallel Preconditioners for Large Scale Partial Difference Equation Systems

Parallelization of incomplete factorization preconditioning of block tridiagonal linear systems with 1-D domain decomposition

ILU Smoothers for Low Mach Navier-Stokes Pressure Solvers

Scalable Hierarchical Parallel Algorithm for the Solution of Super Large-Scale Sparse Linear Equations

Diagonalization-Based Parallel-in-Time Preconditioners for Instationary Fluid Flow Control Problems

Parallel Incomplete Factorization Preconditioning of Block Tridiagonal Linear Systems with 2-D Domain Decomposition