共查询到20条相似文献,搜索用时 140 毫秒
1.
2.
本文采用直接解N-S方程的方法求解振荡叶栅的非定常流场.数值格式采用MC-LB显隐格式,时间步长可以取得较大,故总计算时间并不比用显格式直接解Euler方程长很多.计算结果与无粘结果基本一致.但来流时均静压比定常来流静压略有提高,这是与无粘计算不同之处. 相似文献
3.
4.
根据对流迎风分裂(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格式在计算运动接触间断时,速度存在很大振荡,而新格式不存在振荡. 相似文献
5.
在一个单步格式的基础上构造了一个新的隐式系数矩阵分裂法来数值求解可压Navier-St-okes方程。对方程中的无粘项部分利用守恒型方程中流通量向量为一齐次函数这一特性,根据Jacobian矩阵特征值的符号而将流通向量分裂成两部分。在此基础上据风向而构造逼近于无粘项的差分格式。对方程中的粘性项部分利用算子附加修正的方法来改进计算的收敛过程。所建立的差分格式被用来数值求解Couette流以考查这一方法。 相似文献
6.
7.
为高精度捕捉激波等流场结构,引入一种Power限制器,对NND格式和WNND格式进行改进,分别得到二阶PNND(Power NND)格式和三阶PWNND(Power WNND)格式.该Power类型格式通过Power限制器对相邻待选模板上的一阶导数进行限制,改善了NND格式和WNND格式在间断附近的耗散效应.对各种格式的分析表明,在间断附近采用Power限制器的格式比原格式的表现要好,耗散小且捕捉间断精度高,其中PNND格式虽然只有二阶精度,但在所有算例中与三阶WNND格式的计算结果比较接近,在个别算例中甚至优于WNND格式.最后将PWNND格式应用到二维NACA0012翼型的强迫俯仰振动的数值模拟,计算结果与实验值、参考计算值吻合较好. 相似文献
8.
Ghost Fluid方法与双介质可压缩流动计算 总被引:1,自引:1,他引:0
应用带有Isobaric修正的GhostFluid方法配合LevelSet方法计算可压缩双介质无粘流动.该方法可以消除计算流体界面时所产生的数值跳动和耗散,且编程上比界面跟踪法简单.应用WENO格式数值求解欧拉方程和LevelSet方程,对由刚性气体状态方程所支配的一二维双介质流动进行数值计算,得到了分辨率较高的计算结果. 相似文献
9.
10.
空间-时间守恒(STC)格式是近年来发展出的一种计算格式,在现有的STC格式构造过程中,流动变量在解元中的分布都用其一阶Taylor展开式来表示.STC格式的精度与所采用的Taylor展开式的阶数有关.该文采用流动变量的二阶Taylor展开式来表示其在解元上的分布、构造出了求解一维Euler方程的STC格式.用该格式对几个问题进行了计算,将计算结果与精确解进行了比较,比较表明该格式有较高的精度. 相似文献
11.
主要研究了气动噪声对光传输和成像的影响。流场求解采用法福尔平均N-S方程,并分别利用基于MUSCL格式的AUSMPW+分裂和二阶中心格式对流通量和粘性通量进行离散处理;用四步四阶Runge-Kutta格式对方程进行时间推进求解,采用J-B两方程模型并经Sarkar方法进行可压缩修正后用于湍流的数值模拟。采用基于光学传递函数的物理光学方法分析平均流动和湍流对光质的影响。研究表明:涡脉冲在运动过程中会引起光束沿涡脉冲运动方向产生较大的整体偏折,偏折幅度随时间而变,最大偏折幅度可以达到1.7 μrad;在与涡脉冲运动的垂直方向,整体偏折基本可以忽略不记。在涡脉冲横穿光束的过程中,平均流动会导致焦平面光强峰值有两次下降,平均流动对光强重新分布的影响比湍流脉动影响更强。涡脉冲在运动过程中,先导致光束能量在与涡脉冲运动垂直的方向上扩散,随后使之在涡脉冲运动方向上向外转移。因能量向外扩散,导致Strehl比在0.9 s时刻降低到0.60左右。 相似文献
12.
对一维双曲型守恒律,给出了一种具有较小数值耗散的三阶半离散中心迎风格式.该格式以Liu和Tadmor提出的三阶无振荡重构为基础,同时考虑了波传播的单侧局部速度.时间离散用保持强稳定性的三阶Runge-Kutta方法.由于不需用Riemann解算器,避免了特征分解过程,保持了中心格式简单的优点.数值算例验证本方法可进一步减小数值耗散,提高分辨率. 相似文献
13.
14.
A scheme for third-order ghost interference with thermal light is proposed. The visibility and resolution of the interference fringe related to the bandwidth of the spatial frequency spectrum of the source are analysed. The results show that the visibility of the third-order ghost interference fringe is much higher than that of the second-order one. 相似文献
15.
频率锁定多载波光源是实现太比特系统传输的关键器件之一, 目前已成为光纤通信领域的研究热点.基于循环频移技术产生的频率锁定 多载波光源因具有驱动电压低、串扰小、平坦度高的优点而受到重视, 但是其输出质量容易受到调制器固有的高阶谐波串扰的影响,特别是三阶谐波串扰. 为了减小三阶串扰的影响,本文提出了在I/Q调制器两支路上再加载另一频率为3fm的射频信号的方案, 用其产生的一阶信号来抑制三阶串扰.通过理论分析与数值仿真,得到了载波数为24、 载波间最大功率差小于0.1 dB的多载波光源输出.与未采用三阶串扰抑制方案的结果比较, 光源输出有效光信噪比提高了2 dB. 结果表明,采用本文的改进方案,可以有效地提高多载波光源的输出平坦度和有效光信噪比. 相似文献
16.
A new third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme for scalar and vector hyperbolic equations with piecewise continuous initial conditions is developed. The new scheme is proven to be linearly stable in the energy norm for both continuous and discontinuous solutions. In contrast to the existing high-resolution shock-capturing schemes, no assumption that the reconstruction should be total variation bounded (TVB) is explicitly required to prove stability of the new scheme. We also present new weight functions which drastically improve the accuracy of the third-order ESWENO scheme. Based on a truncation error analysis, we show that the ESWENO scheme is design-order accurate for smooth solutions with any number of vanishing derivatives, if its tuning parameters satisfy certain constraints. Numerical results show that the new ESWENO scheme is stable and significantly outperforms the conventional third-order WENO scheme of Jiang and Shu in terms of accuracy, while providing essentially non-oscillatory solutions near strong discontinuities. 相似文献
17.
By numerical simulations, we show that picosecond Gaussian optical pulses with a precise optimal frequency chirping can transmit stably in full-dispersion compensation optical fiber links, with not only second- but also third-order dispersion compensation, using dispersion shift fibers with opposite dispersion sign. The optimal pre-chirp is determined principally by the second-order dispersion scheme and scarcely affected by third-order dispersion scheme. It demonstrates that, to a high bit rate transmission system, the pre-chirping technology and higher-order dispersion compensation are two very efficient measures in improving performance of system. 相似文献
18.
In the paper, we provide sufficient conditions for the oscillatory and asymptotic behavior of a new type of third-order nonlinear dynamic equations with mixed nonlinear neutral terms. Our theorems not only improve and extend existing theorems in the literature but also provide a new approach as far as the nonlinear neutral terms are concerned. The main results are illustrated by some particular examples. 相似文献
19.
《中国物理 B》2015,(5)
This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time(SRT) lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations,it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version. 相似文献
20.
Shravan K. Veerapaneni Denis Gueyffier George Biros Denis Zorin 《Journal of computational physics》2009,228(19):7233-7249
We extend [Shravan K. Veerapaneni, Denis Gueyffier, Denis Zorin, George Biros, A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D, Journal of Computational Physics 228(7) (2009) 2334–2353] to the case of three-dimensional axisymmetric vesicles of spherical or toroidal topology immersed in viscous flows. Although the main components of the algorithm are similar in spirit to the 2D case—spectral approximation in space, semi-implicit time-stepping scheme—the main differences are that the bending and viscous force require new analysis, the linearization for the semi-implicit schemes must be rederived, a fully implicit scheme must be used for the toroidal topology to eliminate a CFL-type restriction and a novel numerical scheme for the evaluation of the 3D Stokes single layer potential on an axisymmetric surface is necessary to speed up the calculations. By introducing these novel components, we obtain a time-scheme that experimentally is unconditionally stable, has low cost per time step, and is third-order accurate in time. We present numerical results to analyze the cost and convergence rates of the scheme. To verify the solver, we compare it to a constrained variational approach to compute equilibrium shapes that does not involve interactions with a viscous fluid. To illustrate the applicability of method, we consider a few vesicle-flow interaction problems: the sedimentation of a vesicle, interactions of one and three vesicles with a background Poiseuille flow. 相似文献