Spatio-temporal dynamics of a reaction–diffusion–advection food-limited system with nonlocal delayed competition and Dirichlet boundary condition

Spatio-temporal dynamics of a reaction–diffusion–advection food-limited population model with nonlocal delayed competition and Dirichlet boundary condition are considered. Existence and stability of the positive spatially nonhomogeneous steady state solution are shown. Existence and direction of the spatially nonhomogeneous steady-state-Hopf bifurcation are proved. Stable spatio-temporal patterns near the steady-state-Hopf bifurcation point are numerically obtained. We also investigate the joint influences of some important parameters including advection rate, food-limited parameter and nonlocal delayed competition on the dynamics. It is found that the effect of advection on Hopf bifurcation is opposite with the corresponding no-flux system. The theoretical results provide some interesting highlights in ecological protection in streams or rivers.     

Chaotic Dynamics of a Delayed Discrete-Time Hopfield Network of Two Nonidentical Neurons with no Self-Connections

A bifurcation analysis is undertaken for a discrete-time Hopfield neural network of two neurons with two delays, two internal decays and no self-connections, choosing the product of the interconnection coefficients as the characteristic parameter for the system. The stability domain of the null solution is found, the values of the characteristic parameter for which bifurcations occur at the origin are identified, and the existence of Fold/Cusp, Neimark–Sacker and Flip bifurcations is proved. All these bifurcations are analyzed by applying the center manifold theorem and the normal form theory. It is shown that the dynamics in a neighborhood of the null solution become more and more complex as the characteristic parameter grows in magnitude and passes through the bifurcation values. Under certain conditions, it is proved that if the magnitudes of the interconnection coefficients are large enough, the neural network exhibits Marotto’s chaotic behavior.      

Algebraic decay and variable speeds in wavefront solutions of a scalar reaction-diffusion equation

Wavefront solutions of scalar reaction-diffusion equations havebeen intensively studied for many years. There are two basiccases, typified by quadratic and cubic kinetics. An intermediatecase is considered in this paper, namely, u_l = u_xx + u²(1 –u). It is shown that there is a unique travelling-wave solution,with a speed 1/2, for which the decay to zero ahead of the waveis exponential with x. Moreover, numerical evidence is presentedwhich suggests that initial conditions with such exponentialdecay evolve to this travelling-wave solution, independentlyof the half-life of the initial decay. It is then shown thatfor all speeds greater than 1/2 there is also a travelling-wavesolution, but that these faster waves decay to zero algebraically,in proportion to 1/x. The numerical evidence suggests that initialconditions with such a decay rate evolve to one of these fasterwaves; the particular speed depends in a simple way on the detailsof the initial condition. Finally, initial conditions decayingalgebraically for a power law other than 1/x are considered.It is shown numerically that such initial conditions evolveeither to an algebraically decaying travelling wave, or in somecases to a wavefront whose shape and speed vary as a functionof time. This variation is monotonic and can be quite pronounced,and the speed is a function of u as well as of time. Using asimple linearization argument, an approximate formula is derivedfor the wave speed which compares extremely well with the numericalresults. Finally, the extension of the results to the more generalcase of u_l = u_xx + u^m(1 – u), with m > 1, is discussed.     

Advection-mediated competition in general environments

We consider a reaction–diffusion–advection system of two competing species with one of the species dispersing by random diffusion as well as a biased movement upward along resource gradient, while the other species by random diffusion only. It has been shown that, under some non-degeneracy conditions on the environment function, the two species always coexist when the advection is strong. In this paper, we show that for general smooth environment function, in contrast to what is known, there can be competitive exclusion when the advection is strong, and, we give a sharp criterion for coexistence that includes all previously considered cases. Moreover, when the domain is one-dimensional, we derive in the strong advection limit a system of two equations defined on different domains. Uniqueness of steady states of this non-standard system is obtained when one of the diffusion rates is large.     

Analytical solution of a spatially variable coefficient advection–diffusion equation in up to three dimensions

Analytical solutions are provided for the two- and three-dimensional advection–diffusion equation with spatially variable velocity and diffusion coefficients. We assume that the velocity component is proportional to the distance and that the diffusion coefficient is proportional to the square of the corresponding velocity component. There is a simple transformation which reduces the spatially variable equation to a constant coefficient problem for which there are available a large number of known analytical solutions for general initial and boundary conditions. These solutions are also solutions to the spatially variable advection–diffusion equation. The special form of the spatial coefficients has practical relevance and for divergent free flow represent corner or straining flow. Unlike many other analytical solutions, we use the transformation to obtain solutions of the spatially variable coefficient advection–diffusion equation in two and three dimensions. The analytical solutions, which are simple to evaluate, can be used to validate numerical models for solving the advection–diffusion equation with spatially variable coefficients. For numerical schemes which cannot handle flow stagnation points, we provide analytical solution to the spatially variable coefficient advection–diffusion equation for two-dimensional corner flow which contains an impermeable flow boundary. The impermeable flow boundary coincides with a streamline along which the fluid velocity is finite but the concentration vanishes. This example is useful for validating numerical schemes designed to predict transport around a curved boundary.     

Equilibrium attractivity of Runge-Kutta methods

For dissipative differential equations y' = f (y) it is knownthat contractivity of the exact solution is reproduced by algebraicallystable Runge–Kutta methods. In this paper we investigatewhether a different property of the exact solution also holdsfor Runge–Kutta solutions. This property, called equilibriumattractivity, means that the norm of the righthand side f neverincreases. It is a property dual to algebraic stability sinceneither is sufficient for the other, in general. We derive sufficientalgebraic conditions for Runge–Kutta methods and proveequilibrium attractivity of the high-order algebraically stableRadau-IIA and Lobatto-IIIC methods and the Lobatto-IIIA collocationmethods (which are not algebraically stable). No smoothnessassumptions on f and no stepsize restrictions are required but,except for some simple cases, f has to satisfy certain additionalproperties which are generalizations of the simple one-sidedLipschitz condition using more than two argument points. Thesemultipoint conditions are discussed in detail.     

The Transient Motion of a Partially Immersed Rolling Strip in Water

An infinitely long thin strip is partly immersed in deep fluidand is fixed along its submerged end. It is given a small disturbancefrom its equilibrium position firstly by an applied rotationalforce and then by an initial angular displacement, the fluidcoming to rest before the strip, subject to a linear restoringforce, is released. Under the assumptions that viscosity andsurface tension may be neglected and the equations of motionlinearized, the ensuing two-dimensional displacement is describedby functions involving a Fourier integral. These are computedand the transient motion found. It is shown that after an initialstage the behaviour of the body is approximated by a dampedharmonic oscillatory motion which cannot be represented by asecond order differential equation with constant coefficients.This conflicts with the theory used by engineers and naval architectsin ship hydrodynamics. Ultimately the decay of motion is monotonic,decaying like t^–7 when an initial force is applied andt^–6 when there is an initial displacement. Comparisonis made with the corresponding behaviour of the undamped system.     

Hopf bifurcation in a delayed reaction–diffusion–advection population model

In this paper, we investigate a reaction–diffusion–advection model with time delay effect. The stability/instability of the spatially nonhomogeneous positive steady state and the associated Hopf bifurcation are investigated when the given parameter of the model is near the principle eigenvalue of an elliptic operator. Our results imply that time delay can make the spatially nonhomogeneous positive steady state unstable for a reaction–diffusion–advection model, and the model can exhibit oscillatory pattern through Hopf bifurcation. The effect of advection on Hopf bifurcation values is also considered, and our results suggest that Hopf bifurcation is more likely to occur when the advection rate increases.     

Pattern formation in a 2D simple chemical system with general orders of autocatalysis and decay

In this paper, we investigate pattern formation in a coupledsystem of reaction–diffusion equations in two spatialdimensions. These equations arise as a model of isothermal chemicalautocatalysis with termination in which the orders of autocatalysisand termination, m and n, respectively, are such that 1 <n < m. We build on the preliminary work by Leach & Wei(2003, Physica D, 180, 185–209) for this coupled systemin one spatial dimension, by presenting rigorous stability analysisand detailed numerical simulations for the coupled system intwo spatial dimensions. We demonstrate that spotty patternsare observed over a wide parameter range.     

On optimal convergence rates for higher-order Navier-Stokes approximations. I. Error estimates for the spatial discretization

^** Email: bause{at}am.uni-erlangen.de Due to the increasing use of higher-order methods in computationalfluid dynamics, the question of optimal approximability of theNavier–Stokes equations under realistic assumptions onthe data has become important. It is well known that the regularitycustomarily hypothesized in the error analysis for parabolicproblems cannot be assumed for the Navier–Stokes equations,as it depends on non-local compatibility conditions for thedata at time t = 0, which cannot be verified in practice. Takinginto account this loss of regularity at t = 0, improved convergenceof the order (min{h^(5/2)–,h³/t^(1/4)+}), for any >0, is shown locally in time for the spatial discretization ofthe velocity field by (non-)conforming finite elements of third-orderapproximability properties. The error estimate itself is provedby energy methods, but it is based on sharp a priori estimatesfor the Navier–Stokes solution in fractional-order spacesthat are derived by semigroup methods and complex interpolationtheory and reflect the optimal regularity of the solution ast 0.     

A WKB analysis of the buckling condition for a cylindrical shell of arbitrary thickness subjected to an external pressure

We consider the plane-strain buckling of a cylindrical shellof arbitrary thickness which is made of a Varga material andis subjected to an external hydrostatic pressure on its outersurface. The WKB method is used to solve the eigenvalue problemthat results from the linear bifurcation analysis. We show thatthe circular cross-section buckles into a non-circular shapeat a value of µ₁ which depends on A₁/A₂ and a mode number,where A₁ and A₂ are the undeformed inner and outer radii, andµ₁ is the ratio of the deformed inner radius to A₁. Inthe large mode number limit, we find that the dependence ofµ₁ on A₁/A₂ has a boundary layer structure: it is constantover almost the entire region of 0 < A₁/A₂ < 1 and decreasessharply from this constant value to unity as A₁/A₂ tends tounity. Our asymptotic results for A₁ – 1 = O(1) and A₁– 1 = O(1/n) are shown to agree with the numerical resultsobtained by using the compound matrix method.     

Stability and bifurcation in a two-species reaction–diffusion–advection competition model with time delay

In this paper, we are concerned with the dynamics of a class of two-species reaction–diffusion–advection competition models with time delay subject to the homogeneous Dirichlet boundary condition or no-flux boundary condition in a bounded domain. The existence of steady state solution is investigated by means of the Lyapunov–Schmidt reduction method. The stability and Hopf bifurcation at the spatially nonhomogeneous steady-state are obtained by analyzing the distribution of the associated eigenvalues. Finally, the effect of advection on Hopf bifurcation is explored, which shows that with the increase of convection rate, the Hopf bifurcation phenomenon is more likely to emerge.     

Centralizers in Domains of Gelfand-Kirillov Dimension 2

Given an affine domain of Gelfand–Kirillov dimension 2over an algebraically closed field, it is shown that the centralizerof any non-scalar element of this domain is a commutative domainof Gelfand–Kirillov dimension 1 whenever the domain isnot polynomial identity. It is shown that the maximal subfieldsof the quotient division ring of a finitely graded Goldie algebraof Gelfand–Kirillov dimension 2 over a field F all havetranscendence degree 1 over F. Finally, centralizers of elementsin a finitely graded Goldie domain of Gelfand–Kirillovdimension 2 over an algebraically closed field are considered.In this case, it is shown that the centralizer of a non-scalarelement is an affine commutative domain of Gelfand–Kirillovdimension 1. 2000 Mathematics Subject Classification 16P90.     

Weakly non-linear bifurcation analysis of pattern formation in strained alloy film growth

We consider a model for the growth of alloy films that includessurface diffusion, the effect of stresses due to misfit andstresses due to composition gradients when the alloy componentshave different sizes. Linear stability theory predicts a bifurcationfrom the planar homogeneous film to a non-planar compositionallymodulated film at a critical deposition rate. In this paper,we perform a weakly non-linear bifurcation analysis of hexagonaland band patterns using an asymptotic analysis of the systemclose to its critical state. A novel feature of the analysisis that the formulation of the adjoint problem involves thesolution of the composition-driven elasticity problem in thepresence of surface diffusion and requires multiple scales inthe growth direction. Our results characterize the transcriticalbifurcation to hexagons and the pitchfork bifurcation to bandsnear threshold. Finally, we apply our results to the growthof Si_1–XGe_X films on Si_0.5Ge_0.5 substrates and describehow the amplitude of surface undulations and the amplitude ofcompositional modulations corresponding to hexagons and bandsdepend on Ge composition.     

Wild Recurrent Critical Points

It is conjectured that a rational map whose coefficients arealgebraic over Q_p has no wandering components of the Fatou set.Benedetto has shown that any counterexample to this conjecturemust have a wild recurrent critical point. We provide the firstexamples of rational maps whose coefficients are algebraic over Q_p and that have a (wild) recurrent critical point. In fact,it is shown that there is such a rational map in every one-parameterfamily of rational maps that is defined over a finite extensionof Q_p and that has a Misiurewicz bifurcation.     

The streamline-diffusion method for nonconforming Qrot1 elements on rectangular tensor-product meshes

When the streamline–diffusion finite element method isapplied to convection–diffusion problems using nonconformingtrial spaces, it has previously been observed that stabilityand convergence problems may occur. It has consequently beenproposed that certain jump terms should be added to the bilinearform to obtain the same stability and convergence behaviouras in the conforming case. The analysis in this paper showsthat for the Q^rot₁ 1 element on rectangular shape-regular tensor-productmeshes, no jump terms are needed to stabilize the method. Inthis case moreover, for smooth solutions we derive in the streamline–diffusionnorm convergence of order h^3/2 (uniformly in the diffusion coefficientof the problem), where h is the mesh diameter. (This estimateis already known for the conforming case.) Our analysis alsoshows that similar stability and convergence results fail tohold true for analogous piecewise linear nonconforming elements.     

Latin Squares and the Hall-Paige Conjecture

The Hall–Paige conjecture deals with conditions underwhich a finite group G will possess a complete mapping, or equivalentlya Latin square based on the Cayley table of G will possess atransversal. Two necessary conditions are known to be: (i) thatthe Sylow 2-subgroups of G are trivial or non-cyclic, and (ii)that there is some ordering of the elements of G which yieldsa trivial product. These two conditions are known to be equivalent,but the first direct, elementary proof that (i) implies (ii)is given here. It is also shown that the Hall–Paige conjecture impliesthe existence of a duplex in every group table, thereby provinga special case of Rodney's conjecture that every Latin squarecontains a duplex. A duplex is a ‘double transversal’,that is, a set of 2n entries in a Latin square of order n suchthat each row, column and symbol is represented exactly twice.2000 Mathematics Subject Classification 05B15, 20D60.     

Solution to a Class of Coupled Linear Partial Differential Equations

The solution to a coupled system of partial differential equationsinvolving a general linear time-independent operator L is presented.Examples of these equations include coupled diffusion equationsor coupled convection–dispersion equations. The solutionconsists of a convolution of the Green's function appropriatefor the operator L and a function independent of the operatorL. The method enables one to write software to calculate thesolution to a wide range of problems. The change of solutionupon changing the problem often only involves a substitutionof the Green's function. A specific example of physical significanceis given.     

Notes on the Numerical Solution of the Biharmonic Equation

It is shown that the solution of the differenced form of theDirichlet problem for the biharmonic equation by iteration isachieved in O(h^–1)log h² steps of iteration where h isa step-size parameter from the mesh.     

Estimates of the indeterminacy set for elliptic boundary value problems with uncertain data

Estimates are derived for the so-called indeterminacy set formed by solutions to an elliptic type boundary value problem with not fully determined coefficients. A two-sided estimate for the diameter of the indeterminacy set is obtained in the energy norm. It is shown that this estimate depends on parameters defining the variability range of the coefficients. The analysis is based on functional a posteriori estimates that provide guaranteed bounds of the difference between an approximate solution and any admissible function in the energy space. The estimates are obtained for the diffusion equation. However, the proposed tools can be used for other classes of partial differential equations if functional a posteriori estimates are established. Bibliography: 2 titles. __________ Translated from Problemy Matematicheskogo Analiza, No. 36, 2007, pp. 77–80.