Mg,N掺杂β-Ga2O3光电性质的第一性原理计算

通过基于密度泛函理论的第一性原理计算,研究了Mg单掺杂、N单掺杂和不同浓度的Mg-N共掺杂β-Ga₂O₃的结构性质、电子性质和光学性质,以期获得性能比较优异的p型β-Ga₂O₃材料。建立了五种模型:Mg单掺杂、N单掺杂、1个Mg-N共掺杂、2个Mg-N共掺杂和3个Mg-N共掺杂β-Ga₂O₃。经过计算,3个Mg-N共掺杂β-Ga₂O₃体系的结构最稳定。此外,在5种模型中,3个Mg-N共掺杂β-Ga₂O₃体系的禁带宽度是最小的,并且N 2p和Mg 3s贡献的占据态抑制了氧空位的形成,从而增加了空穴浓度。因此,3个Mg-N共掺杂β-Ga₂O₃体系表现出优异的p型性质。3个Mg-N共掺杂体系的吸收峰出现明显红移,在太阳盲区的光吸收系数较大,这归因于导带Ga 4s、Ga 4p、Mg 3s向价带O 2p、N 2p的带间电子跃迁。本工作将为p型β-Ga₂O₃日盲光电材料的研究和应用提供理论指导。     

Tire and soil effects on power loss: Measurement and comparison with finite element model results

In the present study, the effect of vertical load, tire inflation pressure and soil moisture content on power loss in tire under controlled soil bin conditions were investigated. Also a finite element model of tire-soil interaction in order to achieve a suitable model for predicting power loss in tire was created. Increasing the vertical load on the tire had a noteworthy impact on increasing the tire contact volume with the soil, reducing the percentage of slip, and increasing the rolling resistance; although, reducing the load on the tire had the opposite effect. At a constant inflation pressure, by increasing the vertical load on the tire, the amount of power loss due to the rolling resistance and the total power loss in the tire increased. Increase in soil moisture content increased the power loss caused by slip. Increasing the inflation pressure at a constant vertical load, also increasing the soil moisture content, led to an increase in the power loss caused by rolling resistance, and increase total power loss. The obtained error for estimating power loss of rolling resistance and total power loss was satisfactory and confirmed the acceptability of the model for power loss estimation.     

Acoustoelastic effect simulation by time–space finite element formulation based on quadratic interpolation of the acceleration

Acoustoelastic effect describes the change of ultrasound velocity due to the initial stress. Its simulation involves a numerical analysis of nonlinear elastodynamics and requires high accuracy in the time domain. A time–space finite element formulation, derived from the quadratic interpolation of the acceleration within a time segment, is proposed for an accurate simulation of the acoustoelastic effect in the present study. Ten different integration schemes are generated based on this formulation and nine of them are found to be conditionally stable. Among the nine stable schemes, one is found to obtain a spectral radius of one when the normalized step ratio is less than 5.477, indicating no numerical dissipation or numerical divergence. Compared with integration schemes from previous studies, this integration scheme demonstrates better performance in calculation accuracy and energy conservation. A two-stage approach, namely the static stage and the dynamic stage, has been employed in the simulation of the acoustoelastic effect. The former stage is adopted to obtain the initial stress and the latter stage, where the proposed integration scheme is implemented, is adopted to simulate the ultrasound propagation in an initial stress state. The simulation results of the dynamic stage show that the ultrasound velocity increases in a compression stress state and decreases in a tension stress state for aluminum alloy, which is in good agreement with previous experimental studies. Together with the simulation result of the static stage, it is conjectured that the acoustoelastic effect results from the stress-dependent elastic modulus.     

基于SWT方法的钢绞线索微动疲劳特性分析

为得到钢绞线索丝间接触区的应力场分布并预测微动疲劳裂纹萌生位置和微动疲劳寿命,本文利用参数化方法建立了精细化的钢绞线拉索有限元模型,包括整索模型和不同层丝间接触区域的局部精细化子模型．分析了钢绞线索在两种交变荷载工况下的应力场变化情况,并基于多轴疲劳SWT(Smith-Watson-Topper)临界平面法进行了疲劳特性分析和疲劳寿命预测．主要结论如下：钢绞线索内接触区边缘处的微动幅值较大,中心处几乎没有相对滑动,微动疲劳的初始裂纹萌生点位于接触区域边缘;经不同区域子模型分析比较,在轴向循环荷载作用下,外层钢丝的接触区域比内层钢丝更易发生微动疲劳损伤;在横向位移循环荷载作用下,同层钢丝因位置角度不同而产生了较大的疲劳特性差异,且相比轴向循环拉伸,该工况下最不利单丝的微动疲劳寿命更低;与非接触区域相比,接触区的疲劳寿命大幅降低,微动现象对钢绞线索的抗疲劳性能有明显降低作用．     

A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

颗粒毛细效应影响因素的离散元分析

颗粒毛细效应是指将一根细管插入填充有颗粒物质的容器中并对管施加竖直振动时颗粒在管内上升并最终达到一个稳定的高度的现象, 该现象为颗粒物料的逆重力输运提供了一种潜在的技术途径. 为探究颗粒毛细效应的影响因素, 采用离散元方法, 模拟再现了颗粒毛细效应过程,展示了不同管径下颗粒竖直方向速度演变特性, 考察了不同容器宽度和振动条件下颗粒最终毛细上升高度随管径的演变规律. 结果表明, 在容器宽度与粒径比为40、管振幅与粒径比为14.33、管振动频率为12 Hz情况下, 管径与粒径比$D/d = 3.33$时, 管内颗粒堵塞严重, 使得颗粒上升缓慢,并造成颗粒柱中断; $D/d = 8.33$时, 起初毛细上升高度增加迅速, 随后毛细上升高度的增大逐渐减缓, 管内颗粒在管径方向几乎不存在速度梯度; $D/d =15$时, 随着颗粒毛细上升高度的增大, 管内颗粒柱分离为速度截然不同的两层, 上层颗粒在管径方向几乎不存在速度梯度, 而下层颗粒存在明显的速度梯度.研究还发现, 在毛细效应能够发生的管径范围内, 存在一个对应于颗粒最终毛细上升高度最大值的临界管径, 当管径小于临界管径时, 颗粒最终毛细上升高度随管径的增大而增大, 当管径大于临界管径时, 颗粒最终毛细上升高度随管径的增大而趋于减小; 增大容器宽度,临界管径有所增大; 增大振幅、适当提高频率能够有效促进临界管径的增大.      

基于特征正交分解的一类瞬态非线性热传导问题的新型快速分析方法

提出了一种基于特征正交分解(POD)和有限元法的瞬态非线性热传导问题的模型降阶快速分析方法, 建立了导热系数随温度变化的一类瞬态非线性热传导问题有限元格式的POD降阶模型. 在隐式时间推进方法的基础上有效结合单元预转换方法和多级线性化方法发展了一种加速求解瞬态非线性热传导降阶模型的新型计算方法,并通过二维和三维算例验证了该方法的准确性和高效性. 研究结果表明: (1)降阶模型解的均方根误差在经过初始时段轻微的脉动后稳定于0.01%以下, 而其计算效率比有限元全阶模型提高2$\sim $3个数量级, 并且自由度数量(DOFs)愈大提高的幅度也愈加显著; (2)新型算法解决了常规算法在计算非线性降阶模型时加速性能差的问题, 即使是在DOFs比较小的时候也能够明显提高计算效率; (3)常数边界条件下得到的POD模态可以用来建立相同求解域在各种复杂时变边界条件下的瞬态非线性热传导降阶模型, 并对其传热过程和温度场进行快速准确的分析与预测, 具有很好的工程应用价值.      

A hybrid numerical approach to predict the vibrational responses of panels excited by a turbulent boundary layer

In this work, a hybrid numerical approach to predict the vibrational responses of planar structures excited by a turbulent boundary layer is presented. The approach combines an uncorrelated wall plane wave technique with the finite element method. The wall pressure field induced by a turbulent boundary layer is obtained as a set of uncorrelated wall pressure plane waves. The amplitude of these plane waves are determined from the cross spectrum density function of the wall pressure field given either by empirical models from literature or from experimental data. The response of the planar structure subject to a turbulent boundary layer excitation is then obtained from an ensemble average of the different realizations. The numerical technique is computationally efficient as it rapidly converges using a small number of realizations. To demonstrate the method, the vibrational responses of two panels with simply supported or clamped boundary conditions and excited by a turbulent flow are considered. In the case study comprising a plate with simply supported boundary conditions, an analytical solution is employed for verification of the method. For both cases studies, numerical results from the hybrid approach are compared with experimental data measured in two different anechoic wind tunnels.     

Investigation on the effect of density ratio on the convergence behavior of partitioned method for fluid–structure interaction simulation

In order to investigate the effect of density ratio of fluid and solid on the convergence behavior of partitioned FSI algorithm, three strong-coupling partitioned algorithms (fixed-point method with a constant under-relaxation parameter, Aitken’s method and Quasi-Newton inverse least squares (QN-ILS) method) have been considered in the context of finite element method. We have employed the incompressible Navier–Stokes equations for a Newtonian fluid domain and the total Lagrangian formulation for a non-linear motion of solid domain. Linear-elastic (hyper-elastic) model has been employed for solid material with small (large) deformation. A pulsatile inlet-flow interacting with a 2D circular channel of linear-elastic material and a pressure wave propagation in a 3D flexible vessel have been simulated. Both linear-elastic and hyper-elastic (Mooney–Rivlin) models have been adopted for the 3D flexible vessel. From the present numerical experiments, we have found that QN-ILS outperforms the others leading to a robust convergence regardless of the density ratio for both linear-elastic and hyper-elastic models. On the other hand, the performances of the fixed-point method with a constant under-relaxation parameter and the Aitken’s method depend strongly on the density ratio, relaxation parameter selected for coupling iteration, and degree of deformation. Although the QN-ILS of this work is still slower than a monolithic method for serial computation, it has an advantage of easier parallelization due to the modularity of the partitioned FSI algorithm.     

Stochastic sensitivity analysis of coupled acoustic-structural systems under non-stationary random excitations

This paper presents a new sensitivity analysis method for coupled acoustic–structural systems subjected to non-stationary random excitations. The integral of the response power spectrum density (PSD) of the coupled system is taken as the objective function. The thickness of each structural element is used as a design variable. A time-domain algorithm integrating the pseudo excitation method (PEM), direct differentiation method (DDM) and high precision direct (HPD) integration method is proposed for the sensitivity analysis of the objective function with respect to design variables. Firstly, the PEM is adopted to transform the sensitivity analysis under non-stationary random excitations into the sensitivity analysis under pseudo transient excitations. Then, the sensitivity analysis equation of the coupled system under pseudo transient excitations is derived based on the DDM. Moreover, the HPD integration method is used to efficiently solve the sensitivity analysis equation under pseudo transient excitations in a reduced-order modal space. Numerical examples are presented to demonstrate the validity of the proposed method.