An almost fourth order uniformly convergent difference scheme for a semilinear singularly perturbed reaction-diffusion problem

Summary. This paper is concerned with a high order convergent discretization for the semilinear reaction-diffusion problem: , for , subject to , where . We assume that on , which guarantees uniqueness of a solution to the problem. Asymptotic properties of this solution are discussed. We consider a polynomial-based three-point difference scheme on a simple piecewise equidistant mesh of Shishkin type. Existence and local uniqueness of a solution to the scheme are analysed. We prove that the scheme is almost fourth order accurate in the discrete maximum norm, uniformly in the perturbation parameter . We present numerical results in support of this result. Received February 25, 1994     

Thermal degradation of two liquid fuels and detonation tests for pulse detonation engine studies

The use of liquid fuels such as kerosene is of interest for the pulse detonation engine (PDE). Within this context, the aim of this work, which is a preliminary study, was to show the feasibility to initiate a detonation in air with liquid-fuel pyrolysis products, using energies and dimensions of test facility similars to those of PDEs. Therefore, two liquids fuels have been compared, JP10, which is a synthesis fuel generally used in the field of missile applications, and decane, which is one of the major components of standard kerosenes (F-34, Jet A1, ...). The thermal degradation of these fuels was studied with two pyrolysis processes, a batch reactor and a flow reactor. The temperatures varied from 600°C to 1,000°C and residence times for the batch reactor and the flow reactor were, respectively, between 10–30 s and 0.1–2 s. Subsequently, the detonability of synthetic gaseous mixtures, which was a schematisation of the decomposition state after the pyrolysis process, has been studied. The detonability study, regarding nitrogen dilution and equivalence ratio, was investigated in a 50 mm-diameter, 2.5 m-long detonation tube. These dimensions are compatible with applications in the aircraft industry and, more particularly, in PDEs. Therefore, JP10 and decane were compared to choose the best candidate for liquid-fuel PDE studies. This paper was based on work that was presented at the 20th International Colloquium on the Dynamics of Explosions and Reactive Systems, Montreal, Canada, July 31 – August 5, 2005.     

A robust and efficient adaptive reweighted estimator of multivariate location and scatter

This article proposes a reweighted estimator of multivariate location and scatter, with weights adaptively computed from the data. Its breakdown point and asymptotic behavior under elliptical distributions are established. This adaptive estimator is able to attain simultaneously the maximum possible breakdown point for affine equivariant estimators and full asymptotic efficiency at the multivariate normal distribution. For the special case of hard-rejection weights and the MCD as initial estimator, it is shown to be more efficient than its non-adaptive counterpart for a broad range of heavy-tailed elliptical distributions. A Monte Carlo study shows that the adaptive estimator is as robust as its non-adaptive relative for several types of bias-inducing contaminations, while it is remarkably more efficient under normality for sample sizes as small as 200.     

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.     

Existence of positive solutions for nonlinear three-point problems on time scales

In this paper, by using fixed point theorems in cones, we study the existence of at least one, two and three positive solutions of a nonlinear second-order three-point boundary value problem for dynamic equations on time scales. As an application, we also give some examples to demonstrate our results.     

Semi-implicit complementary volume scheme for solving level set like equations in image processing and curve evolution

Summary. We introduce linear semi-implicit complementary volume numerical scheme for solving level set like nonlinear degenerate diffusion equations arising in image processing and curve evolution problems. We study discretization of image selective smoothing equation of mean curvature flow type given by Alvarez, Lions and Morel ([3]). Solution of the level set equation of Osher and Sethian ([26], \[30]) is also included in the study. We prove and estimates for the proposed scheme and give existence of its (generalized) solution in every discrete time-scale step. Efficiency of the scheme is given by its linearity and stability. Preconditioned iterative solvers are used for computing arising linear systems. We present computational results related to image processing and plane curve evolution. Received April 25, 2000 / Revised version received June 11, 2001 / Published online November 15, 2001     

HAC estimation and strong linearity testing in weak ARMA models

In the framework of ARMA models, we consider testing the reliability of the standard asymptotic covariance matrix (ACM) of the least-squares estimator. The standard formula for this ACM is derived under the assumption that the errors are independent and identically distributed, and is in general invalid when the errors are only uncorrelated. The test statistic is based on the difference between a conventional estimator of the ACM of the least-squares estimator of the ARMA coefficients and its robust HAC-type version. The asymptotic distribution of the HAC estimator is established under the null hypothesis of independence, and under a large class of alternatives. The asymptotic distribution of the proposed statistic is shown to be a standard χ² under the null, and a noncentral χ² under the alternatives. The choice of the HAC estimator is discussed through asymptotic power comparisons. The finite sample properties of the test are analyzed via Monte Carlo simulation.     

Runge–Kutta convolution quadrature methods for well-posed equations with memory

Runge–Kutta based convolution quadrature methods for abstract, well-posed, linear, and homogeneous Volterra equations, non necessarily of sectorial type, are developed. A general representation of the numerical solution in terms of the continuous one is given. The error and stability analysis is based on this representation, which, for the particular case of the backward Euler method, also shows that the numerical solution inherits some interesting qualitative properties, such as positivity, of the exact solution. Numerical illustrations are provided.     

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.