Construction and Analysis of Lattice Boltzmann Methods Applied to a 1D Convection-Diffusion Equation

To solve the 1D (linear) convection-diffusion equation, we construct and we analyze two LBM schemes built on the D1Q2 lattice. We obtain these LBM schemes by showing that the 1D convection-diffusion equation is the fluid limit of a discrete velocity kinetic system. Then, we show in the periodic case that these LBM schemes are equivalent to a finite difference type scheme named LFCCDF scheme. This allows us, firstly, to prove the convergence in L ^∞ of these schemes, and to obtain discrete maximum principles for any time step in the case of the 1D diffusion equation with different boundary conditions. Secondly, this allows us to obtain most of these results for the Du Fort-Frankel scheme for a particular choice of the first iterate. We also underline that these LBM schemes can be applied to the (linear) advection equation and we obtain a stability result in L ^∞ under a classical CFL condition. Moreover, by proposing a probabilistic interpretation of these LBM schemes, we also obtain Monte-Carlo algorithms which approach the 1D (linear) diffusion equation. At last, we present numerical applications justifying these results.     

High‐order numerical methods for one‐dimensional parabolic singularly perturbed problems with regular layers

In this work we construct and analyze some finite difference schemes used to solve a class of time‐dependent one‐dimensional convection‐diffusion problems, which present only regular layers in their solution. We use the implicit Euler or the Crank‐Nicolson method to discretize the time variable and a HODIE finite difference scheme, defined on a piecewise uniform Shishkin mesh, to discretize the spatial variable. In both cases we prove that the numerical method is uniformly convergent with respect to the diffusion parameter, having order near two in space and order one or 3/2, depending on the method used, in time. We show some numerical examples which illustrate the theoretical results, in the case of using the Euler implicit method, and give better numerical behaviour than that predicted theoretically, showing order two in time and order N^?2log²N in space, if the Crank‐Nicolson scheme is used to discretize the time variable. Finally, we construct a numerical algorithm by combining a third order A‐stable SDIRK with two stages and a third‐order HODIE difference scheme, showing its uniformly convergent behavior, reaching order three, up to a logarithmic factor. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Extension of high‐order compact schemes to time‐dependent problems

In this article, we present an extension of our previous approaches for steady‐state higher‐order compact (HOC) difference methods to time‐dependent problems. The formulation also provides a framework for similar treatment of other HOC spatial schemes. A stability analysis is provided for transient convection‐diffusion in 1D and transient diffusion in 2D. Supporting numerical experiments are included to illustrate stability and accuracy as well as oscillatory and dissipative behavior. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 657–672, 2001     

The method of fundamental solutions with eigenfunction expansion method for nonhomogeneous diffusion equation

In this article we describe a numerical method to solve a nonhomogeneous diffusion equation with arbitrary geometry by combining the method of fundamental solutions (MFS), the method of particular solutions (MPS), and the eigenfunction expansion method (EEM). This forms a meshless numerical scheme of the MFS‐MPS‐EEM model to solve nonhomogeneous diffusion equations with time‐independent source terms and boundary conditions for any time and any shape. Nonhomogeneous diffusion equation with complex domain can be separated into a Poisson equation and a homogeneous diffusion equation using this model. The Poisson equation is solved by the MFS‐MPS model, in which the compactly supported radial basis functions are adopted for the MPS. On the other hand, utilizing the EEM the diffusion equation is first translated to a Helmholtz equation, which is then solved by the MFS together with the technique of the singular value decomposition (SVD). Since the present meshless method does not need mesh generation, nodal connectivity, or numerical integration, the computational effort and memory storage required are minimal as compared with other numerical schemes. Test results for two 2D diffusion problems show good comparability with the analytical solutions. The proposed algorithm is then extended to solve a problem with irregular domain and the results compare very well with solutions of a finite element scheme. Therefore, the present scheme has been proved to be very promising as a meshfree numerical method to solve nonhomogeneous diffusion equations with time‐independent source terms of any time frame, and for any arbitrary geometry. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A general mixed covolume framework for constructing conservative schemes for elliptic problems

We present a general framework for the finite volume or covolume schemes developed for second order elliptic problems in mixed form, i.e., written as first order systems. We connect these schemes to standard mixed finite element methods via a one-to-one transfer operator between trial and test spaces. In the nonsymmetric case (convection-diffusion equation) we show one-half order convergence rate for the flux variable which is approximated either by the lowest order Raviart-Thomas space or by its image in the space of discontinuous piecewise constants. In the symmetric case (diffusion equation) a first order convergence rate is obtained for both the state variable (e.g., concentration) and its flux. Numerical experiments are included.

     

The fast scalar auxiliary variable approach with unconditional energy stability for nonlocal Cahn–Hilliard equation

Comparing with the classical local gradient flow and phase field models, the nonlocal models such as nonlocal Cahn–Hilliard equations equipped with nonlocal diffusion operator can describe more practical phenomena for modeling phase transitions. In this paper, we construct an accurate and efficient scalar auxiliary variable approach for the nonlocal Cahn–Hilliard equation with general nonlinear potential. The first contribution is that we have proved the unconditional energy stability for nonlocal Cahn–Hilliard model and its semi‐discrete schemes carefully and rigorously. Second, what we need to focus on is that the nonlocality of the nonlocal diffusion term will lead the stiffness matrix to be almost full matrix which generates huge computational work and memory requirement. For spatial discretizaion by finite difference method, we find that the discretizaition for nonlocal operator will lead to a block‐Toeplitz–Toeplitz‐block matrix by applying four transformation operators. Based on this special structure, we present a fast procedure to reduce the computational work and memory requirement. Finally, several numerical simulations are demonstrated to verify the accuracy and efficiency of our proposed schemes.     

Upper bound for ruin probabilities under optimal investment and proportional reinsurance

In this paper, we consider the optimal investment and reinsurance from an insurer's point of view to maximize the adjustment coefficient. We obtain the explicit expressions for the optimal results in the diffusion approximation (D‐A) case as well as in the jump‐diffusion (J‐D) case. Furthermore, we derive a sharper bound on the ruin probability, from which we conclude that the case with investment is always better than the case without investment. Some numerical examples are presented to show that the ruin probability in the D‐A case sometimes underestimates the ruin probability in the J‐D case. Copyright © 2007 John Wiley & Sons, Ltd.     

A numerical scheme to the McKendrick–von Foerster equation with diffusion in age

In this paper a numerical scheme for McKendrick–von Foerster equation with diffusion in age (MV‐D) is proposed. First, we discretize the time variable to get a second‐order ordinary differential equation (ODE). At each time level, well‐posedness of this ODE is established using classical methods. Stability estimates for this semidiscrete scheme are derived. Later we construct piecewise linear (in time) functions using the solutions of the semidiscrete problems to approximate the solution to MV‐D and establish the convergence result. Numerical results are presented in some cases and compared with the corresponding analytic solutions where the latter is known explicitly.     

High‐order difference schemes for 2D elliptic and parabolic problems with mixed derivatives

We propose a 9‐point fourth‐order finite difference scheme for 2D elliptic problems with a mixed derivative and variable coefficients. The same approach is extended to derive a class of two‐level high‐order compact schemes with weighted time discretization for solving 2D parabolic problems with a mixed derivative. The schemes are fourth‐order accurate in space and second‐ or lower‐order accurate in time depending on the choice of a weighted average parameter μ. Unconditional stability is proved for 0.5 ≤ μ ≤ 1, and numerical experiments supporting our theoretical analysis and confirming the high‐order accuracy of the schemes are presented. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 366–378, 2007     

High‐order asymptotic‐preserving schemes for linear systems: Application to the Goldstein–Taylor equations

In this work, we address the numerical approximation of linear systems with possibly stiff source terms which induce an asymptotic diffusion limit. More precisely, we are interested in the design of high‐order asymptotic‐preserving schemes. Our approach is based on a very simple modification of the numerical flux associated with the usual HLL scheme. This alteration can be understood as a numerical diffusion reduction technique and allows to capture the correct asymptotic behavior in the diffusion limit and to consider uniformly high‐order extensions. We more specifically consider the case of the Goldstein–Taylor model but the overall approach is shown to be easily adapted to more general systems.     

High-order finite difference schemes for the solution of second-order BVPs

We introduce new methods in the class of boundary value methods (BVMs) to solve boundary value problems (BVPs) for a second-order ODE. These formulae correspond to the high-order generalizations of classical finite difference schemes for the first and second derivatives. In this research, we carry out the analysis of the conditioning and of the time-reversal symmetry of the discrete solution for a linear convection–diffusion ODE problem. We present numerical examples emphasizing the good convergence behavior of the new schemes. Finally, we show how these methods can be applied in several space dimensions on a uniform mesh.     

Further results on the impulsive synchronization of uncertain complex‐variable chaotic delayed systems

Synchronization and control of nonlinear dynamical systems with complex variables has attracted much more attention in various fields of science and engineering. In this article, we investigate the problem of impulsive synchronization for the complex‐variable delayed chaotic systems with parameters perturbation and unknown parameters in which the time delay is also included in the impulsive moment. Based on the theories of adaptive control and impulsive control, synchronization schemes are designed to make a class of complex‐variable chaotic delayed systems asymptotically synchronized, and unknown parameters are identified simultaneously in the process of synchronization. Sufficient conditions are derived to synchronize the complex‐variable chaotic systems include delayed impulses. To illustrate the effectiveness of the proposed schemes, several numerical examples are given. © 2014 Wiley Periodicals, Inc. Complexity 21: 131–142, 2016     

On the stability of alternating‐direction explicit methods for advection‐diffusion equations

Alternating‐Direction Explicit (A.D.E.) finite‐difference methods make use of two approximations that are implemented for computations proceeding in alternating directions, e.g., from left to right and from right to left, with each approximation being explicit in its respective direction of computation. Stable A.D.E. schemes for solving the linear parabolic partial differential equations that model heat diffusion are well‐known, as are stable A.D.E. schemes for solving the first‐order equations of fluid advection. Several of these are combined here to derive A.D.E. schemes for solving time‐dependent advection‐diffusion equations, and their stability characteristics are discussed. In each case, it is found that it is the advection term that limits the stability of the scheme. The most stable of the combinations presented comprises an unconditionally stable approximation for computations carried out in the direction of advection of the system, from left to right in this case, and a conditionally stable approximation for computations proceeding in the opposite direction. To illustrate the application of the methods and verify the stability conditions, they are applied to some quasi‐linear one‐dimensional advection‐diffusion problems. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Discrete artificial boundary conditions for the linearized Korteweg–de Vries equation

We consider the derivation of continuous and fully discrete artificial boundary conditions for the linearized Korteweg–de Vries equation. We show that we can obtain them for any constant velocities and any dispersion. The discrete artificial boundary conditions are provided for two different numerical schemes. In both continuous and discrete case, the boundary conditions are nonlocal with respect to time variable. We propose fast evaluations of discrete convolutions. We present various numerical tests which show the effectiveness of the artificial boundary conditions.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1455–1484, 2016     

Local discontinuous Galerkin methods with explicit Runge‐Kutta time marching for nonlinear carburizing model

A fully discrete local discontinuous Galerkin (LDG) method coupled with 3 total variation diminishing Runge‐Kutta time‐marching schemes, for solving a nonlinear carburizing model, will be analyzed and implemented in this paper. On the basis of a suitable numerical flux setting in the LDG method, we obtain the optimal error estimate for the Runge‐Kutta–LDG schemes by energy analysis, under the condition τ ≤ λh², where h and τ are mesh size and time step, respectively, λ is a positive constant independent of h. Numerical experiments are presented to verify the accuracy and capability of the proposed schemes. For the carburizing diffusion processes of steel and the diffusion simulation for Cu‐Ni system, the numerical results show good agreement with the experimental results.     

Stability of solutions of a class of nonlinear fractional diffusion equations with respect to a pseudo‐differential operator

In this paper, we study the forward and the backward in time problems for a class of nonlinear diffusion equations with respect to the pseudo‐differential operator. Herein, we investigate the stability of the solution of the forward problem in relationship with parameters of the pseudo‐differential operator and initial data. Besides, as known, the backward in time problem is instability. Hence, we give a method to regularize the solution of the backward problem in the case of the parameters are perturbed.     

Novel numerical techniques for the finite moment log stable computational model for European call option

Option pricing models are often used to describe the dynamic characteristics of prices in financial markets. Unlike the classical Black–Scholes (BS) model, the finite moment log stable (FMLS) model can explain large movements of prices during small time steps. In the FMLS, the second-order spatial derivative of the BS model is replaced by a fractional operator of order α which generates an α-stable Lévy process. In this paper, we consider the finite difference method to approximate the FMLS model. We present two numerical schemes for this approximation: the implicit numerical scheme and the Crank–Nicolson scheme. We carry out convergence and stability analyses for the proposed schemes. Since the fractional operator routinely generates dense matrices which often require high computational cost and storage memory, we explore three methods for solving the approximation schemes: the Gaussian elimination method, the bi-conjugate gradient stabilized method (Bi-CGSTAB) and the fast Bi-CGSTAB (FBi-CGSTAB) in order to compare the cost of calculations. Finally, two numerical examples with exact solutions are presented where we also use extrapolation techniques to achieve higher-order convergence. The results suggest that the proposed schemes are unconditionally stable and convergent, and the FMLS model is useful for pricing options.     

Arbitrary high-order linearly implicit energy-preserving algorithms for Hamiltonian PDEs

In this paper, we present a novel strategy to systematically construct linearly implicit energy-preserving schemes with arbitrary order of accuracy for Hamiltonian PDEs. Such a novel strategy is based on the newly developed exponential scalar auxiliary variable (ESAV) approach that can remove the bounded-from-below restriction of nonlinear terms in the Hamiltonian functional and provides a totally explicit discretization of the auxiliary variable without computing extra inner products. So it is more effective and applicable than the traditional scalar auxiliary variable (SAV) approach. To achieve arbitrary high-order accuracy and energy preservation, we utilize symplectic Runge-Kutta methods for both solution variables and the auxiliary variable, where the values of the internal stages in nonlinear terms are explicitly derived via an extrapolation from numerical solutions already obtained in the preceding calculation. A prediction-correction strategy is proposed to further improve the accuracy. Fourier pseudo-spectral method is then employed to obtain fully discrete schemes. Compared with the SAV schemes, the solution variables and the auxiliary variable in these ESAV schemes are now decoupled. Moreover, when the linear terms have constant coefficients, the solution variables can be explicitly solved by using the fast Fourier transform. Numerical experiments are carried out for three Hamiltonian PDEs to demonstrate the efficiency and conservation of the ESAV schemes.

     

A reduced‐order extrapolated natural boundary element method based on POD for the 2D hyperbolic equation in unbounded domain

In this article, we primarily focuses to study the order‐reduction for the classical natural boundary element (NBE) method for the two‐dimensional (2D) hyperbolic equation in unbounded domain. To this end, we first build a semi‐discretized format about time for the hyperbolic equation and discuss the existence, stability, and convergence of the time semi‐discretized solutions. We then establish the classical fully discretized NBE format from the time semi‐discretized one and analyze the existence, stability, and convergence of the classical NBE solutions. Next, using proper orthogonal decomposition method, we build a reduced‐order extrapolated NBE (ROENBE) format containing very few unknowns but having adequately high accuracy, and we also discuss the existence, stability, and convergence of the ROENBE solutions. Finally, we use some numerical examples to show that the ROENBE method is far superior to the classical NBE one. It shows that the ROENBE method is reliable and effective for solving the 2D hyperbolic equation with the unbounded domain.     

Supersonic turbulence as an agent of structure formation in space

Supersonic turbulence is a main agent for structure formation in space, from stellar winds to wind-blown bubbles, supernovae remnants, -ray bursts, and star-formation in molecular clouds. Dissipation is mostly due to shocks. Under isothermal conditions, the density follows a log-normal distribution whose variance grows with Mach number. The exact scaling is under debate. We focus on problems with not strictly periodic boundary conditions and present 2D and 3D numerical simulation of supersonic turbulence in plan-parallel slabs. From dimensional analysis we derive scaling laws for the mean quantities under isothermal conditions. Simulations including spatially extended radiative cooling layers reveal a complex interplay between turbulence and RM- and RT-unstable interfaces. Turbulent Mach numbers here are significantly reduced compared to the isothermal case and are of order 1. However, unstable interfaces between the different thermal component of the gas are equally able to produce structure. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)