基于OpenFOAM的投影法磁流体求解器开发与验证

在开源计算流体力学C++工具包OpenFOAM环境下开发了低磁雷诺数条件下的磁流体求解器,并进行了验证。采用投影算法求解动量方程和压力泊松方程;采用非结构网格同位相容守恒算法求解电势泊松方程、感应电流和洛伦兹力;采用边界耦合方法求解流固耦合电势场。通过对均匀磁场下导电方管和导电圆管内的完全发展磁流体层流的数值模拟和解析解的对比,对求解器进行了验证。进一步对非均匀强磁场作用下导电方管和导电圆管内完全发展磁流体层流进行了数值模拟,并与ALEX实验结果进行了比较。数值解和实验结果吻合良好。所开发的求解器可用于复杂结构强磁场作用下磁流体的数值模拟研究。     

基于OpenFOAM的投影法磁流体求解器开发与验证

在开源计算流体力学C++工具包OpenFOAM环境下开发了低磁雷诺数条件下的磁流体求解器,并进行了验证。采用投影算法求解动量方程和压力泊松方程;采用非结构网格同位相容守恒算法求解电势泊松方程、感应电流和洛伦兹力;采用边界耦合方法求解流固耦合电势场。通过对均匀磁场下导电方管和导电圆管内的完全发展磁流体层流的数值模拟和解析解的对比,对求解器进行了验证。进一步对非均匀强磁场作用下导电方管和导电圆管内完全发展磁流体层流进行了数值模拟,并与ALEX实验结果进行了比较。数值解和实验结果吻合良好。所开发的求解器可用于复杂结构强磁场作用下磁流体的数值模拟研究。     

基于浸入边界-多松弛时间格子玻尔兹曼通量求解法的流固耦合算法研究

针对流固耦合问题,发展了基于浸入边界-多松弛时间格子玻尔兹曼通量求解法(immersed boundary method multi-relaxation-time lattice Boltzmann flux solver,IB-MRT-LBFS)的弱耦合算法.依据多尺度Chapman-Enskog展开,建立不可压宏观方程状态变量和通量与格子玻尔兹曼方程中粒子密度分布函数之间的关系;采用强制浸入边界法处理流固界面使固壁表面满足无滑移边界条件,根据修正的速度求解动量方程力源项;结构运动方程采用四阶龙格-库塔法求解.格子模型与浸入边界法的引入使流固耦合计算可以在笛卡尔网格下进行,无需生成贴体网格及运用动网格技术,简化了计算过程.数值模拟了单圆柱横向涡激振动、单圆柱及串列双圆柱双自由度涡激振动问题.结果表明,IB-MRT-LBFS能够准确预测圆柱涡激振动的锁定区间、振动响应、受力情况以及捕捉尾流场结构形态,验证了该算法在求解流固耦合问题的有效性和可行性.     

流动/运动松耦合与紧耦合计算方法及稳定性分析

为了更加精确地模拟流动/运动耦合问题, 建立了耦合动态混合网格生成非定常流场计算和六自由度运动方程求解的一体化计算方法, 并在统一框架内同时实现了松耦合与紧耦合方法.通过圆柱涡致自激振荡(vortex induced vibration, VIV)的模拟, 对不同时间精度的松耦合和紧耦合算法的优劣及适用范围进行了评估和分析; 通过引入附加质量的概念, 对耦合算法的稳定性进行了理论分析.研究表明:在流体的密度与物体的密度接近时, 松耦合方法是不稳定的, 必须采用紧耦合方法.最后利用耦合算法对二维鱼体的自主游动和钝锥三自由度自由飞过程进行了数值模拟, 证实了理论分析的结论.      

基于CFD/CSD的舵面控制气动弹性数值模拟方法

在动态网格上通过耦合求解流动控制方程和结构动力学方程, 发展了一种舵面控制下飞行器运动响应过程中气动弹性数值模拟研究方法.流动控制方程采用N-S方程, 结构动力学采用线性模态叠加方法, 其中流动控制方程空间离散采用基于非结构网格的有限体积方法, 对流通量采用计算HLLC格式, 非定常时间离散采用基于LU-SGS的双时间步长方法.模拟中, 气动运动和结构变形在双时间步长方法推进过程中采用改进松耦合方法, 气动网格与结构网格之间信息交换采用无限平板样条法实现, 飞行器的运动和变形采用基于重叠网格和Delaunay图映射变形网格相结合的方法进行处理.采用多个考核算例对发展的数值方法进行考核验证, 结果表明该方法可以高效精确模拟舵面开环控制下飞行器运动响应过程中的气动弹性特性.      

基于有限体积方法的复合材料固化过程模拟

传热与固化过程的耦合作用在复合材料成型工艺中具有重要影响。本文基于非结构化网格有限体积法提出一种解决结构形式复杂的复合材料成型过程中各向异性非稳态导热与固化反应动力学耦合问题的算法。各向异性复合材料的导热控制方程空间项采用非结构化网格有限体积法离散,时间项采用欧拉隐格式求解。固化动力学方程采用了欧拉显格式计算,并通过固化反应速率计算出能量方程的热源,同时根据温度分布及其变化更新固化反应速率,实现导热与固化反应的双向耦合。最后,通过文献中的几个经典算例验证了本文方法的正确性和可靠性,同时也表明本文方法可以处理复杂结构形式的复合材料成型过程的传热与化学反应的耦合模拟问题。     

耦合辐射的三维热化学非平衡流数值模拟

发展耦合辐射的三维热化学非平衡流场计算方法,可用于非结构网格.采用Jameson有限体积法求解耦合辐射源项的三维N-S方程.辐射源项通过求解辐射输运方程(Radiative Transport Equation RTE)获得.在空间方向上离散后,采用有限体积法求解辐射输运方程.化学模型包含11个组元,20个化学反应.采用该数值方法计算MUSES-C模型在速度为11.6 km·s^-1时的绕流流场及前驻点处的辐射热流密度.并通过对比,分析热辐射对流场的影响.     

谐振管内非线性驻波的有限体积数值算法

为了扩展谐振管内非线性驻波在工程中的应用, 以及克服现有数值计算方法仅局限于求解直圆柱形和指数形谐振管内非线性驻波的问题. 根据变截面的非稳态可压缩热黏性流体Navier-Stokes方程和空间守恒方程, 并基于求解压力速度耦合方程的半隐式算法和交错网格技术, 构建一种能够计算任意形状轴对称谐振管受活塞驱动时内部非线性驻波的有限体积算法. 分别对圆柱形、指数形和圆锥形谐振管内的非线性驻波进行仿真计算. 通过与现有试验结果以及数值仿真结果的对比, 验证了该方法的正确性.并获得除驻波声压之外的另外一些新的物理结果, 包括速度、密度、温度的瞬时变化.在直圆柱形谐振管内产生冲击声压波, 速度波形中出现钉状结构.而在指数形和圆锥形谐振管内产生高声压幅值的驻波, 没有出现冲击波, 速度波形中均未发现钉状结构. 计算结果表明谐振管内非线性驻波的物理属性与谐振管形状之间有密切关系.     

基于非结构网格的不可压流动耦合求解器

采用非结构化网格有限容积法求解了不可压N-S方程组,对流项采用GAMMA格式,扩散项采用二阶中心差分格式建立离散方程,用SOAR算法处理压力与速度的耦合关系,得到了一种求解不可压N-S方程的非结构网格耦合求解器。通过方腔顶盖驱动流、后台阶绕流以及方腔自然对流等几个典型的算例,考察了求解器的计算精度及收敛特性,并与SIMPLE算法进行了比较,结果表明该求解器是有效可行的。     

三维机翼跨音速颤振的分布式并行计算

以多台微机组成分布式计算系统,利用无限插值方法(TFI)生成三维多块贴体运动网格,以Navier Stokes方程为控制方程,求解三维机翼的跨音速非定常气动力,并与颤振方程耦合迭代计算,求解机翼广义位移响应的时间历程,根据广义位移的时间历程的衰减、等幅和发散振荡等情况确定机翼跨音速颤振临界条件.通过算例验证,计算结果与实验结果和理论分析相吻合.     

A coupled finite volume solver for the solution of incompressible flows on unstructured grids

This paper reports on a newly developed fully coupled pressure-based algorithm for the solution of laminar incompressible flow problems on collocated unstructured grids. The implicit pressure-velocity coupling is accomplished by deriving a pressure equation in a procedure similar to a segregated SIMPLE algorithm using the Rhie–Chow interpolation technique and assembling the coefficients of the momentum and continuity equations into one diagonally dominant matrix. The extended systems of continuity and momentum equations are solved simultaneously and their convergence is accelerated by using an algebraic multigrid solver. The performance of the coupled approach as compared to the segregated approach, exemplified by SIMPLE, is tested by solving five laminar flow problems using both methodologies and comparing their computational costs. Results indicate that the number of iterations needed by the coupled solver for the solution to converge to a desired level on both structured and unstructured meshes is grid independent. For relatively coarse meshes, the CPU time required by the coupled solver on structured grid is lower than the CPU time required on unstructured grid. On dense meshes however, this is no longer true. For low and moderate values of the grid aspect ratio, the number of iterations required by the coupled solver remains unchanged, while the computational cost slightly increases. For structured and unstructured grid systems, the required number of iterations is almost independent of the grid size at any value of the grid expansion ratio. Recorded CPU time values show that the coupled approach substantially reduces the computational cost as compared to the segregated approach with the reduction rate increasing as the grid size increases.     

应用笛卡尔非结构切割网格进行外挂物投放的数值模拟

描述了一种新的网格生成技术,即笛卡尔非结构切割网格技术,采用叉树数据结构,完成了几种单段和多段翼型以及三维机翼的网格生成.应用中心有限体积法,对其绕流问题进行Euler方程数值模拟,并将计算结果与实验数据进行对比.在机翼绕流数值模拟的基础上,求解出机翼带外挂物的分离投放的流场计算问题.     

Mesh refinement algorithms in an unstructured solver for multiphase flow simulation using discrete particles

This study developed spray-adaptive mesh refinement algorithms with directional sensitivity in an unstructured solver to improve spray simulation for internal combustion engine application. Inadequate spatial resolution is often found to cause inaccuracies in spray simulation using the Lagrangian–Eulerian approach due to the over-estimated diffusion and inappropriate liquid–gas phase coupling. Dynamic mesh refinement algorithms adaptive to fuel sprays and vapor gradients were developed in order to increase the grid resolution in the spray region to improve simulation accuracy. The local refinement introduced the coarse-fine face interface that requires advanced numerical schemes for flux calculation and grid rezoning with moving boundaries. To resolve the issue in flux calculation, this work implemented the refinement/coarsening algorithms into a collocated solver to avoid tedious interpolations in solving the momentum equations. A pressure correction method was applied to address unphysical pressure oscillations due to the collocation of pressure and velocity. An edge-based algorithm was used to evaluate the edge-centered quantities in order to account for the contributions from all the cells around an edge at the coarse-fine interface. A quasi-second-order upwind scheme with strong monotonicity was also modified to accommodate the coarse-fine interface for convective fluxes. To resolve the issue related to grid rezoning, rezoning was applied to the initial baseline mesh only and the new locations of the refined grids were obtained by interpolating the updated baseline mesh. The time step constraints were also re-evaluated to account for the change resulting from mesh refinement. The present refinement algorithm was used in simulating fuel sprays in an engine combustion chamber. It was found that the present approach could produce the same level of results as those using the uniformly fine mesh with substantially reduced computer time. Results also showed that this approach could alleviate the artifacts related to the Lagrangian discrete modeling of spray drops due to insufficient spatial resolution.     

液舱横荡非结构网格高阶格式数值模拟

提出一种液舱横荡数值模拟的方法,将气液两相交界面视为物理间断,通过高阶精度离散格式捕捉间断.根据NVD(Normalized Variable Diagram)实现非结构化网格上高精度离散格式,建立固定网格上自由表面运动模拟方法.在开发的非结构网格有限体积法求解器GTEA(General Transport Equation Analyzer)基础上,实现上述方法.首先对经典的溃坝过程进行模拟,并与文献结果对比验证方法和程序的可信度.对二维矩形液舱在不同激振频率时的横荡进行数值计算,并与实验以及商业软件CFX计算结果进行比较.结果表明方法和软件可以模拟自由面的翻卷、破碎运动现象,对距自由面较深点处流体载荷的计算结果与实验值符合较好,与商业软件CFX相比,在相同计算网格下,算法可以更好的计算次峰值,验证方法正确可行.     

CGNS API和FVM在非结构混合网格计算中的应用

用CGNS API（CFD General Notation System Application Programming Interface）作为非结构混合网格求解器的前处理和后处理，用FVM（Finite Volume Method）作为偏微分方程求解方法．在前处理过程中，用hash表法对内部网格面和边界网格面进行编号，并计算出相应的几何信息，以满足FVM求解器的需要．从FVM求解器计算出来的各种场信息可以写入原来的CGNS文件，该文件可以被许多专业商业后处理软件（如Tecplot，Fluent，CFX等）读取和进行可视化；对于求解器，用基于网格中心的FVM及SIMPLEC（Semi Implicit Methodfor Pressure Linked Equation Consistent）方法求解压力速度耦合．最后给出两个说明算例．     

Spectral difference method for compressible flow on unstructured grids with mixed elements

This paper presents the development of a 2D solver for inviscid and viscous compressible flows using the spectral difference (SD) method for unstructured grids with mixed elements. A mixed quadrilateral and triangular grid is first refined using one-level h-refinement to generate a quadrilateral grid while keeping the curvature of boundary edges. The SD method designed for quadrilateral meshes can subsequently be applied for the refined unstructured grid. Results obtained with the SD method for both inviscid and viscous compressible flows compare well with analytical solutions and other published results.     

不可压缩流基于块预处理的并行有限元计算

发展了一个模拟非定常不可压缩粘性流的并行有限元求解器,时间离散使用具有二阶精度的隐式中点格式,基于三维非结构四面体网格剖分,使用高阶混合有限元离散速度场（P2）和压力场（P1）.全离散格式产生的代数方程组是大型、稀疏、非对称和病态的,基于修正的压力对流扩散预处理（PCD）和精心设计的子问题迭代执行策略,采用预处理的GMRES迭代法来高效求解线性方程组.利用相同的子问题迭代策略,同时给出基于最小二乘交换子（LSC）预处理的并行效率对比.大量数值算例验证了算法的精度、可扩展性和可靠性.三维驱动方腔流模拟结果（Re=3200.0）清晰地显示了方腔流中主涡（PE）、下游二次涡（DSE）、上游二次涡（USE）、侧壁涡（EWV）和TGL涡的存在.     

A sweeping preconditioner for time-harmonic Maxwell’s equations with finite elements

This paper is concerned with preconditioning the stiffness matrix resulting from finite element discretizations of Maxwell’s equations in the high frequency regime. The moving PML sweeping preconditioner, first introduced for the Helmholtz equation on a Cartesian finite difference grid, is generalized to an unstructured mesh with finite elements. The method dramatically reduces the number of GMRES iterations necessary for convergence, resulting in an almost linear complexity solver. Numerical examples including electromagnetic cloaking simulations are presented to demonstrate the efficiency of the proposed method.     

A balanced force refined level set grid method for two-phase flows on unstructured flow solver grids

A balanced force refined level set grid method for two-phase flows on structured and unstructured flow solver grids is presented. To accurately track the phase interface location, an auxiliary, high-resolution equidistant Cartesian grid is introduced. In conjunction with a dual-layer narrow band approach, this refined level set grid method allows for parallel, efficient grid convergence and error estimation studies of the interface tracking method. The Navier–Stokes equations are solved on an unstructured flow solver grid with a novel balanced force algorithm for level set methods based on the recently proposed method by Francois et al. [M.M. Francois, S.J. Cummins, E.D. Dendy, D.B. Kothe, J.M. Sicilian, M.W. Williams, A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework, J. Comput. Phys. 213 (2006) 141–173] for volume of fluid methods on structured grids. To minimize spurious currents, a second order converging curvature evaluation technique for level set methods is presented. The results of several different test cases demonstrate the effectiveness of the proposed method, showing good mass conservation properties and second order converging spurious current magnitudes.