[论文解读] A Nonlinear ODE System for the Unsteady Hydrodynamic Force – A New Approach

Osama A. Marzouk|arXiv (Cornell University)|Oct 3, 2010

Fluid Dynamics Simulations and Interactions被引用 4

一句话总结

本文提出了一种新颖的两自由度非线性常微分方程系统，通过直接模拟随时间演变的合力系数（CT）及其方向（β），精确建模了圆柱体上的非定常总水动力力。该模型通过使用从模拟中获得的雷诺数相关参数，避免求解完整的纳维-斯托克斯偏微分方程，从而在计算成本极低的情况下，实现了在雷诺数Re = 100–500范围内对力的平均分量和振荡分量的高保真度预测。

ABSTRACT

We propose a reduced-order model for the instantaneous hydrodynamic force on a cylinder. The model consists of a system of two ordinary differential equations (ODEs), which can be integrated in time to yield very accurate histories of the resultant force and its direction. In contrast to several existing models, the proposed model considers the actual (total) hydrodynamic force rather than its<br> perpendicular or parallel projection (the lift and drag), and captures the complete force rather than the oscillatory part only. We study and provide descriptions of the relationship between the model parameters, evaluated utilizing results from numerical simulations, and the Reynolds number so that the model can be used at any arbitrary value within the considered range of 100 to 500 to provide accurate representation of the force without the need to perform time-consuming simulations and solving the partial differential equations (PDEs) governing the flow field.

研究动机与目标

• 开发一种降阶模型，直接模拟总水动力力（CT 和 β），而非其升力与阻力分量的投影。
• 消除对耗时的纳维-斯托克斯模拟及速度/压力场后处理的依赖。
• 通过参数化常微分方程系统，在Re = 100–500范围内任意雷诺数下实现精确的力预测。
• 为流固耦合研究提供一种具有物理解释性且可解析处理的模型。

提出的方法

• 构建两方程系统：一个具有三次非线性的自激振子方程用于建模力的方向β，一个受迫振子方程用于建模总力系数CT。
• 采用具有软化刚度（负三次项）的杜芬型方程来模拟β的动力学。
• 通过非线性阻尼和耦合系数将CT与β方程耦合。
• 利用Re = 100至500范围内的高保真数值模拟校准模型参数。
• 推导出每个参数与Re/100之间连续的多项式函数关系（R² > 0.9927），以实现任意Re值的插值。
• 在Re = 100、300和500下，将模型与完整的纳维-斯托克斯模拟进行验证，结果显示CT和β随时间的历史变化具有极佳的一致性。

实验结果

研究问题

• RQ1 是否能够通过降阶常微分方程系统，在不求解纳维-斯托克斯方程的前提下，准确捕捉圆柱体上总水动力力（大小与方向）的完整时间历史？
• RQ2 在尾迹不稳定性范围（Re = 100–500）内，模型参数（阻尼、刚度、耦合）如何随雷诺数变化？
• RQ3 该模型在多大程度上能够再现CT和β的关键频谱特征，例如CT的2St频率和β的基频St频率？
• RQ4 该模型是否保持了时间平均CT非零而平均β为零的物理特性？
• RQ5 该模型能否作为流固耦合研究中全尺寸CFD模拟的可靠、参数化替代方案？

主要发现

• 该模型在Re = 100、300和500下，对CT和β的时间历史具有高度一致的再现能力，与纳维-斯托克斯模拟结果几乎完全一致。
• 总力系数CT在2St处具有主导谐波，更高次谐波可忽略不计；而β主要为基频St，更高次谐波（如a3β）小一个数量级。
• 线性频率ω在Re = 400处达到峰值后趋于稳定，而杜芬项系数始终保持负值（软化特性），并在Re = 350前持续增大。
• 阻尼与耦合系数在Re = 100至150间急剧减小，随后呈现振荡行为，且参数-雷诺数关系的多项式拟合决定系数R²高达0.9927–0.9996。
• 即使在高涡脱频率条件下（如Re = 500），该模型仍保持高精度，CT中仅存在微小偏差，归因于a3β效应的增强。
• 所提出的模型可直接、高效地模拟完整力矢量，无需分解为升力与阻力，从而保持物理保真度，并可与结构动力学实现耦合。