Soret and Dufour effects on Double-Diffusive Free Convection From a Corrugated Vertical Surface in a Non-Darcy Porous Medium

Combined heat and mass transfer process by natural convection from a wavy vertical surface immersed in a fluid-saturated semi-infinite porous medium due to Soret and Dufour effects for Forchheimer extended non-Darcy model has been analyzed. A similarity transformation followed by a wavy to flat surface transformation is applied to the governing coupled non-linear partial differential equations, and they are reduced to boundary layer equations. The obtained boundary layer equations are solved by finite difference scheme based on the Keller-Box approach in conjunction with block-tridiagonal solver. Detailed simulations are carried out for a wide range of parameters like Groshof number (Gr*), Lewis number (Le), Buoyancy ratio (B), Wavy wall amplitude (a), Soret number (S _r ), and Dufour number (D _f ). Comparison tables local and average Nusselt (Nu) number, local and average Sherwood (Sh) number plots are presented.     

Spatial characteristics of dynamic properties of modulated wavy vortex flow in a rotating Couette system

Time-dependent velocity profiles V _z (z, t) have been measured in a cylindrical Couette system for a range of Reynolds numbers corresponding to wavy vortex flow (WVF) and modulated wavy vortex flow (MWV). The spatial distribution of the power spectrum was obtained and the behaviour of the dynamic properties of WVF and MWV was investigated with respect to Reynolds number (Re) at four corners of cells (inside-outside and Inflow-Outflow). All modes of f ₁, f ₂ and B were found to be equally strong and broad at different corners for lower Re number. At large Re, however, f ₂ remains broad only at the outside-Inflow corner, while B is strong at Outflow. It was also found that B disappears at the outside-Inflow corner for R ^* > 18. These results suggest that the dynamic properties of WVF and MWV would have a spatial dependence within a cell structure scale.     

Experimental study of vortex-induced vibrations of a cylinder near a rigid plane boundary in steady flow

In this study, the vortex-induced vibrations of a cylinder near a rigid plane boundary in a steady flow are studied experimentally. The phenomenon of vortex-induced vibrations of the cylinder near the rigid plane boundary is reproduced in the flume. The vortex shedding frequency and mode are also measured by the methods of hot film velocimeter and hydrogen bubbles. A parametric study is carded out to investigate the influences of reduced velocity, gap-to-diameter ratio, stability parameter and mass ratio on the amplitude and frequency responses of the cylinder. Experimental results indicate： （1） the Strouhal number （St） is around 0.2 for the stationary cylinder near a plane boundary in the sub-criti- cal flow regime; （2） with increasing gap-to-diameter ratio （eo/D）, the amplitude ratio （A/D） gets larger but frequency ratio （f/fn） has a slight variation for the case of larger values of eo/D（eo/D 〉 0.66 in this study）; （3） there is a clear difference of amplitude and frequency responses of the cylin- derbetween the larger gap-to-diameter ratios （e0/D 〉 0.66） and the smaller ones （e0/D 〈 0.3）; （4） the vibration of the cylinder is easier to occur and the range of vibration in terms of Vr number becomes more extensive with decrease of the stability parameter, but the frequency response is affected slightly by the stability parameter; （5） with decreasing mass ratio, the width of the lock-in ranges in terms of Vr and the frequency ratio （f/fn） become larger.     

The effect of suction or injection on the boundary layer flows induced by continuous surfaces stretched with prescribed skin friction

The boundary layer flow and heat transfer on a stretched surface moving with prescribed skin friction is studied for permeable surface. Three major cases are studied for isothermal surface (n=0) stretched corresponding to different dimensional skin friction boundary conditions namely; skin friction at the surface scales as (x ^?1/2) at m=0, constant skin friction at m=1/3 and skin friction scales as (x) at m=1. The constants m and n are the indices of the power law velocity and temperature exponent respectively. Similarity solutions are obtained for the boundary layer equations subject to power law temperature and velocity variation. The effect of various governing parameters, such as Prandtl number Pr, suction/injection parameter f _w , m and n are studied. The results show that for isothermal surface increasing m enhances the dimensionless heat transfer coefficient for fixed f _w at the suction case and the reverse is true at the injection case. Furthermore, for fixed m, as f _w increases the dimensionless heat transfer coefficient increases. Large enhancements are observed in the heat transfer coefficient as the temperature boundary condition along the surface changes from uniform to linear where the dimensional skin friction is of order (x) at m=1. This enhancement decreases as the suction increases.     

Mixed convection boundary layer flow over a horizontal circular cylinder with Newtonian heating

The steady mixed convection boundary layer flow over a horizontal circular cylinder, generated by Newtonian heating in which the heat transfer from the surface is proportional to the local surface temperature, is considered in this study. The governing boundary layer equations are first transformed into a system of non-dimensional equations via the non-dimensional variables, and then into non-similar equations before they are solved numerically using a numerical scheme known as the Keller-box method. Numerical solutions are obtained for the skin friction coefficient Re ^1/2 C _f and the local wall temperature θ _w (x) as well as the velocity and temperature profiles with two parameters, namely the mixed convection parameter λ and the Prandtl number Pr.     

A note on the heat transfer in convective fins

In this paper a generalized approach to the problem of heat transfer through convective fins is given. The proper dimensionless variables, which specify the general problem are identified, and upper bounds of the values of the dimensionless number N_r defined as “the ratio of the heat transferred by the fin to that of the corresponding bare surface” are derived. It was shwon that these limiting values of the N_r are 1/√B₁ and √2/B₁ for longitudinal fins and spines respectively, where B₁ is the Biot number hb/k, while for annular fins of constant thickness and hyperbolic profile, N_r? K(β)/√B_i, where K(β) is a number determined by the profile of the fin and the ratio β=x₂/x₁ of the outside to the inside radii. It was also shown that for longitudinal fins and spinces the possible adverse insulating effect by the use of the fin is avoided, if one selects the value of √hA/KC < 1, which is a rather stricter criterion than the one reported in the literature, namely that of hA/kC < 1 [2–5]. An example is given to show how one may utilize the appropriate value of N_r and the fin effectiveness e, to obtain the dimensions of the fin.     

Low-thrust Lambert transfer based on two-stage constant-vector thrust control method

The low-thrust Lambert transfer refers to that the spacecraft achieves the orbital transfer whose boundary conditions are represented by two sets of orbital elements at initial and final time by the low-thrust propulsion system. The modulus and direction of the low-thrust solutions in previous methods change with time, which leads to high control requirements for the engine. In this paper, to reduce the requirements of the engine, a practical two-stage constant-vector thrust control method is proposed, in which the magnitude and direction of the thrust are deemed as segmental constant value in TNH frame, where three components of the thrust are f_t, f_n, and f_h. First, the mathematical model of the two-stage constant-vector thrust is formulated, and a rapid algorithm is presented to obtain the solution based on the linearized sensitivity matrix, which describes the relationship between the constant-vector thrust and the change of the orbital elements approximately. Furthermore, two low-thrust Lambert strategies based on the two-stage constant-vector thrust are presented for cases of short-time transfer and long-time transfer. A sequence of numerical simulations demonstrated the efficiency of the proposed approaches. The proposed control strategies are solved rapidly, and they are also suitable for different types of orbits with J₂ perturbation, which are practical options for engineering applications.

     

Heat transfer and skin-friction in a turbulent boundary layer under a non-equilibrium longitudinal adverse pressure gradient

The experimental data on the effect of weak and moderate non-equilibrium adverse pressure gradients (APG) on the parameters of dynamic and thermal boundary layers are presented. The Reynolds number based on the momentum thickness at the beginning of the APG region was Re^** = 5500. The APG region was a slot channel with upper wall expansion angles from 0 to 14°. The profiles of the mean and fluctuation velocity components were measured using a single-component hot-wire anemometer. The friction coefficients were determined using two methods, namely, the indirect Clauser method and the direct method of weighting the lower wall region on a single-component strain-gage balance. The heat transfer coefficients were determined by a transient method using an IR camera. It is noticed that in the pressure gradient range realized the universal logarithmic region in the boundary layer profile is conserved. The values of the relative (divided by the parameters in zero gradient flow at the same value of Re^**) friction and heat transfer coefficients, together with the Reynolds analogy factor, are determined as functions of the longitudinal pressure gradient. The values of the relative friction coefficient reduced to c_f/c_f0 = 0.7 and those of the heat transfer to St/St₀ = 0.9. A maximum value of the Reynolds analogy factor (St/St₀)/(c_f/c_f0) = 1.16 was reached for the pressure gradient parameter β = 2.9.     

Reynolds' analogy for the thermal convection driven by nonisothermal stretching surfaces

In the present paper the steady boundary-layer flows induced by permeable stretching surfaces with variable temperature distribution are investigated under the aspect of Reynolds' analogy r = St _x /C _f(x). It is shown that for certain stretching velocities and wall temperature distributions, “Reynolds' function”r, i.e. the ratio of the local Stanton number St _x and the skin friction coefficient C _f(x) equals −1/2 for any value of the Prandtl number Pr and of the dimensionless suction/injection velocity f _w. In all of these cases, the dimensionless temperature field ϑ is connected to the dimensionless downstream velocity f ^′ by the simple relationship ϑ=(f ^′)^Pr. It is also shown that in the general case, Reynolds' function r may possess several singularities in f _w. The largest of them represents a critical value, so that for f _w<f _w,crit the solutions of the energy equation (although they still satisfy all the boundary conditions) become nonphysical.     

The asymptotic expansions of singularly perturbed boundary value problems

In this paper we study the singularly perturbed boundary value problem:εy″=f(t,y,ε),y(0)=ξ(ε),y(1)=η(ε).where εis a positive small parameter.In the conditions:f_(?)(0,y,0)≥m_0 ,f_(?)(l,y,0)≥m_0 and f_(?)f(t,y,ε)≥0 ,we prove the existences,and uniformly valid asymptotic expansions of solutions for the given boundary value problems,and hence we improve the existing results.     

Computational modeling of biomagnetic micropolar blood flow and heat transfer in a two-dimensional non-Darcian porous medium

We study theoretically and computationally the incompressible, non-conducting, micropolar, biomagnetic (blood) flow and heat transfer through a two-dimensional square porous medium in an (x,y) coordinate system, bound by impermeable walls. The magnetic field acting on the fluid is generated by an electrical current flowing normal to the x–y plane, at a distance l beneath the base side of the square. The flow regime is affected by the magnetization B ₀ and a linear relation is used to define the relationship between magnetization and magnetic field intensity. The steady governing equations for x-direction translational (linear) momentum, y-direction translational (linear) momentum, angular momentum (micro-rotation) and energy (heat) conservation are presented. The energy equation incorporates a special term designating the thermal power per unit volume due to the magnetocaloric effect. The governing equations are non-dimensionalized into a dimensionless (ξ,η) coordinate system using a set of similarity transformations. The resulting two point boundary value problem is shown to be represented by five dependent non-dimensional variables, f _ξ  (velocity), f _η (velocity), g (micro-rotation), E (magnetic field intensity) and θ (temperature) with appropriate boundary conditions at the walls. The thermophysical parameters controlling the flow are the micropolar parameter (R), biomagnetic parameter (N _H ), Darcy number (Da), Forchheimer (Fs), magnetic field strength parameter (Mn), Eckert number (Ec) and Prandtl number (Pr). Numerical solutions are obtained using the finite element method and also the finite difference method for Ec=2.476×10⁻⁶ and Prandtl number Pr=20, which represent realistic biomagnetic hemodynamic and heat transfer scenarios. Temperatures are shown to be considerably increased with Mn values but depressed by a rise in biomagnetic parameter (N _H ) and also a rise in micropolarity (R). Translational velocity components are found to decrease substantially with micropolarity (R), a trend consistent with Newtonian blood flows. Micro-rotation values are shown to increase considerably with a rise in R values but are reduced with a rise in biomagnetic parameter (N _H ). Both translational velocities are boosted with a rise in Darcy number as is micro-rotation. Forchheimer number is also shown to decrease translational velocities but increase micro-rotation. Excellent agreement is demonstrated between both numerical solutions. The mathematical model finds applications in blood flow control devices, hemodynamics in porous biomaterials and also biomagnetic flows in highly perfused skeletal tissue. Dedicated to Professor Y.C. Fung (1919-), Emeritus Professor of Biomechanics, Bioengineering Department, University of California at San Diego, USA for his seminal contributions to biomechanics and physiological fluid mechanics over four decades and his excellent encouragement to the authors, in particular OAB, with computational biofluid dynamics research.     

Radiative and Thermal Dispersion Effects on Non-Darcy Natural Convection with Lateral Mass Flux for Non-Newtonian Fluid from a Vertical Flat Plate in a Saturated Porous Medium

The problem of natural convective heat transfer for a non-Newtonian fluid from an impermeable vertical plate embedded in a fluid-saturated porous medium has been analyzed. Non-Darcian, radiative and thermal dispersion effects have been considered in the present analysis. The governing boundary layer equations and boundary conditions are cast into a dimensionless form and simplified by using a similarity transformation. The resulting system of equations is solved by using a double shooting Runge–Kutta method. The effect of viscosity index n, the conduction–radiation parameter R, the non-Darcy parameter Gr*, the thermal dispersion parameter Ds and the suction/injection parameter f_w on the fluid velocities, temperatures and the local Nusselt number are discussed.     

Forced convection heat transfer at the boundary layer of a packed bed

Numerical investigations of the nature of the fluid flow pattern and heat transfer at the boundary layer of a packed bed are reported. A volume averaged Navier-Stokes equation is used to predict the fluid flow and a volume averaged heat balance equation the heat transfer. A variable porosity in the packing is assumed in the region near the wall. Simulations are performed using a modified penalty Galerkin finite element method. The case of fully developed hydrodynamic flow and developing thermal flow is studied. The Nusselt number is found to depend on the Reynolds number, Graetz number and ratio of thermal conductivity of the solid and fluid phases. Comparison is made to some experimental literature values.Nomenclature A constant - [A] Navier-Stokes type matrix - B constant - [B] divergence matrix - C _p constant pressure heat capacity - d characteristic length - D _p particle diameter - D _t tube diameter - {F} solicitation vector - Gz Graetz number, z D _t ^–1 Pr _f Re _p - k permeability term - k _f Thermal conductivity of the fluid phase - k _s Thermal conductivity of the solid phase - [K] coefficient matrix for temperature equation - n normal vector - P pressure - Pr _f Prandtl number for the fluid _f C _p k _f ^-1 - r radial coordinate - R _t tube radius - R residual - R _m residual - Re _p Reynolds number for particle, - t tortuosity factor - T temperature - interstitial velocity - z axial coordinate - effective thermal conductivity - penalty parameter - boundary of solution domain - porosity - viscosity - density - test function - solution domain - test function     

Nusselt number and friction factor in thermally stratified turbulent channel flow under Non-Oberbeck–Boussinesq conditions

In stably stratified turbulence, computations under Oberbeck–Boussinesq (OB) hypothesis of temperature-independent fluid properties may lead to inaccurate representation of the flow field and to wrong estimates of momentum/heat transfer coefficients. This is clearly assessed here comparing direct numerical simulations of stratified turbulence under OB conditions to simulations under NOB (Non-Oberbeck–Boussinesq) conditions of temperature-dependent fluid viscosity and thermal expansion coefficient. Compared to the OB case, NOB conditions may induce local flow relaminarization with significant variations (up to 30%) of heat and momentum transfer coefficients. Together with DNS results, we propose a phenomenological model (based on turbulent bursts) for heat transfer prediction in stratified turbulence under OB and NOB conditions. Implications of NOB assumptions on mixing efficiency (i.e. flux Richardson number Ri_f) and turbulent Prandtl number (Pr_t) are also discussed. These results are of specific importance in RANS modelling, where the condition Pr_t = 1 is usually assumed (Reynolds analogy). Although this assumption is valid in some situations (i.e. boundary layer, pipe flow) there is uncertainty about its validity for stably-stratified turbulence. We demonstrate that this assumption is inaccurate when NOB effects become significant.     

Singular perturbation of boundary value problem for a vector fourth order nonlinear differential equation

We study the vector boundary value problem with boundary perturbations: ε~2y~((4))=f(x,y,y″,ε, μ) ( μ<χ<1-μ) y(χ,ε,μ)l_(χ-μ)= A_1(ε,μ), y(χ,ε,μ)l_(χ-1-μ)=B_1(ε,μ) y″(χ,ε,μ)l_(χ-μ)=A_2(ε,μ),y″(χ,ε,μ)l_(χ-1-μ)=B_2(ε,μ)where yf, A_j and B_j (j=1,2) are n-dimensional vector functions and ε,μ are two small positive parameters. This vector boundary value problem does not appear to have been studied, although the scalar boundary value problem has been treated. Under appropriate assumptions, using the method of differential inequalities we find a solution of the vector boundary value problem and obtain the uniformly valid asymptotic expansions.     

Transient behavior of heat removal from a cylindrical heat storage vessel packed with spherical porous particles

The transient heat transfer behavior in the case of heat removal from a cylindrical heat storage vessel packed with spherical particles was investigated experimentally for various factors (flow rate, diameter of spherical particles packed, temperature difference between flowing cold air and spherical particles accumulating heat, and physical properties of spherical particles). The experiments were covered in ranges of Reynolds number based on the mean diameter of spherical particles packed Re_d = 10.3–2200, porosity?=0.310 to 0.475, ratio of spherical particle diameter to cylinder diameterd/D = 0.0075–0.177 and ratio of length of the cylinder to cylinder diameterL/D=2.5–10. It was found that especially the flow rate and the dimension of spherical particles played an important role in estimating the transient local heat transfer characteristics near the wall of the cylindrical vessel in the present heat storage system. As flow rate and diameter of spherical particles were increased under a given diameter of the cylinder heat storage vessel, the mean heat transfer coefficient between the flow cold air and the hot spherical particles increased and the time period to finish removing heat from the vessel reduced. In addition, the useful experimental correlation equations of mean heat transfer coefficient between both phases and the time period to finish removing heat from the vessel were derived with the functional relationship of Nusselt numberNu _d=f [modified Prandtl numberPr ^* (d/D), Re_d) and Fourier numberFo = f(d/D, L/D, Pr^*, Re_d).     

Effects of Knudsen number and geometry on gaseous flow and heat transfer in a constricted microchannel

A flow and heat transfer numerical simulation is performed for a 2D laminar incompressible gas flow through a constricted microchannel in the slip regime with constant wall temperature. The effects of rarefaction, creeping flow, first order slip boundary conditions and hydrodynamically/thermally developing flow are assumed. The effects of Knudsen number and geometry on thermal and hydrodynamic characteristics of flow in a constricted microchannel are explored. SIMPLE algorithm in curvilinear coordinate is used to solve the governing equations including continuity, energy and momentum with the temperature jump and velocity slip conditions at the solid walls in discretized form. The resulting velocity and temperature profiles are then utilized to obtain the microchannel C _f Re and Nusselt number as a function of Knudsen number and geometry. The results show that Knudsen number has declining effect on the C _f Re and Nusselt number in the constricted microchannel. In addition, the temperature jump on wall and slip velocity increase with increasing Knudsen number. Moreover, by decreasing the throttle area, the fluid flow characteristics experience more intense variations in the constricted region. To verify the code a comparison is carried out with available results and good agreement is achieved.     

Mixed convection and conduction heat transfer in open cavities

A numerical study has been carried out on mixed convection and conduction in open cavities. The study covers the Rayleigh number from 10 to 10⁶, the Reynolds number from 1 to 10³ andk _r =k _w /k _a from 1 to 100 forPr=0.72 (air) and cavity aspect ratioB from 0.8 to 1.3. The flow is assumed to be laminar and two-dimensional. The density variation is taken into account by the Boussinesq approximation. The controlvolume approach is used for solving the governing equations of conjugate heat transfer involving conduction in the walls. Streamlines and isotherms in the system are produced and the results are represented in terms of the Nusselt number as function of other parameters. Correlations are derived to calculate heat transfer through the cavity openings.     

Transient conjugated mixed-convective heat transfer in a vertical plate channel with one wall heated discretely

 This study presents a numerical solution of the unsteady conjugated mixed-convection heat transfer in a vertical plate channel with one wall suddenly subjected to either isoflux or isothermal discrete heat sources. The effects of the dimensionless heat source length H ₁, the dimensionless spacing between heat sources H ₂, the dimensionless channel length L, the dimensionless heated-plate thickness B _l, the wall-to-fluid conductivity ratio K and the ratio of Grashof number to Reynolds number Gr/Re on the interface heat flux, Nusselt number and bulk fluid temperature are discussed in detail. Results show that the discrete heating can cause the heat transfer direction conversely from the fluid to the heated plate during the transient period, which is more significant for the cases with larger L and H ₂. For the system with isoflux discrete heat sources, the time required to reach the steady-state is shorter for larger H ₂. While the trend is reverse for system with isothermal discrete heat sources. Additionally, a higher ratio of the input energy is axially conducted through the plate wall from heated sections to unheated regions for a larger H ₂ and B _l or smaller L. Received on 9 November 1998