Numerical study of a diffusive epidemic model of influenza with variable transmission coefficient

A diffusive epidemic model for H1N1 influenza is formulated with a view to gain basic understanding of the virus behavior. All newborns are assumed to be susceptible. Mortality rate for infective individuals in the population is assumed to be greater than natural mortality rate. Latent, infectious and immune periods are assumed to be constants throughout this study. The numerical solutions of this model are carried out under three different initial populations distribution. In order to investigate the effect of the disease transmission coefficient on the spread of disease, β is taken to be constant as well as a function of seasonally varying time t and a function of spatial variable x  . The threshold quantity (_R0)$({\mathfrak{R}}_0)$ that governs the disease dynamics is derived. Numerical simulation shows that the system supports the existence of sustained and damped oscillations depending on initial populations distribution, the disease transmission rate and diffusion.     

Asymptotic behavior of a differential operator with discontinuities at two points

In this paper, we study a Sturm–Liouville operator with eigenparameter‐dependent boundary conditions and transmission conditions at two interior points. By establishing a new operator A associated with the problem, we prove that the operator A is self‐adjoint in an appropriate space H, discuss completeness of its eigenfunctions in H, and obtain its Green function. Copyright © 2010 John Wiley & Sons, Ltd.     

Nonstationary Stokes system in weighted Sobolev spaces

We consider an initial‐boundary value problem for nonstationary Stokes system in a bounded domain Omega??³ with slip boundary conditions. We assume that Ω is crossed by an axis L. Let us introduce the following weighted Sobolev spaces with finite norms: and where ?(x) = dist{x, L}. We proved the result. Given the external force f∈L_{2, ?µ}(Ω^T), initial velocity v₀∈H(Ω), µ∈?₊\? there exist velocity v∈H(Ω^T) and the pressure p, ?p∈L_{2, ?µ}(Ω^T) and a constant c, independent of v, p, f, such that As we consider the Stokes system in weighted Sobolev spaces the following two things must be used:

• 1. the slip boundary condition and
• 2. the Helmholtz–Weyl decomposition.

Copyright © 2010 John Wiley & Sons, Ltd.     

Homogenization of the Signorini boundary‐value problem in a thick junction and boundary integral equations for the homogenized problem

We consider a mixed boundary‐value problem for the Poisson equation in a thick junction Ω_ε which is the union of a domain Ω₀ and a large number of ε—periodically situated thin cylinders. The non‐uniform Signorini conditions are given on the lateral surfaces of the cylinders. The asymptotic analysis of this problem is done as ε→0, i.e. when the number of the thin cylinders infinitely increases and their thickness tends to zero. We prove a convergence theorem and show that the non‐uniform Signorini boundary conditions are transformed in the limiting variational inequalities in the region that is filled up by the thin cylinders as ε→0. The convergence of the energy integrals is proved as well. The existence and uniqueness of the solution to this non‐standard limit problem is established. This solution can be constructed by using a penalty formulation and successive iteration. For some subclass, these problems can be reduced to an obstacle problem in Ω₀ and an appropriate postprocessing. The equations in Ω₀ finally are also treated with boundary integral equations. Copyright © 2010 John Wiley & Sons, Ltd.     

Three‐point bending tests—Part I: Mathematical study and asymptotic analysis

The goal of this work is to study the static behaviour of a three‐dimensional elastic beam when is subjected to a three‐point bending test. In the first part, under suitable compatibility conditions on the applied forces and on the geometry of the beam, we will prove the existence of a unique solution for the associated contact elastic problem; these conditions of compatibility on the data come from the absence of a Dirichlet condition on the beam boundary. In the second part, we will study the asymptotic behaviour of this problem; in particular, we will deduce the one‐dimensional models associated with the displacement components, and we will give the existence and uniqueness of solution for them. Moreover, we will give an expression for the normal axial stress in the beam which is related to the modulus of rupture of brittle materials. In the final part of the work, we will deal with the regularity of the solution for the bending problem and we will prove some properties of the coincidence set. Copyright © 2011 John Wiley & Sons, Ltd.     

Reproducing Kernel method for the solution of nonlinear hyperbolic telegraph equation with an integral condition

In this article, an iterative method is proposed for solving nonlinear hyperbolic telegraph equation with an integral condition. Its exact solution is represented in the form of series in the reproducing kernel space. In the mean time, the n‐term approximation u_n(x, t) of the exact solution u(x, t) is obtained and is proved to converge to the exact solution. Moreover, the partial derivatives of u_n(x, t) are also convergent to the partial derivatives of u(x, t). Some numerical examples have been studied to demonstrate the accuracy of the present method. Results obtained by the method have been compared with the exact solution of each example and are found to be in good agreement with each other. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 867–886, 2011     

On Two-point Boundary Value Problems for Second-order Difference Equation

In this paper, we aim to investigate the difference equation \begin{align*} \Delta^{2}y(t-1)+|y(t)|=0, \ \ \ \ \ t\in[1,T]_{\mathbb{Z}} \end{align*} with different boundary conditions $y(0)=0$ or $\Delta y(0)=0$ and $y(T+1)=B$ or $\Delta y(T)=B$,\ where\ $T\geq 1$ is an integer and $B\in\mathbb{R}$. We will show that how the values of $T$ and $B$ influence the existence and uniqueness of the solutions to the about problem. In details, for the different problems, the $TB$-plane explicitly divided into different parts according to the number of the solutions to the above problems. These parts of $TB$-plane for the value of $T$ and $B$ guarantee the uniqueness, the existence and the nonexistence of solutions respectively.     

非对称碰撞约束下形状记忆合金梁的随机响应与分岔

形状记忆合金(SMA)是二十一世纪具有形状记忆效应的新型智能材料.针对具有非对称约束的SMA梁,本文构造了碰撞振动系统.在无碰撞和有碰撞两种情况下,利用随机平均法给出了近似解析结果.数值模拟作为验证解析结果的工具.结果表明,系统能量的概率响应曲线具有非光滑特性.当约束位置发生变化时,系统会出现随机P分岔和D分岔.     

随机结构非线性动力响应的概率密度演化分析

提出了随机结构非线性动力响应分析的概率密度演化方法．根据结构动力响应的随机状态方程，利用概率守恒原理，建立了随机结构非线性动力响应的概率密度演化方程．结合Newmark-Beta时程积分方法与Lax-Wendroff差分格式，提出了概率密度演化方程的数值分析方法．通过与Monte Carlo分析方法对比，表明所给出的概率密度演化方法具有良好的计算精度和较小的计算工作量．研究表明：随机结构非线性动力响应概率密度具有典型的演化特征，随着时间增长，概率密度曲线分布趋于复杂．