Finite-difference analysis of MHD Stokes problem for a vertical infinite plate in dissipative fluid

Finite-difference solution of MHD flow past an impulsively started vertical infinite plate in an electrically conducting fluid has been presented on taking into account the viscous dissipative heat. Results for velocity and temperature are shown graphically whereas the numerical values of the skin-friction and the rate of heat transfer are entered in the table. The results are discussed in terms of the parameters M (the Hartmann number), G (the Grashof number, G>0, cooling of the plate by free convection, G<0, heating of the plate by free convection currents), E (the Eckert number) and P (the Prandtl number).Nomenclature B ₀ applied magnetic field - c _p specific heat at constant pressure - g acceleration due to gravity - k thermal conductivity - t time - T temperature of the fluid near the plate - T temperature of the fluid far away from the plate - U ₀ velocity of the plate - u velocity of the fluid - coefficient of volume expansion - kinematic viscosity - scalar electrical conductivity - coefficient of viscosity - density of the fluid     

Oscillatory hydromagnetic disc flow

Summary This paper presents the solution to the problem of determining the flow field and the fluctuating torque necessary to sustain the motion of a torsionally oscillating plate in a viscous conducting fluid subjected to a uniform axial field under the assumptions that the amplitude of the oscillation is small and that the magnetic Prandtl number Pr _m () is small enough to justify the neglect of induced fields. The analysis reveals that the field decreases the flow velocities, but increases the magnitude of the fluctuating torque.     

Transient free convection flow of an incompressible viscous dissipative fluid

A finite-difference solution of transient free convection flow of a viscous dissipative fluid past an infinite vertical plate, on taking into account viscous dissipative heat is presented. Velocity profiles, temperature profiles are shown for different values of Pr, the Prandtl number and E, the Eckert number. The numerical values of the skin-friction and the rate of heat transfer are entered in a Table. It is observed that greater viscous dissipative heat causes a rise in the velocity, temperature and the skin-friction and a fall in the rate of heat transfer. An increase in Pr leads to a fall in the velocity, temperature and the skin-friction but the rate of heat transfer increases with increasing Pr.     

Heat transfer in the flow of a viscoelastic fluid over a stretching surface

An analysis is made of heat transfer in the boundary layer of a viscoelastic fluid flowing over a stretching surface. The velocity of the surface varies linearly with the distance x from a fixed point and the surface is held at a uniform temperature T _w higher than the temperature T _∞ of the ambient fluid. An exact analytical solution for the temperature distribution is found by solving the energy equation after taking into account strain energy stored in the fluid (due to its elastic property) and viscous dissipation. It is shown that the temperature profiles are nonsimilar in marked contrast with the case when these profiles are found to be similar in the absence of viscous dissipation and strain energy. It is also found that temperature at a point increases due to the combined influence of these two effects in comparison with its corresponding value in the absence of these two effects. A novel result of this analysis is that for small values of x, heat flows from the surface to the fluid while for moderate and large values of x, heat flows from the fluid to the surface even when T _w >T _∞. Temperature distribution and the surface heat flux are determined for various values of the Prandtl number P, the elastic parameter K ₁ and the viscous dissipation parameter a. Numerical solutions are also obtained through a fourth-order accurate compact finite difference scheme. Received on 14 October 1997     

Unsteady mixed convection on the stagnation-point flow adjacent to a vertical plate with a magnetic field

An analysis is performed to study the unsteady combined forced and free convection flow (mixed convection flow) of a viscous incompressible electrically conducting fluid in the vicinity of an axisymmetric stagnation point adjacent to a heated vertical surface. The unsteadiness in the flow and temperature fields is due to the free stream velocity, which varies arbitrarily with time. Both constant wall temperature and constant heat flux conditions are considered in this analysis. By using suitable transformations, the Navier–Stokes and energy equations with four independent variables (x, y, z, t) are reduced to a system of partial differential equations with two independent variables (, ). These transformations also uncouple the momentum and energy equations resulting in a primary axisymmetric flow, in an energy equation dependent on the primary flow and in a buoyancy-induced secondary flow dependent on both primary flow and energy. The resulting system of partial differential equations has been solved numerically by using both implicit finite-difference scheme and differential-difference method. An interesting result is that for a decelerating free stream velocity, flow reversal occurs in the primary flow after certain instant of time and the magnetic field delays or prevents the flow reversal. The surface heat transfer and the surface shear stress in the primary flow increase with the magnetic field, but the surface shear stress in the buoyancy-induced secondary flow decreases. Further the heat transfer increases with the Prandtl number, but the surface shear stress in the secondary flow decreases.     

Heat transfer from a slowly rotating sphere

The problem of steady state forced convection heat transfer in a viscous incompressible fluid occupying the annular region between two concentric spheres is considered. The inner sphere is maintained at a constant temperatureT ₀ and rotates slowly around an axis through the centre. The outer sphere is at rest and the temperature of its surface is prescribed as a function of the spherical coordinates and. It is shown that, when viscous dissipation is small, the overall rate of heat transfer from the rotating sphere into the fluid is unaffected by convection from the sphere surface, in case of a slow rotation, where the Stokes solution holds.     

Heat transfer in a laminar wall jet over a curved surface

The steady state heat transfer characteristics of the wall jet over a curved surface are obtained for constant wall temperature and constant wall heat flux boundary conditions. Both concave and convex curvatures have been considered. Numerical results for the temperature distribution are obtained and solutions for the wall values of the temperature functions have been tabulated for Prandtl number ranging from 0.01 to 100 while the curvature parameter was varied from –0.03 to 0.07.Nomenclature f velocity profile function - h heat transfer coefficient - K thermal conductivity - Nu Nusselt number - Pr Prandtl number - q _w heat flux at the wall - Re Reynolds number - R ₀ surface radius of curvature - T temperature - U characteristic velocity - u velocity component in x direction - v velocity component in y direction - x distance parallel to the surface - y distance normal to the surface - curvature parameter - dimensionless coordinate - dimensionless temperature - dynamic viscosity - kinematic viscosity - fluid density - shear stress - w conditions at the wall - conditions far away from the surface     

Velocity and temperature distributions on a semi-infinite flat plate embedded in a saturated porous medium

In this paper the velocity and temperature distributions on a semi-infinite flat plate embedded in a saturated porous medium are obtained for the governing equations (Kaviany [7]) following the technique adopted by Chandrashekara [2] which are concerned with the interesting situations of the existence of transverse, velocity and thermal boundary layers. Here the pressure gradient is just balanced by the first and second order solid matrix resistances for small permeability and observed that by increasing of the flow resistance the asymptotic value for the heat transfer rate increases. Further we concluded that the transverse boundary layers are thicker than that of axial boundary layers. Hence we evaluated the expressions for the boundary layer thickness, the shear stress at the semi-infinite plate and _T (the ratio of the thicknesses of the thermal boundary layer and momentum boundary layer). The variations of these quantities for different values of the porous parameterB and the flow resistanceF have been discussed in detail with the help of tables. The curves for velocity and temperature distributions have been plotted for different values ofB andF.Lastly we have evaluated the heat fluxq(x) and found that it depends entirely upon the Reynolds numberRe, Prandtl numberPr,B andF.     

Unsteady Couette flow in a rotating system

We have studied the unsteady Couette flow of a viscous incompressible fluid confined between parallel plates, rotating with an uniform angular velocity about an axis normal to the plates. The flow is induced by the motion of the upper plate and the fluid and plates rotate in unison with the same constant angular velocity. An exact solution of the governing equations have been obtained for small and large time $\tau$ by applying Laplace transform technique. It is found that the primary velocity decreases with increase in rotation parameter for small as well as large time. It is interesting to note that a back flow occurs in the region $0.0⩽\eta ⩽0.7$ for large time with increase in K when $K=4$ and 5. The secondary velocity increases in magnitude for small time with increase in rotation parameter. It is observed that the secondary velocity increases in magnitude for small values of rotation parameter. On the other hand, for large values of rotation parameter $K^2$, it decreases near the stationary plate and increases near the moving plate. The shear stress due to primary flow decreases with increase in rotation parameter $K^2$. On the other hand, it increases due to secondary flow with increase in rotation parameter for small time. It is noticed that for large time there exists separation in the primary and secondary flows due to high rotation.     

Temperature-dependent heat sources or sinks in a stagnation point flow

Summary An analytical study has been made to determine the heat transfer characteristics of a stagnation point flow in which there are temperature-dependent heat sources or sinks. Results have been obtained for both strong and weak sources or sinks for a Prandtl number of 0.7. An analytical method, applicable to all Prandtl numbers, was utilized which circumvented the need for extensive numerical solutions and which, at the same time, provided a closed-form representation for the heat transfer. A few numerical solutions were carried out to verify the method.Nomenclature a _i constants depending on Prandtl number - c _p specific heat at constant pressure - f dimensionless velocity variable - g function defined by equation (13) - g _n functions of (n=1, 2, 3,...) - k thermal conductivity - Pr Prandtl number, c _p /k - q heat transfer rate per unit area at surface - Q heat flux parameter, q/k(u ₁/)^1/2 - S rate of heat generation or removal per unit volume (divided by c _p ) - T static temperature; T _w , wall temperature; T , free-stream temperature - u ₁ proportionality constant for free-stream velocity - U free-stream velocity - v normal velocity component - x coordinate measuring distance along surface from stagnation point - y coordinate measuring distance normal to surface - heat generation parameter, equation (3) - dimensionless normal coordinate, - dimensionless temperature - _n functions of (n=1, 2, 3,...) - absolute viscosity - kinematic viscosity - density     

Laminar natural convection heat transfer from vertical plate to power-law fluid

Summary Approximate solutions for laminar natural convection heat transfer between a vertical plate and a power-law fluid with high Prandtl number were obtained using an integral method for cases with various types of boundary conditions. The results were found in good agreement with available experimental evidence.Nomenclature a exponent defined by equations (28) and (29) - A, B, C, D, E constants defined by equations (15) to (19) - C ₁, C ₂, M ₁, M ₂ coefficients for Nusselt number expression defined by (32b), (33b) - f temperature difference, equal to T _s–T - f ⁺ dimensionless temperature difference - g gravitational acceleration - Gr Grashof number defined by (25), (50) and (66), respectively - H heat flux at plate surface - h _x local heat transfer coefficient - K consistency index for Power-law fluid - k thermal conductivity of fluid - K ₁, K ₂ constants defined by (50) and (51) - L height of plate - n flow behavior index for Power-law fluid - P a quantity defined by (54a) - T temperature - T _s plate temperature - T temperature of the bulk of fluid - s constant given by (35) - u velocity component along x-direction - u _x maximum velocity induced by natural convection current, (10) - v velocity component along y-direction - x distance measured along direction parallel to that of gravitational force - x ⁺ dimensionless quantity, defined as x/L - y distance measured away from plate - Nu _x local Nusselt number - Nu _av average Nusselt number - Pr Prandtl number defined by (24) - T temperature difference according to boundary conditions - thermal diffusivity of fluid - coefficient of thermal expression of fluid - boundary layer thickness - ⁺ dimensionless boundary layer thickness - dimensionless velocity profile - dimensionless variable, defined as y/ - dimensionless temperature difference     

On heat transfer effects of a viscous fluid squeezed and extruded between two parallel plates

The unsteady squeezing and extrusion of a viscous fluid between two parallel plates of constant temperature is examined. The dimensionless extrusion parameter,=U/V, is introduced to represent the effects of the extrusion on the squeezing velocities. The squeezing parameter=VH/, represents the effect of the inertial forces on heat and fluid flow characteristics. It is found that increasing the extrusion parameter will increase both the velocity and the heat transfer rates to the viscous fluid. Increasing the squeezing parameter had also decreased the fluid velocity and enhanced heat transfer rates. Increasing the viscous effects or the Eckert number E=U²/c_p (T_E – T_s) heated the fluid and consequently decreased the heat transfer rates. Different velocity profiles, temperature profiles, and Nusselt numbers against various dimensionless groups are drawn.     

Flow and Heat Transfer of Viscous Fluids Saturated in Porous Media over a Permeable Non-isothermal Stretching Sheet

The present paper deals with the flow and heat transfer of a viscous fluid saturated in a porous medium past a permeable and non-isothermal stretching sheet with internal heat generation or absorption and radiation. Closed-form solutions to steady, two dimensional momentum equations with neglecting quadratic inertia terms and heat transfer equation are found using a similarity transformation. Asymptotic expressions of the temperature functions are also presented valid for both very large and very small modified Prandtl numbers. Attention is focused on the effects of porous parameter K, suction parameter R, radiation parameter Nr, viscosity ratio Λ, internal heat parameter α and Prandtl number P to the characteristics of flow and heat transfer.     

Mixed convection flow over a vertical plate with localized heating (cooling), magnetic field and suction (injection)

An analysis is carried out to study the effects of localized heating (cooling), suction (injection), buoyancy forces and magnetic field for the mixed convection flow on a heated vertical plate. The localized heating or cooling introduces a finite discontinuity in the mathematical formulation of the problem and increases its complexity. In order to overcome this difficulty, a non-uniform distribution of wall temperature is taken at finite sections of the plate. The nonlinear coupled parabolic partial differential equations governing the flow have been solved by using an implicit finite-difference scheme. The effect of the localized heating or cooling is found to be very significant on the heat transfer, but its effect on the skin friction is comparatively small. The buoyancy, magnetic and suction parameters increase the skin friction and heat transfer. The positive buoyancy force (beyond a certain value) causes an overshoot in the velocity profiles.A mass transfer constant - B magnetic field - C_fx skin friction coefficient in the x-direction - C_p specific heat at constant pressure, kJ.kg^–1.K - C_v specific heat at constant volume, kJ.kg^–1.K^–1 - E electric field - g acceleration due to gravity, 9.81 m.s^–2 - Gr Grashof number - h heat transfer coefficient, W.m².K^–1 - Ha Hartmann number - k thermal conductivity, W.m^–1.K - L characteristic length, m - M magnetic parameter - Nu_x local Nusselt number - p pressure, Pa, N.m^–2 - Pr Prandtl number - q heat flux, W.m^–2 - Re Reynolds number - Re_m magnetic Reynolds number - T temperature, K - T_o constant plate temperature, K - u,v velocity components, m.s^–1 - V characteristic velocity, m.s^–1 - x,y Cartesian coordinates - thermal diffusivity, m².s^–1 - coefficient of thermal expansion, K^–1 - , transformed similarity variables - dynamic viscosity, kg.m^–1.s^–1 - ₀ magnetic permeability - kinematic viscosity, m².s^–1 - density, kg.m^–3 - buoyancy parameter - electrical conductivity - stream function, m².s^–1 - dimensionless constant - dimensionless temperature, K - w, conditions at the wall and at infinity     

Analysing flow and heat transfer of a viscoelastic fluid over a semi-infinite horizontal moving flat plate

This paper presents a numerical study of the flow and heat transfer of an incompressible homogeneous second grade type fluid above a flat plate moving with constant velocity U. Such a viscoelastic fluid is at rest and the motion is created by the sheet. The effects of the non-Newtonian nature of the fluid are governed by the local Deborah number K (the ratio between the relaxation time of the fluid and the characteristic time of the flow). When , a new analytical solution for this flow is presented and the effects of fluid's elasticity on flow characteristics, dimensionless stream function and its derivatives are analysed in a wide domain of K. A novel result of the analysis is that a change in the flow solution's behaviour occurs when the dimensionless stream function at the edge of the boundary layer, f_∞, equals 1.0. It is found that velocity at a point decreases with increase in the elasticity of the fluid and, as expected, the amount of fluid entrained diminishes when the effects of fluid's elasticity are augmented. In our heat transfer analyses we assume that the surface temperature has a power-law variation. Two cases are studied, namely, (i) the sheet with prescribed surface temperature (PST case) and (ii) the sheet with prescribed heat flux (PHF case). Local similarity heat-transfer solutions are given for PST case when s=2 (the wall temperature parameter) whereas when a similarity solution takes place in the case of prescribed wall heat flux. The numerical results obtained are fairly in good agreement with the aforementioned analytical ones.     

Transient temperature response of a porous matrix entering a planetary atmosphere

Summary An analysis is made of the transient temperature behavior of a transpiration-cooled porous matrix entering a planetary atmosphere with constant velocity and negative entry angle. The analysis is based on one dimensional heat conduction in a porous plate subjected to a time dependent heat flux at one side and cooled internally by mass injection from a constant temperature reservoir at the opposite side. An exact closed-form solution is obtained and temperature charts are presented for a wide range of Fourier number and coolant flow parameters.Nomenclature A surface area, ft² - C constant, 17,600 Btu/ft^3/2-sec - C _c constant pressure specific heat of coolant, Btu/lb_m-°F - g local gravitational acceleration, ft/sec² - g _c coolant flow parameter, defined by equation (15) - h height of entry above planet surface, ft - K ₁ ratio of local heat flux to stagnation point heat flux - K thermal conductivity of plate material, Btu/sec-ft-°F - L plate thickness, ft - m constant, 3.15 - m _c coolant mass flow rate, lb_m/sec - M _n roots of equation (33) - n constant, 0.50 - N defined by equation (30) - P porosity - q surface heat flux, Btu/ft²-sec - q ₀ surface heat flux at t=0, Btu/ft²-sec, defined by (6) - r distance from planet center, ft - R radius of curvature at stagnation point, ft - t time, sec - T temperature, °F - T _c coolant supply temperature, °F - V velocity, ft/sec - x normal coordinate through plate, ft - y altitude, ft - thermal diffusivity of plate, ft²/sec - atmospheric density decay parameter, 1/23500 ft^–1 - flight path angle relative to local horizontal direction, positive for climbing and negative for descent, deg - dimensionless temperature parameter, defined by (12) - dimensionless distance, defined by (13) - free stream atmospheric density, slug/ft³ - ₀ atmospheric density at reference state, slug/ft³ - Fourier number, defined by (14) - 1/sec, defined by (7) - flight entry parameter, defined by (16)     

Thermal Capacity Effect on Transient Free Convection Adjacent to a Vertical Surface in a Porous Medium

In this paper we analyse how the presence of the thermal capacity of a vertical flat plate of finite thickness, which is embedded in a porous medium affects the transient free convection boundary-layer flow. At the time t = 0, the plate is suddenly loaded internally with a constant heat flux rate q, so that a transient boundary-layer flow is initiated adjacent to the plate. Initially, the transient effects due to the imposition of the uniform heat flux rate at the plate are confined to a thin fluid region near to the surface and are described by a small time solution. These effects continue to penetrate outwards and eventually evolve into a new steady state flow. Analytical solutions have been derived for these transient (small time) and steady state (large time) flow regimes, which are then matched by a numerical solution of the full boundary-layer equations. It has been found that the non-dimensional fluid temperature (or fluid velocity) profiles are reduced when the thermal capacity effects, described by a parameter Q ^*, are reduced. For small values of Q ^*, the approach of these profiles to their steady state values is monotonic. However, for large values of Q ^*, the temperature profiles are observed to locally exceed (pass through a maximum value) the final steady state values at certain distances from the plate. In general, the maxima in the temperature profiles increase in size as Q ^* increases and the time taken to approach the steady state solutions increases significantly.     

Unified treatment of free convection adjacent to a vertical plate with three thermal boundary conditions

Free convection adjacent to a vertical plate is considered with three different boundary conditions, namely, the plate is subjected to a prescribed temperature, a prescribed heat flux or a prescribed heat transfer coefficient. By a unified treatment of similarity analysis, the governing equations of free convection are reduced to identical system of ordinary differential equations for all the cases. It is shown that these equations are invariant under a certain transformation group and solution for one case can be used to obtain the solutions for the other two cases by a simple method. The critical cases are found for which the solutions for all the three cases are identical.

---

Einheitliche Behandlung der freien Konvektion längs einer vertikalen Platte für drei verschiedene Randbedingungen

Zusammenfassung Die freie Konvektion längs einer vertikalen Platte wird für drei verschiedene Randbedingungen betrachtet: einmal wird die Platte auf einer bestimmten Temperatur gehalten, dann einem bestimmten Wärmestrom ausgesetzt und letztens mit einem bestimmten Wärmeübergangskoeffizienten untersucht. Bei einer einheitlichen Behandlung der Ähnlichkeitsanalyse werden für alle Fälle die bestimmenden Gleichungen für freie Konvektion zu einem identischen System von gewöhnlichen Differentialgleichungen reduziert. Es wird beweisen, daß diese Gleichungen unter bestimmten Transformationsgruppen invariant sind und Lösungen von einem Fall dazu benützt werden können, um die Lösungen der restlichen Fälle zu erhalten. Es werden die kritischen Fälle gefunden, für die alle Lösungen der drei entsprechenden Randbedingungen identisch sind.

---

Nomenclature c _o positive constant - c _p specific heat of the convective fluid - C function defined by Eq. (6a) - E energy convection by the boundary layer - f dimensionless stream function - g gravitational acceleration - G function defined by Eq. (6b) - k thermal conductivity - N heat transfer coefficient - PHF prescribed heat flux - PHTC prescribed heat transfer coefficient - Pr Prandtl number - PT prescribed temperature - q _w surface heat flux - T temperature - T _e ambient temperature - u velocity component inx direction - v velocity component iny direction - x coordinate in the vertical direction - y coordinate in the horizontal direction - coefficient of thermal expansion - dimensionless similarity variable - kinematic viscosity - dimensionless temperature - density of the convective fluid - stream function     

Combined natural and forced convection in a laminar wall jet along a vertical plate with uniform surface heat flux

The steady state flow and heat transfer characteristics of the combined natural and forced convection in a two dimensional, laminar, incompressible wall jet over a vertical wall are obtained for constant wall heat flux boundary condition. The velocity and temperature distribution are assumed to be power series, where the zeroth term corresponds to that for a plane wall jet in the absence of buoyancy effects. Numerical results for the momentum and thermal series functions are presented for a Prandtl number of 0.73. Wall values of the momentum and thermal series functions are presented for Prandtl numbers ranging from 0.01 to 1000.Nomenclature Gr* modified Grashof number - k thermal conductivity - Nu Nusselt number - Pr Prandtl number - q _w heat flux at the wall - Re Reynolds number - T temperature - u velocity component in x-direction - v velocity component in y-direction - x co-ordinate along the plane wall - y co-ordinate normal to the wall - () gamma function - non-dimensional co-ordinate defined in (6) - non-dimensional temperature - dynamic viscosity - kinematic viscosity - non-dimensional co-ordinate defined in (6) - density - w values at the wall - values at large distances away from the wall     

Mass transfer effects on flow past a vertical oscillating plate with variable temperature

An exact solution to the flow of a viscous incompressible fluid past an infinite vertical oscillating plate, in the presence of a foreign mass has been derived by the Laplace-transform technique when the plate temperature is linearly varying as time. The velocity profiles are shown on graphs and the numerical values of the skin-friction are listed in a table. It is observed that the skin-friction increases with increasingSc, t orGr but decreases with increasingGm ort, whereSc (the Schmidt number), (frequency),t (time),Gr (the Grashof number) andGm is (the modified Grashof number) andt.

---

Einfluß von Stoffübergang auf die Strömung entlang einer senkrechten oszillierenden Platte veränderlicher Temperaturen

Zusammenfassung Mit Hilfe der Laplace-Transformation wird eine exakte Lösung für die Strömung einer zähen, inkompressiblen Flüssigkeit entlang einer unendlich ausgedehnten, senkrechten, oszillierenden Platte gewonnen, wobei die Einwirkung eines Feststoffs Berücksichtigung findet und die Plattentemperatur linear mit der Zeit veränderlich sein soll. Die Geschwindigkeitsprofile sind in Diagrammen dargestellt und die numerischen Werte der Reibungsschubspannung in einer Tabelle. Letztere wächst mit der Schmidt-ZahlSc, der Grashof-ZahlGr und dem Produkt aus Frequenz und Zeitt; sie nimmt ab, wennGm (die modifizierte Grashof-Zahl) undt zunehmen.

---

Nomenclature C species concentration in the fluid near the plate - C ^– species concentration in the fluid away from the plate - C _W species concentration at the plate - C _p specific heat at constant pressure - D chemical molecular diffusivity - Gm modified Grashof number - Gr Grashof number - g acceleration due to gravity - K thermal conductivity - P Prandtl number - Sc Schmidt number - T Temperature of the fluid near the plate - T _W temperature of the plate - T temperature of the fluid far away from the plate - t time - u velocity of the fluid in the upward direction - U ₀ amplitude of oscillation - x coordinate axis along the plate in the vertically upward direction - y coordinate axis normal to the plate Greek symbols skin-friction - viscosity - coefficient of volume expansion - ^* coefficient of species expansion - density - frequency