Falkner-Skan equation for flow past a moving wedge with suction or injection

The characteristics of steady two-dimensional laminar boundary layer flow of a viscous and incompressible fluid past a moving wedge with suction or injection are theoretically investigated. The transformed boundary layer equations are solved numerically using an implicit finite-difference scheme known as the Keller-box method. The effects of Falkner-Skan power-law parameter (m), suction/injection parameter (f0) and the ratio of free stream velocity to boundary velocity parameter (λ) are discussed in detail. The numerical results for velocity distribution and skin friction coefficient are given for several values of these parameters. Comparisons with the existing results obtained by other researchers under certain conditions are made. The critical values off ₀,m and λ are obtained numerically and their significance on the skin friction and velocity profiles is discussed. The numerical evidence would seem to indicate the onset of reverse flow as it has been found by Riley and Weidman in 1989 for the Falkner-Skan equation for flow past an impermeable stretching boundary.     

Falkner-Skan equation for flow past a stretching surface with suction or blowing: Analytical solutions

The simultaneous effects of suction and injection on tangential movement of a nonlinear power-law stretching surface governed by laminar boundary layer flow of a viscous and incompressible fluid beneath a non-uniform free with stream pressure gradient is considered. The self-similar flow is governed by Falkner-Skan equation, with transpiration parameter γ, wall slip velocity λ and stretching sheet (or pressure gradient) parameter β. The exact solution for β = −1 and three closed form asymptotic solutions for β large, large suction γ, and λ → 1 have also been presented. Dual solutions are found for β = −1 for each value of the transpiration parameter, including the non-permeable surface, for each prescribed value of the wall slip velocity λ. The large β asymptotic solution also dual with respect to wall slip velocity λ, but do not depend on suction and blowing. The critical values of γ, β and λ are obtained and their significance on the skin friction and velocity profiles is discussed. An approximate solution by integral method for a trial velocity profile is presented and results are compared with the exact solutions.     

On a certain boundary value problem arising in shrinking sheet flows

This paper presents an analysis of the boundary value problem resulting from the magnetohydrodynamic (MHD) viscous flow influenced by a shrinking sheet with suction for the cases of two-dimensional (m = 1) and axisymmetric (m = 2) shrinking. The influences of the parameter m as well as the effects of suction parameter s and Hartmann number M² on similar entrainment velocity f(∞) and flow characteristics are studied. To this purpose, the resulting nonlinear ordinary differential equation is solved numerically using the 4th order Runge-Kutta method in combination with a shooting procedure. The obtained results elucidate reliability and efficiency of the technique from which interesting features between the skin friction coefficient f″(0) and the entrainment velocity f(∞) as function of the mass transfer parameter s can also be obtained.     

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.     

Construction of nonnegative symmetric matrices with given spectrum

Let σ = (λ₁, … , λ_n) be the spectrum of a nonnegative symmetric matrix A with the Perron eigenvalue λ₁, a diagonal entry c and let τ = (μ₁, … , μ_m) be the spectrum of a nonnegative symmetric matrix B with the Perron eigenvalue μ₁. We show how to construct a nonnegative symmetric matrix C with the spectrum

(λ₁+max{0,μ₁-c},λ₂,…,λ_n,μ₂,…,μ_m).     

Normal Forms for Nonautonomous Differential Equations

We extend Henry Poincarés normal form theory for autonomous differential equations x=f(x) to nonautonomous differential equations x=f(t, x). Poincarés nonresonance condition λ_j−∑ⁿ_i=1 ?_iλ_i≠0 for eigenvalues is generalized to the new nonresonance condition λ_j∩∑ⁿ_i=1 ?_iλ_i=∅ for spectral intervals.     

On singular integrals depending on three parameters

In this paper we will prove the pointwise convergence of L(f; x, y, λ) to f(x₀, y₀), as (x, y, λ) tends to (x₀, y₀, λ₀) in the space L_2π, by the three parameter family of singular operators. In contrast to previous works, the kernel function is radial.     

Bernstein operators for extended Chebyshev systems

Let U_n ⊂ Cⁿ[a, b] be an extended Chebyshev space of dimension n + 1. Suppose that f₀ ∈ U_n is strictly positive and f₁ ∈ U_n has the property that f₁/f₀ is strictly increasing. We search for conditions ensuring the existence of points t₀, …, t_n ∈ [a, b] and positive coefficients α₀, …, α_n such that for all f ∈ C[a, b], the operator B_n:C[a, b] → U_n defined by satisfies B_nf₀ = f₀ and B_nf₁ = f₁. Here it is assumed that p_n,k, k = 0, …, n, is a Bernstein basis, defined by the property that each p_n,k has a zero of order k at a and a zero of order n − k at b.     

Solving second order initial value problems by a hybrid multistep method without predictors

A three-step seventh order hybrid linear multistep method (HLMM) with three non-step points is proposed for the direct solution of the special second order initial value problems (IVPs) of the form y″ = f(x, y) with an extension to y″ = f(x, y, y′). The main method and additional methods are obtained from the same continuous scheme derived via interpolation and collocation procedures.The stability properties of the methods are discussed by expressing them as a one-step method in higher dimension. The methods are then applied in block form as simultaneous numerical integrators over non-overlapping intervals. Numerical results obtained using the proposed block form reveal that it is highly competitive with existing methods in the literature.     

Equivalence of inverse Sturm-Liouville problems with boundary conditions rationally dependent on the eigenparameter

Three inverse problems for a Sturm-Liouville boundary value problem −y″+qy=λy, y(0)cosα=y′(0)sinα and y′(1)=f(λ)y(1) are considered for rational f. It is shown that the Weyl m-function uniquely determines α, f, and q, and is in turn uniquely determined by either two spectra from different values of α or by the Prüfer angle. For this it is necessary to produce direct results, of independent interest, on asymptotics and oscillation.     

Analyzing the longitudinal effect of hypersonic flow past a conical cone via the perturbation method

To analyze the hypersonic flow past a conical cone, the variations of gasdynamic properties subjected to the longitudinal curvature effect by using the perturbation method. An outer perturbation expansion has been carried out by recent researchers, but a problem occurred, the outer expansion solutions are not uniformly valid in the shock layer, however, the outcome near the conical body surface called vortical layer remains deflective. This study intends to discover uniformly valid analytical solutions in the shock layer by applying the inner perturbation expansions matching with the out expansions to analyze the characteristics in the whole region including shock layer and vortical layer. Starting from the zero-order approximate solutions for hypersonic conical flow is then applied as the basic solutions for the outer perturbation expansions of a flow field. The governing equations and boundary conditions are also expanded via outer perturbations. Using an approximate analytical scheme in the derivation process, first-order perturbation equations can be simplified and the approximate closed-form solutions are obtained; furthermore, the various flow field quantities, including the normal force coefficient on the cone surface, have been calculated. According to the variations of gasdynamic properties, the longitudinal curvature effect for the hypersonic flow past a conical cone can be determined. Thicknesses of shock layer and vortical layer can be predicted as well. The physical phenomena inside both layers can be investigated carefully, the conditions for an elliptic cone with longitudinal curvature, m = 1 and n = 2 and other conditions of parameters; the perturbation parameter, ε_m2 = 0.1, semi-vertex angle of the unperturbed cone, δ = 10°, and hypersonic similarity parameter, K_δ = M_∞δ = 1.0, the thickness of vortical layer, η_VL, can be calculated at the position angle of conical cone body, ? = 30° was demonstrated here. Results show how very thin the vortical layer is approximately only 10% of the shock layer close to the body, the pressure in the whole shock layer is verified to be uniformly valid which agrees with previous studies. Large gradient changes in entropy and density are found when the flow approaches the cone surface, the most important is, this method provides a benchmark solution to the hypersonic flow past a conical cone and to assist the grids and numerics for numerical computation should be fashioned to accommodate the whole flow field region including the vortical layer of rapid adjustment, and let the analysis become more effective and low cost.     

On single and double Soules matrices

Until now the concept of a Soules basis matrix of sign patternN consisted of an orthogonal matrix R∈R^n,n, generated in a certain way from a positive n-vector, which has the property that for any diagonal matrix Λ = diag(λ₁, … , λ_n), with λ₁ ? ? ? λ_n ? 0, the symmetric matrix A = RΛR^T has nonnegative entries only. In the present paper we introduce the notion of a pair of double Soules basis matrices of sign patternN which is a pair of matrices (P, Q), each in R^n,n, which are not necessarily orthogonal and which are generated in a certain way from two positive vectors, but such that PQ^T = I and such that for any of the aforementioned diagonal matrices Λ, the matrix A = PΛQ^T (also) has nonnegative entries only. We investigate the interesting properties which such matrices A have.As a preamble to the above investigation we show that the iterates, , generated in the course of the QR-algorithm when it is applied to A = RΛR^T, where R is a Soules basis matrix of sign pattern N, are again symmetric matrices generated by the Soules basis matrices R_k of sign pattern N which are themselves modified as the algorithm progresses.Our work here extends earlier works by Soules and Elsner et al.     

Pointwise Estimate for Linear Combinations of Bernstein-Kantorovich Operators

For linear combinations of Bernstein-Kantorovich operators K_n, r(f, x), we give an equivalent theorem with ω^2r_?^λ(f, t). The theorem unites the corresponding results of classical and Ditzian-Totik moduli of smoothness.     

The domination numbers of cylindrical grid graphs

Let γ(P_m □ C_n) denote the domination number of the cylindrical grid graph formed by the Cartesian product of the graphs P_m, the path of length m, m ? 2 and the graph C_n, the cycle of length n, n ? 3. In this paper, methods to find the domination numbers of graphs of the form P_m □ C_n with n ? 3 and m = 2, 3 and 4 are proposed. Moreover, bounds on domination numbers of the graphs P₅ □ C_n, n ? 3 are found. The methods that are used to prove that results readily lead to algorithms for finding minimum dominating sets of the above mentioned graphs.     

Laplacian energy of a graph

Let G be a graph with n vertices and m edges. Let λ₁, λ₂, … , λ_n be the eigenvalues of the adjacency matrix of G, and let μ₁, μ₂, … , μ_n be the eigenvalues of the Laplacian matrix of G. An earlier much studied quantity is the energy of the graph G. We now define and investigate the Laplacian energy as . There is a great deal of analogy between the properties of E(G) and LE(G), but also some significant differences.     

Similarity solutions for mixed convection boundary layer flow over a permeable horizontal flat plate

The effects of suction and injection on steady laminar mixed convection boundary layer flow over a permeable horizontal flat plate in a viscous and incompressible fluid is investigated in this paper. The similarity solutions of the governing boundary layer equations are obtained for some values of the suction and injection parameter f₀, the constant exponent n of the wall temperature as well as the mixed convection parameter λ. The resulting system of nonlinear ordinary differential equations is solved numerically for both assisting and opposing flow regimes using a finite-difference scheme known as the Keller-box method. Numerical results for the reduced skin friction coefficient, the reduced local Nusselt number, and the velocity and temperature profiles are obtained for various values of the parameters considered. Dual solutions are found to exist for the opposing flow.     

Inequalities between ∥f(A + B)∥ and ∥f(A) + f(B)∥

The conjecture posed by Aujla and Silva [J.S. Aujla, F.C. Silva, Weak majorization inequalities and convex functions, Linear Algebra Appl. 369 (2003) 217-233] is proved. It is shown that for any m-tuple of positive-semidefinite n × n complex matrices A_j and for any non-negative convex function f on [0, ∞) with f(0) = 0 the inequality ?f(A₁) + f(A₂) + ? + f(A_m)? ? ? f(A₁ + A₂ + ? + A_m)? holds for any unitarily invariant norm ? · ?. It is also proved that ?f(A₁) + f(A₂) + ? + f(A_m)? ? f(?A₁ + A₂ + ? + A_m?), where f is a non-negative concave function on [0, ∞) and ? · ? is normalized.     

Meromorphic Functions That Share Fixed-Points

In this paper, we prove the following result: Let f(z) and g(z) be two nonconstant meromorphic(entire) functions, n ≥ 11(n ≥ 6) a positive integer. If fⁿ(z)f′(z) and gⁿ(z)g′(z) have the same fixed-points, then either f(z) = c₁e^cz2, g(z) = c₂e^− cz2, where c₁, c₂, and c are three constants satisfying 4(c₁c₂)^n + 1c² = −1, or f(z) ≡ tg(z) for a constant t such that t^n + 1 = 1.     

Bounding the gap between extremal Laplacian eigenvalues of graphs

Let G be a graph whose Laplacian eigenvalues are 0 = λ₁ ? λ₂ ? ? ? λ_n. We investigate the gap (expressed either as a difference or as a ratio) between the extremal non-trivial Laplacian eigenvalues of a connected graph (that is λ_n and λ₂). This gap is closely related to the average density of cuts in a graph. We focus here on the problem of bounding the gap from below.