多段翼型位流的准确解

本文给出了求解多段翼型位流准确解的一种方法。采用求解拉普拉斯方程把多段翼型绕流区域变换成一个有限矩形域,然后把此矩形域变换成多个圆,从而通过求绕多个圆的准确位流解来求得绕多段翼型位流的准确解。本文提供的方法能广泛应用到任意外形的多段翼型位流计算。采用这种方法计算所得的位流压力分布与用其它方法求得的压力分布非常吻合。     

振荡机翼跨音速非定常粘性绕流的数值计算

以非定常Ｎ－Ｓ方程为主管方程，采用ＬＵ－ＮＮＤ混合差分格式，Ｃ和Ｃ－Ｈ型贴体运动网格，Ｂ－Ｌ双层代数紊流模型，求解绕振荡翼型和三维机翼的跨音速非定常粘性流场，分别计算了ＮＡＣＡ００１２翼型和Ｍ６机翼作俯仰振荡时跨音速非定常粘性绕流流场。研究了非定常绕流的气动特性，部分计算结果和风洞实验值作了比较。     

带副翼的翼身组合体绕流的Euler和N-S方程解

将对接分区网格与分区求解算法结合,有效地求解了带副翼偏转的翼身组合体绕流的N S方程.数值方法中选用VanLeer分裂格式离散无粘通量项,采用中心差分格式来离散粘性通量项.分区交界面采用了一种满足通量守恒的内边界耦合条件.数值算例表明该方法是求解带操纵面偏转的翼身组合体绕流的有效方法.     

BGK方法在非结构网格上的应用

采用旋转局部坐标的方法,发展了一种针对非结构网格的BGK计算方法.该方法属于有限体积法,大致分为两个步骤:①空间离散:②通量计算及时间推进.在第①步中,采用基于最小二乘法的高阶ENO格式来获得宏观物理量的高阶导数;在第②步中,采用旋转坐标轴的方法来计算非结构网格单元各边的通量.并得出了后台阶绕流(Backward Facing Step)及翼型绕流(Flow Over an Airfoil)两个算例的计算结果.     

Euler及N-S方程的二维四边形/叉树混合网格自适应算法

描述了一种基于直角叉树网格的Euler和N-S方程自适应算法.由于考虑了粘性的作用,提出并使用了四边形叉树混合网格的方法,在几何表面附近生成贴体的四边形网格,外流场使用直角叉树网格.采用中心有限体积法,对Euler及N-S方程进行数值求解,对N-S方程的计算中加入了B-L代数湍流模型.在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离等特性的捕捉和分辨能力.采用上述方法,数值分析了单段和多段翼型的绕流问题.     

隐式格式求解拟压缩性非定常不可压Navier-Stokes方程

采用Rogers发展的双时间步拟压缩方法,数值求解不可压非定常问题.数值通量分别采用三阶精度Roe格式和二阶精度Harten-Yee的TVD格式离散.为了加快收敛,提高求解效率,试验了几种隐式格式(ADI-LU,LGS,LU-SGS).针对经典的低雷诺数(Re=200)圆柱绕流问题,比较了不同隐式方法的计算结果和求解效率,以及两种数值离散格式计算结果的异同.最后采用Roe格式数值求解了两种典型的低速非定常流动问题:绕转动圆柱(ω=1)低雷诺数流动;NACA0015翼型等速拉起数值模拟.     

多块网格技术及其在操纵面绕流计算中的应用

采用中心有限体积法求解Navier Stokes方程组,利用点对点搭接及面搭接多块网格技术成功地进行了多段翼型和存在外形间断(剪刀差)的带操纵面三维机翼的网格生成及粘性绕流流场数值计算,使用重叠面积加权插值及Ramshaw算法实现了面搭接的多块网格之间的信息传递并保证了通量守恒,并将流场计算结果与实验数据进行了对比.     

不规则对接网格的分区求解计算

采用一种保持通量守恒的不规则对接网格分区求解中交界面耦合条件的计算方法, 结合有限体积法求解了Euler 方程, 无粘通量取用Van Leer 分裂格式, 构造了一种限制器以实现格式的二阶精度和TVD 性质, 并给出数值算例。     

二维翼型跨音速绕流的数值模拟

本文结合高精度TVD格式的数值通量和时间进展多步法给出了一种求解定常流问题的数值方法。同时给出一些特殊处理来加快数值解的收敛速度。本文用以上方法计算了翼型跨音速绕流问题,结果表明此方法具有分辨率高,收敛速度较快之性质。     

用涡方法数值模拟运动物体非定常流场

采用涡方法对静止圆柱绕流、简谐振动圆柱、静止NACA0012翼型绕流和运动扑翼的非定常流场进行了数值模拟和数值验证,对比验证结果较好.控制方程采用二维不可压缩N-S方程,根据算子分裂法,对流和扩散分别采用不同的时间步长.涡方法基于Lagrange坐标系的特点使其方便模拟运动物体的非定常流场.本文的工作为今后的翼型动态失...     

���·��̵ķ��Ѹ�ʽ�����ڽ�յ������廷��ײ�о��е�Ӧ��

从Lagrange观点出发,采用分裂格式法求解一维辐射流体力学方程组中的三温方程,用于紧凑等离子体环与靶碰撞产生的X射线辐射的数值模拟研究。提出了求解电子、离子、辐射场三温相脱离的能量方程的隐式差分格式,介绍了能量交换项与压力做功及热传导项分开计算的分裂格式方法,数值模拟得到了其对温度变化的影响。     

A conservative numerical method for the Cahn–Hilliard equation in complex domains

We propose an efficient finite difference scheme for solving the Cahn–Hilliard equation with a variable mobility in complex domains. Our method employs a type of unconditionally gradient stable splitting discretization. We also extend the scheme to compute the Cahn–Hilliard equation in arbitrarily shaped domains. We prove the mass conservation property of the proposed discrete scheme for complex domains. The resulting discretized equations are solved using a multigrid method. Numerical simulations are presented to demonstrate that the proposed scheme can deal with complex geometries robustly. Furthermore, the multigrid efficiency is retained even if the embedded domain is present.     

An unconditionally stable rotational velocity-correction scheme for incompressible flows

We present an unconditionally stable splitting scheme for incompressible Navier–Stokes equations based on the rotational velocity-correction formulation. The main advantages of the scheme are: (i) it allows the use of time step sizes considerably larger than the widely-used semi-implicit type schemes: the time step size is only constrained by accuracy; (ii) it does not require the velocity and pressure approximation spaces to satisfy the usual inf–sup condition: in particular, the equal-order finite element/spectral element approximation spaces can be used; (iii) it only requires solving a pressure Poisson equation and a linear convection–diffusion equation at each time step. Numerical tests indicate that the computational cost of the new scheme for each time step, under identical time step sizes, is even less expensive than the semi-implicit scheme with low element orders. Therefore, the total computational cost of the new scheme can be significantly less than the usual semi-implicit scheme.     

Woods-Saxon势下的赝自旋对称性

通过求解具有Woods-Saxon径向标量势与矢量势的Dirac方程, 分别分析了原子核中赝自旋双重态能级劈裂和波函数下分量劈裂随着Woods-Saxon势参数的变化关系, 其中Wood-Saxon势中的参数a以及R对于赝自旋能级劈裂和波函数劈裂的影响最大, 而势阱深度V₀只产生微小的影响. 随着平均场的变化, 赝自旋双重态会出现能级反转现象. 具有较大n或l的赝自旋双重态能级劈裂对于参数的变化更敏感. 赝自旋波函数下分量劈裂随着参数a和R的增大向核表面扩散, 并且在核表面附近达到最大. 赝自旋波函数劈裂随着参数a的增大而增大, 但是随着参数R或V₀的增大它却在减小的. 由于在特定的同位素链中, 参数R和V₀都与核子数有关系, 而参数a又是描述核的表面性质的, 所以以这些参数为变量对于赝自旋劈裂的研究是有意义的, 研究的结果至少在定性上可以应用到大部分原子核中.     

Splitting a qudit state via Greenberger-Horne-Zeilinger states of qubits

We present a tripartite quantum information splitting scheme which splits a qutrit state via two GHZ states. The scheme is then generalized to splitting a qudit state among any number of receivers. We show that this scheme is also applicable to splitting any multi-qudit entangled states.     

求解多维欧拉方程的二阶非结构网格混合旋转Riemann求解器

将基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法推广到非结构网格,采用基于网格中心的有限体积法,梯度的计算采用基于节点的方法引入更多的控制体模板,限制器的构造采用与非结构化网格相适应的形式.在求解Riemann问题时,沿具有一定物理意义的两个迎风方向,即控制体界面两侧速度差矢量方向及与之正交的方向.能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或"红斑"现象.为减小计算量,采用HLL和Roe FDS混合旋转格式.     

Explicit Multi-Symplectic Splitting Methods for the Nonlinear Dirac Equation

In this paper, we propose two new explicit multi-symplectic splitting methods for the nonlinear Dirac (NLD) equation. Based on its multi-symplectic formulation, the NLD equation is split into one linear multi-symplectic system and one nonlinear infinite Hamiltonian system. Then multi-symplectic Fourier pseudospectral method and multi-symplectic Preissmann scheme are employed to discretize the linear subproblem, respectively. And the nonlinear subsystem is solved by a symplectic scheme. Finally, a composition method is applied to obtain the final schemes for the NLD equation. We find that the two proposed schemes preserve the total symplecticity and can be solved explicitly. Numerical experiments are presented to show the effectiveness of the proposed methods.     

Algebraic Splitting for Incompressible Navier–Stokes Equations

Fully discretized incompressible Navier–Stokes equations are solved by splitting the algebraic system with an approximate factorization. This splitting affects the temporal convergence order of velocity and pressure. The splitting error is proportional to the pressure variable, and a simple analysis shows that the original convergence order of the time-integration scheme can be retained by solving for incremental pressure. The combination of splitting and incremental pressure is shown to be equivalent to an error-correcting method using the full pressure. In numerical experiments employing a third-order time-integration scheme and various orders for the pressure increment, the splitting error is shown to control the convergence order, and the full order of the scheme is recaptured for both velocity and pressure. The difference between perturbing the momentum or the continuity equation is also explored.     

A generalized tripartite scheme for splitting an arbitrary 2-qubit state with three 2-qubit partially entangled states and the Kraus measurement

We put forward a generalized tripartite scheme for splitting an arbitrary 2-qubit pure state with three 2-qubit non-maximally en-tangled states as quantum channels.The scheme for the first time incorporates the Kraus measurement into quantum information splitting scheme.In contrast to the similar scheme using the same quantum channels and the ancilla-entangled measurement,our scheme is superior in terms of operation and complexity,success probability,resource consumption and effciency.     

一种基于AUSM思想的通量分裂方法

根据对流迎风分裂(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格式在计算运动接触间断时,速度存在很大振荡,而新格式不存在振荡.