Simulation of anisotropic diffusion processes in fluids with smoothed particle hydrodynamics

Anisotropic diffusion phenomenon in fluids is simulated using smoothed particle hydrodynamics (SPH). A new SPH approximation for diffusion operator, named anisotropic SPH approximation for anisotropic diffusion (ASPHAD), is derived. Basic idea of the derivation is that anisotropic diffusion operator is first approximated by an integral in a coordinate system in which it is isotropic. The coordinate transformation is a combination of a coordinate rotation and a scaling in accordance with diffusion tensor. Then, inverse coordinate transformation and particle discretization are applied to the integral to achieve ASPHAD. Noting that weight function used in the integral approximation has anisotropic smoothing length, which becomes isotropic under the inverse transformation. ASPHAD is general and unique for both isotropic and anisotropic diffusions with either constant or variable diffusing coefficients. ASPHAD was numerically examined in some cases of isotropic and anisotropic diffusions of a contaminant in fluid, and the simulation results are very consistent with corresponding analytical solutions. Copyright © 2016 John Wiley & Sons, Ltd.     

Discontinuous Galerkin finite element discretization of a strongly anisotropic diffusion operator

The discretization of a diffusion equation with a strong anisotropy by a discontinuous Galerkin finite element method is investigated. This diffusion term is implemented in the tracer equation of an ocean model, thanks to a symmetric tensor that is composed of diapycnal and isopycnal diffusions. The strong anisotropy comes from the difference of magnitude order between both diffusions. As the ocean model uses interior penalty terms to ensure numerical stability, a new penalty factor is required in order to correctly deal with the anisotropy of this diffusion. Two penalty factors from the literature are improved and established from the coercivity property. One of them takes into account the diffusion in the direction normal to the interface between the elements. After comparison, the latter is better because the spurious numerical diffusion is weaker than with the penalty factor proposed in the literature. It is computed with a transformed coordinate system in which the diffusivity tensor is diagonal, using its eigenvalue decomposition. Furthermore, this numerical scheme is validated with the method of manufactured solutions. It is finally applied to simulate the evolution of temperature and salinity due to turbulent processes in an idealized Arctic Ocean. Copyright © 2014 John Wiley & Sons, Ltd.     

多体系统Lagrange方程数值算法的研究进展

Lagrange方法是建立多体系统动力学方程的普遍方法之一, 其方程的形式为常微分方程组或微分-代数方程组,数值计算与数 值分析是研究多体系统动力学特性的重要方法。本文简要介绍了多 体系统动 力学方程的第一、二类Lagrange方程和修正的Lagrange方 程的基本形式及这些方程的正则形式,着重介绍了正则方程在数值 计算中的特点,就多体系统Lagrange方程的隐式算法、辛算法和多 体系统动力学特性的数值分析方法(包括数值仿真、 Poincarè映射 和Lyapunov指数的计算方法)的研究现状进行了综述。     

多体系统Lagrange方程数值算法的研究进展

Lagrange方法是建立多体系统动力学方程的普遍方法之一,其方程的形式为常微分方程组或微分 - 代数方程组,数值计算与数值分析是研究多体系统动力学特性的重要方法.本文简要介绍了多体系统动力学方程的第一、二类Lagrange方程和修正的Lagrange方程的基本形式及这些方程的正则形式,着重介绍了正则方程在数值计算中的特点,就多体系统Lagrange方程的隐式算法、辛算法和多体系统动力学特性的数值分析方法(包括数值仿真、Poincar'e映射和Lyapunov指数的计算方法)的研究现状进行了综述.     

Complexity analysis of a Cournot–Bertrand duopoly game with different expectations

In this paper, we consider a Cournot–Bertrand mixed duopoly model with different expectations, where the market has linear demand and the firms have fixed marginal cost functions. Two firms choose output and price as decision variables, respectively, under the assumption that there is a certain degree of differentiation between the products offered by firms to avoid the whole market is occupied by the one that applies a lower price. The two players are considered to have bounded rational and static expectations. The existence and local stability condition of Nash equilibrium is investigated. We find the stability region of Cournot–Bertrand system is bigger than that of Cournot or Bertrand system under the same conditions. Furthermore, there are two different kinds of bifurcations when the parameters pass through the different boundary curves of the stability region, which is different from the Cournot or Bertrand model. Numerical simulation method is used to display the dynamic behaviors of the dynamical system, such as periodic cycles, bifurcation diagrams and strange attractors of the systems. The economic explanations of the complex dynamic behaviors are also given.     

竖置管流中液固两相脉动特性和颗粒浓度分布

利用激光多普勒分相测量技术,考察了液固两相自下而上通过竖置矩形管时,固、液两相的时均速度、流向及横向的脉动强度和颗粒相的相对浓度分布,证实了颗粒浓度的横向分布主要取决于颗粒的横向脉动强度分布（即npvp′2^-=常数）的分析结果。     

Asymptotic Behavior of the Principal Eigenvalue for Cooperative Elliptic Systems and Applications

The asymptotic behavior of the principal eigenvalue for general linear cooperative elliptic systems with small diffusion rates is determined. As an application, we show that if a cooperative system of ordinary differential equations has a unique positive equilibrium which is globally asymptotically stable, then the corresponding reaction-diffusion system with either the Neumann boundary condition or the Robin boundary condition also has a unique positive steady state which is globally asymptotically stable, provided that the diffusion coefficients are sufficiently small. Moreover, as the diffusion coefficients approach zero, the positive steady state of the reaction-diffusion system converges uniformly to the equilibrium of the corresponding kinetic system.     

一类混沌系统中的簇发振荡及其延迟叉形分岔行为

由于多时间尺度问题在实际工程系统中广泛存在,关于其复杂动力学行为及其产生机制的研究已成为当前国内外的热点课题之一.簇发振荡是多时间尺度系统复杂动力学行为的典型代表,而分岔延迟又是簇发振荡中的常见现象.本文为探讨非线性系统中分岔延迟所引发的簇发振荡的分岔机制,在一个三维混沌系统中引入参数激励,当激励频率远小于系统的固有频率时,系统产生了两时间尺度簇发振荡.将整个激励项看做慢变参数,激励系统转化为广义自治系统也即快子系统,分析快子系统平衡点的稳定性以及分岔条件,并运用快慢分析法和转换相图揭示了簇发振荡的动力学机理.文中考察了4组参数条件下系统的动力学行为,研究发现当慢变激励项周期性地通过分岔点时,系统产生了明显的超临界叉形分岔延迟行为,随着参数激励振幅的增大,分岔延迟的时间也逐渐延长,当这种延迟的动态行为终止于不同的参数区域时,导致系统轨线围绕不同稳定吸引子(平衡点,极限环)运动,从而得到了不同的簇发振荡行为.      

Chaotic attitude motion of gyrostat statellites in a central force field

The nonlinear attitude motion of gyrostat satellites in a central force field is investigated, with particular emphasis on their long-time dynamic behavior for a wide range of parameters. The numbers of equilibrium solutions, as well as their stability, vary with the rotor speed, and bifurcation diagrams have been obtained. Various dynamic behaviors of gyrostat satcllites, e.g. periodic, quasiperiodic, and chaotic, are studied via the Poincaré map technique. It is shown that the rotor speed has a significant effect on the dynamic behavior of gyrostat satellites.     

COMPLICATED BURSTING BEHAVIORS AS WELL AS THE MECHANISM OF A THREE DIMENSIONAL NONLINEAR SYSTEM 1)

由多时间尺度耦合效应引起的簇发振荡行为是非线性动力学研究的重要课题之一.本文针对一类参数激励下的三维非线性电机系统(该系统可以描述两种自激同极发电机系统的动力学行为,两种系统在数学上等效),研究了当参数激励频率远小于系统自然频率时的各种复杂簇发振荡行为及其产生机理.通过快慢分析方法, 将参数激励作为慢变参数,得到了非自治系统对应的广义自治系统及快子系统和慢变量,并给出了快子系统的稳定性和分岔条件以及系统关于典型参数的单参数分岔图.借助转换相图与分岔图的叠加, 分析了对称式delayed subHopf/fold cycle簇发振荡的产生机理及其动力学转迁, 即delayed subHopf/fold cycle簇发振荡、焦点/焦点型对称式叉形分岔滞后簇发振荡和焦点/焦点型叉形分岔滞后簇发振荡.研究结果表明, 系统会出现两种不同的分岔滞后形式, 一种是亚临界Hopf分岔滞后,另一种是叉形分岔滞后,而且控制参数显著影响平衡点的稳定性和分岔滞后区间的宽度.同时初始点的选取则会影响系统动力学行为的对称性.本文的研究进一步加深了对由分岔滞后引起的簇发振荡的认识和理解.     

Attractivity of Equilibrium Sets of Systems with Dry Friction

The dynamics of mechanical systems with dry friction elements, modelledby set-valued force laws, can be described by differential inclusions.An equilibrium set of such a differential inclusion corresponds to astationary mode for which the friction elements are sticking. Theattractivity properties of the equilibrium set are of major importancefor the overall dynamic behaviour of this type of systems. Conditionsfor the attractivity of the equilibrium set of MDOF mechanical systemswith multiple friction elements are presented. These results areobtained by application of a generalisation of LaSalle's principle fordifferential inclusions of Filippov-type. Besides passive systems, alsosystems with negative viscous damping are considered. For such systems,only local attractivity of the equilibrium set can be assured undercertain conditions. Moreover, an estimate for the region of attractionis given for these cases. The effectiveness of the results isillustrated by means of both 1DOF and MDOF examples.     

一类三维非线性系统的复杂簇发振荡行为及其机理

由多时间尺度耦合效应引起的簇发振荡行为是非线性动力学研究的重要课题之一.本文针对一类参数激励下的三维非线性电机系统(该系统可以描述两种自激同极发电机系统的动力学行为,两种系统在数学上等效),研究了当参数激励频率远小于系统自然频率时的各种复杂簇发振荡行为及其产生机理.通过快慢分析方法, 将参数激励作为慢变参数,得到了非自治系统对应的广义自治系统及快子系统和慢变量,并给出了快子系统的稳定性和分岔条件以及系统关于典型参数的单参数分岔图.借助转换相图与分岔图的叠加, 分析了对称式delayed subHopf/fold cycle簇发振荡的产生机理及其动力学转迁, 即delayed subHopf/fold cycle簇发振荡、焦点/焦点型对称式叉形分岔滞后簇发振荡和焦点/焦点型叉形分岔滞后簇发振荡.研究结果表明, 系统会出现两种不同的分岔滞后形式, 一种是亚临界Hopf分岔滞后,另一种是叉形分岔滞后,而且控制参数显著影响平衡点的稳定性和分岔滞后区间的宽度.同时初始点的选取则会影响系统动力学行为的对称性.本文的研究进一步加深了对由分岔滞后引起的簇发振荡的认识和理解.      

Global analysis of Ivlev's type predator-prey dynamic systems

Consider a class of Ivlev's type predator-prey dynamic systems with prey and predator both having linear density restricts. By using the qualitative methods of ODE, the global stability of positive equilibrium and existence and uniqueness of non-small amplitude stable limit cycle are obtained. Especially under certain conditions, it shows that existence and uniqueness of non-small amplitude stable limit cycle is equivalent to the local un-stability of positive equilibrium and the local stability of positive equilibrium implies its global stability. That is to say, the global dynamic of the system is entirely determined by the local stability of the positive equilibrium.     

Kinematic simulation for stably stratified and rotating turbulence

The properties of one-particle and particle-pair diffusion in rotating and stratified turbulence are studied by applying the rapid distortion theory (RDT) to a kinematic simulation (KS) of the Boussinesq equation with a Coriolis term.Scalings for one- and two-particle horizontal and vertical diffusions in purely rotating turbulence are proposed for small Rossby numbers.Particular attention is given to the locality-in-scale hypothesis for two-particle diffusion in purely rotating turbulence both in the horizontal and the vertical directions. It is observed that both rotation and stratification decrease the pair diffusivity and improve the validity of the locality-in-scale hypothesis. In the case of stratification the range of scales over which the locality-in-scale hypothesis is observed is increased.It is found that rotation decreases the diffusion in the horizontal direction as well as, though to a much lesser extent, in the vertical direction.     

Global analysis of Ivlev＇s type predator-prey dynamic systems

Consider a class of Ivlev＇s type predator-prey dynamic systems with prey and predator both having linear density restricts. By using the qualitative methods of ODE, the global stability of positive equilibrium and existence and uniqueness of non-small amplitude stable limit cycle are obtained. Especially under certain conditions, it shows that existence and uniqueness of non-small amplitude stable limit cycle is equivalent to the local un-stability of positive equilibrium and the local stability of positive equilibrium implies its global stability. That is to say, the global dynamic of the system is entirely determined by the local stability of the positive equilibrium.     

Constructing a chaotic system with any number of equilibria

In the chaotic Lorenz system, Chen system and Rössler system, their equilibria are unstable and the number of the equilibria are no more than three. This paper shows how to construct some simple chaotic systems that can have any preassigned number of equilibria. First, a chaotic system with no equilibrium is presented and discussed. Then a methodology is presented by adding symmetry to a new chaotic system with only one stable equilibrium, to show that chaotic systems with any preassigned number of equilibria can be generated. By adjusting the only parameter in these systems, one can further control the stability of their equilibria. This result reveals an intrinsic relationship of the global dynamical behaviors with the number and stability of the equilibria of a chaotic system.     

Analytical and computational studies on the vacuum performance of a chevron ejector

The effects of chevrons on the performance of a supersonic vacuum ejector-diffuser system are investigated numerically and evaluated theoretically in this work. A three-dimensional geometrical domain is numerically solved using a fully implicit finite volume scheme based on the unsteady Reynolds stress model. A one-dimensional mathematical model provides a useful tool to reveal the steady flow physics inside the vacuum ejector-diffuser system. The effects of the chevron nozzle on the generation of recirculation regions and Reynolds stress behaviors are studied and compared with those of a conventional convergent nozzle. The present performance parameters obtained from the simulated results and the mathematical results are validated with existing experimental data and show good agreement. Primary results show that the duration of the transient period and the secondary chamber pressure at a dynamic equilibrium state depend strongly on the primary jet conditions, such as inlet pressure and primary nozzle shape. Complicated oscillatory flow, generated by the unsteady movement of recirculation, finally settles into a dynamic equilibrium state. As a vortex generator, the chevron demonstrated its strong entrainment capacity to accelerate the starting transient flows to a certain extent and reduce the dynamic equilibrium pressure of the secondary chamber significantly.     

Thickness Effect of Pulse Shaper on Dynamic Stress Equilibrium and Dynamic Deformation Behavior in the Polycarbonate Using SHPB Technique

0Introduction Thehighstrainratestress strainresponsesofpolymersandpolymericcompositematerialshave receivedincreasedscientificandindustrialattentioninrecentyears.Polymericmaterialsaresubjected todynamicloadingandhighstrainratedeformationinavarietyofimporta…     

对流扩散方程的摄动有限体积(PFV)方法及讨论

在有限体积(FV)方法的重构近似中,引入数值摄动处理,即把界面数值通量摄动展开成网格间距的幂级数,并利用积分方程自身的性质求出幂级数的系数,同时获得高精度迎风和中心型摄动有限体积(PFV)格式.对标量输运方程给出积分近似为二阶、重构近似为二、三和四阶迎风和中心型PFV格式,这些PFV格式的结构形式及使用基点数与一阶迎风格式完全一致,迎风PFV格式满足对流有界准则;二阶和四阶中心PFV格式对网格Peclet数的任意值均为正型格式,比常用的二阶中心格式优越.用一维标量输运和方腔流动算例说明PFV格式的优良性能,并把PFV方法与性质相近的摄动有限差分(PFD)方法及相关的高精度方法作了对比分析.     

Characterizing dynamic load propagation in cohesionless granular packing using force chain

When dynamic load is applied on a granular assembly, the time-dependent dynamic load and initial static load (such as gravity stress) act together on individual particles. In order to better understand how dynamic load triggers the micro-structure's evolution and furtherly the ensemble behavior of a granular assembly, we propose a criterion to recognize the major propagation path of dynamic load in 2D granular materials, called the “dynamic force chain”. Two steps are involved in recognizing dynamic force chains: (1) pick out particles with dynamic load larger than the threshold stress, where the attenuation of dynamic stress with distance is considered; (2) among which quasi-linear arrangement of three or more particles are identified as a force chain. The spatial distribution of dynamic force chains in indentation of granular materials provides a direct measure of dynamic load diffusion. The statistical evolution of dynamic force chains shows strong correlation with the indentation behaviors.