两邻边铰支两邻边固支von-Karman矩形板的研究

本文运用双重Ｆｏｕｒｉｅｒ级数分析了受横向均布载荷作用下两邻边铰支两邻边固支矩形板的非线性弯曲问题。利用该级数正交性,可以把非线性偏微分控制方程化为含级数系数的非线性代数方程组并通过迭代解出所有系数。最后给出了位移和应力的计算结果。     

高压输电网设计中的二阶非线性方程的次谐波分枝

本文用Melnikov方法研究高压输电网设计中提出的非线性二阶受迫振动方程的次谐波分枝,分析了系统可能出现次谐解的条件和参数区域,以此为工程应用提供参考。     

旋转对管内流体流动的影响

本文描述了旋转管流的基本运动方程。在笛卡尔坐标系和柱坐标系中,针对两种基本的几何流动,推导出旋转情况下管内的流动方程,定义了影响旋转管流的流动和换热的两个无量纲准则数J和Ra,分析了二次流产生的原因,从而说明旋转对流动的影响以及旋转管流的运动特征。     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究,将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元,采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元,更好地表述了真实物面特征;物面边界相邻体网格单元相应构造为高阶体网格单元,其余体网格单元采用八节点六面体单元,以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解,验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外,由于采用了隐式计算方法,仅需较少的时间步就能迭代收敛。     

Automatic differentiation based discrete adjoint method for aerodynamic design optimization on unstructured meshes

A systematic methodology for formulating, implementing, solving and verifying discrete adjoint of the compressible Reynolds-averaged Navier-Stokes (RANS) equations for aerodynamic design optimization on unstructured meshes is proposed. First, a general adjoint formulation is con-structed for the entire optimization problem, including parameterization, mesh deformation, flow solution and computation of the objective function, which is followed by detailed formulations of matrix-vector products arising in the adjoint model. According to this formulation, procedural components of implementing the required matrix-vector products are generated by means of auto-matic differentiation (AD) in a structured and modular manner. Furthermore, a duality-preserving iterative algorithm is employed to solve flow adjoint equations arising in the adjoint model, ensur-ing identical convergence rates for the tangent and the adjoint models. A three-step strategy is adopted to verify the adjoint computation. The proposed method has several remarkable features:the use of AD techniques avoids tedious and error-prone manual derivation and programming;duality is strictly preserved so that consistent and highly accurate discrete sensitivities can be obtained; and comparable efficiency to hand-coded implementation can be achieved. Upon the cur-rent discrete adjoint method, a gradient-based optimization framework has been developed and applied to a drag reduction problem.     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

内埋式弹舱舱门气动载荷计算分析研究

采用中心有限体积法进行空间离散和四步Runge-Kutta显式推进方法进行时间推进数值求解了非定常欧拉方程;采用Delaunay方法生成内埋式弹舱舱盖打开时的非结构网格.数值模拟了内埋式轰炸机弹舱单开舱门在分别打开8°和45°时的复杂流场.通过对上述复杂流场的分析,发现作用于舱门的气动载荷具有使其趋向于关闭的特性,且有随打开角度的增大而减小的趋势.     

固体火箭发动机燃烧室三维流动数值计算

采用ｋε双方程湍流模型,以ＳＩＭＰＬＥ 计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流ＮＳ 方程．在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构．计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡．     

用非结构直角网格和欧拉方程计算飞机绕流

 采用八叉树结构 ,生成复杂外形绕流计算的非结构直角网格。物面附近用投影方法 ,使网格贴体。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。对歼击机模型的绕流计算表明 ,网格生成的机时花费很少 ,总体质量好 ,因而欧拉方程解算的收敛质量也好。