Time integration methods for solving initial value problems are an important component of many scientific and engineering simulations. Implicit time integrators are desirable for their stability properties, significantly relaxing restrictions on timestep size. However, implicit methods require solutions to one or more systems of nonlinear equations at each timestep, which for large simulations can be prohibitively expensive. This paper introduces a new family of linearly implicit multistep methods (LIMM), which only requires the solution of one linear system per timestep. Order conditions and stability theory for these methods are presented, as well as design and implementation considerations. Practical methods of order up to five are developed that have similar error coefficients, but improved stability regions, when compared to the widely used BDF methods. Numerical testing of a self-starting variable stepsize and variable order implementation of the new LIMM methods shows measurable performance improvement over a similar BDF implementation.
翻译:解决初始价值问题的时间整合方法是许多科学和工程模拟的一个重要部分,隐含的时间集成器对于其稳定性是可取的,大大放宽了对时间步骤大小的限制,然而,隐含的方法要求每个时间步骤解决一个或多个非线性方程式系统,而对于大型模拟来说,这种系统的费用太高,令人望而却步;本文件引入了一套新的线性隐含的多步法(LIMM),它只要求每个时间步骤解决一条线性系统;提出了这些方法的秩序条件和稳定性理论,以及设计和执行方面的考虑;制定了多达五种实际的秩序方法,这些方法有类似的误差系数,但与广泛使用的BDF方法相比,稳定性区域有所改善;对自启动的可变步骤和新的LIMM方法的可变顺序实施进行数值测试,表明在类似的BDF实施过程中可以衡量的绩效改进。