A gas-kinetic scheme for shallow-water equations with source terms

In this paper, the Kinetic Flux Vector Splitting (KFVS) scheme is extended to solving the shallow water equations with source terms. To develop a well-balanced scheme between the source term and the flow convection, the source term effect is accounted in the flux evaluation across cell interfaces. This leads to a modified gas-kinetic scheme with particular application to the shallow water equations with bottom topography. Numerical experiments show better resolution of the unsteady solution than conventional finite difference method and KFVS method with little additional cost. Moreover, some positivity properties of the gas-kinetic scheme is established.     

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.     

Fonctions qui opèrent sur les espaces de Besov

Sans résumé

---

     

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.     

A p-version finite element method for nonlinear elliptic variational inequalities in 2D

This article introduces and analyzes a p-version FEM for variational inequalities resulting from obstacle problems for some quasi-linear elliptic partial differential operators. We approximate the solution by controlling the obstacle condition in images of the Gauss–Lobatto points. We show existence and uniqueness for the discrete solution u _p from the p-version for the obstacle problem. We prove the convergence of u _p towards the solution with respect to the energy norm, and assuming some additional regularity for the solution we derive an a priori error estimate. In numerical experiments the p-version turns out to be superior to the h-version concerning the convergence rate and the number of unknowns needed to achieve a certain exactness of the approximation.     

Asymptotic and boundary integral formulae for the computational reconstruction of small conductivity imperfections

In this paper we investigate the electrostatic problem of determining conductivity profiles from the knowledge of boundary currents and voltages. We obtain an improved estimate for the voltage potential of a two-dimensional conductor having finitely many circular inclusions and piecewise constant conductivity profile. We derive an asymptotic expansion for the voltage potential in terms of the reference voltage potential and the location, size, and conductivity of the inhomogeneities. This representation is used to formulate the nonlinear least squares problem for estimating the location and size of the inhomogeneities. Required boundary data for the voltage potential are generated numerically by solving a system of integral equations. Computational experiments are presented to demonstrate the effectiveness of our identification procedure.     

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