Asymptotic preserving Implicit-Explicit Runge-Kutta methods for non linear kinetic equations

We discuss Implicit-Explicit (IMEX) Runge Kutta methods which are particularly adapted to stiff kinetic equations of Boltzmann type. We consider both the case of easy invertible collision operators and the challenging case of Boltzmann collision operators. We give sufficient conditions in order that such methods are asymptotic preserving and asymptotically accurate. Their monotonicity properties are also studied. In the case of the Boltzmann operator, the methods are based on the introduction of a penalization technique for the collision integral. This reformulation of the collision operator permits to construct penalized IMEX schemes which work uniformly for a wide range of relaxation times avoiding the expensive implicit resolution of the collision operator. Finally we show some numerical results which confirm the theoretical analysis

Asymptotic-Preserving Monte Carlo methods for transport equations in the diffusive limit

We develop a new Monte Carlo method that solves hyperbolic transport equations with stiff terms, characterized by a (small) scaling parameter. In particular, we focus on systems which lead to a reduced problem of parabolic type in the limit when the scaling parameter tends to zero. Classical Monte Carlo methods suffer of severe time step limitations in these situations, due to the fact that the characteristic speeds go to infinity in the diffusion limit. This makes the problem a real challenge, since the scaling parameter may differ by several orders of magnitude in the domain. To circumvent these time step limitations, we construct a new, asymptotic-preserving Monte Carlo method that is stable independently of the scaling parameter and degenerates to a standard probabilistic approach for solving the limiting equation in the diffusion limit. The method uses an implicit time discretization to formulate a modified equation in which the characteristic speeds do not grow indefinitely when the scaling factor tends to zero. The resulting modified equation can readily be discretized by a Monte Carlo scheme, in which the particles combine a finite propagation speed with a time-step dependent diffusion term. We show the performance of the method by comparing it with standard (deterministic) approaches in the literature

Optimal control of epidemic spreading in presence of social heterogeneity

The spread of COVID-19 has been thwarted in most countries through non-pharmaceutical interventions. In particular, the most effective measures in this direction have been the stay-at-home and closure strategies of businesses and schools. However, population-wide lockdowns are far from being optimal carrying heavy economic consequences. Therefore, there is nowadays a strong interest in designing more efficient restrictions. In this work, starting from a recent kinetic-type model which takes into account the heterogeneity described by the social contact of individuals, we analyze the effects of introducing an optimal control strategy into the system, to limit selectively the mean number of contacts and reduce consequently the number of infected cases. Thanks to a data-driven approach, we show that this new mathematical model permits to assess the effects of the social limitations. Finally, using the model introduced here and starting from the available data, we show the effectivity of the proposed selective measures to dampen the epidemic trends

Analyses of shuttle orbiter approach and landing conditions

A study of one shuttle orbiter approach and landing conditions are summarized. Causes of observed PIO like flight deficiencies are identified and potential cures are examined. Closed loop pilot/vehicle analyses are described and path/attitude stability boundaries defined. The latter novel technique proved of great value in delineating and illustrating the basic causes of this multiloop pilot control problem. The analytical results are shown to be consistent with flight test and fixed base simulation. Conclusions are drawn relating to possible improvements of the shuttle orbiter/digital flight control system

The Moment Guided Monte Carlo method for the Boltzmann equation

In this work we propose a generalization of the Moment Guided Monte Carlo method developed in [11]. This approach permits to reduce the variance of the particle methods through a matching with a set of suitable macroscopic moment equations. In order to guarantee that the moment equations provide the correct solutions, they are coupled to the kinetic equation through a non equilibrium term. Here, at the contrary to the previous work in which we considered the simplified BGK operator, we deal with the full Boltzmann operator. Moreover, we introduce an hybrid setting which permits to entirely remove the resolution of the kinetic equation in the limit of infinite number of collisions and to consider only the solution of the compressible Euler equation. This modification additionally reduce the statistical error with respect to our previous work and permits to perform simulations of non equilibrium gases using only a few number of particles. We show at the end of the paper several numerical tests which prove the efficiency and the low level of numerical noise of the method.Comment: arXiv admin note: text overlap with arXiv:0908.026

Kinetic models for epidemic dynamics with social heterogeneity

We introduce a mathematical description of the impact of sociality in the spread of infectious diseases by integrating an epidemiological dynamics with a kinetic modeling of population-based contacts. The kinetic description leads to study the evolution over time of Boltzmann-type equations describing the number densities of social contacts of susceptible, infected and recovered individuals, whose proportions are driven by a classical SIR-type compartmental model in epidemiology. Explicit calculations show that the spread of the disease is closely related to moments of the contact distribution. Furthermore, the kinetic model allows to clarify how a selective control can be assumed to achieve a minimal lockdown strategy by only reducing individuals undergoing a very large number of daily contacts. We conduct numerical simulations which confirm the ability of the model to describe different phenomena characteristic of the rapid spread of an epidemic. Motivated by the COVID-19 pandemic, a last part is dedicated to fit numerical solutions of the proposed model with infection data coming from different European countries

A Multilevel Monte Carlo Asymptotic-Preserving Particle Method for Kinetic Equations in the Diffusion Limit

We propose a multilevel Monte Carlo method for a particle-based asymptotic-preserving scheme for kinetic equations. Kinetic equations model transport and collision of particles in a position-velocity phase-space. With a diffusive scaling, the kinetic equation converges to an advection-diffusion equation in the limit of zero mean free path. Classical particle-based techniques suffer from a strict time-step restriction to maintain stability in this limit. Asymptotic-preserving schemes provide a solution to this time step restriction, but introduce a first-order error in the time step size. We demonstrate how the multilevel Monte Carlo method can be used as a bias reduction technique to perform accurate simulations in the diffusive regime, while leveraging the reduced simulation cost given by the asymptotic-preserving scheme. We describe how to achieve the necessary correlation between simulation paths at different levels and demonstrate the potential of the approach via numerical experiments.Comment: 20 pages, 7 figures, published in Monte Carlo and Quasi-Monte Carlo Methods 2018, correction of minor typographical error

Asymptotically complexity diminishing schemes (ACDS) for kinetic equations in the diffusive scaling

In this work, we develop a new class of numerical schemes for collisional kinetic equations in the diffusive regime. The first step consists in reformulating the problem by decomposing the solution in the time evolution of an equilibrium state plus a perturbation. Then, the scheme combines a Monte Carlo solver for the perturbation with an Eulerian method for the equilibrium part, and is designed in such a way to be uniformly stable with respect to the diffusive scaling and to be consistent with the asymptotic diffusion equation. Moreover, since particles are only used to describe the perturbation part of the solution, the scheme becomes computationally less expensive – and is thus an asymptotically complexity diminishing scheme (ACDS) – as the solution approaches the equilibrium state due to the fact that the number of particles diminishes accordingly. This contrasts with standard methods for kinetic equations where the computational cost increases (or at least does not decrease) with the number of interactions. At the same time, the statistical error due to the Monte Carlo part of the solution decreases as the system approaches the equilibrium state: the method automatically degenerates to a solution of the macroscopic diffusion equation in the limit of infinite number of interactions. After a detailed description of the method, we perform several numerical tests and compare this new approach with classical numerical methods on various problems up to the full three dimensional case

A convenient approach to characterizing model uncertainty with application to early dark energy solutions of the Hubble tension

Despite increasingly precise observations and sophisticated theoretical models, the discrepancy between measurements of H0 from the cosmic microwave background or from Baryon Acoustic Oscillations combined with Big-Bang Nucleosynthesis versus those from local distance ladder probes -- commonly known as the $H_0$ tension -- continues to perplex the scientific community. To address this tension, Early Dark Energy (EDE) models have been proposed as alternatives to $\Lambda$CDM, as they can change the observed sound horizon and the inferred Hubble constant from measurements based on this. In this paper, we investigate the use of Bayesian Model Averaging (BMA) to evaluate EDE as a solution to the H0 tension. BMA consists of assigning a prior to the model and deriving a posterior as for any other unknown parameter in a Bayesian analysis. BMA can be computationally challenging in that one must approximate the joint posterior of both model and parameters. Here we present a computational strategy for BMA that exploits existing MCMC software and combines model-specific posteriors post-hoc. In application to a comprehensive analysis of cosmological datasets, we quantify the impact of EDE on the H0 discrepancy. We find an EDE model probability of $\sim$90% whenever we include the H0 measurement from Type Ia Supernovae in the analysis, whereas the other data show a strong preference for the standard cosmological model. We finally present constraints on common parameters marginalized over both cosmological models. For reasonable priors on models with and without EDE, the H0 tension is reduced by at least 20%