The primary impediments impeding the implementation of high-order methods in simulating viscous flow over complex configurations are robustness and convergence.These challenges impose significant constraints on comput...The primary impediments impeding the implementation of high-order methods in simulating viscous flow over complex configurations are robustness and convergence.These challenges impose significant constraints on computational efficiency,particularly in the domain of engineering applications.To address these concerns,this paper proposes a robust implicit high-order discontinuous Galerkin(DG)method for solving compressible Navier-Stokes(NS)equations on arbitrary grids.The method achieves a favorable equilibrium between computational stability and efficiency.To solve the linear system,an exact Jacobian matrix solving strategy is employed for preconditioning and matrix-vector generation in the generalized minimal residual(GMRES)method.This approach mitigates numerical errors in Jacobian solution during implicit calculations and facilitates the implementation of an adaptive Courant-Friedrichs-Lewy(CFL)number increasing strategy,with the aim of improving convergence and robustness.To further enhance the applicability of the proposed method for intricate grid distortions,all simulations are performed in the reference domain.This practice significantly improves the reversibility of the mass matrix in implicit calculations.A comprehensive analysis of various parameters influencing computational stability and efficiency is conducted,including CFL number,Krylov subspace size,and GMRES convergence criteria.The computed results from a series of numerical test cases demonstrate the promising results achieved by combining the DG method,GMRES solver,exact Jacobian matrix,adaptive CFL number,and reference domain calculations in terms of robustness,convergence,and accuracy.These analysis results can serve as a reference for implicit computation in high-order calculations.展开更多
This paper proposes a vector-borne plant disease model with discontinuous treatment strategies. Constructing Lyapunov function and applying non-smooth theory to analyze discontinuous differential equations, the basic ...This paper proposes a vector-borne plant disease model with discontinuous treatment strategies. Constructing Lyapunov function and applying non-smooth theory to analyze discontinuous differential equations, the basic reproductive number R0 is proved, which determines whether the plant disease will be extinct or not. If R0 R0 > 1 , there exists a unique endemic equilibrium which is globally stable. The numerical simulations are provided to verify our theoretical results, which indicate that after infective individuals reach some level, strengthening treatment measures is proved to be beneficial in controlling disease transmission.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.12102247)the Technology Development Program(Grant No.JCKY2022110C119).
文摘The primary impediments impeding the implementation of high-order methods in simulating viscous flow over complex configurations are robustness and convergence.These challenges impose significant constraints on computational efficiency,particularly in the domain of engineering applications.To address these concerns,this paper proposes a robust implicit high-order discontinuous Galerkin(DG)method for solving compressible Navier-Stokes(NS)equations on arbitrary grids.The method achieves a favorable equilibrium between computational stability and efficiency.To solve the linear system,an exact Jacobian matrix solving strategy is employed for preconditioning and matrix-vector generation in the generalized minimal residual(GMRES)method.This approach mitigates numerical errors in Jacobian solution during implicit calculations and facilitates the implementation of an adaptive Courant-Friedrichs-Lewy(CFL)number increasing strategy,with the aim of improving convergence and robustness.To further enhance the applicability of the proposed method for intricate grid distortions,all simulations are performed in the reference domain.This practice significantly improves the reversibility of the mass matrix in implicit calculations.A comprehensive analysis of various parameters influencing computational stability and efficiency is conducted,including CFL number,Krylov subspace size,and GMRES convergence criteria.The computed results from a series of numerical test cases demonstrate the promising results achieved by combining the DG method,GMRES solver,exact Jacobian matrix,adaptive CFL number,and reference domain calculations in terms of robustness,convergence,and accuracy.These analysis results can serve as a reference for implicit computation in high-order calculations.
文摘This paper proposes a vector-borne plant disease model with discontinuous treatment strategies. Constructing Lyapunov function and applying non-smooth theory to analyze discontinuous differential equations, the basic reproductive number R0 is proved, which determines whether the plant disease will be extinct or not. If R0 R0 > 1 , there exists a unique endemic equilibrium which is globally stable. The numerical simulations are provided to verify our theoretical results, which indicate that after infective individuals reach some level, strengthening treatment measures is proved to be beneficial in controlling disease transmission.
