Viscous compressible barotropic symmetric flows with free boundary under general mass force Part I: Uniform‐in‐time bounds and stabilization

We consider symmetric flows of a viscous compressible barotropic fluid with a free boundary, under a general mass force depending both on the Eulerian and Lagrangian co‐ordinates, with arbitrarily large initial data. For a general non‐monotone state function p, we prove uniform‐in‐time energy bound and the uniform bounds for the density ρ, together with the stabilization as t → ∞ of the kinetic and potential energies. We also obtain H¹‐stabilization of the velocity v to zero provided that the second viscosity is zero. For either increasing or non‐decreasing p, we study the L^λ‐stabilization of ρ and the stabilization of the free boundary together with the corresponding ω‐limit set in the general case of non‐unique stationary solution possibly with zones of vacuum. In the case of increasing p and stationary densities ρ_S separated from zero, we establish the uniform‐in‐time H¹‐bounds and the uniform stabilization for ρ and v. All these results are stated and mainly proved in the Eulerian co‐ordinates. They are supplemented with the corresponding stabilization results in the Lagrangian co‐ordinates in the case of ρ_S separated from zero. Copyright © 2005 John Wiley & Sons, Ltd.     

Modeling λ‐invariants by p‐adic random matrices

How often is the p‐adic λ‐invariant of an imaginary quadratic field equal to m? We model this by the statistics of random p‐adic matrices, and test these predictions numerically. © 2011 Wiley Periodicals, Inc.     

Methods for the numerical solution of the Benjamin‐Bona‐Mahony‐Burgers equation

L^∞‐error estimates for finite element for Galerkin solutions for the Benjamin‐Bona‐Mahony‐Burgers (BBMB) equation are considered. A priori bound and the semidiscrete Galerkin scheme are studied using appropriate projections. For fully discrete Galerkin schemes, we consider the backward Euler method and analyze the corresponding error estimates. For a second order accuracy in time, we propose a three‐level backward method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

On cubic non‐Cayley vertex‐transitive graphs

In 1983, the second author [D. Maru?i?, Ars Combinatoria 16B (1983), 297–302] asked for which positive integers n there exists a non‐Cayley vertex‐transitive graph on n vertices. (The term non‐Cayley numbers has later been given to such integers.) Motivated by this problem, Feng [Discrete Math 248 (2002), 265–269] asked to determine the smallest valency ?(n) among valencies of non‐Cayley vertex‐transitive graphs of order n. As cycles are clearly Cayley graphs, ?(n)?3 for any non‐Cayley number n. In this paper a goal is set to determine those non‐Cayley numbers n for which ?(n) = 3, and among the latter to determine those for which the generalized Petersen graphs are the only non‐Cayley vertex‐transitive graphs of order n. It is known that for a prime p every vertex‐transitive graph of order p, p² or p³ is a Cayley graph, and that, with the exception of the Coxeter graph, every cubic non‐Cayley vertex‐transitive graph of order 2p, 4p or 2p² is a generalized Petersen graph. In this paper the next natural step is taken by proving that every cubic non‐Cayley vertex‐transitive graph of order 4p², p>7 a prime, is a generalized Petersen graph. In addition, cubic non‐Cayley vertex‐transitive graphs of order 2p^k, where p>7 is a prime and k?p, are characterized. © 2011 Wiley Periodicals, Inc. J Graph Theory 69: 77–95, 2012     

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.     

Non‐Cayley‐isomorphic self‐complementary circulant graphs

Non‐CI self‐complementary circulant graphs of prime‐squared order are constructed and enumerated. It is shown that for prime p, there exists a self‐complementary circulant graph of order p² not Cayley isomorphic to its complement if and only if p ≡ 1 (mod 8). Such graphs are also enumerated. © 2000 John Wiley & Sons, Inc. J Graph Theory 34: 128–141, 2000     

Asymptotic and numerical solutions of three‐dimensional boundary‐layer flow past a moving wedge

We consider a laminar boundary‐layer flow of a viscous and incompressible fluid past a moving wedge in which the wedge is moving either in the direction of the mainstream flow or opposite to it. The mainstream flows outside the boundary layer are approximated by a power of the distance from the leading boundary layer. The variable pressure gradient is imposed on the boundary layer so that the system admits similarity solutions. The model is described using 3‐dimensional boundary‐layer equations that contains 2 physical parameters: pressure gradient (β) and shear‐to‐strain‐rate ratio parameter (α). Two methods are used: a linear asymptotic analysis in the neighborhood of the edge of the boundary layer and the Keller‐box numerical method for the full nonlinear system. The results show that the flow field is divided into near‐field region (mainly dominated by viscous forces) and far‐field region (mainstream flows); the velocity profiles form through an interaction between 2 regions. Also, all simulations show that the subsequent dynamics involving overshoot and undershoot of the solutions for varying parameter characterizing 3‐dimensional flows. The pressure gradient (favorable) has a tendency of decreasing the boundary‐layer thickness in which the velocity profiles are benign. The wall shear stresses increase unboundedly for increasing α when the wedge is moving in the x‐direction, while the case is different when it is moving in the y‐direction. Further, both analysis show that 3‐dimensional boundary‐layer solutions exist in the range −1<α<∞. These are some interesting results linked to an important class of boundary‐layer flows.     

Entropy flux-splittings for hyperbolic conservation laws part I: General framework

A general framework is proposed for the derivation and analysis of flux-splittings and the corresponding flux-splitting schemes for systems of conservation laws endowed with a strictly convex entropy. The approach leads to several new properties of the existing flux-splittings and to a method for the construction of entropy flux-splittings for general situations. A large family of genuine entropy flux-splittings is derived for several significant examples: the scalar conservation laws, the p-system, and the Euler system of isentropic gas dynamics. In particular, for the isentropic Euler system, we obtain a family of splittings that satisfy the entropy inequality associated with the mechanical energy. For this system, it is proved that there exists a unique genuine entropy flux-splitting that satisfies all of the entropy inequalities, which is also the unique diagonalizable splitting. This splitting can be also derived by the so-called kinetic formulation. Simple and useful difference schemes are derived from the flux-splittings for hyperbolic systems. Such entropy flux-splitting schemes are shown to satisfy a discrete cell entropy inequality. For the diagonalizable splitting schemes, an a priori L^∞ estimate is provided by applying the principle of bounded invariant regions. The convergence of entropy flux-splitting schemes is proved for the 2 × 2 systems of conservation laws and the isentropic Euler system. ©1995 John Wiley & Sons, Inc.     

Tetravalent vertex‐transitive graphs of order 4p

A graph is vertex‐transitive if its automorphism group acts transitively on vertices of the graph. A vertex‐transitive graph is a Cayley graph if its automorphism group contains a subgroup acting regularly on its vertices. In this article, the tetravalent vertex‐transitive non‐Cayley graphs of order 4p are classified for each prime p. As a result, there are one sporadic and five infinite families of such graphs, of which the sporadic one has order 20, and one infinite family exists for every prime p>3, two families exist if and only if p≡1 (mod 8) and the other two families exist if and only if p≡1 (mod 4). For each family there is a unique graph for a given order. © 2011 Wiley Periodicals, Inc.     

G‐invariant persistent homology

Classical persistent homology is a powerful mathematical tool for shape comparison. Unfortunately, it is not tailored to study the action of transformation groups that are different from the group Homeo(X) of all self‐homeomorphisms of a topological space X. This fact restricts its use in applications. In order to obtain better lower bounds for the natural pseudo‐distance d_G associated with a group G ? Homeo(X), we need to adapt persistent homology and consider G‐invariant persistent homology. Roughly speaking, the main idea consists in defining persistent homology by means of a set of chains that is invariant under the action of G. In this paper, we formalize this idea and prove the stability of the persistent Betti number functions in G‐invariant persistent homology with respect to the natural pseudo‐distance d_G. We also show how G‐invariant persistent homology could be used in applications concerning shape comparison, when the invariance group is a proper subgroup of the group of all self‐homeomorphisms of a topological space. In this paper, we will assume that the space X is triangulable, in order to guarantee that the persistent Betti number functions are finite without using any tameness assumption. Copyright © 2014 John Wiley & Sons, Ltd.     

3‐Factor‐Criticality of Vertex‐Transitive Graphs

A graph of order n is p ‐factor‐critical, where p is an integer of the same parity as n, if the removal of any set of p vertices results in a graph with a perfect matching. 1‐factor‐critical graphs and 2‐factor‐critical graphs are factor‐critical graphs and bicritical graphs, respectively. It is well known that every connected vertex‐transitive graph of odd order is factor‐critical and every connected nonbipartite vertex‐transitive graph of even order is bicritical. In this article, we show that a simple connected vertex‐transitive graph of odd order at least five is 3‐factor‐critical if and only if it is not a cycle.     

Blow‐up phenomena for a class of fourth‐order nonlinear wave equations with a viscous damping term

This paper deals with the blow‐up phenomena for a class of fourth‐order nonlinear wave equations with a viscous damping term with Ω = (0,1) and α > 0. Here, f(s) is a given nonlinear smooth function. For 0 < α < p – 1, we prove that the blow‐up occurs in finite time for arbitrary positive initial energy and suitable initial data. This result extends the recent results obtained by Xu et al. (Applicable Analysis)(2013) and Chen and Lu (J. Math. Anal. Appl.)(2009). Copyright © 2015 John Wiley & Sons, Ltd.     

Vertex‐Transitive Cubic Graphs of Square‐Free Order

A classification of connected vertex‐transitive cubic graphs of square‐free order is provided. It is shown that such graphs are well‐characterized metacirculants (including dihedrants, generalized Petersen graphs, Möbius bands), or Tutte's 8‐cage, or graphs arisen from simple groups PSL(2, p).     

Single‐cell compact finite‐difference discretization of order two and four for multidimensional triharmonic problems

In this article, we discuss finite‐difference methods of order two and four for the solution of two‐and three‐dimensional triharmonic equations, where the values of u,(?²u/?n²) and (?⁴u/?n⁴) are prescribed on the boundary. For 2D case, we use 9‐ and for 3D case, we use 19‐ uniform grid points and a single computational cell. We introduce new ideas to handle the boundary conditions and do not require to discretize the boundary conditions at the boundary. The Laplacian and the biharmonic of the solution are obtained as byproduct of the methods. The resulting matrix system is solved by using the appropriate block iterative methods. Computational results are provided to demonstrate the fourth‐order accuracy of the proposed methods. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Stability and convergence of two‐grid Crank‐Nicolson extrapolation scheme for the time‐dependent natural convection equations

In this paper, we consider the Crank‐Nicolson extrapolation scheme for the 2D/3D unsteady natural convection problem. Our numerical scheme includes the implicit Crank‐Nicolson scheme for linear terms and the recursive linear method for nonlinear terms. Standard Galerkin finite element method is used to approximate the spatial discretization. Stability and optimal error estimates are provided for the numerical solutions. Furthermore, a fully discrete two‐grid Crank‐Nicolson extrapolation scheme is developed, the corresponding stability and convergence results are derived for the approximate solutions. Comparison from aspects of the theoretical results and computational efficiency, the two‐grid Crank‐Nicolson extrapolation scheme has the same order as the one grid method for velocity and temperature in H¹‐norm and for pressure in L²‐norm. However, the two‐grid scheme involves much less work than one grid method. Finally, some numerical examples are provided to verify the established theoretical results and illustrate the performances of the developed numerical schemes.     

Aubry‐Mather theory and periodic solutions of the forced Burgers equation

Consider a Hamiltonian system with Hamiltonian of the form H(x, t, p) where H is convex in p and periodic in x, and t and x ∈ ℝ¹. It is well‐known that its smooth invariant curves correspond to smooth Z²‐periodic solutions of the PDE u_t + H(x, t, u)_x = 0. In this paper, we establish a connection between the Aubry‐Mather theory of invariant sets of the Hamiltonian system and Z²‐periodic weak solutions of this PDE by realizing the Aubry‐Mather sets as closed subsets of the graphs of these weak solutions. We show that the complement of the Aubry‐Mather set on the graph can be viewed as a subset of the generalized unstable manifold of the Aubry‐Mather set, defined in (2.24). The graph itself is a backward‐invariant set of the Hamiltonian system. The basic idea is to embed the globally minimizing orbits used in the Aubry‐Mather theory into the characteristic fields of the above PDE. This is done by making use of one‐ and two‐sided minimizers, a notion introduced in [12] and inspired by the work of Morse on geodesics of type A [26]. The asymptotic slope of the minimizers, also known as the rotation number, is given by the derivative of the homogenized Hamiltonian, defined in [21]. As an application, we prove that the Z²‐periodic weak solution of the above PDE with given irrational asymptotic slope is unique. A similar connection also exists in multidimensional problems with the convex Hamiltonian, except that in higher dimensions, two‐sided minimizers with a specified asymptotic slope may not exist. © 1999 John Wiley & Sons, Inc.     

Stability of small‐amplitude shock profiles of symmetric hyperbolic‐parabolic systems

Combining pointwise Green's function bounds obtained in a companion paper [36] with earlier, spectral stability results obtained in [16], we establish nonlinear orbital stability of small‐amplitude Lax‐type viscous shock profiles for the class of dissipative symmetric hyperbolic‐parabolic systems identified by Kawashima [20], notably including compressible Navier‐Stokes equations and the equations of magnetohydrodynamics, obtaining sharp rates of decay in L^p with respect to small L¹ ∩ H³ perturbations, 2 ≤ p ≤ ∞. Our analysis extends and somewhat refines the approach introduced in [35] to treat stability of relaxation profiles. © 2004 Wiley Periodicals, Inc.     

The second‐order upwind finite difference fractional steps method for moving boundary value problem of oil‐water percolation

The research of the three‐dimensional (3D) compressible miscible (oil and water) displacement problem with moving boundary values is of great value to the history of oil‐gas transport and accumulation in basin evolution, as well as to the rational evaluation in prospecting and exploiting oil‐gas resources, and numerical simulation of seawater intrusion. The mathematical model can be described as a 3D‐coupled system of nonlinear partial differential equations with moving boundary values. For a generic case of 3D‐bounded region, a kind of second‐order upwind finite difference fractional steps schemes applicable to parallel arithmetic is put forward. Some techniques, such as the change of variables, calculus of variations, and the theory of a priori estimates, are adopted. Optimal order estimates in l² norm are derived for the errors in approximate solutions. The research is important both theoretically and practically for model analysis in the field, for model numerical method and for software development. Thus, the well‐known problem has been solved.Copyright © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1103–1129, 2014     

A convolution formula for time‐invariant linear communication systems

For every linear and time‐invariant time‐discrete (communication) system T : l^∞ → l^∞, formally, the following convolution formula can be derived: \input amssym.def $$(Tf)(n)=\sum_{k \in {\Bbb Z}} h(n‐k) f(k),\quad n \in {\Bbb Z}, f \in l^{\infty},$$ where h = Tδ is the delta impulse response. This paper is concerned with the question under which assumptions linear and time‐invariant time‐discrete systems T : l^∞ → l^∞ can be characterized by this formula. For this purpose we derive a convolution formula in a more general situation which also leads to a well‐known convolution formula in the time‐continuous case. Copyright © 1999 John Wiley & Sons, Ltd.