Quasiequilibrium lattice Boltzmann models with tunable bulk viscosity for enhancing stability

Taking advantage of a closed-form generalized Maxwell distribution function [ P. Asinari and I. V. Karlin Phys. Rev. E 79 036703 (2009)] and splitting the relaxation to the equilibrium in two steps, an entropic quasiequilibrium (EQE) kinetic model is proposed for the simulation of low Mach number flows, which enjoys both the H theorem and a free-tunable parameter for controlling the bulk viscosity in such a way as to enhance numerical stability in the incompressible flow limit. Moreover, the proposed model admits a simplification based on a proper expansion in the low Mach number limit (LQE model). The lattice Boltzmann implementation of both the EQE and LQE is as simple as that of the standard lattice Bhatnagar-Gross-Krook (LBGK) method, and practical details are reported. Extensive numerical testing with the lid driven cavity flow in two dimensions is presented in order to verify the enhancement of the stability region. The proposed models achieve the same accuracy as the LBGK method with much rougher meshes, leading to an effective computational speed-up of almost three times for EQE and of more than four times for the LQE. Three-dimensional extension of EQE and LQE is also discussed

Aerothermal modeling program, phase 2

The main objectives of the Aerothermal Modeling Program, Phase 2 are: to develop an improved numerical scheme for incorporation in a 3-D combustor flow model; to conduct a benchmark quality experiment to study the interaction of a primary jet with a confined swirling crossflow and to assess current and advanced turbulence and scalar transport models; and to conduct experimental evaluation of the air swirler interaction with fuel injectors, assessments of current two-phase models, and verification the improved spray evaporation/dispersion models

Improved numerical methods for turbulent viscous flows aerothermal modeling program, phase 2

The details of a study to develop accurate and efficient numerical schemes to predict complex flows are described. In this program, several discretization schemes were evaluated using simple test cases. This assessment led to the selection of three schemes for an in-depth evaluation based on two-dimensional flows. The scheme with the superior overall performance was incorporated in a computer program for three-dimensional flows. To improve the computational efficiency, the selected discretization scheme was combined with a direct solution approach in which the fluid flow equations are solved simultaneously rather than sequentially

Screening effects in flow through rough channels

A surprising similarity is found between the distribution of hydrodynamic stress on the wall of an irregular channel and the distribution of flux from a purely Laplacian field on the same geometry. This finding is a direct outcome from numerical simulations of the Navier-Stokes equations for flow at low Reynolds numbers in two-dimensional channels with rough walls presenting either deterministic or random self-similar geometries. For high Reynolds numbers, when inertial effects become relevant, the distribution of wall stresses on deterministic and random fractal rough channels becomes substantially dependent on the microscopic details of the walls geometry. In addition, we find that, while the permeability of the random channel follows the usual decrease with Reynolds, our results indicate an unexpected permeability increase for the deterministic case, i.e., ``the rougher the better''. We show that this complex behavior is closely related with the presence and relative intensity of recirculation zones in the reentrant regions of the rough channel.Comment: 4 pages, 5 figure

Autocatalytic plume pinch-off

A localized source of buoyancy flux in a non-reactive fluid medium creates a plume. The flux can be provided by either heat, a compositional difference between the fluid comprising the plume and its surroundings, or a combination of both. For autocatalytic plumes produced by the iodate-arsenous acid reaction, however, buoyancy is produced along the entire reacting interface between the plume and its surroundings. Buoyancy production at the moving interface drives fluid motion, which in turn generates flow that advects the reaction front. As a consequence of this interplay between fluid flow and chemical reaction, autocatalytic plumes exhibit a rich dynamics during their ascent through the reactant medium. One of the more interesting dynamical features is the production of an accelerating vortical plume head that in certain cases pinches-off and detaches from the upwelling conduit. After pinch-off, a new plume head forms in the conduit below, and this can lead to multiple generations of plume heads for a single plume initiation. We investigated the pinch-off process using both experimentation and simulation. Experiments were performed using various concentrations of glycerol, in which it was found that repeated pinch-off occurs exclusively in a specific concentration range. Autocatalytic plume simulations revealed that pinch-off is triggered by the appearance of accelerating flow in the plume conduit.Comment: 10 figures. Accepted for publication in Phys Rev E. See also http://www.physics.utoronto.ca/nonlinear/papers_chemwave.htm

Numerical investigation of oil injection in a Roots blower operated as expander

The adoption of positive displacement machines as expanders in Organic Rankine Cycles (ORCs) is increasingly common, especially in the low to medium power range. At the same time, these devices often serve as compressor in Vapor-Compression Refrigeration Systems. In both cases, the application of Computation Fluid Dynamics (CFD) to optimize such machines has become an integrated tool in the design process. As a consequence, several challenges associated with the numerical simulation have to be taken into account. For example, the modeling of the gap represents a challenge for the stability of the numerical analysis. The dynamic of the process, combined with deformations of the clearances and of the working chamber has to be considered with extra care. To raise the efficiency of the machine, oil is typically injected. Its numerical modeling imply an extra challenge in the simulation of the actual operation of the machine. The present work is mainly focused on the multi-phase nature of the flow, with a particular analysis of the lubricant oil injected. In this work, a two-lobe Roots blower operated as expander has been simulated with the open-source software OpenFOAM-v1812, using the SCORG-V5.2.2. This analysis highlights the areas that are affected the most by the oil presence in order to highlight the sealing effect it provides

Fluctuating hydrodynamic modelling of fluids at the nanoscale

A good representation of mesoscopic fluids is required to combine with molecular simulations at larger length and time scales (De Fabritiis {\it et. al}, Phys. Rev. Lett. 97, 134501 (2006)). However, accurate computational models of the hydrodynamics of nanoscale molecular assemblies are lacking, at least in part because of the stochastic character of the underlying fluctuating hydrodynamic equations. Here we derive a finite volume discretization of the compressible isothermal fluctuating hydrodynamic equations over a regular grid in the Eulerian reference system. We apply it to fluids such as argon at arbitrary densities and water under ambient conditions. To that end, molecular dynamics simulations are used to derive the required fluid properties. The equilibrium state of the model is shown to be thermodynamically consistent and correctly reproduces linear hydrodynamics including relaxation of sound and shear modes. We also consider non-equilibrium states involving diffusion and convection in cavities with no-slip boundary conditions