On the application of slip boundary condition on curved boundaries

Hydrodynamic simulations of sloshing phenomena often involve the application of slip boundary condition at the wetted surfaces. If these surfaces are curved, the ambiguous nature of the normal vector in the discretized problem can interfere with the application of such a boundary condition. Even the use of consistent normal vectors, preferred from the point of view of conservation, does not assure good approximation of the continuum slip condition in the discrete problem, and non‐physical recirculating flow fields may be observed. As a remedy, we consider the Navier slip condition, and more successfully, the so‐called BC‐free boundary condition. Copyright © 2004 John Wiley & Sons, Ltd.     

A consistent and fast weakly compressible smoothed particle hydrodynamics with a new wall boundary condition

A modified weakly compressible smoothed particle hydrodynamics (WCSPH) is presented, which utilizes consistent discretization schemes for spatial derivatives in the flow equations. Here, each SPH particle is considered as a computational point that represents a specific part of the fluid. To overcome non‐physical oscillations that usually arise in standard WCSPH, we modified the mass conservation equation by using a numerical filter. This modification is based on the difference between two discretization schemes used for the term . Furthermore, a new implementation of wall boundary condition in SPH is introduced. This condition is imposed on the pressure of wall boundary particles to ensure that the acceleration of each boundary particle in normal direction to the wall is zero. Thus, no penetration through walls will occur. To examine the performance of the modified method, we solved a series of two‐dimensional incompressible internal flow benchmark problems. By comparing the result with analytical solutions and the results of the standard WCSPH, we show that the use of consistent schemes in conjunction with the proposed numerical filter improves both accuracy and speed of the numerical method. Copyright © 2011 John Wiley & Sons, Ltd.     

NUMERICAL SOLUTION OF THE SINGULARLY PERTURBED PROBLEM WITH NONLOCAL BOUNDARY CONDITION

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｉｓｐａｐｅｒｉｓｃｏｎｃｅｒｎｅｄｗｉｔｈｔｈｅｎｕｍｅｒｉｃａｌｓｏｌｕｔｉｏｎ ,ｂｙｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄ ,ｏｆｔｈｅｆｏｌｌｏｗｉｎｇｓｉｎｇｕｌａｒｌｙｐｅｒｔｕｒｂｅｄｎｏｎｌｏｃａｌｂｏｕｎｄａｒｙｖａｌｕｅｐｒｏｂｌｅｍ :Ｌｕ≡-εｕ″ +ａ(ｘ)ｕ=ｆ(ｘ)　　 ( 0 <ｘ<ｌ) ,( 1 )Ｌ0 ｕ≡-εｕ′( 0 ) +γｕ( 0 ) =μ0 , ( 2 )Ｌ1ｕ≡ｕ(ｌ) -δｕ(ｄ) =μｌ　　 ( 0 <ｄ<ｌ) ,( 3 )ｗｈｅｒｅεｉｓａｓｍａｌｌｐｏｓｉｔｉｖｅｐａｒａｍｅｔｅｒ,…     

Flow of a fluid-solid mixture: Normal stress differences and slip boundary condition

In this paper, using mixture theory we study the flow of a dense suspension, composed of solid particles and a fluid; the emphasis is on the influence of the slip boundary condition and the effect of normal stress differences. Very little work has been done considering both the slip at the walls and the normal stress effects in the frame of a two-component flow. In this paper, the stress tensor for the solid component is modeled as a nonlinear fluid which not only includes the viscous effects but also the normal stress effects; the fluid constituent is modeled as a viscous fluid. We look at the flow between two flat plates.     

Stokes flow of micro-polar fluids by peristaltic pumping through tube with slip boundary condition

This paper studies the Stokes flow of micro-polar fluids by peristaltic pumping through the cylindrical tube under the effect of the slip boundary condition. The motion of the wall is governed by the sinusoidal wave equation. The analytical and numerical solutions for the axial velocity, the micro-polar vector, the stream function, the pressure gradient, the friction force, and the mechanical efficiency are obtained by using the lubrication theory under the low Reynolds number and long wavelength approximations. The impacts of the emerging parameters, such as the coupling number, the micro-polar parameter, the slip parameter on pumping characteristics, the friction force, the velocity profile, the mechanical efficiency, and the trapping phenomenon are depicted graphically. The numerical results infer that large pressure is required for peristaltic pumping when the coupling number is large, while opposite behaviors are found for the micro-polar parameter and the slip parameter. The size of the trapped bolus reduces with the increase in the coupling number and the micro-polar parameter, whereas it blows up with the increase in the slip parameter.     

Electroosmotic flow of Eyring fluid in slit microchannel with slip boundary condition

In consideration of the electroosmotic flow in a slit microchannel, the con-stitutive relationship of the Eyring fluid model is utilized. Navier＇s slip condition is used as the boundary condition. The governing equations are solved analytically, yielding the velocity distribution. The approximate expressions of the velocity distribution are also given and discussed. Furthermore, the effects of the dimensionless parameters, the electrokinetic parameter, and the slip length on the flow are studied numerically, and appropriate conclusions are drawn.     

Time‐dependent flow across a step: the slip with friction boundary condition

The paper studies numerically the slip with friction boundary condition in the time‐dependent incompressible Navier–Stokes equations. Numerical tests on two‐ and three‐dimensional channel flows across a step using this boundary condition on the bottom wall are performed. The influence of the friction parameter on the flow field is studied and the results are explained according to the physics of the flow. Due to the stretching and tilting of vortices, the three‐dimensional results differ in many respects from the two‐dimensional ones. Copyright © 2005 John Wiley & Sons, Ltd.     

Meshless Galerkin analysis of Stokes slip flow with boundary integral equations

This paper presents a novel meshless Galerkin scheme for modeling incompressible slip Stokes flows in 2D. The boundary value problem is reformulated as boundary integral equations of the first kind which is then converted into an equivalent variational problem with constraint. We introduce a Lagrangian multiplier to incorporate the constraint and apply the moving least‐squares approximations to generate trial and test functions. In this boundary‐type meshless method, boundary conditions can be implemented exactly and system matrices are symmetric. Unlike the domain‐type method, this Galerkin scheme requires only a nodal structure on the bounding surface of a body for approximation of boundary unknowns. The convergence and abstract error estimates of this new approach are given. Numerical examples are also presented to show the efficiency of the method. Copyright © 2009 John Wiley & Sons, Ltd.     

Implementation of slip boundary conditions in the finite volume method: new techniques

Two different techniques for the implementation of the linear and nonlinear slip boundary conditions into a finite volume method based numerical code are presented. For the linear Navier slip boundary condition, an implicit implementation in the system of equations is carried out for which there is no need for any relaxation, especially when handling high slip coefficients. For three different nonlinear slip boundary conditions, two different methods are devised, one based on solving a transcendental equation for the boundary and the other on the linearization of the slip law. For assessment purposes, comparison is made between these new methods and the usual iterative process. With these new methods, the convergence difficulties, typical of the iterative procedure, are eliminated, and for some of the test cases, the convergence rate even increased with the slip velocity. The details of these implementations are given first for a simple geometry using orthogonal meshes and Cartesian coordinates followed by their generalization to non‐Cartesian coordinates and nonorthogonal meshes. The developed code was tested in the benchmark slip‐stick and 4:1 contraction flows, evidencing the robustness of the proposed procedures. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical study of particle deposition in a turbulent channel flow with transverse roughness elements on one wall

A numerical study is presented for the effect of wall roughness on the deposition of solid spherical particles in a fully developed turbulent channel flow based on large eddy simulation combined with a Lagrangian particle-tracking scheme. The interest is focused on particles with response times in wall units in the range of 2.5 ≤ τ_p⁺ ≤ 600 depositing onto a vertical rough surface consisting of two-dimensional transverse square bars separated by a rectangular cavity. Predictions of particle deposition rates are obtained for several values of the cavity width to roughness element height ratio and particle response time. It is shown that the accumulation of particles in the near wall region and their preferential concentration in flow areas of low streamwise fluid velocity that occur in turbulent flows at flat channels are significantly affected by the roughness elements. Particle deposition onto the rough wall is considerably increased, exhibiting a subtle dependence on the particle inertia and the spacing between the bars. The observed augmentation of deposition coefficient can be attributed to the flow modifications induced by the roughness elements and to the inertial impaction of particles onto the frontal deposition area of the protruding square bars.     

Numerical study of bubble and particle motion in a turbulent boundary layer using proper orthogonal decomposition

We have studied the motion of bubbles and particles in the near-wall region of a turbulent boundary layer, to investigate the influence of the unsteady turbulent structure. The velocity field was computed using Proper Orthogonal Decomposition (POD), and the trajectories of bubbles and particle have been computed by integrating their equation of motion. We have used this to investigate the roles, and the relative importance, of the different forces acting on bubbles and particles, We find that the unsteady turbulent structure plays an important role in the preferential accumulation of bubbles and particles. The accumulation of bubbles depends on a rather complicated interaction between the pressure gradient and the lift force; neither is sufficient, acting on its own, to explain the strong accumulation observed when they act together.     

DISTRIBUTION OF RANDOM ELEMENTS SUBJECTED TO A FLEXIBLE BOUNDARY CONDITION

The probability distribution function of n random elements subjected to the flexible boundary condition is derived. The probability density is a descending curve and converges to a delta function as n tends to infinity. The distribution of the minimum value is discussed in context of ordered statistics.     

Flow of fluids with pressure dependent viscosities in an orthogonal rheometer subject to slip boundary conditions

We consider slow steady flows of Navier–Stokes-like fluids with pressure dependent viscosities between rotating infinite parallel plates with Navier slip boundary conditions. We derive exact solutions which correspond to flows in orthogonal and torsional rheometers, and investigate the effect of the slip coefficient and the material parameters on the solutions. We find that even when inertial effects are ignored, vorticity boundary layers develop at the upper boundary due to the pressure dependence of the viscosity. These boundary layers diminish and eventually disappear with increased slippage.     

A high‐order scheme for the incompressible Navier–Stokes equations with open boundary condition

We present a numerical scheme to solve the incompressible Navier–Stokes equations with open boundary condition. After replacing the incompressibility constraint by the pressure Poisson equation, the key is how to give an appropriate boundary condition for the pressure Poisson equation. We propose a new boundary condition for the pressure on the open boundary. Some numerical experiments are presented to verify the accuracy and stability of scheme. Copyright © 2013 John Wiley & Sons, Ltd.     

Creeping motion of a slip spherical particle in a circular cylindrical pore

A combined analytical–numerical study for the creeping flow caused by a spherical fluid or solid particle with a slip-flow surface translating in a viscous fluid along the centerline of a circular cylindrical pore is presented. To solve the axisymmetric Stokes equations for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the pore wall by the Fourier transforms and then on the particle surface by a collocation technique. Numerical results for the hydrodynamic drag force acting on the particle are obtained with good convergence for various values of the relative viscosity or slip coefficient of the particle, the slip parameter of the pore wall, and the ratio of radii of the particle and pore. For the motion of a fluid sphere along the axis of a cylindrical pore, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force for all cases is a monotonic increasing function of the ratio of particle-to-pore radii, and approaches infinity in the limit. Except for the case that the cylindrical pore is hardly slip and the value of the ratio of particle-to-pore radii is close to unity, the drag force exerted on the particle increases monotonically with an increase in its relative viscosity or with a decrease in its slip coefficient for a constant ratio of radii. In a comparison for the pore shape effect on the axial translation of a slip sphere, it is found that the particle in a circular cylindrical pore in general acquires a lower hydrodynamic drag than in a spherical cavity, but this trend can be reversed for the case of highly slippery particles and pore walls.     

Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

A projection scheme for Navier-Stokes with variable viscosity and natural boundary condition

Combiningthe Navier-Stokes systems with Neumann (or natural) boundary condition to characterize a fluid flow is frequent. The popular projection (or pressure correction) methods inspired by Chorin and Temam are not well adapted to such boundary condition, which translate in loss of accuracy. If some alternative projection methods have been proposed to reduce the accuracy loss due to the Neumann condition in case of Newtonian fluids, little has been proposed for generalized Newtonian fluids. In this work, we propose two methods derived from the incremental pressure correction projection that can be used for fluids with inhomogeneous or variable viscosity with natural boundary condition. Both time and space accuracy of the methods will be illustrated using a manufactured solution.     

Improvement of boundary conditions for non-planar boundaries represented by polygons with an initial particle arrangement technique

The boundary conditions represented by polygons in moving particle semi-implicit (MPS) method (Koshizuka and Oka, Nuclear Science and Engineering, 1996 Koshizuka, S., and Y. Oka. 1996. “Moving-Particle Semi-Implicit Method for Fragmentation of Incompressible Fluid.” Nuclear Science and Engineering 123: 421–434.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) have been widely used in the industry simulations since it can simply simulate complex geometry with high efficiency. However, the inaccurate particle number density near non-planar wall boundaries dramatically affects the accuracy of simulations. In this paper, we propose an initial boundary particle arrangement technique coupled with the wall weight function method (Zhang et al. Transaction of JSCES, 2015 Zhang, T.G., S. Koshizuka, K. Shibata, K. Murotani, and E. Ishii. 2015. “Improved Wall Weight Function With Polygon Boundary in Moving Particle Semi-Implicit Method.” Transaction of JSCES. No. 20150012. [Google Scholar]) to improve the particle number density near slopes and curved surfaces with boundary conditions represented by polygons in three dimensions. Two uniform grids are utilized in the proposed technique. The grid points in the first uniform grid are used to construct boundary particles, and the second uniform grid stores the same information as in the work by Zhang et al. The wall weight functions of the grid points in the second uniform grid are calculated by newly constructed boundary particles. The wall weight functions of the fluid particles are interpolated from the values stored on the grid points in the second uniform grid. Because boundary particles are located on the polygons, complex geometries can be accurately represented. The proposed method can dramatically improve the particle number density and maintain the high efficiency. The performance of the previously proposed wall weight function (Zhang et al.) with the boundary particle arrangement technique is verified in comparison with the wall weight function without boundary particle arrangement by investigating two example geometries. The simulations of a water tank with a wedge and a complex geometry show the general applicability of the boundary particle arrangement technique to complex geometries and demonstrate its improvement of the wall weight function near the slopes and curved surfaces.