A BGK‐penalization‐based asymptotic‐preserving scheme for the multispecies Boltzmann equation

An asymptotic‐preserving (AP) scheme is efficient in solving multiscale problems where kinetic and hydrodynamic regimes coexist. In this article, we extend the BGK‐penalization‐based AP scheme, originally introduced by Filbet and Jin for the single species Boltzmann equation (Filbet and Jin, J Comput Phys 229 (2010) 7625–7648), to its multispecies counterpart. For the multispecies Boltzmann equation, the new difficulties arise due to: (1) the breaking down of the conservation laws for each species and (2) different convergence rates to equilibria for different species in disparate masses systems. To resolve these issues, we find a suitable penalty function—the local Maxwellian that is based on the mean velocity and mean temperature and justify various asymptotic properties of this method. This AP scheme does not contain any nonlinear nonlocal implicit solver, yet it can capture the fluid dynamic limit with time step and mesh size independent of the Knudsen number. Numerical examples demonstrate the correct asymptotic‐behavior of the scheme. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Linear stability of WENO schemes coupled with explicit Runge–Kutta schemes

This paper focuses on the results of the linear stability analysis of the finite‐difference weighted essentially non‐oscillatory (WENO) schemes with optimal weights. The standard WENO schemes between the third and 11th order, the order‐optimised WENO schemes of the sixth and eighth order and the bandwidth‐optimised WENO schemes of the third and fourth order are considered. Several explicit Runge–Kutta schemes including the recently published strong stability‐preserving explicit Runge–Kutta schemes are considered for time discretisation. The stability limits as well as dissipation and dispersion properties dependent on the Courant–Friedrichs–Lewy number are presented for a hyperbolic model equation. The different combinations of space and time discretisation schemes are compared in terms of their accuracy and efficiency. For a parabolic model equation, the viscous term is discretised with high‐order central differences. The stability limits for the parabolic problem are presented as well. Numerical results of linear test cases are shown. Copyright © 2011 John Wiley & Sons, Ltd.     

Higher-order time integration with a local Lax–Wendroff procedure for a central scheme on an overlapping grid

We have implemented a high-order Lax–Wendroff type time integration for a central scheme on an overlapping grid for conservation law problems. Using a local iterative approach presented by Dumbser et al. (JCP, 2008) [12], we extend a local high-order spatial reconstruction on each cell to a local higher-order space–time polynomial on the cell. We rewrite the central scheme in a fully discrete form to avoid volume integration in the space–time domain. The fluxes at cell interfaces are calculated directly via integrating a higher-order space–time reconstruction of the flux. We compare this approach with the corresponding multi-stage Runge–Kutta time integration (RK). Numerical results show that the new time integration is more cost-effective.     

A thermodynamic and dynamic subgrid closure model for two‐material cells

A general and robust subgrid closure model for two‐material cells is proposed. The conservative quantities of the entire cell are apportioned between two materials, and then, pressure and velocity are fully or partially equilibrated by modeling subgrid wave interactions. An unconditionally stable and entropy‐satisfying solution of the processes has been successfully found. The solution is valid for arbitrary level of relaxation. The model is numerically designed with care for general materials and is computationally efficient without recourse to subgrid iterations or subcycling in time. The model is implemented and tested in the Lagrange‐remap framework. Two interesting results are observed in 1D tests. First, on the basis of the closure model without any pressure and velocity relaxation, a material interface can be resolved without creating numerical oscillations and/or large nonphysical jumps in the problem of the modified Sod shock tube. Second, the overheating problem seen near the wall surface can be solved by the present entropy‐satisfying closure model. The generality, robustness, and efficiency of the model make it useful in principle in algorithms, such as ALE methods, volume of fluid methods, and even some mixture models, for compressible two‐phase flow computations. Copyright © 2013 John Wiley & Sons, Ltd.     

A numerical analysis of a class of generalized Boussinesq‐type equations using continuous/discontinuous FEM

An improved class of Boussinesq systems of an arbitrary order using a wave surface elevation and velocity potential formulation is derived. Dissipative effects and wave generation due to a time‐dependent varying seabed are included. Thus, high‐order source functions are considered. For the reduction of the system order and maintenance of some dispersive characteristics of the higher‐order models, an extra O(μ²ⁿ⁺²) term (n ∈ ?) is included in the velocity potential expansion. We introduce a nonlocal continuous/discontinuous Galerkin FEM with inner penalty terms to calculate the numerical solutions of the improved fourth‐order models. The discretization of the spatial variables is made using continuous P₂ Lagrange elements. A predictor‐corrector scheme with an initialization given by an explicit Runge–Kutta method is also used for the time‐variable integration. Moreover, a CFL‐type condition is deduced for the linear problem with a constant bathymetry. To demonstrate the applicability of the model, we considered several test cases. Improved stability is achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

WENO schemes applied to the quasi-relativistic Vlasov–Maxwell model for laser–plasma interaction

In this paper we focus on WENO-based methods for the simulation of the 1D Quasi-Relativistic Vlasov–Maxwell (QRVM) model used to describe how a laser wave interacts with and heats a plasma by penetrating into it. We propose several non-oscillatory methods based on either Runge–Kutta (explicit) or Time-Splitting (implicit) time discretizations. We then show preliminary numerical experiments.     

Automorphisms of Hilbert schemes of points on a generic projective K3 surface

We study automorphisms of the Hilbert scheme of n points on a generic projective K3 surface S, for any . We show that is either trivial or generated by a non‐symplectic involution and we determine numerical and divisorial conditions which allow us to distinguish between the two cases. As an application of these results we prove that, for any , there exist infinitely many admissible degrees for the polarization of the K3 surface S such that admits a non‐natural involution. This provides a generalization of the results of [7] for .