Influence of hydrodynamic slip on convective transport in flow past a circular cylinder

The presence of a finite tangential velocity on a hydrodynamically slipping surface is known to reduce vorticity production in bluff body flows substantially while at the same time enhancing its convection downstream and into the wake. Here, we investigate the effect of hydrodynamic slippage on the convective heat transfer (scalar transport) from a heated isothermal circular cylinder placed in a uniform cross-flow of an incompressible fluid through analytical and simulation techniques. At low Reynolds (\({\textit{Re}}\ll 1\)) and high Péclet (\({\textit{Pe}}\gg 1\)) numbers, our theoretical analysis based on Oseen and thermal boundary layer equations allows for an explicit determination of the dependence of the thermal transport on the non-dimensional slip length \(l_s\). In this case, the surface-averaged Nusselt number, Nu transitions gradually between the asymptotic limits of \(Nu \sim {\textit{Pe}}^{1/3}\) and \(Nu \sim {\textit{Pe}}^{1/2}\) for no-slip (\(l_s \rightarrow 0\)) and shear-free (\(l_s \rightarrow \infty \)) boundaries, respectively. Boundary layer analysis also shows that the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) holds for a shear-free cylinder surface in the asymptotic limit of \({\textit{Re}}\gg 1\) so that the corresponding heat transfer rate becomes independent of the fluid viscosity. At finite \({\textit{Re}}\), results from our two-dimensional simulations confirm the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) for a shear-free boundary over the range \(0.1 \le {\textit{Re}}\le 10^3\) and \(0.1\le {\textit{Pr}}\le 10\). A gradual transition from the lower asymptotic limit corresponding to a no-slip surface, to the upper limit for a shear-free boundary, with \(l_s\), is observed in both the maximum slip velocity and the Nu. The local time-averaged Nusselt number \(Nu_{\theta }\) for a shear-free surface exceeds the one for a no-slip surface all along the cylinder boundary except over the downstream portion where unsteady separation and flow reversal lead to an appreciable rise in the local heat transfer rates, especially at high \({\textit{Re}}\) and Pr. At a Reynolds number of \(10^3\), the formation of secondary recirculating eddy pairs results in appearance of additional local maxima in \(Nu_{\theta }\) at locations that are in close proximity to the mean secondary stagnation points. As a consequence, Nu exhibits a non-monotonic variation with \(l_s\) increasing initially from its lowermost value for a no-slip surface and then decreasing before rising gradually toward the upper asymptotic limit for a shear-free cylinder. A non-monotonic dependence of the spanwise-averaged Nu on \(l_s\) is observed in three dimensions as well with the three-dimensional wake instabilities that appear at sufficiently low \(l_s\), strongly influencing the convective thermal transport from the cylinder. The analogy between heat transfer and single-component mass transfer implies that our results can directly be applied to determine the dependency of convective mass transfer of a single solute on hydrodynamic slip length in similar configurations through straightforward replacement of Nu and \({\textit{Pr}}\) with Sherwood and Schmidt numbers, respectively.     

Linearized Stability for Semiflows Generated by a Class of Neutral Equations,with Applications to State-Dependent Delays

We prove a principle of linearized stability for semiflows generated by neutral functional differential equations of the form x′(t) = g(? x _t , x _t ). The state space is a closed subset in a manifold of C ²-functions. Applications include equations with state-dependent delay, as for example x′(t) = a x′(t + d(x(t))) + f (x(t + r(x(t)))) with \({a\in\mathbb{R}, d:\mathbb{R}\to(-h,0), f:\mathbb{R}\to\mathbb{R}, r:\mathbb{R}\to[-h,0]}\).     

Parametric study on heat transfer enhancement and pressure drop of an internal blade tip-wall with pin-fin arrays

One way to cool gas turbine tips is to design serpentine passages with 180° turns inside the blades to fully utilize the coolant potential. It is therefore a desire to improve the cooling of the blade tips to ensure a long durability and safe operation. In the present work, a two-pass channel with a 180° turn and various arrays of pin-fins mounted internally on the tip-cap is considered. The effects of pin-fin height, diameter and pitches on the heat transfer enhancement and pressure drop are investigated numerically. The nominal ratio of height to diameter (H/D) of the pin-fins is 2, and the ratio of tip clearance to pin-fin height is about 10. The inlet Reynolds numbers based on hydraulic diameter are ranging from 100,000 to 600,000. Details of the three dimensional fluid flow and heat transfer over the pin-finned tips are presented. The overall performances of various tips are compared. It is found that due to the combination of turning, impingement and pin-fin crossflow, the heat transfer coefficient of the pin-finned tips is up to a factor of 2.1 higher than that of the smooth tip. This augmentation is achieved at the expense of a penalty of pressure drop around 30%. Results show that the magnitude of the heat transfer enhancement depends upon pin-fin configuration and arrangement. It is suggested that pin-fins are suitable to enhance the blade tip heat transfer and thus to improve the tip cooling.     

Miscible Thermo-Viscous Fingering Instability in Porous Media. Part 1: Linear Stability Analysis

The development of the thermo-viscous fingering instability of miscible displacements in homogeneous porous media is examined. In this first part of the study dealing with stability analysis, the basic equations and the parameters governing the problem in a rectilinear geometry are developed. An exponential dependence of viscosity on temperature and concentration is represented by two parameters, thermal mobility ratio β _T and a solutal mobility ratio β _C , respectively. Other parameters involved are the Lewis number Le and a thermal-lag coefficient λ. The governing equations are linearized and solved to obtain instability characteristics using either a quasi-steady-state approximation (QSSA) or initial value calculations (IVC). Exact analytical solutions are also obtained for very weakly diffusing systems. Using the QSSA approach, it was found that an increase in thermal mobility ratio β _T is seen to enhance the instability for fixed β _C , Le and λ. For fixed β _C and β _T , a decrease in the thermal-lag coefficient and/or an increase in the Lewis number always decrease the instability. Moreover, strong thermal diffusion at large Le as well as enhanced redistribution of heat between the solid and fluid phases at small λ is seen to alleviate the destabilizing effects of positive β _T . Consequently, the instability gets strictly dominated by the solutal front. The linear stability analysis using IVC approach leads to conclusions similar to the QSSA approach except for the case of large Le and unity λ flow where the instability is seen to get even less pronounced than in the case of a reference isothermal flow of the same β _C , but β _T  = 0. At practically, small value of λ, however, the instability ultimately approaches that due to β _C only.     

Inviscid Models Generalizing the Two-dimensional Euler and the Surface Quasi-geostrophic Equations

Any classical solution of the two-dimensional incompressible Euler equation is global in time. However, it remains an outstanding open problem whether classical solutions of the surface quasi-geostrophic (SQG) equation preserve their regularity for all time. This paper studies solutions of a family of active scalar equations in which each component u _j of the velocity field u is determined by the scalar θ through \({u_j =\mathcal{R}\Lambda^{-1}P(\Lambda) \theta}\) , where \({\mathcal{R}}\) is a Riesz transform and Λ = (?Δ)^1/2. The two-dimensional Euler vorticity equation corresponds to the special case P(Λ) = I while the SQG equation corresponds to the case P(Λ) = Λ. We develop tools to bound \({\|\nabla u||_{L^\infty}}\) for a general class of operators P and establish the global regularity for the Loglog-Euler equation for which P(Λ) = (log(I + log(I ? Δ))) ^γ with 0 ≦ γ ≦ 1. In addition, a regularity criterion for the model corresponding to P(Λ) = Λ ^β with 0 ≦ β ≦ 1 is also obtained.     

Families of rational solutions of the <Emphasis Type="Italic">y</Emphasis>-nonlocal Davey–Stewartson II equation

The y-nonlocal Davey–Stewartson II equation is an extension of the usual DS II equation involving a partially parity-time-symmetric potential only with respect to the spatial variable y. By using the Hirota bilinear method, families of n-order rational solutions are obtained, which include lumps in the (x, y)-plane and the (y, t)-plane, growing-and-decaying line waves in the (x, t)-plane, and hybrid solutions of interacting line rogue waves and lumps in the (x, y)-plane.     

Three-Dimensional Problem of Radiative Gasdynamics of the Apollo-4 Command Module during Superorbital Atmospheric Entry

The three-dimensional problem of radiative gasdynamics of the superorbital entry of the Apollo-4 command module into the dense terrestrial atmosphere at an angle of attack of 25° is numerically solved. The flow conditions corresponding to the flight velocity V_∞ = 10.5 km/s at an altitude H = 67.3 km are considered in detail. The distributions of the densities of convective and radiative heat fluxes along the surface in a flow are obtained. The spectral composition of the thermal radiation attaining the surface is studied. The results of the calculations are successfully compared with the data of two-dimensional calculations.     

Optimization of a stepped circular pin-fin array to enhance heat transfer performance

A Stepped circular pin-fin array is formulated numerically and optimized with Kriging metamodeling technique to enhance heat transfer performance. The problem is defined by two non-dimensional geometric design variables composed of height of the channel, height of smaller diameter part of the pin-fins, and smaller diameter of the pin-fins, to maximize heat transfer rate compromising with friction loss. Ten designs generated by Latin hypercube sampling were evaluated by three-dimensional Reynolds-averaged Navier–Stokes solver and the evaluated objectives were used to construct the surrogate model. The predictions of objective function by Kriging model at optimum point show reasonable accuracy in comparison with the values calculated by RANS analysis. Optimum shape of pin-fins strongly depends on the weighting factor which measures importance of the friction loss term in the objective function. The thermal performances are much higher than that of the straight pin-fin at sampling optimum points with different weighting factors.     

Multipulse Phases in k-Mixtures of Bose–Einstein Condensates

For the system

$-\Delta U_i+ U_i=U_i^3-\beta U_i\sum_{j\neq i}U_j^2,\quad i=1,\dots,k,$

(with k ≧ 3), we prove the existence for β large of positive radial solutions on \({\mathbb R^N}\) . We show that as β →  + ∞, the profile of each component U _i separates, in many pulses, from the others. Moreover, we can prescribe the location of such pulses in terms of the oscillations of the changing-sign solutions of the scalar equation  ? ΔW  +  W  =  W³. Within an Hartree–Fock approximation, this provides a theoretical indication of phase separation into many nodal domains for the k-mixtures of Bose–Einstein condensates.

     

Mechanical Properties of Hard W-C Coating on Steel Substrate Deduced from Nanoindentation and Finite Element Modeling

Mechanical properties of a hard and stiff W-C coating on steel substrate have been investigated using nanoindentation combined with finite element modeling (FEM) and extended FEM (XFEM). The significant pile-up observed around the indents in steel substrate caused an overestimation of hardness and indentation modulus. A simple geometrical model, considering the additional contact surfaces due to pile-up, has been proposed to reduce this overestimation. The presence of W-C coating suppressed the pile-up in the steel substrate and a transition to sink-in behavior occurred. The FEM simulations adequately reproduced the surface topography of the indents in the substrate and coating/substrate systems as well. The maximum principal stresses of the indented W-C/steel coated system were tensile; they were always located in the coating and evolved in 3 stages. Cohesive cracking occurred during loading in the sink-in zone (stage III) when the ultimate tensile strength (σ _max ) of the coating was reached. The obtained hardness (H _c ), indentation modulus (E _c ), yield stress (Y) and strength (σ _max ) of the W-C coating were H _c ? =?20 GPa, E _c ? =?250 GPa, Y?=?9.0 GPa and σ _max ? =?9.35 GPa, respectively. XFEM resulted in fracture energy of the W-C coating of G?=?38.1 J?·?m^-2 and fracture toughness of K _IC ? =?3.5 MPa?·?m^0.5.     

Numerical Simulation of Two Phase Turbulent Flow of Nanofluids in Confined Slot Impinging Jet

In this article, a numerical investigation is performed on flow and heat transfer of confined impinging slot jet, with a mixture of water and Al₂O₃ nanoparticles as the working fluid. Two-dimensional turbulent flow is considered and a constant temperature is applied on the impingement surface. The k ? ω turbulence model is used for the turbulence computations. Two-phase mixture model is implemented to study such a flow field. The governing equations are solved using the finite volume method. In order to consider the effect of obstacle angle on temperature fields in the channel, the numerical simulations were performed for different obstacle angles of 0° ? 60°. Also different geometrical parameters, volume fractions and Reynolds numbers have been considered to study the behavior of the system in terms of stagnation point, average and local Nusselt number and stream function contours. The results showed that the intensity and size of the vortex structures depend on jet- impingement surface distance ratio (H/W) and volume fraction. The maximum Nusselt number occurs at the stagnation point with the highest values at about H/W = 1. Increasing obstacle angle, from 15° to 60°, enhances the heat transfer rate. It was also revealed that the minimum value of average Nusselt number occurs in higher H/W ratios with decreasing the channel length.     

Rheology and structural changes of plasticized zeins in the molten state

Zeins, storage proteins from maize, are suitable for making biobased thermoplastic materials. The rheological behavior of a commercial zein plasticized with 20 w% glycerol was studied in the molten state by steady-state flow experiments in extrusion conditions and oscillatory rheometry. For low residence times, a shear-thinning viscoelastic behavior was observed, with G″ exceeding G′. After 300 s at 130 °C, the complex viscosity |η ^?|?=?7?×?10³ ω ^?0.46 was found to be similar to that of thermoplastic polymer melts used in fused deposition modeling. However, the ratio between the exponents of the power laws describing G′(ω) and G″(ω) did not meet the typical value of 2 for entangled polymer melts. Moreover, for longer residence times, the viscosity increased and a gelation phenomenon was observed with a crossing over of G′(ω) and G″(ω). Gel times ranged from 6000 s at 120 °C to 1700 s at 150 °C. The evolution of the macromolecular structure assessed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and high-performance size exclusion chromatography suggested that this gelation phenomenon involves various types of covalent and non-covalent cross-links. Disulfide bonds played a significant role in gelation kinetics despite a very low cysteine residue content in the protein primary structure (about 1 mol%). These results suggested that plasticized zeins initially behave like a low-viscosity non-entangled polymer melt, before cross-linking progressively led to a continuous network.     

Growths of mechanical elasticity and electrical conductance of graphene nanoplatelet/poly(lactic acid) composites under strong electric field: correlation with time evolution of higher order structure of graphene nanoplatelets

For a composite of poly(lactic acid) containing graphene nanoplatelets (GNPs) at a low concentration (0.34 vol%), this study examined growths of mechanical and electrical properties under an alternating current (AC) electric field, focusing on field-induced GNP structures governing those properties. The composite was subjected to the AC field (60 Hz) of various intensities E for various times t _E at 190 °C. A fraction of randomly oriented GNPs was aligned by the field and then connected into columns, as suggested from optical microscopy. This structural evolution led to qualitatively similar growths of low-frequency storage modulus and static electrical conductivity. The key quantity for understanding this growth was a time t _E* for occurrence of short circuit that detected formation of GNP columns conductively bridging the electrodes. The growths of both modulus and conductivity for various E were summarized as functions of a reduced variable, t _E/t _E*, confirming the growths commonly reflected the evolution of the GNP columns. However, the modulus grew fast and leveled off by t _E/t _E* ~ 1, whereas the conductivity kept growing gradually even at t _E/t _E* > 1. This difference was discussed in relation to the matrix chains and leftover GNPs out the column.     

Planar Vector Fields with a Given Set of Orbits

We determine all the \({\mathcal{C}^1}\) planar vector fields with a given set of orbits of the form y ? y(x) = 0 satisfying convenient assumptions. The case when these orbits are branches of an algebraic curve is also study. We show that if a quadratic vector field admits a unique irreducible invariant algebraic curve \({g(x, y) = \sum_{j=0}^S a_j(x) y^{S-j}= 0}\) with S branches with respect to the variable y, then the degree of the polynomial g is at most 4S.     

Two-dimensional interaction between an incident shock and a turbulent boundary layer in the presence of an entropy layer

The results of an experimental and numerical investigation of flow and heat transfer in the region of the interaction between an incident oblique shock and turbulent boundary layers on sharp and blunt plates are presented for the Mach numbers M_∞ = 5 and 6 and the Reynolds numbers Re_∞L = 27×10⁶ and 14×10⁶. The plate bluntness and the incident shock position were varied. It is shown that the maximum Stanton number St _m in the shock incidence zone decreases with increase in the plate bluntness radius r to a certain value and then varies only slightly with further increase in r. In the case of a turbulent undisturbed boundary layer heat transfer is diminished with increase in r more slowly than in the case of a laminar undisturbed flow. In the presence of an incident shock the bluntness of the leading edge of the flat plate results in a greater decrease in the Stanton number than in the absence of the shock. With increase in the bluntness of the leading edge of the plate the separation zone first sharply lengthens and then decreases in size or remains constant.     

Regularity and Singularities of Optimal Convex Shapes in the Plane

We focus here on the analysis of the regularity or singularity of solutions Ω ₀ to shape optimization problems among convex planar sets, namely:

$J(\Omega_{0})={\rm min} \{J(\Omega), \Omega \quad {\rm convex},\Omega \in \mathcal{S}_{\rm ad}\},$

where \({\mathcal{S}_{\rm ad}}\) is a set of 2-dimensional admissible shapes and \({J:\mathcal{S}_{\rm ad}\rightarrow\mathbb{R}}\) is a shape functional. Our main goal is to obtain qualitative properties of these optimal shapes by using first and second order optimality conditions, including the infinite dimensional Lagrange multiplier due to the convexity constraint. We prove two types of results:

1. i)
   under a suitable convexity property of the functional J, we prove that Ω ₀ is a W ^2,p -set, \({p\in[1, \infty]}\). This result applies, for instance, with p = ∞ when the shape functional can be written as J(Ω) = R(Ω) + P(Ω), where R(Ω) = F(|Ω|, E _f (Ω), λ₁(Ω)) involves the area |Ω|, the Dirichlet energy E _f (Ω) or the first eigenvalue of the Laplace–Dirichlet operator λ₁(Ω), and P(Ω) is the perimeter of Ω;
    

1. ii)
   under a suitable concavity assumption on the functional J, we prove that Ω ₀ is a polygon. This result applies, for instance, when the functional is now written as J(Ω) = R(Ω) ? P(Ω), with the same notations as above.
    

     

Regularity of Optimal Maps on the Sphere: the Quadratic Cost and the Reflector Antenna

Building on the results of Ma et al. (in Arch. Rational Mech. Anal. 177(2), 151–183 (2005)), and of the author Loeper (in Acta Math., to appear), we study two problems of optimal transportation on the sphere: the first corresponds to the cost function d ²(x, y), where d(·, ·) is the Riemannian distance of the round sphere; the second corresponds to the cost function ?log |x ? y|, known as the reflector antenna problem. We show that in both cases, the cost-sectional curvature is uniformly positive, and establish the geometrical properties so that the results of Loeper (in Acta Math., to appear) and Ma et al. (in Arch. Rational Mech. Anal. 177(2), 151–183 (2005)) can apply: global smooth solutions exist for arbitrary smooth positive data and optimal maps are Hölder continuous under weak assumptions on the data.     

Intensity of convective mass/heat transfer in a rotary regenerator rotor with the transverse needle-fins

A forced convective mass transfer coefficient was electrochemically measured for a cylindrical bundle of transverse needle-fins ?1 × 10.9, applied as the rotor porous matrix of a rotary heat regenerator. The baffle inside the rotor was present. The technique based on the ferricyanide–ferrocyanide redox reaction controlled at the cathode, in the presence of a sodium hydroxide based electrolyte, was used in this experiment. A set of the six neighbouring fins, connected in parallel, was the cathode. The distribution of the mass transfer coefficient according to different static rotor angle position and the mean mass transfer Chilton–Colburn coefficient correlation j _M  = j _M (Re) for rotation numbers, Ro: 0, 0.8, 1.6 and 2.0 were stated in the mean Reynolds number, Re, range 180–985. The comparison was made between the convective heat fluxes of the pin-fins and the sheet rotor, for Ro = 0.     

Simultaneous impacts of MHD and variable wall temperature on transient mixed Casson nanofluid flow in the stagnation point of rotating sphere

A numerical analysis is provided to scrutinize time-dependent magnetohydrodynamics(MHD) free and forced convection of an electrically conducting non-Newtonian Casson nanofluid flow in the forward stagnation point region of an impulsively rotating sphere with variable wall temperature. A single-phase flow of nanofluid model is reflected with a number of experimental formulae for both effective viscosity and thermal conductivity of nanofluid. Exceedingly nonlinear governing partial differential equations(PDEs)subject to their compatible boundary conditions are mutated into a system of nonlinear ordinary differential equations(ODEs). The derived nonlinear system is solved numerically with implementation of an implicit finite difference procedure merging with a technique of quasi-linearization. The controlled parameter impacts are clarified by a parametric study of the entire flow regime. It is depicted that from all the exhibited nanoparticles,Cu possesses the best convection. The surface heat transfer and surface shear stresses in the x-and z-directions are boosted with maximizing the values of nanoparticle solid volume fraction ? and rotation λ. Besides, as both the surface temperature exponent n and the Casson parameter γ upgrade, an enhancement of the Nusselt number is given.     

Monte Carlo Assessment of the Impact of Oscillatory and Pulsating Boundary Conditions on the Flow Through Porous Media

The linear stability analysis of vertical throughflow of power law fluid for double-diffusive convection with Soret effect in a porous channel is investigated in this study. The upper and lower boundaries are assumed to be permeable, isothermal and isosolutal. The linear stability of vertical through flow is influenced by the interactions among the non-Newtonian Rayleigh number (Ra), Buoyancy ratio (N), Lewis number (Le), Péclet number (Pe), Soret parameter (Sr) and power law index (n). The results indicate that the Soret parameter has a significant influence on convective instability of power law fluid. It has also been noticed that buoyancy ratio has a dual effect on the instability of fluid flow. Further, it is noticed that the basic temperature and concentration profiles have singularities at \(Pe = 0\) and \(Le = 1\), the convective instability is looked into for the limiting case of \(Pe\rightarrow 0\) and \(Le \rightarrow 1\). For the case of pure thermal convection with no vertical throughflow, the present numerical results coincide with the solution of standard Horton–Rogers–Lapwood problem. The present results for critical Rayleigh number obtained using bvp4c and two-term Galerkin approximation are compared with those available in the literature and are tabulated.