Well-posedness of second order evolution equation on discrete time

We characterize the well-posedness for second order discrete evolution equations in unconditional martingale difference spaces by means of Fourier multipliers and R-boundedness properties of the resolvent operator which defines the equation. Applications to semilinear problems are given.     

Error estimates for a finite element-finite volume discretization of convection-diffusion equations

We consider a time-dependent linear convection-diffusion equation. This equation is approximated by a combined finite element-finite volume method: the diffusion term is discretized by Crouzeix-Raviart piecewise linear finite elements, and the convection term by upwind barycentric finite volumes on a triangular grid. An implicit Euler approach is used for time discretization. It is shown that the error associated with this scheme, measured by a discrete L^∞-L²- and L²-H¹-norm, respectively, decays linearly with the mesh size and the time step. This result holds without any link between mesh size and time step. The dependence of the corresponding error bound on the diffusion coefficient is completely explicit.     

Construction par des éléments finis p1 d'une approximation convergente de problèmes variationnels

We consider some variational problems with nonlinear boundary conditions. We approximate them by some discrete problems using P₁ finite elements. We know that the discrete problems converge to the continuous ones when we consider linear boundary conditions (see [1]). We are interested by the nonlinear case. Mainly, we prove convergence for this kind of problems.     

On the complexity of inverse convex ordered 1-median problem on the plane and on tree networks

An ordered median function is used in location theory to generalize a class of problems, including median and center problems. In this paper we consider the complexity of inverse ordered 1-median problems on the plane and on trees, where the multipliers are sorted nondecreasingly. Based on the convexity of the objective function, we prove that the problems with variable weights or variable coordinates on the line are NP-hard. Then we can directly get the NP-hardness result for the corresponding problem on the plane. We finally develop a cubic time algorithm that solves the inverse convex ordered 1-median problem on trees with relaxation on modification bounds.     

Two-time scales in spatially structured models of population dynamics: A semigroup approach

The aim of this work is to provide a unified approach to the treatment of a class of spatially structured population dynamics models whose evolution processes occur at two different time scales. In the setting of the C₀-semigroup theory, we will consider a general formulation of some semilinear evolution problems defined on a Banach space in which the two-time scales are represented by a parameter ε>0 small enough, that mathematically gives rise to a singular perturbation problem. Applying the so-called aggregation of variables method, a simplified model called the aggregated model is constructed. A nontrivial mathematical task consists of comparing the asymptotic behaviour of solutions to both problems when ε→₀₊, under the assumption that the aggregated model has a compact attractor. Applications of the method to a class of two-time reaction-diffusion models of spatially structured population dynamics and to models with discrete spatial structure are given.     

On Fictitious Domain Formulations for Maxwell''s Equations

We consider fictitious domain-Lagrange multiplier formulations for variational problems in the space H(curl: Ω{\bf)} derived from Maxwell's equations. Boundary conditions and the divergence constraint are imposed weakly by using Lagrange multipliers. Both the time dependent and time harmonic formulations of the Maxwell's equations are considered, and we derive well-posed formulations for both cases. The variational problem that arises can be discretized by functions that do not satisfy an a-priori divergence constraint.     

Local Solutions of Weakly Parabolic Quasilinear Differential Equations

We consider a quasilinear parabolic boundary value problem, the elliptic part of which degenerates near the boundary. In order to solve this problem, we approximate it by a system of linear degenerate elliptic boundary value problems by means of semidiscretization with respect to time. We use the theory of degenerate elliptic operators and weighted Sobolev spaces to find a priori estimates for the solutions of the approximating problems. These solutions converge to a local solution, if the step size of the time-discretization goes to zero. It is worth pointing out that we do not require any growth conditions on the nonlinear coefficients and right-hand side, since we lire able to prove L∞ - estimates.     

Fourier decomposition of properly elliptic boundary value problems in a half-plane

For a general class of second-order elliptic boundary value problems in the lower half-plane, we show that the existence and uniqueness of solutions in L ^p Sobolev spaces is reduced to the invertibility of the ordinary differential operators obtained by Fourier decomposition. This terminology refers to the partial Fourier series expansion in the case of horizontally periodic solutions and to the partial Fourier transform otherwise. The problem is straightforward when p = 2 and, in the periodic case, the same question on a strip with finite width can also be quickly settled by indirect arguments irrespective of ${p \in (1, \infty )}$ . However, in the half-plane, the infinite depth raises serious difficulties when p ≠ 2. These difficulties are overcome by writing the problem as a first-order system and using existing abstract results about operator valued Fourier multipliers. In that approach, the randomized boundedness of the resolvent becomes the central issue.     

Nonlinear Fourier Transforms and the mKdV Equation in the Quarter Plane

The unified transform method introduced by Fokas can be used to analyze initial‐boundary value problems for integrable evolution equations. The method involves several steps, including the definition of spectral functions via nonlinear Fourier transforms and the formulation of a Riemann‐Hilbert problem. We provide a rigorous implementation of these steps in the case of the mKdV equation in the quarter plane under limited regularity and decay assumptions. We give detailed estimates for the relevant nonlinear Fourier transforms. Using the theory of L²‐RH problems, we consider the construction of quarter plane solutions which are C¹ in time and C³ in space.     

Analyse de Fourier,multiplicateurs de Schur sur Sp et ensembles Λ(p)cb non commutatifs

This work deals with various questions concerning Fourier multipliers on L^p, Schur multipliers on the Schatten class S^p as well as their completely bounded versions when L^p and S^p are viewed as operator spaces. We use for this aim subsets ofenjoying the Λ(p)_cb-property which is much stronger than the usual Λ(p)-property. We start by studying the notion of Λ(p)_cb-sets in the general case of an arbitrary discrete group before turning to .     

Birkhoff normal form for splitting methods applied to semilinear Hamiltonian PDEs. Part II. Abstract splitting

We consider Hamiltonian PDEs that can be split into a linear unbounded operator and a regular non linear part. We consider abstract splitting methods associated with this decomposition where no discretization in space is made. We prove a normal form result for the corresponding discrete flow under generic non resonance conditions on the frequencies of the linear operator and on the step size, and under a condition of zero momentum on the nonlinearity. This result implies the conservation of the regularity of the numerical solution associated with the splitting method over arbitrary long time, provided the initial data is small enough. This result holds for rounded numerical schemes avoiding at each step possible high frequency energy drift. We apply these results to nonlinear Schrödinger equations as well as the nonlinear wave equation.     

A theoretical analysis of a new second order finite volume approximation based on a low-order scheme using general admissible spatial meshes for the one dimensional wave equation

This paper details our note [6] and it is an extension of our previous works  and  which dealt with first order (both in time and space) and second order time accurate (second order in time and first order in space) implicit finite volume schemes for second order hyperbolic equations with Dirichlet boundary conditions on general nonconforming multidimensional spatial meshes introduced recently in [14]. We aim in this work (and some forthcoming studies) to get higher order (both in time and space) finite volume approximations for the exact solution of hyperbolic equations using the class of spatial generic meshes introduced recently in [14] on low order schemes from which the matrices used to compute the discrete solutions are sparse. We focus in the present contribution on the one dimensional wave equation and on one of its implicit finite volume schemes described in [4]. The implicit finite volume scheme approximating the one dimensional wave equation we consider (hereafter referred to as the basic finite volume scheme) yields linear systems to be solved successively. The matrices involved in these linear systems are tridiagonal, symmetric and definite positive. The finite volume approximate solution of the basic finite volume scheme is of order h+k$h+k$, where h (resp. k  ) is the mesh size of the spatial (resp. time) discretization. We construct a new finite volume approximation of order (h+k)²${(h+k)}^2$ in several discrete norms which allow us to get approximations of order two for the exact solution and its first derivatives. This new high-order approximation can be computed using linear systems whose matrices are the same ones used to compute the discrete solution of the basic finite volume scheme while the right hand sides are corrected. The construction of these right hand sides includes the approximation of some high order spatial derivatives of the exact solution. The computation of the approximation of these high order spatial derivatives can be performed using the same matrices stated above with another two tridiagonal matrices. The manner by which this new high-order approximation is constructed can be repeated to compute successively finite volume approximations of arbitrary order using the same matrices stated above. These high-order approximations can be obtained on any one dimensional admissible finite volume mesh in the sense of [13] without any condition. To reach the above results, a theoretical framework is developed and some numerical examples supporting the theory are presented. Some of the tools of this framework are new and interesting and they are stated in the one space dimension but they can be extended to several space dimensions. In particular a new and useful a prior estimate for a suitable discrete problem is developed and proved. The proof of this a prior estimate result is based essentially on the decomposition of the solution of the discrete problem into the solutions of two suitable discrete problems. A new technique is used in order to get a convenient finite volume approximation whose discrete time derivatives of order up to order two are also converging towards the solution of the wave equation and their corresponding time derivatives.     

Asymptotics of odd solutions for cubic nonlinear Schrödinger equations

We consider the Cauchy problem for a cubic nonlinear Schrödinger equation in the case of an odd initial data from H²∩H0,2. We prove the global existence in time of solutions to the Cauchy problem and construct the modified asymptotics for large values of time.     

Asymptotics of odd solutions of quadratic nonlinear Schrödinger equations

We consider the Cauchy problem for a quadratic nonlinear Schrödinger equation in the case of odd initial data from H²∩H0,2. We prove the global existence in time of solutions to the Cauchy problem and construct the modified asymptotics for large values of time.     

A class of Fourier multipliers on starlike Lipschitz surfaces

In this paper, we consider a class of Fourier multipliers whose symbols are controlled by a polynomial on starlike Lipschitz surfaces and get the L² boundedness of these operators on Sobolev spaces and their endpoint estimates.     

A Jacobi-collocation method for solving second kind Fredholm integral equations with weakly singular kernels

In this work, we propose a Jacobi-collocation method to solve the second kind linear Fredholm integral equations with weakly singular kernels. Particularly, we consider the case when the underlying solutions are sufficiently smooth. In this case, the proposed method leads to a fully discrete linear system. We show that the fully discrete integral operator is stable in both infinite and weighted square norms. Furthermore, we establish that the approximate solution arrives at an optimal convergence order under the two norms. Finally, we give some numerical examples, which confirm the theoretical prediction of the exponential rate of convergence.     

On the solution of discrete bottleneck problems

We consider a class of discrete bottleneck problems which includes bottleneck integral flow problems. It is demonstrated that a problem in this class, with k discrete variables, can be optimized by solving atmost k problems with real valued variables.     

A generalization of Chernoff's product formula for time-dependent operators

In this article we provide a set of sufficient conditions that allow a natural extension of Chernoff's product formula to the case of certain one-parameter family of functions taking values in the algebra L(B) of all bounded linear operators defined on a complex Banach space B. Those functions need not be contraction-valued and are intimately related to certain evolution operators U(t,s)_0?s?t?T on B. The most direct consequences of our main result are new formulae of the Trotter-Kato type which involve either semigroups with time-dependent generators, or the resolvent operators associated with these generators. In the general case we can apply such formulae to evolution problems of parabolic type, as well as to Schrödinger evolution equations albeit in some very special cases. The formulae we prove may also be relevant to the numerical analysis of non-autonomous ordinary and partial differential equations.     

An efficient numerical approach to solve Schrödinger equations with space fractional derivative

We design and analyze an efficient numerical approach to solve the coupled Schrödinger equations with space‐fractional derivative. The numerical scheme is based on leap‐frog in time direction and Fourier method in spatial direction. The advantage of the numerical scheme is that only a linear equation needs to be solved for each time step size, and we proved that the energy and mass of space‐fractional coupled Schrödinger equations (SFCSEs) are conserved in the case of full‐discrete scheme. Moreover, we also analyze the error estimate of the numerical scheme, and numerical solutions converge with the order in L² norm. Numerical examples are illustrated to verify the theoretical results.     

Fourier spectral projection method and nonlinear convergence analysis for Navier-Stokes equations

In this paper, we investigate the convergence rate of the Fourier spectral projection methods for the periodic problem of n-dimensional Navier-Stokes equations. Based on some alternative formulations of the Navier-Stokes equations and the related projection methods, the error estimates are carried out by a global nonlinear error analysis. It simplifies the analysis, relaxes the restriction on the time step size, weakens the regularity requirements on the genuine solution, and leads to some improved convergence results. A new correction technique is proposed for improving the accuracy of the numerical pressure.