Algebraic isomorphisms and strongly double triangle subspace lattices

Let D={{0},K,L,M,X} be a strongly double triangle subspace lattice on a non-zero complex reflexive Banach space X, which means that at least one of three sums K + L, L + M and M + K is closed. It is proved that a non-zero element S of AlgD is single in the sense that for any A,B∈AlgD, either AS = 0 or SB = 0 whenever ASB = 0, if and only if S is of rank two. We also show that every algebraic isomorphism between two strongly double triangle subspace lattice algebras is quasi-spatial.     

The representations of the Drazin inverse of differences of two matrices

In this paper, we give explicit representations of (P ± Q)^D, (P ± PQ)^D and (PQ)^# of two matrices P and Q, as a function of P, Q, P^D and Q^D, under the conditions P³Q = QP and Q³P = PQ.     

Falkner-Skan boundary layer flow of a power-law fluid past a stretching wedge

The steady two-dimensional laminar boundary layer flow of a power-law fluid past a permeable stretching wedge beneath a variable free stream is studied in this paper. Using appropriate similarity variables, the governing equations are reduced to a single third order highly nonlinear ordinary differential equation in the dimensionless stream function, which is solved numerically using the Runge-Kutta scheme coupled with a conventional shooting procedure. The flow is governed by the wedge velocity parameter λ, the transpiration parameter f₀, the fluid power-law index n, and the computed wall shear stress is f″(0). It is found that dual solutions exist for each value of f₀, m and n considered in λ − f″(0) parameter space. A stability analysis for this self-similar flow reveals that for each value of f₀, m and n, lower solution branches are unstable while upper solution branches are stable. Very good agreements are found between the results of the present paper and that of Weidman et al. [28] for n = 1 (Newtonian fluid) and m = 0 (Blasius problem [31]).     

Nonexistence of Isochronous Centers in Planar Polynomial Hamiltonian Systems of Degree Four

In this paper we study centers of planar polynomial Hamiltonian systems and we are interested in the isochronous ones. We prove that every center of a polynomial Hamiltonian system of degree four (that is, with its homogeneous part of degree four not identically zero) is nonisochronous. The proof uses the geometric properties of the period annulus and it requires the study of the Hamiltonian systems associated to a Hamiltonian function of the form H(x, y)=A(x)+B(x) y+C(x) y²+D(x) y³.     

Approximation algorithms for constructing some required structures in digraphs

We consider a new problem of constructing some required structures in digraphs, where all arcs installed in such required structures are supposed to be cut from some pieces of a specific material of length L. Formally, we consider the model: a digraph D = (V, A; w), a structure S and a specific material of length L, where w: A → R⁺, we are asked to construct a subdigraph D′ from D, having the structure S, such that each arc in D′ is constructed by a part of a piece or/and some whole pieces of such a specific material, the objective is to minimize the number of pieces of such a specific material to construct all arcs in D′.     

The common Re-nnd and Re-pd solutions to the matrix equations AX = C and XB = D

In this article, we consider common Re-nnd and Re-pd solutions of the matrix equations AX = C and XB = D with respect to X, where A, B, C and D are given matrices. We give necessary and sufficient conditions for the existence of common Re-nnd and Re-pd solutions to the pair of the matrix equations and derive a representation of the common Re-nnd and Re-pd solutions to these two equations when they exist. The presented examples show the advantage of the proposed approach.     

Thermal dispersion driven by the spontaneous imbibition process

In this work, we have theoretically analyzed the thermal dispersion process under the influence of the spontaneous imbibition of a liquid trapped in a capillary element, considering the presence of a uniform temperature gradient. The capillary element is represented by a porous medium which is initially found at temperature T₀ and pressure P₀. Suddenly, the lower part of the porous medium touches a liquid reservoir at temperature T_l and pressure P₀. This contact between both phases, in turn causes spontaneously the imbibition process. Using a one-dimensional formulation of the average conservation laws, we derive the corresponding nondimensional momentum and energy equations. The numerical solutions permit us to evaluate the position and velocity of the imbibition front as well as the temperature profiles and Nusselt numbers. The above results are shown by taking into account the influence of three dimensionless parameters: the ratio of the characteristic thermal time to the characteristic imbibition time, β, the ratio of the hydrostatic head of the imbibed liquid to the characteristic pressure difference for the imbibition front, α, and the ratio of the dispersive thermal diffusivity to the effective thermal diffusivity of the medium, Ω. The predictions show that temperature profiles and the heat transfer process are strongly dependent on thermal dispersion effects, indicating a clear deviation in comparison with the case of Ω = 0 that represents the absence of the thermal dispersion.     

New results on the energy of integral circulant graphs

Circulant graphs are an important class of interconnection networks in parallel and distributed computing. Integral circulant graphs play an important role in modeling quantum spin networks supporting the perfect state transfer as well. The integral circulant graph ICG_n(D) has the vertex set Z_n = {0, 1, 2, … , n − 1} and vertices a and b are adjacent if gcd(a − b, n) ∈ D, where D ⊆ {d : d∣n, 1 ? d < n}. These graphs are highly symmetric, have integral spectra and some remarkable properties connecting chemical graph theory and number theory. The energy of a graph was first defined by Gutman, as the sum of the absolute values of the eigenvalues of the adjacency matrix. Recently, there was a vast research for the pairs and families of non-cospectral graphs having equal energies. Following Bapat and Pati [R.B. Bapat, S. Pati, Energy of a graph is never an odd integer, Bull. Kerala Math. Assoc. 1 (2004) 129-132], we characterize the energy of integral circulant graph modulo 4. Furthermore, we establish some general closed form expressions for the energy of integral circulant graphs and generalize some results from Ili? [A. Ili?, The energy of unitary Cayley graphs, Linear Algebra Appl. 431 (2009), 1881-1889]. We close the paper by proposing some open problems and characterizing extremal graphs with minimal energy among integral circulant graphs with n vertices, provided n is even.     

Numerical computation of end plate effect on Taylor vortices between rotating conical cylinders

End plate effect on Taylor vortices between rotating conical cylinders is studied by numerical method in this paper. We suppose that the inner cone rotates together with the end plate at the top and the outer one as well as the end plate at the bottom remains at rest. It is found that the instability sets in at a critical Reynolds number about Re = 80. Increase Re to about Re = 200 the first single clockwise vortex is formed near the top of the flow system. Further increase Re to about Re = 440 another clockwise vortex is formed under the first one. At about Re = 448 the third vortex is formed which rotates in counterclockwise direction between the first two vortices. With increasing of Re the process continues. Finally, a configuration is obtained with an odd number of vortices in the annulus at about Re = 700, which confirms the experimental observation. Moreover, the local extreme values of pressure and velocity are achieved at the adjacent lines between neighboring vortices or at the medium lines of vortices. The effect of gap size on vortices is also considered. It is shown that the number of vortices increases with decreasing of the gap size and the end plates play an important role in the parity of the number of the vortices.     

Representations for the Drazin inverses of 2 × 2 block matrices

Let M denote a 2 × 2 block complex matrix , where A and D are square matrices, not necessarily with the same orders. In this paper explicit representations for the Drazin inverse of M are presented under the condition that BDⁱC = 0 for i = 0, 1, … , n − 1, where n is the order of D.     

Classification of small (0, 1) matrices

Denote by A_n the set of square (0, 1) matrices of order n. The set A_n, n ? 8, is partitioned into row/column permutation equivalence classes enabling derivation of various facts by simple counting. For example, the number of regular (0, 1) matrices of order 8 is 10160459763342013440. Let D_n, S_n denote the set of absolute determinant values and Smith normal forms of matrices from A_n. Denote by a_n the smallest integer not in D_n. The sets D₉ and S₉ are obtained; especially, a₉ = 103. The lower bounds for a_n, 10 ? n ? 19 (exceeding the known lower bound a_n ? 2f_n − 1, where f_n is nth Fibonacci number) are obtained. Row/permutation equivalence classes of A_n correspond to bipartite graphs with n black and n white vertices, and so the other applications of the classification are possible.     

Viscous flow and temperature distribution in rotating annulus with suction and injection

Part A.—Steady viscous incompressible flow in a rotating coaxial cylindrical annulus with suction and injection is studied. The unsteady flow is also considered. Part B.—An exact solution for temperature distribution at different constant wall temperatures is obtained. It is assumed that the rate of injection at the inner wall equals the rate of suction at the outer wall.     

Compositions of Polynomials with Coefficients in a Given Field

Let F ⊂ K be fields of characteristic 0, and let K[x] denote the ring of polynomials with coefficients in K. Let p(x) = ∑_k = 0ⁿa_kx^k ∈ K[x], a_n ≠ 0. For p ∈ K[x]\F[x], define D_F(p), the F deficit of p, to equal n − max{0 ≤ k ≤ n : a_k∉F}. For p ∈ F[x], define D_F(p) = n. Let p(x) = ∑_k = 0ⁿa_kx^k and let q(x) = ∑_j = 0^mb_jx^j, with a_n ≠ 0, b_m ≠ 0, a_n, b_m ∈ F, b_j∉F for some j ≥ 1. Suppose that p ∈ K[x], q ∈ K[x]\F[x], p, not constant. Our main result is that p ° q ∉ F[x] and D_F(p ° q) = D_F(q). With only the assumption that a_nb_m ∈ F, we prove the inequality D_F(p ° q) ≥ D_F(q). This inequality also holds if F and K are only rings. Similar results are proven for fields of finite characteristic with the additional assumption that the characteristic of the field does not divide the degree of p. Finally we extend our results to polynomials in two variables and compositions of the form p(q(x, y)), where p is a polynomial in one variable.     

Eigenvector-free solutions to AX = B with PX = XP and X = sX constraints

The matrix equation AX = B with PX = XP and X^H = sX constraints is considered, where P is a given Hermitian involutory matrix and s = ±1. By an eigenvalue decomposition of P, we equivalently transform the constrained problem to two well-known constrained problems and represent the solutions in terms of the eigenvectors of P. Using Moore-Penrose generalized inverses of the products generated by matrices A, B and P, the involved eigenvectors can be released and eigenvector-free formulas of the general solutions are presented. Similar strategy is applied to the equations AX = B, XC = D with the same constraints.     

Fundamental solutions of singular SPDEs

This paper deals with some models of mathematical physics, where random fluctuations are modeled by white noise or other singular Gaussian generalized processes. White noise, as the distributional derivative od Brownian motion, which is the most important case of a Lévy process, is defined in the framework of Hida distribution spaces. The Fourier transformation in the framework of singular generalized stochastic processes is introduced and its applications to solving stochastic differential equations involving Wick products and singularities such as the Dirac delta distribution are presented. Explicit solutions are obtained in form of a chaos expansion in the Kondratiev white noise space, while the coefficients of the expansion are tempered distributions. Stochastic differential equations of the form P(ω, D) ◊ u(x, ω) = A(x, ω) are considered, where A is a singular generalized stochastic process and P(ω, D) is a partial differential operator with random coefficients. We introduce the Wick-convolution operator which enables us to express the solution as u = sA ◊ I^◊(−1), where s denotes the fundamental solution and I is the unit random variable. In particular, the stochastic Helmholtz equation is solved, which in physical interpretation describes waves propagating with a random speed from randomly appearing point sources.     

Numerical study of melting in an enclosure with discrete protruding heat sources

In order to explore the capability of a solid–liquid phase change material (PCM) for cooling electronic or heat storage applications, melting of a PCM in a vertical rectangular enclosure was studied. Three protruding generating heat sources are attached on one of the vertical walls of the enclosure, and generating heat at a constant and uniform volumetric rate. The horizontal walls are adiabatic. The power generated in heat sources is dissipated in PCM (n-eicosane with the melting temperature, T_m = 36 °C) that filled the rectangular enclosure. The advantage of using PCM is that it is able to absorb high amount of heat generated by heat sources due to its relatively high energy density. To investigate the thermal behaviour and thermal performance of the proposed system, a mathematical model based on the mass, momentum and energy conservation equations was developed. The governing equations are next discretised using a control volume approach in a staggered mesh and a pressure correction equation method is employed for the pressure–velocity coupling. The PCM energy equation is solved using the enthalpy method. The solid regions (wall and heat sources) are treated as fluid regions with infinite viscosity and the thermal coupling between solid and fluid regions is taken into account using the harmonic mean of the thermal conductivity method. The dimensionless independent parameters that govern the thermal behaviour of the system were next identified. After validating the proposed mathematical model against experimental data, a numerical investigation was next conducted in order to examine the thermal behaviour of the system by analyzing the flow structure and the heat transfer during the melting process, for a given values of governing parameters.     

Mean matrices and infinite divisibility

We consider matrices M with entries m_ij = m(λ_i, λ_j) where λ₁, … ,λ_n are positive numbers and m is a binary mean dominated by the geometric mean, and matrices W with entries w_ij = 1/m (λ_i, λ_j) where m is a binary mean that dominates the geometric mean. We show that these matrices are infinitely divisible for several much-studied classes of means.     

A non-Newtonian fluid flow model for blood flow through a catheterized artery—Steady flow

In this paper the effects catheterization and non-Newtonian nature of blood in small arteries of diameter less than 100 μm, on velocity, flow resistance and wall shear stress are analyzed mathematically by modeling blood as a Herschel–Bulkley fluid with parameters n and θ and the artery and catheter by coaxial rigid circular cylinders. The influence of the catheter radius and the yield stress of the fluid on the yield plane locations, velocity distributions, flow rate, wall shear stress and frictional resistance are investigated assuming the flow to be steady. It is shown that the velocity decreases as the yield stress increases for given values of other parameters. The frictional resistance as well as the wall shear stress increases with increasing yield stress, whereas the frictional resistance increases and the wall shear stress decreases with increasing catheter radius ratio k (catheter radius to vessel radius). For the range of catheter radius ratio 0.3–0.6, in smaller arteries where blood is modeled by Herschel–Bulkley fluid with yield stress θ = 0.1, the resistance increases by a factor 3.98–21.12 for n = 0.95 and by a factor 4.35–25.09 for n = 1.05. When θ = 0.3, these factors are 7.47–124.6 when n = 0.95 and 8.97–247.76 when n = 1.05.     

Dirac-bracket structure in multidimensional mode conversion

The intersection of two (2n − 1)-dimensional dispersion manifolds D_a and D_b in the 2n-dimensional ray phase space P yields a (2n − 2)-dimensional conversion manifold M≡D_a∩D_b that naturally possesses a Dirac-bracket structure that is inherited from the canonical Poisson bracket on ray phase space. The canonical symplectic two-form Ω ≡ Ω_∥ + Ω_⊥, defined on the 2n-dimensional tangent plane T_z0P≡T_z0M⊕(T_z0M)_⊥, can thus be decomposed into the Dirac two-form Ω_∥ on the (2n − 2)-dimensional tangent plane T_z0M at a conversion point z₀∈M, and the symplectic two-form Ω_⊥ on its orthogonal 2-dimensional complement (T_z0M)_⊥. These two symplectic two-forms are introduced in our analysis of multidimensional mode conversion, where their respective geometrical roles are defined. We note that since the Dirac-bracket structure Ω_∥ vanishes identically when n = 1, it represents a new structure in multidimensional (n > 1) mode conversion theory.     

Automatic Controls for Nonlinear Dynamical Systems with Lipschitzian Trajectories

We would like to investigate on the solution to the automatic control problem given by the differential equation y′(t) = f(t, y(t), w(t)) for a given initial function x in the initial domain D(x, ω, Y) for almost all t in the interval I, with controls given by w(t) = g(t, y(t), T(y)(t)), where T is a nonanticipating and Lipschitzian operator. The result will be generalized for a dynamical system y′(t) = f(t, y(t), T(y), u(t)).