Development of accurate, robust liquid equations of state for multi-phase CFD simulations with a modified AUSM+-up scheme

Eight equations of state (EOS) have been evaluated for the simulation of compressible liquid water properties, based on empirical correlations, the principle of corresponding states and thermodynamic relations. The IAPWS-IF97 EOS for water was employed as the reference case. These EOSs were coupled to a modified AUSM+-up convective flux solver to determine flow profiles for three test cases of differing flow conditions. The impact of the non-viscous interaction term discretisation scheme, interfacial pressure method and selection of low-Mach number diffusion were also compared. It was shown that a consistent discretisation scheme using the AUSM+-up solver for both the convective flux and the non-viscous interfacial term demonstrated both robustness and accuracy whilst facilitating a computationally cheaper solution than discretisation of the interfacial term independently by a central scheme. The simple empirical correlations gave excellent results in comparison to the reference IAPWS-IF97 EOS and were recommended for developmental work involving water as a cheaper and more accurate EOS than the more commonly used stiffened-gas model. The correlations based on the principles of corresponding-states and the modified Peng-Robinson cubic EOS also demonstrated a high degree of accuracy, which is promising for future work with generic fluids. Further work will encompass extension of the solver to multiple dimensions and to account for other source terms such as surface tension, along with the incorporation of phase changes. © 2013.

CFD investigation of turbulent premixed flame response to transverse forcing

Screech is a high frequency oscillation that is usually characterized by instabilities caused by large-scale coherent flow structures in the wake of bluff-body flameholders and shear layers. Such oscillations can lead to changes in flame surface area which can cause the flame to burn unsteadily, but also couple with the acoustic modes and inherent fluid-mechanical instabilities that are present in the system. In this study, the flame response to hydrodynamic oscillations is analyzed in a controlled manner using high-fidelity Computational Fluid Dynamics (CFD) with an unsteady Reynolds-averaged Navier-Stokes approach. The response of a premixed flame with and without transverse velocity forcing is analyzed. When unforced, the flame is shown to exhibit a self-excitation that is attributed to the anti-symmetric shedding of vortices in the wake of the flameholder. The flame is also forced using two different kinds of low-amplitude out-of-phase inlet velocity forcing signals. The first forcing method is harmonic forcing with a single characteristic frequency, while the second forcing method involves a broadband forcing signal with frequencies in the range of 500 - 1000 Hz. For the harmonic forcing method, the flame is perturbed only lightly about its mean position and exhibits a limit cycle oscillation that is characteristic of the forcing frequency. For the broadband forcing method, larger changes in the flame surface area and detachment of the flame sheet can be seen. Transition to a complicated trajectory in the phase space is observed. When analyzed systematically with system identification methods, the CFD results, expressed in the form of the Flame Transfer Function (FTF) are capable of elucidating the flame response to the imposed perturbation. The FTF also serves to identify, both spatially and temporally, regions where the flame responds linearly and nonlinearly. Locking-in between the flame's natural self-excited frequency and the subharmonic frequencies of the broadband forcing signal is found to alter the dynamical behaviour of the flame. Copyright © 2013 by ASME.

数值摄动算法及其CFD格式

数值摄动算法将流体动力学效应耦合进NS方程组和对流扩散(CD)方程离散的数学基本格式(MBS),特别是耦合进最简单的一阶迎风和二阶中心格式之中,由此构建成一系列新的摄动格式(PS).构建PS的主要步骤是将MBS中的通量重构为步长的幂级数,利用空间分裂和导出的高阶流体动力学线性关系式,并引入下游不影响上游的对流运动规律,通过消除重构格式修正微分方程的截断误差诸项求出幂级数的待定系数,由此获得非线性PS.PS的项是MBS中对应项与R△x(及λR△x)之简单多项式的乘积,R△x和λ分别是网格Reynolds数和网格CFL数.PS和MBS使用相同结点,简单性彼此相当,但PS精度高,稳定范围大,例如PS包含了许多绝对稳定高阶迎风和中心有限差分(FD)格式和绝对正型有限体积(FV)格式,这些格式对网格Reynolds数的任意值均为不振荡格式.数值摄动算法因此是构建高精度不振荡CFD格式的新方法.PS用于计算不可压缩流,可压缩流,液滴萃取传质,微通道两相流等,均获得良好数值结果或与已有Benchmark解一致的数值结果.已有文献称数值摄动算法为新型高精度方法和高算法,文中也讨论了一些值得进一步研究的课题

基于并行CFD方法的颤振特性计算

将流体和结构运动方程分别构造为含子迭代的计算格式,发展一种紧耦合气动弹性分析方法。其中流体计算的空间离散采用改进的HLLEW(Harten-Lax-Van Leer-Einfeldt-Wada)格式,同时采用MPI通信实现并行计算。弹性体网格变形法用于生成随结构变形的四面体非结构网格。利用所发展的方法,对一标模445.6翼型进行了跨声速气动颤振特性研究,分析其颤振特性并对比实验结果。同时还运用所发展的方法对某机带控制舵垂尾的颤振特性进行了研究,确立了所发展方法对实际工程颤振问题的适应性。

CFD analysis of thermodynamic cycles in a pulse tube refrigerator

The objectives of this paper are to study the thermodynamic cycles in an inertance tube pulse tube refrigerator (ITPTR) by means of CFD method The simulation results show that gas parcels working in different parts of ITPTR undergo different thermodynamic cycles The net effects of those thermodynamic cycles are pumping heat from the low temperature part to the high temperature part of the system The simulation results also show that under different frequencies of piston movement the gas parcels working in the same part of the system will undergo the same type of thermodynamic cycles The simulated thermal cycles are compared with those thermodynamic analysis results from a reference Comparisons show that both CFD simulations and theoretical analysis predict the same type of thermal cycles at the same location However only CFD simulation can give the quantitative results while the thermodynamic analysis is still remaining in quality (C) 2010 Elsevier Ltd All rights reserved

Sand storms: CFD analysis of Reynolds stress and collision stress of particles near sand bed

Sand storm is a serious environmental threat to humans. Sand particles are transported by saltation and suspension, causing soil erosion in one place and deposition in another. In order to prevent and predict sand storms, the causes and the manners of particle motions must be studied in detail. In this paper a standard k-epsilon model is used for the gas phase simulation and the discrete element method (DEM) is used to predict the movements of particles using an in-house procedure. The data are summarized in an Eulerian-Eulerian regime after simulation to get the statistical particle Reynolds stress and particle collision stress. The results show that for the current case the Reynolds stress and the air shear stress predominate in the region 20-250 mm above the initial sand bed surface. However, in the region below 3 mm, the collision stress must be taken into account in predicting particle movement. (C) 2010 Chinese Society of Particuology and Institute of Process Engineering, Chinese Academy of Sciences. Published by Elsevier B.V. All rights reserved.