共查询到20条相似文献,搜索用时 171 毫秒
1.
流动稳定性问题常常归结于巨型非对称矩阵特征值问题。多数求解巨型非对称矩阵特征问题的算法均是经基本的Arnoldi算法演化而来。首先简述基本的Arnoldi算法; 其次简述基于Arnoldi算法的几类变体, 如显式重启Arnoldi算法,隐式重启Arnoldi算法与多重隐式重启Arnoldi算法; 最后基于Arnoldi算法及其变体结合谱位移技术求解计算流动稳定性问题, 并通过数值实验比较可知结合谱位移技术的多重隐式重启Arnoldi算法的求解效率最高。 相似文献
2.
本文对谱方法用于周期性非定常流动的隐式求解方法进行了探讨,分析了影响计算稳定性和收敛速度的因素.提出了结合多重网格的隐式求解方法并对算法进行了验证,初步计算表明本文算法具有良好的稳定性和收敛速度.对于周期性非定常流动,结合本文提出的隐式求解的时域谱方法可以达到很高的精度且具有良好的计算效率. 相似文献
3.
树形多体Hamilton系统辛算法 总被引:3,自引:0,他引:3
研究了树形多体Hamilton系统的隐式辛算法。用矩阵形式给出了系统的正则方程及其右端函数的Jacobi矩阵,并给出该矩阵的分块算法,可提高计算效率。隐式辛Runge-Kutta算法被采用,数值结果表明给出的算法计算效率高,并可保持长期数值计算的稳定性。 相似文献
4.
5.
6.
7.
针对现有可见光隐式成像通信调制算法隐式效果差、算法复杂的问题,提出了一种基于蓝色通道拉普拉斯最高层金字塔的隐式信息调制算法.该算法利用人眼视觉系统对黄蓝分量不敏感的特性,仅对蓝色通道图像进行拉普拉斯金字塔分解,并在最高层嵌入隐式信息.此外,考虑到图像纹理差异对隐式效果的影响,采用图像熵作为衡量图像间纹理复杂程度的指标,... 相似文献
8.
光伏电池组件非线性输出特性的物理建模及其优化参数的准确提取是光伏发电系统设计计算、性能评估及优化控制的重要前提.相对于传统的隐式单二极管模型,该文在光伏电池显式单二极管模型的基础上利用Lambert W函数推导了光伏组件的显式单二极管模型,提出一种基于重启边界约束Nelder-Mead单纯形算法的参数提取方法rbcNM,并利用两种典型光伏电池组件的实测数据对隐式、显式单二极管模型的准确性进行了对比测试和验证.结果表明:rbcNM算法可以快速准确的提取隐式、显式单二极管模型的优化参数,计算结果与实测数据具有很好的一致性,相对于已有文献在准确度上取得了大幅度的提升;显式单二极管模型的准确性显著高于隐式单二极管模型,对光伏电池组件的电流-电压和功率-电压特性曲线具有更高的拟合精度. 相似文献
9.
10.
一种能克服两方程湍流模型刚性的流热一体化算法 总被引:1,自引:0,他引:1
使用流热一体化算法解决传统耦合算法难于处理的流场与固体温度场耦合计算问题.为保证其计算效率,提出一个简单有效的方法用以克服两方程湍流模型刚性.给出该方法的理论证明,并验证其配合多重网格提高LUSGS隐式时间推进法的计算效率. 相似文献
11.
A new high order finite-difference method utilizing the idea of Harten ENO subcell resolution method is proposed for chemical reactive flows and combustion. In reaction problems, when the reaction time scale is very small, e.g., orders of magnitude smaller than the fluid dynamics time scales, the governing equations will become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present proposed method is a modified fractional step method which solves the convection step and reaction step separately. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with the computed flow variables in the shock region modified by the Harten subcell resolution idea. For numerical experiments, a fifth-order finite-difference WENO scheme and its anti-diffusion WENO variant are considered. A wide range of 1D and 2D scalar and Euler system test cases are investigated. Studies indicate that for the considered test cases, the new method maintains high order accuracy in space for smooth flows, and for stiff source terms with discontinuities, it can capture the correct propagation speed of discontinuities in very coarse meshes with reasonable CFL numbers. 相似文献
12.
13.
14.
用分步傅里叶变换法求解二阶孤子传输的非线性薛定谔方程, 得到了在此条件下孤子传输的数值图形, 发现二阶孤子在传输中被压缩, 幅值振荡变化。2个二阶孤子在传输过程中没有出现象2个一阶孤子那样周期性碰撞, 但2个二阶孤子时间间隔较小时, 随传输距离在2个二阶孤子中间周期性地衍生出第3个孤子。研究证明:二阶孤子的传输具有与一阶孤子明显不同的特征。 相似文献
15.
I. K. Larin M. L. Kuskov 《Russian Journal of Physical Chemistry B, Focus on Physics》2013,7(4):509-513
The characteristics of the chain processes of stratospheric ozone depletion are considered. It is shown that the basic feature of these processes is the presence of a limiting stage, which determines the rate of chain propagation in this process and, accordingly, the rate of ozone depletion. It is also shown that the practice, so far widespread in the stratospheric chemistry, of defining the rate-limiting step of the chain process as a single reaction with the lowest rate throughout the stratosphere does not enable to correctly determine the rare of the chain process and leads to a significant overestimation of the latter. Methods for correctly calculating the rate of the limiting step for an arbitrary number of chain propagation reactions and for determining the termination rate and the chain length are for the first time proposed. Application of the developed method is demonstrated by the example of the nitrogen oxide cycle. 相似文献
16.
17.
A pressure based, iterative finite volume method is developed for calculation of compressible, viscous, heat conductive gas flows at all speeds. The method does not need the use of under-relaxation coefficient in order to ensure a convergence of the iterative process. The method is derived from a general form of system of equations describing the motion of compressible, viscous gas. An emphasis is done on the calculation of gaseous microfluidic problems. A fast transient process of gas wave propagation in a two-dimensional microchannel is used as a benchmark problem. The results obtained by using the new method are compared with the numerical solution obtained by using SIMPLE (iterative) and PISO (non-iterative) methods. It is shown that the new iterative method is faster than SIMPLE. For the considered problem the new method is slightly faster than PISO as well. Calculated are also some typical microfluidic subsonic and supersonic flows, and the Rayleigh–Bénard convection of a rarefied gas in continuum limit. The numerical results are compared with other analytical and numerical solutions. 相似文献
18.
A new fractional-step method is proposed for the numerical solution of high speed reacting flows, where the chemical time scales are often much smaller than the fluid dynamical time scales. When the problem is stiff, because of insufficient spatial/temporal resolution, a well-known spurious numerical phenomenon occurs in standard finite volume schemes: the incorrect calculation of the speed of propagation of discontinuities. The new method is first illustrated considering a one-dimensional scalar hyperbolic advection/reaction equation with stiff source term, which may be considered as a model problem to under-resolved detonations. During the reaction step, the proposed scheme replaces the cell average representation with a two-value reconstruction, which allows us to locate the discontinuity position inside the cell during the computation of the source term. This results in the correct propagation of discontinuities even in the stiff case. The method is proved to be second-order accurate for smooth solutions of scalar equations and is applied successfully to the solution of the one-dimensional reactive Euler equations for Chapman–Jouguet detonations. 相似文献
19.
20.
Tindle CT O'Driscoll LM Higham CJ 《The Journal of the Acoustical Society of America》2000,108(1):76-83
The conventional coupled mode solution is combined with perturbation theory to give a fast, accurate range-dependent normal mode solution for deep water acoustic propagation. Perturbation theory is used to calculate the new normal modes at each range step. The new modes are obtained as a linear combination of the modes for the previous step without requiring a numerical solution of the depth-separated wave equation. The process may be repeated for many steps and yields normal modes and eigenvalues which are sufficiently accurate for solution of practical problems in deep water. The method is applied to long-range propagation through oceanic fronts. 相似文献