船舶力学

对比结果	本文结果	文献1^[19]	文献2^[20]
St	0.188	0.192	0.190

对比结果	本文结果	文献1^[19]	文献2^[20]
St	0.188	0.192	0.190

参数类型	管束数	管束直径	管束分布圆直径	节距
符号/单位	N_t	D_t /m	D_d /m	P /m
参数值	7	0.032	0.128	0.048

参数类型	管束数	管束直径	管束分布圆直径	节距
符号/单位	N_t	D_t /m	D_d /m	P /m
参数值	7	0.032	0.128	0.048

工况	网格方案		评判标准参数
Case 1	120	0.00084	179436	0.1971
Case 2	160	0.00063	220685	0.2383
Case 3	200	0.00050	271341	0.2532
Case 4	200	0.00050	335071	0.2536

工况	网格方案		评判标准参数
Case 1	120	0.00084	179436	0.1971
Case 2	160	0.00063	220685	0.2383
Case 3	200	0.00050	271341	0.2532
Case 4	200	0.00050	335071	0.2536

模态参数	Mode 0	Mode 1	Mode 2	Mode 3	Mode 4	Mode 5
振幅	2656591	34533	4552	1179	678	461
增长率/（×10⁻⁵）	0.00	−1.82	−8.98	−22.70	−22.17	−24.32
频率/Hz	0.00	16.59	33.19	49.78	66.37	82.97

模态参数	Mode 0	Mode 1	Mode 2	Mode 3	Mode 4	Mode 5
振幅	2656591	34533	4552	1179	678	461
增长率/（×10⁻⁵）	0.00	−1.82	−8.98	−22.70	−22.17	−24.32
频率/Hz	0.00	16.59	33.19	49.78	66.37	82.97

错列管束绕流的动力学模态分解参数研究

PDF下载

韩旺 ¹ , 曹贻鹏 ¹ , 刘晨 ¹ , 孙凯浪 ² , 孙玉东 ² , 赵晓臣 ¹

船舶力学 | 流体力学 2026,30(3): 387-398

收起

船舶力学 |流体力学 2026 , 30 (3) : 387 -398

错列管束绕流的动力学模态分解参数研究

全屏

韩旺¹, 曹贻鹏¹, 刘晨¹, 孙凯浪², 孙玉东², 赵晓臣¹

作者信息

^1.哈尔滨工程大学动力与能源工程学院，哈尔滨 150001

^2.中国船舶科学研究中心，江苏无锡 214082

通讯作者:

曹贻鹏（1980–），男，副教授，通讯作者，E-mail：yipengcao06@163.com

作者简介:

韩　旺（1997–），男，博士研究生，E-mail：han_wang@hrbeu.edu.cn

The study of dynamic mode decomposition parameters for the flow around staggered tube bundles

Wang HAN¹, Yi-peng CAO¹, Chen LIU¹, Kai-lang SUN², Yu-dong SUN², Xiao-chen ZHAO¹

Affiliations

^1.Harbin Engineering University, Harbin 150001, China

^2.China Ship Scientific Research Center, Wuxi 214082, China

出版时间: 2026-03-15 doi: 10.3969/j.issn.1007-7294.2026.03.005

文章导航

摘要

收起

本文对错列换热管束绕流场进行了非定常数值仿真，使用动力学模态分解（Dynamic Mode Decomposition, DMD）方法对数值计算得到的管束绕流涡量数据进行了分析及流场重构，研究了快照数、奇异值的截断阶数对流场重构的影响。结果表明：基于DMD方法能够准确有效地获取错列管束绕流非定常流场中涡量数据的模态结构和对应模态的频率、增长率，基于主要模态重构的流场与数值结果吻合较好，验证了DMD方法在错列管束绕流分析中的有效性；根据无量纲周期数，提出了最小快照数准则，能够直接得到有效的低误差快照数，极大减少了计算成本；从涵盖能量的程度出发，提出了最优奇异值截断阶数准则，有效避免过大的奇异值截断阶数而带来的误差。该研究可为复杂换热器管束绕流场的DMD研究提供参数选取依据。

关键词

错列管束 / 动力学模态分解 / 流场重构 / 快照数 / 奇异值截断阶数

Abstract

收起

Unsteady numerical simulations were systematically conducted to investigate the transient flow characteristics around a staggered heat exchanger tube bundle. The vorticity distribution derived from numerical simulations was analyzed using Dynamic Mode Decomposition (DMD) for modal extraction and flow field reconstruction. The effects of the number of snapshots and singular value truncation order on the accuracy of reconstructed flow fields were studied. The results indicate that DMD enables accurate extraction of both vortical mode structures and their associated characteristic parameters, including frequencies and growth rates, from transient flow fields surrounding a staggered tube bundle. Flow field reconstructions obtained from dominant coherent modes exhibit strong agreement with high-fidelity numerical benchmarks. This confirms the effectiveness of the DMD methodology in resolving vortex-driven flow interactions in staggered tube bundle configurations. A minimum snapshot criterion is proposed based on the dimensionless period number, allowing direct determination of an appropriate number of snapshots with minimal computational error, which in turn substantially reduces computational costs. In terms of the energy retention ratio, a criterion for determining the optimal singular value truncation order is established, alleviating errors caused by overly aggressive truncation thresholds. This study provides theoretical insights into parameters selection for DMD analysis of complex flow fields around heat exchanger tube bundles.

Key words

staggered tube bundle / Dynamic Mode Decomposition / flow reconstruction / number of snapshots / singular value truncation order

引用本文

韩旺, 曹贻鹏, 刘晨, 孙凯浪, 孙玉东, 赵晓臣. 错列管束绕流的动力学模态分解参数研究. 船舶力学, 2026 , 30 (3) : 387 -398 . DOI: 10.3969/j.issn.1007-7294.2026.03.005

Wang HAN, Yi-peng CAO, Chen LIU, Kai-lang SUN, Yu-dong SUN, Xiao-chen ZHAO. The study of dynamic mode decomposition parameters for the flow around staggered tube bundles[J]. Journal of Ship Mechanics, 2026 , 30 (3) : 387 -398 . DOI: 10.3969/j.issn.1007-7294.2026.03.005

正文

收起

0　引　言

收起

管壳式换热器在船舶动力系统等领域中有着极为广泛的应用，其壳程内部的换热管按一定结构、尺寸和排列形式组成了不同类型的换热管束，壳程介质在流经错列管束时，会出现剪切层分离、失稳及转捩等现象，并在尾流中形成交替脱落的旋涡^[1]。管束在湍流涡结构影响下的噪声和振动问题长期存在，不仅会造成换热器结构破坏、性能降低，还将消耗大量维护费用。因此，研究换热器壳程错列管束绕流场特性至关重要。

随着高精度计算流体力学的迅猛发展，已能够实现换热器壳程管束绕流的数值模拟研究。然而，面对低分辨率且高维度的庞大流场数据，如何从复杂的非线性管束绕流场数据中提取低维空间下的主要流动特征，仍是一个需要解决的问题^[2]。本征正交分解（Proper Orthogonal Decomposition, POD）和动力学模态分解（Dynamic Mode Decomposition, DMD）是复杂流动分析中建立降阶模型较为典型的两种方法^[3]。POD将高阶流场矩阵投影到若干正交的低阶向量上，实现了空间层面上对流场信息的分析，但该方法无法在时间层面上对流场的信息进行分析，可能出现模态重叠的现象^[4]。DMD方法对时空耦合结构数据建立降阶模型（Reduced Order Model, ROM），并结合主要模态和频率等信息对流场结构进行分析。Schmid^[5]首先提出了基于相似矩阵的DMD，并将其作为系统线性稳定性的分析方法。Rowley^[6]又提出了基于伴随矩阵的DMD，将其视为计算近似Koopman分解的特定数值算法。之后，学者们将DMD方法推广到了不同的流动分析中。例如，Sakai^[7]、叶坤^[8]、Bai^[9]、孙婉荣^[10]、袁猛^[11]将DMD分别应用于单圆柱与双圆柱绕流，分析其稳定性以及尾流流态的流动特性；韩忠华^[12]、Jonathan^[13]、Sampath^[14]、Seena^[15]、Rowley^[16]、丁杰^[17]将DMD分别应用于机翼等钝体绕流及其他流动，验证了DMD分析的有效性并获得了其流场结构特征，上述研究充分说明了DMD方法对流动问题分析的可靠性。学者们还针对DMD参数对流场重构误差的影响开展了研究，如快照数、时间步长和奇异值截断阶数等，得出了快照数和奇异值截断阶数对重构误差影响显著的结论。

相较于单圆柱和双圆柱绕流，错列管束绕流的流动机制更为复杂。然而，在现有的管束绕流研究中，动态模式分解方法的适用性研究仍存在不足。此外，现有研究通常采用枚举不同参数值的方式，通过对比不同参数设置下的重构效果来确定各参数的最优值。尽管这种方法能够通过多次试验实现较低的重构误差，但其对操作者经验的依赖性较高，缺乏统一的动态模式分解参数选取标准，且尚未针对错列管束扰流问题建立专门的DMD参数准则。基于此，本文开展错列管束绕流的动力学模态分解研究，在验证动力学模态分解错列管束绕流问题的适用性的同时，探究不同快照数对流场重构精度的影响，讨论不同奇异值截断阶数下的预测效果，提出满足错列管束绕流低误差的最小快照数准则与最优奇异值截断阶数准则。

1　计算方法

收起

1.1　流场控制方程

本文采用SST k-ω湍流模型来求解流体控制方程，其中k和ω的输运方程为

$ \frac{\partial }{\partial t}\left(\rho k\right)+\frac{\partial }{\partial {x}_{i}}\left(\rho k{u}_{i}\right)={P}_{k}-{\beta }^{*}\rho \omega k+\frac{\partial }{\partial {x}_{j}}\left[\left(\mu +{\sigma }_{k}{\mu }_{t}\right)\frac{\partial k}{{x}_{j}}\right] $

(1)

$ \frac{\partial }{\partial t}\left(\rho \omega \right)+\frac{\partial }{\partial {x}_{j}}\left(\rho k{u}_{j}\right)=\frac{\partial }{\partial {x}_{j}}\left[\left(\mu +{\sigma }_{\omega }{\mu }_{t}\right)\frac{\partial \omega }{{x}_{j}}\right]+2\left(1-{F}_{1}\right)\frac{\rho {\sigma }_{\omega 2}}{\omega }\frac{\partial k}{\partial {x}_{j}}\frac{\partial \omega }{\partial {x}_{j}}+\frac{\gamma }{{v}_{t}}P-\beta \rho {\omega }^{2} $

(2)

在方程（1）、（2）中，湍动能生成项

${P}_{k} $

、湍流粘度

${\mu }_{t} $

、混合函数

${F}_{1} $

、

${F}_{2} $

、以及应变率张量

$S $

的表达式如下

$ \begin{cases} {P}_{k}=\min \left(P,10{\beta }^{*}\rho k\omega \right)\\{\mu }_{t}=\dfrac{\rho {a}_{1}k}{\max \left({a}_{1}\omega ,S{F}_{2}\right)}\\{F}_{1}=\tanh \left(\arg _{1}^{4}\right),{\arg }_{1}=\min \left\{\max \left(\dfrac{\sqrt{k}}{{\beta }^{*}\omega y},\dfrac{500\mu }{\rho {y}^{2}\omega }\right),\dfrac{4\rho {\sigma }_{\omega 2}k}{C{D}_{k\omega }{y}^{2}}\right\}，C{D}_{k\omega }=\max \left(2\rho {\sigma }_{\omega 2}\dfrac{1}{\omega }\dfrac{\partial k}{\partial {x}_{j}}\dfrac{\partial \omega }{\partial {x}_{j}},{10}^{-10}\right)\\{F}_{2}=\tanh \left(\arg _{2}^{2}\right),{\arg }_{2}=\max \left(\dfrac{2\sqrt{k}}{0.09\omega y},\dfrac{500\mu }{\rho {y}^{2}\omega }\right)\\S=\sqrt{2{S}_{ij}{S}_{ij}},{S}_{ij}=\dfrac{1}{2}\left(\dfrac{\partial {u}_{i}}{\partial {x}_{j}}-\dfrac{\partial {u}_{j}}{\partial {x}_{i}}\right)\end{cases} $

(3)

式中：

$y $

为距最近壁面的法向近距离，

${\beta }^{*} $

、

${\beta }_{1} $

、

${\beta }_{2} $

、

${a}_{1} $

、

${\sigma }_{\omega 1} $

、

${\sigma }_{\omega 2} $

、

$ {\sigma }_{k1} $

、

${\sigma }_{k2} $

均为模型系数。

需要说明的是，

${F}_{1} $

并非方程（1）、（2）的未知量，而是SST k-ω模型的混合函数，其一用于控制ω方程交叉扩散项中的系数

$ 2\left(1-{F}_{1}\right) $

，其二用于模型常数的加权混合。式（2）中的

$P $

为由平均速度梯度计算得到的原始湍动能生成项，定义为

$P=2{\mu }_{t}{S}^{2} $

，为避免在强剪切区域的数值不稳定性，对湍动能方程中的生成项进行约束，即

${P}_{k}=\min \left(P,10{\beta }^{*}\rho k\omega \right) $

。式（2）中的

${v}_{t} $

为湍流运动粘度，定义为

${v}_{t}={\mu }_{t}/\rho $

。

1.2　动力学模态分解方法

DMD方法将流动过程分解为若干流体模态及其特征参数，通过选取合适的模态建立ROM来近似还原物理场特征。获取动力学ROM可以采取两种方法，分别是伴随矩阵法和相似矩阵法，而后者比前者具有更好的鲁棒性^[18]。因此，本文采用相似矩阵法对错列管束绕流场进行DMD求解。

设在初始时刻流场中某物理量的数值为一个列向量

$ {x}_{1} $

，其为初始时刻的快照，列向量的元素中包含了系统不同空间位置的该物理量的值，之后每隔相等的时间步长

$ {\Delta }t $

进行一次数据采样，直到第

$ k\times {\Delta }t $

时刻停止，得到共计

$ k $

组快照构成的数据集合{

$ {x}_{1},{x}_{2},{x}_{3},\cdots ,{x}_{k} $

}。

基于Koopman算子可将非线性的数据快照近似为线性动力学系统，即快照

$ {x}_{i+1} $

可通过快照

$ {x}_{i} $

的线性映射表示为

$ {x}_{i+1}={\boldsymbol{A}}{x}_{i} $

，

$ {\boldsymbol{A}} $

包含了系统的时间演化信息。对于1到k组流场快照，有如下关系

$ {\boldsymbol{X}}_{2}^{\left(k\right)}={\boldsymbol{AX}}_{1}^{\left(k-1\right)} $

(4)

式中：

$ {\boldsymbol{X}}_{1}^{\left(k-1\right)} $

$ {x}_{1},{x}_{2},{x}_{3},\cdots ,{x}_{k-1} $

]，

$ {\boldsymbol{X}}_{2}^{\left(k\right)} $

$ {x}_{2},{x}_{3},{x}_{4},\cdots ,{x}_{k} $

]。

由于矩阵

$ {\boldsymbol{A}} $

的维度很高且矩阵

$ {\boldsymbol{X}}_{1}^{\left(k-1\right)} $

、

$ {\boldsymbol{X}}_{2}^{\left(k\right)} $

并不是方阵，因此无法直接求解，需要寻找

$ {\boldsymbol{A}} $

的近似矩阵

$ {\boldsymbol{A}}' $

来获取系统的动力学特征。首先，对矩阵

$ {\boldsymbol{X}}_{1}^{\left(k-1\right)} $

进行奇异值分解

$ {\boldsymbol{X}}_{1}^{\left(k-1\right)}={\boldsymbol{USV}}^{H} $

(5)

式中：

$ {\boldsymbol{U}} $

为左奇异矩阵，

$ {\boldsymbol{S}} $

为奇异值矩阵，

$ {\boldsymbol{V}} $

为右奇异矩阵。用

$ {\boldsymbol{C}} $

表示矩阵，在右上角加数字表示矩阵的维度，则

$ {\boldsymbol{U}}\in {{\boldsymbol{C}}}^{n\times r} $

，

$ {\boldsymbol{S}}\in {{\boldsymbol{C}}}^{r\times r} $

，

$ {\boldsymbol{V}}\in {{\boldsymbol{C}}}^{\left(k-1\right)\times r} $

。

$ r $

是指奇异值分解之后截断的阶数。

这样矩阵

$ {\boldsymbol{A}} $

便可用矩阵

$ {\boldsymbol{X}}_{1}^{\left(k-1\right)} $

的左奇异矩阵

$ {\boldsymbol{U}} $

表示

$ {\boldsymbol{A}}\approx {\boldsymbol{UA}}'{{\boldsymbol{U}}}^{H} $

(6)

将式（5）和（6）代入式（4）中，可得到近似矩阵

$ {\boldsymbol{A}}' $

的表达式为

$ {\boldsymbol{A}}'={{\boldsymbol{U}}}^{H}{\boldsymbol{X}}_{2}^{\left(k\right)}{\boldsymbol{VS}}^{-1} $

(7)

对近似矩阵

$ {\boldsymbol{A}}' $

进行特征分解，得到其特征值和特征向量

$ {\boldsymbol{A}}'{\boldsymbol{W}}={\boldsymbol{W}}{{\boldsymbol{\varLambda}} } $

(8)

式中：

$ {\boldsymbol{W}} $

为

$ {\boldsymbol{A}}' $

的特征向量构成的矩阵，

$ {{\boldsymbol{\varLambda}} } $

为对角矩阵，其对角线元素由

$ {\boldsymbol{A}}' $

的特征值构成。使用

$ {\boldsymbol{A}}' $

的特征值和特征向量可以得到DMD模态

$ {{{\varPhi}} }={{\boldsymbol{X}}}_{2}^{\left(\textit{k}\right)}{\boldsymbol{VS}}^{-1}{\boldsymbol{W}} $

(9)

式中：

$ {\varPhi } $

的列向量

$ {\varphi }_{i} $

即为DMD模态，模态对应的特征值就是

$ {{\boldsymbol{\varLambda}} } $

对角线上的元素

$ {\lambda }_{i} $

，

$ \ln \left({\lambda }_{i}\right)/{\Delta }t $

的实部是模态的增长率，虚部是模态的频率。根据上述结果，便可重构

$ t+{\Delta }t $

时刻的数据快照

$ \textit{x}\left(\textit{t}+{\Delta }t\right)=\sum\limits_{i=1}^{n}{\varphi }_{i}\exp \left(\ln {\lambda }_{i}t/{\Delta }t\right){b}_{i}={\Phi }{{{\boldsymbol{\varLambda}} }}^{t/{\Delta }t}b $

(10)

式中：

$ {b}_{i} $

为第

$ i $

阶模态的振幅，

$ b $

为由元素

$ {b}_{i} $

构成的列向量，

$ b $

可以根据

$ {\varPhi }b={{\boldsymbol{X}}}_{1} $

得到。

2　非定常流场数值计算

收起

2.1　验证模型

在对错列管束绕流模型进行数值计算之前，有必要对数值计算方法及参数进行验证。本文选择了一个二维圆柱绕流模型进行验证^[19]。圆直径为5.12 mm，Re=200，流动介质为水。计算域网格划分为结构化网格，确保圆柱边界层网格y⁺＜1。建立的二维圆柱流计算网格如图1所示。

对二维圆柱绕流非定常流动进行数值计算，并取圆柱斯特劳哈尔数St与文献结果进行对比，如表1所示。本文结果与文献结果的最大误差不足3%，充分说明了数值方法和计算参数的有效性。

2.2　错列管束模型

工业领域中实际应用的管壳式换热器管束数量众多且内部结构复杂，为合理地分析管束绕流特性，忽略壳程支撑结构，将管束绕流场简化为由7根转置正三角形排列的换热管束组成的最小“单元”。计算模型取为二维矩形流域，最小管束“单元”位于流域内部，流体从左侧流入、右侧流出，二维错列管束绕流模型如图2所示，表2显示了模型详细的几何参数。

对二维错列管束模型划分结构化网格。其中，管束壁面处第一层网格高度设置为1.0×10⁻⁵ m，增长率为1.1，确保y⁺<1。考虑网格密度对数值计算的影响，共讨论4组网格密度下计算结果的差异，相应网格参数列于表3中。网格密度的控制参数包括管束圆周向节点数

$ {N}_{{\mathrm{C}}} $

、无因次参数

$ {\Delta }_{{\mathrm{C}}}/{N}_{{\mathrm{C}}} $

（

$ {\Delta }_{{\mathrm{C}}} $

为管束的周向长度）以及总单元数

$ N $

^[21]。以中心管束上的斯特劳哈尔数St作为评判网格无关性的标准参数，不同网格密度算例的计算结果如表3所示。结果表明，Case 3、Case 4的St基本一致，相对误差在1%以内，出于节约计算资源的目的，本文选取Case 3（

$ {N}_{{\mathrm{C}}} $

=200，

$ {\Delta }_{{\mathrm{C}}}/{N}_{{\mathrm{C}}} $

=0.00050，

$ N $

=271341）的网格参数用于该研究的数值模拟工作。

模型入口定义为速度边界，流速为3 m/s；出口设置为压力边界，表压为0 Pa；流场介质为水，对应的物性参数均按照0 ℃、一个大气压力下选取。为加快流场计算进程，先使用Standard

$ k-\varepsilon $

湍流模型对二维错列管束模型进行定常流动计算，待参数波动稳定后，以定常流场作为流场初始值进行非定常流动计算。压力速度耦合选用SIMPLE算法，压力选用二阶离散格式，动量、湍动能及湍流耗散率均采用二阶迎风离散格式。瞬态计算时间步长为

$ {\Delta }t $

=0.01 s，总计算时长

$ T $

=60 s。

在计算的同时监测7根管束壁面上的升力和阻力系数，为减少瞬态时间步长的影响，横坐标取为无量纲时间步长（

$ {V}_{0}T/{D}_{{\mathrm{t}}} $

），如图3所示；稳态、瞬态管束绕流场的涡量分布如图4所示。从图3可以看出，经过定常流场计算后，非定常流动很快便达到了稳定阶段，此时各管束升力、阻力系数随无量纲时间步长增加，在经历不稳定阶段后出现了周期性的变化特征（

$ {V}_{0}T/{D}_{{\mathrm{t}}}\geq $

937.5）；从图4可以看出，在当前圆柱间距和参数设置下，流场的两个阶段中管束下游产生的涡街方向均一致，属于同步同相模式。其中，在非定常流动中，尾流中出现了交替产生的涡，这一特征与各管束升力、阻力系数波动趋势一致。

3　DMD计算分析

收起

3.1　动力学模态分析

在图3中的动态稳定阶段，取937.5~2237.5无量纲时间步跨度的流场涡量数据作为DMD分析的快照数据，共计1399个快照。图5为各模态特征值的实部和虚部在单位圆的分布，位于特征圆上的点共有27个（红色标识），分别对应13对共轭模态和一个主导模态，其均为周期性模态，位于特征圆外部的均为不稳定模态，其余特征值均位于单位元内部，属于稳定模态。但是受计算机舍入误差的影响，到圆心距离严格为1的点并不存在，因此基于特征圆的分析方法仅能定性分析模态稳定性^[11]。

根据

$ {\boldsymbol{X}}_{1}^{\left(k-1\right)} $

奇异值的变化，使用硬阈值作为判据以截断奇异值^[22]。在保留前19阶奇异值后，按照振幅从大到小对前6阶模态排序，其模态参数如表4所示、振幅和增长率的关系如图6所示。

从表4、图6中可以看出，Mode 0点对应的模态系数幅值最大，模态增长率为0，因此该模态的幅值基本不变且对各时刻流场的作用基本相同，因此对应的是流场的平均状态。Mode 1点和Mode 2点对应的模态增长率较小且为负，幅值没有继续增大的趋势，反应了流场随时间的变化特征。Mode 3点、Mode 4点和Mode 5点对应的模态增长率较大且为负，是稳定模态。以上计算完成了对二维错列管束绕流场的动力学模态分解，建立了错列管束绕流场的涡量降阶模型，在此基础上，研究建立的降阶模型在涡量场重构中的效果，重构过程如图7所示，重构错列管束涡量场与真实错列管束涡量场对比如图8所示。

从图8中可以看出，对于所取的快照空间时刻内的涡量场，使用DMD方法已经能够捕捉错列管束绕流中关键的流场特征并进行流场重构，重构流场与原始流场高度吻合，表明重构效果良好，丢失的流场特征较少。但由于实际的流场涡量特征更为复杂，所用1399个快照和截断前19阶模态得到的近似涡量场仍与真实涡量场间存在5%的偏差。

3.2　输入快照数对流场重构效果的影响

本文提出采用无量纲周期数

$ {N}_{{\mathrm{T}}} $

对快照数的影响进行衡量。对于类圆柱绕流，其流动具有特征频率，因此使用斯特劳哈尔数St对其流动进行描述

$ {S}{t}=fL/{V} $

(11)

式中：

$ f $

为涡脱落频率，

$ L $

为特征长度，

$ {V} $

为流体速度。其中，

$ f $

又可以写成

$ f=1/T $

(12)

式中：

$ T $

为涡脱落周期。在DMD快照取样时，应当至少包含一个完整周期性的数据特征，即快照数量需要满足

$ N\geq T/{\Delta }t $

(13)

式中：

$ N $

为快照数，

$ {\Delta }t $

为流场瞬态时间步长。考虑到类圆柱绕流结构的相似性，以St作为相似性准则，把错列管束绕流等效为圆柱绕流，由式（11）、（12）、（13）便可以得到关于错列管束绕流的最小快照数

$ {N}_{{\mathrm{T}}}={C}_{1}{D}_{{\mathrm{d}}}/\left({C}_{2}{S}{t}{V}{\Delta }t\right) $

(14)

式中：

$ {N}_{{\mathrm{T}}} $

是无量纲周期快照数，

$ {C}_{1} $

是管束分布系数，

$ {C}_{1}={A}_{{\mathrm{t}}}/{A}_{{\mathrm{d}}} $

，

$ {A}_{{\mathrm{t}}} $

、

$ {A}_{{\mathrm{d}}} $

是管束和分布圆面积，

$ {C}_{2} $

是频率系数，

$ {C}_{2}={f}_{{\mathrm{feature}}}/\left(n{f}_{\max }\right) $

，

$ {f}_{{\mathrm{feature}}} $

是特征频率，

$ {f}_{\max } $

是最大分析频率，

$ n $

是频率放大系数（

$ n $

=1,2,3,4,5,6），

$ {V} $

是入口流速。本文

$ {C}_{1}=0.4375 $

，

$ {C}_{2} $

考虑了不同频率放大系数，依次取为：0.110、0.055、0.036、0.027、0.022和0.018，对应快照数

$ {N}_{{\mathrm{T}}} $

分别为533、1066、1599、2132、2665和3198。

同时，为了更全面地考虑所有时刻下全局尺度上的重构误差，采用均方根误差

$ {x}_{{\mathrm{RMS}}} $

对DMD模态重构的涡量场误差进行全面评估，

$ {x}_{{\mathrm{RMS}}} $

定义为

$ {x}_{{\mathrm{RMS}}}=\sqrt{\frac{1}{N}\sum\limits_{i=1}^{N}{\left| {x}_{{\mathrm{real}}}\left(i\right)-{x}_{{\mathrm{DMD}}}\left(i\right)\right| }^{2}} $

(15)

式中：

$ N $

是快照数，

$ {x}_{{\mathrm{real}}}\left(i\right) $

与

$ {x}_{{\mathrm{DMD}}}\left(i\right) $

分别是真实流场与DMD模态重构的第

$ i $

个涡量场快照。

从图3中动态稳定阶段第937.5个无量纲时间步开始，取不同

$ {N}_{{\mathrm{T}}} $

下管束绕流场涡量数据作为输入快照，图9为不同快照数下重构涡量场与真实涡量场的误差对比。

从图中可以看出，重构涡量场与真实涡量场的误差整体较小，主要分布在管束7的正下游及流域的尾流位置，其中管束7的正下游处是最大误差位置，DMD可以较精确地捕捉低频大尺度的旋涡，但在捕捉高频小尺度涡上存在一定误差。随快照数的增加，重构涡量场与真实涡量场的差异逐渐减小，当快照数不少于

$ {N}_{{\mathrm{T}}}=1599 $

（

$ n $

=3）时，整个流场内的均方根误差均小于10 s⁻¹，与真实涡量场之间的误差比不足3.3%，说明此时快照数对于DMD重构结果的影响可以忽略不计。

综合以上结果，将

$ n $

=3代入式（14），得到了适用于错列管束绕流问题最小快照数

$ {N}_{\min } $

的表达式为

$ {N}_{\min }=\frac{3{A}_{{\mathrm{t}}}{f}_{\max }{D}_{{\mathrm{d}}}}{{A}_{{\mathrm{d}}}{f}_{{\mathrm{feature}}}{S}{t}{V}{\Delta }t} $

(16)

以快照数作为横坐标，以重构涡量场的最大均方根误差作为纵坐标进行作图，结果如图10所示。从图中可以直观地看出，涡量最大均方根误差随快照数增加逐渐减小，减小速率随快照数增加而变小，在最小快照数

$ {N}_{\min } $

位置涡量最大均方根误差不足10 s⁻¹，之后随快照数增加，最大均方根误差基本保持不变。说明当满足最小快照数

$ {N}_{\min } $

时，快照数对建立的降阶模型重构误差的影响可以忽略不计。

根据不同快照数下涡量场重构均方根误差的计算与分析，本文提出了最小快照数准则的表达式。

3.3　奇异值截断阶数对流场重构效果的影响

从能量的角度来看，奇异值截断的阶数决定了重构流场的能量保留程度。通过保留能量占比大的奇异值，可以保留流场的大部分能量，减少计算量和存储需求，同时控制重构误差。因此，在最小快照数

$ {N}_{\min } $

（

$ {N}_{{\mathrm{T}}}=1599 $

）下，本文从建立的ROM模型与原始流场总能量占比的角度出发，讨论奇异值的截断阶数对错列管束绕流场重构精度的影响。

在降阶模型与原始模型的一致性研究方面，常以95%能量保留和99%能量保留作为能量占比规则，95%的能量保留已被证明能够在保证精度的同时，显著降低模型的复杂性，而99%的能量保留可以更好地保证降阶后的模型与原始模型的相似度。因此为保证重构流场精度，本文分别取截断模态总能量在原始流场总能量的占比为95%、99%、99.50%和99.99%时的奇异值的截断阶数进行讨论。截断模态总能量在原始流场总能量的表达式如式（17）所示。

$ {E}_{{\mathrm{DMD}}}=\sum\limits_{i=1}^{r}{\varphi }_{i}{b}_{i}{Q}_{i}, {E}_{{\mathrm{real}}}=\sum\limits_{i=1}^{N}{\varphi }_{i}{b}_{i}{Q}_{\mathrm{i}}, {r}=\mathrm{ar}\mathrm{gmi}\mathrm{n}\left\{\left({E}_{{\mathrm{DMD}}}\right)/\left({E}_{{\mathrm{real}}}\right)\geq C\right\} $

(17)

式中：

$ {E}_{{{\mathrm{DMD}}{\mathrm{}}}} $

是降阶模型总能量，

$ {E}_{\text{real}} $

是原始流场总能量，以DMD分解得到的全部模态能量总和表示，

$ Q $

是由特征值组成的范德蒙矩阵，

$ Q=\left[{\boldsymbol{\lambda }}^{0},{\boldsymbol{\lambda }}^{1},{\boldsymbol{\lambda }}^{2},\cdots ,{\boldsymbol{\lambda }}^{N-2}\right] $

，

$ \boldsymbol{\lambda } $

是由特征值组成的列向量，

$ \boldsymbol{\lambda }\mathbf{=}{\left[{\lambda }_{1},{\lambda }_{2},{\lambda }_{3},\cdots ,{\lambda }_{r}\right]}^{{\mathrm{T}}} $

。

将

$ {E}_{{\mathrm{DMD}}}/{E}_{{\mathrm{real}}} $

定义为能量占比。通过式（17）可求得，当能量占比为95%、99%、99.50%和99.99%时，相应奇异值的截断阶数分别为11、21、85、217。下面对不同截断阶数下流场涡量值差异进行分析。图11为位于管束4正下游方向，距离管束4圆心0.25 m处的中心采样点的涡量值变化的重构值与真实值对比。结合最大误差分布位置，图12为位于管束7正下游方向，距离管束7圆心0.05 m处的最大误差采样点的涡量值变化的重构值与真实值的对比。

从图11、图12可以看出，当能量占比为95%（r=11）时，两采样点处涡量重构值与真实值差距较大，基于DMD 方法的重构结果与真实值的变化周期基本一致，但是同步性较差且幅值存在偏差，尤其在最大误差采样点处的重构值幅值差异极大。当能量占比为99%（r=21）时，虽然两采样点处涡量重构结果在幅值的准确性上差异较小，且两曲线的变化规律完全同步，此时通过DMD方法可以用最少计算资源实现较好的流场重构效果。当能量占比为99.5%（r=85）时，虽然用于重构流场的模态个数增加，但增加的模态中同时包含了一部分不稳定模态，对重构结果引入了不必要的干扰，导致重构涡量曲线和真实结果的幅值与同步性偏差较大。当能量占比为99.99%（r=217）时，此时所建立的错列管束ROM模型与真实流场基本无异，两采样点处的重构值与原始值在1500个快照时刻以内都高度吻合，相比于能量占比为99%（r=21）时流场重构精度更高，但是将消耗更大的计算资源。

根据不同奇异值截断阶数下涡量重构值的计算与分析，本文提出，最优奇异值截断阶数应满足能量占比为99%的条件。

4　结　论

收起

本文利用DMD方法对数值计算得到的错列管束绕流场的涡量数据进行了分析，建立了降阶模型并对流场进行了重构，验证了DMD方法对错列管束绕流问题的适用性，在此基础上研究了不同快照数、不同奇异值截断阶数对DMD重构误差的影响，进而得出以下主要结论：

（1）基于DMD方法，针对错列管束绕流场，通过提取流场主要特征和各主导模态，进而重构流场。得到的结果与数值结果吻合较好，验证了本文工作的有效性。

（2）本文使用无量纲周期数来判别快照数量，进而提出了最小快照数准则。通过该准则，能够直接得到有效的低误差快照数，这避免在计算中可能设置的冗余快照数，以及为寻找合适的低误差快照数而开展的枚举试算工作，进而极大地减少了计算成本。

（3）本文从能量占比角度出发，根据降阶模型能量占流场总能量的比值来得到最优奇异值截断阶数。进而通过核对涡量重构值的准确性得到了最优奇异值截断阶数准则。这能够有效避免过大的奇异值截断阶数带来的误差。

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

杜晓庆, 林芝强, 吴葛菲. 串列双圆柱绕流的模态特征及流场重构[J]. 哈尔滨工业大学学报, 2022, 54(10): 49‒57.

Du X Q, Du Z Q, Wu G F. Mode characteristics and reconstruction of flow field around two tandem circular cylinders[J]. Journal of Harbin Institute of Technology, 2022, 54(10): 49‒57. (in Chinese)

黄瑾, 明平剑, 孙婉荣. 流体横掠管束的动模态分析研究[J]. 工程热物理学报, 2023, 44(8): 2279‒2284.

Huang J, Ming P J, Sun W R. Dynamic modal analysis of fluid sweeping rods bundle[J]. Journal of Engineering Thermophysics, 2023, 44(8): 2279‒2284. (in Chinese)

Taira K, Hemati M S, Brunton S L, et al. Modal analysis of fluid flows: Applications and outlook[J]. American Institute of Aeronautics and Astronautics Journal, 2020, 58(3): 998.

Berkooz G, Holmes P, Lumley J L. The proper orthogonal decomposition in the analysis of turbulent flows[J]. Annual Review of Fluid Mechanics, 1993, 25: 539‒575.

Schmid P J. Dynamic mode decomposition of numerical and experimental data[J]. Journal of Fluid Mechanics, 2010, 656: 5‒28.

Nathan J K, Rowley C W, Steven L B, et al. On dynamic mode decomposition: theory and applications[J]. Journal of Computational Dynamics, 2015, 1(2): 391‒421.

Sakai M, Sunada Y, Imamura T, et al. Experimental and numerical studies on flow behind a circular cylinder based on POD and DMD[J]. Transactions of the Japan Society for Aeronautical and Space Sciences, 2015, 58(2): 100.

叶坤, 武洁, 叶正寅, 等. 动力学模态分解和本征正交分解对圆柱绕流稳定性的分析[J]. 西北工业大学学报, 2017, 35(4): 599.

Ye K, Wu J, Ye Z Y, et al. Analysis circular cylinder flow using dynamic mode and proper orthogonal decomposition[J]. Chinese Journal of Northwestern Polytechnical University, 2017, 35(4): 599. (in Chinese)

Bai X D, Zhang W, Wang Y. Deflected oscillatory wake pattern behind two side-by-side circular cylinders[J]. Ocean Engineering, 2020, 197: 106847.

孙婉荣, 明平剑, 刘勇. 二维圆柱绕流的动模态分析研究[J]. 工程热物理学报, 2021, 42(5): 1275‒1281.

Sun W R, Ming P J, Liu Y. Study on dynamic mode analysis of two-dimensional cylinder flow[J]. Journal of Engineering Thermophysics, 2021, 42(5): 1275‒1281. (in Chinese)

袁猛, 张新玉, 柳贡民, 等. 并列双圆柱绕流的动力学模态分解[J]. 船舶力学, 2021, 25(5): 527‒534.

Yuan M, Zhang X Y, Liu G M, et al. Dynamic mode decomposition of flow past parallel twin cylinders[J]. Journal of Ship Mechanics, 2021, 25(5): 527‒534. (in Chinese)

韩忠华, 王绍楠, 韩莉, 等. 一种基于动模态分解的翼型流动转捩预测新方法[J]. 航空学报, 2017, 31(8): 120034.

Han Z H, Wang S N, Han L, et al. A novel method for automatic transition prediction of flows over airfoils based on dynamic mode decomposition[J]. Journal of Aeronautics Northwestern Polytechnical University, 2017, 31(8): 120034. (in Chinese)

Jonathan H T, Rowley C W, Dirk M L, et al. On dynamic mode decomposition: Theory and applications[J]. Journal of Computational Dynamics, 2014, 1(2): 391‒421.

Sampath R, Chakravarthy S R. Proper orthogonal and dynamic mode decompositions of time-resolved PIV of confined backward-facing step flow[J]. Experiments in Fluids, 2014, 55(9): 1792.

Seena A, Sung H J. Dynamic mode decomposition of turbulent cavity flows for self-sustained oscillations[J]. International Journal of Heat and Fluid Flow, 2011, 32(6): 1098‒1110.

Rowley C W, Mezic I, Bagheri S, et al. Spectral analysis of nonlinear flows[J]. Journal of Fluid Mechanics, 2009, 641: 115‒127.

丁杰, 胡晨星, 刘鹏寅, 等. 无叶扩压器中非定常流动的DMD模态分析[J]. 动力工程学报, 2017, 37(6): 447‒453.

Ding J, Hu C X, Liu P Y, et al. Dynamic mode decomposition of the unsteady flow in a vaneless diffuser[J]. Journal of Chinese Society of Power Engineering, 2017, 37(6): 447‒453. (in Chinese)

张琦. 湍流载荷激励下换热器直管束随机振动特性研究[D]. 哈尔滨: 哈尔滨工程大学, 2022.

Zhang Q. Study on random vibration characteristics of straight tube bundle in shell-and-tube heat exchangers[D]. Harbin: Harbin Engineering University, 2022. (in Chinese)

顾罡. 二维单圆柱、双圆柱绕流问题和三维垂荡板运动的数值模拟[D]. 上海: 上海交通大学, 2007.

Gu G. Numerical simulation of two-dimensional flow around single & double cylinders and three-dimensional oscillation of heave plates[D]. Shanghai: Shanghai Jiao Tong University, 2007. (in Chinese)

叶春明, 吴文权. 数值模拟圆柱绕流旋涡运动及尾流不稳定性分析[J]. 工程热物理学报, 1997, 18(2): 169‒172.

Ye C M, Wu W Q. Numerical simulation of vortex motion and wake instability analysis around a cylinder without restraint[J]. Journal of Engineering Thermophysics, 1997, 18(2): 169‒172. (in Chinese)

王帅, 陈正寿, 杜炳鑫, 等. 柔性管异形表面结构对其涡激振动影响的数值研究[J]. 船舶力学, 2024, 28(3): 450‒464.

Wang S, Chen Z S, Du B X, et al. Numerical study on the effect of flexible pipe special surface structures on vortex-induced vibration[J]. Journal of Ship Mechanics, 2024, 28(3): 450‒464. (in Chinese)

Kutz J N, Brunton S L, Brunton B W, et al. Dynamic mode decomposition: Data-driven modeling of complex systems[M]. Society for Industrial and Applied Mathematics, Washington State University, 2016.

2026年第30卷第3期

PDF下载

107

引用本文

BibTeX

文章信息

doi: 10.3969/j.issn.1007-7294.2026.03.005

接收时间：2025-06-14
首发时间：2026-07-07
出版时间：2026-03-15

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2025-06-14

基金

作者信息

^1.哈尔滨工程大学动力与能源工程学院，哈尔滨 150001

^2.中国船舶科学研究中心，江苏无锡 214082

通讯作者:

曹贻鹏（1980–），男，副教授，通讯作者，E-mail：yipengcao06@163.com

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/cblx/CN/10.3969/j.issn.1007-7294.2026.03.005

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

对比结果	本文结果	文献1^[19]	文献2^[20]
St	0.188	0.192	0.190

对比结果

本文结果

文献1^[19]

文献2^[20]

0.188

0.192

0.190

参数类型	管束数	管束直径	管束分布圆直径	节距
符号/单位	N_t	D_t /m	D_d /m	P /m
参数值	7	0.032	0.128	0.048

参数类型

管束数

管束直径

管束分布圆直径

节距

符号/单位

N_t

D_t /m

D_d /m

P /m

参数值

0.032

0.128

0.048

工况	网格方案		评判标准参数
Case 1	120	0.00084	179436	0.1971
Case 2	160	0.00063	220685	0.2383
Case 3	200	0.00050	271341	0.2532
Case 4	200	0.00050	335071	0.2536

工况

网格方案

评判标准参数

$ {N}_{{\mathrm{C}}} $

$ {\Delta }_{{\mathrm{C}}}/{N}_{{\mathrm{C}}} $

$ N $

Case 1

120

0.00084

179436

0.1971

Case 2

160

0.00063

220685

0.2383

Case 3

200

0.00050

271341

0.2532

Case 4

200

0.00050

335071

0.2536

模态参数	Mode 0	Mode 1	Mode 2	Mode 3	Mode 4	Mode 5
振幅	2656591	34533	4552	1179	678	461
增长率/（×10⁻⁵）	0.00	−1.82	−8.98	−22.70	−22.17	−24.32
频率/Hz	0.00	16.59	33.19	49.78	66.37	82.97

模态参数

Mode 0

Mode 1

Mode 2

Mode 3

Mode 4

Mode 5

振幅

2656591

34533

4552

1179

678

461

增长率/（×10⁻⁵）

0.00

−1.82

−8.98

−22.70

−22.17

−24.32

频率/Hz

0.00

16.59

33.19

49.78

66.37

82.97