Two splitting schemes for nonstationary convection-diffusion problems on tetrahedral meshes

Two splitting schemes are proposed for the numerical solution of three-dimensional nonstationary convection-diffusion problems on unstructured meshes in the case of a full diffusion tensor. An advantage of the first scheme is that splitting is generated by the properties of the approximation spaces and does not reduce the order of accuracy. An advantage of the second scheme is that the resulting numerical solutions are nonnegative. A numerical study is conducted to compare the splitting schemes with classical methods, such as finite elements and mixed finite elements. The numerical results show that the splitting schemes are characterized by low dissipation, high-order accuracy, and versatility.     

Mortar element methods for parabolic problems

In this article a standard mortar finite element method and a mortar element method with Lagrange multiplier are used for spatial discretization of a class of parabolic initial‐boundary value problems. Optimal error estimates in L^∞(L²) and L^∞(H¹)‐norms for semidiscrete methods for both the cases are established. The key feature that we have adopted here is to introduce a modified elliptic projection. In the standard mortar element method, a completely discrete scheme using backward Euler scheme is discussed and optimal error estimates are derived. The results of numerical experiments support the theoretical results obtained in this article. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

A new difference scheme with high accuracy and absolute stability for solving convection–diffusion equations

In this paper, we use a semi-discrete and a padé approximation method to propose a new difference scheme for solving convection–diffusion problems. The truncation error of the difference scheme is O(h⁴+τ⁵). It is shown through analysis that the scheme is unconditionally stable. Numerical experiments are conducted to test its high accuracy and to compare it with Crank–Nicolson method.     

An efficient ghost fluid method for compressible multifluids in Lagrangian coordinate

In this paper, we present an extended ghost fluid method (GFM) for computations of inviscid compressible multifluids in Lagrangian coordinate. That is, we capture the appropriate interface conditions by defining a fluid that has the velocity and the pressure of the real fluid at each point, but the entropy or the internal energy of some other fluid. Meanwhile, a single-fluid solver, CWENO-type central-upwind scheme, is developed in Lagrangian coordinate. The high resolution and the non-oscillatory quality of the scheme can be verified by solving several numerical experiments.     

On the reservoir technique convergence for nonlinear hyperbolic conservation laws - Part I

This paper is devoted to the analysis of flux schemes coupled with the reservoir technique for approximating hyperbolic equations and linear hyperbolic systems of conservation laws [F. Alouges, F. De Vuyst, G. Le Coq, E. Lorin, The reservoir scheme for systems of conservation laws, in: Finite Volumes for Complex Applications, III, Porquerolles, 2002, Lab. Anal. Topol. Probab. CNRS, Marseille, 2002, pp. 247-254 (electronic); F. Alouges, F. De Vuyst, G. Le Coq, E. Lorin, Un procédé de réduction de la diffusion numérique des schémas à différence de flux d'ordre un pour les systèmes hyperboliques non linéaires, C. R. Math. Acad. Sci. Paris 335 (7) (2002) 627-632; F. Alouges, F. De Vuyst, G. Le Coq, E. Lorin, The reservoir technique: A way to make Godunov-type schemes zero or very low diffusive. Application to Colella-Glaz, Eur. J. Mech. B Fluids 27 (6) (2008)]. We prove the long time convergence of the reservoir technique and its TVD property for some specific but still general configurations. Proofs are based on a precise study of the treatment by the reservoir technique of shock and rarefaction waves.     

曲线的一类非线性多分辨率算法及其收敛性

在曲线的多分辨率分析基础上,构造了一种新的非线性三分多分辨率算法.并研究这个正则三分多分辨率算法的收敛性和稳定性,进一步,证明了小波参数的收敛性精密地依靠这个基本的多分辨率细分算法的收敛性.     

Euler schemes and large deviations for stochastic Volterra equations with singular kernels

In this paper, an Euler type approximation is constructed for stochastic Volterra equation with singular kernels, which provides an algorithm for numerical calculation. Then, the large deviation estimates of small perturbation to equations of this type are obtained. We finally apply them to SDEs with the kernel of fractional Brownian motion with Hurst parameter H∈(0,1).     

Building nonredundant adaptive wavelets by update lifting

In this paper, we propose a technique for building adaptive wavelets by means of an extension of the lifting scheme. Our scheme comprises an adaptive update lifting step and a fixed prediction lifting step. The adaptivity consists hereof that the system can choose between two different update filters, and that this choice is triggered by the local gradient of the original signal. If the gradient is large (in some seminorm sense) it chooses one filter, if it is small the other. We derive necessary and sufficient conditions for the invertibility of such an adaptive system for various scenarios. Furthermore, we present some examples to illustrate our theoretical results.     

Convergence analysis of implicit iterative algorithms for asymptotically nonexpansive mappings

In this paper, we consider the weak and strong convergence of an implicit iterative process with errors for two finite families of asymptotically nonexpansive mappings in the framework of Banach spaces. Our results presented in this paper improve and extend the recent ones announced by many others.     

RCM‐TVD hybrid scheme for hyperbolic conservation laws

We describe a hybrid method for the solution of hyperbolic conservation laws. A third‐order total variation diminishing (TVD) finite difference scheme is conjugated with a random choice method (RCM) in a grid‐based adaptive way. An efficient multi‐resolution technique is used to detect the high gradient regions of the numerical solution in order to capture the shock with RCM while the smooth regions are computed with the more efficient TVD scheme. The hybrid scheme captures correctly the discontinuities of the solution and saves CPU time. Numerical experiments with one‐ and two‐dimensional problems are presented. Copyright © 2007 John Wiley & Sons, Ltd.