[论文解读] A fully nonlinear iterative solution method for self-similar potential flows with a free boundary

Алессандро Яфрати|arXiv (Cornell University)|Dec 30, 2012

Fluid Dynamics Simulations and Interactions参考文献 4被引用 3

一句话总结

本文提出了一种用于求解具有自由表面的自相似势流的全非线性迭代边界积分方法，采用伪时间推进法并引入修正速度势以简化自由表面边界条件。该方法在水下入水与水下撞击问题中准确捕捉了细射流和复杂流动特征，通过多个基准案例（包括楔形体入水、平板入水及倾斜平板水下撞击）的定量验证，展示了其高精度与计算效率。

ABSTRACT

An iterative solution method for fully nonlinear boundary value problems governing self-similar flows with a free boundary is presented. Specifically, the method is developed for application to water entry problems, which can be studied under the assumptions of an ideal and incompressible fluid with negligible gravity and surface tension effects. The approach is based on a pseudo time stepping procedure, which uses a boundary integral equation method for the solution of the Laplace problem governing the velocity potential at each iteration. In order to demonstrate the flexibility and the capabilities of the approach, several applications are presented: the classical wedge entry problem, which is also used for a validation of the approach, the block sliding along an inclined sea bed, the vertical water entry of a flat plate and the ditching of an inclined plate. The solution procedure is also applied to cases in which the body surface is either porous or perforated. Comparisons with numerical or experimental data available in literature are presented for the purpose of validation.

研究动机与目标

• 开发一种鲁棒的全非线性求解方法，用于自相似自由表面势流且边界未知的问题。
• 解决水下入水问题中细射流与强梯度区域的挑战，且不依赖渐近近似方法。
• 将边界积分方法的应用范围扩展至复杂几何形状与边界条件，包括多孔或穿孔体。
• 在多种构型（如楔形体入水、垂直平板入水及倾斜平板水下撞击）下，将方法与已知实验与数值数据进行验证。
• 实现在自相似公式中无需时间域积分即可准确计算水动力载荷与射流动力学。

提出的方法

• 采用伪时间推进法迭代求解全非线性边值问题，将自由表面视为移动边界。
• 引入修正速度势，以简化自由表面上动力与运动边界条件的施加。
• 在每次迭代中，使用常数面板离散化的边界积分方程法求解速度势的拉普拉斯方程。
• 在自由表面上施加狄利克雷条件，而在物体表面施加诺伊曼条件，通过边缘条件确保自由表面的切向脱离。
• 采用拉伸的自相似变量将时间依赖问题转化为稳态边值问题，降低计算复杂度。
• 在细射流区域应用浅水模型，以降低计算成本，同时保持精度。

实验结果

研究问题

• RQ1 如何构建一种全非线性迭代方法，以求解具有未知边界的自相似自由表面势流？
• RQ2 该方法能否准确捕捉水下入水与水下撞击问题中的细射流与复杂流动结构？
• RQ3 使用修正速度势在多大程度上提升了边界积分解法的稳定性与效率？
• RQ4 该方法在多大程度上可扩展以包含多孔或穿孔体边界？
• RQ5 与实验和数值基准相比，该方法的精度与计算效率如何？

主要发现

• 该方法成功捕捉了经典楔形体入水问题的自相似流场，结果与现有文献高度一致。
• 在垂直平板入水问题中，该方法预测的射流尖端速度约为平板速度的三倍，符合物理预期。
• 在10°倾斜平板水下撞击问题中，当V/U = 0.03时，射流根部以1.5倍平板水平速度移动，且压力峰值位于边缘下游0.4976倍平板速度处。
• 倾斜平板上的水动力载荷计算结果为F ≈ 0.156ρU³t，由压力分布的数值积分获得。
• 自由表面形成飞溅，内角约为61度，两侧分别与静水位成52°和67°倾角。
• 由于高效浅水模型与自适应离散化，该方法在高流速与高压梯度区域（尤其是细射流区域）仍保持高精度。