Navier Slip Boundary Condition Research Articles

A 3D Non-Stationary Micropolar Fluids Equations with Navier Slip Boundary Conditions

Micropolar fluids are fluids with microstructure and belong to a class of fluids with asymmetric stress tensor that called Polar fluids, and include, as a special case, the well-established Navier–Stokes model. In this work we study a 3D micropolar fluids model with Navier boundary conditions without friction for the velocity field and homogeneous Dirichlet boundary conditions for the angular velocity. Using the Galerkin method, we prove the existence of weak solutions and establish a Prodi–Serrin regularity type result which allow us to obtain global-in-time strong solutions at finite time.

Axisymmetric lattice Boltzmann model with slip boundary conditions for liquid flows in microtube

Lattice Boltzmann (LB) method, which is a mesoscopic numerical method, has been considered as a powerful tool to study microscale gas and liquid flows. Boundary slip phenomenon, which is a significant feature of both microscale gas and liquid flows, has been extensively studied by the LB method in the past decade. However, most of the previous works have focused on the microchannel flows and studies on the microtube flows are very rare. In this paper, we investigate the widely used slip boundary conditions, i.e., combined bounce-back and specular-reflection (BS) scheme, combined Maxwell-diffusion and specular-reflection (MS) scheme, and combined Maxwell-diffusion and bounce-back (MB) scheme, for the axisymmetric LB model with multi-relaxation-time (MRT) in detail. In order to realize the Navier slip boundary condition for liquid flows, we put forward to a reasonable strategy for determining the combination coefficients and the relaxation time. The proposed boundary schemes are validated by some numerical tests including the Hagen–Poiseuille flow, axisymmetric Womersley flow, Poiseuille flow in a circular annulus, and Womersley flow in a circular annulus. Numerical results are consistent with the analytical solutions, which demonstrate that the proposed boundary schemes are suitable for microscale liquid flows in microtube.

Numerical investigation of flow around a structure using Navier-slip boundary conditions

PurposeComplicated motion of vortex is frequently observed in the wake of islands. These kinds of swirling fluid cause the trap of sediments or pollutants, subsequently inducing the dead zone, odor or poor water quality. Therefore, the understanding of flow past a circular cylinder is significant in predicting water quality and positioning the immersed structures. This study aims to investigate the flow properties around a structure using Navier-slip boundary conditions.Design/methodology/approachBoundary conditions are a major factor affecting the flow pattern because the magnitude of flow detachment on a surface can redistribute the tangential stress on the wall. Therefore, the authors performed an analysis of laminar flow passing through a circular structure to investigate the effect of boundary conditions on the flow pattern.FindingsThe authors examined the relationship between the partial-slip boundary conditions and the flow behavior at low Reynolds number past a circular cylinder considering velocity and vorticity distributions behind the cylinder, lift coefficient and Strouhal number. The amplitude of lift coefficient by the partial slip condition had relatively small value compared with that of no-slip condition, as the wall shear stress acting on the cylinder became smaller by the velocity along the cylinder surface. The frequency of the asymmetrical vortex formation with partial slip velocity was increased compared with no-slip case due to the intrinsic inertial effect of Navier-slip condition.Originality/valueThe ability to engineer slip could have dramatic influences on flow, as the viscous dominated motion can lead to large pressure drops and large axial dispersion. By the slip length control, no-slip, partial-slip and free-slip boundary conditions are tunable, and the velocity distributions at the wall, vortex formation and wake pattern including the amplitude of lift coefficient and frequency were significantly affected by slip length parameter.

Liouville theorem of D-solutions to the stationary magnetohydrodynamics system in a slab

In this paper, we study Liouville theorems of D-solutions to the stationary magnetohydrodynamic system in a slab. We will prove trivialness of the velocity and the magnetic field with various boundary conditions. In some boundary conditions, the additional assumption that the horizontal angular component(s) of the velocity or (and) the magnetic field is (are) axially symmetric is needed. More precisely, five types of boundary conditions will be considered: the vertical periodic boundary condition for the velocity and the magnetic field, the Navier-slip boundary condition for the velocity, the perfectly conducting or insulating boundary condition for the magnetic field, the non-slip boundary condition for the velocity, and the perfectly conducting or insulating boundary condition for the magnetic field. One of our innovations is that we do not impose finite Dirichlet integral assumption on the magnetic field compared with previous works.

Initial boundary value problems for two-dimensional isentropic compressible Navier-Stokes equations with rotating effect terms

We consider the initial boundary value problem for two-dimensional isentropic compressible Navier-Stokes equations with rotating effect terms. This problem is related to the motion of the compressible viscous flow surrounding a rotating obstacle. We show the global existence and uniqueness of strong solutions with initial data close to some approximate solution, provided that the friction coefficient α attached to the Navier-slip boundary condition is small. The large time behavior of the solutions is also studied.

Hydrodynamic instability and wave formation of a viscous film flowing down a slippery inclined substrate: Effect of odd-viscosity

Hydrodynamic instability and formation of waves for thin viscous film flowing down a slippery inclined substrate with broken time-reversal-symmetry have been discussed in this present study. The effect of slip on the substrate is modeled using the Navier slip boundary condition. We have derived a nonlinear evolution model in the framework of the long-wave expansion technique. The one equation model can track the free surface evolution and involve the viscosity, gravity, surface tension, and the effect of the slip of the wall. The odd part of the Cauchy stress tensor with an odd-viscosity coefficient significantly modifies the characteristic of the film flow. We performed the linear stability analysis with the Orr Sommerfeld technique and classical temporal analysis from the model. Comparing the results shows a satisfactory agreement between Orr Sommerfeld and the model when the velocity scale is chosen as twice the free surface velocity. Analyzing the traveling wave solutions and studying different bifurcation analyzes and phase diagrams, we focused on and discussed two significant wave families, namely γ1 (low-wave-speed) and γ2 (high-wave-speed). The spatio-temporal study of the model has been performed numerically for different odd-viscosity and slip parameters. Results show that the odd-viscosity plays a stabilizing role while the slip on the substrate gives a destabilizing effect within the parameter range of our consideration.

Correlation in thermal fluctuations induced by phase-locked hydrodynamic modes.

Thermal fluctuations constitute a fundamental equilibrium phenomenon whose spatial and temporal correlations are governed by the relevant scales of molecular collisions. From the continuum point of view, thermal fluctuations in a fluid can be regarded as comprising a multitude of hydrodynamic modes (HMs) with random phases, each one having one degree of freedom. We show that in a two-dimensional fluid channel with the Navier slip boundary condition, in which the HMs are represented by periodic arrays of vortex and antivortex pairs, periodic modulation of the slip boundary condition can selectively suppress noncommensurate HMs while phase lock the remaining eigenmodes. As a result, thermal fluctuations would exhibit mesoscopic-scale spatial correlations, manifest as a spatially varying diffusion constant when evaluated from the fluctuation-dissipation theorem. Good agreement is shown with the molecular dynamics results. Such manifestation of equilibrium collective motion implies that instead of just being an alternative mathematical basis for expressing thermal fluctuations, in mesoscopic systems the HMs may be manipulated to have physical consequences very different from those expected in bulk fluid.

A cut finite element method for non-Newtonian free surface flows in 2D - application to glacier modelling

In ice sheet and glacier modelling, the Finite Element Method is rapidly gaining popularity. However, constructing and updating meshes for ice sheets and glaciers is a non-trivial and computationally demanding task due to their thin, irregular, and time dependent geometry. In this paper we introduce a novel approach to ice dynamics computations based on the unfitted Finite Element Method CutFEM, which lets the domain boundary cut through elements. By employing CutFEM, complex meshing and remeshing is avoided as the glacier can be immersed in a simple background mesh without loss of accuracy. The ice is modelled as a non-Newtonian, shear-thinning fluid obeying the p-Stokes (full Stokes) equations with the ice atmosphere interface as a moving free surface. A Navier slip boundary condition applies at the glacier base allowing both bedrock and subglacial lakes to be represented. Within the CutFEM framework we develop a strategy for handling non-linear viscosities and thin domains and show how glacier deformation can be modelled using a level set function. In numerical experiments we show that the expected order of accuracy is achieved and that the method is robust with respect to penalty parameters. As an application we compute the velocity field of the Swiss mountain glacier Haut Glacier d'Arolla in 2D with and without an underlying subglacial lake, and simulate the glacier deformation from year 1930 to 1932, with and without surface accumulation and basal melt.

Global Existence and Large Time Behavior of Strong Solutions for 3D Nonhomogeneous Heat-Conducting Magnetohydrodynamic Equations

We are concerned with an initial boundary value problem of nonhomogeneous heat-conducting magnetohydrodynamic equations in a bounded simply connected smooth domain $$\Omega \subseteq {\mathbb {R}}^3$$ , with Navier-slip boundary conditions for the velocity and magnetic field and Neumann boundary condition for the temperature. We prove the global existence of a unique strong solution provided that $$\left( \Vert \sqrt{\rho _0}{\mathbf {u}}_0\Vert _{L^2}^2 +\Vert {\mathbf {b}}_0\Vert _{L^2}^2\right) $$ $$\left( \Vert {{\,\mathrm{curl}\,}}{\mathbf {u}}_0\Vert _{L^2}^2 +\Vert {{\,\mathrm{curl}\,}}{\mathbf {b}}_0\Vert _{L^2}^2\right) $$ is suitably small. Moreover, we also obtain large time decay rates of the solution.

Thermocapillary instability on a film falling down a non-uniformly heated slippery incline

A gravity-driven, thin, incompressible liquid film flow on a non-uniformly heated, slippery inclined plane is considered within the framework of the long-wave approximation method. A mathematical model incorporating variation in surface tension with temperature has been formulated by coupling the Navier–Stokes equation, governing the flow, with the equation of energy. For the slippery substrate, the Navier slip boundary condition is used at the solid–liquid interface. An evolution equation is formed in terms of the free surface, which includes the effects of inertia, thermocapillary as well as slip length. Using the normal mode approach, linear stability analysis is carried out and a critical Reynolds number is obtained, which reflects its dependence on the Marangoni number Mn as well as slip length δ. This depicts that δ and Mn both have the destabilization effect on the flow field. The linear study also reveals that the inertia force has a negligible effect compare to the thermocapillary or slip. In addition, the study highlights a critical Marangoni number at which the instability sets in when the thermocapillary stress attains a critical value. The method of multiple scales is used to investigate the weakly nonlinear stability analysis of the flow. The study interprets that the variation of Mn and δ have substantial effects on different stable/unstable zones. It also shows that within a considered parametric domain, the unconditional stable zone completely vanishes for any value of Mn, when the slip length δ attains a critical value. The study also divulges that in the subcritical unstable (supercritical stable) zone the threshold amplitude (ζa) decreases (increases) with the increment of Mn and δ. Further, we discussed the spatial uniform solution of the complex Ginzburg–Landau equation for sideband disturbances. Employing the Crank–Nicolson method, the nonlinear evolution equation of the free surface is solved numerically in a periodic domain, considering the sinusoidal initial perturbation of small amplitude. The nonlinear simulations are found to be in good agreement with the linear and weakly nonlinear stability analysis. The evolution of the maximum hmax and minimum hmin thickness amplifies, for small change of Mn and δ. Further, it shows that the influence of the thermocapillary force amplifies the destabilizing nature of δ. The traveling wave solution confirms the existence of a fixed point for the considered parametric domain, chosen from the experimental data. Finally, the Hopf bifurcation of the fixed point exhibits that the nonlinear wave speed at the transcritical point increases as δ increases but decreases as Mn increases. The noteworthy result which arises from the study is that a transcritical Hopf bifurcation exists if the slip length δ>max16Mn−13,12Mn−23−Mn.

Well-posedness of stochastic third grade fluid equation

In this paper, we establish the well-posedness for the third grade fluid equation perturbed by a multiplicative white noise. This equation describes the motion of a non-Newtonian fluid of differential type with relevant viscoelastic properties. We are faced with a strongly nonlinear stochastic partial differential equation supplemented with a Navier slip boundary condition. Taking the initial condition in the Sobolev space H2, we show the existence and the uniqueness of the strong (in the probability sense) solution in a two dimensional and non axisymmetric bounded domain.

Breakup dynamics of capillary bridges on hydrophobic stripes

The breakup dynamics of a capillary bridge on a hydrophobic stripe between two hydrophilic stripes is studied experimentally and numerically using direct numerical simulations. The capillary bridge is formed from an evaporating water droplet wetting three neighboring stripes of a chemically patterned surface. By considering the breakup process in a phase space representation, the breakup dynamics can be evaluated without the uncertainty in determining the precise breakup time. The simulations are based on the Volume-of-Fluid (VOF) method implemented in Free Surface 3D (FS3D). In order to construct physically realistic initial data for the VOF simulation, Surface Evolver is employed to calculate an initial configuration consistent with experiments. Numerical instabilities at the contact line are reduced by a novel discretization of the Navier-slip boundary condition on staggered grids. The breakup of the capillary bridge cannot be characterized by a unique scaling relationship. Instead, at different stages of the breakup process different scaling exponents apply, and the structure of the bridge undergoes a qualitative change. In the final stage of breakup, the capillary bridge forms a liquid thread that breaks up consistently with the Rayleigh-Plateau instability.

Approximations for Steady Unidirectional Slip Flows in Elliptic Microchannels

Abstract Consider steady flows in a channel whose cross section Ω is an ellipse, flows with the Navier slip boundary condition. Denote the volume flow rate by Q. We apply to elliptic cross section a recent simple approximation, a rigorous lower bound R on Q, requiring, along with the channel's area and perimeter, the calculation of just the torsional rigidity and two other domain functionals. This avoids the need for solving the partial differential equation repeatedly for differing values of the slip parameter.

The immersed boundary-lattice Boltzmann method for solving solid-fluid interaction problem with Navier-slip boundary condition

To solve the solid-fluid interaction with Navier-slip boundary condition, we propose a new method in this paper. In this approach, the solid body is modeled by the Newton’s equations and the fluid is assumed to satisfy the discrete lattice Boltzmann equations. The Navier-slip boundary condition is applied as the velocity boundary condition for the solid-fluid interface. We discretize the Newton’s equations by a prediction/correction approach and the discrete lattice Boltzmann equations by a second-order finite difference scheme in space together with a second-order Range–Kutta scheme in time. The Navier-slip boundary condition and the immersed boundary body force are almost implicitly updated by using the implicit velocity correction-based immersed boundary-lattice Boltzmann method. Numerical simulations are carried out to verify the new presented method.

RBF collocation and hybrid-LHI methods for Stokes systems and its application to controllability problems

The purpose of this article is to introduce radial basis function (RBF) methods for solving both direct Stokes equations and controllability problems for the Stokes system with few internal scalar controls. In both cases, Dirichlet or Navier-slip boundary conditions are considered. We introduce two radial basis function solvers, one global and the other local, to solve Stokes equations. These methods are used to discretize the primal and adjoint systems related to the controllability problems. Both techniques are based on divergence-free global RBFs. A global colocation technique based on Div-free inverse multi-quadrics is formulated and analyzed. A generalization of scalar hybrid kernels to a vector divergence-free hybrid RBFs setting is defined. Based on these kernels, the local Hermite interpolation (LHI) method in vector form is introduced. Due to the properties of the hybrid kernel, we show that due to the properties of the hybrid kernel this local method, can reduce up to double precision, the value of the condition number of the local Gram matrices. Simultaneously, it is proved that the real components of the eigenvalues corresponding to the global LHI matrix are negative and that consequently backward difference formulas are stable for time integration. The conjugate gradient algorithm is adapted to the radial basis function setting to solve the controllability problems. Several benchmarks problems in two dimensions with a non-convex domain (a star shape) are numerically solved by these RBFs methods to display and compare their feasibility. The solutions to these problems are also implemented by finite element techniques to study their relative performance.

Local null controllability of a fluid–rigid body interaction problem with Navier slip boundary conditions

The aim of this work is to show the local null controllability of a fluid–solid interaction system by using a distributed control located in the fluid. The fluid is modeled by the incompressible Navier–Stokes system with Navier slip boundary conditions and the rigid body is governed by the Newton laws. Our main result yields that we can drive the velocities of the fluid and of the structure to 0 and we can control exactly the position of the rigid body. One important ingredient consists in a new Carleman estimate for a linear fluid–rigid body system with Navier boundary conditions. This work is done without imposing any geometrical conditions on the rigid body.

Asymptotic analysis of an optimal control problem for a viscous incompressible fluid with Navier slip boundary conditions

We consider an optimal control problem for the Navier–Stokes system with Navier slip boundary conditions. We denote by α the friction coefficient and we analyze the asymptotic behavior of such a problem as α → ∞. More precisely, we prove that if we take an optimal control for each α, then there exists a sequence of optimal controls converging to an optimal control of the same optimal control problem for the Navier–Stokes system with the Dirichlet boundary condition. We also show the convergence of the corresponding direct and adjoint states.

Global well-posedness of the Navier-Stokes equations with Navier-slip boundary conditions in a strip domain

<p style='text-indent:20px;'>This paper is concerned with the existence and uniqueness of the strong solution to the incompressible Navier-Stokes equations with Navier-slip boundary conditions in a two-dimensional strip domain where the slip coefficients may not have defined sign. In the meantime, we also establish a number of Gagliardo-Nirenberg inequalities in the corresponding Sobolev spaces which will be applicable to other similar situations.</p>

A partitioned numerical scheme for fluid–structure interaction with slip

We present a loosely coupled, partitioned scheme for solving fluid–structure interaction (FSI) problems with the Navier slip boundary condition. The fluid flow is modeled by the Navier–Stokes equations for an incompressible, viscous fluid, interacting with a thin elastic structure modeled by the membrane or Koiter shell type equations. The fluid and structure are coupled via two sets of coupling conditions: a dynamic coupling condition describing balance of forces, and a kinematic coupling condition describing fluid slipping tangentially to the moving fluid–structure interface, with no penetration in the normal direction. Problems of this type arise in, e.g., FSI with hydrophobic structures or surfaces treated with a no-stick coating, and in biologic FSI involving rough surfaces of elastic tissues or tissue scaffolds. We propose a novel, efficient partitioned scheme where the fluid sub-problem is solved separately from the structure sub-problem, and there is no need for sub-iterations at every time step to achieve stability, convergence, and its first-order accuracy. We derive energy estimates, which prove that the proposed scheme is unconditionally stable for the corresponding linear problem. Moreover, we present convergence analysis and show that under a time-step condition, the method is first-order accurate in time and optimally convergent in space for a Finite Element Method-based spatial discretization. The theoretical rates of convergence in time are confirmed numerically on an example with an explicit solution using the method of manufactured solutions, and on a benchmark problem describing propagation of a pressure pulse in a two-dimensional channel. The effects of the slip rate and fluid viscosity on the FSI solution are numerically investigated in two additional examples: a 2D cylindrical FSI example for which an exact Navier slip Poiseuille-type solution is found and used for comparison, and a squeezed ketchup bottle example with gravity enhanced flow. We show that the Navier-slip boundary condition increases the outflow mass flow rate by 21% for a bottle angled at 45 degrees pointing downward, in the direction of gravity.

Navier–Stokes Equations in a Curved Thin Domain, Part II: Global Existence of a Strong Solution

We consider the Navier-Stokes equations in a three-dimensional curved thin domain around a given closed surface under Navier's slip boundary conditions. When the thickness of the thin domain is sufficiently small, we establish the global existence of a strong solution for large data. We also show several estimates for the strong solution with constants explicitly depending on the thickness of the thin domain. The proofs of these results are based on a standard energy method and a good product estimate for the convection and viscous terms following from a detailed study of average operators in the thin direction. We use the average operators to decompose a three-dimensional vector field on the thin domain into the almost two-dimensional average part and the residual part, and derive good estimates for them which play an important role in the proof of the product estimate.