Step size strategies for the numerical integration of systems of differential equations

In this study, the step size strategies are obtained such that the local error is smaller than the desired error level in the numerical integration of a type of nonlinear equation system in interval [_t0,T]$[t_0,T]$. The algorithms are given for calculating step sizes and numerical solutions according to these strategies. The algorithms are supported by the numerical examples.     

A criterion for the strong ellipticity of the equilibrium equations of an isotropic non-linearly elastic material

A new criterion for satisfying the condition of strong ellipticity of the equilibrium equations of a homogeneous isotropic elastic medium is obtained, consisting of twelve unconditional and the same number of conditional inequalities. The criterion differs from previously known criteria in the fact that it is more effective for the analytical verification of the condition of strong ellipticity. A number of sufficient criteria of the strong ellipticity of the equations of the non-linear statics of an elastic body is presented. It is shown that the method used in this paper is also effective in the case of an incompressible material. A new constitutive inequality in the non-linear theory of elasticity is proposed that, in particular, ensures the strong ellipticity of the equilibrium equations and the stability of an elastic body in a state of hydrostatic compression. A mechanical interpretation of some implications of the above mentioned inequality is given. Using the criterion obtained, the condition of strong ellipticity is verified for a number of specific models of compressible elastic materials and examples of the violation of this condition are also considered.     

On the rate of convergence to equilibrium in the Becker-Döring equations

We provide an explicit rate of convergence to equilibrium for solutions of the Becker-Döring equations using the energy/energy-dissipation relation. The main difficulty is the structure of equilibria of the Becker-Döring equations, which do not correspond to a Gaussian measure, such that a logarithmic Sobolev-inequality is not available. We prove a weaker inequality which still implies for fast decaying data that the solution converges to equilibrium as e^−ct1/3.     

具有垂直传染和接触传染的传染病模型的稳定性研究

本文研究了一类具有垂直传染和接触传染的传染病模型.利用常微分方程定性与稳定性方法,分析了该模型非负平衡点的存在性及其局部稳定性.同时,利用LaSalle不变性原理和通过构造适当的Lyapunov函数,获得了平凡平衡点、无病平衡点和地方病平衡点全局渐近稳定的充分条件.结果表明当基本再生数小于等于1时,所有种群趋于灭绝;当基本再生数大于1和病毒主导再生数小于1时,病毒很快被清除;当基本再生数大于1和病毒主导再生数大于1以及满足一定条件时,病毒持续流行并将成为一种地方病.     

Linearized stability for nonlinear evolution equations

We present a general principle of linearized stability at an equilibrium point for the Cauchy problem , for an -accretive, possibly multivalued, operator in a Banach space , that has a linear 'resolvent-derivative' . The result is applied to derive linearized stability results for the case of under 'minimal' differentiability assumptions on the operators and at the equilibrium point, as well as for partial differential delay equations. RID="h1" ID="h1"To the memory of Philippe Bénilan     

Existence of solutions for the equations

We introduce a concept of weak solution for a boundary value problem modelling the interactive motion of a coupled system consisting in a rigid body immersed in a viscous fluid. The fluid, and the solid are contained in a fixed open bounded set of R3. The motion of the fluid is governed by the incompresible Navier-Stokes equations and the standard conservation's laws of linear, and angular momentum rules the dynamics of the rigid body. The time variation of the fluid's domain (due to the motion of the rigid body) is not known apriori, so we deal with a free boundary value problem. Our main theorem asserts the existence of at least one weak solution for this problem. The result is global in time provided that the rigid body does not touch the boundary     

Two-grid methods for finite volume element approximations of nonlinear parabolic equations

Two-grid methods are studied for solving a two dimensional nonlinear parabolic equation using finite volume element method. The methods are based on one coarse-grid space and one fine-grid space. The nonsymmetric and nonlinear iterations are only executed on the coarse grid and the fine-grid solution can be obtained in a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. The two-grid methods achieve asymptotically optimal approximation as long as the mesh sizes satisfy h=O(H³|lnH|)$h=O\left(H^3|lnH|\right)$. As a result, solving such a large class of nonlinear parabolic equations will not be much more difficult than solving one single linearized equation.     

Global stability of a virus infection model with two delays and two types of target cells

In this paper, incorporating the delay of viral cytopathicity within target cells, we first presented a basic model of viral infection with delay, and then extended it into a model with two delays and two types of target cells. For the models proposed here, both their basic reproduction numbers are found. By constructing Lyapunov functionals, necessary and sufficient conditions ensuring the global stability of the models with delays are given. The obtained results show that, when the basic reproduction number is not greater than one, the infection-free equilibrium is globally stable in the feasible region, which implies that the virus infection goes extinct eventually; when it is greater than one, the infection equilibrium is globally stable in the feasible region, which implies that the virus infection persists in the body of host.     

An iterative substructuring algorithm for equilibrium equations

Summary The topic of iterative substructuring methods, and more generally domain decomposition methods, has been extensively studied over the past few years, and the topic is well advanced with respect to first and second order elliptic problems. However, relatively little work has been done on more general constrained least squares problems (or equivalent formulations) involving equilibrium equations such as those arising, for example, in realistic structural analysis applications. The potential is good for effective use of iterative algorithms on these problems, but such methods are still far from being competitive with direct methods in industrial codes. The purpose of this paper is to investigate an order reducing, preconditioned conjugate gradient method proposed by Barlow, Nichols and Plemmons for solving problems of this type. The relationships between this method and nullspace methods, such as the force method for structures and the dual variable method for fluids, are examined. Convergence properties are discussed in relation to recent optimality results for Varga's theory ofp-cyclic SOR. We suggest a mixed approach for solving equilibrium equations, consisting of both direct reduction in the substructures and the conjugate gradient iterative algorithm to complete the computations.Dedicated to R. S. Varga on the occasion of his 60th birthdayResearch completed while pursuing graduate studies sponsored by the Department of Mathematical Sciences, US Air Force Academy, CO, and funded by the Air Force Institute of Technology, WPAFB, OHResearch supported by the Air Force under grant no. AFOSR-88-0285 and by the National Science Foundation under grant no. DMS-89-02121     

Two-grid methods for characteristic finite volume element solution of semilinear convection-diffusion equations

Two-grid methods for characteristic finite volume element solutions are presented for a kind of semilinear convection-dominated diffusion equations. The methods are based on the method of characteristics, two-grid method and the finite volume element method. The nonsymmetric and nonlinear iterations are only executed on the coarse grid (with grid size H). And the fine-grid solution (with grid size h) can be obtained by a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. The two-grid methods achieve asymptotically optimal approximation as long as the mesh sizes satisfy H = O(h^1/3).     

具有免疫应答和吸收效应的病毒感染模型分析

研究了具有免疫应答和吸收效应的病毒动力学模型的动力学行为.通过构造适当的Lyapunov泛函,使用LaSalle不变性原理,证明了基本再生数、CTL免疫再生数、抗体免疫再生数、CTL免疫竞争再生数和抗体免疫竞争再生数决定了模型的全局性态.若基本再生数小于等于1,病毒在体内清除.若基本再生数大于1,正解在满足条件max{...     

Solutions of sturm-liouville boundary value problems for higher-order differential equations

The existence of solutions of a class of two-point boundary value problems for higher order differential equations is studied. Sufficient conditions for the existence of at least one solution are established. It is of interest that the nonlinearityf in the equation depends on all lower derivatives, and the growth conditions imposed onf are allowed to be super-linear (the degrees of phases variables are allowed to be greater than 1 if it is a polynomial). The results are different from known ones since we don’t apply the Green’s functions of the corresponding problem and the method to obtain a priori bound of solutions are different enough from known ones. Examples that can not be solved by known results are given to illustrate our theorems.     

Boundary controllability for the semilinear Schrödinger equations on Riemannian manifolds

We study the boundary exact controllability for the semilinear Schrödinger equation defined on an open, bounded, connected set Ω of a complete, n-dimensional, Riemannian manifold M with metric g. We prove the locally exact controllability around the equilibria under some checkable geometrical conditions. Our results show that exact controllability is geometrical characters of a Riemannian metric, given by the coefficients and equilibria of the semilinear Schrödinger equation. We then establish the globally exact controllability in such a way that the state of the semilinear Schrödinger equation moves from an equilibrium in one location to an equilibrium in another location.     

The stability of generalized Volterra equations

In this paper, a particular type of a system of generalized Volterra equations [1], whose solutions are assured to be nonnegative for arbitrary nonnegative initial values, is considered. The extended stability theorem of LaSalle is used for deriving conditions for a nonnegative equilibrium point to be stable with respect to a certain subset of the Euclidean space. The obtained stability theorem has a close relation with Lyapunov's stability condition for linear systems with constant coefficients and is generally less restrictive than conditions known so far.     

Integral invariants of the Hamilton equations

Conditions are found for the existence of integral invariants of Hamiltonian systems. For two-degrees-of-freedom systems these conditions are intimately related to the existence of nontrivial symmetry fields and multivalued integrals. Any integral invariant of a geodesic flow on an analytic surface of genus greater than 1 is shown to be a constant multiple of the Poincaré-Cartan invariant. Poincaré's conjecture that there are no additional integral invariants in the restricted three-body problem is proved. Translated fromMatematicheskie Zametki, Vol. 58, No. 3, pp. 379–393, September, 1995. The work was financially supported by the Russian Foundation for Basic Research (grant No. 242 93-013-16244), International Science Foundation (grant No. MCY 000), and INTAS (grant No. 93-339).     

Characterization of equilibrium measures for critical reversible Nearest Particle Systems

We show that for critical reversible attractive Nearest Particle Systems all equilibrium measures are convex combinations of the upper invariant equilibrium measure and the point mass at all zeros, provided the underlying renewal sequence possesses moments of order strictly greater than and obeys some natural regularity conditions. The research of this author is partially supported by the FNS Grant #200021-107475/1     

Boundary controllability for the quasi-linear wave equations coupled in parallel

We study the boundary exact controllability for a system of two quasi-linear wave equations coupled in parallel with springs and viscous terms. We prove the locally exact controllability around superposition equilibria under some checkable geometrical conditions. We then establish the globally exact controllability in such a way that the state of the coupled quasi-linear system moves from a superposition equilibrium in one location to a superposition equilibrium in another location. Our results show that exact controllability is geometrical characters of a Riemannian metric, given by the coefficients and superposition equilibria of the system.     

Least-squares finite element methods for the Navier-Stokes equations for generalized Newtonian fluids

In this contribution we present the least-squares finite element method (LSFEM) for the incompressible Navier-Stokes equations. In detail, we consider a non-Newtonian fluid flow, which is described by a power-law model, see [1]. The second-order problem is reformulated by introducing a first-order div-grad system consisting of the equilibrium condition, the incompressibility condition and the constitutive equation, which are written in residual forms, see [2]. Here, higher-order finite elements which are an important aspect regarding accuracy for the present formulation are investigated. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence to equilibria for compressible navier-stokes equations with large data

We show that for any sequence t_n→∞ and any global weak solution (ρ(t), u(t)) of barotropic compressible Navier-Stokes equations in 2 and 3 space dimensions with potential force there exists a subsequence {s_n} such that ϱ(s_n) in L _ω ^r (Ω) of spatially periodic functions, where ϱ_∞ is an equilibrium density, r>1 suitable. If the equilibrium is unique then the convergence holds for all t→∞. No smallness of data is assumed. Entrata in Redazione il 22 febbraio 1999. Ricevuta nuova versione il 24 agosto 1999. This work was started and completed during the stay of the second author at the University of Toulon in March 1997 and April 1998, which is acknowledged for the financial support, and also partially supported by the Grant Agency of the Czech Republic, Grant No. 201/96/0313.