Heat transfer of nanofluids considering nanoparticle migration and second-order slip velocity

The heat transfer of a magnetohydrodynamics nanofluid inside an annulus considering the second-order slip condition and nanoparticle migration is theoret-ically investigated. A second-order slip condition, which appropriately represents the non-equilibrium region near the interface, is prescribed rather than the no-slip condition and the linear Navier slip condition. To impose different temperature gradients, the outer wall is subjected to q2, the inner wall is subjected to q1, and q1 〉 q2. A modified two-component four-equation non-homogeneous equilibrium model is employed for the nanofiuid, which have been reduced to two-point ordinary boundary value differential equations in the consideration of the thermally and hydrodynamically fully developed flow. The homotopy analysis method （HAM） is employed to solve the equations, and the h-curves are plotted to verify the accuracy and efficiency of the solutions. Moreover, the effects of the physical factors on the flow and heat transfer are discussed in detail, and the semi-analytical relation between NUB and NBT is obtained.

Effect of boundary slip on electroosmotic flow in a curved rectangular microchannel

The aim of this study is to numerically investigate the impact of boundary slip on electroosmotic flow(EOF) in curved rectangular microchannels. Navier slip boundary conditions were employed at the curved microchannel walls. The electric potential distribution was governed by the Poisson–Boltzmann equation, whereas the velocity distribution was determined by the Navier–Stokes equation. The finite-difference method was employed to solve these two equations. The detailed discussion focuses on the impact of the curvature ratio, electrokinetic width, aspect ratio and slip length on the velocity. The results indicate that the present problem is strongly dependent on these parameters. The results demonstrate that by varying the dimensionless slip length from 0.001 to 0.01 while maintaining a curvature ratio of 0.5 there is a twofold increase in the maximum velocity. Moreover, this increase becomes more pronounced at higher curvature ratios. In addition, the velocity difference between the inner and outer radial regions increases with increasing slip length. Therefore, the incorporation of the slip boundary condition results in an augmented velocity and a more non-uniform velocity distribution. The findings presented here offer valuable insights into the design and optimization of EOF performance in curved hydrophobic microchannels featuring rectangular cross-sections.

POSITIVE SOLUTIONS OF BOUNDARY VALUE PROBLEMS FOR SECOND-ORDER SINGULAR NONLINEAR DIFFERENTIAL EQUATIONS

New existence results are presented for the singular second-order nonlinear boundary value problems u ' + g(t)f(u) = 0, 0 < t < 1, au(0) - betau ' (0) = 0, gammau(1) + deltau ' (1) = 0 under the conditions 0 less than or equal to f(0)(+) < M-1, m(1) < f(infinity)(-)less than or equal to infinity or 0 less than or equal to f(infinity)(+)< M-1, m(1) < f (-)(0)less than or equal to infinity where f(0)(+) = lim(u -->0)f(u)/u, f(infinity)(-)= lim(u --> infinity)f(u)/u, f(0)(-)= lim(u -->0)f(u)/u, f(infinity)(+) = lim(u --> infinity)f(u)/u, g may be singular at t = 0 and/or t = 1. The proof uses a fixed point theorem in cone theory.

Second-Order Adjoint Sensitivity Analysis Methodology for Computing Exactly Response Sensitivities to Uncertain Parameters and Boundaries of Linear Systems: Mathematical Framework

This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2nd-CASAM)” for the efficient and exact computation of 1st- and 2nd-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (i.e., model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2nd-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1st-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2nd-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (N2/2 + 3 N/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2nd-LASS requires very little additional effort beyond the construction of the 1st-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1st-LASS and 2nd-LASS requires the same computational solvers as needed for solving (i.e., “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1st-LASS and the 2nd-LASS. Since neither the 1st-LASS nor the 2nd-LASS involves any differentials of the operators underlying the original system, the 1st-LASS is designated as a “first-level” (as opposed to a “first-order”) adjoint sensitivity system, while the 2nd-LASS is designated as a “second-level” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2nd-LASS that involve the imprecisely known boundary parameters. Notably, the 2nd-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.

Three-dimensional mixed convection stagnation-point fow past a vertical surface with second-order slip velocity

This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is the first study presenting this very interesting analysis.Nonlinear partial differential equations for the flow problem are transformed into nonlinear ordinary differential equations(ODEs)by using appropriate similarity transformation.These ODEs with the corresponding boundary conditions are numerically solved by utilizing the bvp4c solver in MATLAB programming language.The effects of the governing parameters on the non-dimensional velocity profiles,temperature profiles,skin friction coefficients,and the local Nusselt number are presented in detail through a series of graphs and tables.Interestingly,it is reported that the reduced skin friction coefficient decreases for the assisting flow situation and increases for the opposing flow situation.The numerical computations of the present work are compared with those from other research available in specific situations,and an excellent consensus is observed.Another exciting feature for this work is the existence of dual solutions.An important remark is that the dual solutions exist for both assisting and opposing flows.A linear stability analysis is performed showing that one solution is stable and the other solution is not stable.We notice that the mixed convection and velocity slip parameters have strong effects on the flow characteristics.These effects are depicted in graphs and discussed in this paper.The obtained results show that the first-order and second-order slip parameters have a considerable effect on the flow,as well as on the heat transfer characteristics.

Positive Solutions for a Class of Second-order m-point Boundary Value Problems

Using a fixed point theorem in cones, the paper consider the existence of positive solutions for a class of second-order m-point boundary value problem. Sufficient conditions to ensure the existence of double positive solutions are obtained. The associated Green function of this problem is also given.

Positive solutions of second-order differential systems with nonhomogeneous boundary conditions

This paper is to investigate the positive solutions of the systems of second-order ordinary differential equations with nonhomogeneous multi-point boundary conditions. By the lower and upper solutions method, Schauder fixed point theorem and fixed point index theory, under certain conditions, it is proved that there exist appropriate regions of parameters in which the problem has at least two, at least one or no positive solution.

Statistical second-order two-scale analysis and computation for heat conduction problem with radiation boundary condition in porous materials

This paper discusses a statistical second-order two-scale（SSOTS） analysis and computation for a heat conduction problem with a radiation boundary condition in random porous materials.Firstly,the microscopic configuration for the structure with random distribution is briefly characterized.Secondly,the SSOTS formulae for computing the heat transfer problem are derived successively by means of the construction way for each cell.Then,the statistical prediction algorithm based on the proposed two-scale model is described in detail.Finally,some numerical experiments are proposed,which show that the SSOTS method developed in this paper is effective for predicting the heat transfer performance of porous materials and demonstrating its significant applications in actual engineering computation.

ON SOLUTION TO THE NAVIER-STOKES EQUATIONS WITH NAVIER SLIP BOUNDARY CONDITION FOR THREE DIMENSIONAL INCOMPRESSIBLE FLUID

In this article, we prove the existence and uniqueness of solutions of the NavierStokes equations with Navier slip boundary condition for incompressible fluid in a bounded domain of R^3. The results are established by the Galerkin approximation method and improved the existing results.

THE NAVIER-STOKES EQUATIONS WITH THE KINEMATIC AND VORTICITY BOUNDARY CONDITIONS ON NON-FLAT BOUNDARIES

We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition （without the incompressibility condition）, which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in R^n（n ≥ 3） with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.

Effect of Pressure on Boundary Slip of Thin Film Lubrication Using Atomistic Simulation

Mechanical systems on all length scales may be subjected to nanoscale thin film lubrication(TFL). Molecular dynamics(MD) simulations were conducted to investigate the lubrication mechanism and boundary slip of squalane confined in nanogap at 293 K with two different film thicknesses and a wide range of pressures. The molecular distribution, density and velocity profiles of squalane were analyzed. The results show that the lubricant atoms tend to form layers parallel to the wall, but the lubricant molecules orient randomly throughout the film in the directions both parallel and perpendicular to the wall. Most squalane molecules appear twisted and folded, and extend to several atomic layers so that there are no slips between lubricant layers. The distances between the lubricant layers are irregular rather than broadening far away from the walls. The boundary slip at the interface of bcc Fe(001) and squalane only occurs at high pressure because of the strong nonbond interactions between lubricant atoms and wall atoms. The tendency of boundary slip is more obvious for films with thinner film thickness. According to the simulations, the relationship between the slip length and the pressure is given.

Scaling group transformation for MHD boundary layer flow over permeable stretching sheet in presence of slip flow with Newtonian heating effects

Taking into account the slip flow effects, Newtonian heating, and thermal radiation, two-dimensional magnetohydrodynamic （MHD） flows and heat transfer past a permeable stretching sheet are investigated numerically. We use one parameter group transformation to develop similarity transformation. By using the similarity transformation, we transform the governing boundary layer equations along with the boundary conditions into ordinary differential equations with relevant boundary conditions. The obtained ordinary differential equations are solved with the fourth-fifth order Runge-Kutta- Fehlberg method using MAPLE 13. The present paper is compared with a published one. Good agreement is obtained. Numerical results for dimensionless velocity, temperature distributions, skin friction factor, and heat transfer rates are discussed for various values of controlling parameters.

Effects of hysteresis of static contact angle (HSCA) and boundary slip on the hydrodynamics of water striders

It is known that contact lines keep relatively still on solids until static contact angles exceed an interval of hysteresis of static contact angle (HSCA), and contact angles keep changing as contact lines relatively slide on the solid. Here, the effects of HSCA and boundary slip were first distinguished on the micro-curvature force (MCF) on the seta. Hence, the total MCF is partitioned into static and dynamic MCFs correspondingly. The static MCF was found proportional to the HSCA and related with the asymmetry of the micro-meniscus near the seta. The dynamic MCF, exerting on the relatively sliding contact line, is aroused by the boundary slip. Based on the Blake-Haynes mechanism, the dynamic MCF was proved important for water walking insects with legs slower than the minimum wave speed . As insects brush the water by laterally swinging legs backwards, setae on the front side of the leg are pulled and the ones on the back side are pushed to cooperatively propel bodies forward. If they pierce the water surface by vertically swinging legs downwards, setae on the upside of the legs are pulled, and the ones on the downside are pushed to cooperatively obtain a jumping force. Based on the dependency between the slip length and shear rate, the dynamic MCF was found correlated with the leg speed U, as , where and are determined by the dimple depth. Discrete points on this curve could give fitted relations as (Suter et al., J. Exp. Biol. 200, 2523-2538, 1997). Finally, the axial torque on the inclined and partially submerged seta was found determined by the surface tension, contact angle, HSCA, seta width, and tilt angle. The torque direction coincides with the orientation of the spiral grooves of the seta, which encourages us to surmise it is a mechanical incentive for the formation of the spiral morphology of the setae of water striders.

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.

An improved convolution perfectly matched layer for elastic second-order wave equation

A convolution perfectly matched layer(CPML)can efficiently absorb boundary reflection in numerical simulation.However,the CPML is suitable for the first-order elastic wave equation and is difficult to apply directly to the second-order elastic wave equation.In view of this,based on the first-order CPML absorbing boundary condition,we propose a new CPML(NCPML)boundary which can be directly applied to the second-order wave equation.We first systematically extend the first-order CPML technique into second-order wave equations,neglecting the space-varying characteristics of the partial damping coefficient in the complex-frequency domain,avoiding the generation of convolution in the time domain.We then transform the technique back to the time domain through the inverse Fourier transform.Numerical simulation indicates that the space-varying characteristics of the attenuation factor have little influence on the absorption effect and increase the memory at the same time.A number of numerical examples show that the NCPML proposed in this study is effective in simulating elastic wave propagation,and this algorithm is more efficient and requires less memory allocation than the conventional PML absorbing boundary.

SQUEEZE FILM FLOW WITH NONLINEAR BOUNDARY SLIP

A nonlinear boundary slip model consisting of an initial slip length and a critical shear rate was used to study the nonlinear boundary slip of squeeze fluid film confined between two approaching spheres. It is found that the initial slip length controls the slip behavior at small shear rate, but the critical shear rate controls the boundary slip at high shear rate. The boundary slip at the squeeze fluid film of spherical surfaces is a strongly nonlinear function of the radius coordinate. At the center or far from the center of the squeeze film, the slip length equals the initial slip length due to the small shear rate. However, in the high shear rate regime the slip length increases very much. The hydrodynamic force of the spherical squeeze film decreases with increasing the initial slip length and decreasing the critical shear rate. The effect of initial slip length on the hydrodynamic force seems less than that of the critical shear rate. When the critical shear rate is very small the hydrodynamic force increases very slowly with a decrease in minimum film thickness. The theoretical predictions agree well with the experiment measurements.

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 lu- brication 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.

Effects of boundary slippage on thin-film lubricationbetween two nonparallel plane plates

Hydrodynamic lubrications between two plane plates with an intersection angle θ have been investigated us- ing the boundary slippage theory, and relations are obtained between dimensionless pressures and coordinate x, between bearing capacity, friction force, friction coefficient and di- mensionless slipping size factor. The results show that bear- ing capacity of two plane plates without boundary slippage significantly increases with increasing intersection angle θ when 0 〈 θ 〈 1°, whereas decreases with increasing in- tersection angle 0 when θ 〉 1°. The results also show that negative pressure occurs in fluid entrance region and bearing capacity decreases, and friction force and friction coefficient increase with the increase of dimensionless slipping size fac- tor.

Boundary Layer Stagnation-Point Slip Flow and Heat Transfer towards a Shrinking/Stretching Cylinder over a Permeable Surface

In this paper, the boundary layer stagnation-point slip flow and heat transfer towards a shrinking/stretching cylinder over a permeable surface is considered. The governing equations are first transformed into a system of non-dimensional equations via the non-dimensional variables, and then into self-similar ordinary differential equations before they are solved numerically using the shooting method. Numerical results are obtained for the skin friction coefficient and the local Nusselt number as well as the velocity and temperature profiles for some values of the governing parameters, namely the velocity slip parameter (α), the thermal slip parameter (β), the curvature parameter (γ) and the velocity ratio parameter (c/a). The physical quantities of interest are the skin friction coefficient and the local Nusselt number measured by f’’(0) and –θ’(0), respectively. The numerical results show that the velocity slip parameter α increases the heat transfer rate at the surface, while the thermal slip parameter β decreases it. On the other hand, increasing the velocity slip parameter α causes the decrease in the flow velocity. Further, it is found that the solutions for a shrinking cylinder (c/ac/a>0) case. Finally, it is also found that the values of f’’(0) and –θ’(0) increase as the curvature parameter γ increases.

Local Existence for Boussinesq Equations with Slip Boundary Condition in a Bounded Domain

In this paper, we are interested in the local existence for the Boussinesq equations with the slip boundary conditions. Energy method and Gerlakin approach are employed in this paper to get the main result.