Abstract

ABSTRACTMany researchers have proposed special treatments for outlet boundary conditions owing to lack of information at the outlet. Among them, the simplest method requires a large enough computational domain to prevent or reduce numerical errors at the boundaries. However, an efficient method generally requires special treatment to overcome the problems raised by the outlet boundary condition used. For example, mass flux is not conserved and the fluid field is not divergence-free at the outlet boundary. Overcoming these problems requires additional computational cost. In this paper, we present a simple and efficient outflow boundary condition for the incompressible Navier–Stokes equations, aiming to reduce the computational domain for simulating flow inside a long channel in the streamwise direction. The proposed outflow boundary condition is based on the transparent equation, where a weak formulation is used. The pressure boundary condition is derived by using the Navier–Stokes equations and the outlet flow boundary condition. In the numerical algorithm, a staggered marker-and-cell grid is used and temporal discretization is based on a projection method. The intermediate velocity boundary condition is consistently adopted to handle the velocity–pressure coupling. Characteristic numerical experiments are presented to demonstrate the robustness and accuracy of the proposed numerical scheme. Furthermore, the agreement of computational results from small and large domains suggests that our proposed outflow boundary condition can significantly reduce computational domain sizes.

Highlights

  • Microfluidics, which is mostly in a low Reynolds number flow regime according to its velocity and length scales, offers fundamentally new capabilities for the control of concentrations of molecules in space and time (Whitesides, 2006)

  • We presented a simple and efficient outflow boundary condition for an incompressible Navier–Stokes equation

  • The proposed outflow boundary condition is based on the transparent equation, where a weak formulation is used

Read more

Summary

Introduction

Microfluidics, which is mostly in a low Reynolds number flow regime according to its velocity and length scales, offers fundamentally new capabilities for the control of concentrations of molecules in space and time (Whitesides, 2006). Once mass flux has been kept as a constant, the convective boundary condition becomes a fixed Dirichlet boundary condition Another effective method derived from a wave equation is the nonreflecting boundary method (Hedstrom, 1979; Jin & Braza, 1993; Johansson, 1993; Rudy & Strikwerda, 1980; Sani & Gresho, 1994; Thompson, 1987), which performs well to minimize the spurious artifacts at the outlet. This kind of model is suitable for wake and jet flow with moderate and high Reynolds numbers.

Governing equations and boundary conditions
Velocity boundary conditions
Re uyy and vt uvx
Pressure boundary conditions
Numerical methods
Boundary conditions for velocity and pressure gradient
Boundary condition for intermediate velocity
Mass flux correction algorithm
Parabolic flow
Convergence test
Backward-facing step flow
Comparison with Neumann boundary condition
Efficiency of the proposed method
Three-dimensional flow
Block on the outlet
Findings
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call