MFC: An open-source high-order multi-component, multi-phase, and multi-scale compressible flow solver

MFC is an open-source tool for solving multi-component, multi-phase, and bubbly compressible flows. It is capable of efficiently solving a wide range of flows, including droplet atomization, shock–bubble interaction, and bubble dynamics. We present the 5- and 6-equation thermodynamically-consistent diffuse-interface models we use to handle such flows, which are coupled to high-order interface-capturing methods, HLL-type Riemann solvers, and TVD time-integration schemes that are capable of simulating unsteady flows with strong shocks. The numerical methods are implemented in a flexible, modular framework that is amenable to future development. The methods we employ are validated via comparisons to experimental results for shock–bubble, shock–droplet, and shock–water-cylinder interaction problems and verified to be free of spurious oscillations for material-interface advection and gas–liquid Riemann problems. For smooth solutions, such as the advection of an isentropic vortex, the methods are verified to be high-order accurate. Illustrative examples involving shock–bubble-vessel-wall and acoustic–bubble-net interactions are used to demonstrate the full capabilities of MFC

A limiter-based well-balanced discontinuous Galerkin method for shallow-water flows with wetting and drying: Triangular grids

A novel wetting and drying treatment for second-order Runge-Kutta discontinuous Galerkin (RKDG2) methods solving the non-linear shallow water equations is proposed. It is developed for general conforming two-dimensional triangular meshes and utilizes a slope limiting strategy to accurately model inundation. The method features a non-destructive limiter, which concurrently meets the requirements for linear stability and wetting and drying. It further combines existing approaches for positivity preservation and well-balancing with an innovative velocity-based limiting of the momentum. This limiting controls spurious velocities in the vicinity of the wet/dry interface. It leads to a computationally stable and robust scheme -- even on unstructured grids -- and allows for large time steps in combination with explicit time integrators. The scheme comprises only one free parameter, to which it is not sensitive in terms of stability. A number of numerical test cases, ranging from analytical tests to near-realistic laboratory benchmarks, demonstrate the performance of the method for inundation applications. In particular, super-linear convergence, mass-conservation, well-balancedness, and stability are verified

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

Thermoacoustic instability - a dynamical system and time domain analysis

This study focuses on the Rijke tube problem, which includes features relevant to the modeling of thermoacoustic coupling in reactive flows: a compact acoustic source, an empirical model for the heat source, and nonlinearities. This thermo-acoustic system features a complex dynamical behavior. In order to synthesize accurate time-series, we tackle this problem from a numerical point-of-view, and start by proposing a dedicated solver designed for dealing with the underlying stiffness, in particular, the retarded time and the discontinuity at the location of the heat source. Stability analysis is performed on the limit of low-amplitude disturbances by means of the projection method proposed by Jarlebring (2008), which alleviates the linearization with respect to the retarded time. The results are then compared to the analytical solution of the undamped system, and to Galerkin projection methods commonly used in this setting. This analysis provides insight into the consequences of the various assumptions and simplifications that justify the use of Galerkin expansions based on the eigenmodes of the unheated resonator. We illustrate that due to the presence of a discontinuity in the spatial domain, the eigenmodes in the heated case, predicted by using Galerkin expansion, show spurious oscillations resulting from the Gibbs phenomenon. By comparing the modes of the linear to that of the nonlinear regime, we are able to illustrate the mean-flow modulation and frequency switching. Finally, time-series in the fully nonlinear regime, where a limit cycle is established, are analyzed and dominant modes are extracted. The analysis of the saturated limit cycles shows the presence of higher frequency modes, which are linearly stable but become significant through nonlinear growth of the signal. This bimodal effect is not captured when the coupling between different frequencies is not accounted for.Comment: Submitted to Journal of Fluid Mechanic

Nonlinear physics of electrical wave propagation in the heart: a review

The beating of the heart is a synchronized contraction of muscle cells (myocytes) that are triggered by a periodic sequence of electrical waves (action potentials) originating in the sino-atrial node and propagating over the atria and the ventricles. Cardiac arrhythmias like atrial and ventricular fibrillation (AF,VF) or ventricular tachycardia (VT) are caused by disruptions and instabilities of these electrical excitations, that lead to the emergence of rotating waves (VT) and turbulent wave patterns (AF,VF). Numerous simulation and experimental studies during the last 20 years have addressed these topics. In this review we focus on the nonlinear dynamics of wave propagation in the heart with an emphasis on the theory of pulses, spirals and scroll waves and their instabilities in excitable media and their application to cardiac modeling. After an introduction into electrophysiological models for action potential propagation, the modeling and analysis of spatiotemporal alternans, spiral and scroll meandering, spiral breakup and scroll wave instabilities like negative line tension and sproing are reviewed in depth and discussed with emphasis on their impact in cardiac arrhythmias.Peer ReviewedPreprin

Out of equilibrium dynamics of coherent non-abelian gauge fields

We study out-of-equilibrium dynamics of intense non-abelian gauge fields. Generalizing the well-known Nielsen-Olesen instabilities for constant initial color-magnetic fields, we investigate the impact of temporal modulations and fluctuations in the initial conditions. This leads to a remarkable coexistence of the original Nielsen-Olesen instability and the subdominant phenomenon of parametric resonance. Taking into account that the fields may be correlated only over a limited transverse size, we model characteristic aspects of the dynamics of color flux tubes relevant in the context of heavy-ion collisions.Comment: 12 pages, 10 figures; PRD version, minor change

Can a microscopic stochastic model explain the emergence of pain cycles in patients?

A stochastic model is here introduced to investigate the molecular mechanisms which trigger the perception of pain. The action of analgesic drug compounds is discussed in a dynamical context, where the competition with inactive species is explicitly accounted for. Finite size effects inevitably perturb the mean-field dynamics: Oscillations in the amount of bound receptors spontaneously manifest, driven by the noise which is intrinsic to the system under scrutiny. These effects are investigated both numerically, via stochastic simulations and analytically, through a large-size expansion. The claim that our findings could provide a consistent interpretative framework to explain the emergence of cyclic behaviors in response to analgesic treatments, is substantiated.Comment: J. Stat. Mech. (Proceedings UPON2008