Importance of curvature evaluation scale for predictive simulations of dynamic gas-liquid interfaces

JOURNAL OF COMPUTATIONAL PHYSICS

卷:365

Article

Owkes, Mark¹  Cauble, Eric¹  Senecal, Jacob¹  Currie, Robert A.¹ 

[1] Montana State Univ, Mech & Ind Engn, Bozeman, MT 59717 USA

关键词: Adjustable Curvature Evaluation Scale (ACES);    Surface tension;    Volume of fluid;    Marker particle;    Height function;    Discretization scale;   

DOI  :  10.1016/j.jcp.2018.03.018

来源: Elsevier

【 摘 要 】

The effect of the scale used to compute the interfacial curvature on the prediction of dynamic gas-liquid interfaces is investigated. A new interface curvature calculation methodology referred to herein as the Adjustable Curvature Evaluation Scale (ACES) is proposed. ACES leverages a weighted least squares regression to fit a polynomial through points computed on the volume-of-fluid representation of the gas-liquid interface. The interface curvature is evaluated from this polynomial. Varying the least squares weight with distance from the location where the curvature is being computed, adjusts the scale the curvature is evaluated on. ACES is verified using canonical static test cases and compared against second- and fourth-order height function methods. Simulations of dynamic interfaces, including a standing wave and oscillating droplet, are performed to assess the impact of the curvature evaluation scale for predicting interface motions. ACES and the height function methods are combined with two different unsplit geometric volume-of-fluid (VoF) schemes that define the interface on meshes with different levels of refinement. We find that the results depend significantly on curvature evaluation scale. Particularly, the ACES scheme with a properly chosen weight function is accurate, but fails when the scale is too small or large. Surprisingly, the second-order height function method is more accurate than the fourth-order variant for the dynamic tests even though the fourth-order method performs better for static interfaces. Comparing the curvature evaluation scale of the second- and fourth-order height function methods, we find the second-order method is closer to the optimum scale identified with ACES. This result suggests that the curvature scale is driving the accuracy of the dynamics. This work highlights the importance of studying numerical methods with realistic (dynamic) test cases and that the interactions of the various discretizations is as important as the accuracy of one part of the discretization. (C) 2018 Elsevier Inc. All rights reserved.

【 授权许可】

Free