Moment-based formulation of Navier–Maxwell slip boundary conditions for lattice Boltzmann simulations of rarefied flows in microchannels

We present an implementation of first-order Navier–Maxwell slip boundary conditions for simulating near-continuum rarefied flows in microchannels with the lattice Boltzmann method. Rather than imposing boundary conditions directly on the particle velocity distribution functions, following the existing discrete analogs of the specular and diffuse reflection conditions from continuous kinetic theory, we use a moment-based method to impose the Navier–Maxwell slip boundary conditions that relate the velocity and the strain rate at the boundary. We use these conditions to solve for the unknown distribution functions that propagate into the\ud domain across the boundary. We achieve second-order accuracy by reformulating these conditions for the second set of distribution functions that arise in the derivation of the lattice Boltzmann method by an integration along characteristics. The results are in excellent agreement with asymptotic solutions of the compressible Navier-Stokes equations for microchannel flows in the slip regime. Our moment formalism is also valuable for analysing the existing boundary conditions, and explains the origin of numerical slip in the bounce-back and other common boundary conditions that impose explicit conditions on the higher moments instead of on the local tangential velocity

A random projection method for sharp phase boundaries in lattice Boltzmann simulations

Existing lattice Boltzmann models that have been designed to recover a macroscopic description of immiscible liquids are only able to make predictions that are quantitatively correct when the interface that exists between the fluids is smeared over several nodal points. Attempts to minimise the thickness of this interface generally leads to a phenomenon known as lattice pinning, the precise cause of which is not well understood. This spurious behaviour is remarkably similar to that associated with the numerical simulation of hyperbolic partial differential equations coupled with a stiff source term. Inspired by the seminal work in this field, we derive a lattice Boltzmann implementation of a model equation used to investigate such peculiarities. This implementation is extended to different spacial discretisations in one and two dimensions. We shown that the inclusion of a quasi-random threshold dramatically delays the onset of pinning and facetting

Phase mixing vs. nonlinear advection in drift-kinetic plasma turbulence

A scaling theory of long-wavelength electrostatic turbulence in a magnetised, weakly collisional plasma (e.g., ITG turbulence) is proposed, with account taken both of the nonlinear advection of the perturbed particle distribution by fluctuating ExB flows and of its phase mixing, which is caused by the streaming of the particles along the mean magnetic field and, in a linear problem, would lead to Landau damping. It is found that it is possible to construct a consistent theory in which very little free energy leaks into high velocity moments of the distribution function, rendering the turbulent cascade in the energetically relevant part of the wave-number space essentially fluid-like. The velocity-space spectra of free energy expressed in terms of Hermite-moment orders are steep power laws and so the free-energy content of the phase space does not diverge at infinitesimal collisionality (while it does for a linear problem); collisional heating due to long-wavelength perturbations vanishes in this limit (also in contrast with the linear problem, in which it occurs at the finite rate equal to the Landau-damping rate). The ability of the free energy to stay in the low velocity moments of the distribution function is facilitated by the "anti-phase-mixing" effect, whose presence in the nonlinear system is due to the stochastic version of the plasma echo (the advecting velocity couples the phase-mixing and anti-phase-mixing perturbations). The partitioning of the wave-number space between the (energetically dominant) region where this is the case and the region where linear phase mixing wins its competition with nonlinear advection is governed by the "critical balance" between linear and nonlinear timescales (which for high Hermite moments splits into two thresholds, one demarcating the wave-number region where phase mixing predominates, the other where plasma echo does).Comment: 45 pages (single-column), 3 figures, replaced with version published in JP

Three-dimensional central-moments-based lattice Boltzmann method with external forcing: A consistent, concise and universal formulation

The cascaded or central-moments-based lattice Boltzmann method (CM-LBM) is a robust alternative to the more conventional BGK-LBM for the simulation of high-Reynolds number flows. Unfortunately, its original formulation makes its extension to a broader range of physics quite difficult. To tackle this issue, a recent work [A. De Rosis, Phys. Rev. E 95, 013310 (2017)] proposed a more generic way to derive concise and efficient three-dimensional CM-LBMs. Knowing the original model also relies on central moments that are derived in an adhoc manner, i.e., by mimicking those of the Maxwell-Boltzmann distribution to ensure their Galilean invariance a posteriori, a very recent effort [A. De Rosis and K. H. Luo, Phys. Rev. E 99, 013301 (2019)] was proposed to further generalize their derivation. The latter has shown that one could derive Galilean invariant CMs in a systematic and a priori manner by taking into account high-order Hermite polynomials in the derivation of the discrete equilibrium state. Combining these two approaches, a compact and mathematically sound formulation of the CM-LBM with external forcing is proposed. More specifically, the proposed formalism fully takes advantage of the D3Q27 discretization by relying on the corresponding set of 27 Hermite polynomials (up to the sixth order) for the derivation of both the discrete equilibrium state and the forcing term. The present methodology is more consistent than previous approaches, as it properly explains how to derive Galilean invariant CMs of the forcing term in an a priori manner. Furthermore, while keeping the numerical properties of the original CM-LBM, the present work leads to a compact and simple algorithm, representing a universal methodology based on CMs and external forcing within the lattice Boltzmann framework.Comment: Published in Phys. Fluids as Editor's Pic

A volume-preserving sharpening approach for the propagation of sharp phase boundaries in multiphase lattice Boltzmann simulations

Lattice Boltzmann models that recover a macroscopic description of multiphase flow of immiscible liquids typically represent the boundaries between phases using a scalar function, the phase field, that varies smoothly over several grid points. Attempts to tune the model parameters to minimise the thicknesses of these interfaces typically lead to the interfaces becoming fixed to the underlying grid instead of advecting with the fluid velocity. This phenomenon, known as lattice pinning, is strikingly similar to that associated with the numerical simulation of conservation laws coupled to stiff algebraic source terms. We present a lattice Boltzmann formulation of the model problem proposed by LeVeque and Yee [J. Comput. Phys. 86, 187] to study the latter phenomenon in the context of computational combustion, and offer a volume-conserving extension in multiple space dimensions. Inspired by the random projection method of Bao and Jin [J. Comput. Phys. 163, 216] we further generalise this formulation by introducing a uniformly distributed quasi-random variable into the term responsible for the sharpening of phase boundaries. This method is mass conserving and the statistical average of this method is shown to significantly delay the onset of pinning

Brewing of filter coffee

We report progress on mathematical modelling of coffee grounds in a drip filter coffee machine. The report focuses on the evolution of the shape of the bed of coffee grounds during extraction with some work also carried out on the chemistry of extraction. This work was sponsored by Philips who are interested in understanding an observed correlation between the final shape of the coffee grounds and the quality of the coffee. We used experimental data gathered by Philips and ourselves to identify regimes in the coffee brewing process and relevant regions of parameter space. Our work makes it clear that a number of separate processes define the shape of the coffee bed depending on the values of the parameters involved e.g. the size of the grains and the speed of fluid flow during extraction. We began work on constructing mathematical models of the redistribution of the coffee grounds specialised to each region and on a model of extraction. A variety of analytic and numerical tools were used. Furthermore our research has progressed far enough to allow us to begin to exploit connections between this problem and other areas of science, in particular the areas of sedimentology and geomorphology, where the processes we have observed in coffee brewing have been studied

A critical perspective on second-order empathy in understanding psychopathology: phenomenology and ethics

The centenary of Karl Jaspers’ General Psychopathology was recognised in 2013 with the publication of a volume of essays dedicated to his work (edited by Stanghellini and Fuchs). Leading phenomenological-psychopathologists and philosophers of psychiatry examined Jaspers notion of empathic understanding and his declaration that certain schizophrenic phenomena are ‘un-understandable’. The consensus reached by the authors was that Jaspers operated with a narrow conception of phenomenology and empathy and that schizophrenic phenomena can be understood through what they variously called second-order and radical empathy. This article offers a critical examination of the second-order empathic stance along phenomenological and ethical lines. It asks: (1) Is second-order empathy (phenomenologically) possible? (2) Is the second-order empathic stance an ethically acceptable attitude towards persons diagnosed with schizophrenia? I argue that second-order empathy is an incoherent method that cannot be realised. Further, the attitude promoted by this method is ethically problematic insofar as the emphasis placed on radical otherness disinvests persons diagnosed with schizophrenia from a fair chance to participate in the public construction of their identity and, hence, to redress traditional symbolic injustices

Nonlinear dynamics of magnetohydrodynamic flows of heavy fluid over an arbitrary surface in shallow water approximation

The magnetohydrodynamic equations system for heavy fluid over an arbitrary surface in shallow water approximation is studied in the present paper. It is shown that simple wave solutions exist only for underlying surfaces that are slopes of constant inclination. All self-similar discontinuous and continuous solutions are found. The exact explicit solutions of initial discontinuity decay problem over a flat plane and a slope are found. It is shown that the initial discontinuity decay solution is represented by one of five possible wave configurations. For each configuration the necessary and sufficient conditions for its realization are found. The change of dependent and independent variables transforming the initial equations over a slope to those over a flat plane is found.Comment: 43 pages, submitted to Theoretical and Computational Fluid Dynamic