skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

Abstract

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes–Brinkman–Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes–Brinkman–Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies the effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition atmore » an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.« less

Authors:
 [1];  [2];  [1]
  1. Irstea, Antony Regional Centre, HBAN, 1 rue Pierre-Gilles de Gennes CS 10030, 92761 Antony cedex (France)
  2. CNRS (UMR 7608), Laboratoire FAST, Batiment 502, Campus University, 91405 Orsay (France)
Publication Date:
OSTI Identifier:
22622277
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 335; Other Information: Copyright (c) 2017 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY; APPROXIMATIONS; BOLTZMANN EQUATION; BOUNDARY CONDITIONS; BOUNDARY LAYERS; COMPUTERIZED SIMULATION; CORRECTIONS; DEFECTS; DIRICHLET PROBLEM; EXACT SOLUTIONS; FINITE ELEMENT METHOD; FLOW MODELS; INTERFERENCE; LAPLACIAN; OSCILLATIONS; PERIODICITY; PERMEABILITY; POROUS MATERIALS; RELAXATION TIME; RESOLUTION; VISCOSITY

Citation Formats

Silva, Goncalo, E-mail: goncalo.nuno.silva@gmail.com, Talon, Laurent, E-mail: talon@fast.u-psud.fr, and Ginzburg, Irina, E-mail: irina.ginzburg@irstea.fr. Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes. United States: N. p., 2017. Web. doi:10.1016/J.JCP.2017.01.023.
Silva, Goncalo, E-mail: goncalo.nuno.silva@gmail.com, Talon, Laurent, E-mail: talon@fast.u-psud.fr, & Ginzburg, Irina, E-mail: irina.ginzburg@irstea.fr. Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes. United States. doi:10.1016/J.JCP.2017.01.023.
Silva, Goncalo, E-mail: goncalo.nuno.silva@gmail.com, Talon, Laurent, E-mail: talon@fast.u-psud.fr, and Ginzburg, Irina, E-mail: irina.ginzburg@irstea.fr. Sat . "Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes". United States. doi:10.1016/J.JCP.2017.01.023.
@article{osti_22622277,
title = {Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes},
author = {Silva, Goncalo, E-mail: goncalo.nuno.silva@gmail.com and Talon, Laurent, E-mail: talon@fast.u-psud.fr and Ginzburg, Irina, E-mail: irina.ginzburg@irstea.fr},
abstractNote = {The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes–Brinkman–Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes–Brinkman–Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies the effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.},
doi = {10.1016/J.JCP.2017.01.023},
journal = {Journal of Computational Physics},
number = ,
volume = 335,
place = {United States},
year = {Sat Apr 15 00:00:00 EDT 2017},
month = {Sat Apr 15 00:00:00 EDT 2017}
}
  • Experimental results for natural convection in a horizontal porous cavity of aspect ratio A = 5 and heated from below are reported. A wide range of governing parameters are covered by careful selection of bead size, solid material, and fluid. These results fully support the effects of fluid-flow parameters (Rayleigh and Prandtl numbers), porous matrix-structure parameters (Darcy and Forchheimer numbers), and the conductivity ratio as predicted by the formulation based on the Darcy-Brinkman-Forchheimer (DBF) equations of motion. The DBF flow model, with variable porosity and variable thermal conductivity in the wall regions, predicts reasonably well in comparison with the experimentalmore » data. However, the difference between the predictions and the measurements increases as the ratio of solid-to-fluid thermal conductivity becomes very large. 32 refs.« less
  • A new immersed boundary-lattice Boltzmann method (IB-LBM) is presented for fully resolved simulations of incompressible viscous flows laden with rigid particles. The immersed boundary method (IBM) recently developed by Breugem (2012) [19] is adopted in the present method, development including the retraction technique, the multi-direct forcing method and the direct account of the inertia of the fluid contained within the particles. The present IB-LBM is, however, formulated with further improvement with the implementation of the high-order Runge–Kutta schemes in the coupled fluid–particle interaction. The major challenge to implement high-order Runge–Kutta schemes in the LBM is that the flow information suchmore » as density and velocity cannot be directly obtained at a fractional time step from the LBM since the LBM only provides the flow information at an integer time step. This challenge can be, however, overcome as given in the present IB-LBM by extrapolating the flow field around particles from the known flow field at the previous integer time step. The newly calculated fluid–particle interactions from the previous fractional time steps of the current integer time step are also accounted for in the extrapolation. The IB-LBM with high-order Runge–Kutta schemes developed in this study is validated by several benchmark applications. It is demonstrated, for the first time, that the IB-LBM has the capacity to resolve the translational and rotational motion of particles with the second-order accuracy. The optimal retraction distances for spheres and tubes that help the method achieve the second-order accuracy are found to be around 0.30 and −0.47 times of the lattice spacing, respectively. Simulations of the Stokes flow through a simple cubic lattice of rotational spheres indicate that the lift force produced by the Magnus effect can be very significant in view of the magnitude of the drag force when the practical rotating speed of the spheres is encountered. This finding may lead to more comprehensive studies of the effect of the particle rotation on fluid–solid drag laws. It is also demonstrated that, when the third-order or the fourth-order Runge–Kutta scheme is used, the numerical stability of the present IB-LBM is better than that of all methods in the literature, including the previous IB-LBMs and also the methods with the combination of the IBM and the traditional incompressible Navier–Stokes solver. - Highlights: • The IBM is embedded in the LBM using Runge–Kutta time schemes. • The effectiveness of the present IB-LBM is validated by benchmark applications. • For the first time, the IB-LBM achieves the second-order accuracy. • The numerical stability of the present IB-LBM is better than previous methods.« less
  • Analysis of flow and convective heat transfer in volumetrically heated porous layer has become a separate topic for research in the last twenty five years in view of its importance in various engineering applications, such as heat removal from nuclear fuel debris, heat transfer associated with storage of nuclear waste, exothermic reaction in packed-bed reactors, heat recovery from geothermal systems and particularly in the field of large storage systems of agricultural products. Here, a pressure-velocity solution for natural convection for fluid saturated heat generating porous medium in a square enclosure is analyzed by finite element method. The numerical solutions obtainedmore » for wide range of fluid Rayleigh number, Ra{sub f}, Darcy number, Da, and heat generating number, Q{sub d}. The justification for taking these non-dimensional parameters independently is to establish the effect of individual parameters on flow patterns. It has been observed that peak temperature occurs at the top central part and weaker velocity prevails near the vertical walls of the enclosure due to the heat generation parameter alone. On comparison, the modified Rayleigh number used by the earlier investigators, can not explain explicitly the effect of heat generation parameter on natural convection within an enclosure having differentially heated vertical walls. At higher Darcy number, the peak temperature and peak velocity are comparatively more, resulting in better enhancement of heat transfer rate.« less
  • Rarefied gas flow behavior is usually described by the Boltzmann equation, the Navier-Stokes system being valid when the gas is less rarefied. Slip boundary conditions for the Navier-Stokes equations are derived in a rigorous and systematic way from the boundary condition at the kinetic level (Boltzmann equation). These slip conditions are explicitly written in terms of asymptotic behavior of some linear half-space problems. The validity of this analysis is established in the simple case of the Couette flow, for which is it proved that the right boundary conditions are obtained.
  • We present a systematic method for constructing boundary conditions (numerical and physical) of the required accuracy, for compact (Pade-like) high-order finite-difference schemes for hyperbolic systems. First a proper summation-by-parts formula is found for the approximate derivative. A [open quotes]simultaneous approximation term[close quotes] is then introduced to treat the boundary conditions. This procedure leads to time-stable schemes even in the system case. An explicit construction of the fourth-order compact case is given. Numerical studies are presented to verify the efficacy of the approach. 8 refs., 5 tabs.