Residual equilibrium schemes for time dependent partial differential equations

Many applications involve partial differential equations which admits nontrivial steady state solutions. The design of schemes which are able to describe correctly these equilibrium states may be challenging for numerical methods, in particular for high order ones. In this paper, inspired by micro-macro decomposition methods for kinetic equations, we present a class of schemes which are capable to preserve the steady state solution and achieve high order accuracy for a class of time dependent partial differential equations including nonlinear diffusion equations and kinetic equations. Extension to systems of conservation laws with source terms are also discussed.Comment: 23 pages, 12 figure

Asymptotic solutions of the nonlinear Boltzmann equation for dissipative systems

Analytic solutions $F(v,t)$ of the nonlinear Boltzmann equation in $d$-dimensions are studied for a new class of dissipative models, called inelastic repulsive scatterers, interacting through pseudo-power law repulsions, characterized by a strength parameter $\nu$, and embedding inelastic hard spheres ($\nu=1$) and inelastic Maxwell models ($\nu=0$). The systems are either freely cooling without energy input or driven by thermostats, e.g. white noise, and approach stable nonequilibrium steady states, or marginally stable homogeneous cooling states, where the data, $v^d_0(t) F(v,t)$ plotted versus $c=v/v_0(t)$, collapse on a scaling or similarity solution $f(c)$, where $v_0(t)$ is the r.m.s. velocity. The dissipative interactions generate overpopulated high energy tails, described generically by stretched Gaussians, $f(c) \sim \exp[-\beta c^b]$ with $0 < b < 2$, where $b=\nu$ with $\nu>0$ in free cooling, and $b=1+{1/2} \nu$ with $\nu \geq 0$ when driven by white noise. Power law tails, $f(c) \sim 1/c^{a+d}$, are only found in marginal cases, where the exponent $a$ is the root of a transcendental equation. The stability threshold depend on the type of thermostat, and is for the case of free cooling located at $\nu=0$. Moreover we analyze an inelastic BGK-type kinetic equation with an energy dependent collision frequency coupled to a thermostat, that captures all qualitative properties of the velocity distribution function in Maxwell models, as predicted by the full nonlinear Boltzmann equation, but fails for harder interactions with $\nu>0$.Comment: Submitted to: "Granular Gas Dynamics", T. Poeschel, N. Brilliantov (eds.), Lecture Notes in Physics, Vol. LNP 624, Springer-Verlag, Berlin-Heidelberg-New York, 200

Whitham Averaged Equations and Modulational Stability of Periodic Traveling Waves of a Hyperbolic-Parabolic Balance Law

In this note, we report on recent findings concerning the spectral and nonlinear stability of periodic traveling wave solutions of hyperbolic-parabolic systems of balance laws, as applied to the St. Venant equations of shallow water flow down an incline. We begin by introducing a natural set of spectral stability assumptions, motivated by considerations from the Whitham averaged equations, and outline the recent proof yielding nonlinear stability under these conditions. We then turn to an analytical and numerical investigation of the verification of these spectral stability assumptions. While spectral instability is shown analytically to hold in both the Hopf and homoclinic limits, our numerical studies indicates spectrally stable periodic solutions of intermediate period. A mechanism for this moderate-amplitude stabilization is proposed in terms of numerically observed "metastability" of the the limiting homoclinic orbits.Comment: 27 pages, 5 figures. Minor changes throughou

Asymptotic Preserving time-discretization of optimal control problems for the Goldstein-Taylor model

We consider the development of implicit-explicit time integration schemes for optimal control problems governed by the Goldstein-Taylor model. In the diffusive scaling this model is a hyperbolic approximation to the heat equation. We investigate the relation of time integration schemes and the formal Chapman-Enskog type limiting procedure. For the class of stiffly accurate implicit-explicit Runge-Kutta methods (IMEX) the discrete optimality system also provides a stable numerical method for optimal control problems governed by the heat equation. Numerical examples illustrate the expected behavior