Screening effect of plasma flow on the resonant magnetic perturbation penetration in tokamaks based on two-fluid model

Numerical simulation on the resonant magnetic perturbation penetration is carried out by the newly-updated initial value code MDC(MHD@Dalian Code).Based on a set of two-fluid fourfield equations,the bootstrap current,parallel,and perpendicular transport effects are included appropriately.Taking into account the bootstrap current,a mode penetration-like phenomenon is found,which is essentially different from the classical tearing mode model.To reveal the influence of the plasma flow on the mode penetration process,E×B drift flow and diamagnetic drift flow are separately applied to compare their effects.Numerical results show that a sufficiently large diamagnetic drift flow can drive a strong stabilizing effect on the neoclassical tearing mode.Furthermore,an oscillation phenomenon of island width is discovered.By analyzing it in depth,it is found that this oscillation phenomenon is due to the negative feedback regulation of pressure on the magnetic island.This physical mechanism is verified again by key parameter scanning.

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.

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.

A kind of extended Korteweg-de Vries equation and solitary wave solutions for interfacial waves in a two-fluid system

This paper considers interfacial waves propagating along the interface between a two-dimensional two-fluid with a flat bottom and a rigid upper boundary. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. It just focuses on the weakly non-linear small amplitude waves by introducing two small independent parameters： the nonlinearity ratio ε, represented by the ratio of amplitude to depth, and the dispersion ratio μ, represented by the square of the ratio of depth to wave length, which quantify the relative importance of nonlinearity and dispersion. It derives an extended KdV equation of the interfacial waves using the method adopted by Dullin et al in the study of the surface waves when considering the order up to O（μ^2）. As expected, the equation derived from the present work includes, as special cases, those obtained by Dullin et al for surface waves when the surface tension is neglected. The equation derived using an alternative method here is the same as the equation presented by Choi and Camassa. Also it solves the equation by borrowing the method presented by Marchant used for surface waves, and obtains its asymptotic solitary wave solutions when the weakly nonlinear and weakly dispersive terms are balanced in the extended KdV equation.

Two-fluid oscillatory flow in a channel

The validity of Navier’s partial slip condition is investigated by studying the oscillatory flow in a coated channel.The two-fluid model is used to solve the unsteady viscous equations exactly.Partial slip is experienced by the core fluid.It is found that Naviers condition does not hold for an unsteady core fluid.

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.

An Interacting Two-Fluid Scenario for Quintom Dark Energy

The Quintom dark energy is a proposal that explains the recent observations that mildly favor the equation of state of dark energy w crossing -1 near the past. The Quintom model is often constructed by two scalar fields, where one is the quintessence field and another is the phantom field. The cosmological implication of the coupling of the two fields of the dark energy is out of question worth investigating. However, the consideration of the coupling in the field scenario is somewhat complex thus we propose an interacting two-fluid Quintom scenario for simplicity. The interaction between the two components is parametrized by a constant η in this scenario. The cosmological implications of this parametrization are investigated in detail in this paper. Also, a diagnostic for this model is performed by using the statefinder pairs {s, r} and {q, r}.

Higher-order Boussinesq-type equations for interfacial waves in a two-fluid system

Interracial waves propagating along the interface between a three-dimensional two-fluid system with a rigid upper boundary and an uneven bottom are considered. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. A set of higher-order Boussinesq-type equations in terms of the depth-averaged velocities accounting for stronger nonlinearity are derived. When the small parameter measuring frequency dispersion keeping up to lower-order and full nonlinearity are considered, the equations include the Choi and Camassa＇s results （1999）. The enhanced equations in terms of the depth-averaged velocities are obtained by applying the enhancement technique introduced by Madsen et al. （1991） and Schaffer and Madsen （1995a）. It is noted that the equations derived from the present study include, as special cases, those obtained by Madsen and Schaffer （1998）. By comparison with the dispersion relation of the linear Stokes waves, we found that the dispersion relation is more improved than Choi and Camassa＇s （1999） results, and the applicable scope of water depth is deeper.

Validation of an improved two-fluid model with particle rotation for gas-solid fluidized bed

Gas-solid fluidized beds are widely applied in chemical and process engineering.It is of significance to establish a reasonable and effective mathematical model to explore the hydrodynamics of gas-particle system for industrial applications.As a less computationally demanding alternative to the discrete descriptions,two-fluid model considering kinetic theory of granular flow is often adopted to describe the fluidized behaviors of particles,but it cannot characterize the rotation of particles and its influence on the fluidized behaviors.In this study,to address the rotation effect of the fluidized particles,a two-fluid model combining the classical fluid and micropolar fluid is established,namely CMTFM.In the CMTFM,classical fluid is used to describe the motion of gas phase,while micropolar fluid is adopted to describe the motion of particle phase,and the rotation of particles and its influence on the hydrodynamics of the gas-particle system are characterized by the degree of freedom of microrotation and the improved drag force based on micropolar viscosities.In the calculation of the gas-solid bubbling fluidized bed,we investigated the influence of the microstructure parameters,particle-particle collision restitution coefficient and inlet velocity,and the results are compared to those from TFM model and experiments.Through the analysis,it manifests that pressure drop and expansion height of the fluidized bed under the consideration of the microrotation effect are closer to the experiments,which demonstrates the feasibility and advantage of the classical-micropolar two-fluid model.

压裂复杂裂缝中支撑剂输送数值模拟研究

在非常规储层水力压裂施工过程中,水力裂缝容易与天然裂缝交汇形成复杂裂缝,支撑剂的运移和铺置形态直接决定了储层的增产效果。为了研究支撑剂在复杂裂缝中的铺置规律,基于欧拉-欧拉方法建立了支撑剂-压裂液两相流数学模型,并利用复杂裂缝室内实验装置与数值模拟的砂堤铺置形态对比的方式进行了模型验证,结果表明,选取的模型可用于研究滑溜水在复杂裂缝内的携带支撑剂的运移铺置规律。以相似准则为基础,建立了简化的复杂裂缝平板模型,并将砂堤形态参数、砂堤面积进行归一化处理,获得了压裂液排量、黏度、支撑剂粒径、裂缝宽度以及裂缝形态对支撑剂在复杂裂缝中铺置的影响规律。结果表明:增大排量、黏度和减小粒径均有利于支撑剂向裂缝深处运移,并且排量对支缝中的砂堤形态影响最明显,但增大排量不利于支撑剂填充近井地带。缝宽的影响体现在壁面效应,在注入时间相同的情况下,当压裂液黏度从1 mPa·s增加到5 mPa·s时,主缝砂堤的归一化长度平均增加了37.9%,归一化高度平均降低了61.4%。裂缝结构越复杂,所有支缝中的支撑剂铺置面积占比越高,分流作用越大。随着复杂裂缝的支缝数量、级数以及支缝延伸长度的增大,支缝中的砂堤高度与长度均有所降低;相对于T1型裂缝,T3和T5型裂缝中的砂堤长度、高度减小最多。

Effects of diamagnetic drift on nonlinear interaction between multi-helicity neoclassical tearing modes

A numerical study of the diamagnetic drift effect on the nonlinear interaction between multi-helicity neoclassical tearing modes(NTMs) is carried out using a set of four-field equations including two-fluid effects.The results show that,in contrast to the single-fluid case,5/3 NTM cannot be completely suppressed by 3/2 NTM with diamagnetic drift flow.Both modes exhibit oscillation and coexist in the saturated phase.To better understand the effect of the diamagnetic drift flow on multiple-helicity NTMs,the influence of typical relevant parameters is investigated.It is found that the average saturated magnetic island width increases with increasing bootstrap current fraction f_(b) but decreases with the ion skin depth δ.In addition,as the ratio of parallel to perpendicular transport coefficients χ_(‖)/χ_(⊥) increases,the average saturated magnetic island widths of the 3/2 and 5/3 NTMs increase.The underlying mechanisms behind these observations are discussed in detail.

径向油孔结构对环下润滑高速轴承内部流动特性的影响研究

为了研究径向油孔结构对环下润滑高速轴承内部流动特性的影响规律,根据高速轴承内部流动特征建立了数值计算方法,针对包含输油通道、径向油孔和轴承组件的物理仿真模型开展了油气两相流动计算,对比并讨论了不同径向油孔结构下轴承内部的油气分布和黏性摩擦损失。数值模拟结果表明,径向油孔孔径增大后轴承内部的平均滑油体积分数单调增加且滑油分布更均匀。径向油孔布设在周向油槽的同侧有利于提升轴承内部的平均滑油体积分数,同时滑油沿周向分布的均匀性较好。轴承组件表面扩展参数受径向油孔结构的影响,其变化趋势与平均滑油体积分数的变化相似。全部采用经验公式预测得到的黏性摩擦损失整体偏高,且不能反映径向油孔结构参数的影响,数值模拟和经验公式相结合计算得到的黏性摩擦损失与直接采用数值模拟获得的结果表现出较好的一致性。

声共振混合器高黏度流体混合的功耗特性研究

声共振混合利用机械共振产生高加速度振动,从而促进流体流动,其功耗特性对于其设计及应用具有重要作用。为研究声共振混合器的功耗特性,基于CFD建立了声共振混合过程仿真模型,分析了高黏度流体声共振混合过程壁面对物料的作用力和做功功率,探究了黏度和振动参数改变对混合器功耗特性的影响,并建立了声共振混合器混合功率的预测函数。研究结果表明,在混合过程中,壁面对液相做功的瞬时功率呈现先减少后稳定波动的趋势,而有效功率呈现先增加后稳定波动的趋势,这种不同的变化趋势是由于两者的相位差发生变化所导致的。增加振幅、频率或等加速度下低频大振幅都能够增加瞬时功率和有效功率,并减少液相进入稳定流动阶段所需吸收的外界能量。

煤岩体突水通道骨料灌注运移堆积机制试验

为揭示突水通道中骨料运移堆积的演化规律,开展了动水条件下骨料灌注试验,得到了骨料的初始沉积位置、堆积体形成方向以及堆积体极限堆积高度与水流速、骨料灌注速度等影响因素的内在关联。结果表明:在1.0 m/s高流速动水环境中,骨料堆积体形成方向为先逆水流方向生长至投料孔底后转为顺水流方向生长;在0.3 m/s低流速条件下,采用口径36 mm漏斗投料同口径18 mm相比,骨料灌注总量相同时,其堆积体的极限高度提升了41%。基于数值力学模型——计算流体动力学-离散元方法(computational fluid dynamics and distinct element method,CFD-DEM),研究了骨料灌注期间突水通道内流场的时空分布特征以及骨料堆积体的过水阻力作用对突水通道内静水压力的影响机制。

喷油润滑球轴承流动特性的数值研究

基于流体体积(VOF)方法建立了喷油润滑球轴承的数值仿真模型.分别采用滑移网格法和多重参考系(MRF)方法对轴承油气两相流动进行了数值计算,得到了轴承内部流场分布和搅油力矩.搭建了轴承流场可视化试验台,对比分析了仿真流场与高速摄影实测流场.结果表明:滑移网格下的轴承内部油液体积分数大于用MRF方法的计算结果,且当转速逐渐增大时,轴承内部油液体积分数均逐渐减小,且在离心力作用下,润滑油更多地分布在轴承外圈滚道.滑移网格下的轴承内部流场油液分布更接近实测流场.仿真计算的搅油力矩随转速增大而增大,MRF方法下搅油力矩大于滑移网格下的搅油力矩,平均误差约为12.3%.两种方法对搅油力矩的计算均具有较好的一致性,但MRF时间成本更低.

PEMFC船形堵块阴极流场的性能

建立流体体积(VOF)两相流模型,研究船形堵块流道的排水性能。建立不同开孔率的船形堵块流道三维模型,研究船形堵块及开孔率对质子交换膜燃料电池(PEMFC)性能的影响。船形堵块流道峰值净功率密度相较于传统直流道可提升9.4%,相较于同堵塞率下的梯形堵块流道可提高2.9%,具有较好的强化传质作用。通过提高船形堵块流道的开孔率,PEMFC可以获得更好的性能。在相同开孔率下,船形堵块流道的氧气摩尔浓度较直流道和梯形堵块流道分别提高28.6%和14.0%。结果表明,相较于传统直流道,船形堵块流道可降低排水周期和流道内平均水含量,具有更好的排水性能。

基于多孔介质方法的钠冷快堆冷却剂沸腾现象模拟

基于两流体六方程模型针对钠的气液两相分别构建守恒方程,采用蒸发冷凝模型表征两相质量交换,分别使用显式和隐式处理方法对蒸发冷凝模型进行计算,同时考虑了Sobolev阻力模型、两相对流换热模型以及相间动量交换等本构关系,开发了适用于模拟钠冷快堆冷却剂沸腾的多孔介质分析方法,利用KNS-37失流实验L22工况数据进行了对比验证,并利用L29工况流量数据验证模型的适用性。结果表明,所建立的钠沸腾多孔介质分析方法可以较好地模拟钠冷快堆沸腾现象,预测沸腾发生时间在6.3 s左右,与实验相差0.2 s,温度和流量的总体变化趋势与实验数据吻合较好。

恶劣海况下船舶砰击颤振响应特性数值计算与试验研究

[目的]针对恶劣海况下船舶所受砰击颤振响应现象,探究船舶非线性波浪载荷与瞬态高幅值砰击载荷的耦合作用。[方法]采用计算流体动力学与有限元方法(CFD-FEM)相结合的双向流固耦合方法对S175集装箱船进行数值仿真计算,并与试验结果及切片理论计算结果进行对比验证;采用分段变截面弹性龙骨梁模型开展船舶的砰击颤振特性模型试验,基于CFD-FEM双向流固耦合方法开展船艏砰击载荷及高频非线性砰击颤振响应特性分析,并与模型试验结果进行对比验证。[结果]结果显示,波浪砰击载荷对船艏颤振响应的影响不可忽视,6级海况下由砰击颤振诱发的二阶高频成分分量占低频波浪弯矩的59.86%。[结论]采用基于CFD-FEM的双向流固耦合方法可准确计算船首砰击颤振响应;在高海况下船舶所受非线性波浪载荷及结构动态响应易受船首瞬态砰击载荷的影响,在船舶结构设计与安全评估中需考虑高频砰击颤振的情况。

Two-fluid modeling of Geldart A particles in gas-fluidized beds

We have investigated the effect of cohesion and drag models on the bed hydrodynamics of Geldart A particles based on the two-fluid （TF） model. For a high gas velocity U0 = 0.03 m/s, we found a transition from the homogeneous fluidization to bubbling fluidization with an increase of the coefficient C1, which is used to account for the contribution of cohesion to the excess compressibility. Thus cohesion can play a role in the bed expansion of Geldart A particles. Apart from cohesion, we have also investigated the influence of the drag models. When using the Wen and Yu drag correlation with an exponent n = 4.65, we find an under-prediction of the bed expansion at low gas velocities （U0 = 0.009 m/s）. When using a larger exponent （n = 9.6）, as reported in experimental studies of gas-fluidization, a much better agreement with the experimental bed expansion is obtained. These findings suggest that at low gas velocity, a scale-down of the commonly used drag model is required. On the other hand, a scale-up of the commonly used drag model is necessary at high gas velocity （U0 = 0.2 and 0.06 m/s）. We therefore conclude that scaling the drag force represent only an ad hoc way of repairing the deficiencies of the TF model, and that a far more detailed study is required into the origin of the failure of the TF model for simulating fluidized beds of fine powders.