The Crank–Nicolson finite spectral element method and numerical simulations for 2D non-stationary Navier–Stokes equations

In this paper, we first build a semi-discretized Crank–Nicolson (CN) model about time for the two-dimensional (2D) non-stationary Navier–Stokes equations about vorticity–stream functions and discuss the existence, stability, and convergence of the time semi-discretized CN solutions. And then, we build a fully discretized finite spectral element CN (FSECN) model based on the bilinear trigonometric basic functions on quadrilateral elements for the 2D non-stationary Navier–Stokes equations about the vorticity–stream functions and discuss the existence, stability, and convergence of the FSECN solutions. Finally, we utilize two sets of numerical experiments to check out the correctness of theoretical consequences.     

Optimal control and stability analysis of an online game addiction model with two stages

In this paper, we establish a mathematical model of online game addiction with two stages to research the dynamic properties of it. The existence of all equilibria is obtained, and the basic reproduction number is calculated by the method of next-generation matrix. The global asymptotic stability of disease-free equilibrium (DFE) is proved by comparison principle, and the global asymptotic stability of endemic equilibrium (EE) is proved by constructing an appropriate Lyapunov function. Then we use the Pontryagin's maximum principle to find the optimal solution of the model, so that the number of infected people can be minimized. In numerical simulation, firstly, we validate the global stability of DFE and EE. Secondly, we consider three kind of control measures (treatment, isolation, and education) and divide them into four cases. The models with control and without control are solved numerically by using forward and backward sweep Runge-Kutta method. In order to achieve the best control effect, we suggest that three kind of measures should be used simultaneously according to the optimal control strategy.     

Global stability of a multistrain SIS model with superinfection and patch structure

We study the global stability of a multistrain SIS model with superinfection and patch structure. We establish an iterative procedure to obtain a sequence of threshold parameters. By a repeated application of a result by Takeuchi et al. [Nonlinear Anal Real World Appl. 2006;7:235–247], we show that these parameters completely determine the global dynamics of the system: for any number of patches and strains with different infectivities, any subset of the strains can stably coexist depending on the particular choice of the parameters.     

过载作用下聚乙烯疲劳裂纹扩展行为研究

论文针对中密度聚乙烯材料(MDPE),采用平板试样进行了I型疲劳裂纹扩展和单次过载下裂纹扩展试验．发现与金属材料类似,单次拉伸过载对聚乙烯(PE)的疲劳裂纹扩展有明显的迟滞作用,降低了裂纹扩展速率．试验还通过变载荷刻线法获取疲劳裂纹扩展前缘的实际形貌和变化规律,对常规变载荷刻线方法进行了调整和验证,其修正方法对高分子材料的疲劳裂纹扩展前缘刻线具有较好的效果．通过观察发现含楔形塑性区的裂尖钝化是裂纹迟滞的主要原因．过载引入的塑性区内残余应力对裂纹迟滞也起了重要作用．论文利用Dugdale模型计算了塑性区尺寸,使用基于残余应力的Wheeler模型对过载迟滞进行了很好的拟合．     

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

A proper orthogonal decomposition variational multiscale meshless interpolating element-free Galerkin method for incompressible magnetohydrodynamics flow

In the recent decade, the meshless methods have been handled for solving most of PDEs due to easiness of the meshless methods. One of the popular meshless methods is the element-free Galerkin (EFG) method that was first proposed for solving some problems in the solid mechanics. The test and trial functions of the EFG are based on the special basis. Recently, some modifications have been developed to improve the EFG method. One of these improvements is the variational multiscale EFG procedure. In the current article, the shape functions of interpolation moving least squares approximation have been applied to the variational multiscale EFG technique for solving the incompressible magnetohydrodynamics flow. In order to reduce the elapsed CPU time of simulation, we employ a reduced-order model based on the proper orthogonal decomposition technique. The current combination can be referred to as the reduced-order variational multiscale EFG technique. To illustrate the reduction in CPU time used as well as the efficiency of the proposed method, we applied it for the two-dimensional cases.     

无轴承旋翼柔性梁载荷特性研究

应用Hamilton原理建立了双路传力的无轴承旋翼运动方程。采用均匀入流模型,基于直升机飞行平衡条件,建立了无轴承旋翼柔性梁载荷的计算模型,并通过算例验证了模型的精度。利用该模型,研究了全机重心位置、机身气动阻力以及平尾安装角对柔性梁载荷特性的影响,给出了各因素对柔性梁载荷的影响趋势,得出了降低柔性梁载荷的方法。数值结果表明:2cm左右重心位置的变化能够引起9%~11%的柔性梁载荷变化量,15%气动阻力的增加会导致约9%的柔性梁载荷的增大;2°平尾安装角的变化引起约10%柔性梁载荷的变化量,3°平尾安装角的变化引起约26%柔性梁载荷的变化。     

桉树杂交种与其亲本的近红外光谱判别

研究桉树控制授粉后目标性状的基因作用方式是探索其基因重组规律的重要内容。常规的数量统计分析精度往往不高,而DNA分析的专业要求高,且费时费力。该研究利用近红外光谱(NIRs)研究不同基因型桉树杂交种、亲本及杂交种与亲本间近红外光谱信息的关系,探索NIRs用于桉树杂交种与其亲本判别的可行性和准确性。以控制授粉的桉树亲本及其杂交F1代材料为对象,每种基因型从各自田间试验分别选取10个单株,采集树冠中上部新鲜健康叶片。用手持式近红外仪Phazir Rx(1624)采集桉树杂交种与其亲本叶片的NIRs信息。每单株选10片完全生理成熟的健康叶片,避开叶脉扫描其正面光谱5次,以50条NIRs信息的均值代表单个叶片的NIRs信息,最终每个基因型获得10条NIRs信息。对原始NIRs采用二阶多项式S.G一阶导数预处理。预处理后的NIRs用于多元统计分析,首先对桉树杂交亲本和子代样本进行主成分分析(PCA),直观展示不同基因型的分类情况。然后运用簇类独立软模式(SIMCA)和偏最小二乘判别分析(PLS-DA)两种有监督的判别模式验证NIRs用于桉树杂交种与其亲本树种的分类判别效果。PCA结果显示,不同的亲本间、杂交种间及杂交种与亲本间样本的主因子得分可以清晰地将各基因型分开。SIMCA模式判别分析中,桉树杂交种样本到亲本PCA模型的样本距离显示,待判别样本能够形成单独的聚类,且能直观反映两者的遗传相似。PLS-DA判别结果显示,桉树杂交亲本的PLS模型能通过预测其杂交子代的响应变量将其与亲本准确分开。结果表明,桉树叶片的NIRs信息可以准确地反映桉树杂交子代遗传信息的传递规律,NIRs判别模型可以准确地将各种基因型予以区分。因此,NIRs信息不仅可用于桉树杂交种和纯种的定性判别,还可以分析桉树基因重组过程中加性遗传效应的大小,从而为桉树遗传基础分析及其育种改良研究提供理论支撑。     

AlGaN/GaN高电子迁移率器件外部边缘电容的物理模型

AlGaN/GaN HEMT外部边缘电容Cofd是由栅极垂直侧壁与二维电子气水平壁之间的电场构成的等效电容.本文基于保角映射法对Cofd进行物理建模,考虑沟道长度调制效应,研究外部偏置、阈值电压漂移和温度变化对Cofd的影响:随着漏源偏压从零开始增加,Cofd先保持不变再开始衰减,其衰减速率随栅源偏压的增加而减缓;AlGaN势垒层中施主杂质浓度的减小和Al组分的减小都可引起阈值电压的正向漂移,正向阈值漂移会加强沟道长度调制效应对Cofd的影响,导致Cofd呈线性衰减.在大漏极偏压工作情况下,Cofd对器件工作温度的变化更加敏感.