No abstract uploaded!

In this paper, we present entropy satisfying schemes for solving an integro-differential equation that describes the evolution of a population structured with respect to a continuous trait. In [P.-E. Jabin and G. Raoul, J. Math. Biol., 63 (2011), pp. 493–517] solutions are shown to converge toward the so-called evolutionary stable distribution (ESD) as time becomes large, using the relative entropy. At the discrete level, the ESD is shown to be the solution to a quadratic programming problem and can be computed by any well-established nonlinear programing algorithm. The schemes are then shown to satisfy the entropy dissipation inequality on the set where initial data are positive and the numerical solutions tend toward the discrete ESD in time. An alternative algorithm is presented to capture the global ESD for nonnegative initial data, which is made possible due to the mutation mechanism built into the modified scheme. A series of numerical tests are given to confirm both accuracy and the entropy satisfying property and to underline the efficiency of capturing the large time asymptotic behavior of numerical solutions in various settings.

We propose efficient Eulerian methods for approximating the finite-time Lyapunov exponent (FTLE). The idea is to compute the related flow map using the Level Set Method and the Liouville equation. There are several advantages of the proposed approach. Unlike the usual Lagrangian-type computations, the resulting method requires the velocity field defined only at discrete locations. No interpolation of the velocity field is needed. Also, the method automatically stops a particle trajectory in the case when the ray hits the boundary of the computational domain. The computational complexity of the algorithm is <i>O</i>(<i>x</i>^(<i>d</i>+1)) with <i>d</i> the dimension of the physical space. Since there are the same number of mesh points in the <i>x</i><i>t</i> space, the computational complexity of the proposed Eulerian approach is optimal in the sense that each grid point is visited for only <i>O</i>(1) time. We also extend the algorithm to compute the FTLE

No abstract uploaded!

A novel discontinuous Galerkin (DG) method is developed to solve timedependent bi-harmonic type equations involving fourth derivatives in one and multiple space dimensions. We present the spatial DG discretization based on a mixed formulation and central interface numerical fluxes so that the resulting semi-discrete schemes are L^2 stable even without interior penalty. For time discretization, we use Crank–Nicolson so that the resulting scheme is unconditionally stable and second order in time. We present the optimal L^2 error estimate of O(h^{k+1}) for polynomials of degree k for semi-discrete DG schemes, and the L2 error of O(h^{k+1} + (\Delta t)^2) for fully discrete DG schemes. Extensions to more general fourth order partial differential equations and cases with non-homogeneous boundary conditions are provided. Numerical results are presented to verify the stability and accuracy of the schemes. Finally, an application to the one-dimensional Swift–Hohenberg equation endowed with a decay free energy is presented.