A second‐order linearized difference scheme for a strongly coupled reaction‐diffusion system

This article deals with the numerical solution to some models described by the system of strongly coupled reaction–diffusion equations with the Neumann boundary value conditions. A linearized three‐level scheme is derived by the method of reduction of order. The uniquely solvability and second‐order convergence in L₂‐norm are proved by the energy method. A numerical example is presented to demonstrate the accuracy and efficiency of the proposed method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Finite difference discretization of the Benjamin‐Bona‐Mahony‐Burgers equation

Numerical solutions of the Benjamin‐Bona‐Mahony‐Burgers equation in one space dimension are considered using Crank‐Nicolson‐type finite difference method. Existence of solutions is shown by using the Brower's fixed point theorem. The stability and uniqueness of the corresponding methods are proved by the means of the discrete energy method. The convergence in L^∞‐norm of the difference solution is obtained. A conservative difference scheme is presented for the Benjamin‐Bona‐Mahony equation. Some numerical experiments have been conducted in order to validate the theoretical results.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

The recognition of symmetric latin squares

A latin square S is isotopic to another latin square S′ if S′ can be obtained from S by permuting the row indices, the column indices and the symbols in S. Because the three permutations used above may all be different, a latin square which is isotopic to a symmetric latin square need not be symmetric. We call the problem of determining whether a latin square is isotopic to a symmetric latin square the symmetry recognition problem. It is the purpose of this article to give a solution to this problem. For this purpose we will introduce a cocycle corresponding to a latin square which transforms very simply under isotopy, and we show this cocycle contains all the information needed to determine whether a latin square is isotopic to a symmetric latin square. Our results relate to 1‐factorizations of the complete graph on n + 1 vertices, K_{n + 1}. There is a well known construction which can be used to make an n × n latin square from a 1‐factorization on n + 1 vertices. The symmetric idempotent latin squares are exactly the latin squares that result from this construction. The idempotent recognition problem is simple for symmetric latin squares, so our results enable us to recognize exactly which latin squares arise from 1‐factorizations of K_{n + 1}. As an example we show that the patterned starter 1‐factorization for the group G gives rise to a latin square which is in the main class of the Cayley latin square for G if and only if G is abelian. Hence, every non‐abelian group gives rise to two latin squares in different main classes. © 2007 Wiley Periodicals, Inc. J Combin Designs 16: 291–300, 2008     

Constructing orthogonal pandiagonal Latin squares and panmagic squares from modular n‐queens solutions

In this article, we show how to construct pairs of orthogonal pandiagonal Latin squares and panmagic squares from certain types of modular n‐queens solutions. We prove that when these modular n‐queens solutions are symmetric, the panmagic squares thus constructed will be associative, where for an n × n associative magic square A = (a_ij), for all i and j it holds that a_ij + an?i?1,n?j?1 = c for a fixed c. We further show how to construct orthogonal Latin squares whose modular difference diagonals are Latin from any modular n‐queens solution. As well, we analyze constructing orthogonal pandiagonal Latin squares from particular classes of non‐linear modular n‐queens solutions. These pandiagonal Latin squares are not row cyclic, giving a partial solution to a problem of Hedayat. © 2007 Wiley Periodicals, Inc. J Combin Designs 15: 221–234, 2007     

Technical note: A fourth‐order finite difference scheme for a system of a 2D nonlinear elliptic partial differential equations

In this article we present a fourth‐order finite difference scheme, for a system of two‐dimensional, second‐order, nonlinear elliptic partial differential equations with mixed spatial derivative terms, using 13‐point stencils with a uniform mesh size h on a square region R subject to Dirichlet boundary conditions. The scheme of order h⁴ is derived using the local solution of the system on a single stencil. The resulting system of algebraic equations can be solved by iterative methods. The difference scheme can be easily modified to obtain formulae for grid points near the boundary. Computational results are given to demonstrate the performance of the scheme on some problems including Navier‐Stokes equations. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 43–53, 2001     

An almost L‐stable BDF‐type method for the numerical solution of stiff ODEs arising from the method of lines

A new BDF‐type scheme is proposed for the numerical integration of the system of ordinary differential equations that arises in the Method of Lines solution of time‐dependent partial differential equations. This system is usually stiff, so it is desirable for the numerical method to solve it to have good properties concerning stability. The method proposed in this article is almost L‐stable and of algebraic order three. Numerical experiments illustrate the performance of the new method on different stiff systems of ODEs after discretizing in the space variable some PDE problems. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Complex dynamics of a discrete fractional‐order Leslie‐Gower predator‐prey model

A proposed discretized form of fractional‐order prey‐predator model is investigated. A sufficient condition for the solution of the discrete system to exist and to be unique is determined. Jury stability test is applied for studying stability of equilibrium points of the discretized system. Then, the effects of varying fractional order and other parameters of the systems on its dynamics are examined. The system undergoes Neimark‐Sacker and flip bifurcation under certain conditions. We observe that the model exhibits chaotic dynamics following stable states as the memory parameter α decreases and step size h increases. Theoretical results illustrate the rich dynamics and complexity of the model. Numerical simulation validates theoretical results and demonstrates the presence of rich dynamical behaviors include S‐asymptotically bounded periodic orbits, quasi‐periodicity, and chaos. The system exhibits a wide range of dynamical behaviors for fractional‐order α key parameter.     

A two‐grid characteristic finite volume element method for semilinear advection‐dominated diffusion equations

A two‐grid finite volume element method, combined with the modified method of characteristics, is presented and analyzed for semilinear time‐dependent advection‐dominated diffusion equations in two space dimensions. The solution of a nonlinear system on the fine‐grid space (with grid size h) is reduced to the solution of two small (one linear and one nonlinear) systems on the coarse‐grid space (with grid size H) and a linear system on the fine‐grid space. An optimal error estimate in H¹ ‐norm is obtained for the two‐grid method. It shows that the two‐grid method achieves asymptotically optimal approximation, as long as the mesh sizes satisfy h = O(H²). Numerical example is presented to validate the usefulness and efficiency of the method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Orthogonal spline collocation methods for the stream function‐vorticity formulation of the Navier–Stokes equations

An orthogonal spline collocation (OSC) spatial discretization is proposed for the solution of the fully coupled stream function‐vorticity formulation of the Navier–Stokes equations in two dimensions. For the time‐stepping, a three‐level leapfrog scheme is employed. This method is algebraically linear, and, at each time step, gives rise to a system of linear equations of the form arising in the OSC approximation of the biharmonic Dirichlet problem and can be solved by a fast direct method. Error estimates in the H^l–norm in space, l = 1,2, are derived for the semi‐discrete method and the fully‐discrete leapfrog scheme which is also shown to be second order accurate in time. Numerical results are presented which confirm the theoretical analysis and exhibit superconvergence phenomena, which provide superconvergent approximations to the components of the velocity. © John Wiley & Sons, Inc. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

H∞ robust control based on event‐triggered sampling for hybrid systems with singular Markovian jump

The problem of H_∞ robust control based on event‐triggered sampling for a class of singular hybrid systems with Markovian jump is considered in this paper. The primary object of this paper here is to design the event‐triggered sampling controller for a class of uncertain singular Markovian systems, and two fundamental issues on mean square exponential admissibility and H_∞ robust performance are fully addressed. By making use of a suitable Lyapunov functional, in combination with both infinitesimal operator and linear matrices inequalities(LMIs), the sufficient criteria are derived to guarantee the controlled singular hybrid system with Markovian jump is robustly exponentially mean‐square admissible and has a prescribed H_∞ performance γ. Finally, a typical RLC circuit system is given to show the effectiveness of the proposed control method.     

The conjugate gradient method solves the neutron transport equation h‐optimally

We report a serendipitous finding that the conjugate gradient method is able to obtain the solution to the linear system derived by the discontinuous Galerkin method for the discrete‐ordinate approximation of the mono‐energetic time‐independent neutron transport equation to a specified accuracy in a number of iterations which is bounded independently of the order of the system. A theoretical explanation for this phenomenon is given. The results of an experiment to test the theory are presented. Copyright © 2007 John Wiley & Sons, Ltd.     

Lorenz‐gauged vector potential formulations for the time‐harmonic eddy‐current problem with L∞‐regularity of material properties

In this paper, we consider some Lorenz‐gauged vector potential formulations of the eddy‐current problem for the time‐harmonic Maxwell equations with material properties having only L^∞‐regularity. We prove that there exists a unique solution of these problems, and we show the convergence of a suitable finite element approximation scheme. Moreover, we show that some previously proposed Lorenz‐gauged formulations are indeed formulations in terms of the modified magnetic vector potential, for which the electric scalar potential is vanishing. Copyright © 2007 John Wiley & Sons, Ltd.     

Regularity and derivative bounds for a convection-diffusion problem with a Neumann outflow condition

A convection-diffusion problem is considered on the unit square. The convective direction is parallel to two of the square's sides. A Neumann condition is imposed on the outflow boundary, with Dirichlet conditions on the other three sides. The precise relationship between the regularity of the solution and the global smoothness and corner compatibility of the data is elucidated. Pointwise bounds on derivatives of the solution are obtained; their dependence on the data regularity and compatibility and on the small diffusion parameter is made explicit. The analysis uses Fourier transforms and Mikhlin multipliers to sharpen regularity results previously published for certain subproblems in a decomposition of the solution.     

Mixed two‐grid finite difference methods for solving one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations

The aim of this paper is to propose mixed two‐grid finite difference methods to obtain the numerical solution of the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. The finite difference equations at all interior grid points form a large‐sparse linear system, which needs to be solved efficiently. The solution cost of this sparse linear system usually dominates the total cost of solving the discretized partial differential equation. The proposed method is based on applying a family of finite difference methods for discretizing the spatial and time derivatives. The obtained system has been solved by two‐grid method, where the two‐grid method is used for solving the large‐sparse linear systems. Also, in the proposed method, the spectral radius with local Fourier analysis is calculated for different values of h and Δt. The numerical examples show the efficiency of this algorithm for solving the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. Copyright © 2016 John Wiley & Sons, Ltd.     

Non‐homogeneous Navier–Stokes systems with order‐parameter‐dependent stresses

We consider the Navier–Stokes system with variable density and variable viscosity coupled to a transport equation for an order‐parameter c. Moreover, an extra stress depending on c and ?c, which describes surface tension like effects, is included in the Navier–Stokes system. Such a system arises, e.g. for certain models of granular flows and as a diffuse interface model for a two‐phase flow of viscous incompressible fluids. The so‐called density‐dependent Navier–Stokes system is also a special case of our system. We prove short‐time existence of strong solution in L^q‐Sobolev spaces with q>d. We consider the case of a bounded domain and an asymptotically flat layer with a combination of a Dirichlet boundary condition and a free surface boundary condition. The result is based on a maximal regularity result for the linearized system. Copyright © 2010 John Wiley & Sons, Ltd.     

On the spectrum of critical sets in latin squares of order 2n

Suppose that L is a latin square of order m and P ? L is a partial latin square. If L is the only latin square of order m which contains P, and no proper subset of P has this property, then P is a critical set of L. The critical set spectrum problem is to determine, for a given m, the set of integers t for which there exists a latin square of order m with a critical set of size t. We outline a partial solution to the critical set spectrum problem for latin squares of order 2ⁿ. The back circulant latin square of even order m has a well‐known critical set of size m²/4, and this is the smallest known critical set for a latin square of order m. The abelian 2‐group of order 2ⁿ has a critical set of size 4ⁿ‐3ⁿ, and this is the largest known critical set for a latin square of order 2ⁿ. We construct a set of latin squares with associated critical sets which are intermediate between the back circulant latin square of order 2ⁿ and the abelian 2‐group of order 2ⁿ. © 2007 Wiley Periodicals, Inc. J Combin Designs 16: 25–43, 2008     

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.     

A remark on the decay of solutions to the 3‐D Navier–Stokes equations

In this paper we derive a decay rate of the L²‐norm of the solution to the 3‐D Navier–Stokes equations. Although the result which is proved by Fourier splitting method is well known, our method is new, concise and direct. Moreover, it turns out that the new method established here has a wide application on other equations. Copyright © 2007 John Wiley & Sons, Ltd.     

The Vlasov‐Poisson‐Boltzmann system near Maxwellians

The dynamics of dilute electrons can be modeled by the Vlasov‐Poisson‐Boltz‐mann system, where electrons interact with themselves through collisions and with their self‐consistent electric field. It is shown that any smooth, periodic initial perturbation of a given global Maxwellian that preserves the same mass, momentum, and total energy (including both kinetic and electric energy), leads to a unique global‐in‐time classical solution. The construction of global solutions is based on an energy method with a new estimate of dissipation from the collision: ∫₀^t〈Lf(s), f(s)〉ds is positive definite for solution f(t,x,v) with small amplitude to the Vlasov‐Poisson‐Boltzmann system (1.4). © 2002 Wiley Periodicals, Inc.     

Global weak solutions for a periodic two‐component μ‐Camassa–Holm system

In this paper, we consider the global existence of weak solutions for a two‐component μ‐Camassa–Holm system in the periodic setting. Global existence for strong solutions to the system with smooth approximate initial value is derived. Then, we show that the limit of approximate solutions is a global‐in‐time weak solution of the two‐component μ‐Camassa–Holm system. Copyright © 2013 John Wiley & Sons, Ltd.