A robust and efficient parallel SVD solver based on restarted Lanczos bidiagonalization.
We study a depth-averaged model of gravity-driven flows made of solid grains and fluid, moving over variable basal surface. In particular, we are interested in applications to geophysical flows such as avalanches and debris flows, which typically contain both solid material and interstitial fluid. The model system consists of mass and momentum balance equations for the solid and fluid components, coupled together by both conservative and non-conservative terms involving the derivatives of the...
This paper is concerned with the problem of determining the typical features of a curve when it is observed with noise. It has been shown that one can characterize the Lipschitz singularities of a signal by following the propagation across scales of the modulus maxima of its continuous wavelet transform. A nonparametric approach, based on appropriate thresholding of the empirical wavelet coefficients, is proposed to estimate the wavelet maxima of a signal observed with noise at various scales. In...
This paper is concerned with the problem of determining the typical features of a curve when it is observed with noise. It has been shown that one can characterize the Lipschitz singularities of a signal by following the propagation across scales of the modulus maxima of its continuous wavelet transform. A nonparametric approach, based on appropriate thresholding of the empirical wavelet coefficients, is proposed to estimate the wavelet maxima of a signal observed with noise at various scales....
We build a non-dissipative second order algorithm for the approximate resolution of the one-dimensional Euler system of compressible gas dynamics with two components. The considered model was proposed in [1]. The algorithm is based on [8] which deals with a non-dissipative first order resolution in Lagrange-remap formalism. In the present paper we describe, in the same framework, an algorithm that is second order accurate in time and space, and that...
In this paper, by use of affine biquadratic elements, we construct and analyze a finite volume element scheme for elliptic equations on quadrilateral meshes. The scheme is shown to be of second-order in -norm, provided that each quadrilateral in partition is almost a parallelogram. Numerical experiments are presented to confirm the usefulness and efficiency of the method.
In this paper, by use of affine biquadratic elements, we construct and analyze a finite volume element scheme for elliptic equations on quadrilateral meshes. The scheme is shown to be of second-order in H1-norm, provided that each quadrilateral in partition is almost a parallelogram. Numerical experiments are presented to confirm the usefulness and efficiency of the method.
The aim of this paper is to present a method using both the ideas of sectional approach and moment methods in order to accurately simulate evaporation phenomena in gas-droplets flows. Using the underlying kinetic interpretation of the sectional method [Y. Tambour, Combust. Flame 60 (1985) 15–28] exposed in [F. Laurent and M. Massot, Combust. Theory Model. 5 (2001) 537–572], we propose an extension of this approach based on a more accurate representation of the droplet size number density in each...
The aim of this paper is to present a method using both the ideas of sectional approach and moment methods in order to accurately simulate evaporation phenomena in gas-droplets flows. Using the underlying kinetic interpretation of the sectional method [Y. Tambour, Combust. Flame60 (1985) 15–28] exposed in [F. Laurent and M. Massot, Combust. Theory Model.5 (2001) 537–572], we propose an extension of this approach based on a more accurate representation of the droplet size number density in each...
By using some NCP functions, we reformulate the extended linear complementarity problem as a nonsmooth equation. Then we propose a self-adaptive trust region algorithm for solving this nonsmooth equation. The novelty of this method is that the trust region radius is controlled by the objective function value which can be adjusted automatically according to the algorithm. The global convergence is obtained under mild conditions and the local superlinear convergence rate is also established under...