带副翼的翼身组合体绕流的Euler和N-S方程解

将对接分区网格与分区求解算法结合,有效地求解了带副翼偏转的翼身组合体绕流的N S方程.数值方法中选用VanLeer分裂格式离散无粘通量项,采用中心差分格式来离散粘性通量项.分区交界面采用了一种满足通量守恒的内边界耦合条件.数值算例表明该方法是求解带操纵面偏转的翼身组合体绕流的有效方法.     

一种基于AUSM思想的通量分裂方法

根据对流迎风分裂(AUSM)思想提出一种通量分裂方法,称为K-CUSP格式.它与传统H-CUSP和E-CUSP格式的最大差异在于总能量的分裂:K-CUSP格式将无粘守恒通量中所有的运动学量分裂到对流项,所有的热力学量分裂到压力项,即总能量被分裂成动能和静焓.对于压力项的数值通量,采用一种新的界面构造方法.数值测试表明:①K-CUSP格式继承了FVS格式的简单性和稳健性.在激波后不易出现压力过冲,在膨胀区域没有振荡,优于AUSM和WPS格式;②K-CUSP格式继承了FDS格式的分辨率.激波间断的分辨率和H-CUSP、E-CUSP格式基本相同,接触间断的分辨率高于FVS格式,低于Roe、AUSM和WPS格式.AUSM和WPS格式在计算运动接触间断时,速度存在很大振荡,而新格式不存在振荡.     

一种基于AUSM分裂的真正多维HLL格式

文章给出了一种真正多维的HLL Riemann解算器.采用AUSM分裂将通量分解成为对流通量和压力通量, 其中对流通量的计算采用迎风格式, 压力通量的计算采用HLL格式, 且将HLL格式的耗散项中的密度差用压力差代替, 从而使得格式能够分辨接触间断.为了实现数值格式真正多维的特性, 分别计算了网格界面中点和角点上的数值通量, 并且采用Simpson公式加权组合中点和角点上的数值通量得到网格界面的数值通量.为了减少重构角点处状态时的模板宽度, 计算中采用基于SDWLS梯度的线性重构获得2阶空间精度, 而时间离散采用2阶保强稳Runge-Kutta方法.数值实验表明, 相比于传统的一维HLL格式, 文章的真正多维HLL格式具有能够分辨接触间断, 以及更大的时间步长等优点.与其他能够分辨接触间断的格式(例如HLLC格式)不同, 真正多维的HLL格式在计算二维问题时不会出现激波不稳定现象.      

槽道湍流的展向振荡电磁力壁面减阻

采用直接数值模拟方法,对槽道湍流的展向振荡电磁力的减阻效果和减阻机理进行了研究,讨论了电磁力强度和振荡频率对湍流猝发事件以及壁面减阻率的影响.结果表明,电磁力强度或振荡频率变化时,湍流猝发频率和猝发强度的变化趋势是相反的,所以存在最优参数使得减阻效果最好.等价壁面展向速度可以很好地描述电磁力强度和振荡频率的变化对减阻效果的综合效应。     

抛物方程显隐修正迎风区域分裂差分法

给出抛物方程一种有效的区域分裂差分格式,提高了计算效率.对一阶项采用二阶迎风差分格式,内边界点和各子区域分别采用显隐差分格式.在较弱的稳定性条件下,得到离散l2模误差估计结果.最后给出具体的数值算例,以验证方法的实用性.     

气相爆轰高阶中心差分-WENO组合格式自适应网格方法

研究一种高阶中心差分-WENO组合格式,并采用自适应网格方法进行二维和三维气相爆轰波的数值模拟.采用ZND爆轰模型的控制方程为包含化学反应源项的Euler方程组.组合格式在大梯度区采用WENO格式捕捉间断,在光滑区采用高阶中心差分格式提高计算效率.采用一种基于流场结构特征的自适应网格.计算结果,表明这种方法同时具有高精度、高分辨率和高效率的特点.     

求解Burgers方程的一种交替分段隐式算法

首先提出一个新的求解Burgers方程的差分格式,然后在此差分格式的基础上构造了便于并行计算的交替分段隐格式,并作了线性化稳定性分析.数值结果表明,本方法具有较高的精度,尤适于扩散项系数较小时的计算,且有效避免了数值结果的非物理振荡.     

任意马赫数非定常流动数值模拟的统一算法

发展适用于从低速到高速任意马赫数非定常流动数值模拟的统一算法.通过引入一个伪时间导数项和一个新的预处理矩阵,得到双时间非定常预处理可压缩Navier-Stokes方程.方程的对流项采用三阶Roe通量近似差分格式离散,粘性项采用二阶中心差分格式离散.基于数值通量的线性化技术,实现伪时间步的隐式ADI-LU格式迭代,进而获得物理时间步的二阶推进精度.重点以低马赫数流动为例,求解了圆柱绕流和NACA0015翼型等速上仰动态失速问题.计算结果表明该统一算法能够较好地模拟低马赫数乃至任意马赫数非定常流动.     

多车种LWR交通流模型的半离散中心迎风格式

对多车种LWR交通流模型,给出一种半离散中心迎风格式,该格式以五阶WENO-Z重构和半离散中心迎风数值通量为基础.WENO-Z重构方法的引入提高了格式的精度,并保证格式具有基本无振荡的性质.时间的离散采用保持强稳定性的Runge-Kutta方法.通过数值算例验证了格式的有效性.     

一种新的可计算可压缩流动的预处理方法

针对使用可压缩流动数值方法求解不可压缩流动存在的刚性问题,基于虚拟压缩法思想,构造了一种以Mach数、速度、密度、温度等变量为元素的预处理矩阵,改变了控制方程组的特征根并使其量级更接近.通过理论推导与分析,证明新方法相比Weiss, Pletcher, Dailey和Choi的方法而言,不仅能降低方程组的刚性,提高了数值求解效率,而且拥有更好的稳定性,此外还能实现低速流动和高速流动之间的光滑过渡.采用有限差分格式进行离散,对流项的Roe格式作为基本加权无振荡(WENO)格式的求解器,黏性项则使用中心型紧致差分格式来计算,与预处理矩阵相结合展开数值实验,结果表明新预处理方法可以实现对无黏和有黏不可压缩流动问题的高精度模拟,且拥有比Weiss和Pletcher等提出的方法更好的收敛性和稳定性.     

一种基于多矩的有限体积浸入边界法

基于多矩VSIAM3格式及浸入边界法,提出一套在复杂计算区域内求解不可压缩流动的数值格式.不可压N-S方程使用VSIAM3格式进行离散,引入浸入边界法处理复杂、移动边界,使用虚拟网格方法计算动量方程修正项,同时还考虑了对连续方程的修正.使用标准算例对数值模式进行验证.     

磁各向异性色散介质散射的Padé时域有限差分方法分析

根据矩阵Padé逼近理论，把磁化色散介质的相对磁导率张量表示成以jω为自变量的矩阵函数形式，用/t代替jω后过渡到时域，再引入离散时域移位算子代替时间微分算子.进而导出磁化色散介质中的磁感应强度B和磁场强度H在离散时域的色散关系，并将其具体应用于旋磁介质，得到了这种磁化色散介质的Padé时域有限差分方法的递推表达式.作为验证，用这种方法计算了磁化铁氧体球的后向雷达散射截面，所得结果与文献结果一致.理论推导及算例表明，该方法是正确和有效的. 关键词： 各向异性介质 色散介质 矩阵Padé逼近 时域有限差分方法     

A deferred correction coupling strategy for low Mach number flow with complex chemistry

We present a new thermodynamic coupling strategy for complex reacting flow in a low Mach number framework. In such flows, the advection, diffusion and reaction processes span a broad range of time scales. In order to reduce splitting errors inherent in Strang splitting approaches, we couple the processes with a multi-implicit spectral deferred correction strategy. Our iterative scheme uses a series of relatively simple correction equations to reduce the error in the solution. The new method retains the efficiencies of Strang splitting compared to a traditional method-of-lines approach in that each process is discretised sequentially using a numerical method well suited for its particular time scale. We demonstrate that the overall scheme is second-order accurate and provides increased accuracy with less computational work compared to Strang splitting for terrestrial and astrophysical flames. The overall framework also sets the stage for higher-order coupling strategies.     

Modelling two-phase flow in porous media at the pore scale using the volume-of-fluid method

We present a stable numerical scheme for modelling multiphase flow in porous media, where the characteristic size of the flow domain is of the order of microns to millimetres. The numerical method is developed for efficient modelling of multiphase flow in porous media with complex interface motion and irregular solid boundaries. The Navier–Stokes equations are discretised using a finite volume approach, while the volume-of-fluid method is used to capture the location of interfaces. Capillary forces are computed using a semi-sharp surface force model, in which the transition area for capillary pressure is effectively limited to one grid block. This new formulation along with two new filtering methods, developed for correcting capillary forces, permits simulations at very low capillary numbers and avoids non-physical velocities. Capillary forces are implemented using a semi-implicit formulation, which allows larger time step sizes at low capillary numbers. We verify the accuracy and stability of the numerical method on several test cases, which indicate the potential of the method to predict multiphase flow processes.     

Analysis of Excitations for a Groove Guide Resonator at 10 GHz by means of the FDTD Method

Various excitations of a new groove guide resonator working in the X-band (8 GHz – 12 GHz) are investigated by means of numerical simulations. For the numerical simulations the Finite-Difference Time-Domain method is used. The groove guide resonator, modelled both with and without excitation structures, is discretised in space. The results in the time domain are then transformed into the frequency domain in order to obtain the resonance frequency spectrum. Comparison between simulations with and without excitations shows the effect of the excitations on the resonance frequency spectrum. The results are compared with those of previous analytical methods.     

A robust numerical scheme for highly compressible magnetohydrodynamics: Nonlinear stability,implementation and tests

The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate that the theoretical stability properties are valid in practice with negligible compromises to accuracy. The result is a highly robust scheme with state-of-the-art efficiency. The scheme’s robustness is due to entropy stability, positivity and properly discretised Powell terms. The implementation takes the form of a modification of the MHD module in the FLASH code, an adaptive mesh refinement code. We compare the new scheme with the standard FLASH implementation for MHD. Results show comparable accuracy to standard FLASH with the Roe solver, but highly improved efficiency and stability, particularly for high Mach number flows and low plasma β. The tests include 1D shock tubes, 2D instabilities and highly supersonic, 3D turbulence. We consider turbulent flows with RMS sonic Mach numbers up to 10, typical of gas flows in the interstellar medium. We investigate both strong initial magnetic fields and magnetic field amplification by the turbulent dynamo from extremely high plasma β. The energy spectra show a reasonable decrease in dissipation with grid refinement, and at a resolution of 512³ grid cells we identify a narrow inertial range with the expected power law scaling. The turbulent dynamo exhibits exponential growth of magnetic pressure, with the growth rate higher from solenoidal forcing than from compressive forcing. Two versions of the new scheme are presented, using relaxation-based 3-wave and 5-wave approximate Riemann solvers, respectively. The 5-wave solver is more accurate in some cases, and its computational cost is close to the 3-wave solver.     

Large-scale turbulence modelling via α-regularisation for atmospheric simulations

We study the possibility of obtaining a computational turbulence model by means of non-dissipative regularisation of the compressible atmospheric equations for climate-type applications. We use an α-regularisation (Lagrangian averaging) of the atmospheric equations. For the hydrostatic and compressible atmospheric equations discretised using a finite volume method on unstructured grids, deterministic and non-deterministic numerical experiments are conducted to compare the individual solutions and the statistics of the regularised equations to those of the original model. The impact of the regularisation parameter is investigated. Our results confirm the principal compatibility of α-regularisation with atmospheric dynamics and encourage further investigations within atmospheric model including complex physical parametrisations.     

Wave-induced bed flows by a Lagrangian vortex scheme

Nominally 2-dimensional viscous flow induced by gravity waves over a spatially periodic bed is simulated by a Lagrangian vortex scheme. A vortex sheet is introduced on the surface at each time step to satisfy the zero velocity conditions. The sheet is discretised; the vortex-in-cell method is used to convect vorticity and random walks are added to effect viscous diffusion. Good agreement with analytical theory is obtained for velocity profiles in uniform sinusoidal flow and for mass transport due to linear waves. Mass transport for finite amplitude waves is also obtained. For separated flow over rippled beds, which is still liminar, a vortex decay factor is required to produce agreement with experiment and is thought to compensate for large scale 3-dimensional effects.