Reliability of a consecutive (r, s)-out-of-(m, n):F lattice system with conditions on the number of failed components in the system

A consecutive(r, s)-out-of-(m, n):F lattice system which is defined as a two-dimensional version of a consecutive k-out-of-n:F system is used as a reliability evaluation model for a sensor system, an X-ray diagnostic system, a pattern search system, etc. This system consists of m × n components arranged like an (m, n) matrix and fails iff the system has an (r, s) submatrix that contains all failed components. In this paper we deal a combined model of a k-out-of-mn:F and a consecutive (r, s)-out-of-(m, n):F lattice system. Namely, the system has one more condition of system down, that is the total number of failed components, in addition to that of a consecutive (r, s)-out-of-(m, n):F lattice system. We present a method to obtain reliability of the system. The proposed method obtains the reliability by using a combinatorial equation that does not depend on the system size. Some numerical examples are presented to show the relationship between component reliability and system reliability.     

On the Hermitian positive definite solutions of nonlinear matrix equation X + A∗XA = Q

Nonlinear matrix equation X^s + A∗X^−tA = Q, where A, Q are n × n complex matrices with Q Hermitian positive definite, has widely applied background. In this paper, we consider the Hermitian positive definite solutions of this matrix equation with two cases: s ? 1, 0 < t ? 1 and 0 < s ? 1, t ? 1. We derive necessary conditions and sufficient conditions for the existence of Hermitian positive definite solutions for the matrix equation and obtain some properties of the solutions. We also propose iterative methods for obtaining the extremal Hermitian positive definite solution of the matrix equation. Finally, we give some numerical examples to show the efficiency of the proposed iterative methods.     

Nordsieck methods with computationally verified algebraic stability

We describe the search for algebraically stable Nordsieck methods of order p = s and stage order q = p, where s is the number of stages. This search is based on the theoretical criteria for algebraic stability proposed recently by Hill, and Hewitt and Hill, for general linear methods for ordinary differential equations. These criteria, which are expressed in terms of the non-negativity of the eigenvalues of a Hermitian matrix on the unit circle, are then verified computationally for the derived Nordsieck methods of order p ? 2.     

Lost-sales inventory systems with a service level criterion

Competitive retail environments are characterized by service levels and lost sales in case of excess demand. We contribute to research on lost-sales models with a service level criterion in multiple ways. First, we study the optimal replenishment policy for this type of inventory system as well as base-stock policies and (R, s, S) policies. Furthermore, we derive lower and upper bounds on the order-up-to level, and we propose efficient approximation procedures to determine the order-up-to level. The procedures find values of the inventory control variables that are close to the best (R, s, S) policy and comply to the service level restriction for most of the instances, with an average cost increase of 2.3% and 1.2% for the case without and with fixed order costs, respectively.     

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.     

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.     

Three-term asymptotic expansions for the moments of the random walk with triangular distributed interference of chance

In this study, a semi-Markovian random walk with a discrete interference of chance (X(t)) is considered and under some weak assumptions the ergodicity of this process is discussed. The exact formulas for the first four moments of ergodic distribution of the process X(t) are obtained when the random variable ζ₁, which is describing a discrete interference of chance, has a triangular distribution in the interval [s, S] with center (S + s)/2. Based on these results, the asymptotic expansions with three-term are obtained for the first four moments of the ergodic distribution of X(t), as a ≡ (S − s)/2 → ∞. Furthermore, the asymptotic expansions for the variance, skewness and kurtosis of the ergodic distribution of the process X(t) are established. Finally, by using Monte Carlo experiments it is shown that the given approximating formulas provide high accuracy even for small values of parameter a.     

A simple approach to asymptotic expansions for Fourier integrals of singular functions

In this work, we are concerned with the derivation of full asymptotic expansions for Fourier integrals as s → ∞, where s is real positive, [a, b] is a finite interval, and the functions f(x) may have different types of algebraic and logarithmic singularities at x = a and x = b. This problem has been treated in the literature by techniques involving neutralizers and Mellin transforms. Here, we derive the relevant asymptotic expansions by a method that employs simpler and less sophisticated tools.     

Distributive equations of implications based on nilpotent triangular norms

In this paper, we explore the distributive equations of implications, both independently and along with other equations. In detail, we consider three classes of equations. (1) By means of the section of I, we give out the sufficient and necessary conditions of solutions for the distributive equation of implication I(x, T(y, z)) = T(I(x, y), (x, z)) based on a nilpotent triangular norm T and an unknown function I, which indicates that there are no continuous solutions satisfying the boundary conditions of implications. Under the assumptions that I is continuous except the vertical section I(0, y), y ∈ [0, 1), we get its complete characterizations. (2) We prove that there are no solutions for the functional equations I(x, T(y, z)) = T(I(x, y), I(x, z)), I(x, I(y, z)) = I(T(x, y), z). (3) We obtain the sufficient and necessary conditions on T and I to be solutions of the functional equations I(x, T(y, z)) = T(I(x, y), I(x, z)), I(x, y) = I(N(y), N(x)).     

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.     

Rooms with a view: A novel approach to iterated multidimensional wave conversion

Ray tracing in the presence of linear mode conversion leads to a ‘splitting’ of an incoming ray into two outgoing rays. When the rays are confined to a cavity, the rays can re-enter the conversion region many times, leading to iterated conversion. In this paper, we present new methods for the analysis of this problem. These involve a shift from local to global methods of analysis, and a shift in emphasis from the study of ray evolution in the dispersion surface to the study of the iterated dynamics of rays crossing the conversion surface. The analytical methods are quite general and can be applied in phase spaces of arbitrary dimension. In two spatial dimensions, (x, y), i.e. with a four-dimensional ray space, (x, y, k_x, k_y), rays are confined to three-dimensional regions called rooms, with one room for each wave type. In these rooms the rays do not cross, but when they intersect the conversion surface a family of converted rays is produced in the other room. The use of rooms allows a full view of the phase space dynamics of the iterated conversion of ray families. A simple two-dimensional model, inspired by the Budden resonance model, is presented as an example of these ideas.     

Critical Exponents of Quasilinear Parabolic Equations

In this paper we study the critical exponents of the Cauchy problem in Rⁿ of the quasilinear singular parabolic equations: u_t = div(|∇u|^m − 1∇u) + t^s|x|^σu^p, with non-negative initial data. Here s ≥ 0, (n − 1)/(n + 1) < m < 1, p > 1 and σ > n(1 − m) − (1 + m + 2s). We prove that p_c ≡ m + (1 + m + 2s + σ)/n > 1 is the critical exponent. That is, if 1 < p ≤ p_c then every non-trivial solution blows up in finite time, but for p > p_c, a small positive global solution exists.     

Composition operators from F(p, q, s) spaces to the nth weighted-type spaces on the unit disc

In this note we characterize the boundedness and compactness of the composition operator from the general function space F(p, q, s) to the nth weighted-type space on the unit disk, where the nth weighted-type space has been recently introduced by Stevo Stevi?.     

Modeling fully developed laminar flow in a helical duct with rectangular cross section and finite pitch

The mass and momentum transport equations are written in an orthogonal coordinate system using Germano’s transformation to model a laminar flow in a helical duct with a rectangular cross section and finite pitch. The system of governing equations are discretized and solved by the finite-volume numerical method. The three dimensional domain is reduced to a two dimensional slab of cells, orthogonal to the main flow direction, enforcing the fully developed state for 2π/(τ · d_h) >> 1 where τ and d_h representing the duct’s centerline torsion and its hydraulic diameter. This approximation and the use of an orthogonal grid allow a great simplification on the numerical procedure. Comparisons of the numerical solution against experimental data are drawn to assess the accuracy of the approximation.     

OMC for scalar multiples of unitaries

We establish the following case of the Determinantal Conjecture of Marcus [M. Marcus, Derivations, Plücker relations and the numerical range, Indiana Univ. Math. J. 22 (1973) 1137-1149] and de Oliveira [G.N. de Oliveira, Research problem: Normal matrices, Linear and Multilinear Algebra 12 (1982) 153-154]. Let A and B be unitary n × n matrices with prescribed eigenvalues a₁, … , a_n and b₁, … , b_n, respectively. Then for any scalars t and s

     

Exact solutions of nonlinear dispersive K(m, n) model with variable coefficients

The variable-coefficient Korteweg-de Vries (KdV) equation with additional terms contributed from the inhomogeneity in the axial direction and the strong transverse confinement of the condense was presented to describe the dynamics of nonlinear excitations in trapped quasi-one-dimensional Bose-Einstein condensates with repulsive atom-atom interactions. To understand the role of nonlinear dispersion in this variable-coefficient model, we introduce and study a new variable-coefficient KdV with nonlinear dispersion (called vc-K(m, n) equation). With the aid of symbolic computation, we obtain its compacton-like solutions and solitary pattern-like solutions. Moreover, we also present some conservation laws for both vc-K⁺(n, n) equation and vc-K⁻(n, n) equation.     

Rules for selecting the parameters of Oustaloup recursive approximation for the simulation of linear feedback systems containing PID controller

Oustaloup recursive approximation (ORA) is widely used to find a rational integer-order approximation for fractional-order integrators and differentiators of the form s^v, v ∈ (−1, 1). In this method the lower bound, the upper bound and the order of approximation should be determined beforehand, which is currently performed by trial and error and may be inefficient in some cases. The aim of this paper is to provide efficient rules for determining the suitable value of these parameters when a fractional-order PID controller is used in a stable linear feedback system. Two numerical examples are also presented to confirm the effectiveness of the proposed formulas.     

Fractal oscillations of self-adjoint and damped linear differential equations of second-order

For a prescribed real number s ∈ [1, 2), we give some sufficient conditions on the coefficients p(x) and q(x) such that every solution y = y(x), y ∈ C²((0, T]) of the linear differential equation (p(x)y′)′ + q(x)y = 0 on (0, T], is bounded and fractal oscillatory near x = 0 with the fractal dimension equal to s. This means that y oscillates near x = 0 and the fractal (box-counting) dimension of the graph Γ(y) of y is equal to s as well as the s dimensional upper Minkowski content (generalized length) of Γ(y) is finite and strictly positive. It verifies that y admits similar kind of the fractal geometric asymptotic behaviour near x = 0 like the chirp function y_ch(x) = a(x)S(φ(x)), which often occurs in the time-frequency analysis and its various applications. Furthermore, this kind of oscillations is established for the Bessel, chirp and other types of damped linear differential equations given in the form y″ + (μ/x)y′ + g(x)y = 0, x ∈ (0, T]. In order to prove the main results, we state a new criterion for fractal oscillations near x = 0 of real continuous functions which essentially improves related one presented in [1].     

A new off-step high order approximation for the solution of three-space dimensional nonlinear wave equations

In this paper, we propose a new high accuracy numerical method of O(k² + k²h² + h⁴) based on off-step discretization for the solution of 3-space dimensional non-linear wave equation of the form u_tt = A(x,y,z,t)u_xx + B(x,y,z,t)u_yy + C(x,y,z,t)u_zz + g(x,y,z,t,u,u_x,u_y,u_z,u_t), 0 < x,y,z < 1,t > 0 subject to given appropriate initial and Dirichlet boundary conditions, where k > 0 and h > 0 are mesh sizes in time and space directions respectively. We use only seven evaluations of the function g as compared to nine evaluations of the same function discussed in  and . We describe the derivation procedure in details of the algorithm. The proposed numerical algorithm is directly applicable to wave equation in polar coordinates and we do not require any fictitious points to discretize the differential equation. The proposed method when applied to a telegraphic equation is also shown to be unconditionally stable. Comparative numerical results are provided to justify the usefulness of the proposed method.     

Generalized Gaussian quadrature rules over two-dimensional regions with linear sides

This paper presents a generalized Gaussian quadrature method for numerical integration over triangular, parallelogram and quadrilateral elements with linear sides. In order to derive the quadrature rule, a general transformation of the regions, R₁ = {(x, y)∣a ? x ? b, g(x) ? y ? h(x)} and R₂ = {(x, y)∣a ? y ? b, g(y) ? x ? h(y)}, where g(x), h(x), g(y) and h(y) are linear functions, is given from (x, y) space to a square in (ξ, η) space, S: {(ξ, η)∣0 ? ξ ? 1, 0 ? η ? 1}. Generlized Gaussian quadrature nodes and weights introduced by Ma et.al. in 1997 are used in the product formula presented in this paper to evaluate the integral over S, as it is proved to give more accurate results than the classical Gauss Legendre nodes and weights. The method can be used to integrate a wide class of functions including smooth functions and functions with end-point singularities, over any two-dimensional region, bounded by linear sides. The performance of the method is illustrated for different functions over different two-dimensional regions with numerical examples.