This paper is mainly concerned with the global asymptotic behaviour of the unique solution to a class of singular Dirichlet problems − Δu = b(x)g(u), u > 0, x ∈ Ω, u|∂Ω = 0, where Ω is a bounded smooth domain in ℝn, g ∈ C1(0, ∞) is positive and decreasing in (0, ∞) with $\lim _{s\rightarrow 0^+}g(s)=\infty$, b ∈ Cα(Ω) for some α ∈ (0, 1), which is positive in Ω, but may vanish or blow up on the boundary properly. Moreover, we reveal the asymptotic behaviour of such a solution when the parameters on b tend to the corresponding critical values.

The effect of signals on stability, stable throughput region, and delay in a two-user slotted ALOHA-based random-access system with collisions is considered. This work gives rise to the development of random access G-networks, which can model security attacks, expiration of deadlines, or other malfunctions, and introduce load balancing among highly interacting queues. The users are equipped with infinite capacity buffers accepting external bursty arrivals. We consider both negative and triggering signals. Negative signals delete a packet from a user queue, while triggering signals cause the instantaneous transfer of packets among user queues. We obtain the exact stability region, and show that the stable throughput region is a subset of it. Moreover, we perform a compact mathematical analysis to obtain exact expressions for the queueing delay by solving a non-homogeneous Riemann boundary value problem. A computationally efficient way to obtain explicit bounds for the expected number of buffered packets at user queues is also presented. The theoretical findings are numerically evaluated and insights regarding the system performance are derived.

The present work is part of a series of studies conducted by the authors on analytical models of avascular tumour growth that exhibit both geometrical anisotropy and physical inhomogeneity. In particular, we consider a tumour structure formed in distinct ellipsoidal regions occupied by cell populations at a certain stage of their biological cycle. The cancer cells receive nutrient by diffusion from an inhomogeneous supply and they are subject to also an inhomogeneous pressure field imposed by the tumour microenvironment. It is proved that the lack of symmetry is strongly connected to a special condition that should hold between the data imposed by the tumour’s surrounding medium, in order for the ellipsoidal growth to be realizable, a feature already present in other non-symmetrical yet more degenerate models. The nutrient and the inhibitor concentration, as well as the pressure field, are provided in analytical fashion via closed-form series solutions in terms of ellipsoidal eigenfunctions, while their behaviour is demonstrated by indicative plots. The evolution equation of all the tumour’s ellipsoidal interfaces is postulated in ellipsoidal terms and a numerical implementation is provided in view of its solution. From the mathematical point of view, the ellipsoidal system is the most general coordinate system that the Laplace operator, which dominates the mathematical models of avascular growth, enjoys spectral decomposition. Therefore, we consider the ellipsoidal model presented in this work, as the most general analytic model describing the avascular growth in inhomogeneous environment. Additionally, due to the intrinsic degrees of freedom inherited to the model by the ellipsoidal geometry, the ellipsoidal model presented can be adapted to a very populous class of avascular tumours, varying in figure and in orientation.

We propose a hybrid spectral element method for fractional two-point boundary value problem (FBVPs) involving both Caputo and Riemann-Liouville (RL) fractional derivatives. We first formulate these FBVPs as a second kind Volterra integral equation (VIEs) with weakly singular kernel, following a similar procedure in [16]. We then design a hybrid spectral element method with generalized Jacobi functions and Legendre polynomials as basis functions. The use of generalized Jacobi functions allow us to deal with the usual singularity of solutions at t = 0. We establish the existence and uniqueness of the numerical solution, and derive a hptype error estimates under L2(I)-norm for the transformed VIEs. Numerical results are provided to show the effectiveness of the proposed methods.

We consider a single server system accepting two types of retrial customers, which arrive according to two independent Poisson streams. The service station can handle at most one customer, and in case of blocking, type i customer, i=1, 2, is routed to a separate type i orbit queue of infinite capacity. Customers from the orbits try to access the server according to the constant retrial policy. We consider coupled orbit queues, and thus, when both orbit queues are non-empty, the orbit queue i tries to re-dispatch a blocked customer of type i to the main service station after an exponentially distributed time with rate μi. If an orbit queue empties, the other orbit queue changes its re-dispatch rate from μi to $\mu_{i}^{\ast}$. We consider both exponential and arbitrary distributed service requirements, and show that the probability generating function of the joint stationary orbit queue length distribution can be determined using the theory of Riemann (–Hilbert) boundary value problems. For exponential service requirements, we also investigate the exact tail asymptotic behavior of the stationary joint probability distribution of the two orbits with either an idle or a busy server by using the kernel method. Performance metrics are obtained, computational issues are discussed and a simple numerical example is presented.

In this paper, Chebyshev series and rigorous numerics are combined to compute solutions of the Euler-Lagrange equations for the one-dimensional Ginzburg-Landau model of superconductivity. The idea is to recast solutions as fixed points of a Newton-like operator defined on a Banach space of rapidly decaying Chebyshev coefficients. Analytic estimates, the radii polynomials and the contraction mapping theorem are combined to show existence of solutions near numerical approximations. Coexistence of as many as seven nontrivial solutions is proved.

We present the formal geometric derivation of a non-equilibrium growth model that takes the form of a parabolic partial differential equation. Subsequently, we study its stationary radial solutions by means of variational techniques. Our results depend on the size of a parameter that plays the role of the strength of forcing. For small forcing we prove the existence and multiplicity of solutions to the elliptic problem. We discuss our results in the context of non-equilibrium statistical mechanics.

We define generalized Atiyah-Patodi-Singer boundary conditions of product type for Dirac operators associated to C*-vector bundles on the product of a compact manifold with boundary and a closed manifold. We prove a product formula for the K-theoretic index classes, which we use to generalize the product formula for the topological signature to higher signatures.

This paper is concerned with constructing polynomial solutions to ordinary boundary value problems. A semi-analytic technique using two-point Hermite interpolation is compared with conventional methods via a series of examples and is shown to be generally superior, particularly for problems involving nonlinear equations and/or boundary conditions.

We study the geometry of the set of closed extensions of index 0 of an elliptic differential cone operator and its model operator in connection with the spectra of the extensions, and we give a necessary and sufficient condition for the existence of rays of minimal growth for such operators.

We study the theoretical and numericalcoupling of two hyperbolic systems of conservation laws at a fixed interface. As already proven in the scalar case, the couplingpreserves in a weak sense the continuity of the solution at the interfacewithout imposing the overall conservativity of the coupled model. We develop a detailed analysis of the coupling inthe linear case. In the nonlinear case, we either use a linearized approach or a coupling method based on the solution of a Riemann problem. We discuss both approaches in the case of the coupling of two fluid models at a material contact discontinuity, the models being the usual gas dynamics equations with different equations ofstate. We also study the coupling of two-temperature plasma fluid models and illustrate the approach by numericalsimulations.

We derive expressions for the generating function of the equilibrium queue length probability distribution in a single server queue with general service times and independent Poisson arrival streams of both ordinary, positive customers and negative customers which eliminate a positive customer if present. For the case of first come first served queueing discipline for the positive customers, we compare the killing strategies in which either the last customer in the queue or the one in service is removed by a negative customer. We then consider preemptive-restart with resampling last come first served queueing discipline for the positive customers, combined with the elimination of the customer in service by a negative customer—the case of elimination of the last customer yields an analysis similar to first come first served discipline for positive customers. The results show different generating functions in contrast to the case where service times are exponentially distributed. This is also reflected in the stability conditions. Incidently, this leads to a full study of the preemptive-restart with resampling last come first served case without negative customers. Finally, approaches to solving the Fredholm integral equation of the first kind which arises, for instance, in the first case are considered as well as an alternative iterative solution method.

The Laplace transform of the probability distribution of the end-to-end delay in tandem networks is obtained where the first and/or second queue are G-queues, i.e. they have negative arrivals. For the most general case the method is based on the solution of a boundary value problem on a closed contour in the complex plane, which itself reduces to the solution of a Fredholm integral equation of the second kind. We also consider the dependence or independence of the sojourn times at each queue in the two special cases where only one of the queues is a G-queue, the other having no negative arrivals.