Finite volume methods for unidirectional dispersive wave models

We extend the framework of the finite volume method to dispersive unidirectional water wave propagation in one space dimension. In particular we consider a KdV-BBM type equation. Explicit and IMEX Runge-Kutta type methods are used for time discretizations. The fully discrete schemes are validated by direct comparisons to analytic solutions. Invariants conservation properties are also studied. Main applications include important nonlinear phenomena such as dispersive shock wave formation, solitary waves and their various interactions.Comment: 25 pages, 12 figures, 51 references. Other authors papers can be downloaded at http://www.lama.univ-savoie.fr/~dutykh

Decoupled and unidirectional asymptotic models for the propagation of internal waves

We study the relevance of various scalar equations, such as inviscid Burgers', Korteweg-de Vries (KdV), extended KdV, and higher order equations (of Camassa-Holm type), as asymptotic models for the propagation of internal waves in a two-fluid system. These scalar evolution equations may be justified with two approaches. The first method consists in approximating the flow with two decoupled, counterpropagating waves, each one satisfying such an equation. One also recovers homologous equations when focusing on a given direction of propagation, and seeking unidirectional approximate solutions. This second justification is more restrictive as for the admissible initial data, but yields greater accuracy. Additionally, we present several new coupled asymptotic models: a Green-Naghdi type model, its simplified version in the so-called Camassa-Holm regime, and a weakly decoupled model. All of the models are rigorously justified in the sense of consistency

Global regularity of two-dimensional flocking hydrodynamics

We study the systems of Euler equations which arise from agent-based dynamics driven by velocity \emph{alignment}. It is known that smooth solutions of such systems must flock, namely -- the large time behavior of the velocity field approaches a limiting "flocking" velocity. To address the question of global regularity, we derive sharp critical thresholds in the phase space of initial configuration which characterize the global regularity and hence flocking behavior of such \emph{two-dimensional} systems. Specifically, we prove for that a large class of \emph{sub-critical} initial conditions such that the initial divergence is "not too negative" and the initial spectral gap is "not too large", global regularity persists for all time

A Space-time Smooth Artificial Viscosity Method For Nonlinear Conservation Laws

We introduce a new methodology for adding localized, space-time smooth, artificial viscosity to nonlinear systems of conservation laws which propagate shock waves, rarefactions, and contact discontinuities, which we call the $C$-method. We shall focus our attention on the compressible Euler equations in one space dimension. The novel feature of our approach involves the coupling of a linear scalar reaction-diffusion equation to our system of conservation laws, whose solution $C(x,t)$ is the coefficient to an additional (and artificial) term added to the flux, which determines the location, localization, and strength of the artificial viscosity. Near shock discontinuities, $C(x,t)$ is large and localized, and transitions smoothly in space-time to zero away from discontinuities. Our approach is a provably convergent, spacetime-regularized variant of the original idea of Richtmeyer and Von Neumann, and is provided at the level of the PDE, thus allowing a host of numerical discretization schemes to be employed. We demonstrate the effectiveness of the $C$-method with three different numerical implementations and apply these to a collection of classical problems: the Sod shock-tube, the Osher-Shu shock-tube, the Woodward-Colella blast wave and the Leblanc shock-tube. First, we use a classical continuous finite-element implementation using second-order discretization in both space and time, FEM-C. Second, we use a simplified WENO scheme within our $C$-method framework, WENO-C. Third, we use WENO with the Lax-Friedrichs flux together with the $C$-equation, and call this WENO-LF-C. All three schemes yield higher-order discretization strategies, which provide sharp shock resolution with minimal overshoot and noise, and compare well with higher-order WENO schemes that employ approximate Riemann solvers, outperforming them for the difficult Leblanc shock tube experiment.Comment: 34 pages, 27 figure

Steady and Stable: Numerical Investigations of Nonlinear Partial Differential Equations

Excerpt: Mathematics is a language which can describe patterns in everyday life as well as abstract concepts existing only in our minds. Patterns exist in data, functions, and sets constructed around a common theme, but the most tangible patterns are visual. Visual demonstrations can help undergraduate students connect to abstract concepts in advanced mathematical courses. The study of partial differential equations, in particular, benefits from numerical analysis and simulation

Oscillating instanton solutions in curved space

We investigate oscillating instanton solutions of a self-gravitating scalar field between degenerate vacua. We show that there exist O(4)-symmetric oscillating solutions in a de Sitter background. The geometry of this solution is finite and preserves the $Z_{2}$ symmetry. The nontrivial solution corresponding to tunneling is possible only if the effect of gravity is taken into account. We present numerical solutions of this instanton, including the phase diagram of solutions in terms of the parameters of the present work and the variation of energy densities. Our solutions can be interpreted as solutions describing an instanton-induced domain wall or braneworld-like object rather than a kink-induced domain wall or braneworld. The oscillating instanton solutions have a thick wall and the solutions can be interpreted as a mechanism providing nucleation of the thick wall for topological inflation. We remark that $Z_{2}$ invariant solutions also exist in a flat and anti-de Sitter background, though the physical significance is not clear.Comment: 25 pages, 11 figues. Some typos corrected, references added, and Ch3. modified according to referee's comment

Convergence of the Crank-Nicolson-Galerkin finite element method for a class of nonlocal parabolic systems with moving boundaries

The aim of this paper is to establish the convergence and error bounds to the fully discrete solution for a class of nonlinear systems of reaction-diffusion nonlocal type with moving boundaries, using a linearized Crank-Nicolson-Galerkin finite element method with polynomial approximations of any degree. A coordinate transformation which fixes the boundaries is used. Some numerical tests to compare our Matlab code with some existing moving finite elements methods are investigated