An adaptive homotopy method for computing bifurcations of nonlinear parametric systems

Wenrui Hao,Chunyue Zheng
DOI: https://doi.org/10.48550/arXiv.2002.03460
2020-02-13
Abstract:In this paper, we present an adaptive step-size homotopy tracking method for computing bifurcation points of nonlinear systems. There are four components in this new method: 1) an adaptive tracking technique is developed near bifurcation points; 2) an inflation technique is backed up when the adaptive tracking fails; 3) Puiseux series interpolation is used to compute bifurcation points; and 4)the tangent cone structure of the bifurcation point is approximated numerically to compute solutions on different branches. Various numerical examples of nonlinear systems are given to illustrate the efficiency of this new approach. This new adaptive homotopy tracking method is also applied to a system of nonlinear PDEs and shows robustness and efficiency for large-scale nonlinear discretized systems.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: the numerical challenges faced when calculating bifurcation points in large - scale nonlinear parametric systems. Specifically, when parameter changes lead to changes in the solution structure, that is, the bifurcation phenomenon occurs. For large - scale nonlinear discretized systems (such as systems from nonlinear partial differential equations), the existing numerical methods are less efficient and computationally costly when dealing with these bifurcation points. Therefore, an efficient and adaptive homotopy tracking method is required to accurately calculate these bifurcation points. ### Main contributions of the paper The paper proposes an Adaptive Homotopy Tracking Method, which combines techniques in numerical algebraic geometry and scientific computing to improve the efficiency and accuracy of calculating bifurcation points in large - scale nonlinear systems. This method contains the following key components: 1. **Adaptive tracking technique**: Use adaptive step - size control when approaching bifurcation points to ensure the stability and accuracy of the tracking process. 2. **Inflation technique**: When the adaptive tracking fails, use the inflation technique to overcome ill - conditioned problems. 3. **Puiseux series interpolation**: Used for approximate calculation of bifurcation points. 4. **Numerical approximation of tangent cone structure**: Used to determine solutions on different branches. ### Mathematical formulas The mathematical formulas involved in the paper include, but are not limited to, the following: - Expression of the nonlinear parametric system: \[ F(u, p)=0 \] where \(u\in\mathbb{R}^n\) is the variable vector and \(p\in\mathbb{R}\) is the parameter. - Expression of the homotopy extended system: \[ \tilde{F}(u, p)= \begin{pmatrix} F(u, p)\\ gv^T(u - u_0)(1 - s)+s(p - p_0)-h \end{pmatrix} \] where \(g = \text{sign}\left(-v^T F_u(u_0, p_0)^{-1} F_p(u_0, p_0)\right)/\| F_u(\tilde{u}, \tilde{p})^{-1} F_p(\tilde{u}, \tilde{p})\|\), \(s = |\lambda_{\min}|/|\tilde{\lambda}_{\min}|\), \(\lambda_{\min}\) is the real part of the minimum eigenvalue of \(F_u\) at \((u_0, p_0)\), and \(v\) is the corresponding eigenvector. - Puiseux series expansion: \[ u(\lambda)=u_b+\sum_{j = 1}^{\infty}a_i\lambda^{j/c_1},\quad p(\lambda)=p_b+\sum_{j = 1}^{\infty}b_i\lambda^{j/c_2} \] Through these methods and techniques, the paper aims to provide an efficient numerical tool for studying and understanding bifurcation behaviors in complex large - scale nonlinear systems.