共查询到20条相似文献,搜索用时 78 毫秒
1.
2.
新的解跨音速Euler方程的隐式杂交方法 总被引:5,自引:0,他引:5
本文提出一个新的杂交格式,它由Jameson-Turkel LU分解法和Harten TVD(total variation diminishing)格式的修正数值通量耦合而成。LU分解用于Euler方程隐式部分的离散化而Harten的数值通量控制所谓的稳态解残差。精心构造LU算子,可使矩阵算子L和U很容易逐点推进求逆;以Harten的数值通量计算杂交格式右边的通量,这样做在a=1/2的条件下,可保证收敛的稳态解具有二阶精度并能高分辨率地捕获激波。数值实验表明,本文格式对于稳态计算是有生命力的并能保证较高的激波分辨率。 相似文献
3.
计算三维粘性/无粘性干扰流动的准联立分区算法 总被引:1,自引:1,他引:1
应用摄动理论导出了粘性/无粘性流动的匹配条件,与边界层方程联立形成能计及粘性/无粘性干扰效应的模型方程。同时引用等效吮吸理论对该模型方程与全位势方程进行交替迭代求解,建立了计算高雷诺数流动的一种快速收敛的粘性/无粘性分区算法。计算结果表明该算法适用于包括小尺度分离的流动;与传统的半逆解法相比,计算效率明显提高 相似文献
4.
介绍了一种基于经验模式分解(Empirical Mode Decomposition,EMD)与粒子群优化算法相结合的飞机结构模态参数辨识方法。一个复杂的脉冲响应信号利用EMD方法使得耦合在一起的多阶模态响应信号分解为与各单阶模态响应信号一一对应的分量,得到前几阶主要的内禀模式函数(Intrinsic Mode Function,IMF),再对分解得到的每一单阶IMF利用粒子群算法辨识得到各阶模态参数。试验仿真结果表明该方法有较高的计算精度,可应用于结构运行模态分析,为飞机等结构设计、运行检测提供有力保障。 相似文献
5.
基于谱方法的管路充填过程仿真 总被引:3,自引:1,他引:3
基于一维瞬变流理论,建立了推进剂供应系统管路充填过程数学模型.提出采用Fourier谱方法求解充填管道中液体瞬变流控制方程的新方法.利用坐标变换方法,很好地处理了解算管内气液移动边界时,非线性迭代收敛速度较慢这一难题.对由贮箱、隔离阀、管道组成的系统充填过程进行了仿真计算,并将计算结果与已发表的采用特征线方法获得的结果进行对比分析,表明该方法和仿真结果可信.进一步研究了预存气体压力、贮箱压力、多变指数等影响因素对最大充填压力峰的影响规律. 相似文献
6.
7.
本文给出了一种求解二维可压Navier-Stokes(简称N-S)方程的边界元法。使用线化技术把控制方程变成变系数的线性偏微分方程,其基本解由组合代数的方法构造,从而给出了N-S方程的积分表达式。本文计算了几例低R_e数的跨声速翼型绕流,结果表明边界元法可以应用到可压粘流中。 相似文献
8.
基于控制理论的Euler方程翼型减阻优化设计 总被引:3,自引:0,他引:3
在Jameson的控制论气动优化设计思想下,应用Euler方程研究了翼型的反设计和减阻问题。设计过程中的梯度是通过求解一个伴随偏微分方程而得到的,每个设计迭代只需求解一次Euler方程和一次伴随方,一与设计变量数无关。在具体数值执行中:①利用分部积分将目标函数变分的计算转化成了类似于计算Euler方程残值的形式,节省了机器时间;②根据控制理论的要求,将减阻问题转化成了相应地修改伴随方程的物面边界条件和目标函数的变分表达,使问题得到了简化;③利用特征不变量分析法,处理了伴承方程的物面和远场边界条件。设计算例证明了本文方法可靠性好、收敛快、大大节约设计时间。 相似文献
9.
本文从机场跑道不平度准则出发,基于功率谱概念,介绍了测定跑道不平度随机变量的基本方法.对跑道不平度功率谱曲线的拟合和飞机地面滑行动力响应进行初步分析。 相似文献
10.
低速叶型气动反问题设计方法 总被引:1,自引:1,他引:1
低马赫数不可压流动中声速与流速大小差别巨大,采用基于可压缩流动控制方程的计算格式求解流场时,由于数值黏性的污染,解的精度低且收敛性差,通常可使用时间预处理技术来解决这一问题。在基于控制理论的优化方法中,共轭方程的Jacobian矩阵和流动方程的系数矩阵相似,因而在低流动马赫数下,求解共轭方程存在着与求解流动方程相同的数值污染和数值刚性问题。首先推导了带有预处理的Roe格式,然后发展了适合全速度流动的共轭方程求解方法,最后选取翼型和叶栅两个典型算例进行了验证。计算结果表明所发展的方法可很好地用于低马赫数时的气动反问题设计。 相似文献
11.
12.
一种计算非定常二维流动的无网格算法 总被引:1,自引:0,他引:1
主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。 相似文献
13.
提出了一种三维动态非结构重叠网格Navier-Stokes(N-S)方程的并行计算方法。N-S方程的空间离散采用格点有限体积方法,时间离散采用隐式的双时间步长方法。应用一方程Spalart-Allmaras(S-A)模型来计算湍流黏性。并行计算采用动态的区域分裂方法,在每一物理时间步利用METIS网格分区系统对网格进行分区。为了实现各CPU之间的负载均衡,每块网格都按CPU个数进行分区并对活动节点和非活动节点进行了加权处理。最后,通过对外挂物投放无黏流动的数值模拟和内埋武器弹舱开启黏性流动的数值模拟,验证了该并行程序的准确性、高性能并行计算以及处理复杂几何外形的能力。 相似文献
14.
The very high accuracy of the Doppler and range measurements between the two low-flying and co-orbiting spacecraft of the
GRACE mission, which will be at the μm/sec and ≈10 μm levels respectively, requires that special procedures be applied in
the processing of these data. Parts of the existing orbit determination and gravity field parameters retrieval methods and
software must be modified in order to fully benefit from the capabilities of this mission. This is being done in the following
areas: (i) numerical integration of the equations of motion (summed form, accuracy of the predictor-corrector loop, Encke's
formulation): (ii) special inter-satellite dynamical parameterization for very short arcs; (iii) accurate solution of large
least-squares problems (normal equations vs. orthogonal decomposition of observation equations); (iv) handling the observation
equations with high accuracy. Theoretical concepts and first tests of some of the newly implemented algorithms are presented.
This revised version was published online in August 2006 with corrections to the Cover Date. 相似文献
15.
16.
提出了近壁湍流对称单个相干结构的理论模型,采用紧致有限差分和Fourier谱展开相结合的数值方法研究它的演化问题,得到雷诺应力分布及高剪切层的形成过程与实验结果一致。 相似文献
17.
Godunov方法在复杂流场数值模拟中的应用 总被引:5,自引:3,他引:2
本文根据Van Leer提出的方法给出了求解三维非定常Euler方程和三维定常Euler方程的二阶Godunov方法。方法要点为:在各体积离散单元中,用流动变量的分段线性分布代替一阶方法中的分段常量分布,在确定线性分布的斜率中引入单调性限制条件;在离散单元界面构造Riemann问题,由Riemann间断解确定单元侧面的通量;用预测-修正两步法对离散方程推进求解,用上述方法结合单区网格或多区重迭网格 相似文献
18.
19.
The ability to generate rational models of time series plays an important role in such applications as adaptive filtering, spectral estimation, digital control, array processing, and forecasting. A method for effecting an autoregressive moving average (ARMA) model estimate is presented which possesses a number of admirable properties: 1) it has an elegant algebraic structure, 2) its modeling performance in spectral estimation applications has been empirically found to typically exceed that of such contemporary techniques as the periodogram, the Burg method, and the Box-Jenkins method on a variety of problems, 3) it is implementable by computationally efficient algorithms, and 4) it is based on pseudomaximum likelihood concepts. Taken in combination, these properties mark this method as being an effective tool in challenging applications requiring high modeling performance in a real time setting. 相似文献