No Arabic abstract
We propose a multi-resolution strategy that is compatible with the lattice Greens function (LGF) technique for solving viscous, incompressible flows on unbounded domains. The LGF method exploits the regularity of a finite-volume scheme on a formally unbounded Cartesian mesh to yield robust and computationally efficient solutions. The original method is spatially adaptive, but challenging to integrate with embedded mesh refinement as the underlying LGF is only defined for a fixed resolution. We present an ansatz for adaptive mesh refinement, where the solutions to the pressure Poisson equation are approximated using the LGF technique on a composite mesh constructed from a series of infinite lattices of differing resolution. To solve the incompressible Navier-Stokes equations, this is further combined with an integrating factor for the viscous terms and an appropriate Runge-Kutta scheme for the resulting differential-algebraic equations. The parallelized algorithm is verified through with numerical simulations of vortex rings, and the collision of vortex rings at high Reynolds number is simulated to demonstrate the reduction in computational cells achievable with both spatial and refinement adaptivity.
In this paper, a multiple-distribution-function lattice Boltzmann method (MDF-LBM) with multiple-relaxation-time model is proposed for incompressible Navier-Stokes equations (NSEs) which are considered as the coupled convection-diffusion equations (CDEs). Through direct Taylor expansion analysis, we show that the Navier-Stokes equations can be recovered correctly from the present MDF-LBM, and additionally, it is also found that the velocity and pressure can be directly computed through the zero and first-order moments of distribution function. Then in the framework of present MDF-LBM, we develop a locally computational scheme for the velocity gradient where the first-order moment of the non-equilibrium distribution is used, this scheme is also extended to calculate the velocity divergence, strain rate tensor, shear stress and vorticity. Finally, we also conduct some simulations to test the MDF-LBM, and find that the numerical results not only agree with some available analytical and numerical solutions, but also have a second-order convergence rate in space.
In this article, a coupled Two-relaxation-time Lattice Boltzmann-Volume penalization (TRT-LBM-VP) method is presented to simulate flows past obstacles. Two relaxation times are used in the collision operator, of which one is related to the fluid viscosity and the other one is related to the numerical stability and accuracy. The volume penalization method is introduced into the TRT-LBM by an external forcing term. In the procedure of the TRT-LBM-VP, the processes of interpolating velocities on the boundaries points and distributing the force density to the Eulerian points are unneeded. Performing the TRT-LBM-VP on a certain point, only the variables of this point are needed. As a consequence, the TRT-LBM-VP can be conducted parallelly. From the comparison between the result of the cylindrical Couette flow solved by the TRT-LBM-VP and that solved by the Single-relaxation-time LBM-VP (SRT-LBM-VP), the accuracy of the TRT-LBM-VP is higher than that of the SRT-LBM-VP. Flows past a single circular cylinder, a pair of cylinders in tandem and side-by-side arrangements, two counter-rotating cylinders and a NACA-0012 airfoil are chosen as numerical experiments to verify the present method further. Good agreements between the present results and those in the previous literatures are achieved.
In this brief report, a thermal lattice-Boltzmann (LB) model is presented for axisymmetric thermal flows in the incompressible limit. The model is based on the double-distribution-function LB method, which has attracted much attention since its emergence for its excellent numerical stability. Compared with the existing axisymmetric thermal LB models, the present model is simpler and retains the inherent features of the standard LB method. Numerical simulations are carried out for the thermally developing laminar flows in circular ducts and the natural convection in an annulus between two coaxial vertical cylinders. The Nusselt number obtained from the simulations agrees well with the analytical solutions and/or the results reported in previous studies.
This paper addresses how two time integration schemes, the Heuns scheme for explicit time integration and the second-order Crank-Nicolson scheme for implicit time integration, can be coupled spatially. This coupling is the prerequisite to perform a coupled Large Eddy Simulation / Reynolds Averaged Navier-Stokes computation in an industrial context, using the implicit time procedure for the boundary layer (RANS) and the explicit time integration procedure in the LES region. The coupling procedure is designed in order to switch from explicit to implicit time integrations as fast as possible, while maintaining stability. After introducing the different schemes, the paper presents the initial coupling procedure adapted from a published reference and shows that it can amplify some numerical waves. An alternative procedure, studied in a coupled time/space framework, is shown to be stable and with spectral properties in agreement with the requirements of industrial applications. The coupling technique is validated with standard test cases, ranging from one-dimensional to three-dimensional flows.
In this paper, a coupling lattice Boltzmann (LB) model for simulating thermal flows on the standard D2Q9 lattice is developed in the framework of the double-distribution-function (DDF) approach in which the viscous heat dissipation and compression work are considered. In the model, a density distribution function is used to simulate the flow field, while a total energy distribution function is employed to simulate the temperature field. The discrete equilibrium density and total energy distribution functions are obtained from the Hermite expansions of the corresponding continuous equilibrium distribution functions. The pressure given by the equation of state of perfect gases is recovered in the macroscopic momentum and energy equations. The coupling between the momentum and energy transports makes the model applicable for general thermal flows such as non-Boussinesq flows, while the existing DDF LB models on standard lattices are usually limited to Boussinesq flows in which the temperature variation is small. Meanwhile, the simple structure and basic advantages of the DDF LB approach are retained. The model is tested by numerical simulations of thermal Couette flow, attenuation-driven acoustic streaming, and natural convection in a square cavity with small and large temperature differences. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.