压力修正算法的并行实现方法研究

基于分块非结构化网格上的SIMPLE算法、区域分解算法及MPI并行编程方法,给出了一种并行计算流体力学实施方案。提出了"串行粗粒度,并行细粒度"的区域分解及网格划分方法;依据内边界网格界面上的通量守恒原则,推导出了子区域间内界面上的变量传递关系式。在魔方计算机上使用500处理器核完成了两类流动问题计算。研究表明,并行计算结果、串行计算结果与基准解吻合很好,通信时间、cache命中率对并行加速比有显著影响。     

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

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

非结构网格在平面叶栅内湍流流动数值模拟中的应用

1控制方程和离散方程对二维非定常可压缩粘性流体，雷诺时均Navierstokes方程可写成如下积分形式：式中，0为计算区域；＊D为区域边界；Q为守恒变量；F为通量，包括无粘通量F’和粘性通量F”。计算域采用三角形单元离散，三角形的边作为控制体边界，物理量置于单元中心。对每个三角形控制体，方程（1）可写成：·其中，Q；表示单元的守恒变量值；aC和Vi分别表示控制体的边界和面积，为计算式（2中右端值，将单元所有边界通量相加，最后可得如下方程：式中k（i）为单元边界，只，j为通过单元界面的通量，凸马为界面长度。2无粘通量本文采…     

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

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

辛格式在强激光场一维模型问题计算中的应用

 将辛算法推广到复辛空间，指出了辛算法保定态Schr-dinger方程的Wronskian守恒。将辛算法应用于强场一维模型的计算中，并与Runge-Kutta法作了比较。结果显示，辛算法保持定态Schr-dinger方程的Wronskian守恒，适合于在充分 远空间上计算线性无关解，是计算强激光场一维模型的合理的数值方法。     

双曲型守恒律的一种高精度TVD差分格式

构造了一维双曲型守恒律方程的一个高精度高分辨率的守恒型TVD差分格式.其主要思想是:首先将计算区域划分为互不重叠的小单元,且每个小单元再根据希望的精度阶数分为细小单元;其次,根据流动方向将通量分裂为正、负通量,并通过小单元上的高阶插值逼近得到了细小单元边界上的正、负数值通量,为避免由高阶插值产生的数值振荡,进一步根据流向对其进行TVD校正;再利用高阶Runge KuttaTVD离散方法对时间进行离散,得到了高阶全离散方法.进一步推广到一维方程组情形.最后对一维欧拉方程组计算了几个算例.     

LED显示屏色域边界的快速计算

为了快速计算出LED显示屏的色域边界,提出一种迭代求解映射线与色域边界交点进而拟合每个色相面内色域边界的算法。首先根据LED显示的色度学原理给出了颜色是否在LED显示屏色域内的判断方法。然后在CIELAB均匀颜色空间中,分析了迭代算法的原理,确定了迭代计算初始值,并进行了计算精度和速度分析。最后以D65为参照,在等色相面内进行了LED显示屏色域边界的拟合实验及色差分析。实验结果表明:11次迭代运算后,拟合色域边界已非常光滑,每条映射线上的色域边界点真实值和计算值最大色差仅为0.23。与插值类计算方法相比,最大色差值和计算速度分别降低和提高了一个数量级。本文算法在计算精度和速度上能够满足对LED显示屏进行色域分析和与其他显示设备进行色域映射的要求。     

强场一维模型问题的保Wronskian守恒算法

提出辛算法是保Wronskian守恒的算法,把辛算法应用于强激光场一维模型的计算中,结果显示,Wronskian保持不变,与理论分析一致。     

多介质流体力学计算的守恒型高分辨率格式

应用Lagrange坐标系下的守恒型格式计算多介质流体力学问题,在物质交界面附近采用一阶格式的数值通量,而在其余部分采用高分辨率格式的数值通量,不仅保持了高分辨率的良好性质,而且消除了一般的守恒型格式在界面附近所产生的震荡.     

ALE中的重映算法

介绍了二维非结构网格上的守恒重映算法，重点是基于SFB／DC思想的通量重映算法。用统一的公式表示不同的单元量重映算法，包括原始的贡献网格法、Barth—Jespersen方法、最小二乘法，不同算法间的区别体现为梯度求法的差异。对于交错网格上速度的重映，介绍了SALE和HIS算法。此外，为保证重映算法的有界性，引入了修补方法。     

湍流模型在内流流场数值模拟中的应用

本文将Spalart-Allmaras湍流模型方程推导为守恒形式,并应用于内流湍流流场计算.采用此-方程湍流模型分别对扩压器流场以及NASA 67转子流场进行了数值模拟.结果表明常用在外流计算中的Spalart-Allmaras湍流模型在内流计算中仍具有较高的精度、数值稳定性和效率.     

On the HLLC Riemann solver for interface interaction in compressible multi-fluid flow

In this work, the HLLC Riemann solver, which is much more robust, simpler and faster than iterative Riemann solvers, is extended to obtain interface conditions in sharp-interface methods for compressible multi-fluid flows. For interactions with general equations of state and material interfaces, a new generalized Roe average is proposed. For single-phase interactions, this new Roe average does not introduce artificial states and satisfies the U-property exactly. For interactions at material interfaces, the U-property is satisfied by introducing ghost states for the internal energy. A number of numerical tests suggest that the proposed Riemann solver is suitable for general equations of state and has an accuracy comparable to iterative Riemann solvers, while being significantly more robust and efficient.     

湿蒸汽两相流动问题的耦合求解

控制湿蒸汽两相流动(汽相和液相)的方程个数多,其数值模拟通常采用将两相分开、两相间用参数传递的方法来处理。本文首次尝试将湿蒸汽两相流动汽相和液相的8个控制方程统一求解,推导了这一方程组的特征值及对应的特征向量,结合控制体积方法,采用TVD格式和时间推进法数值求解这一耦合方程组,研发了相应计算软件。对二维喷管中的自发凝结和非均质凝结算例进行的数值模拟和分析,表明了本文方法和计算软件的有效性和可靠性。     

Modeling Three-Dimensional Multiphase Flow Using a Level Contour Reconstruction Method for Front Tracking without Connectivity

Three-dimensional multiphase flow and flow with phase change are simulated using a simplified method of tracking and reconstructing the phase interface. The new level contour reconstruction technique presented here enables front tracking methods to naturally, automatically, and robustly model the merging and breakup of interfaces in three-dimensional flows. The method is designed so that the phase surface is treated as a collection of physically linked but not logically connected surface elements. Eliminating the need to bookkeep logical connections between neighboring surface elements greatly simplifies the Lagrangian tracking of interfaces, particularly for 3D flows exhibiting topology change. The motivation for this new method is the modeling of complex three-dimensional boiling flows where repeated merging and breakup are inherent features of the interface dynamics. Results of 3D film boiling simulations with multiple interacting bubbles are presented. The capabilities of the new interface reconstruction method are also tested in a variety of two-phase flows without phase change. Three-dimensional simulations of bubble merging and droplet collision, coalescence, and breakup demonstrate the new method's ability to easily handle topology change by film rupture or filamentary breakup. Validation tests are conducted for drop oscillation and bubble rise. The susceptibility of the numerical method to parasitic currents is also thoroughly assessed.     

多排叶片内部流动以及叶片排与排气系统相互作用研究

本文采用高分辨率格式和多块多网格方法求解雷诺平均Navier-Stokes方程,数值模拟多排叶片内的三维粘性流动;使用数值激盘模拟叶片排的存在,研究叶片排与非轴对称排气部件之间的相互作用以及复杂的内部流动。文中描述了数值方法,给出了NASA透平导叶和单级透平内部流场的数值结果及其与相应实验结果的对比,也给出了多级透平内部流场的数值结果,以及透平与排气部件之间的耦合流场的结果。     

Flows with suspended and floating particles

The evolution of the configuration of a set of particles dispersed in a flowing liquid is crucial in many applications such as sedimentation, slurry transport, rheology and structured arrays of micro- and nano-particles. Direct simulation based on what is called fictitious domain method coupled with finite element method has been used to study particulate flows and sedimentation process. Here we extend the previously proposed formulations to naturally include buoyancy force and the capillary driven attraction or repulsion of particles located at fluid interfaces. The set of differential equations is discretized using a fully implicit-fully coupled fictitious domain/finite element approach, avoiding numerical instabilities that may arise from explicit integration. The proposed formulation and implementation are validated by comparing the predictions of simple 2D flows to available numerical or analytical solutions. The method is then used to analyze the flotation of 2D particles and capillary driven aggregation at fluid interfaces.     

Numerical simulation of drop impact on a liquid–liquid interface with a multiple marker front-capturing method

The gravity-driven motion of a droplet impacting on a liquid–liquid interface is studied. The full Navier–Stokes equations are solved on a fixed, uniform grid using a finite difference/front-capturing method. For the representation of fluid–fluid interfaces, a coupled Level-Set/Volume-Of-Fluid method [M. Sussman, E.G. Puckett, A coupled Level-Set and Volume-of-Fluid method for computing 3D and axisymmetric incompressible two-phase flows, J. Comp. Phys. 162 (2000) 301–337] is used, in which we introduce the novel approach of describing separate interfaces with different marker functions. As a consequence, we prevent numerical coalescence of the droplet and the liquid–liquid interface without excessive (local) grid refinement. To validate our method, numerical simulations of the drop impact event are compared with experiments [Z. Mohamed-Kassim, E.K. Longmire, Drop impact on a liquid–liquid interface, Phys. Fluids 15 (2003) 3263–3273]. Furthermore, a comparison is made with the numerical results of [A. Esmaeeli, G. Tryggvason, Direct numerical simulations of bubbly flows. Part 2. Moderate Reynolds number arrays, J. Fluid Mech. 385 (1999) 325–358] for an array of rising bubbles. The investigation shows that the multiple marker approach successfully prevents numerical coalescence of interfaces and adequately captures the effect of surface tension.     

An all-speed relaxation scheme for interface flows with surface tension

We consider interface flows where compressibility and capillary forces (surface tension) are significant. These flows are described by a non-conservative, unconditionally hyperbolic multiphase model. The numerical approximation is based on finite-volume method for unstructured grids. At the discrete level, the surface tension is approximated by a volume force (CSF formulation). The interface physical properties are recovered by designing an appropriate linearized Riemann solver (Relaxation scheme) that prevents spurious oscillations near material interfaces. For low-speed flows, a preconditioning linearization is proposed and the low Mach asymptotic is formally recovered. Numerical computations, for a bubble equilibrium, converge to the required Laplace law and the dynamic of a drop, falling under gravity, is in agreement with experimental observations.     

贴体坐标系下模型加力室的大涡模拟

本文利用贴体网格对带V形槽稳定器模型加力燃烧室素流化学反应流流动进行大涡模拟的研究。采用区域法生成模型加力燃烧室的二维贴体网格,并采用多区域耦合法进行区域之间的数据传递,求解加力室整体流场。采用k方程亚网格尺度模型和亚网格EBU燃烧模型分别估算其亚网格紊流粘性和化学反应速率,用热通量辐射模型估算辐射通量,并用交错网格下SIMPLE算法和混合差分格式求解离散方程,壁面函数处理固壁边界条件。计算结果显示了稳定器后面的回流区气流结构,所得的热态流场模拟结果与实验比较吻合,表明采用贴体网格对模型加力燃烧室进行大涡模拟能真实反映流体流动及燃烧过程。