刘瑜^*, 邓家钰^*, 王成恩^,*,2), 苏红星^†*上海交通大学机械与动力工程学院, 上海 200240
^†大连理工大学航空航天学院, 辽宁大连 116024

A MONOLITHIC METHOD FOR SIMULATING CONJUGATE HEAT TRANSFER VIA QUASI-IMPLICIT SCHEME OF CHARACTERISTIC-BASED SPLIT FINITE ELEMENT¹⁾

Liu Yu^*, Deng Jiayu^*, Wang Chengen^,*,2), Su Hongxin^{† *}School of Mechanical Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
^†School of Aeronautics and Astronautics, Dalian University of Technology, Dalian 116024, Liaoning, China

通讯作者: 2)王成恩, 教授, 主要研究方向: 航空发动机数值传热分析技术. E-mail:c.wang@sjtu.edu.cn

收稿日期:2020-08-24接受日期:2021-01-20网络出版日期:2021-04-08

基金资助:

1)装备预研共用技术项目.41423010103

Received:2020-08-24Accepted:2021-01-20Online:2021-04-08
作者简介 About authors


摘要
共轭传热现象在科学和工程领域中大量存在. 随着计算能力的发展, 对共轭传热现象进行准确有效的数值模拟, 成为科学研究和工程设计上的重要挑战.共轭传热数值模拟的方法可以分为两大类: 分区耦合和整体耦合.本文采用有限元法对共轭传热问题进行整体耦合模拟. 固体传热求解采用标准的伽辽金有限元方法.流动求解采用基于特征分裂的有限元方法. 该方法是一种重要的求解流动问题的有限元方法, 可以使用等阶有限元. 该方法的准隐格式与其他格式相比, 具有时间步长大的特点. 将稳定项中的时间步长与全局时间步长分开, 改进了准隐格式的稳定性. 基于改进的特征分裂有限元方法的准隐格式, 发展了一种层流共轭传热数值模拟的整体耦合方法. 采用这种方法可以将流体计算域和固体计算域作为一个整体划分有限元网格, 并且所有变量都可以采用相同的插值函数, 从而有利于程序的实现. 通过对典型问题的模拟, 验证了这种方法的准确性. 本工作还研究了固体区域时间步长对定常共轭传热问题数值模拟收敛性的影响.
关键词： 共轭传热;整体耦合;不可压缩流动;有限元;数值模拟

Abstract
Conjugate heat transfer is widely present in the fields of science and engineering. With the development of computing power, the accurate and effective numerical simulation of conjugate heat transfer has become a major challenge in scientific research and engineering design. The method of numerical simulation of conjugate heat transfer can be divided into two main categories: partitioned method and monolithic method. Each of these methods has its pros and cons. We have developed a monolithic method for simulating the conjugate heat transfer between solid and incompressible laminar flows with the finite element method. Heat conduction in solid is solved by the standard Galerkin finite element method. The flow solution adopts the characteristic-based split finite element method (CBS). This method is an important finite element method for solving flow problems, and equal-order finite elements can be used. Compared with semi-implicit and CBS-AC schemes, the quasi-implicit scheme of this method can adopt a larger time-step. The stability of the quasi-implicit scheme is improved by distinguishing the time step in the stabilization item from the global time step. Based on the quasi-implicit scheme of the improved CBS method, a monolithic method of conjugate heat transfer numerical simulation has been developed. In this way, the fluid part and solid part of the computational domain can be divided into finite element meshes as a whole, and the equal-order interpolation functions can be used for all variables, thus facilitating the realization of the program. The accuracy of this method is validated by simulating the benchmark problems. The effect of the time step for the solid domain on the convergence of steady conjugate heat transfer simulation has also been studied.
Keywords：conjugate heat transfer;monolithic;incompressible flow;finite element;numerical simulation


PDF (3578KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
刘瑜, 邓家钰, 王成恩, 苏红星. 基于特征分裂有限元准隐格式的共轭传热整体耦合数值模拟方法¹⁾. 力学学报[J], 2021, 53(4): 986-997 DOI:10.6052/0459-1879-20-299
Liu Yu, Deng Jiayu, Wang Chengen, Su Hongxin. A MONOLITHIC METHOD FOR SIMULATING CONJUGATE HEAT TRANSFER VIA QUASI-IMPLICIT SCHEME OF CHARACTERISTIC-BASED SPLIT FINITE ELEMENT¹⁾. Chinese Journal of Theoretical and Applied Mechanics[J], 2021, 53(4): 986-997 DOI:10.6052/0459-1879-20-299

引言

共轭传热现象在科学和工程领域中大量存在^[1], 比如核反应系统^[2], 燃气轮机涡轮叶片冷却^[3-4], 航天器热防护^[5-6], 化工处理^[7], 生物医学^[8-9], 电子设备散热^[10], MEMS设备中的微槽道对流换热^[11]等.

共轭传热问题除了实验研究以外, 可以通过理论分析或数值模拟的方法求解.理论分析方法对简单的模型(如平板, 圆管等)进行建模, 通过简化得到精确或近似的理论解^[12]. 理论方法对共轭传热问题的理解很重要, 但只适用于简单的问题.数值模拟将流动的数值解与固体传热的数值解耦合起来, 主要有两种耦合方法: 分区耦合(partitioned)^[13]和整体耦合(monolithic)^[14]方法. 分区耦合方法对流场和固体热分析采用各自相应的成熟的求解器, 然后通过迭代方法满足流热耦合的界面条件, 即在界面上温度和热流应是连续的.分区耦合方法的优点是可以采用现有的流动和固体热分析求解程序, 但是为了满足界面条件需要进行迭代求解, 在不同的求解器之间传递边界条件, 算法实现复杂, 并且还需要考虑分区耦合算法在流固耦合界面上的稳定性条件^[15-20]. 整体耦合方法将流动方程和固体传热方程统一离散, 然后求解单一的方程.整体耦合方法的优点是不需要处理流固耦合界面条件, 缺点是如果流动特征时间和固体传热特征时间差别很大, 非定常共轭传热的计算量将是巨大的^[21-22].

目前存在大量的数值方法模拟共轭传热问题,如有限差分法^[23], 有限体积法^[6,24-25], 浸入式边界法^[26], 无网格法^[27], 格子玻尔兹曼方法^[28-30], 边界元法^[31], 有限元法等.有限元法在计算固体传热问题上具有优势, 如果流场求解也采用有限元求解, 可以降低程序实现的复杂性. 采用有限元法求解共轭传热的文献相对较少. Misa和Sarkar^[32]采用整体耦合方法, 利用满足inf-sup条件的经典伽辽金(Galerkin)有限元方法对方腔自然对流的共轭传热进行了模拟. Malatip 等^[33]采用Rice和Schnipke^[34]提出的交错有限元方法和整体耦合方法求解定常共轭传热问题. Malatip等将Choi和Yoo^[35]提出的二阶时间精度分裂有限元方法推广到求解共轭传热问题^[36]和流-热-结构耦合问题^[37]. Verstraete和Scholl^[19]以及Scholl等^[20]采用压力稳定有限元法计算流动, 定常有限元法计算固体传热, 对分区耦合方法计算定常共轭传热的稳定性问题进行了研究. Long等^[38]将光滑有限元方法(smoothed finite element method)和改进的光滑粒子流体动力学方法耦合起来求解热-流-结构相互作用问题.

利用伽辽金有限元法求解不可压缩流动, 需要有限单元满足所谓的inf-sup条件, 即速度和压力的插值函数不能具有相同的阶数^[39]. 这增加了算法实现的难度, 特别是对于三维问题更是如此. 为了采用等阶元计算, Hughes等^[40]提出了用PSPG (Pressure-Stabilizing/Petrov-Galerkin)稳定项来规避inf-sup条件, 为了实现对流主导流动的稳定计算还需要添加SUPG(Streamline-Upwind/Petrov-Galerkin)稳定项^[41]. Zienkiewicz等^[42]在分数步法的基础上, 提出了基于特征分裂的有限元法(CBS). CBS方法有半隐格式(semi-implicit), 准隐格式(quasi-implicit)和基于人工压缩性的格式(CBS-AC)^[42]. 准隐格式的时间步长只取决于对流稳定性, 计算效率高于半隐格式和基于人工压缩性的格式^[43]. 在作者阅读的有限文献范围内, 还没有看到将CBS方法的准隐格式用于共轭传热数值模拟的例子. CBS方法既适用于不可压缩流动的求解, 也适用于亚声速、超声速流动的求解; 并且同样可以采用等阶插值函数; CBS方法是分数步方法, 每一个分裂步得到的方程是标量方程, 对存储量要求小, 方程的求解也相对容易.本文的主要目的是采用CBS方法的准隐格式, 发展一种整体耦合层流共轭传热数值模拟方法.

本文接下来给出共轭传热的控制方程, 并介绍流场计算采用的基于特征分裂的有限元方法, 随后给出整体耦合共轭传热数值模拟算法.为了验证本文所采用算法的准确性, 对不可压缩流动问题和共轭传热问题进行了模拟.

1 计算方法

1.1 控制方程

非守恒形式的不可压缩Navier-Stokes 方程为

(1)$\frac{\partial u_i}{\partial t} = -u_j \frac{\partial u_i}{\partial x_j} - \frac{\partial p}{\partial x_i} + \frac{1}{Re} \frac{\partial^2 u_i}{\partial x_j^2} $
(2)$\frac{\partial u_i}{\partial x_i} = 0 $
上式是无量纲化的方程, 其中$u_i,(i=1,2,3)$ 是流体速度分量, $p$是流体静压, $Re$为雷诺数.方程中变量下标满足爱因斯坦求和规则.

假设流体的密度是常量, 且其他流体性质也与温度无关.将能量方程表示为温度的形式, 即

(3)$\begin{eqnarray} \frac{\partial {T}}{\partial {t}} + {u}_j \frac{\partial {T}}{\partial {x_j}} = \frac{1}{Re Pr} \frac{\partial^2 {T}}{\partial x_j^2} + \frac{1}{Re Pr} {Q}_{f} \label{eq3} \end{eqnarray} $
其中, $T$为温度, $Q_{f}$为流体中的热源项, $Pr$为普朗特数.在这种情形下,能量方程与连续性方程和动量方程解耦.温度只是简单的输运变量, 在得到了速度场后, 计算温度方程就可以得到温度场的分布.

固体传热方程为

(4)$\begin{eqnarray} cr \frac{\partial {T}}{\partial {t}} = \kappa_{sf} \frac{1}{Re Pr} \frac{\partial^2 {T}}{\partial {x_j^2}} + \frac{1}{Re Pr} {Q}_{s} \label{eq:temperature-solid} \end{eqnarray} $
其中$\kappa_{sf} = \kappa_{s} / \kappa_{f}$, $cr = \rho_{s} c_{ps}/\rho_{f} c_{pf}$. $\rho_{f}$, $\rho_{s}$分别是流体和固体的密度, $\kappa_{f},\kappa_{s}$分别是流体和固体的导热系数, $c_{ps}$和$c_{pf}$分别表示固体和流体的比热. 推导以上无量纲方程采用的参考量见附录.

1.2 求解不可压缩流动的分数步法

方程(1)的半离散形式为

(5)$\begin{eqnarray} \frac{u_i^{n+1}-u_i^n}{\Delta t} = -u_j^{n+\theta_1} \frac{\partial u_i^{n+\theta_1}}{\partial x_j} - \frac{\partial p^{n+\theta_2}}{\partial x_i} + \frac{1}{Re} \frac{\partial^2 u_i^{n+\theta_3}}{\partial x_j^2} \label{eq:semi-implicit-ins} \end{eqnarray} $
其中$\Delta t$是时间步长, $\theta_1$, $\theta_2$和$\theta_3$分别是控制对流项、压力项和黏性项时间离散的参数, 其取值范围为$[0,1]$.为了避免求解非线性方程, 只考虑$\theta_1=0$的情形. 在推导CBS格式之前, 首先引入辅助变量$\Delta u^{*}$和$\Delta u^{**}$, 且令

(6)$\begin{eqnarray} u_i^{n+1} - u_i^n &=& \Delta u_i^* + \Delta u_i^{**} \end{eqnarray} $
采用基于Helmholtz-Hodge分解的2级分数步方法, 将方程分裂为

(7)$\frac{\Delta u_i^*}{\Delta t} = -u_j^{n+\theta_1} \frac{\partial u_i^{n+\theta_1}}{\partial x_j} + \frac{1}{Re} \frac{\partial^2 u_i^{n+\theta_3}}{\partial x_j^2} $
(8)$\frac{\Delta u_i^{**}}{\Delta t} = -\frac{\partial p^{n+\theta_2}}{\partial x_i} $
对方程(8)取散度

(9)$\begin{eqnarray} \frac{\partial \Delta u_i^{**}}{\partial x_i} = -\Delta t \frac{\partial p^{n+\theta_2}}{\partial x_i} \end{eqnarray} $
由于$p^{n+\theta_2} = p^n + \theta_2(p^{n+1} - p^n) = p^n + \theta_2 \Delta p$, 代入式(9)得到

(10)$\begin{eqnarray} \frac{\partial \Delta u_i^{**}}{\partial x_i} = -\Delta t \frac{\partial p^{n}}{\partial x_i} -\Delta t\theta_2 \frac{\partial^2\Delta p}{\partial x_i^2} \label{eq10} \end{eqnarray} $
对式(2)在时间层$n+\theta_4$上展开, 其中$\theta_4$是控制连续性方程时间离散的参数,

(11)$\begin{eqnarray} u_i^{n+\theta_4} = u_i^n + \theta_4 (u_i^{n+1} - u_i^n) \end{eqnarray}$
故

(12)$\begin{eqnarray} \left(\frac{\partial u_i}{\partial x_i}\right)^{n+\theta_4} = \left(\frac{\partial u_i}{\partial x_i}\right)^{n} + \theta_4 \frac{\partial \Delta u_i^*}{\partial x_i} + \theta_4 \frac{\partial \Delta u_i^{**}}{\partial x_i} \end{eqnarray} $
令式(12)等于0, 并代入式(10)得到

(13)$\begin{eqnarray} \Delta t \theta_4 \theta_2 \frac{\partial^2 \Delta p}{\partial x_i^2} = \frac{\partial u_i^n}{\partial x_i} + \theta_4 \frac{\partial \Delta u_i^*}{\partial x_i} -\Delta t \theta_4 \frac{\partial^2 p^n}{\partial x_i^2} \label{eq13} \end{eqnarray} $
方程的求解顺序为: 首先求解式(7)得到$\Delta u_i^*$, 然后将其代入式(13)得到$\Delta p$, 随后计算式(8)得到$\Delta u_i^{**}$. 在计算得到速度场后, 代入方程(3)即可计算温度场.

1.3 对流稳定项

对于对流占优问题, CBS格式需要添加对流稳定项^[43]. 本文在半离散方程(5)中添加显式的稳定项, 即

(14)$\begin{eqnarray} &&\frac{u_i^{n+1} - u_i^n}{\Delta t} = -u_j^{n+\theta_1}\frac{\partial u_i^{n+\theta_1}}{\partial x_j} - \frac{\partial p^{n+\theta_2}}{\partial x_i} + \frac{1}{Re} \frac{\partial^2 u_i^{n+\theta_3}}{\partial x_j^2} +\\&& \frac{\Delta t_{s}}{2}u_k \frac{\partial}{\partial x_k} \left(-u_j^n\frac{\partial u_i^n}{\partial x_j} - \frac{\partial p^n}{\partial x_i} + \frac{1}{Re}\frac{\partial^2 u_i^n}{\partial x_j^2} \right) \end{eqnarray} $
其中$\Delta t_{s}$为对流稳定项需要的时间步长, 可以认为是稳定性参数, 与SUPG中的稳定参数$\tau$^[41]的作用相似. 在分数步方法中, 稳定项也要进行分裂, 在预测步式(7)中添加的稳定项为

(15)$\begin{eqnarray} \frac{\Delta t_{s}}{2}u_k \frac{\partial}{\partial x_k} \left(-u_j^n\frac{\partial u_i^n}{\partial x_j} + \frac{1}{Re}\frac{\partial^2 u_i^n}{\partial x_j^2} \right) \end{eqnarray} $
即不考虑压力项. 在修正步式(8)中, 稳定项仅包含压力项, 即

(16)$\begin{eqnarray} \frac{\Delta t_{s}}{2}u_k \frac{\partial}{\partial x_k} \left(-\frac{\partial p^n}{\partial x_i} \right) \end{eqnarray} $
类似可以得到带稳定项的能量方程半离散形式, 即

(17)$\begin{eqnarray} &&\frac{T^{n+1}-T^n}{\Delta t} = -u_j^{n+\theta_1}\frac{\partial T^{n+\theta_1}}{\partial x_j} + \frac{1}{Re Pr} \frac{\partial^2 T^{n+\theta_3}}{\partial x_j^2} +\\ &&\frac{1}{Re Pr} Q_{f}^n + \frac{\Delta t_{s}}{2}u_k \frac{\partial}{\partial x_k} \Bigg(-u_j^{n}\frac{\partial T^{n}}{\partial x_j} + \frac{1}{Re Pr} \frac{\partial^2 T^{n}}{\partial x_j^2} +\frac{1}{Re Pr} Q_{f}^n \Bigg) \end{eqnarray} $
其中最后一项是对流稳定项.

1.4 基于特征分裂有限元方法的矩阵形式

对半离散方程采用伽辽金有限元方法进行空间离散, 根据参数的不同, 可以得到不同的格式^[43].

1.4.1 半隐格式

令$1/2 \leqslant \theta_2 \leqslant 1$, 对黏性项进行显式处理, 即$\theta_3 = 0$. 半隐格式的完全离散矩阵形式为:

步骤1

(18)$\begin{eqnarray} &&{M}_{{u,i}} \Delta {u}_i^* = - \Delta t {C} {u}_i^n - \Delta t {K} {u}_i^n + \Delta t {F}_i^n + \Delta t \Delta t_{s} {f}_{{s},i}^n \ \ (i=1,2,3) \end{eqnarray} $
其中${M}_{u,i}$, ${C}$, ${K}$分别表示速度质量矩阵, 离散对流算子和黏性算子. ${F}_i$表示黏性边界积分, ${f}_{{s},i}$是稳定项.

步骤2

(19)$\begin{eqnarray} \theta_4 \theta_2 {L} \Delta {p} = \frac{1}{\Delta t} \left({D}_k {u}_k^n + \theta_4 {D}_k \Delta {u}_k^* - \Delta t \theta_4 {L} {p}^n\right) \end{eqnarray} $
步骤3

(20)$\begin{eqnarray} {M}_{u,i} \Delta {u}_i = {M}_{u,i} \Delta {u}_i^{*} - \Delta t{G}_i({p}^n + \theta_2 \Delta {p}) + \Delta t \Delta t_{s} {f }_{p,i}^n \end{eqnarray} $
其中${L},{D},{G}$分别表示离散拉普拉斯算子, 离散散度算子和离散梯度算子, ${f}_{p,i}$是稳定项.半隐格式是条件稳定的, 稳定性由对流项和黏性项确定.

在计算得到了速度场后, 即可以将速度场代入到能量方程计算温度场.

步骤4

(21)$\begin{eqnarray} &&{M}_{T} \Delta {T}^{n+1} = - \Delta t {C} {T}^n - \Delta t {K}_{T} {T}^n + \Delta t {F}_{T}^n +\Delta t {Q}_{T}^n +\Delta t \Delta t_{s} {f}_{T}^n \end{eqnarray} $
其中, ${M}_{T}$, ${K}_{T}$分别是温度方程质量矩阵和温度扩散算子, ${F}_{T}$表示温度扩散项边界积分, ${Q}_{T}$为热源项, ${f}_{T}$为稳定项.

半隐格式的单元稳定时间步长为

(22)$\Delta t^e = \min\left( \Delta t_{con}^e,\Delta t_{vis}^e\right) $
(23)$\Delta t_{con}^e = \frac{h}{|{{u}}|}, \ \ \Delta t_{vis} = \frac{h^2 Re}{2} $
其中上标$e$表示单元, $h$是单元的网格尺度, $\Delta t_{con}$, $\Delta t_{ vis}$分别是对流时间步长和黏性时间步长.全局时间步长取所有单元稳定时间步长的最小值, 即

(24)$\begin{eqnarray} \Delta t = \min_{e} \Delta t^e \label{eq24} \end{eqnarray} $
黏性稳定性在一定条件下会对时间步长产生严重的限制.为了克服这一缺点, 提出了准隐格式.

1.4.2 准隐格式

与半隐格式不同, 准隐格式对黏性项进行隐式处理^[43], 即$1/2 \leqslant \theta_3 \leqslant 1$.

步骤1

(25)$\begin{eqnarray} &&\left({M}_{u,i} + \Delta t \theta_3 {K} \right) {u}_i^* = {M}_{u,i} {u}_i^n - \Delta t {C} {u}_i^n + (\theta_3 - 1)\Delta t {K} {u}_i^n+\\&& \Delta t {F}_i^n + \Delta t \Delta t_{s} {f }_{{s},i}^n \ \ (i=1,2,3) \end{eqnarray} $
黏性项的边界积分仍然作显式处理. 步骤2和步骤3与半隐格式相同.

步骤3还可以采用非投影形式, 即对于速度修正步

(26)$\begin{eqnarray} {M}{u}_i^{n+1} = {M} {u}_i^* - \Delta t{G}_i({p}^n + \theta_2 \Delta {p}) + \Delta t \Delta t_{s} {f }_{p,i}^n \end{eqnarray} $
步骤4

(27)$\begin{eqnarray} &&\left({M}_{T} + \Delta t \theta_3 {K}_{T} \right) {T}^{n+1} = {M}_{T} {T}^n - \Delta t {C} {T}^n+\\ && (\theta_3 - 1)\Delta t {K}_{T} {T}^n + \Delta t {F}_{T}^n + \Delta t {Q}_{T}^n + \Delta t \Delta t_{s} {f }_{T}^n \label{eq27} \end{eqnarray} $
准隐格式的时间步长只受对流项限制, 即

(28)$\begin{eqnarray} \Delta t = \min_{e} \Delta t_{con}^e \label{eq28} \end{eqnarray} $
以上矩阵和源项的具体形式见附录.

说明1: 半隐格式的时间步长受黏性稳定性条件影响;而基于人工可压缩性的格式(CBS-AC), 根据文献[43]报道, 准隐格式的计算效率比其高3 $\sim$ 100倍.因此本文计算采用准隐格式.

说明2: CBS方法允许等阶插值.本文选择等阶有限元, 即速度、压力、温度都采用相同的有限元基函数, 这大大方便了程序的实现.

说明3: 准隐格式求解的四步方程的矩阵是对称正定的, 故可以采用共轭梯度法求解.预处理矩阵本文采用雅克比方法计算.

说明4: 本文计算所采用的参数为$\theta_1=0$, $\theta_2=\theta_3=\theta_4=1$.

1.5 共轭传热耦合算法

在进行流固耦合模拟时, 需要划分流场网格和固体网格. 本文设计了一种流固耦合网格和边界条件表达方式. 出于简单考虑, 流固耦合界面是适配的, 即在流固耦合界面上流场网格点和固体网格点是一致的. 为了方便, 流体区域和固体区域作为一个统一的计算域划分单一的网格. 为了区分流体和固体网格, 需要用一个指标集来表示网格单元的材料属性, 考虑用一个材料属性文件来存储网格单元的属性, 即mateiral.txt.

对于温度场而言, 流固耦合边界属于内面, 因此在流固耦合边界上只需要指定流场变量边界条件.而在固体的其他边界上只需要指定温度边界条件, 不需要指定流场变量的边界条件.

如果采用前面提及的流固耦合网格及边界条件, 可以得到共轭传热耦合算法如Algorithm 1所示.


Algorithm 1
Algorithm 1共轭传热耦合算法

新窗口打开|下载CSV

固体传热方程离散后的形式为

(29)$\begin{eqnarray} &&\left({M}_{T} + \Delta t \theta_3 {K}_{T} \right) {T}^{n+1} = {M}_{T} {T}^n +\\&& \Delta t \left[(\theta_3 - 1) {K}_{T} {T}^n + {Q}_{T}^n + {F}_{T}^n \right] \label{eq29} \end{eqnarray} $
温度方程在流体区域按照式(27)离散, 在固体区域按照式(29)离散, 最后将温度方程装配为统一的方程组.

该算法是一种整体耦合算法, 与分区算法相比, 不需要通过迭代就可以隐式满足流热耦合界面条件.

1.6 时间步长

1.6.1 流体时间步长

Bevan等^[43]研究了CBS方法中的半隐格式、准隐格式和基于人工压缩性的格式的计算效率, 发现准隐格式由于只需要考虑对流时间步长, 相比于其他两种格式具有较高的计算效率.文献[43]中没有对时间推进步长和稳定项中的时间步长进行区分, 在稳定项中采用的也是时间推进步长, 即$\Delta t_{s} = \Delta t$, 并且建议采用半隐格式和准隐格式时, 时间步长应该选择全局时间步长.

作者认为稳定项中的时间步长$\Delta t_{s}$可以视为稳定性参数, 而稳定性参数应具有当地性质. 当计算采用全局时间步长时, 如果整个计算区域的稳定性时间步长都取同样的值, 有可能不能较好地抑制数值不稳定性. 因此本文采用准隐方法进行计算, 对物理推进时间步长和稳定项中的时间步长做了区分. 物理推进时间步长$\Delta t$按式(28)计算,取全场的最小时间步长作为全局时间步长. 稳定项中的时间步长为局部时间步长$\Delta t_{s}$,按照式(22)计算.

1.6.2 固体时间步长

在固体传热方程(4)中, 固体密度与比热的乘积与流体密度与比热的乘积之比, $cr$, 在一般情况下有$cr \gg 1$. 这导致传热方程时间尺度远大于流动方程的时间尺度. He等^[44] 估计了典型共轭传热问题的固体和流体的稳定时间步长之比为$10^3$ $\sim$ $10^5$. 当对共轭传热问题的流体温度方程和固体温度方程统一求解时, 如果采用流体部分的稳定时间步长进行计算, 则固体部分的温度场需要迭代很多步才能收敛. 对于定常问题, 为了提高收敛速度, 在固体部分需要采用比流体时间步长更大的时间步长

(30)$\begin{eqnarray} \Delta t_{sol} = n \cdot \Delta t_{flu} \end{eqnarray} $
其中$\Delta t_{flu}$和$\Delta t_{sol}$分别表示流动计算和固体传热计算采用的时间步长, $n$为放大因子.

式(30)中如果选择$n = cr$, 则相当于在方程(4)中令$cr = 1$, $\Delta t=\Delta t_{flu}=\Delta t_{sol}$. 这是文献[23, 45]采用的做法.本文定义名义上的$cr$值为$cr'$, 即

(31)$cr' = cr/n $
计算时用$cr'$代替$cr$. 名义上固体时间步长与流体时间步长相同, 但是通过调整$cr'$实际上改变了固体时间步长. $cr'$可以足够小, 在满足稳定性的前提下提高收敛速度.这将在下面的算例中得到验证.

2 算法验证

2.1 稳定项时间步长选取方法对计算结果的影响

首先对稳定项中的时间步长的选择方式对计算结果的影响进行比较. 采用3种方式计算稳定项中的时间步长, 分别是:与全局时间步长相同, 采用局部时间步长和取零值(即不添加稳定项). 计算采用的例子是方腔顶盖驱动流.方腔的顶盖按指定的速度运动, 方腔的其他边界固定.在顶盖的作用下, 方腔中流体将发生运动. 计算域$\varOmega$为$(0,1)\times(0,1)$, 边界条件都为Dirichlet条件. 顶盖的运动速度, 按照文献[46], 给定如下

(32)$\begin{eqnarray} u_0 = \frac{1}{2}\left[-\tanh\left(\frac{|x-0.5|-0.495}{0.001} \right) + 1 \right] \end{eqnarray} $
当网格足够稠密时, 稳定项的作用不明显, 为此在这里特意选择比较稀疏的网格进行计算. 网格包含800个线性三角形单元和441个节点. 网格如图1(a)所示. 图1(b) $\sim\!$图1(d)分别给出了由3种方法计算的$Re=5000$时$x$方向上的速度分量$u$的云图. 图2给出了在$y=0.8$上$u$的分布. 可见如果不加入稳定项, 流场变量将出现很大的伪振荡;采用文献[43]中的方法, 由于稳定项中的时间步长全场都是相同的, 不能完全抑制伪振荡的出现; 而采用本文的方法流场变量分布光滑, 抑制了伪振荡的出现.

图1

新窗口打开|下载原图ZIP|生成PPT
图1顶盖驱动流的计算网格和速度云图

Fig.1Computational mesh and contours of $x$-component of velocity

图2

新窗口打开|下载原图ZIP|生成PPT
图2由不同稳定项计算方法得到的速度的$x$分量在$y=0.8$ 上的分布

Fig.2Velocity profile at $y=0.8$ predicted by different stabilization schemes

2.2 层流计算验证

2.2.1 方腔顶盖驱动流

方腔顶盖驱动流是用来检验层流计算准确性的典型算例. 计算和边界条件与上节相同. 采用一阶三角形单元计算, 单元数为3200, 节点数为1681, 在边界附近进行了加密. 分别对$Re=400$, 1000, 5000, 10 000时的情形进行了计算. 图3比较了$x=0.5$上速度分量$v$和$y=0.5$上速度分量$u$的分布与Ghia等^[47]的计算结果, 可见符合很好.

图3

新窗口打开|下载原图ZIP|生成PPT
图3不同雷诺数条件下的速度型面

Fig.3Velocity profiles in the $x$ and $y$ directions



2.2.2 层流后台阶流动

对于这个问题有实验数据可供比较, 因而为很多作者所引用. 本文选择雷诺数为229的情形进行计算. 问题的计算域和边界条件如图4所示. 进口速度$u$的分布的表达式为

(33)$\begin{eqnarray} &&u(y) = 0.662 4y^6-7.554 7y^5+33.9y^4-75.283y^3+ 83.368y^2-37.793y+2.695 9 \end{eqnarray} $
在壁面处按无滑移边界条件处理, 在出口速度指定Neumann条件, 压力指定为0, 在其余边界上压力指定为Neumann条件. 计算网格包含15 200个三角形单元, 7841个节点.

图4

新窗口打开|下载原图ZIP|生成PPT
图4层流后台阶流动计算域和边界条件

Fig.4Domain and boundary conditions of flow over a backward facing step



图5分别给出了压力云图和速度云图. 图(6)给出了不同截面上速度分布与实验数据^[48]的比较, 可见计算结果与实验符合很好.

图5

新窗口打开|下载原图ZIP|生成PPT
图5层流后台阶流动的压力和速度云图

Fig.5Pressure and velocity contours of flow over a backward facing step

图6

新窗口打开|下载原图ZIP|生成PPT
图6层流后台阶流动速度型面与实验的比较

Fig.6Comparison of present predication of velocity profiles at different locations with experimental data

2.3 共轭传热数值模拟

本文对若干共轭传热问题进行模拟, 以检验算法和程序的准确性. 其中, 反向流动热交换器在模拟中考虑了$cr'$的取值对收敛速度的影响, 其余算例都设置$cr'=1$.

2.3.1 共轭传热库埃特流

共轭传热库埃特流^[33]的示意图为图7.上壁以固定速度运动, 温度为$T_2 = 0$; 下壁固定具有厚度, 需要考虑热传导, 下壁底部温度为$T_1 = 1$; 流固耦合界面上温度的精确解为

(34)$\begin{eqnarray} T_{w} = \frac{\kappa_{f} h_1 T_2 + \kappa_{s} h_2 T_1}{\kappa_{f} h_1 + \kappa_{s} h_2} \end{eqnarray} $
其中$h_1$, $h_2$分别为固体和流体区域的厚度. 流体域中的温度和固体域中的温度在$y$方向上为线性分布.

图7

新窗口打开|下载原图ZIP|生成PPT
图7平行平板间的共轭传热库埃特流

Fig.7Conjugate Couette flow problem in parallel plate channel



图8给出了取不同的$\kappa_{sf}$值时, 由一阶三角形单元计算的温度沿$y$方向上的分布曲线, 网格单元数为1200. 可见计算值与理论值十分吻合.

图8

新窗口打开|下载原图ZIP|生成PPT
图8不同条件下共轭传热库埃特流温度分布曲线与理论解的比较

Fig.8Comparison of benchmark solutions for Couette flow condition



2.3.2 二维反向流动热交换器

二维反向流动热交换器的示意图如图9所示. 在热交换器中分布两个平行流道, 流道被一金属平板分隔开. 两个流道的厚度分别为$\delta_1$, $\delta_3$, 隔板的厚度为$\delta_2$. 流道的外壁是绝热的, 流道入口处的速度和温度被认为是均匀的. 计算采用的参数如下: 几何参数$\delta_1=\delta_2=\delta_3=0.1$ m, $L=1$ m; 流动参数$U_1=0.2$ m/s, $T_1=800$ K, $U_3=0.1$ m/s, $T_3=300$ K, $\rho_{f} = 1000$ kg/m$^3$, $\kappa_{f} = 10$ W/(m$\cdot$K), $c_{pf} = 25$ J/(kg$\cdot$K), $\mu = 0.15$ kg(m$\cdot$s); 金属板的物性$\rho_{s} = 8000$ kg/m$^3$, $\kappa_{s} = 50$ W/(m$\cdot$K), $c_{ps} = 500$ J/(kg$\cdot$K).

图9

新窗口打开|下载原图ZIP|生成PPT
图9二维反向流动热交换器

Fig.9A conjugate counter flow heat exchanger



计算网格采用线性三角形单元, 单元数为3360, 节点数为1763.

在本例中$cr$的准确值为160. 为了考虑$cr'$的取值对收敛速率的影响, 分别对$cr'$采用准确值和$cr'=1$, 0.001的情形进行了模拟, 图10给出了温度残值的收敛曲线. 可见$cr'$如果采用准确的值, 收敛将非常缓慢. 如果令$cr'=1$, 收敛速度将大幅增大, 继续减小$cr'$的值,收敛速度会继续增大. $cr'$的取值不影响最终的收敛解.

图10

新窗口打开|下载原图ZIP|生成PPT
图10采用不同的$cr'$ 的收敛速率比较

Fig.10Comparison of convergence rate for different $cr'$ values



计算的速度场和温度场分别如图11(a)和图11(b)所示. 本文还计算了不同的固体和流体热传导系数之比条件下的共轭传热. 在$x=L/2$上温度的分布曲线如图12所示, 并与Malatip等^[37]计算的结果进行比较, 可见与文献给出的结果符合良好.

图11

新窗口打开|下载原图ZIP|生成PPT
图11二维反向流动热交换器 $x$ 方向上的速度云图和温度云图

Fig.11Predicted $x$-component of velocity and temperature contours for a counter flow heat exchanger

图12

新窗口打开|下载原图ZIP|生成PPT
图12分别取$\kappa_{sf}=0.1$, 1, 5, 10, 二维反向流动热交换器的温度在 $x = L/2$ 上的分布与文献[37] 的比较

Fig.12Temperature distributions at $x = L/2$ for $\kappa_{sf}=0.1$, 1, 5, 10 of a counter flow heat exchanger compared with results from Ref.[37]



2.3.3 加热固体的强迫对流冷却

第3个测试例子是流动流经两块平行平板间含有3个矩形加热固体的强迫对流冷却问题. 该问题首次由Davalath和Bayazitoglu^[23] 提出用来模拟集成电路元件的冷却. 该问题的计算域和边界条件如图13所示, 矩形块尺寸相同, 且等间距分布. 流体以完全发展的速度型面从左侧进入, 从右侧流出. 3个固体块中的生成的热可以假设为具有均匀分布的热源项, 本文中假设热源项$Q_s=8$. 计算时取$Pr=0.71$, $\kappa_{sf}=10$, 分别计算了雷诺数为$Re= 100$, 500, 1000的情形. 计算网格如图14所示, 网格中含有个10 900个线性三角形单元, 5708个网格节点. 图15给出了雷诺数为100时的压力、速度和温度云图. 图16给出了固体-流体界面上温度的分布, 并与文献[37]的结果进行了比较, 可见本文计算结果与文献符合很好.

图13

新窗口打开|下载原图ZIP|生成PPT
图13加热固体的强迫对流冷却问题的计算域和边界条件

Fig.13Domain and boundary conditions of forced convection cooling across rectangular blocks

图14

新窗口打开|下载原图ZIP|生成PPT
图14加热固体的强迫对流冷却问题的计算网格

Fig.14Computational mesh of forced convection cooling across rectangular blocks

图15

新窗口打开|下载原图ZIP|生成PPT
图15雷诺数为 100 时, 加热固体的强迫对流冷却的压力, 速度幅度和温度云图

Fig.15Pressure, velocity magnitude and temperature contours for $Re=100$

图16

新窗口打开|下载原图ZIP|生成PPT
图16雷诺数为 100, 500, 1000 时固体-流体界面上温度的分布与文献[37]的比较

Fig.16Comparison of wall temperature distributions along solid-fluid interface of rectangular blocks for $Re=100$, 500, 1000 with Ref.[37]

3 结论

本文通过将物理推进时间步长与稳定项中的时间步长进行区别, 改进了基于特征分裂的有限元法的准隐格式, 使其具有了更好的稳定性. 在此基础上, 利用改进了的准隐格式和整体耦合方法实现了不可压缩流体和固体间共轭传热的数值模拟并通过与理论、 实验和文献中的数值解进行比较, 验证了本文提出方法的正确性. 还研究了固体区域时间步长对定常共轭传热问题数值模拟收敛性的影响, 发现选择合适的固体时间步长可以显著提高收敛速率.

本文的一个不足之处是没有对稳定时间步长对基于特征分裂的有限元法准隐格式的稳定性的影响进行理论分析, 在未来的工作中将对此进行研究; 同时还可以研究物理推进时间步长采用局部时间步长对定常共轭传热问题计算收敛性的影响. 未来还将采用本文提出的方法对自然对流和混合对流问题进行研究, 进一步探索该方法的适用范围.

附录

控制方程的无量纲化

有量纲形式的控制方程为:

连续性方程

(A1)$\frac{\partial u_{i}}{\partial x_{i}}=0$
动量方程

(A2)$\frac{\partial \rho_{\mathrm{f}} u_{i}}{\partial t}+\rho_{\mathrm{f}} u_{j} \frac{\partial u_{i}}{\partial x_{j}}=-\frac{\partial p}{\partial x_{i}}+\mu \frac{\partial^{2} u_{i}}{\partial x_{j}^{2}}$
能量方程

(A3)$\frac{\partial \rho_{\mathrm{f}} c_{p, f} T}{\partial t}+\rho_{\mathrm{f}} c_{p, f} u_{j} \frac{\partial T}{\partial x_{j}}=\kappa_{\mathrm{f}} \frac{\partial^{2} T}{\partial x_{j}^{2}}+Q_{\mathrm{f}}$
固体中的传热方程

(A4)$\frac{\partial \rho_{\mathrm{s}} c_{p, s} T}{\partial t}=\kappa_{\mathrm{s}} \frac{\partial^{2} T}{\partial x_{j}^{2}}+Q_{\mathrm{s}}$
选取参考量为$\rho_{ref}=\rho_{f}$, $V_{ref}$, $L_{ref}$, $p_{ref} = \rho_{ref} V_{ref}^2$, $\mu_{ref}=\mu$, $c_{p,ref}=c_{pf}$, $\kappa_{ref}=\kappa_{f}$, $t_{ref}=L_{ref}/V_{ref}$, $T_{ref} = {\mu_{ref} V_{ref}^2}/{\kappa_{ref}}$, $Q_{ref} = {\mu_{ref} V_{ref}}/{L_{ref}^2}$, 对方程进行无量纲化,即得到本文采用的无量纲化方程.

CBS 方程中的矩阵和源项

(A5)$\left.\begin{array}{l} \boldsymbol{M}_{u, i}=\int_{\Omega} \boldsymbol{N}^{\mathrm{T}} \boldsymbol{N} \mathrm{d} \Omega, \boldsymbol{K}=\int_{\Omega} \frac{1}{R e} \frac{\partial \boldsymbol{N}^{\mathrm{T}}}{\partial x_{j}} \frac{\partial \boldsymbol{N}}{\partial x_{j}} \mathrm{~d} \Omega \\ \boldsymbol{C}=\int_{\Omega} \boldsymbol{N}^{\mathrm{T}} u_{j} \frac{\partial \boldsymbol{N}}{\partial x_{j}} \mathrm{~d} \Omega, \boldsymbol{F}_{i}^{n}=\int_{\Gamma} \boldsymbol{N}^{\mathrm{T}} \boldsymbol{t}_{i} \mathrm{~d} \Gamma \end{array}\right\}$
其中

(A6)$\begin{eqnarray} {t} ={\sigma} \cdot {n} \end{eqnarray}$
${\sigma}$ 是应力张量, ${n}$ 是边界外法向.

(A7)$\boldsymbol{M}_{\mathrm{T}}=\int_{\Omega} c r^{\prime} \boldsymbol{N}^{\mathrm{T}} \boldsymbol{N} \mathrm{d} \Omega, \boldsymbol{K}_{\mathrm{T}}=\int_{\Omega} \frac{\kappa_{\mathrm{sf}}}{\operatorname{RePr}} \frac{\partial \boldsymbol{N}^{\mathrm{T}}}{\partial x_{j}} \frac{\partial \boldsymbol{N}}{\partial x_{j}} \mathrm{~d} \Omega$
(A8)$\boldsymbol{Q}_{\mathrm{T}}=\int_{\Omega} \frac{\kappa_{\mathrm{sf}}}{\operatorname{RePr}} N^{\mathrm{T}} Q \mathrm{~d} \Omega$
在流体域中$cr'=1$, $\kappa_{sf}=1$, $Q=Q_{f}$, 在固体域中$Q=Q_{s}$.

(A9)$\boldsymbol{L}=-\int_{\Omega} \frac{\partial \boldsymbol{N}^{\mathrm{T}}}{\partial x_{j}} \frac{\partial \boldsymbol{N}}{\partial x_{j}} \mathrm{~d} \Omega$
(A10)$\boldsymbol{D}_{i}=\int_{\Omega} \boldsymbol{N}^{\mathrm{T}} \frac{\partial \boldsymbol{N}}{\partial x_{i}} \mathrm{~d} \Omega$
(A11)$\boldsymbol{G}_{i}=\int_{\Omega} \boldsymbol{N}^{\mathrm{T}} \frac{\partial \boldsymbol{N}}{\partial x_{i}} \mathrm{~d} \Omega$
稳定项

(A12)$f_{\mathrm{s}, i}^{n}=\frac{1}{2} \int_{\Omega} u_{k} \frac{\partial \boldsymbol{N}^{\mathrm{T}}}{\partial x_{k}}\left(-u_{j} \frac{\partial u_{i}}{\partial x_{j}}+\frac{1}{R e} \frac{\partial^{2} u_{i}}{\partial x_{j}^{2}}+\boldsymbol{S}_{i}\right) \mathrm{d} \Omega$
(A13)$f_{p, i}^{n}=\frac{1}{2} \int_{\Omega} u_{k} \frac{\partial N^{\mathrm{T}}}{\partial x_{k}}\left(-\frac{\partial p}{\partial x_{i}}\right) \mathrm{d} \Omega$
(A14)$f_{T, s}^{n}=\frac{1}{2} \int_{\Omega} u_{k} \frac{\partial N^{\mathrm{T}}}{\partial x_{k}}\left(-u_{j} \frac{\partial T}{\partial x_{j}}+\frac{1}{\operatorname{RePr}} \frac{\partial^{2} T}{\partial x_{j}^{2}}+Q_{\mathrm{f}}\right) \mathrm{d} \Omega$
其中${S}$ 是体积力项, ${N}$ 是单元形状函数.

参考文献 View Option 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

[1] Afzal A, Samee ADM, Razak RKA, et al. Steady and transient state analyses on conjugate laminar forced convection heat transfer
Archives of Computational Methods in Engineering, 2020,27:135-170
[本文引用: 1]

[2] 冯琳娜, 李权, 黄永忠 等. 基于流固共轭传热的两相流动数值模拟及燃料子通道CHF 预测研究
原子能科学技术, 2020,54(6):1092-1098
[本文引用: 1]

(Feng Linna, Li Quan, Huang Yongzhong, et al. Subcooled boiling simulation and prediction of critical heat flux in fuel subchannel based on conjugate heat transfer
Atomic Energy Science and Technology, 2020,54(6):1092-1098 (in Chinese))
[本文引用: 1]

[3] 周源远, 范小军, 李亮 等. 燃气轮机叶片前缘旋流冷却的热流固耦合数值研究
西安交通大学学报, 2020,54(1):135-142
[本文引用: 1]

(Zhou Yuanyuan, Fan Xiaojun, Li Liang, et al. Numerical research on the swirl cooling using thermal-fluid-structure coupled method for turbine blade leading edge
Journal of Xi'an Jiaotong University, 2019,45(9):1700-1712 (in Chinese))
[本文引用: 1]

[4] Wang X, Xu HZ, Wang JH, et al. Multi-objective optimization of discrete film hole arrangement on a high pressure turbine end-wall with conjugate heat transfer simulations
International Journal of Heat and Fluid Flow, 2019,78:108428
[本文引用: 1]

[5] 桂业伟, 刘磊, 代光月 等. 高超声速飞行器流-热-固耦合研究现状与软件开发
航空学报, 2017,38(7):87-105.
[本文引用: 1]

(Gui Yewei, Liu Lei, Dai Guangyue, et al. Research status of hypersonic vehicle fluid-thermal-solid coupling and software development
Acta Aeronautica, 2017,38(7):87-105(in Chinese))
[本文引用: 1]

[6] H?tte F, Haupt C. Transient 3d conjugate heat transfer simulation of a rectangular GOX-GCH4 rocket combustion chamber and validation
Aerospace Science and Technology, 2020,105:106043
[本文引用: 2]

[7] Das S, Panda A, Deen NG, et al. A sharp-interface immersed boundary method to simulate convective and conjugate heat transfer through highly complex periodic porous structures
Chemical Engineering Science, 2018,191:1-18
[本文引用: 1]

[8] Saxena A, Saha V, Ng EYK. Skin temperature maps as a measure of carotid artery stenosis
Computers in Biology and Medicine, 2020,116:103548
DOIURLPMID [本文引用: 1]
In this study, the effect of carotid artery stenosis on the neck skin temperature maps was investigated. With the presence of stenosis, alterations in the carotid artery hemodynamics bring about changes in the heat transfer to the surrounding tissue. This is expected to be captured in the resulting temperature map over the external neck skin surface; possibly it correlates to the presence of stenosis. A total of twenty carotid artery samples, from ten patients with both sides normal (0% stenosis), stenosis (>50%) on one side, and stenosis (>50%) on both sides, were studied. Duplex Ultrasound and infrared (IR) thermography examinations were performed. A computational study, on an ideal 3-dimensional (3D) carotid artery and jugular vein model encapsulated with a solid neck tissue phantom resembling the human neck, was carried out. Incorporating the patient-specific geometrical (depth of artery and stenosis) and flow (peak systolic and end diastolic inlet velocity) boundary conditions, conjugate bio-heat transfer was studied using a finite volume numerical scheme. Simulation results and in-vivo thermal maps show that the average temperature on the external neck skin surface is significantly higher for normal patients (32.82 +/- 0.53 degrees C versus 32.00 +/- 0.37 degrees C, p < 0.001). Furthermore, the thermal region of interests (TROIs) were extracted from the in-vivo thermal images, which both qualitatively and quantitatively distinguish the normal and diseased cases. This study suggests the potential of thermal feature-based screening of patients with carotid artery stenosis.

[9] Alsabery AI, Naganthran K, Azizul FM, et al. Numerical study of conjugate natural convection heat transfer of a blood filled horizontal concentric annulus
International Communications in Heat and Mass Transfer, 2020,114:104568
[本文引用: 1]

[10] Dbouk T. A new technology for CPU chip cooling by concentrated suspension flow of non-colloidal particles
Applied Thermal Engineering, 2019,146:664-673
[本文引用: 1]

[11] Prajapati YK. Influence of fin height on heat transfer and fluid flow characteristics of rectangular microchannel heat sink
International Journal of Heat and Mass Transfer, 2019,137:1041-1052
[本文引用: 1]

[12] Dorfman AS. Conjugate Problems in Convective Heat Transfer. Boca Raton: CRC Press, 2010
[本文引用: 1]

[13] Meng F, Banks JW, Henshaw WD, et al. A stable and accurate partitioned algorithm for conjugate heat transfer
Journal of Computational Physics, 2017,344:51-85
[本文引用: 1]

[14] Pan XM, Lee CH, Choi JI. Efficient monolithic projection method for time-dependent conjugate heat transfer problems
Journal of Computational Physics, 2018,369:191-208
[本文引用: 1]

[15] Errera MP, Duchaine F. Comparative study of coupling coefficients in dirichlet —— robin procedure for fluid —— structure aerothermal simulations
Journal of Computational Physics, 2016,312:218-234
[本文引用: 1]

[16] Errera MP, Moretti R, Salem R, et al. A single stable scheme for steady conjugate heat transfer problems
Journal of Computational Physics, 2019,394:491-502


[17] Giles MB. Stability analysis of numerical interface conditions in fluidstructure thermal analysis
International Journal for Numerical Methods in Fluids, 1997,25(4):421-436


[18] Henshaw WD, Chand KK. A composite grid solver for conjugate heat transfer in fluid-structure systems
Journal of Computational Physics, 2009,228(10):3708-3741


[19] Verstraete T, Scholl S. Stability analysis of partitioned methods for predicting conjugate heat transfer
International Journal of Heat and Mass Transfer, 2016,101:852-869
[本文引用: 1]

[20] Scholl S, Janssens B, Verstraete T. Stability of static conjugate heat transfer coupling approaches using robin interface conditions
Computers & Fluids, 2018,172:209-225
[本文引用: 2]

[21] He L, Oldfield MLG. Unsteady conjugate heat transfer modeling
Journal of Turbomachinery, 2011,133(3):031022
[本文引用: 1]

[22] He L. Fourier spectral modelling for multi-scale aero-thermal analysis
International Journal of Computational Fluid Dynamics, 2013,27(2):118-129
[本文引用: 1]

[23] Davalath J, Bayazitoglu Y. Forced convection cooling across rectangular blocks
Journal of Heat Transfer, 1987,109(2):321-328
[本文引用: 3]

[24] Konar D, Sultan MA, Roy S. Numerical analysis of 2-d laminar natural convection heat transfer from solid horizontal cylinders with longitudinal fins
International Journal of Thermal Sciences, 2020,154:106391
[本文引用: 1]

[25] Costa R, Nobrega JM, Clain S, et al. Very high-order accurate polygonal mesh finite volume scheme for conjugate heat transfer problems with curved interfaces and imperfect contacts
Computer Methods in Applied Mechanics and Engineering, 2019,357:112560
[本文引用: 1]

[26] Lee S, Hwang W. Development of an efficient immersed-boundary method with sub grid-scale models for conjugate heat transfer analysis using large eddy simulation
International Journal of Heat and Mass Transfer, 2019,134:198-208
[本文引用: 1]

[27] Wijayanta AT, Pranowo. A localized meshless approach using radial basis functions for conjugate heat transfer problems in a heat exchanger
International Journal of Refrigeration, 2020,110:38-46
[本文引用: 1]

[28] 李桥忠, 陈木凤, 李游 等. 浸没边界-简化热格子Boltzmann 方法研究及其应用
力学学报, 2019,51(2):392-404
[本文引用: 1]

(Li Qiaozhong, Chen Mufeng, Li You, et al. Non-orthogonal multiple-relaxation-time lattice Boltzmann method for numerical simulation of thermal coupling with porous square cavity flow containing internal heat source
Chinese Journal of Theoretical and Applied Mechanics, 2019,51(2):392-404 (in Chinese))
[本文引用: 1]

[29] Mohebbi R, Lakzayi H, Rasam H. Numerical simulation of conjugate heat transfer in a square cavity consisting the conducting partitions by utilizing lattice Boltzmann method
Physica A: Statistical Mechanics and its Applications, 2020,546:123050


[30] Spinelli GG, Celik B. Leveling out interface temperature for conjugate heat transfer problems
Computers & Fluids, 2020,210:104652
[本文引用: 1]

[31] Blobner J, Hriber?ek M, Kuhn G. Dual reciprocity BEM-BDIM technique for conjugate heat transfer computations
Computer Methods in Applied Mechanics and Engineering, 2000,190(8):1105-1116
[本文引用: 1]

[32] Misra D, Sarkar A. Finite element analysis of conjugate natural convection in a square enclosure with a conducting vertical wall
Computer Methods in Applied Mechanics and Engineering, 1997,141(3):205-219
[本文引用: 1]

[33] Malatip A, Wansophark N, Dechaumphai P. Combined streamline upwind petrov galerkin method and segregated finite element algorithm for conjugate heat transfer problems
Journal of Mechanical Science and Technology, 2006,20(10):1741-1752
[本文引用: 2]

[34] Rice JG, Schnipke RJ. An equal-order velocity-pressure formulation that does not exhibit spurious pressure modes
Computer Methods in Applied Mechanics and Engineering, 1986,58(2):135-149
[本文引用: 1]

[35] Choi HG, Yoo JY. Streamline upwind scheme for the segregated formulation of the Navier-Stokes equation
Numerical Heat Transfer, Part B: Fundamentals, 1994,25(2):145-161
[本文引用: 1]

[36] Malatip A, Wansophark N, Dechaumphai P. A second-order time accurate finite element method for analysis of conjugate heat transfer between solid and unsteady viscous flow
Journal of Mechanical Science and Technology, 2009,23(3):775-789
[本文引用: 1]

[37] Malatip A, Wansophark N, Dechaumphai P. Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems
Applied Mathematics and Mechanics, 2012,33(1):99-116
[本文引用: 7]

[38] Long T, Yang PY, Liu MB. A novel coupling approach of smoothed finite element method with SPH for thermal fluid structure interaction problems
International Journal of Mechanical Sciences, 2020,174:105558
[本文引用: 1]

[39] Reddy JN. The Finite Element Method in Heat Transfer and Fluid Dynamics, 3e. Boca Raton: CRC Press, 2010
[本文引用: 1]

[40] Hughes TJ, Franca LP, Balestra M. A new finite element formulation for computational fluid dynamics: V.circumventing the Babu? ka-Brezzi condition: a stable petrov-galerkin formulation of the stokes problem accommodating equal-order interpolations
Computer Methods in Applied Mechanics and Engineering, 1986,59(1):85-99
[本文引用: 1]

[41] Tezduyar T. Stabilized finite element formulations for incompressible flow computations
Advances in Applied Mechanics, 1991,28:1-44
[本文引用: 2]

[42] Zienkiewicz OC, Taylor RL, Zhu JZ. The Finite Element Method: Its Basis and Fundamentals,7e. London: Butterworth-Heinemann Press, 2013
[本文引用: 2]

[43] Bevan RL, Boileau E, van Loon R, et al. A comparative study of fractional step method in its quasi-implicit, semiimplicit and fully-explicit forms for incompressible flows
International Journal of Numerical Methods for Heat & Fluid Flow, 2016,26(3/4):595-623
[本文引用: 8]

[44] He L. Closely coupled fluid-solid interface method with moving-average for LES based conjugate heat transfer solution
International Journal of Heat and Fluid Flow, 2019,79:108440
[本文引用: 1]

[45] Chen X, Han P. A note on the solution of conjugate heat transfer problems using SIMPLE-like algorithms
International Journal of Heat and Fluid Flow, 2000,21(4):463-467
[本文引用: 1]

[46] Li LF. A split-step finite-element method for incompressible Navier-Stokes equations with high-order accuracy up-to the boundary
Journal of Computational Physics, 2020,408:109274
[本文引用: 1]

[47] Ghia U, Ghia K, Shin C. High-re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method
Journal of Computational Physics, 1982,48(3):387-411
[本文引用: 1]

[48] Denhamand MK, Patrick MA. Laminar flow over a down-stream facing step in a two-dimensional flow channel
Transactions of the Institution of Chemical Engineers, 1974,52:361-367
[本文引用: 1]