A lowest order divergence-free finite element on rectangular grids

It is shown that the conforming Q _2,1;1,2-Q′₁ mixed element is stable, and provides optimal order of approximation for the Stokes equations on rectangular grids. Here, Q _2,1;1,2 = Q _2,1 × Q _1,2, and Q _2,1 denotes the space of continuous piecewise-polynomials of degree 2 or less in the x direction but of degree 1 in the y direction. Q′₁ is the space of discontinuous bilinear polynomials, with spurious modes filtered. To be precise, Q′₁ is the divergence of the discrete velocity space Q _2,1;1,2. Therefore, the resulting finite element solution for the velocity is divergence-free pointwise, when solving the Stokes equations. This element is the lowest order one in a family of divergence-free element, similar to the families of the Bernardi-Raugel element and the Raviart-Thomas element.     

Group divisible designs with three unequal groups and larger first index

We show the necessary conditions are sufficient for the existence of group divisible designs (or PBIBDs) with block size k=3 with three groups of size (n,2,1) for any n≥2 and any two indices with λ₁>λ₂.     

Auslander–Bridger Modules

Classically, the Auslander–Bridger transpose finds its best applications in the well-known setting of finitely presented modules over a semiperfect ring. We introduce a class of modules over an arbitrary ring R, which we call Auslander–Bridger modules, with the property that the Auslander–Bridger transpose induces a well-behaved bijection between isomorphism classes of Auslander–Bridger right R-modules and isomorphism classes of Auslander–Bridger left R-modules. Thus we generalize what happens for finitely presented modules over a semiperfect ring. Auslander–Bridger modules are characterized by two invariants (epi-isomorphism class and lower-isomorphism class), which are interchanged by the transpose. Via a suitable duality, we find that kernels of morphisms between injective modules of finite Goldie dimension are also characterized by two invariants (mono-isomorphism class and upper-isomorphism class).     

Generating equations approach for quadratic matrix equations

We show how Van Loan's method for annulling the (2,1) block of skew‐Hamiltonian matrices by symplectic‐orthogonal similarity transformation generalizes to general matrices and provides a numerical algorithm for solving the general quadratic matrix equation: For skew‐Hamiltonian matrices we find their canonical form under a similarity transformation and find the class of all symplectic‐orthogonal similarity transformations for annulling the (2,1) block and simultaneously bringing the (1,1) block to Hessenberg form. We present a structure‐preserving algorithm for the solution of continuous‐time algebraic Riccati equation. Unlike other methods in the literature, the final transformed Hamiltonian matrix is not in Hamiltonian–Schur form. Three applications are presented: (a) for a special system of partial differential equations of second order for a single unknown function, we obtain the matrix of partial derivatives of second order of the unknown function by only algebraic operations and differentiation of functions; (b) for a similar transformation of a complex matrix into a symmetric (and three‐diagonal) one by applying only finite algebraic transformations; and (c) for finite‐step reduction of the eigenvalues–eigenvectors problem of a Hermitian matrix to the eigenvalues– eigenvectors problem of a real symmetric matrix of the same dimension. Copyright © 1999 John Wiley & Sons, Ltd.     

Distributions of eigenvalues and inertias of some block Hermitian matrices consisting of orthogonal projectors

A complex square matrix A is called an orthogonal projector if A ²?=?A?=?A*, where A* is the conjugate transpose of A. In this article, we first give some formulas for calculating the distributions of real eigenvalues of a linear combination of two orthogonal projectors. Then, we establish various expansion formulas for calculating the inertias, ranks and signatures of some 2?×?2 and 3?×?3, as well as k?×?k block Hermitian matrices consisting of two orthogonal projectors. Many applications of the formulas are presented in characterizing interval distributions of numbers of eigenvalues, and nonsingularity of these block Hermitian matrices. In addition, necessary and sufficient conditions are given for various equalities and inequalities of these block Hermitian matrices to hold.     

Lanczos,Householder transformations,and implicit deflation for fast and reliable dominant singular subspace computation

Many applications, such as subspace‐based models in information retrieval and signal processing, require the computation of singular subspaces associated with the k dominant, or largest, singular values of an m×n data matrix A, where k?min(m,n). Frequently, A is sparse or structured, which usually means matrix–vector multiplications involving A and its transpose can be done with much less than ??(mn) flops, and A and its transpose can be stored with much less than ??(mn) storage locations. Many Lanczos‐based algorithms have been proposed through the years because the underlying Lanczos method only accesses A and its transpose through matrix–vector multiplications. We implement a new algorithm, called KSVD, in the Matlab environment for computing approximations to the singular subspaces associated with the k dominant singular values of a real or complex matrix A. KSVD is based upon the Lanczos tridiagonalization method, the WY representation for storing products of Householder transformations, implicit deflation, and the QR factorization. Our Matlab simulations suggest it is a fast and reliable strategy for handling troublesome singular‐value spectra. Copyright © 2001 John Wiley & Sons, Ltd.     

The existence of doubly near resolvable (v,3,2)-BIBDs

The existence of doubly near resolvable (v,2,1)-BIBDs was established by Mullin and Wallis in 1975. In this article, we determine the spectrum of a second class of doubly near resolvable balanced incomplete block designs. We prove the existence of DNR(v,3,2)-BIBDs for v ≡ 1 (mod 3), v ≥ 10 and v ? {34,70,85,88,115,124,133,142}. The main construction is a frame construction, and similar constructions can be used to prove the existence of doubly resolvable (v,3,2)-BIBDs and a class of Kirkman squares with block size 3, KS₃(v,2,4). © 1994 John Wiley & Sons, Inc.     

Regularity and solvability of linear differential operators in Gevrey spaces

We are interested in the following question: when regularity properties of a linear differential operator imply solvability of its transpose in the sense of Gevrey ultradistributions? This question is studied for a class of abstract operators that contains the usual differential operators with real‐analytic coefficients. We obtain a new proof of a global result on compact manifolds (global Gevrey hypoellipticity implying global solvability of the transpose), as well as some results in the non‐compact case by means of the so‐called property of non‐confinement of singularities. We provide applications to Hörmander operators, to operators of constant strength and to locally integrable systems of vector fields. We also analyze a conjecture stated in a recent paper of Malaspina and Nicola, which asserts that, in differential complexes naturally arising from locally integrable structures, local solvability in the sense of ultradistributions implies local solvability in the sense of distributions. We establish the validity of the conjecture when the cotangent structure bundle is spanned by the differential of a single first integral.     

A quadratic functional for a third order linear difference equation

We will define a certain quadratic functional and use it to prove various results for the third order difference equation l₃y(t)=Δ³y(t-1)+p(t)Δy(t)+q(t)y(t)=0. In particular we will define kth order generlized zeroes for solutions of this equations and define (2, 1)- and (1,2)-disconjugacy of l₃y=0 on [a,b+3]. Then we will use our quadratic functional to prove sufficient conditions for (2,1)- and (1,2)-disconjugacy. We will also discuss what we call type I and II solutions of l₃y=0 and give properties of these solutions. These later results give asymptotic behavior of solutions at infinity.     

Least‐squares solutions and least‐rank solutions of the matrix equation AXA* = B and their relations

A Hermitian matrix X is called a least‐squares solution of the inconsistent matrix equation AXA^* = B, where B is Hermitian. A^* denotes the conjugate transpose of A if it minimizes the F‐norm of B ? AXA^*; it is called a least‐rank solution of AXA^* = B if it minimizes the rank of B ? AXA^*. In this paper, we study these two types of solutions by using generalized inverses of matrices and some matrix decompositions. In particular, we derive necessary and sufficient conditions for the two types of solutions to coincide. Copyright © 2012 John Wiley & Sons, Ltd.     

Labeling planar graphs with a condition at distance two

An L(2,1)-labeling of a graph is a mapping c:V(G)→{0,…,K} such that the labels assigned to neighboring vertices differ by at least 2 and the labels of vertices at distance two are different. The smallest K for which an L(2,1)-labeling of a graph G exists is denoted by λ_2,1(G). Griggs and Yeh [J.R. Griggs, R.K. Yeh, Labeling graphs with a condition at distance 2, SIAM J. Discrete Math. 5 (1992) 586–595] conjectured that λ_2,1(G)≤Δ² for every graph G with maximum degree Δ≥2. We prove the conjecture for planar graphs with maximum degree Δ≠3. All our results also generalize to the list-coloring setting.     

The unstable integral homology of the mapping class groups of a surface with boundary

We construct a graph complex calculating the integral homology of the bordered mapping class groups. We compute the homology of the bordered mapping class groups of the surfaces S _1,1, S _1,2 and S _2,1. Using the circle action on this graph complex, we build a double complex and a spectral sequence converging to the homology of the unbordered mapping class groups. We compute the homology of the punctured mapping class groups associated to the surfaces S _1,1 and S _2,1. Finally, we use Miller’s operad to get the first Kudo–Araki and Browder operations on our graph complex. We also consider an unstable version of the higher Kudo–Araki–Dyer–Lashoff operations.     

Density results relative to the Dirichlet energy of mappings into a manifold

Let ?? be a smooth, compact, oriented Riemannian manifold without boundary. Weak limits of graphs of smooth maps u_k:Bⁿ → ?? with an equibounded Dirichlet integral give rise to elements of the space cart^2,1 (Bⁿ × ??). Assume that ?? is 1‐connected and that its 2‐homology group has no torsion. In any dimension n we prove that every element T in cart^2,1 (Bⁿ × ??) with no singular vertical part can be approximated weakly in the sense of currents by a sequence of graphs of smooth maps u_k:Bⁿ → ?? with Dirichlet energies converging to the energy of T. © 2006 Wiley Periodicals, Inc.     

Analytical approximate solutions of Riesz fractional diffusion equation and Riesz fractional advection–dispersion equation involving nonlocal space fractional derivatives

In this paper, we consider the analytical solutions of fractional partial differential equations (PDEs) with Riesz space fractional derivatives on a finite domain. Here we considered two types of fractional PDEs with Riesz space fractional derivatives such as Riesz fractional diffusion equation (RFDE) and Riesz fractional advection–dispersion equation (RFADE). The RFDE is obtained from the standard diffusion equation by replacing the second‐order space derivative with the Riesz fractional derivative of order α∈(1,2]. The RFADE is obtained from the standard advection–dispersion equation by replacing the first‐order and second‐order space derivatives with the Riesz fractional derivatives of order β∈(0,1] and of order α∈(1,2] respectively. Here the analytic solutions of both the RFDE and RFADE are derived by using modified homotopy analysis method with Fourier transform. Then, we analyze the results by numerical simulations, which demonstrate the simplicity and effectiveness of the present method. Here the space fractional derivatives are defined as Riesz fractional derivatives. Copyright © 2015 John Wiley & Sons, Ltd.     

Reduction of the adjoint representation of sl(2,1): generators of primitive ideals

We present a reduction of the adjoint representation of the Lie superalge-bra,sl(2,1) and a study of the quotient algebra B(^c,k)= u/u(C?c)+u(D?kc), where c,k are two complex numbers. Under some additional conditions, we prove that every irreducible infinite dimensional representation of B(^c,k) is faithful, and that B(^C,K) is a primitive algebra. We give explicitly a set of generators of primitive degenerate ideal of infinite codimension. Essentially we prove that any minimal primitive ideal of u(sl(2,1)) is generated, as a 2-sided ideal, by its intersection with the algebra of gg-iuvariants.     

The solution of the equation AX + BX ⋆ = 0

We give a complete solution of the matrix equation AX?+?BX ^??=?0, where A, B?∈?? ^m×n are two given matrices, X?∈?? ^n×n is an unknown matrix, and ? denotes the transpose or the conjugate transpose. We provide a closed formula for the dimension of the solution space of the equation in terms of the Kronecker canonical form of the matrix pencil A?+?λB, and we also provide an expression for the solution X in terms of this canonical form, together with two invertible matrices leading A?+?λB to the canonical form by strict equivalence.     

Procrustes problems and inverse eigenproblems for multilevel block α‐circulants

Let n = (n₁,n₂,…,n_k) and α = (α₁,α₂,…,α_k) be integer k‐tuples with α_i∈{1,2,…,n_i?1} and for all i = 1,2,…,k. Multilevel block α ‐circulants are (k + 1)‐level block matrices, where the first k levels have the block α_i‐circulant structure with orders and the entries in the (k + 1)‐st level are unstructured rectangular matrices with the same size . When k = 1, Trench discussed on his paper "Inverse problems for unilevel block α‐circulants" the Procrustes problems and inverse problems of unilevel block α‐circulants and their approximations. But the results are not perfect for the case gcd( α , n ) > 1 (i.e., gcd(α₁,n₁) > 1). In this paper, we also discuss Procrustes problems for multilevel block α ‐circulants. Our results can further make up for the deficiency when k = 1. When , inverse eigenproblems for this kind of matrices are also solved. By using the related results, we can design an artificial Hopfield neural network system that possesses the prescribed equilibria, where the Jacobian matrix of this system has the constrained multilevel α ‐circulative structure. Finally, some examples are employed to illustrate the effectiveness of the proposed results. Copyright © 2016 John Wiley & Sons, Ltd.     

Asymptotic profile of solutions for strongly damped Klein‐Gordon equations

We consider the Cauchy problem in R ⁿ for strongly damped Klein‐Gordon equations. We derive asymptotic profiles of solutions with weighted L^1,1( R ⁿ) initial data by a simple method introduced by the second author. Furthermore, from the obtained asymptotic profile, we get the optimal decay order of the L²‐norm of solutions. The obtained results show that the wave effect will be relatively weak because of the mass term, especially in the low‐dimensional case (n = 1,2) as compared with the strongly damped wave equations without mass term (m = 0), so the most interesting topic in this paper is the n = 1,2 cases to compare the difference.     

Partial Regularity for Singular Solutions to the Monge‐Ampère Equation

We prove that solutions to the Monge‐Ampère inequality in ?ⁿ are strictly convex away from a singular set of Hausdorff (n‐1)‐dimensional measure zero. Furthermore, we show this is optimal by constructing solutions to det D²u = 1 with singular set of Hausdorff dimension as close as we like to n‐1. As a consequence we obtain W^2,1 regularity for the Monge‐Ampère equation with bounded right‐hand side and unique continuation for the Monge‐Ampère equation with sufficiently regular right‐hand side. © 2015 Wiley Periodicals, Inc.     

The (3,4) Mode Interaction in the GeoFlow Framework

The problem of convection in a self‐gravitating spherical shell of fluid is commonly encountered in sciences like astrophysics and geophysics (earth's liquid core). The GEOFLOW‐experiment is a project of the European Space Agency in order to perform the spherical Rayleigh‐Bénard convection problem on the International Space Station in a micro‐gravity environment: the central force field is simulated by a dielectrophoretic one. Beyond a critical Rayleigh number Ra_c, generically an unique spherical ℓ mode becomes unstable and only stationary or travelling waves solutions are expected near the onset. But, for a critical aspect ratio η_c two consecutive modes (ℓ, ℓ + 1) are unstable. The (1,2) and (2,3) interactions have showed a rich bifurcation diagram, in particular, we have found heteroclinic cycles predicted by the theoretical study. Because of the experiment requirements, only the (3,4) one is possible. So, this paper purposes to analyse this bifurcation in non‐rotating case in the GEOFLOWframework using the theory of bifurcation with the spherical symmetry.