A second-order accurate numerical method for a fractional wave equation

We study a generalized Crank–Nicolson scheme for the time discretization of a fractional wave equation, in combination with a space discretization by linear finite elements. The scheme uses a non-uniform grid in time to compensate for the singular behaviour of the exact solution at t = 0. With appropriate assumptions on the data and assuming that the spatial domain is convex or smooth, we show that the error is of order k ² + h ², where k and h are the parameters for the time and space meshes, respectively.     

Reliability and expectation bounds for coherent systems with exchangeable components

Sharp upper and lower bounds are obtained for the reliability functions and the expectations of lifetimes of coherent systems based on dependent exchangeable absolutely continuous components with a given marginal distribution function, by use of the concept of Samaniego's signature. We first show that the distribution of any coherent system based on exchangeable components with absolutely continuous joint distribution is a convex combination of distributions of order statistics (equivalent to the k-out-of-n systems) with the weights identical with the values of the Samaniego signature of the system. This extends the Samaniego representation valid for the case of independent and identically distributed components. Combining the representation with optimal bounds on linear combinations of distribution functions of order statistics from dependent identically distributed samples, we derive the corresponding reliability and expectation bounds, dependent on the signature of the system and marginal distribution of dependent components. We also present the sequences of exchangeable absolutely continuous joint distributions of components which attain the bounds in limit. As an application, we obtain the reliability bounds for all the coherent systems with three and four exchangeable components, expressed in terms of the parent marginal reliability function and specify the respective expectation bounds for exchangeable exponential components, comparing them with the lifetime expectations of systems with independent and identically distributed exponential components.     

The q-numerical radius of weighted shift operators with periodic weights

We deal with the q-numerical radius of weighted unilateral and bilateral shift operators. In particular, the q-numerical radius of weighted shift operators with periodic weights is discussed and computed.     

On the bidirectional vortex and other similarity solutions in spherical coordinates

The bidirectional vortex refers to the bipolar, coaxial swirling motion that can be triggered, for example, in cyclone separators and some liquid rocket engines with tangential aft-end injectors. In this study, we present an exact solution to describe the corresponding bulk motion in spherical coordinates. To do so, we examine both linear and nonlinear solutions of the momentum and vorticity transport equations in spherical coordinates. The assumption will be that of steady, incompressible, inviscid, rotational, and axisymmetric flow. We further relate the vorticity to some power of the stream function. At the outset, three possible types of similarity solutions are shown to fulfill the momentum equation. While the first type is incapable of satisfying the conditions for the bidirectional vortex, it can be used to accommodate other physical settings such as Hill’s vortex. This case is illustrated in the context of inviscid flow over a sphere. The second leads to a closed-form analytical expression that satisfies the boundary conditions for the bidirectional vortex in a straight cylinder. The third type is more general and provides multiple solutions. The spherical bidirectional vortex is derived using separation of variables and the method of variation of parameters. The three-pronged analysis presented here increases our repertoire of general mean flow solutions that rarely appear in spherical geometry. It is hoped that these special forms will permit extending the current approach to other complex fluid motions that are easier to capture using spherical coordinates.     

A fictitious domain approach to the numerical solution of PDEs in stochastic domains

We present an efficient method for the numerical realization of elliptic PDEs in domains depending on random variables. Domains are bounded, and have finite fluctuations. The key feature is the combination of a fictitious domain approach and a polynomial chaos expansion. The PDE is solved in a larger, fixed domain (the fictitious domain), with the original boundary condition enforced via a Lagrange multiplier acting on a random manifold inside the new domain. A (generalized) Wiener expansion is invoked to convert such a stochastic problem into a deterministic one, depending on an extra set of real variables (the stochastic variables). Discretization is accomplished by standard mixed finite elements in the physical variables and a Galerkin projection method with numerical integration (which coincides with a collocation scheme) in the stochastic variables. A stability and convergence analysis of the method, as well as numerical results, are provided. The convergence is “spectral” in the polynomial chaos order, in any subdomain which does not contain the random boundaries.     

Exponential decay rate of survival probability in a disastrous random environment

Summary. We study the exponential decay rate of the survival probability up to time t>0 of a random walker moving in Zopf; ^d in a temporally and spatially fluctuating random environment. When the random walker has a speed parameter κ>0, we investigate the influence of κ on the exponential decay rate λ(d,κ). In particular we prove that for any fixed d≥1, λ(d,κ) behaves like as logκ as κ↘0. Received: 21 May 1996 / In revised form: 2 February 1997     

Metric marginal problems for set-valued or non-measurable variables

Summary In a separable metric space, if two Borel probability measures (laws) are nearby in a suitable metric, then there exist random variables with those laws which are nearby in probability. Specifically, by a well-known theorem of Strassen, the Prohorov distance between two laws is the infimum of Ky Fan distances of random variables with those laws. The present paper considers possible extensions of Strassen's theorem to two random elements one of which may be (compact) set-valued and/or non-measurable. There are positive results in finite-dimensional spaces, but with factors depending on the dimension. Examples show that such factors cannot entirely be avoided, so that the extension of Strassen's theorem to the present situation fails in infinite dimensions.This research was partially supported by a Guggenheim Fellowship, by National Science Foundation grant DMS 8505550 at MSRI-Berkeley, and other NSF grants     

On order-preservation and positive correlations for multidimensional diffusion processes

Summary As a continuation of the study by Herbst and Pitt (1991), this note presents two criteria. The first one is on the order-preservation for two (may be different) multidimensional diffusion processes. The second one is on the preservation of positive correlations for a diffusion process.Research supported in part by the Ying-Tung Fok Educational Foundation and the National Natural Science Foundation of China     

A diffusion approximation result for two parameter processes

Summary We consider a one-dimensional linear wave equation with a small mean zero dissipative field and with the boundary condition imposed by the so-called Goursat problem. In order to observe the effect of the randomness on the solution we perform a space-time rescaling and we rewrite the problem in a diffusion approximation form for two parameter processes. We prove that the solution converges in distribution toward the solution of a two-parameter stochastic differential equation which we identify. The diffusion approximation results for oneparameter processes are well known and well understood. In fact, the solution of the one-parameter analog of the problem we consider here is immediate. Unfortunately, the situation is much more complicated for two-parameter processes and we believe that our result is the first one of its kind.Partially supported by ONR N00014-91-J-1010