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.     

脉冲波作用下的高超音速非定常流动模拟及影响研究

采用高阶精度有限差分方法模拟了快声波脉冲扰动作用下的高超音速非定常流场,分析了脉冲波与高超音速流场的相互干扰,并应用Fourier频谱分析研究扰动波在边界层的发展．结果表明:来流脉冲扰动波与激波及边界层强烈相互作用,弓形激波明显向内弯曲,激波后扰动波被显著放大;来流扰动波与弓形激波干扰形成的边界层外的扰动波和近壁面内形成的边界层扰动波存在明显分界．钝锥头部参数扰动幅值要远大于其他位置参数扰动幅值．在边界层内的发展阶段,一些扰动模态持续增长,一些扰动模态被过滤掉,不再增长,甚至衰减,而也有一些扰动模态先衰减再增长．总的来说,在钝锥头部低频扰动模态为主导模态,随着扰动从流场上游向下游发展,总扰动模态中的低频模态成份和高频模态成份所占的比例开始转变,高频模态成分显著地增大．     

Dissociation effects in hypersonic flow past a circular cone at an angle of attack

The analytical solution to the steady, compressible, non-viscous, inviscid hypersonic flow past a circular cone at an angle of incidence, with an attached Shockwave, in the presence of dissociation of air in the shock layer, has been obtained here under the assumption of thermal equilibrium. Expression for the velocity, pressure, temperature, density, velocity of air, Mach number, pressure, drag and lift coefficients have been obtained both in the shocklayer outside the vortical layer and on the surface of the cone inside the vortical layer.     

On self-similar solutions of semilinear wave equations in higher space dimensions

In this paper we analyze self-similar solutions of the semilinear wave equation Φ_tt − ΔΦ − Φ^p = 0 for n > 3 space dimensions. We found several classes of analytic solutions labeled by a single parameter, the form of which differ in the vicinity of the light cone. We also propose suitable numerical methods to study them.     

On two perturbation estimates of the extreme solutions to the equations X ± A*XA = Q

Two perturbation estimates for maximal positive definite solutions of equations X + A*X⁻¹A = Q and X − A*X⁻¹A = Q are considered. These estimates are proved in [Hasanov et al., Improved perturbation Estimates for the Matrix Equations X ± A*X⁻¹A = Q, Linear Algebra Appl. 379 (2004) 113-135]. We derive new perturbation estimates under weaker restrictions on coefficient matrices of the equations. The theoretical results are illustrated by numerical examples.     

A smoothing self-adaptive Levenberg-Marquardt algorithm for solving system of nonlinear inequalities

In this paper, we consider the smoothing self-adaptive Levenberg-Marquardt algorithm for the system of nonlinear inequalities. By constructing a new smoothing function, the problem is approximated via a family of parameterized smooth equations H(x) = 0. A smoothing self-adaptive Levenberg-Marquardt algorithm is proposed for solving the system of nonlinear inequalities based on the new smoothing function. The Levenberg-Marquardt parameter μ_k is chosen as the product of μ_k = ∥H^k∥^δ with δ ∈ (0, 2] being a positive constant. We will show that if ∥H^k∥^δ provides a local error bound, which is weaker than the non-singularity, the proposed method converges superlinearly to the solution for δ ∈ (0, 1), while quadratically for δ ∈ [1, 2]. Numerical results show that the new method performs very well for system of inequalities.     

Resolution, in L-spaces, of transmission problems set in an unbounded domains

The goal of this work is to give a complete study of some abstract transmission problems (P^δ), for every δ > 0, set in unbounded domain composed of a half-line ]−∞, 0[ and a thin layer ]0, δ[. Existence and uniqueness results are obtained for strict solutions in UMD Banach spaces, by using essentially the semigroup theory and the Dore-Venni’s Theorem given in [8].     

Inventory models for perishable items with inventory level dependent demand rate

This paper presents inventory models for perishable items with inventory level dependent demand rate. The models with and without backlogging are studied. In the backlogging model, it is assumed that the backlogging rate is dependent on the waiting time and the amount of products already backlogged simultaneously. Two cases that holding inventory is profitable or not are studied, respectively. The smallest shelf space to ensure shortage not occur when holding inventory is not profitable is obtained. In the model without backlogging, it is assumed that the remaining stock at the end of the inventory cycle is disposed of with salvage value. The necessary and sufficient conditions for the existence and uniqueness of the optimal solution of these models are investigated. At last, some numerical examples are presented to illustrate the effectiveness of the proposed model. The model in this paper is generalization of present ones. In particularly, the model is reduced to Padmanabhan and Vrat’s when δ₁ = 0, and Dye and Ouyang’s when δ₂ = 0. If S = s and δ₂ = 0, it is Chang, Goyal and Teng’s model.     

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]).     

Special exact soliton solutions for the K(2, 2) equation with non-zero constant pedestal

Special exact solutions of the K(2, 2) equation, u_t + (u²)_x + (u²)_xxx = 0, are investigated by employing the qualitative theory of differential equations. Our procedure shows that the K(2, 2) equation either has loop soliton, cusped soliton and smooth soliton solutions when sitting on the non-zero constant pedestal lim_x→±∞u = A ≠ 0, or possesses compacton solutions only when lim_x→±∞u = 0. Mathematical analysis and numerical simulations are provided for these soliton solutions of the K(2, 2) equation.     

Krull’s theory for the double gamma function

The Barnes’ G-function G(x) = 1/Γ₂, satisfies the functional equation logG(x + 1) − logG(x) = logΓ(x). We complement W. Krull’s work in Bemerkungen zur Differenzengleichung g(x + 1) − g(x) = φ(x), Math. Nachrichten 1 (1948), 365-376 with additional results that yield a different characterization of the function G, new expansions and sharp bounds for G on x > 0 in terms of Gamma and Digamma functions, a new expansion for the Gamma function and summation formulae with Polygamma functions.     

The impact of smooth W-grids in the numerical solution of singular perturbation two-point boundary value problems

This paper develops a semi-analytic technique for generating smooth nonuniform grids for the numerical solution of singularly perturbed two-point boundary value problems. It is based on the usual idea of mapping a uniform grid to the desired nonuniform grid. We introduce the W-grid, which depends on the perturbation parameter ? ? 1. For problems on [0, 1] with a boundary layer at one end point, the local mesh width h_i = x_i+1 − x_i, with 0 = x₀ < x₁ < ? < x_N = 1, is condensed at either 0 or 1. Two simple 2nd order finite element and finite difference methods are combined with the new mesh, and computational experiments demonstrate the advantages of the smooth W-grid compared to the well-known piecewise uniform Shishkin mesh. For small ?, neither the finite difference method nor the finite element method produces satisfactory results on the Shishkin mesh. By contrast, accuracy is vastly improved on the W-grid, which typically produces the nominal 2nd order behavior in L², for large as well as small values of N, and over a wide range of values of ?. We conclude that the smoothness of the mesh is of crucial importance to accuracy, efficiency and robustness.     

Theory of the hypersonic viscous shock layer at high Reynolds numbers and intensive injection of foreign gases : PMM vol. 38, n 6, 1974, pp. 1015–1024

The hypersonic flow around smooth blunted bodies in the presence of intensive injection from the surface of these is considered. Using the method of external and internal expansions the asymptotics of the Navier-Stokes equations is constructed for high Reynolds numbers determined by parameters of the oncoming stream and of the injected gas. The flow in the shock layer falls into three characteristic regions. In regions adjacent to the body surface and the shock wave the effects associated with molecular transport are insignificant, while in the intermediate region they predominate. In the derivation of solution in the first two regions the surface of contact discontinuity is substituted for the region of molecular transport (external problem). An analytic solution of the external problem is obtained for small values of parameters ₁ = ρ_∞/ρs^* and δ = ρ_ω^*1/2ν_ω^*/ρ_∞1/2ν_∞, in the form of corresponding series expansions in these parameters. Asymptotic formulas are presented for velocity profiles, temperatures, and constituent concentration across the shock layer and, also, the shape of the contact discontinuity and of shock wave separation. The derived solution is compared with numerical solutions obtained by other authors. The flow in the region of molecular transport is defined by equations of the boundary layer with asymptotic conditions at plus and minus infinity, determined by the external solution (internal problem). A numerical solution of the internal problem is obtained taking into consideration multicomponent diffusion and heat exchange. The problem of multicomponent gas flow in the shock layer close to the stagnation line was previously considered in [1] with the use of simplified Navier-6tokes equations.The supersonic flow of a homogeneous inviscid and non-heat-conducting gas around blunted bodies in the presence of subsonic injection was considered in [2–7] using Euler's equations. An analytic solution, based on the classic solution obtained by Hill for a spherical vortex, was derived in [2] for a sphere on the assumption of constant but different densities in the layers between the shock wave and the contact discontinuity and between the latter and the body. Certain results of a numerical solution of the problem of intensive injection at the surface of axisymmetric bodies of various forms, obtained by Godunov's method [3], are presented. Telenin's method was used in [4] for numerical investigation of flow around a sphere; the problem was solved in two formulations: in the first, flow parameters were determined for the whole of the shock layer, while in the second this was done for the sutface of contact discontinuity, which was not known prior to the solution of the problem, with the pressure specified by Newton's formula and flow parameters determined only in the layer of injected gases. The flow with injection over blunted cones was numerically investigated in [5] by the approximate method proposed by Maslen. The flow in the shock layer in the neighborhood of the stagnation line was considered in [6, 8], and intensive injection was investigated by methods of the boundary layer theory in [8–12].     

Eigenvalues and the One-Dimensional p-Laplacian

We consider the boundary value problem (?_p(u′))′ + λF(t, u) = 0, with p > 1, t ∈ (0, 1), u(0) = u(1) = 0, and with λ > 0. The value of λ is chosen so that the boundary value problem has a positive solution. In addition, we derive an explicit interval for λ such that, for any λ in this interval, the existence of a positive solution to the boundary value problem is guaranteed. In addition, the existence of two positive solutions for λ in an appropriate interval is also discussed.     

Factoring and decomposing a class of linear functional systems

Within a constructive homological algebra approach, we study the factorization and decomposition problems for a class of linear functional (determined, over-determined, under-determined) systems. Using the concept of Ore algebras of functional operators (e.g., ordinary/partial differential operators, shift operators, time-delay operators), we first concentrate on the computation of morphisms from a finitely presented left module M over an Ore algebra to another one M′, where M (resp., M′) is a module intrinsically associated with the linear functional system Ry = 0 (resp., R′z = 0). These morphisms define applications sending solutions of the system R′z = 0 to solutions of R y = 0. We explicitly characterize the kernel, image, cokernel and coimage of a general morphism. We then show that the existence of a non-injective endomorphism of the module M is equivalent to the existence of a non-trivial factorization R = R₂R₁ of the system matrix R. The corresponding system can then be integrated “in cascade”. Under certain conditions, we also show that the system Ry = 0 is equivalent to a system R′z = 0, where R′ is a block-triangular matrix of the same size as R. We show that the existence of idempotents of the endomorphism ring of the module M allows us to reduce the integration of the system Ry = 0 to the integration of two independent systems R₁y₁ = 0 and R₂y₂ = 0. Furthermore, we prove that, under certain conditions, idempotents provide decompositions of the system Ry = 0, i.e., they allow us to compute an equivalent system R′z = 0, where R′ is a block-diagonal matrix of the same size as R. Applications of these results in mathematical physics and control theory are given. Finally, the different algorithms of the paper are implemented in the Maple package Morphisms based on the library oremodules.     

Peakon, loop and solitary travelling wave solutions for the general modified CH-DP equation

Travelling wave solutions for the general modified CH-DP equation u_t − u_xxt + αu²u_x − βu_xu_xx = uu_xxx are developed. By using the dynamical system method, a peakon and a dark soliton are found to coexist for the same wave speed. Exact explicit blow-up solutions are given. By using numerical simulation, a loop solution for a special case is discussed.     

Convergence rate of the solution toward boundary layer solution for initial-boundary value problem of the 2-D viscous conservation laws

In this paper, we study the large time behavior of the solution to the initial boundary value problem for 2-D viscous conservation laws in the space x ? bt. The global existence and the asymptotic stability of a stationary solution are proved by Kawashima et al. [1]. Here, we investigate the convergence rate of solution toward the boundary layer solution with the non-degenerate case where f′(u₊) − b < 0. Based on the estimate in the H² Sobolev space and via the weighted energy method, we draw the conclusion that the solution converges to the corresponding boundary layer solution with algebraic or exponential rate in time, under the assumption that the initial perturbation decays with algebraic or exponential in the spatial direction.     

Perturbation of the SVD in the presence of small singular values

This paper gives SVD perturbation bounds and expansions that are of use when an m × n, m ? n matrix A has small singular values. The first part of the paper gives subspace bounds that are closely related to those of Wedin but are stated so as to isolate the effect of any small singular values to the left singular subspace. In the second part first and second order approximations are given for perturbed singular values. The subspace bounds are used to show that all approximations retain accuracy when applied to small singular values. The paper concludes by deriving a subspace bound for multiplicative perturbations and using that bound to give a simple approximation to a singular value perturbed by a multiplicative perturbation.     

Multipartite Moore digraphs

We derive some Moore-like bounds for multipartite digraphs, which extend those of bipartite digraphs, under the assumption that every vertex of a given partite set is adjacent to the same number δ of vertices in each of the other independent sets. We determine when a multipartite Moore digraph is weakly distance-regular. Within this framework, some necessary conditions for the existence of a r-partite Moore digraph with interpartite outdegree δ > 1 and diameter k = 2m are obtained. In the case δ = 1, which corresponds to almost Moore digraphs, a necessary condition in terms of the permutation cycle structure is derived. Additionally, we present some constructions of dense multipartite digraphs of diameter two that are vertex-transitive.