A network partition method for solving large-scale complex nonlinear processes

A numerical framework based on network partition and operator splitting is developed to solve nonlinear differential equations of large-scale dynamic processes encountered in physics, chemistry and biology. Under the assumption that those dynamic processes can be characterized by sparse networks, we minimize the number of splitting for constructing subproblems by network partition. Then the numerical simulation of the original system is simplified by solving a small number of subproblems, with each containing uncorrelated elementary processes. In this way, numerical difficulties of conventional methods encountered in large-scale systems such as numerical instability, negative solutions, and convergence issue are avoided. In addition, parallel simulations for each subproblem can be achieved, which is beneficial for large-scale systems. Examples with complex underlying nonlinear processes, including chemical reactions and reaction-diffusion on networks, demonstrate that this method generates convergent solution in a efficient and robust way.