交错网格上拉氏计算的总能守恒算法与界面做功计算

  • 摘要: 经典交错网格相容离散方法通常不直接给出单元边上的压强,因此难以精确计算内部子区域的界面做功.针对这一问题,本文在节点动量方程可写为隅角压力向量合力的前提下,引入外隅角力,给出一种从动量离散构造总能守恒内能格式的推导方法.该方法在隅角层面构造与动量方程离散相容的外隅角力,使节点周围外隅角力抵消,从而保持全域总能守恒.进一步,本文给出基于外隅角力的闭合子区域边界做功计算公式,证明该公式在给定离散格式意义下与子区域总能平衡完全一致.二维Sedov、Dukowicz、三维 Sedov 以及一个应用问题的计算结果表明,该方法能够精确统计全域和子区域能量守恒误差,并可用于提高复杂拉氏流体计算结果的可信度.

     

    Abstract: Classical compatible discretizations on staggered Lagrangian grids do not explicitly provide the pressure on cell faces, which makes the work done on internal interfaces difficult to evaluate by a face-pressure formula. This paper introduces an external corner-force method. Under the assumption that the nodal momentum update can be decomposed into corner pressure vectors, the method constructs a compatible internal-energy update without changing the original momentum scheme. The external corner forces cancel around each node and therefore preserve total energy. A discrete formula for the work done on closed sub-domain boundaries is then derived and proved to be exactly consistent with the discrete energy balance. Numerical examples demonstrate the effectiveness of the method for global and sub-domain energy-conservation diagnostics.

     

/

返回文章
返回