Dynamics of a Predator-Prey Model with Allee Effect and Herd Behavior

This paper deals with dynamics of a predator-prey model with Allee effect and herd behavior. We first study the stability of non-negative constant solutions for such system. We also establish the existence of Hopf bifurcation solutions for such predator-prey model. The stability and bifurcation direction of Hopf bifurcation solution in the case of spatial homogeneity are further discussed. At the same time, several examples are given by MATLAB. Finally, the numerical simulations of the system are carried out through MATLAB, which intuitively verifies and supplements the theoretical analysis results.     

A discontinuous Galerkin immersed boundary solver for compressible flows: Adaptive local time stepping for artificial viscosity–based shock-capturing on cut cells

We present a higher-order cut cell immersed boundary method (IBM) for the simulation of high Mach number flows. As a novelty on a cut cell grid, we evaluate an adaptive local time stepping (LTS) scheme in combination with an artificial viscosity–based shock-capturing approach. The cut cell grid is optimized by a nonintrusive cell agglomeration strategy in order to avoid problems with small or ill-shaped cut cells. Our approach is based on a discontinuous Galerkin discretization of the compressible Euler equations, where the immersed boundary is implicitly defined by the zero isocontour of a level set function. In flow configurations with high Mach numbers, a numerical shock-capturing mechanism is crucial in order to prevent unphysical oscillations of the polynomial approximation in the vicinity of shocks. We achieve this by means of a viscous smoothing where the artificial viscosity follows from a modal decay sensor that has been adapted to the IBM. The problem of the severe time step restriction caused by the additional second-order diffusive term and small nonagglomerated cut cells is addressed by using an adaptive LTS algorithm. The robustness, stability, and accuracy of our approach are verified for several common test cases. Moreover, the results show that our approach lowers the computational costs drastically, especially for unsteady IBM problems with complex geometries.     

改进的五阶WENO-Z+格式

WENO-ZWENO-Z$+\!$格式的性能提升依赖于新增项的作用,该项的作用是在WENO-Z格式的基础上进一步增大欠光滑子模板上的权重. 系数$\lambda$被设置用来调控该项的作用, 以避免负耗散. 本文指出了WENO-Z$+\!$格式的缺陷,其所采用$\lambda $的取值方式既不能充分发挥格式的潜力, 也未完全消除负耗散;提出$\lambda $的值应随当地流场数据变化,方能充分发挥新增项在降低数值耗散、提高分辨率上的潜力. 基于此,本文重新设计了$\lambda $的计算公式,该公式能自适应地调控新增项的作用: 只在欠光滑子模板上的权重容易过度增大的地方削弱该项的作用,以避免负耗散; 在其他地方则充分发挥新增项的作用,最大限度增大欠光滑子模板上的权重, 提高格式的分辨率.将使用该系数公式的新格式命名为WENO-Z++, 并对其数值性能进行了系统的研究.理论分析表明, 新格式在间断处具有基本无振荡(essentially non-oscillatory,ENO)特性和更低的数值耗散. 对近似色散关系(approximate dispersion relation,ADR)的研究表明,新格式有效地避免了因过度增大欠光滑子模板上的权重而带来的负耗散,其频谱特性也得到了显著提升.本文还推导了使新格式在极值点处也能保持最优阶的精度的参数设置.一系列求解Euler方程的数值试验表明,新格式的激波捕捉能力和对复杂流场结构的分辨率都显著优于原WENO-Z$+\!$格式.}     

DISCUSSION ON THE APPLICATION OF NUMERICAL EXPERIMENT IN EXPERIMENTAL TEACHING OF FLUID MECHANICS

为强化学生对流体力学课程理论知识及实验教学环节的掌握,本文针对传统实验教学过程中受限于硬件设施、模型等问题,将数值模拟技术与流体力学实验教学相结合。以伯努利实验为例,应用数值模拟软件进行物理建模、网格划分、变参数数值模拟等实验内容,通过压力云图、流线等模拟结果对能量转换特性进行了分析,增强了学生对流体力学理论知识及相关参数的物理意义、变化特征等的理解,使学生在实验中产生强烈的“参与感”,有效地提高了教学水平和教学质量。     

薄板弯曲分析的节点积分高阶无网格法

采用高阶无网格法求解薄板弯曲问题,在已发展的线性曲率光顺方案的基础上,通过引入泰勒展开技术,建立了能够精确再现纯弯曲和线性弯曲模式的节点积分方法。与之相比,目前无网格薄板分析主要采用的节点积分方法仅能精确再现纯弯曲模式。数值结果表明,本文方法可精确通过纯弯曲和线性弯曲试验,且能得到光滑、无振荡的弯矩场。与标准的高斯积分方法和目前已存在的节点积分方法相比,本文方法在计算精度、效率以及弯矩分布等方面均展现出显著优势。     

压电材料平面裂纹尖端场的杂交应力有限元分析

基于复势理论和杂交变分原理建立了一种适用于力电耦合分析的杂交应力有限元模 型. 给出了建立刚度矩阵的主要公式和推导过程，单元内的位移场和应力场采用满足平 衡方程的复变函数级数解，假设的复变函数级数解事先精确满足裂纹的无应力和电位移法向 分量为零的条件，单元外边界的位移场假设按抛物线变化, 单元的刚度矩阵采用Gauss积分的方法得出. 通过对力电耦合裂尖场的数值计算验证了程序 的正确性和单元的有效性，同时也用所得结果校验了理论解.     

Geometric Integration Methods that Preserve Lyapunov Functions

We consider ordinary differential equations (ODEs) with a known Lyapunov function V. To ensure that a numerical integrator reflects the correct dynamical behaviour of the system, the numerical integrator should have V as a discrete Lyapunov function. Only second-order geometric integrators of this type are known for arbitrary Lyapunov functions. In this paper we describe projection-based methods of arbitrary order that preserve any given Lyapunov function. AMS subject classification (2000) 65L05, 65L06, 65L20, 65P40     

大型海湾水交换计算中随机游动方法的应用研究

水交换矩阵是反映大型海湾不同子区域之间水交换关系的有效工具,该文对随机游动方法在水交换矩阵计算中的应用进行了研究．通过与对流扩散模型的对比,指出用随机游动方法计算水交换矩阵比用对流扩散方程求解更快更灵活．同时,为了对不同区域之间的水交换特性进行长期快速预测,采用了Thompson提出的预测矩阵方法．通过理论分析得出,影响该矩阵预测效果的因素主要有:粒子数目、分区数和预测矩阵取值时间．以一个矩形海湾为例,将数值模拟结果和预测矩阵预测结果进行对比,分析了以上各因素对预测效果的影响．结果表明,粒子数目越多,预测矩阵取值时间越长,预测矩阵的预测效果越好．在粒子数目足够多的情形下,分区数越多,预测效果越好．相反,粒子数目如果较少,过多的分区数将可能导致预测矩阵效果太差而不可用．     

不同控制参数下的弹簧摆

通过拉格朗日方程得出弹簧摆系统的动力学方程,并以频率比作为控制参数,利用Matlab软件对不同控制参数和初始摆角下的弹簧摆进行了数值模拟,从而直观地研究了弹簧摆的运动.     

The influence of temperature gradients on the kinetic boundary layer problem for a condensing droplet

We extend an earlier method for solving kinetic boundary layer problems to the case of particles moving in aspatially inhomogeneous background. The method is developed for a gas mixture containing a supersaturated vapor and a light carrier gas from which a small droplet condenses. The release of heat of condensation causes a temperature difference between droplet and gas in the quasistationary state; the kinetic equation describing the vapor is the stationary Klein-Kramers equation for Brownian particles diffusing in a temperature gradient. By means of an expansion in Burnett functions, this equation is transformed into a set of coupled algebrodifferential equations. By numerical integration we construct fundamental solutions of this equation that are subsequently combined linearly to fulfill appropriate mesoscopic boundary conditions for particles leaving the droplet surface. In view of the intrinsic numerical instability of the system of equations, a novel procedure is developed to remove the admixture of fast growing solutions to the solutions of interest. The procedure is tested for a few model problems and then applied to a slightly simplified condensation problem with parameters corresponding to the condensation of mercury in a background of neon. The effects of thermal gradients and thermodiffusion on the growth rate of the droplet are small (of the order of 1%), but well outside of the margin of error of the method.