An implicit algorithm based on iterative modified approximate factorization method coupling with characteristic boundary conditions for solving subsonic viscous flows

An algorithm composed of an iterative modified approximate factorization(MAF(k)) method with Navier-Stokes characteristic boundary conditions(NSCBC) is proposed for solving subsonic viscous flows.A transformation on the matrix equation in MAF(k) is made in order to impose the implicit boundary conditions properly.To be in consistent with the implicit solver for the interior domain,an implicit scheme for NSCBC is formulated.The performance of the developed algorithm is investigated using spatially evolving zero pressure gradient boundary layer over a flat plate and a wall jet mixing with a cross flow over a flat plate with a square hole as the test cases.The numerical results are compared to the existing experimental datasets and a number of general correlations,together with other available numerical solutions,which demonstrate that the developed algorithm possesses promising capacity for simulating the subsonic viscous flows with large CFL number.

Revised Three-Dimensional Navier-Stokes Characteristic Boundary Conditions for Intense Reactive Turbulence

The three-dimensional Navier-Stokes characteristic boundary conditions(3D-NSCBC), although physically reasonable and popular in many applications, may encounter the instability problem in simulating complex flows, especially for large Reynolds number reactive turbulence where locally the strong reversed flow appears at the outflow boundary surfaces. In the present work, a revised 3D-NSCBC strategy is proposed based on the kinematic relation in different moving coordinate systems. Following this strategy, a systematic formulation is presented for the outflow surface with local reversed flow and can be easily extended to the coupled edge and corner boundaries. Direct numerical simulation(DNS) tests of flow with different turbulence intensities are carried out. Compared with the conventional 3D-NSCBC, the newly proposed method exhibits satisfactory performance to confine numerical instability in the strong reversed flow region. The results confirm the robustness and effectiveness of this newly proposed algorithm.

Numerical simulation of hydrodynamic characteristics during the diversion closure in a horizontal tunnel

Based on water inrush accident of 1841 working face of Desheng Coal Mine in Wu'an, Hebei province, China, an evaluation model of hydrodynamic characteristics of the project is set up and simulated using Matlab. It is assumed that the pipe flow would transform into seepage flow when the aggregates are plugged into the water inrush channel and the seepage flow would disappear along with grouting process. The simulation results show that the flow velocity will increase with an increase in height of aggregates accumulation body during the aggregates filling process; the maximum seepage velocity occurs on the top of plugging zone; and the water flow decreases with increasing plugging height of water inrush channel. Finally, the field construction results show that the water inrush channel can be plugged effectively by the compacted body prepared with aggregate and cement slurry.

A Characteristic Boundary Condition for Multispeed Lattice Boltzmann Methods

We present the development of a non-reflecting boundary condition,based on the Local One-Dimensional Inviscid(LODI)approach,for Lattice Boltzmann Models working with multi-speed stencils.We test and evaluate the LODI implementation with numerical benchmarks,showing significant accuracy gains with respect to the results produced by a simple zerogradient condition.We also implement a simplified approach,which allows handling the unknown distribution functions spanning several layers of nodes in a unified way,still preserving a comparable level of accuracy with respect to the standard formulation.

A COMPUTATIONAL METHOD FOR THE NAVIER-STOKES EQUATIONS AT ALL SPEEDS

A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using A USMPW scheme in conjunction with the third-order MUSCL scheme with Van Leer limiter. The present method was applied to solve the multidimensional compressible Navier-Stokes equations in curvilinear coordinates. Characteristic boundary conditions based on the eigensystem of the preconditioned equations were employed. In order to examine the performance of present method, driven-cavity flow at various Reynolds numbers and viscous flow through a convergent-divergent nozzle at supersonic were selected to rest this method. The computed results were compared with the experimental data or the other numerical results available in literature and good agreements between them are obtained. The results show that the present method is accurate, self-adaptive and stable for a wide range of flow conditions from low speed to supersonic flows.