Slip-flow heat transfer in a microchannel with viscous dissipation

Forced convection flow in a microchannel with constant wall temperature is studied, including viscous dissipation effect. The slip-flow regime is considered by incorporating both the velocity-slip and the temperature-jump conditions at the surface. The energy equation is solved for the developing temperature field using finite integral transform. To increase β_v Kn is to increase the slip velocity at the wall surface, and hence to decrease the friction factor. Effects of the parameters β_v Kn, β, and Br on the heat transfer results are illustrated and discussed in detail. For a fixed Br, the Nusselt number may be either higher or lower than those of the continuum regime, depending on the competition between the effects of β_v Kn and β. At a given β_v Kn the variation of local Nusselt number becomes more even when β becomes larger, accompanied by a shorter thermal entrance length. The fully developed Nusselt number decreases with increasing β irrelevant to β_v Kn. The increase in Nusselt number due to viscous heating is found to be more pronounced at small β_v Kn.     

On the unsteady rotational flow of a fractional second grade fluid through a circular cylinder

Here the velocity field and the associated tangential stress corresponding to the rotational flow of a generalized second grade fluid within an infinite circular cylinder are determined by means of the Laplace and finite Hankel transforms. At time t=0 the fluid is at rest and the motion is produced by the rotation of the cylinder around its axis. The solutions that have been obtained are presented under series form in terms of the generalized G-functions. The similar solutions for ordinary second grade and Newtonian fluids are obtained from general solution for β→1, respectively, β→1 and α ₁→0. Finally, the influences of the pertinent parameters on the fluid motion, as well as a comparison between models, is underlined by graphical illustrations.     

On the unsteady rotational flow of a generalized Maxwell fluid through a circular cylinder

In this paper, the velocity field and the associated tangential stress corresponding to the rotational flow of a generalized Maxwell fluid within an infinite circular cylinder are determined by means of the Laplace and finite Hankel transforms. Initially, the fluid is at rest, and the motion is produced by the rotation of the cylinder about its axis with a unsteady angular velocity. The solutions that have been obtained are presented under series form in terms of the generalized G _a,b,c (, t)-functions. The similar solutions for the ordinary Maxwell and Newtonian fluids, performing the same motion, are obtained as special cases, when β → 1, respectively β → 1 and λ → 0, from general solutions. Finally, the solutions that have been obtained are compared by graphical illustrations, and the influence of the pertinent parameters on the fluid motion is also underlined by graphical illustrations.     

Flow of power-law fluids over a moving wedge surface with wall mass injection

The boundary layer problem of a power-law fluid flow with fluid injection on a wedge whose surface is moving with a constant velocity in the opposite direction to that of the uniform mainstream is analyzed. The free stream velocity, the injection velocity at the surface, moving velocity of the wedge surface, the wedge angle and the power law index of non-Newtonian fluid are assumed variables. The fourth order Runge–Kutta method modified by Gill is used to solve the non-dimensional boundary layer equations for non-Newtonian flow field. Without fluid injection, for every angle of wedge β, a limiting value for velocity ratio λ _cr (velocity of the wedge surface/velocity of the uniform flow) is found for each power-law index n. The value of λ _cr increases with the increasing wedge angle β. The value of wedge angle also restricts the physical characteristics of the fluid to be used. The effects of the different parameters on velocity profile and on skin friction are studied and the drag reduction is discussed. In case of C = 2.5 and velocity ratio λ = 0.2 for wedge angle β = 0.5 with the fluid with power law-index n = 0.5, 48.8% drag reduction is obtained.     

Growth factors for Marangoni instability in a spherical liquid layer under zero-gravity conditions

The neutral-stability analysis presented by Hoefsloot et al. [3] is completed by computing the growth factorsβ for the normal modes and by showing that the neutral states (Re(β)=0) are stationary (Im(β)=0) rather than oscillatory (Im(β)≠0).     

Localized Thickening of a Compressed Elastic Band

An infinite elastic band is compressed along its unbounded direction, giving rise to a continuous family of homogeneous configurations that is parameterized by the compression rate β < 1 (β = 1 when there is no compression). It is assumed that, for some critical value β ₀, the compression force as a function of β has a strict local extremum and that the linearized equation around the corresponding homogeneous configuration is strongly elliptic. Under these conditions, there are nearby localized deformations that are asymptotically homogeneous. When the compression force reaches a strict local maximum at β ₀, they describe localized thickening and they occur for values of β slightly smaller than β ₀. Since the material is supposed to be hyperelastic, homogeneous and isotropic, the localized deformations are not due to localized imperfections. The method follows the one developed by A. Mielke for an elastic band under traction: interpretation of the nonlinear elliptic system as an infinite dimensional dynamical system in which the unbounded direction plays the role of time, its reduction to a center manifold and the existence of a homoclinic solution to the reduced finite dimensional problem in [A. Mielke, Hamiltonian and Lagrangian fiows on center manifolds, Lecture Notes in Mathematics 1489. Springer, Berlin Heidelberg New York, 1991]. The main difference lies in the fact that Agmon's condition does not hold anymore and therefore the linearized problem cannot be analyzed as in Mielke's work.     

Slip coefficients and macroparameter jumps of a diatomic gas with rotational degrees of freedom on a weakly curved spherical surface

Using the half-space moment method, the problem of the slip of a diatomic gas along a rigid spherical surface is solved within the framework of a model kinetic equation previously proposed which takes into account the rotational degrees of freedom of the gas. Second-order slip coefficients (correctionsC _m ^′ , β _R ^′ , and β _R to the isothermal and thermal slip which are linear with respect to the Knudsen number Kn) are obtained. The gas macroparameter jump coefficientsC _v andC _q, which are of the second order in the Knudsen number and characterize the discontinuity of the normal mass and heat fluxes on the gas-rigid phase interface, are calculated. These coefficients are given as functions of the tangential momentum accommodation coefficient, the translational and rotational energy accommodation coefficients, and the Prandtl number. The coefficients are calculated for certain diatomic gases. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 163–173, January–February, 2000.     

MHD flow and heat transfer for the upper-convected Maxwell fluid over a stretching sheet

In the present work, the effect of MHD flow and heat transfer within a boundary layer flow on an upper-convected Maxwell (UCM) fluid over a stretching sheet is examined. The governing boundary layer equations of motion and heat transfer are non-dimensionalized using suitable similarity variables and the resulting transformed, ordinary differential equations are then solved numerically by shooting technique with fourth order Runge–Kutta method. For a UCM fluid, a thinning of the boundary layer and a drop in wall skin friction coefficient is predicted to occur for higher the elastic number. The objective of the present work is to investigate the effect of Maxwell parameter β, magnetic parameter Mn and Prandtl number Pr on the temperature field above the sheet.     

Further Comments on “Combined Forced and Free Convective Flow in a Vertical Porous Channel: The Effects of Viscous Dissipation and Pressure Work”

This note is concerned with the assertion of Barletta and Nield (2009a) that “a fluid with a thermal expansion coefficient greater than that of a perfect gas (β > β _perfect gas) is of marginal or no interest in the framework of convection in porous media”, and that for a remark of Magyari (Transp. Porous Media, 2009) about the forced convection eigenflow solutions, the circumstance β > β _perfect gas does not represent “a sound physical basis”. Here, it is shown, however, that these assertions are in contradiction with the experimentally measured values of β for important technical fluids as e.g., air, nitrogen, carbon dioxide, and ammonia where, in the temperature range between −20 and +100°C, just the inequality β > β _perfect gas holds.     

Existence of Singular Self-Similar Solutions of the Three-Dimensional Euler Equations in a Bounded Domain

A self-similar solution of the three-dimensional (3d) incompressible Euler equations is defined byu(x,t) =U(y)/t*-t) ^{α, y = x/(t* ~ t)β,α,β>} 0, whereU(y) satisfiesζU + βy. ΔU + U. VU + VP = 0,divU = 0. For α = β = 1/2, which is the limiting case of Leray’s self-similar Navier—Stokes equations, we prove the existence of(U,P) ε H³(Ω,R³ X R) in a smooth bounded domain Ω, with the inflow boundary data of non-zero vorticity. This implies the possibility that solutions of the Euler equations blow up at a timet = t*, t* < +∞.     

Asymptotic Behavior of Structures Made of Straight Rods

This paper is devoted to describe the deformations and the elastic energy for structures made of straight rods of thickness 2δ when δ tends to 0. This analysis relies on the decomposition of the large deformation of a single rod introduced in Blanchard and Griso (Anal. Appl. 7(1):21–71, 2009) and on the extension of this technique to a multi-structure. We characterize the asymptotic behavior of the infimum of the total elastic energy as the minimum of a limit functional for an energy of order δ ^β (2<β≤4).     

Miscible Thermo-Viscous Fingering Instability in Porous Media. Part 2: Numerical Simulations

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.     

Secondary flow patterns and mixing in laminar pulsating flow through a curved pipe

Mixing by secondary flow is studied by particle image velocimetry (PIV) in a developing laminar pulsating flow through a circular curved pipe. The pipe curvature ratio is η = r ₀/r _c  = 0.09, and the curvature angle is 90°. Different secondary flow patterns are formed during an oscillation period due to competition among the centrifugal, inertial, and viscous forces. These different secondary-flow structures lead to different transverse-mixing schemes in the flow. Here, transverse mixing enhancement is investigated by imposing different pulsating conditions (Dean number, velocity ratio, and frequency parameter); favorable pulsating conditions for mixing are introduced. To obviate light-refraction effects during PIV measurements, a T-shaped structure is installed downstream of the curved pipe. Experiments are carried out for the Reynolds numbers range 420 ≤ Re_st ≤ 1,000 (Dean numbers 126.6 ≤ Dn ≤ 301.5) corresponding to non-oscillating flow, velocity component ratios 1 ≤ (β = U _max,osc/U _m,st) ≤ 4 (the ratio of velocity amplitude of oscillations to the mean velocity without oscillations), and frequency parameters 8.37 < (α = r ₀(ω/ν)^0.5) < 24.5, where α² is the ratio of viscous diffusion time over the pipe radius to the characteristic oscillation time. The variations in cross-sectional average values of absolute axial vorticity (|ζ|) and transverse strain rate (|ε|) are analyzed in order to quantify mixing. The effects of each parameter (Re_st, β, and α) on transverse mixing are discussed by comparing the dimensionless vorticities (|ζ _P |/|ζ _S |) and dimensionless transverse strain rates (|ε _P |/|ε _S |) during a complete oscillation period.     

Motion of a sphere in an electrically conducting rotating fluid

The combined effect of rotation and magnetic field is investigated for the axisymmetric flow due to the motion of a sphere in an inviscid, incompressible electrically conducting fluid having uniform rotation far upstream. The steady-state linearized equations contain a single parameter α=1/2βR _m, β being the magnetic pressure number and R _m the magnetic Reynolds number. The complete solution for the flow field and magnetic field is obtained and the distribution of vorticity and current density is found. The induced vorticity is O(α⁴) and the current density is O(R _m) on the sphere.     

Experimental Study of Seepage Properties of Broken Sandstone Under Different Porosities

In coal mining the water flow in broken rock is a very common phenomenon. Study of seepage properties of broken rock is one of the basic subjects required in order to understand the stability of rock surrounding roadways, preventing disasters such as water inrush and gas outbursts and developing underground resources. So far, quantitative studies on the nonlinear seepage properties of broken sandstone under different porosities are not extensive in the research literature. In this article, by means of an electro-hydraulic servo-controlled test system (MTS815.02) and a patent seepage device, the seepage properties under different conditions of porosity were tested on broken sandstone of five different grain sizes. Based on the loading method of controlling the axial compression displacement and steady permeating method, we obtained curves of the relation of pore pressure with time, as well as the relation curves between the pore pressure gradient for steady seepage and velocity. Furthermore, we calculated the permeability k and non-Darcy coefficient β corresponding to different porosities by fitting these curves with the binomial expression. This study indicates that: (1) the seepage properties of broken sandstone are closely related to grain size, load levels, and porosity structure; (2) the permeability k decreases, while the coefficient β increases with a decrease in porosity φ, but both the k − φ and the β − φ curves show some local fluctuations; (3) the permeability k of the broken sandstone has a magnitude of 10⁻¹⁴–10⁻¹² m², while the coefficient β ranges from 10¹⁰ to 10¹² m⁻¹. The results obtained provide some information for further study of the nonlinear seepage behavior of broken rock theoretically.     

A Hierarchy of Plate Models Derived from Nonlinear Elasticity by Gamma-Convergence

We derive a hierarchy of plate models from three-dimensional nonlinear elasticity by Γ-convergence. What distinguishes the different limit models is the scaling of the elastic energy per unit volume ∼h^β, where h is the thickness of the plate. This is in turn related to the strength of the applied force ∼h^α. Membrane theory, derived earlier by Le Dret and Raoult, corresponds to α=β=0, nonlinear bending theory to α=β=2, von Kármán theory to α=3, β=4 and linearized vK theory to α>3. Intermediate values of α lead to certain theories with constraints. A key ingredient in the proof is a generalization to higher derivatives of our rigidity result [29] which states that for maps v:(0,1)³→ℝ³, the L² distance of ∇v from a single rotation is bounded by a multiple of the L² distance from the set SO(3) of all rotations.     

Momentum and heat transfer from an asymmetrically confined circular cylinder in a plane channel

Unsteady momentum and heat transfer from an asymmetrically confined circular cylinder in a plane channel is numerically investigated using FLUENT for the ranges of Reynolds numbers as 10≤Re≤500, of the blockage ratio as 0.1≤β≤0.4, and of the gap ratio as 0.125≤γ≤1 for a constant value of the Prandtl number of 0.744. The transition of the flow from steady to unsteady (characterized by critical Re) is determined as a function of γ and β. The effect of γ on the mean drag and lift coefficients, Strouhal number (St), and Nusselt number (Nu _w ) is studied. Critical Re was found to increase with decreasing γ for all values of β. and St were found to increase with decreasing values of γ for fixed β and Re. The effect of decrease in γ on was found to be negligible for all blockage ratios investigated.     

Stress concentration around a small spherical inhomogeneous inclusion on the axis of a circular cylinder in torsion

Stresses are calculated for a small elastic inclusion bonded (in a way to ensure continuous displacements and tractions across the surface of the inclusion) to a large circular cylindrical shaft. The inclusion considered is inhomogeneous and of spherical shape. The shear modulus of the inclusion is assumed to be μ₀ exp(−1/2β² r ²). Some special cases are also discussed.     

Postbuckling of elastic beam subjected to a concentrated moment within the span length of beam

This paper aims to present the exact closed form solutions and postbuckling behavior of the beam under a concentrated moment within the span length of beam. Two approaches are used in this paper. The nonlinear governing differential equations based on elastica theory are derived and solved analytically for the exact closed form solutions in terms of elliptic integral of the first and second kinds. The results are presented in graphical diagram of equilibrium paths, equilibrium configurations and critical loads. For validation of the results from the first approach, the shooting method is employed to solve a set of nonlinear differential equations with boundary conditions. The set of nonlinear governing differential equations are integrated by using Runge–Kutta method fifth order with adaptive step size scheme. The error norms of the end conditions are minimized within prescribed tolerance (10⁻⁵). The results from both approaches are in good agreement. From the results, it is found that the stability of this type of beam exhibits both stable and unstable configurations. The limit load point existed. The roller support can move through the hinged support in some cases of β and leads to the more complex of the configuration shapes of the beam.     

Effects of drift angle on model ship flow

The effects of drift angle on model ship flow are investigated through towing tank tests for the Series 60 C_B=0.6 cargo/container model ship. Resistance, side force, drift moment, sinkage, trim, and heel data are procured for a range of drift angles β and Froude numbers (Fr) and the model free condition. Detailed free-surface and mean velocity and pressure flow maps are procured for high and low Fr=0.316 and 0.16 and β=5° and 10° (free surface) and β=10° (mean velocity and pressure) for the model fixed condition (i.e. fixed with zero sinkage, trim, and heel). Comparison of results at high and low Fr and previous data for β=0° enables identification of important free-surface and drift effects. Geometry, conditions, data, and uncertainty analysis are documented in sufficient detail so as to be useful as a benchmark for computational fluid dynamics (CFD) validation. The resistance increases linearly with β with same slope for all Fr, whereas the increases in the side force, drift moment, sinkage, trim, and heel with β are quadratic. The wave profile is only affected near the bow, i.e. the bow wave amplitude increases/decreases on the windward/leeward sides, whereas the wave elevations are affected throughout the entire wave field. However, the wave envelope angle on both sides is nearly the same as β=0°, i.e. the near-field wave pattern rotates with the hull and remains within a similar wave envelope as β=0°. The wave amplitudes are significantly increased/decreased on the windward/leeward sides. The wake region is also asymmetric with larger wedge angle on the leeward side. The boundary layer and wake are dominated by the hull vortex system consisting of fore body keel, bilge, and wave-breaking vortices and after body bilge and counter-rotating vortices. The occurrence of a wave-breaking vortex for breaking bow waves has not been previously documented in the literature. The trends for the maximum vorticity, circulation, minimum axial velocity, and trajectories are discussed for each vortex. Received: 16 September 1999/Accepted: 8 November 2001