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.     

Heat transfer characteristics from in-line arrays of free impinging jets

An experimental investigation of the convective heat transfer on a flat surface in a multiple-jet system is described. A thin metal sheet was heated electrically and cooled from one side. On the other black coated side the temperature field was measured using an IR camera. Varied parameters were the jet Reynolds number in the range from 1,400 to 41,400, the normalized distance nozzle to sheet H/d from 1 to 10, and the normalized nozzle spacing S/d from 2 to 10. A geometrical arrangement of nine nozzle in-line arrays was tested. The results show that the multiple-jet system enhances the local and average heat transfer in comparison with that of a single nozzle. A maximum of the heat transfer was found for the normalized spacing S/d = 6.0. The normalized distance H/d has nearly no effect on the heat transfer in the range 2 ≤ H/d ≤ 4. The maximum average Nusselt number was correlated as a function of the jet Reynolds number      

Numerical Study of Natural Convection Heat Transfer in an Inclined Porous Cavity with Time-Periodic Boundary Conditions

Kalabin et al. (Numer. Heat Transfer A 47, 621-631, 2005) studied the unsteady natural convection for the sinusoidal oscillating wall temperature on one side wall and constant average temperature on the opposing side wall. The present article is on the unsteady natural convective heat transfer in an inclined porous cavity with similar temperature boundary conditions as those of Kalabin et al. The inclined angle of the cavity is varied from 0° to 80°. The flow field is modeled with the Brinkman-extended Darcy model. The combined effects of inclination angle of the enclosure and oscillation frequency of wall temperature are studied for Ra* = 10³, Da = 10⁻³, , and Pr=1. Some results are also obtained with the Darcy–Brinkman–Forchheimer model and Darcy’s law and are compared with the present Brinkman-extended Darcy model. The maximal heat transfer rate is attained at the oscillating frequency f = 46.7π and the inclined angle .     

The Asymptotic Finite-dimensional Character of a Spectrally-hyperviscous Model of 3D Turbulent Flow

We obtain attractor and inertial-manifold results for a class of 3D turbulent flow models on a periodic spatial domain in which hyperviscous terms are added spectrally to the standard incompressible Navier–Stokes equations (NSE). Let P _m be the projection onto the first m eigenspaces of A =−Δ, let μ and α be positive constants with α ≥3/2, and let Q _m =I − P _m , then we add to the NSE operators μ A _φ in a general family such that A _φ≥Q _m A ^α in the sense of quadratic forms. The models are motivated by characteristics of spectral eddy-viscosity (SEV) and spectral vanishing viscosity (SVV) models. A distinguished class of our models adds extra hyperviscosity terms only to high wavenumbers past a cutoff λ _m0 where m ₀ ≤ m, so that for large enough m ₀ the inertial-range wavenumbers see only standard NSE viscosity. We first obtain estimates on the Hausdorff and fractal dimensions of the attractor (respectively and ). For a constant K _α on the order of unity we show if μ ≥ ν that and if μ ≤ ν that where ν is the standard viscosity coefficient, l ₀ = λ₁^−1/2 represents characteristic macroscopic length, and is the Kolmogorov length scale, i.e. where is Kolmogorov’s mean rate of dissipation of energy in turbulent flow. All bracketed constants and K _α are dimensionless and scale-invariant. The estimate grows in m due to the term λ _m /λ₁ but at a rate lower than m ^3/5, and the estimate grows in μ as the relative size of ν to μ. The exponent on is significantly less than the Landau–Lifschitz predicted value of 3. If we impose the condition , the estimates become for μ ≥ ν and for μ ≤ ν. This result holds independently of α, with K _α and c _α independent of m. In an SVV example μ ≥ ν, and for μ ≤ ν aspects of SEV theory and observation suggest setting for 1/c within α orders of magnitude of unity, giving the estimate where c _α is within an order of magnitude of unity. These choices give straight-up or nearly straight-up agreement with the Landau–Lifschitz predictions for the number of degrees of freedom in 3D turbulent flow with m so large that (e.g. in the distinguished-class case for m ₀ large enough) we would expect our solutions to be very good if not virtually indistinguishable approximants to standard NSE solutions. We would expect lower choices of λ _m (e.g. with a > 1) to still give good NSE approximation with lower powers on l ₀/l _ε, showing the potential of the model to reduce the number of degrees of freedom needed in practical simulations. For the choice , motivated by the Chapman–Enskog expansion in the case m = 0, the condition becomes , giving agreement with Landau–Lifschitz for smaller values of λ _m then as above but still large enough to suggest good NSE approximation. Our final results establish the existence of a inertial manifold for reasonably wide classes of the above models using the Foias/Sell/Temam theory. The first of these results obtains such an of dimension N > m for the general class of operators A _φ if α > 5/2. The special class of A _φ such that P _m A _φ = 0 and Q _m A _φ ≥ Q _m A ^α has a unique spectral-gap property which we can use whenever α ≥ 3/2 to show that we have an inertial manifold of dimension m if m is large enough. As a corollary, for most of the cases of the operators A _φ in the distinguished-class case that we expect will be typically used in practice we also obtain an , now of dimension m ₀ for m ₀ large enough, though under conditions requiring generally larger m ₀ than the m in the special class. In both cases, for large enough m (respectively m ₀), we have an inertial manifold for a system in which the inertial range essentially behaves according to standard NSE physics, and in particular trajectories on are controlled by essentially NSE dynamics.      

Unconditional nonlinear exponential stability of the motionless conduction-diffusion solution

Nonlinear stability of the motionless state of a heterogeneous fluid with constant temperature-gradient and concentration-gradient is studied for both cases of stress-free and rigid boundary conditions. By introducing new energy functionals we have shown that for τ=P _C /P _T ≤1, the motionless state is always stable and for τ≤1, the sufficient and necessary conditions for stability coincide, whereP _C ,P _T ,C andR are the Schmidt number, Prandtl number, Rayleigh number for solute and heat respectively. Moreover, the criteria guarantees the exponential stability. The Project supported by SRF for ROCS, China Postdoctoral Science Foundation and the National Basic Research Project “Nonlinear Science”     

Forced convection gaseous slip flow in circular porous micro-channels

Laminar forced convection of gaseous slip flow in a circular micro-channel filled with porous media under local thermal equilibrium condition is studied numerically using the finite difference technique. Hydrodynamically fully developed flow is considered and the Darcy–Brinkman–Forchheimer model is used to model the flow inside the porous domain. The present study reports the effect of several operating parameters (Knudsen number (Kn), Darcy number (Da), Forchhiemer number (Γ), and modified Reynolds number ) on the velocity slip and temperature jump at the wall. Results are given in terms of the velocity distribution, temperature distribution, skin friction , and the Nusselt number (Nu). It is found that the skin friction is increased by (1) decreasing Knudsen number, (2) increasing Darcy number, and (3) decreasing Forchheimer number. Heat transfer is found to (1) decrease as the Knudsen number, or Forchheimer number increase, (2) increase as the Peclet number or Darcy number increase.     

Roughness induced forced convective laminar-transitional micropipe flow: energy and exergy analysis

Variable fluid property continuity, Navier–Stokes and energy equations are solved for roughness induced forced convective laminar-transitional flow in a micropipe. Influences of Reynolds number, heat flux and surface roughness, on the momentum-energy transport mechanisms and second-law of thermodynamics, are investigated for the ranges of Re = 1–2,000, Q = 5–100 W/m² and ε = 1–50 μm. Numerical investigations put forward that surface roughness accelerates transition with flatter velocity profiles and increased intermittency values (γ); such that a high roughness of ε = 50 μm resulted in transitional character at Re _tra = 450 with γ = 0.136. Normalized friction coefficient (C _f^*) values showed augmentation with Re, as the evaluated C _f^* are 1.006, 1.028 and 1.088 for Re = 100, 500 and 1,500, respectively, at ε = 1 μm, the corresponding values rise to C _f^* = 1.021, 1.116 and 1.350 at ε = 50 μm. Heat transfer rates are also recorded to rise with Re and ε; moreover the growing influence of ε on Nusselt number with Re is determined by the Nu _ε=50 μm/Nu _ε=1 μm ratios of 1.086, 1.168 and 1.259 at Re = 500, 1,000 and 1,500. Thermal volumetric entropy generation values decrease with Re and ε in heating; however the contrary is recorded for frictional volumetric entropy generation data, where the augmentations in are more considerable when compared with the decrease rates of      

Experimental investigations of controlled transition of a laminar separation bubble in an axisymmetric diffuser

The instability of a pressure-induced laminar separation bubble is examined experimentally on an axisymmetric diffuser for a Reynolds number range 7,800 ≤ ≤ 11,400 for an inlet pipe diameter D ₁ (50 mm) and as mean input flow velocity 4.2 m/s ≤ u _m ≤ 6.1 m/s. A characterization of the base flow shows a wide-spread separation at the smooth diverging contour which gives rise to a massive amplification of instabilities. Controlled disturbances are introduced by means of a slot and a membrane actuator to trigger the transition, and the receptivity of the perturbations to the laminar boundary layer is evaluated. Different axisymmetric and azimuthal disturbances are applied in order to study their influence on the laminar–turbulent transition. The measurements show a clear dependence of the transition scenario and the reattachment length on the actuation mode.     

Transition layer thickness in a fluid-porous medium of multi-sized spherical beads

Momentum and mass transfer at fluid–porous interfaces occur in many technical and natural applications. The vertical extend below a fluid–porous interface within which the free fluid velocity reduces to a constant Darcy velocity in the porous medium is known as Brinkman layer. Recently, the Brinkman layer thickness (δ) has been measured for a porous bed of mono-sized spherical beads, and was found to be in the order of the particle diameter (d). In this study, we investigate a porous medium made of multi-sized spherical beads. The measured averaged interfacial velocity field clearly indicated that, in the case of multi-sized beads, δ is in the order of a characteristic diameter given by with x _i and d _i being the weight fraction and diameter of the component i in the mixture.     

Variable permeability effect on convection in binary mixtures saturating a porous layer

The Darcy Model with the Boussinesq approximation is used to study natural convection in a shallow porous layer, with variable permeability, filled with a binary fluid. The permeability of the medium is assumed to vary exponentially with the depth of the layer. The two horizontal walls of the cavity are subject to constant fluxes of heat and solute while the two vertical ones are impermeable and adiabatic. The governing parameters for the problem are the thermal Rayleigh number, R _T, the Lewis number, Le, the buoyancy ratio, φ, the aspect ratio of the cavity, A, the normalized porosity, ε, the variable permeability constant, c, and parameter a defining double-diffusive convection (a = 0) or Soret induced convection (a = 1). For convection in an infinite layer, an analytical solution of the steady form of the governing equations is obtained on the basis of the parallel flow approximation. The onset of supercritical convection, or subcritical, convection are predicted by the present theory. A linear stability analysis of the parallel flow model is conducted and the critical Rayleigh number for the onset of Hopf’s bifurcation is predicted numerically. Numerical solutions of the full governing equations are found to be in excellent agreement with the analytical predictions.     

A Variational Approach to the Fredholm Alternative for the p-Laplacian near the First Eigenvalue

We are concerned with the existence of a weak solution to the degenerate quasi-linear Dirichlet boundary value problem

Thermal characteristics in the annulus with an inner rotating rib-roughness cylinder

This work experimentally investigates the heat transfer characteristics in the annulus with an inner rotating rib-roughness cylinder, whose flow and thermal behaviors are associated with Taylor number (Ta) and centrifugal buoyancy parameter (Gr _Ω/Ta). The operating range of Ta is from 4.90 × 10² to 5.80 × 10⁵, while the surface of the inner cylinder is heated up with several constant heat fluxes (279, 425 and 597 W/m²) to obtain various values of Gr _Ω/Ta. Besides, three modes of the inner cylinder without/with longitudinal ribs are considered. The end of the annular channel is connected to a side chamber to fit practical applications (such as in the rotary blade coupling of a four-wheel-drive vehicle). The experimental results show that the average Nusselt number was almost constant at low Ta, but increased rapidly with Ta when Ta exceeded some critical value (3,000–5,200 for present study). Additionally, the Gr _Ω/Ta effect on the heat transfer was negligible herein. Furthermore, by comparing with the inner cylinder without longitudinal ribs, stalling ribs on the inner cylinder increases the transport of heat by a factor of 1.22 at 10⁵ < Ta < 10⁶, and embedding cavities into the ribs increases the transport of heat by a factor of 1.16 at 10⁵ < Ta < 10⁶. Finally, the relationships between the and the Ta for various modes of test sections were proposed.     

Effects of Large Scales Motion in Unstable Stratified Shear Flows

Homogeneous turbulence under unstable uniform stratification (N ² < 0) and vertical shear is investigated by using the linear theory (or the so-called rapid distortion theory, RDT) for an initial isotropic turbulence over a range −∞ ≤ R _i =N ²/S ² ≤ 0. The initial potential energy is zero and P _r =1 (i.e. the molecular Prandtl number).One-dimensional (streamwise) k ₁−spectra, especially Θ₃₃(k ₁) (i.e., that of the vertical kinetic energy, are investigated. In agreement with previous experiments, it is found that the unstable stratification affects the turbulence quantities at all scales. A significant increase of the vertical kinetic energy is observed at low wavenumbers k ₁ (i.e. large scales) due to an increase of the stratification . The effect of the shear (S) is appreciable only at high wavenumbers k ₁ (i.e. small scales).Based on the importance of the spectral components with k ₁ = 0, the asymptotic forms of Θ _ij (k ₁ = 0) or equivalently the so-called “two-dimensional” energy components (2DEC) are analyzed in detail. The asymptotic form for the ratio of 2DEC is compared to the long-time limit of the ratio of real energies. In the unstable stratified shearless case (S=0,N ² ≠ 0) the variances and the covariances of the velocity and the density are derived analytically in terms of the Weber functions, while when S ≠ 0 and N ² ≠ 0 they are obtained numerically (−100 ≤ R _i <0 and . The results are discussed in connection to previous experimental results in unstable stratified open channel flows cooled from above by Komori et al. Phy Fluids 25, 1539–1546 (1982).It is shown that the Richardson number dependence of the long-time limit of the ratios of real energies is well described by this “simple” model (i.e. the dependence of the long-time limit of 2DEC on R _i ). For example, the ratio of the potential energy to the kinetic energy (q ²/2), approaches −R _i /(1−R _i ), the ratio of turbulent energy production by buoyancy forces to production by shearing forces (i.e. the flux Richardson number, R _if ), approaches R _i . Also, the Richardson number dependence of the principal angle (β) of the Reynolds stress tensor and the angle (β_ρ) of the scalar flux vector is fairly predicted by this model .On the other hand, it is found that the above ratios are insensitive to viscosity, while the ratios ɛ /q ² and , depend on the viscosity and they evolve asymptotically like t ⁻¹. The turbulent Froude number, F _rt =(L _Oz /L _E )^2/3, where L _Oz and L _E are the Ozmidov length scale and the Ellison length scale, respectively, evolves asymptotically like t ^−1/3.     

Dynamic Contact Behavior of a Golf Ball during an Oblique Impact

The oblique impact between a golf ball and a rigid steel target was studied using a high-speed video camera. Video images recorded before and after the impact were used to determine the inbound velocity v _i, rebound velocity v _r, inbound angle θ_i, rebound angle θ_r, and the coefficient of restitution e. The results showed that θ_r and e decreased as v _i increased. The maximum compression ratio η_c, contact time t _c, average angular velocity , and tangential velocity , along the target were determined from images obtained during the impact. The images demonstrated that η_c increased with v _i while t _c decreased. In addition, and increased almost linearly as v _i increased. A rigid body model was used to estimate the final angular velocity ω* and tangential velocity ν_t* at the end of the impact; these results were then compared with experimental data.     

A study of the upper limit of solid scatters density for gray Lattice Boltzmann Method

The upper limit of the solid scatters density ns （x）, a key parameter for the simulation of flows in porous media with a gray Lattice Boltzmann Method, is studied by an analytical way for the infiltration Poiseuille flow between two infinite parallel plates. Analyses of three different gray Lattice Boltzmann schemes, separately proposed by Gao and Sharma et al., Dardis and McCloskey, and Thorne and Sukop, indicate that the effective domain of Gao and Sharma＇s scheme is restricted to ns 〈 1/2√3≈0.289, Dardis and McCloskey＇s scheme is restricted to ns 〈 （√57-1）/28≈0.234, and that there is no extra restriction on ns（x） with Thorne and Sukop＇s scheme. These results are obtained for the dimensionless relaxation time τ= 1. The above analytical results are verified by our numerical simulations. The use of a gray LBM is further illustrated by simulating the flow at the interface of a porous medium. Simulation results yield velocity profiles which agree very well with Brinkman＇s prediction.     

A New Approach to the Fundamental Theorem of Surface Theory

The fundamental theorem of surface theory classically asserts that, if a field of positive-definite symmetric matrices (a _αβ ) of order two and a field of symmetric matrices (b _αβ ) of order two together satisfy the Gauss and Codazzi-Mainardi equations in a simply connected open subset ω of , then there exists an immersion such that these fields are the first and second fundamental forms of the surface , and this surface is unique up to proper isometries in . The main purpose of this paper is to identify new compatibility conditions, expressed again in terms of the functions a _αβ and b _αβ , that likewise lead to a similar existence and uniqueness theorem. These conditions take the form of the matrix equation

Forced convection heat transfer and friction characteristics in multi-exit -rectangular package with inclined tube bundles

An experiment was carried out to investigate the characteristics of the heat transfer and pressure drop for forced convection airflow over tube bundles that are inclined relative to the on-coming flow in a rectangular package with one outlet and two inlets. The experiments included a wide range of angles of attack and were extended over a Reynolds number range from about 250 to 12,500. Correlations for the Nusselt number and pressure drop factor are reported and discussed. As a result, it was found that at a fixed Re, for the tube bundles with attack angle of 45 ° has the best heat transfer coefficient, followed by 60, 75 and 90 °, respectively. This investigation also introduces the factors which can be used for finding the heat transfer and the pressure drop factor on the tube bundles positioned at different angles to the flow direction. Moreover, no perceptible dependence of Cand C on Re was detected. In addition, flow visualizations were explored to broaden our fundamental understanding of the heat transfer for the present study.     

Non-autonomous 2D Navier–Stokes System with Singularly Oscillating External Force and its Global Attractor

We study the global attractor of the non-autonomous 2D Navier–Stokes (N.–S.) system with singularly oscillating external force of the form . If the functions g ₀(x, t) and g ₁ (z, t) are translation bounded in the corresponding spaces, then it is known that the global attractor is bounded in the space H, however, its norm may be unbounded as since the magnitude of the external force is growing. Assuming that the function g ₁ (z, t) has a divergence representation of the form where the functions (see Section 3), we prove that the global attractors of the N.–S. equations are uniformly bounded with respect to for all . We also consider the “limiting” 2D N.–S. system with external force g ₀(x, t). We have found an estimate for the deviation of a solution of the original N.–S. system from a solution u ⁰(x, t) of the “limiting” N.–S. system with the same initial data. If the function g ₁ (z, t) admits the divergence representation, the functions g ₀(x, t) and g ₁ (z, t) are translation compact in the corresponding spaces, and , then we prove that the global attractors converges to the global attractor of the “limiting” system as in the norm of H. In the last section, we present an estimate for the Hausdorff deviation of from of the form: in the case, when the global attractor is exponential (the Grashof number of the “limiting” 2D N.–S. system is small).      

Pore-Scale Analysis of NAPL Blob Dissolution and Mobilization in Porous Media

A pore-scale analysis of nonaqueous phase liquid (NAPL) blob dissolution and mobilization in porous media was presented. Dissolution kinetics of residual NAPLs in an otherwise water-saturated porous medium was investigated by conducting micromodel experiments. Changes in residual NAPL volume were measured from recorded video images to calculate the mass transfer coefficient, K and the lumped mass transfer rate coefficient, k. The morphological characteristics of the blobs such as specific and intrinsic area were found to be independent of water flow rate except at NAPL saturations below 2%. Dissolution process was also investigated by separating the mass transfer into zones of mobile and immobile water. The fractions of total residual NAPL perimeters in contact with mobile water and immobile water were measured and their relationship to the mass transfer coefficient was discussed. In general, residual NAPLs are removed by dissolution and mobilization. Although these two mechanisms were studied individually by others, their simultaneous occurrence was not considered. Therefore, in this study, mobilization of dissolving NAPL blobs was investigated by an analysis of the forces acting on a trapped NAPL blob. A dimensional analysis was performed to quantify the residual blob mobilization in terms of dimensionless Capillary number (Ca _I). If Ca _I is equal to or greater than the trapping number defined as , then blob mobilization is expected.     

The Toda System and Clustering Interfaces in the Allen–Cahn equation

We consider the Allen–Cahn equation in a bounded, smooth domain Ω in , under zero Neumann boundary conditions, where is a small parameter. Let Γ₀ be a segment contained in Ω, connecting orthogonally the boundary. Under certain nondegeneracy and nonminimality assumptions for Γ₀, satisfied for instance by the short axis in an ellipse, we construct, for any given N ≥ 1, a solution exhibiting N transition layers whose mutual distances are and which collapse onto Γ₀ as . Asymptotic location of these interfaces is governed by a Toda-type system and yields in the limit broken lines with an angle at a common height and at main order cutting orthogonally the boundary.