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 κ600;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     

Coupling for Markov renewal processes and the rate of convergence in ergodic theorems for processes with semi-Markov switchings

Coupling procedures for Markov renewal processes are described. Applications to ergodic theorems for processes with semi-Markov switchings are considered.This paper was partly prepared with the support of NFR Grant F-UP 10257-300.     

Nonstationary curves in Hilbert spaces and their correlation functions II

Partially supportted by Grant MM-22/1991 of MESC     

Asymptotics in the random assignment problem

Summary We show that, in the usual probabilistic model for the random assignment problem, the optimal cost tends to a limit constant in probability and in expectation. The method involves construction of an infinite limit structure, in terms of which the limit constant is defined. But we cannot improve on the known numerical bounds for the limit.Research supported by NSF Grant MCS90-01710     

U-statistics in danach spaces

U-statistics in Banach spaces are considered and thoroughly investigated. The martingale structure, estimates of moments, the law of large numbers, the central limit theorem, the invariance principle, estimates of the rate of convergence, and large deviations are established     

Mass spectral analysis and imaging of tissue by ToF-SIMS—The role of buckminsterfullerene, C60, primary ions

Recent developments in desorption/ionisation mass spectrometry techniques have made their application to biological analysis a realistic and successful proposition. Developments in primary ion source technology, mainly through the advent of polyatomic ion beams, have meant that the technique of secondary ion mass spectrometry (SIMS) can now access the depths of information required to allow biological imaging to be a viable option.Here the role of the primary ion C₆₀⁺ is assessed with regard to molecular imaging of lipids and pharmaceuticals within tissue sections. High secondary ion yields and low surface damage accumulation are demonstrated on both model and real biological samples, indicating the high secondary ion efficiency afforded to the analyst by this primary ion when compared to other cluster ion beams used in imaging. The newly developed 40 keV C₆₀⁺ ion source allows the beam to be focused such that high resolution imaging is demonstrated on a tissue sample, and the greater yields allow the molecular signal from the drug raclopride to be imaged within tissue section following in vivo dosing.The localisation shown for this drug alludes to issues regarding the chemical environment affecting the ionisation probability of the molecule; the importance of this effect is demonstrated with model systems and the possibility of using laser post-ionisation as a method for reducing this consequence of bio-sample complexity is demonstrated and discussed.