三维可压缩 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.     

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

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

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

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

三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

A PARALLEL IMAGINARY EBE METHOD FOR SOLVING POSITIVE DEFINITE LINEAR SYSTEMS

ＡＰＡＲＡＬＬＥＬＩＭＡＧＩＮＡＲＹＥＢＥＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧＰＯＳＩＴＩＶＥＤＥＦＩＮＩＴＥＬＩＮＥＡＲＳＹＳＴＥＭＳＰａｎＸｉａｏｓｕ（ＤｅｐａｒｔｍｅｎｔｏｆＡｐｐｌｉｅｄＭａｔｈｅｍａｔｉｃｓ，ＰｈｙｓｉｃｓａｎｄＭｅｃｈａｎｉｓ...     

角接触球轴承分析模型的数值求解

建立了内圈静止、外圈旋转的角接触球轴承分析模型,考虑了滚动体的离心力、陀螺力矩以及弹性流体动 力润滑作用.建立了非线性方程组未知量的约束条件,解决了Newton-Raphson法求解方程易出现不收敛的问 题.以两个角接触球轴承为算例,分析了轴承内、外圈分别旋转对应力的影响,以及转速和径向载荷对轴承工作 性能的影响.结果表明,轴承外圈转时接触应力随转速增大而显著增大,内圈转时则无明显变化;低转速下外圈 转的接触应力明显低于内圈转;较小的径向栽荷下外圈转的接触应力高于内圈转;轴承微区滑动速度受转速和径向栽荷的影响显著.通过Jones球轴承拟静力学分析程序计算,验证了模型和程序的有效性.     

关于拟线性椭圆型方程在广义Sobolev空间中解的正则性

本文讨论了下面方程Ｄｉｒｉｃｈｌｅｔ问题在广义Ｓｏｂｏｌｅｖ空间中解中的正则性－ｄ／ｄｘ１ａ１（ｘ，ｕ，Ｄｕ）＋ａ（ｘ，ｕ，Ｄｕ）＝０，ｘ∈Ω其中Ω∈Ｒ是有界区域，证明了上述问题在Ｗ（Ω）和Ｗ（Ω）在存在有界广义解。     

一类非线性椭圆型方程非平凡解和分歧点的存在性

本文给出了所谓集中到列紧原理的另一形式,并利用它讨论了一类非线性椭圆型方程广义银和分歧点的顾在性。     

平面Poisson问题的间接变量等价边界积分方程

使用变分法确立平面Laplace外问题的确切形式；在此基础上，为了避免使用无限区域上的Green公式(不成立)，通过构造紧支集上的无穷可微函数，获得满足非齐次控制微分方程的特解的重要结果。使用非解析开拓数学方法，导出满足控制微分方程的位势函数集的等价边界积分表示式。获得Poisson问题的间接变量等价边界积分方程；通过一些例子，详细讨论了传统的外问题和边界积分方程的形式，表明它们不具有普遍适用性。