突扩圆管内液-固两相流固体颗粒运动特性的DPM数值模拟

采用考虑颗粒碰撞的欧拉一拉格朗日数值模拟方法(DPM),对水平突扩圆管中液固两相流固体颗粒的碰撞过程进行了数值计算.在模型中,对液相采用欧拉法建立控制方程,对离散颗粒采用拉格朗日方法模拟.采用硬球模型描述颗粒间的碰撞作用.计算结果表明,该模型可以真实地模拟液固两相流中固体颗粒运动的动态变化过程以及颗粒的非均匀分布特征,从单颗粒层次上提供颗粒的运动信息,这有助于深入研究液固两相流中固体颗粒的运动规律.     

风沙跃移中颗粒与多粒径床面碰撞的数值模拟

采用考虑颗粒碰撞的欧拉-拉格朗日数值模拟方法,对风沙跃移中颗粒冲击多粒径床面的碰撞过程进行了数值计算。在模型中,对气相采用欧拉方法建立控制方程,对离散颗粒采用拉格朗日方法模拟,颗粒间碰撞作用采用软球模型描述。计算结果表明该模型可以模拟风沙运动中颗粒冲击多粒径床面的动态运动过程。而且在多粒径非均匀床面上的颗粒起跳具有较大的随机性。这有助于进一步揭示风沙运动中颗粒碰撞起跳机理。     

数值模拟颗粒旋转对流化床内气固两相流动特性的影响

流化床内颗粒自旋转将影响颗粒相的流动特性.本文运用基于颗粒动理学理论的欧拉-欧拉气固多相流模型,考虑颗粒自旋转流动对颗粒碰撞能量交换和耗散的影响,数值模拟流化床内气体颗粒两相流动特性.计算结果表明颗粒的自旋转使得床内更容易形成气泡,颗粒浓度分布变化增大.颗粒自旋转运动将导致床内非均匀结构更明显.     

直接模拟中不同边界条件的实施及对沉降规律的影响

在牛顿流体中, 对颗粒在4种不同边界的垂直通道中的沉降运动进行了直接数值模拟. 计算结果表明:通过计算区域随颗粒运动而移动构建的无限长通道能准确模拟颗粒自由下落到稳定沉降的发展过程; 周期性边界条件由于流场变化, 对颗粒沉降产生了影响, 不能模拟颗粒的自由沉降过程; 底部封闭边界适合模拟封闭容器内颗粒与固壁的相互作用过程, 若颗粒达到稳定沉降, 也能模拟无限长通道内的沉降过程; 流化边界适合模拟流化床内气固两相流动. 计算结果有助于更好地理解和使用不同边界条件. 关键词： 直接数值模拟 边界条件 沉降 任意拉格朗日-欧拉方法     

三维混合沙输运数值模拟

采用计算流体力学和颗粒离散元耦合的方法模拟了三维混合沙输运过程。采用体平均的Navier-stokes方程来描述气相运动,考虑了气相和颗粒相的相互作用。颗粒运动通过求解牛顿运动方程来求解,采用硬球模型描述颗粒和颗粒及颗粒和壁面的碰撞。本模型中,颗粒运动是三维的而气相运动是二维的。计算结果表明:总输沙率沿高度方向在大于2cm以上按照指数衰减,在2 cm以下则偏大;各粒径颗粒具有不同的输沙率分布,粗粒径颗粒按指数规律衰减,其它粒径颗粒输沙率随高度先指数增加后减少;各粒径颗粒平均水平速度随高度对数函数增加且同高度时随粒径增大而减小,1 cm高度以下则相反;沙粒平均粒径沿高度线性递减,2 cm以下粒径偏大。     

循环床内气固两相流中稠密颗粒间碰撞的数值模拟

循环床内部的流动属于复杂不均匀的稠密气固两相流，稠密颗粒间的相互作用是影响颗粒运动和浓度分布的不可忽视的重要因素。本文采用直接模拟Monte－Carlo算法（DSMC算法）来模拟颗粒间的相互碰撞过程，并与随机轨道模型结合起来综合考虑湍流和颗粒碰撞对颗粒运动和浓度分布的影响，模拟结果预报了床内分层流动结构和颗粒在稀相区的不均匀分布，计算结果与实验定性符合。     

鼓泡流化床埋管磨损量及其分布的数值研究

本文采用离散颗粒单元法对流化床内颗粒运动及其与固定埋管受热面的相互作用进行颗粒直接数值模拟,其中颗粒之间的碰撞采用Tsuji等提出的软球碰撞模型处理,而流场的计算采用大涡模拟,其亚网格应力为Smagorinsky涡黏性模型,流动工况为两维鼓泡流化床.磨损量的估计是基于祝京旭等人的埋管磨损试验研究的结论,并结合本文数值模拟的结果,揭示了流化床埋管磨损量及其分布的若干规律.     

气固鼓泡床的全尺度直接数值模拟

颗粒全尺度直接数值模拟由于不引入任何相间作用力封闭模型、同时又考虑了颗粒相与流体相之间的四向耦合,因而具有极高的计算精度和准度。本文采用内嵌边界多重直接力算法结合软球碰撞模型对实验室尺度的气固鼓泡流化床装置进行了全尺度直接数值模拟。计算结果真实还原了鼓泡流态化的流动情形,且能有效地捕捉装置内部流场的详细运动情况和涡结构。对流化床内颗粒在传统的离散元模型框架下进行曳力统计发现,平均曳力比全尺度模拟结果小约20%~30%。其具体数值因选取比较的曳力模型和统计网格而异。改进传统曳力模型需考虑颗粒群的非均匀性以及颗粒拟温度。     

颗粒与顺列管束磨损的数值模拟

采用直接数值模拟方法模拟了在管道中直径80μm的颗粒对顺列10×10管束的碰撞和磨损.圆管与流场流体之间的相互作用采用内嵌边界方法进行计算并得到流场结构.颗粒的运动采用了拉格朗日跟踪方法计算并与流场之间进行了双向耦合.计算分析了沿流向以及垂直于流向各排管束受到颗粒的碰撞产生的磨损量.计算结果表明沿流向第三排管束以后的各排管束应采取防磨措旌.     

颗粒圆孔射流碰撞并行直接数值模拟算法研究

本文对空间模式发展的颗粒圆孔射流碰撞进行了并行直接数值模拟算法研究。气相采用可压缩的N-S方程直接求解。颗粒相采用Lagrangian方法跟踪实际的颗粒运动。利用并行求解算法,实现了颗粒穿越边界面的模拟以及高效颗粒碰撞算法。考虑了颗粒和流体的双相耦合以及颗粒之间的碰撞。在本文的计算条件下,颗粒的直径远小于网格的间距,平均的Kolmogorov尺度和网格的间距在一个量级。气相和颗粒相的应力与实验的对比研究表明,本文的颗粒并行程序是可信的。     

Two-Level Stabilized Finite Volume Methods for the Stationary Navier-Stokes Equations

In this work, two-level stabilized finite volume formulations for the 2D steady Navier-Stokes equations are considered. These methods are based on the local Gauss integration technique and the lowest equal-order finite element pair. Moreover, the two-level stabilized finite volume methods involve solving one small Navier-Stokes problem on a coarse mesh with mesh size $H$, a large general Stokes problem for the Simple and Oseen two-level stabilized finite volume methods on the fine mesh with mesh size $h$=$\mathcal{O}(H^2)$ or a large general Stokes equations for the Newton two-level stabilized finite volume method on a fine mesh with mesh size $h$=$\mathcal{O}(|\log h|^{1/2}H^3)$. These methods we studied provide an approximate solution $(\widetilde{u}_h^v,\widetilde{p}_h^v)$ with the convergence rate of same order as the standard stabilized finite volume method, which involve solving one large nonlinear problem on a fine mesh with mesh size $h$. Hence, our methods can save a large amount of computational time.     

Numerical model reduction of 2D steady incompressible laminar flows: Application on the flow over a backward-facing step

The numerical solution of most fluid mechanics problems usually needs such a fine mesh that the associated computational times become non-negligible parts in any design process. In order to couple numerical modelling schemes with inversion or control algorithms, the size of such models needs to be highly reduced. The identification method is a way to build low-order models that fit with the original ones. The laminar flow over a backward-facing step is used as a test case. Presented solutions are found to be in good agreement with experimental and numerical results found in the literature.     

An effective explicit pressure gradient scheme implemented in the two-level non-staggered grids for incompressible Navier–Stokes equations

In this paper, an improved two-level method is presented for effectively solving the incompressible Navier–Stokes equations. This proposed method solves a smaller system of nonlinear Navier–Stokes equations on the coarse mesh and needs to solve the Oseen-type linearized equations of motion only once on the fine mesh level. Within the proposed two-level framework, a prolongation operator, which is required to linearize the convective terms at the fine mesh level using the convergent Navier–Stokes solutions computed at the coarse mesh level, is rigorously derived to increase the prediction accuracy. This indispensable prolongation operator can properly communicate the flow velocities between the two mesh levels because it is locally analytic. Solution convergence can therefore be accelerated. For the sake of numerical accuracy, momentum equations are discretized by employing the general solution for the two-dimensional convection–diffusion–reaction model equation. The convective instability problem can be simultaneously eliminated thanks to the proper treatment of convective terms. The converged solution is, thus, very high in accuracy as well as in yielding a quadratic spatial rate of convergence. For the sake of programming simplicity and computational efficiency, pressure gradient terms are rigorously discretized within the explicit framework in the non-staggered grid system. The proposed analytical prolongation operator for the mapping of solutions from the coarse to fine meshes and the explicit pressure gradient discretization scheme, which accommodates the dispersion-relation-preserving property, have been both rigorously justified from the predicted Navier–Stokes solutions.     

基于不连续因子的沸水堆芯三维瞬态数值模拟

基于不连续因子校正的粗网格有限差分法是实现堆芯瞬态三维数值模拟的高效方法之一, 粗网节块的界面不连续因子与边界反照率的计算方法决定了实时数值模拟过程中的精度。在计算不连续因子的过程中, 省去了细网节块计算与粗网均匀化过程, 直接在粗网格划分情况下, 基于节块展开法和非线性迭代策略, 推导了粗网格界面不连续因子比率与边界反照率的计算公式, 并编制了相应的计算程序。沸水堆典型算例的三维瞬态模拟证实该方法可在空间域和时间域两方面, 使静态、瞬态精度均达到与先进节块法相等同的程度, 并且计算效率优于先进节块法, 为核电站全范围模拟机三维堆芯的实时仿真模型开发提供了一种切实可行的选择。     

On the flexibility of agglomeration based physical space discontinuous Galerkin discretizations

In this work we show that the flexibility of the discontinuous Galerkin (dG) discretization can be fruitfully exploited to implement numerical solution strategies based on the use of elements with very general shapes. Thanks to the freedom in defining the mesh topology, we propose a new h-adaptive technique based on agglomeration coarsening of a fine mesh. The possibility to enhance the error distribution over the computational domain is investigated on a Poisson problem with the goal of obtaining a mesh independent discretization.The main building block of our dG method consists of defining discrete polynomial spaces directly on physical frame elements. For this purpose we orthonormalize with respect to the L²-product a set of monomials relocated in a specific element frame and we introduce an easy way to reduce the cost related to numerical integration on agglomerated meshes. To complete the dG formulation for second order problems, two extensions of the BR2 scheme to arbitrary polyhedral grids, including an estimate of the stabilization parameter ensuring the coercivity property, are here proposed.     

周期孔洞区域中热力耦合问题的双尺度有限元计算

复合材料的研究中经常遇到具有周期孔洞结构的材料,由于区域的小周期性及剧烈振荡性,用传统的有限元计算方法来计算这些材料对应的问题时需要大量的计算机存储空间及计算时间.对这类材料的热力耦合问题给出了一种新型的高阶双尺度渐近解,得到了对应的均匀化常数、均匀化方程及对应的有限元算法.数值算例表明,周期单胞的局部结构对局部应力与应变有较大的影响.算法对数值模拟这类材料的力学行为是高效和可行的. 关键词： 双尺度方法 热力耦合 周期孔洞区域 有限元方法     

驻波声场中单分散细颗粒的相互作用特性

利用外加声场促进悬浮在气相中的细颗粒发生相互作用,进而引起颗粒的碰撞和凝并,使得颗粒平均粒径增大、数目浓度降低,是控制细颗粒排放的重要技术途径.为探究驻波声场中单分散细颗粒的相互作用,建立包含曳力、重力、声尾流效应的颗粒相互作用模型,采用四阶经典龙格-库塔算法和二阶隐式亚当斯插值算法对模型进行求解.将数值模拟得到的颗粒声波夹带速度和相互作用过程与相应的解析解和实验结果进行对比,验证模型的准确性.进而研究颗粒初始条件和直径对相互作用特性的影响.结果表明,初始时刻颗粒中心连线越接近声波波动方向、颗粒位置越接近波腹点,颗粒间的声尾流效应就越强,颗粒发生碰撞所需要的时间就越短.研究还发现,颗粒直径对颗粒相互作用的影响取决于初始时刻颗粒中心连线偏离声波波动方向的程度.当偏离较小时,颗粒直径越大,颗粒发生碰撞所需要的时间越短;当偏离很大时,直径较小的颗粒能够发生碰撞,而直径较大的颗粒则无法发生碰撞.     

Electromagnetic Modelling of Fine Features in Photonic Applications

Straightforward application of numerical modelling approaches to Photonic Band-Gap structures demands the use of fine meshes, notably finer that the size of a single scattering element, in order that high accuracy is achieved. However simulations employing sufficiently fine meshing to correctly represent the geometry of the scatterers lead to very long computational run-times and huge memory consumption. Such direct numerical approaches to PBG characterisation are only suited to the analysis of single unit cells or for benchmarking. In practice, the computational overheads significantly increase when one deals with an array of cells or when modelling the global behaviour of a large number of devices integrated on a single substrate. Therefore in order to model the macroscopic response of PBG structures the possibility of employing meshes of much larger size than the individual scattering elements has been explored. A suitable approach, initially developed for Electromagnetic Compatibility predictions, has been modified to permit modelling of Photonic Crystal Waveguides. The method provides second order accuracy and leaves a designer with the flexibility to discretise the entire problem space with an arbitrary number of scatterers per mesh cell. Results are first presented for small clusters of scatterers and subsequently for complete photonic structures.     

辐射流体力学的分子动理学自适应加密方法(英文)

通过构造新的平衡分布函数和结合分区自适应网格加密方法,对不带扩散项的平衡辐射流体力学方程,构造二阶的分子动理学BGK-AMR格式.一方面在关心的计算区域中局部加密计算网格,提高计算精度的同时大大节省了计算网格数量和计算时间;另一方面,不同于已有的参数强耦合平衡分布函数,新构造的平衡分布函数中各参数不相互依赖,简化了辐射流体力学分子动理学格式的计算.一维和二维的数值算例显示了格式的性能.     

定态对流扩散方程的渐近数值方法

本文提出了一种求解定态对流扩散方程的渐近数值方法,在边界层附近不必取很细的网格,对模型问题的数值计算表明,利用中等大小的步长就可得到边界层内的数值解。