Particulate flow modelling in a spiral separator by using the Eulerian multi-fluid VOF approach

The Euler-Euler model is less effective in capturing the free surface of flow film in the spiral separator,and thus a Eulerian multi-fluid volume of fluid(VOF)model was first proposed to describe the particulate flow in spiral separators.In order to improve the applicability of the model in the high solid concentration system,the Bagnold effect was incorporated into the modelling framework.The capability of the proposed model in terms of predicting the flow film shape in a LD9 spiral separator was evaluated via comparison with measured flow film thicknesses reported in literature.Results showed that sharp air–water and air-pulp interfaces can be obtained using the proposed model,and the shapes of the predicted flow films before and after particle addition were reasonably consistent with the observations reported in literature.Furthermore,the experimental and numerical simulation of the separation of quartz and hematite were performed in a laboratory-scale spiral separator.When the Bagnold lift force model was considered,predictions of the grade of iron and solid concentration by mass for different trough lengths were more consistent with experimental data.In the initial development stage,the quartz particles at the bottom of the flow layer were more possible to be lifted due to the Bagnold force.Thus,a better predicted vertical stratification between quartz and hematite particles was obtained,which provided favorable conditions for subsequent radial segregation.

High-order targeted essentially non-oscillatory scheme for two-fluid plasma model

The weakly ionized plasma flows in aerospace are commonly simulated by the single-fluid model,which cannot describe certain nonequilibrium phenomena by finite collisions of particles,decreasing the fidelity of the solution.Based on an alternative formulation of the targeted essentially non-oscillatory(TENO)scheme,a novel high-order numerical scheme is proposed to simulate the two-fluid plasmas problems.The numerical flux is constructed by the TENO interpolation of the solution and its derivatives,instead of being reconstructed from the physical flux.The present scheme is used to solve the two sets of Euler equations coupled with Maxwell's equations.The numerical methods are verified by several classical plasma problems.The results show that compared with the original TENO scheme,the present scheme can suppress the non-physical oscillations and reduce the numerical dissipation.

Thermophoresis effects on gas-particle phases flow behaviors in entrained flow coal gasifier using Eulerian model

A numerical model based on the Eulerian–Eulerian two-fluid approach is used to simulate the gasification of coal char inside an entrained flow gasifier. In this model, effects of thermophoresis of coal char particles are thoroughly investigated. The thermophoresis is due to the gas temperature gradient caused by absorpted heat of coal char gasification. This work, firstly, calculates the gas temperature gradient and thermophoretic force at1100 °C,1200 °C,1300 °C and 1400 °C wall temperatures. Then, the changes of particle volume fraction and velocity in the gasifier are studied in the simulation with thermophoresis or not. The results indicate that considering the particle thermophoresis has some effects on the calculation of particle volume fraction in the gasifier, especially at wall temperature of 1400 °C, and the maximum particle volume fraction variance ratio reaches up to 1.38% on wall surface of the gasifier. These effects are mainly caused by large gas temperature gradient along the radial direction of the gasifier. For the particle velocity, the changes are small but can be observable along radial direction of the gasifier, which has good agreement with the distributions of radial gas temperature gradient and thermophoretic force. These changes above may have certain effects on gasification reaction rates in this Eulerian model. So the change of gasification reaction rates in the simulation with thermophoresis or not is studied finally.

NUMERICAL SIMULATION OF WET STEAM CONDENSING FLOW WITH AN EULERIAN/EULERIAN MODEL

Wet steam condensing flow in low-pressure steam turbine leads to efficiencylosses and blade erosions. In order to investigate this problem by numerical approach, anEulerian/Eulerian model has been developed, in which the wet steam is regarded as mixture comprisingtwo coupled systems: the vapor phase and the liquid phase. These two systems are both described byconservation equations. High resolution TVD scheme is employed to capture condensing phenomena inwet steam flow. This model has been validated by numerical simulations of condensing flows in 1D and2D nozzles. Compared with experimental data, a good agreement is observed. This Eulerian/Eulerianmodel can be extended to 3D calculation of condensing flow.

Sulfur Distribution and Transport Studies in East Asia Using Eulerian Model

A three-dimensional regional Eulenan model of sulfur deposition and transport is developed.Processes treated in the model include emission,transport,diffusion,gas-phase and aqueous-phase chemical process,dry deposition,ramout and washout of sulfur.A 'looking up table' method is provided to deal with the gas-phase chemical process including sulfur transfer Dry-depositon velocity considers the influence of underlying surface,wind,degree of stability by parameterization Model calculated values reasonably agrees with observation.Distribution of sulfur deposition and transport in East Asia are also analyzed in this paper Some amount of sulfur emission of different countries transport across boundaries,but the main origin of sulfur deposition in each country in East Asia is from itself.Furthermore.some transport paths on different layers and outlet or inlet zones are found.According to sulfur balance and budget we concluded that sulfur outlets are bigger than inlets across boundary and emissions are more than deposition in most places of East Asia

A hybrid Lagrangian-Eulerian numerical model for sea-ice dynamics

A hybrid Lagrangian - Eulerian （HLE） method is developed for sea ice dynamics, which combines the high computational efficiency of finite difference method （FDM） with the high numerical accuracy of smoothed particle hydrodynamics （SPH）. In this HLE model, the sea ice cover is represented by a group of Lagrangian ice particles with their own thicknesses and concentrations. These ice variables are interpolated to the Eularian gird nodes using the Gaussian interpolation function. The FDM is used to determine the ice velocities at Eulerian grid nodes, and the velocities of Lagrangian ice particles are interpolated from these grid velocities with the Gaussian function also. The thicknesses and concentrations of ice particles are determined based on their new locations. With the HLE numerical model, the ice ridging process in a rectangular basin is simulated, and the simulated results are validated with the analytical solution. This method is also applied to the simulation of sea ice dynamics in a vortex wind field. At last, this HLE model is applied to the Bohai Sea, and the simulated concentration, thickness and velocity match the satellite images and the field observed data well.

Temperature and current sensitivity extraction of optical superconducting transition-edge sensors based on a two-fluid model

Optical superconducting transition-edge sensor(TES)has been widely used in quantum information,biological imaging,and fluorescence microscopy owing to its high quantum efficiency,low dark count,and photon number resolving capability.The temperature sensitivity(α_(I))and current sensitivity(β_(I))are important parameters for optical TESs,which are generally extracted from the complex impedance.Here we present a method to extractα_(I)andβ_(I)based on a two-fluid model and compare the calculated current-voltage curves,pulse response,and theoretical energy resolution with the measured ones.This method shows qualitative agreement that is suitable for further optimization of optical TESs.

Three-dimensional CFD simulation of inlet structure flow in pumping station based on Eulerian solid- liquid two-phase flow model

Sediment deposition in the pumping station has a huge negative impact on unit operation.The three-dimensional CFD method has been used to simulate inlet structure flow in pumping station based on the Eulerian solid- liquid two-phase flow model. The numerical results of the preliminary scheme show that sediment deposition occurs in the forebay of pumping station because of poor flow pattern therein. In order to improve hydraulic configuration in the forebay,one modified measure reconstructs water diversion weir shape,and another measure sets a water retaining sill in the approach channel. The simulation results of the modified scheme prove that back flow in the forebay has been eliminated and the sediment deposition region has also been reduced.

Anisotropic Plane Symmetric Two-Fluid Cosmological Model with Time-Varying G and A

We investigate a two-fluid anisotropic plane symmetric cosmological model with variable gravitational constant G（t） and cosmological term A（t）. In the two-fluid model, one fluid is chosen to be that of the radiation field modeling the cosmic microwave background and the other one a perfect fluid modeling the material content of the universe. Exact solutions of the field equations are obtained by using a special form for the average scale factor which corresponds to a specific time-varying deceleration parameter. The model obtained presents a cosmological scenario which describes an early acceleration and late-time deceleration. The gravitation constant increases with the cosmic time whereas the cosmological term decreases and asymptotically tends to zero. The physical and kinematical behaviors of the associated fluid parameters are discussed.

Universal critical properties of the Eulerian bond-cubic model

We investigate the Eulerian bond-cubic model on the square lattice by means of Monte Carlo simulations, using an efficient cluster algorithm and a finite-size scaling analysis. The critical points and four critical exponents of the model are determined for several values of n. Two of the exponents are fractal dimensions, which are obtained numerically for the first time. Our results are consistent with the Coulomb gas predictions for the critical O（n） branch for n 〈 2 and the results obtained by previous transfer matrix calculations. For n = 2, we find that the thermal exponent, the magnetic exponent and the fractal dimension of the largest critical Eulerian bond component are different from those of the critical 0（2） loop model. These results confirm that the cubic anisotropy is marginal at n = 2 but irrelevant for n〈2.

Numerical Simulation of Countercurrent Flow in a PWR Hot Leg by Using Two-Fluid Model

In order to evaluate CCFL （countercurrent flow limitation） characteristics in a PWR （pressurized water reactor） hot leg under reflux condensation, numerical simulations have been conducted using a 2F （two-fluid） model and a VOF （volume of fluid） method implemented in the CFD （computational fluid dynamics） software, FLUENT6.3.26. The 2F model gave good agreement with CCFL data in low pressure conditions but did not give good results for high pressure steam-water conditions. In the previous study, the computational grid and schemes were improved in the VOF method to improve calculations in circular tubes, and the calculated CCFL characteristics agreed well with the UPTF （Upper Plenum Test Facility） data at 1.5 MPa. In this study, therefore, using the 2F model and the computational grid previously improved for the VOF calculations, numerical simulations were conducted for steam-water flows at 1.5 MPa under PWR full-scale conditions. In the range of medium gas volumetric fluxes, the calculated CCFL characteristics agreed well with the values calculated by the VOF method and the UPTF data at 1.5 MPa. This indicated that the reference set of the interfacial drag correlations employed in this study could be applied not only to low pressures but also to high pressures.

正十二烷喷雾火焰中碳烟的瞬态生成特性

针对压燃式发动机中碳烟演化过程难以被实验准确捕捉的问题,研究了高温高压环境中正十二烷喷雾燃烧焰中碳烟瞬态生成特性及动力学行为,在实验方面,应用消光辐射联合技术以及燃烧成像测速技术,同步获得高温高压环境中单孔喷油器正十二烷喷雾燃烧火焰中的碳烟体积分数、碳烟温度和碳烟速度场的瞬态分布;数值计算方面,基于OpenFOAM环境,发展欧拉喷雾模型和非稳态火焰面进度变量湍流燃烧模型的耦合模型,喷雾燃烧过程的数值计算结果与实验数据有高度的一致性。研究结果表明:相比于参考工况,低喷油压力工况火焰浮起长度更短,碳烟在当量比Φ>2区间驻留时间更长,推动了碳烟的生成,碳烟峰值质量增加了80%;高温环境工况碳烟在当量比Φ>2区间驻留时间没有明显差异,但更短的火焰浮起长度使得燃料在更高的当量比开始燃烧,导致其碳烟峰值质量增加了88%;富氧工况的燃烧温度最高,碳烟粒子运动速度增大,使得碳烟在当量比Φ>2区间驻留时间明显减少,同时火焰内部当量比较小,因而碳烟峰值质量减少了42%。

基于CEL法熔滴冲击基板的TPF/FSI数值模拟

航空发动机和燃气轮机的高温工作环境会对其热端金属部件造成不可逆损伤,热障涂层(thermal barrier coatings,TBCs)能够承受高温和高压的侵蚀性环境,从而保证金属部件使用的安全性及可靠性。以空心Y_(2)O_(3)-ZrO_(2)(yttrium-stabilized zirconia,YSZ)粒子在等离子焰流中的2种形态,即全熔熔滴和空心熔滴为研究对象对其冲击铺展过程进行了模拟。基于ABAQUS/EXPLICIT耦合欧拉-拉格朗日(coupled Eulerian-Lagrangian,CEL)有限元方法,首先针对ABAQUS缺少相变模型,导致使用CEL方法计算熔滴铺展形貌失真的问题,给出了适用的动力黏度随温度变化的经验公式。此外,考虑周围空气对熔滴铺展过程的影响,提出了“两相流(two phase flow,TPF)/流-固耦合(fluid-structure-interaction,FSI)”2.5D模型,对熔滴冲击基板凝固成型及空气卷入的过程进行了模拟,并揭示了2种熔滴铺展形貌存在较大差异的机理,对制备隔热性能更优的热障涂层具有指导意义。

水平井射孔簇内支撑剂分布数值模拟

为研究水平井射孔簇内支撑剂的运移规律,建立了水平井筒几何模型,并采用欧拉双流体模型分析了不同注入流量、压裂液黏度和注入φ(砂粒)对支撑剂运移规律的影响。结果表明,在水平井筒中,支撑剂的分布向趾部偏移。增加注入流速能促进颗粒的分散,但会降低射孔内的均匀性;提高压裂液黏度能减缓颗粒沉降,有助于支撑剂的均匀分布;增加注入φ(砂粒)会增加井筒堵塞风险,因此选择注入φ(砂粒)<10%。此外,雷诺数与归一化标准差呈自然对数关系,较高的雷诺数有助于提高颗粒在射孔流出的均匀性。

履齿结构对深海沉积物扰动分析

为探究集矿车在行进过程中履带运动对深海沉积物的扰动机理,并建立适用于机械结构与沉积物作用的仿真模型,结合水槽实验和模拟分析,设定对照试验,研究机械结构运动对沉积物的扰动作用和羽状流扩散现象。实验与仿真结果表明:偏置一字形履齿低扰动性能较好,扰动最大速度在1.5 m/s以下;扰动源的下陷深度和转速对水槽中的浊度变化有显著影响,浊度计最大和最小测量值为4 381.3 mg/L和3 073.3 mg/L,相差29.9%,峰值时间相差超过50%。搭建欧拉多相流模型,设定参数和边界条件,仿真与试验结果接近,可用于后续开展小范围内的机械扰动研究。

适用欧拉-拉格朗日方法模拟气液泡状流的气泡破碎模型

欧拉-拉格朗日方法已被广泛应用于模拟鼓泡塔等气-液反应器内的流型、气泡尺寸(或气含率)及其分布。文献中该方法主要基于临界Weber数观点来描述气泡破碎行为,且破碎后的子气泡尺寸由随机数确定。但现有实验和理论研究表明,临界Weber数约束不能体现气体密度等物性参数和泡内气体重分布对气泡破碎行为的影响。针对这些不足,提出了适用欧拉-拉格朗日框架且考虑泡内气体重分布贡献的气泡破碎机理模型,并利用开源软件OpenFOAM开发了基于新破碎模型的求解器。新模型预测结果能较好地吻合实验测量的时均轴向液速、气泡尺寸及其分布等实验数据。特别地,考虑泡内气体重分布现象的破碎机理模型成功预测了实验观测的气泡尺寸双峰分布特征。

云下洗脱-Eulerian模式

在Eulerian模式中,在大气氨浓度大于硝酸的条件下,简要讨论了各种污染物对雨水组成的影响.结果表明,在本文考虑的情况下,氨对SO_2气体洗脱和SO^(2-)形成的影响最大,且能抑制雨水的酸化;过氧化氢的作用略小于氨.气溶胶对雨水SO_4^(2-)浓度的贡献小于气体.云下H_2O_2和O_3氧化能力的大小取决于雨水的pH、降落距离和各自的浓度.在雨水pH高,降落距离短和各自浓度小时,O_3的氧化能力大于H_2O_2,反之则H_2O_2大于O_3.对韶关降水酸化的模拟表明,考虑气体污染物补充机制的Eulerian模式可以较好地模拟现场降水酸化过程.

地—井瞬变电磁三维各向异性响应特征研究

地—井瞬变电磁法由于抗干扰能力强、响应信号明显等优势而应用广泛。在深部矿产资源调查和地质构造勘查等任务中,多个油气储存在层理发育的沉积岩地区大多数情况下表现为电各向异性,但是目前地—井瞬变电磁反演算法均基于各向同性模型,传统反演成像技术不再适用。为了分析地层电导率各向异性对地—井瞬变电磁勘探结果的影响,对各向异性介质模型研究其电磁响应特性。基于各向异性介质的电导率张量,提出了一种基于有限元方法的三维分析模型,采用PARDISO求解器进行求解。针对多个各向异性角度的情况模拟了各向异性介质的瞬变电磁响应,分析介质各向异性对电磁波传播规律的影响,探讨了各向异性瞬变电磁响应特征。研究结果表明:在垂直各向异性介质中,电磁波传播方向垂直于电导率最大的方向时,电磁响应变化最明显;电磁波传播方向平行于电导率最大的方向时,电磁响应变化最弱。在水平各向异性介质中,电磁波垂直于电导率最小的方向时,电磁响应变化最明显;电磁波平行于电导率最小的方向时,响应变化较弱。最后,针对地—井瞬变电磁方位测井,分析了发射源沿不同方向条件下各向异性地层和异常体模型的电磁响应,结果表明其瞬变电磁响应特征与发射源的方向有关,且对于主轴各向异性异常体,发射源方向的选择直接影响对异常体各向异性方向的判断。该研究结论对油气开发地球物理勘探方法的研究具有一定的参考价值。

群体机器人觅食任务Eulerian建模与分析

为了提高觅食任务性能,根据蚂蚁觅食原理引入了觅食概率的概念,机器人按照该概率决定是否离开Home区域,通过改变该概率值来调整环境中机器人的数量。为了研究觅食概率对系统性能的影响,利用Eulerian数学建模方法建立了觅食任务模型,并通过仿真实验分析了觅食概率对系统整体性能的影响。实验结果表明:觅食概率越大,环境中机器人数量越增加,机器人之间的干扰也增加,从而降低系统性能。最后通过实验分析给出了觅食概率与环境中任务数量以及机器人数之间的大致关系。

Clustered张拉整体结构的动力学建模

针对clustered张拉整体结构,提出了一种基于任意拉格朗日-欧拉(arbitrary Lagrangian-Eulerian, ALE)法的多体动力学模型。相较于传统的拉格朗日法模型,该研究中的模型具有更为简单的运动学约束。首先,引入了一种ALE时变长度索单元,其网格节点与物质点可以独立运动,提供了一种自然的方式描述结构中运动的滑轮和滑动绳索;其次,利用达朗贝尔原理,推导了该单元的广义力向量并计算了相应的雅可比矩阵;然后,建立了张拉整体系统的动力学方程,并利用广义α算法对其进行求解,选取节点的全局位置坐标和物质坐标作为广义坐标,其中全局位置坐标可以被不同物体共享,以减少动力学方程的自由度数和消除物体间的约束;最后,展示了一个数值算例,10层可折叠张拉整体塔架,对其折叠过程进行了准静态和动力学仿真,验证了模型的有效性。所提出的模型和算法可为clustered张拉整体结构的设计提供理论指导,具有工程意义。