Neural firing rate model with a steep firing rate function

In this study we justify rigorously the approximation of the steep firing rate functions with a unit step function in a two-population neural firing rate model with steep firing rate functions. We do this justification by exploiting the theory of switching dynamical systems. It has been demonstrated that switching dynamics offer a possibility of simplifying the dynamical system and getting approximations of the solution of the system for any specific choice of parameters. In this approach the phase space of the system is divided into regular and singular domains, where the limit dynamics can be written down explicitly. The advantages of this method are illustrated by a number of numerical examples for different cases of the singular domains (i.e. for black, white and transparent walls) and for specific choices of the parameters involved. General conditions have been formulated on these parameters to give black, white and transparent walls. Further, the existence and stability of regular and singular stationary points have been investigated. It has been shown that the regular stationary points (i.e. stationary points inside the regular domains) are always stable and this property is preserved for smooth and sufficiently steep activation functions. In the most technical part of the paper we have provided conditions on the existence and stability of singular stationary points (i.e. those belonging to the singular domains). For the existence results, the implicit function theorem has been used, whereas the stability of singular stationary points is addressed by applying singular perturbation analysis and the Tikhonov theorem.     

Stationary availability of a semi‐Markov system with random maintenance

We consider a reparable system with a finite state space, evolving in time according to a semi‐Markov process. The system is stopped for it to be preventively maintained at random times for a random duration. Our aim is to find the preventive maintenance policy that optimizes the stationary availability, whenever it exists. The computation of the stationary availability is based on the fact that the above maintained system evolves according to a semi‐regenerative process. As for the optimization, we observe on numerical examples that it is possible to limit the study to the maintenance actions that begin at deterministic times. We demonstrate this result in a particular case and we study the deterministic maintenance policies in that case. In particular, we show that, if the initial system has an increasing failure rate, the maintenance actions improve the stationary availability if and only if they are not too long on the average, compared to the repairs ( a bound for the mean duration of the maintenance actions is provided). On the contrary, if the initial system has a decreasing failure rate, the maintenance policy lowers the stationary availability. A few other cases are studied. Copyright © 2000 John Wiley & Sons, Ltd.     

Wave simulations of Gray‐Scott reaction‐diffusion system

In this work, we study the numerical simulation of the one‐dimensional reaction‐diffusion system known as the Gray‐Scott model. This model is responsible for the spatial pattern formation, which we often meet in nature as the result of some chemical reactions. We have used the trigonometric quartic B‐spline (T4B) functions for space discretization with the Crank‐Nicolson method for time integration to integrate the nonlinear reaction‐diffusion equation into a system of algebraic equations. The solutions of the Gray‐Scott model are presented with different wave simulations. Test problems are chosen from the literature to illustrate the stationary waves, pulse‐splitting waves, and self‐replicating waves.     

Comparison of the bifurcation curves of a two-variable and a three-variable circadian rhythm model

Two dynamical systems describing the circadian fluctuation of two proteins (PER and TIM) in cells are compared. A simplified model with two variables has already been investigated. Detailed study of the possible bifurcation has been carried out. Periodic solutions of the differential equations with 24-h period have been obtained numerically. Here the general, more realistic model having three variables is investigated. The possible phase portraits and local bifurcations are studied in detail. The saddle-node and Hopf-bifurcation curves are determined in the plane of two parameters by using the parametric representation method. Using these curves the number and the type of the stationary points can be determined. The relation of the two bifurcation curves and the Takens–Bogdanov bifurcation points are also studied. The bifurcation curves are compared to those obtained for the simplified two-variable system.     

Superconvergence of a stabilized finite element approximation for the Stokes equations using a local coarse mesh L2 projection

This article first recalls the results of a stabilized finite element method based on a local Gauss integration method for the stationary Stokes equations approximated by low equal‐order elements that do not satisfy the inf‐sup condition. Then, we derive general superconvergence results for this stabilized method by using a local coarse mesh L² projection. These supervergence results have three prominent features. First, they are based on a multiscale method defined for any quasi‐uniform mesh. Second, they are derived on the basis of a large sparse, symmetric positive‐definite system of linear equations for the solution of the stationary Stokes problem. Third, the finite elements used fail to satisfy the inf‐sup condition. This article combines the merits of the new stabilized method with that of the L² projection method. This projection method is of practical importance in scientific computation. Finally, a series of numerical experiments are presented to check the theoretical results obtained. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 28: 115‐126, 2012     

M/C2/c/K queuing model and optimization for a geriatric care center

Queuing systems with finite buffers are reasonable models for many manufacturing, telecommunication, and healthcare systems. Although some approximations exist, the exact analysis of multi‐server and finite‐buffer queues with general service time distribution is unknown. However, the phase‐type assumption for service time is a frequently used approach. Because the Cox distribution, a kind of phase‐type distribution, provides a good representation of data with great variability, it has a vast area of application in modeling service times. The research focus is twofold. First, a theoretical structure of a multi‐server and finite‐buffer queuing system in which the service time is modeled by the two‐phase Cox distribution is studied. It is focused on finding an efficient solution to the stationary probabilities using the matrix‐geometric method. It is shown that the stationary probability vector can be obtained with the matrix‐geometric method by using level‐dependent rate matrices, and the mean queue length is computed. Second, an empirical analysis of the model is presented. The proposed methodology is applied in a case study concerning the geriatric patients. Some numerical calculations and optimizations are performed by using geriatric data. Copyright © 2015 John Wiley & Sons, Ltd.     

带刚性限位的双层隔振系统的离散随机模型

研究由多刚体组成的带刚性限位的双层隔振系统,对其冲击后受到周期性外激励和低强度噪声扰动共同作用下可能会产生的碰撞进行了分析．基于单向约束多体动力学理论,导出了此隔振系统的最大Poincaré映射,建立了其冲击后的零次近似随机离散模型和一次近似随机离散模型．通过对一MTU公司的柴油机隔振系统冲击作用后振动响应的调查指出,由于可能发生间歇性碰撞,该系统呈现复杂的非线性特性．零次近似模型和一次近似模型有较大的区别,低强度的噪声也会对系统产生较大的影响．得到的结果对如何正确设计带刚性限位的双层隔振系统提供了理论参考依据．     

On Multi‐component Ermakov Systems in a Two‐Layer Fluid: Integrable Hamiltonian Structures and Exact Vortex Solutions

By introducing an elliptic vortex ansatz, the 2+1‐dimensional two‐layer fluid system is reduced to a finite‐dimensional nonlinear dynamical system. Time‐modulated variables are then introduced and multicomponent Ermakov systems are isolated. The latter is shown to be also Hamiltonian, thereby admitting general solutions in terms of an elliptic integral representation. In particular, a subclass of vortex solutions is obtained and their behaviors are simulated. Such solutions have recently found applications in oceanic and atmospheric dynamics. Moreover, it is proved that the Hamiltonian system is equivalent to the stationary nonlinear cubic Schrödinger equations coupled with a Steen‐Ermakov‐Pinney equation.     

Stationary solutions to the Falk system on shape memory alloys

We study the Falk model system describing martensitic phase transitions in shape memory alloys. Its physically closed stationary state is formulated as a nonlinear eigenvalue problem with a non‐local term. Then, some results on existence, stability, and bifurcation of the solution are proven. In particular, we prove the existence of dynamically stable nontrivial stationary solutions. Copyright © 2009 John Wiley & Sons, Ltd.     

具有标准发生率的随机SIR传染病模型的建立及平稳分布与疾病灭绝研究

研究了一类具有标准发生率以及考虑随机扰动与系统变量成正比的随机SIR传染病模型.首先,对于任意的正的初值,系统存在唯一的全局正解以及通过构造合适的随机李雅普诺夫函数,得到了模型遍历平稳分布存在的充分条件.其次,给出了疾病灭绝的充分条件,并与模型遍历平稳分布存在的充分条件作对比,得出了在特定条件下随机SIR模型的阈值.最后通过数值模拟验证了结果的正确性.     

An analytical and numerical study of liquid dynamics in a one‐dimensional capillary under entrapped gas action

Capillary dynamics has been and is yet an important field of research, because of its very relevant role played as the core mechanism at the base of many applications. In this context, we are particularly interested in the liquid penetration inspection technique. Because of the obviously needed level of reliability involved with such a non‐destructive test, this paper is devoted to study how the presence of an entrapped gas in a close‐end capillary may affect the inspection outcome. This study is carried out through a one‐dimensional ordinary differential model that despite its simplicity is able to point out quite well the capillary dynamics under the effect of an entrapped gas. The paper is divided into two main parts; the first starts from an introductory historical review of capillary flows modelling, goes on presenting the one‐dimensional second order ordinary differential model, taking into account the presence of an entrapped gas and therefore ends by showing some numerical simulation results. The second part is devoted to the analytical study of the model by separating the initial transitory behaviour from the stationary one. Besides, these solutions are compared with the numerical ones, and finally, an expression is deduced for the threshold radius switching from a fully damped transitory to an oscillatory one. Copyright © 2013 John Wiley & Sons, Ltd.     

Ergodicity and threshold behaviors of a predator‐prey model in stochastic chemostat driven by regime switching

This paper deals with a stochastic predator‐prey model in chemostat which is driven by Markov regime switching. For the asymptotic behaviors of this stochastic system, we establish the sufficient conditions for the existence of the stationary distribution. Then, we investigate, respectively, the extinction of the prey and predator populations. We explore the new critical numbers between survival and extinction for species of the dual‐threshold chemostat model. Numerical simulations are accomplished to confirm our analytical conclusions.     

Transformation of Time-Delay Systems to a Form with Zero Dynamics

The problem of transforming a controlled linear stationary system of differential equations with commensurable time delays into a canonical form with zero dynamics is considered. This problem has been well studied for ODE systems and is closely related to the concept of a relative degree. In this paper, the results are extended to time-delay systems.     

Analysis of a PDE model of the swelling of mitochondria accounting for spatial movement

We analyze existence and asymptotic behavior of a system of semilinear diffusion‐reaction equations that arises in the modeling of the mitochondrial swelling process. The model itself expands previous work in which the mitochondria were assumed to be stationary, whereas now their movement is modeled by linear diffusion. While in the previous model certain formal structural conditions were required for the rate functions describing the swelling process, we show that these are not required in the extended model. Numerical simulations are included to visualize the solutions of the new model and to compare them with the solutions of the previous model.     

Mathematical and Numerical Study of the 3D Atmospheric Boundary Layer Flows Including Pollution Dispersion

The paper presents a mathematical and numerical investigation of the atmospheric boundary layer (ABL) flow over 3D complex terrain part of which is represented by the real topography of the Krkonose mountains located in the Czech Republic. The flow is supposed to be turbulent, non‐stratified, viscous, incompressible and stationary. Two mathematical models have been formulated. The first model is based upon the RANS equations in the conservative form and the second one uses the Boussinesq approximation of RANS equations and takes the non‐conservative form. Also pollution dispersion over the complex 3D terrain has been considered in both models. The problem closure is achieved by an algebraic turbulence model and given boundary conditions.     

Asymptotic dynamics of the one‐dimensional attraction–repulsion Keller–Segel model

The asymptotic behavior of the attraction–repulsion Keller–Segel model in one dimension is studied in this paper. The global existence of classical solutions and nonconstant stationary solutions of the attraction–repulsion Keller–Segel model in one dimension were previously established by Liu and Wang (2012), which, however, only provided a time‐dependent bound for solutions. In this paper, we improve the results of Liu and Wang (2012) by deriving a uniform‐in‐time bound for solutions and furthermore prove that the model possesses a global attractor. For a special case where the attractive and repulsive chemical signals have the same degradation rate, we show that the solution converges to a stationary solution algebraically as time tends to infinity if the attraction dominates. Copyright © 2014 John Wiley & Sons, Ltd.     

具有非抢占优先站点的轮询系统顾客队长

本文讨论轮询系统在系统平稳条件下,对于具有一个非抢占的优先权站点且采用穷尽服务方式下的轮询系统进行理论分析,给出服务员轮询到每个站点时该站点的队长及数学期望。     

Confined steady states of a Vlasov‐Poisson plasma in an infinitely long cylinder

We consider the two‐dimensional Vlasov‐Poisson system to model a two‐component plasma whose distribution function is constant with respect to the third space dimension. First, we show how this two‐dimensional Vlasov‐Poisson system can be derived from the full three‐dimensional model. The existence of compactly supported steady states with vanishing electric potential in a three‐dimensional setting has already been investigated in the literature. We show that these results can easily be adapted to the two‐dimensional system. However, our main result is to prove the existence of compactly supported steady states even with a nontrivial self‐consistent electric potential.     

Nontrivial stationary points of two-species self-structured communities

The two-species model of self-structured stationary biological communities proposed by U. Dieckmann and R. Law is considered. A way of investigating the system of integro-differential equations describing the model equilibrium is developed, nontrivial stationary points are found, and constraints on the model parameter space resulting in similar stationary points are studied. The results are applied to a number of widely known biological scenarios.     

索赔次数为复合Poisson-Geometric过程的风险模型及破产概率

本文引入一类复合Poisson-Geometric分布,这类分布包括两个参数,是普通Poisson分布的一种推广,并在保险中有其实际的应用背景;基于此分布产生一个计数过程,称之为复合Poisson-Geometric过程.本文着重研究了索赔次数为复合Poisson-Geometric过程的风险模型,这种模型是经典风险模型的一个推广.针对此模型,本文给出了破产概率公式及更新方程.作为特例,当索赔额服从指数分布时,给出了破产概率的显式表达式.