In the vicinity of the liquid--vapor critical point, supercritical fluids behave strongly compressibly and, in parallel, thermophysical properties have strong state dependence. These lead to various peculiar phenomena, one of which being the piston effect where a sudden heating induces a mechanical pulse. The coupling between thermal and mechanical processes, in the linear approximation, yields a non-trivially rich thermoacoustics. The numerous applications of supercritical fluids raise the need for reliable yet fast and efficient numerical solution for thermoacoustic time and space dependence in this sensitive domain. Here, we present a second-order accurate, fully explicit staggered space-time grid finite difference method for such coupled linear thermoacoustic problems. Time integration is based on the splitting of the state space vector field representing the interactions that affect the dynamics into reversible and irreversible parts, which splitting procedure leads to decoupled wave and heat equations. The former is a hyperbolic partial differential equation, while the latter is a parabolic one, therefore, different time integration algorithms must be amalgamated to obtain a reliable, dispersion error-free, and dissipation error-free numerical solution. Finally, the thermoacoustic approximation of the supercritical piston effect is investigated via the developed method.
翻译:在液-气临界点附近,超临界流体表现出强烈的可压缩性,同时其热物理性质具有显著的状态依赖性。这些特性导致多种特殊现象,其中之一为活塞效应——即突然加热会引发机械脉冲。在线性近似下,热过程与机械过程的耦合产生了非平凡且丰富的热声学现象。超临界流体的广泛应用催生了对此敏感域内热声时空依赖关系进行可靠、快速且高效数值求解的需求。本文针对此类线性耦合热声问题,提出了一种二阶精度、完全显式的交错时空网格有限差分方法。时间积分基于将表征动力学相互作用的状态空间矢量场分解为可逆与不可逆部分,该分裂过程导致波动方程与热传导方程的解耦。前者为双曲型偏微分方程,后者为抛物型方程,因此需要融合不同的时间积分算法以获得可靠、无频散误差且无耗散误差的数值解。最后,通过所开发的方法对超临界活塞效应的热声近似进行了研究。