Search

In this work, we obtain quantitative estimates of the continuity constant for the $L^p$ maximal regularity of relatively continuous nonautonomous operators $\mathbb{A} : I \longrightarrow \mathcal{L}(D,X)$, where $D \subset X$ densely and compactly. They allow in particular to establish a new general growth condition for the global existence of strong solutions of Cauchy problems for nonlocal quasilinear equations for a certain class of nonlinearities $u \longrightarrow \mathbb{A}(u)$. The estimates obtained rely on the precise asymptotic analysis of the continuity constant with respect to perturbations of the operator of the form $\mathbb{A}(\cdot) + \lambda I$ as $\lambda \longrightarrow \pm \infty$. A complementary work in preparation supplements this abstract inquiry with an application of these results to nonlocal parabolic equations in noncylindrical domains depending on the time variable.

This manuscript deals with the analysis of numerical methods for the full discretization (in time and space) of the linear heat equation with Neumann boundary conditions, and it provides the reader with error estimates that are uniform in time. First, we consider the homogeneous equation with homogeneous Neumann boundary conditions over a finite interval. Using finite differences in space and the Euler method in time, we prove that our method is of order 1 in space, uniformly in time, under a classical CFL condition, and despite its lack of consistency at the boundaries. Second, we consider the nonhomogeneous equation with nonhomogeneous Neumann boundary conditions over a finite interval. Using a tailored similar scheme, we prove that our method is also of order 1 in space, uniformly in time, under a classical CFL condition. We indicate how this numerical method allows for a new way to compute steady states of such equations when they exist. We conclude by several numerical experiments to illustrate the sharpness and relevance of our theoretical results, as well as to examine situations that do not meet the hypotheses of our theoretical results, and to illustrate how our results extend to higher dimensions.

Starting from coupled fluid-kinetic equations for the modeling of laden flows, we derive relevant viscous corrections to be added to asymptotic hydrodynamic systems, by means of Chapman-Enskog expansions and analyse the shock profile structure for such limiting systems. Our main findings can be summarized as follows. Firstly, we consider simplified models, which are intended to reproduce the main difficulties and features of more intricate systems. However, while they are more easily accessible to analysis, such toy-models should be considered with caution since they might lose many important structural properties of the more realistic systems. Secondly, shock profiles can be identified also in such a case, which can be proven to be stable at least in the regime of small amplitude shocks. Last, but not least, regarding at the temperature of the mixture flow as a parameter of the problem, we show that the zero-temperature model admits viscous shock profiles. Numerical results indicate that a similar conclusion should apply in the regime of small positive temperatures.

In this article, we propose and study several discrete versions of homogeneous and inhomogeneous one-dimensional Fokker-Planck equations. In particular, for these discretizations of velocity and space, we prove the exponential convergence to the equilibrium of the solutions, for time-continuous equations as well as for time-discrete equations. Our method uses new types of discrete Poincar\'e inequalities for a "two-direction" discretization of the derivative in velocity. For the inhomogeneous problem, we adapt hypocoercive methods to the discrete cases.

This paper deals with the numerical integration of well-posed multiscale systems of ODEs or evolutionary PDEs. As these systems appear naturally in engineering problems, time-subcycling techniques are widely used every day to improve computational efficiency. These methods rely on a decomposition of the vector field in a fast part and a slow part and take advantage of that decomposition. This way, if an unconditionnally stable (semi-)implicit scheme cannot be easily implemented, one can integrate the fast equations with a much smaller time step than that of the slow equations, instead of having to integrate the whole system with a very small time-step to ensure stability. Then, one can build a numerical integrator using a standard composition method, such as a Lie or a Strang formula for example. Such methods are primarily designed to be convergent in short-time to the solution of the original problems. However, their longtime behavior rises interesting questions, the answers to which are not very well known. In particular, when the solutions of the problems converge in time to an asymptotic equilibrium state, the question of the asymptotic accuracy of the numerical longtime limit of the schemes as well as that of the rate of convergence is certainly of interest. In this context, the asymptotic error is defined as the difference between the exact and numerical asymptotic states. The goal of this paper is to apply that kind of numerical methods based on splitting schemes with subcycling to some simple examples of evolutionary ODEs and PDEs that have attractive equilibrium states, to address the aforementioned questions of asymptotic accuracy, to perform a rigorous analysis, and to compare them with their counterparts without subcycling. Our analysis is developed on simple linear ODE and PDE toy-models and is illustrated with several numerical experiments on these toy-models as well as on more complex systems. Lie and, Comment: IMA Journal of Numerical Analysis, Oxford University Press (OUP): Policy A - Oxford Open Option A, 2015

We present a general, high-order, fully explicit relaxation scheme which can be applied to any system of nonlinear hyperbolic conservation laws in multiple dimensions. The scheme consists of two steps. In a first (relaxation) step, the nonlinear hyperbolic conservation law is approximated by a kinetic equation with stiff BGK source term. Then, this kinetic equation is integrated in time using a projective integration method. After taking a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution, the time derivative is estimated and used in an (outer) Runge-Kutta method of arbitrary order. We show that, with an appropriate choice of inner step size, the time step restriction on the outer time step is similar to the CFL condition for the hyperbolic conservation law. Moreover, the number of inner time steps is also independent of the stiffness of the BGK source term. We discuss stability and consistency, and illustrate with numerical results (linear advection, Burgers' equation and the shallow water and Euler equations) in one and two spatial dimensions., Comment: 26 pages, 11 figures

We investigate a high-order, fully explicit, asymptotic-preserving scheme for a kinetic equation with linear relaxation, both in the hydrodynamic and diffusive scalings in which a hyperbolic, resp. parabolic, limiting equation exists. The scheme first takes a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution and estimate the time derivative of the slow components. These estimated time derivatives are then used in an (outer) Runge-Kutta method of arbitrary order. We show that, with an appropriate choice of inner step size, the time-step restriction on the outer time step is similar to the stability condition for the limiting macroscopic equation. Moreover, the number of inner time steps is also independent of the scaling parameter. We analyse stability and consistency, and illustrate with numerical results., Comment: 32 pages, 8 figures

We analyze numerically a forward-backward diffusion equation with a cubic-like diffusion function, -emerging in the framework of phase transitions modeling- and its "entropy" formulation determined by considering it as the singular limit of a third-order pseudo-parabolic equation. Precisely, we propose schemes for both the second and the third order equations, we discuss the analytical properties of their semi-discrete counter-parts and we compare the numerical results in the case of initial data of Riemann type, showing strengths and flaws of the two approaches, the main emphasis being given to the propagation of transition interfaces.

We investigate a projective integration scheme for a kinetic equation in the limit of vanishing mean free path, in which the kinetic description approaches a diffusion phenomenon. The scheme first takes a few small steps with a simple, explicit method, such as a spatial centered flux/forward Euler time integration, and subsequently projects the results forward in time over a large time step on the diffusion time scale. We show that, with an appropriate choice of the inner step size, the time-step restriction on the outer time step is similar to the stability condition for the diffusion equation, whereas the required number of inner steps does not depend on the mean free path. We also provide a consistency result. The presented method is asymptotic-preserving, in the sense that the method converges to a standard finite volume scheme for the diffusion equation in the limit of vanishing mean free path. The analysis is illustrated with numerical results, and we present an application to the Su-Olson test.