We derive an adjoint method for the Direct Simulation Monte Carlo (DSMC) method for the spatially homogeneous Boltzmann equation with a general collision law. This generalizes our previous results in [Caflisch, R., Silantyev, D. and Yang, Y., 2021. Journal of Computational Physics, 439, p.110404], which was restricted to the case of Maxwell molecules, for which the collision rate is constant. The main difficulty in generalizing the previous results is that a rejection sampling step is required in the DSMC algorithm in order to handle the variable collision rate. We find a new term corresponding to the so-called score function in the adjoint equation and a new adjoint Jacobian matrix capturing the dependence of the collision parameter on the velocities. The new formula works for a much more general class of collision models.
翻译:我们从直接模拟蒙特卡洛(DSMC)法中得出一种联合方法,用于空间均匀的布尔茨曼方程式,并附有一般碰撞法。这概括了我们以前在[Caflisch、R、Silantyev、D和Yang,Y.,2021]中的结果。《计算物理杂志》,439,p.1104404]中的结果,该结果仅限于马克斯韦尔分子的情况,其碰撞率保持不变。推广先前结果的主要困难是,DSMC算法中需要一个拒绝抽样步骤,以便处理可变碰撞率。我们发现一个新的术语与所谓的对等函数相对应,并找到一个新的联合 Jacobian矩阵,它反映了碰撞参数对速度的依赖性。新的公式为更一般性的碰撞模型工作提供了更普遍的分类。