Adaptive multiresolution computations applied to detonations

A space-time adaptive method is presented for the reactive Euler equations describing chemically reacting gas flow where a two species model is used for the chemistry. The governing equations are discretized with a finite volume method and dynamic space adaptivity is introduced using multiresolution analysis. A time splitting method of Strang is applied to be able to consider stiff problems while keeping the method explicit. For time adaptivity an improved Runge--Kutta--Fehlberg scheme is used. Applications deal with detonation problems in one and two space dimensions. A comparison of the adaptive scheme with reference computations on a regular grid allow to assess the accuracy and the computational efficiency, in terms of CPU time and memory requirements.Comment: Zeitschrift f\"ur Physicalische Chemie, accepte

Finite Volume vs.vs. Streaming-based Lattice Boltzmann algorithm for fluid-dynamics simulations: a one-to-one accuracy and performance study

A new finite volume (FV) discretisation method for the Lattice Boltzmann (LB) equation which combines high accuracy with limited computational cost is presented. In order to assess the performance of the FV method we carry out a systematic comparison, focused on accuracy and computational performances, with the standard $streaming$ (ST) Lattice Boltzmann equation algorithm. To our knowledge such a systematic comparison has never been previously reported. In particular we aim at clarifying whether and in which conditions the proposed algorithm, and more generally any FV algorithm, can be taken as the method of choice in fluid-dynamics LB simulations. For this reason the comparative analysis is further extended to the case of realistic flows, in particular thermally driven flows in turbulent conditions. We report the first successful simulation of high-Rayleigh number convective flow performed by a Lattice Boltzmann FV based algorithm with wall grid refinement.Comment: 15 pages, 14 figures (discussion changes, improved figure readability

An adaptive, fully implicit multigrid phase-field model for the quantitative simulation of non-isothermal binary alloy solidification

Using state-of-the-art numerical techniques, such as mesh adaptivity, implicit time-stepping and a non-linear multi-grid solver, the phase-field equations for the non-isothermal solidification of a dilute binary alloy have been solved. Using the quantitative, thin-interface formulation of the problem we have found that at high Lewis number a minimum in the dendrite tip radius is predicted with increasing undercooling, as predicted by marginal stability theory. Over the dimensionless undercooling range 0.2–0.8 the radius selection parameter, σ*, was observed to vary by over a factor of 2 and in a non-monotonic fashion, despite the anisotropy strength being constant