Combined effects of internal heat generation and higher order chemical reaction on the non‐darcian forced convective flow of a viscous incompressible fluid with variable viscosity and thermal conductivity over a stretching surface embedded in a porous medium

In this paper, we study the combined effects of internal heat generation and higher order chemical reaction on a steady two‐dimensional non‐Darcian forced convective flow of a viscous incompressible fluid with variable dynamic viscosity and thermal conductivity in a fluid saturated porous medium passing over a linear stretching sheet. Using similarity transformations, the governing nonlinear‐coupled partial differential equations are made dimensionless and solved numerically for similarity solutions using very robust computer algebra software Maple 8. The non‐dimensional velocity, temperature and concentration distributions are presented graphically for various pertinent parameters such as relative temperature difference parameter, Darcy number, porosity parameter, reaction rate parameter and the order of the chemical reaction. The variations of Prandtl number and Schmidt number within the boundary layer are also displayed graphically when the fluid dynamic viscosity and thermal conductivity are temperature dependent. From the present numerical computations it is found that Prandtl number as well as Schmidt number must be taken as variables within the flow domain when the fluid's dynamic viscosity and thermal conductivity are variable. In the presence of internal heat generation, dynamic viscosity and thermal conductivity of the fluid are found to be higher than when it is absent. Increasing Darcy number reduces dynamic viscosity as well as thermal conductivity whereas increasing pore size reduces the Schmidt number and increases the Prandtl number within the boundary layer. For higher order reaction the rate of increase in mass transfer function is less compared to the rate of increase for the lower order reaction. © 2011 Canadian Society for Chemical Engineering     

MHD flow and heat transfer of a micropolar fluid over a nonlinear stretching surface with variable surface heat flux and heat generation

An analysis has been carried out to study magnetohydrodynamic boundary layer flow and heat transfer of an electrically conducting micropolar fluid over a nonlinear stretching surface with variable wall heat flux in the presence of heat generation/absorption and a non‐uniform transverse magnetic field. The governing system of partial differential equations is first transformed into a system of ordinary differential equations using similarity transformation. The transformed equations are solved numerically. Results for the dimensionless velocity, micro‐rotation, and temperature profiles are displayed graphically delineating the effects of various parameters characterising the flow. The results show that the velocity profile decreases as the magnetic parameter and the velocity exponent increase, while it increases as the material parameter increases. The results show also that the temperature profile increases as the magnetic parameter, the velocity exponent, and the heat generation parameter increase. Furthermore, the temperature profile decreases as the material parameter, the heat absorption parameter, and the Prandtl number increase.     

Combined effects of internal heat generation and buoyancy force on boundary layer flow over a vertical plate with a convective surface boundary condition

This paper considers the effect of buoyancy force and internal heat generation on laminar thermal boundary layer over a vertical plate with a convective surface boundary condition. We assumed that left surface of the plate is in contact with a hot fluid while a stream of cold fluid flows steadily over the right surface with a heat source that decays exponentially. Using a similarity variable, the steady state governing non‐linear partial differential equations have been transformed into a set of coupled non‐linear ordinary differential equations, which are solved numerically by applying shooting iteration technique together with fourth order Runge–Kutta integration scheme. The effects of Prandtl number, local Biot number, the internal heat generation parameter and the local Grashof number on the velocity and temperature profiles are illustrated and interpreted in physical terms. A comparison with previously published results on special case of the problem shows excellent agreement. From our results, an overshoot of fluid velocity within the boundary layer is observed due to combined effect of buoyancy force and internal heat generation, in addition, internal heat generation causes thickening of thermal boundary layer. © 2011 Canadian Society for Chemical Engineering     

UNSTEADY FLUID AND HEAT FLOW INDUCED BY A STRETCHING SHEET WITH MASS TRANSFER AND CHEMICAL REACTION

The effect of chemical reaction on the flow, heat, and mass transfer within a viscous fluid on an unsteady stretching sheet is examined. The stretching rate, temperature and concentration of the sheet, and the chemical reaction rate are assumed to vary with time. The time-dependent boundary layer equations governing the flow are reduced through a convenient similarity transformation to a set of ordinary differential equations, which are numerically solved by applying the fourth-order Runge-Kutta-Fehlberg scheme with the shooting technique. Results for the velocity, temperature, and concentration distributions as well as the wall temperature and concentration gradients are presented graphically for various values of the unsteadiness parameter A, Prandtl number Pr, Schmidt number Sc, and chemical reaction parameter γ.     

Mixed convection flow of non‐Newtonian fluid from a slotted vertical surface with uniform surface heat flux

In the present paper, the combined convection flow of an Ostwald–de Waele type power‐law non‐Newtonian fluid past a vertical slotted surface has been investigated numerically. The boundary condition of uniform surface heat flux is considered. The equations governing the flow and the heat transfer are reduced to local non‐similarity form. The transformed boundary layer equations are solved numerically using implicit finite difference method. Solutions for the heat transfer rate obtained for the rigid surface compare well with those documented in the published literature. From the present analysis, it is observed that, an increase in χ leads to increase in skin friction as well as reduction in heat transfer at the surface. As the power‐law index n increases, the friction factor as well as heat transfer increase.     

Three‐Dimensional Simulation of Bubble Formation Through a Microchannel T‐Junction

Three‐dimensional simulations of bubble formation in Newtonian and non‐Newtonian fluids through a microchannel T‐junction are conducted by the volume‐of‐fluid method. For Newtonian fluids, the critical capillary number Ca for the transition of the bubble breakup mechanism is dependent on the velocity ratio between the two phases and the microchannel dimension. For the power law fluid, the bubble diameter decreases and the generation frequency increases with higher viscosity parameter K and power law index n. For a Bingham fluid, the viscous force plays a more important role in microbubble formation. Due to the yield stress τ_y, a high‐viscous region is developed in the central area of the channel and bubbles deform to a flat ellipsoid shape in this region. The bubble diameter and generation frequency are almost independent of K.     

Mixed convection heat transfer in porous media in the non-darcy regime

In this study mixed convection heat transfer in a homogeneous porous duct of square cross section in a horizontal orientation is examined. Results from a generalized Forchheimer model are compared with that from the Darcy model. The heat transfer rate and the flow behavior depend on the following parameters: Grashof number, Gr = Q'gβKa/kv², an axial flow pressure drop parameter, ζ = (aK/vμ)dp'/dz', an inertial parameter ξ = mK/a, appearing in the Forccheimer model and the Prandtl number, Pr = C_pμ/k. In the Darcy limit, ξ → 0, the role of the axial flow parameter, λ is reduced to a mere scale factor and the flow behavior is determined by a single parameter, λ = Gr · Pr. Both the Darcy and the Forchheimer models exhibit dual solutions and a hysteresis behavior over a certain range of Gr. Such parametric dependence can be used as an additional tool along with carefully designed experiments to determine the importance of inertial and Prandtl number effects on convective heat transfer in porous media.     

Isocyanate trimerization kinetics and heat transfer in structural reaction injection molding

Guidelines are developed for molding large composite parts via structural reaction injection molding using glass preforms and polyisocyanurate resins. These are based on numerical simulations of the simultaneous heat transfer and reaction kinetics of a commercial system during and after mold filling. Premised requirements are that resin does not gel before the mold is filled, yet, reactions are sufficiently vigorous to approach completion. An existing mechanistic kinetic model is used and material parameters found from a chemical kinetics study employing an insulated cup. It is found desirable to use a high mold temperature and a low preform temperature in molding. Nondimensionalization of the governing equations reveals the existence of a Nusselt number (N_u), which describes the relative importance of heat transfer between resin and glass relative to thermal diffusion to the mold wall. With a Nusselt number of about 50 or higher it is possible to use the cooling capacity of the preform to extend gel time. The magnitude of N_u is influenced by part thickness, glass fraction, strand diameter, and flow velocity. Thus, the effect of the preform on extending resin gel time is within control of the molder.     

NUMERICAL MODELING OF MHD CONVECTIVE HEAT AND MASS TRANSFER IN PRESENCE OF FIRST-ORDER CHEMICAL REACTION AND THERMAL RADIATION

An analysis was carried out numerically to study unsteady heat and mass transfer by free convection flow of a viscous, incompressible, electrically conducting Newtonian fluid along a vertical permeable plate under the action of transverse magnetic field taking into account thermal radiation as well as homogeneous chemical reaction of first order. The fluid considered here is an optically thin gray gas, absorbing-emitting radiation, but a non-scattering medium. The porous plate was subjected to a constant suction velocity with variable surface temperature and concentration. The dimensionless governing coupled, nonlinear boundary layer partial differential equations were solved by an efficient, accurate, extensively validated, and unconditionally stable finite difference scheme of the Crank-Nicolson type. The velocity, temperature, and concentration fields were studied for the effects of Hartmann number (M), radiation parameter (R), chemical reaction (K), and Schmidt number (Sc). The local skin friction, Nusselt number, and Sherwood number are also presented and analyzed graphically. It is found that velocity is reduced considerably with a rise in the magnetic body parameter (M), whereas the temperature and concentration are found to be markedly boosted with an increase in the magnetic body parameter (M). An increase in the conduction-radiation parameter (R) is found to escalate the local skin friction (τ), Nusselt number, and concentration, whereas an increase in the conduction-radiation parameter (R) is shown to exert the opposite effect on either velocity or temperature field. Similarly, the local skin friction and the Sherwood number are both considerably increased with an increase in the chemical reaction parameter. Possible applications of the present study include laminar magneto-aerodynamics, materials processing, and MHD propulsion thermo-fluid dynamics.     

A bifurcation study of convective heat transfer in a hele-shaw cell

Steady-state multiplicity characteristics of convective heat transfer within a Hele-Shaw cell are investigated. The Navier-Stokes equations and the energy equation are averaged across the narrow gap, d, of the cell. The resulting two-dimensional, stationary equations depend on the following parameters: (i) the length to height aspect ratio γ, (ii) the tilt anle ? (iii) the Prandtl number Pr, (iv) an inertia parameter ξ = d²/ 12a², and (v) the Grashof number. Gr = Q_gβga⁵/kv². Here a is the height of the cell and Q, is the heat generation rate per unit volume. The complete structure of symmetric and asymmetric stationary solutions are traced using recent algorithms from bifurcation theory. In the double limit of ξ → 0 and Gr → ∞ such that Ra = 4GrPrξ remains finite (where Ra is the Rayleigh number for the Darcy model) the Hele-Shaw model reduces to that of the Darcy model.     

Mixed convective boundary layer flow over a vertical wedge embedded in a porous medium saturated with a nanofluid: Natural Convection Dominated Regime

A boundary layer analysis is presented for the mixed convection past a vertical wedge in a porous medium saturated with a nano fluid. The governing partial differential equations are transformed into a set of non-similar equations and solved numerically by an efficient, implicit, iterative, finite-difference method. A parametric study illustrating the influence of various physical parameters is performed. Numerical results for the velocity, temperature, and nanoparticles volume fraction profiles, as well as the friction factor, surface heat and mass transfer rates have been presented for parametric variations of the buoyancy ratio parameter Nr, Brownian motion parameter Nb, thermophoresis parameter Nt, and Lewis number Le. The dependency of the friction factor, surface heat transfer rate (Nusselt number), and mass transfer rate (Sherwood number) on these parameters has been discussed.     

Shell and tube heat exchanger design using mixed‐integer linear programming

The design of heat exchangers, especially shell and tube heat exchangers was originally proposed as a trial and error procedure where guesses of the heat transfer coefficient were made and then verified after the design was finished. This traditional approach is highly dependent of the experience of a skilled engineer and it usually results in oversizing. Later, optimization techniques were proposed for the automatic generation of the best design alternative. Among these methods, there are heuristic and stochastic approaches as well as mathematical programming. In all cases, the models are mixed integer non‐linear and non‐convex. In the case of mathematical programming solution procedures, all the solution approaches were likely to be trapped in a local optimum solution, unless global optimization is used. In addition, it is very well‐known that local solvers need good initial values or sometimes they do not even find a feasible solution. In this article, we propose to use a robust mixed integer global optimization procedure to obtain the optimal design. Our model is linear thanks to the use of standardized and discrete geometric values of the heat exchanger main mechanical components and a reformulation of integer nonlinear expressions without losing any rigor. © 2016 American Institute of Chemical Engineers AIChE J, 63: 1907–1922, 2017     

Modeling and numerical analysis of nanoliquid (titanium oxide,graphene oxide) flow viscous fluid with second order velocity slip and entropy generation

The prime objective of the present communication is to examine the entropy-optimized second order velocity slip Darcy–Forchheimer hybrid nanofluid flow of viscous material between two rotating disks.Electrical conducting flow is considered and saturated through Darcy–Forchheimer relation. Both the disks are rotating with different angular frequencies and stretches with different rates. Here graphene oxide and titanium dioxide are considered for hybrid nanoparticles and water as a continuous phase liquid. Joule heating, heat generation/absorption and viscous dissipation effects are incorporated in the mathematical modeling of energy expression. Furthermore, binary chemical reaction with activation energy is considered. The total entropy rate is calculated in the presence of heat transfer irreversibility, fluid friction irreversibility,Joule heating irreversibility, porosity irreversibility and chemical reaction irreversibility through thermodynamics second law. The nonlinear governing equations are first converted into ordinary differential equations through implementation of appropriate similarity transformations and then numerical solutions are calculated through Built-in-Shooting method. Characteristics of sundry flow variables on the entropy generation rate, velocity, concentration, Bejan number, temperature are discussed graphically for both graphene oxide and titanium dioxide hybrid nanoparticles. The engineering interest like skin friction coefficient and Nusselt number are computed numerically and presented through tables. It is noticed from the obtained results that entropy generation rate and Bejan number have similar effects versus diffusion parameter. Also entropy generation rate is more against the higher Brinkman number.     

Natural convection in a heat generating fluid bounded by two horizontal concentric cylinders

Convective flow and heat transfer of a Boussinesq fluid contained between two horizontal concentric cylinders is investigated under the effects of two driving mechanisms – an externally-imposed temperature gradient across the annulus, and a uniform internal heat generation. Numerical results for flow field and temperature distribution are obtained in terms of four dimensionless parameters, namely the radius ratio, R, the Prandtl number, Pr, the Rayleigh number, Ra*, and the ratio, S, between the characteristic temperature induced by internal heating and the applied temperature difference between the boundaries. Depending on the value of S, the flow pattern is made up of either one or two vortices in each half cavity, and heat is transferred into or out of the cavity through the hot wall. In particular, for a certain value of the applied temperature difference, the hot wall apparently acts as a thermally-insulated boundary, the internal heat is completely lost through the cold wall, and the fluid undergoes a transition from a bicellular to a unicellular flow regime.     

Estimation of solubility parameter components of solutes and polymers using heat of vaporization and heat of sorption of solutes

In a recent study, a two‐dimensional solubility parameter model was used to correlate the heat of solution for solutes ranging from n‐alkanes to alcohols, dissolved in isotatic polypropylene (PP), poly(ethyl ethylene) (PEE), and poly(dimethylsiloxane) (PDMS). When literature data of solubility parameter components of solutes were used, the correlation had some scattering for solutes with low values of cohesive energy density. In this study, the components of solubility parameters of solutes and polymers were estimated from cohesive energy and heat of sorption of solutes. Good correlation was obtained for the specific heat of sorption (ΔU_sorp/V) for solutes ranging from n‐alkanes to alcohols, and PDMS had a polar component as previously estimated. Free volume effect in solution process may be the source of a small systematic deviation from the model. © 2009 Wiley Periodicals, Inc. J Appl Polym Sci, 2009     

Thermal Diffusion and MHD Effects on Combined Free‐Forced Convection and Mass Transfer of a Viscous Fluid Flow Through a Porous Medium with Heat Generation

The problem of thermal diffusion and magnetic field effects on combined free‐forced convection and mass transfer flow past a vertical porous flat plate, in the presence of heat generation is studied numerically. The governing momentum, energy and concentration equations are converted into a system of nonlinear ordinary differential equations by means of similarity transformations. The resulting system of coupled nonlinear ordinary differential equations is solved numerically by using the Shooting method. Numerical results are presented for velocity, temperature and concentration profiles within the boundary layer for different parameters of the problem including suction parameter, heat generation parameter, Soret number, Dufour number, magnetic parameter, etc. In addition, the effects of the pertinent parameters on the skin friction and the rates of heat and mass transfer are discussed numerically and illustrated graphically.     

圆形自由水射流冲击换热及喷嘴布置

引言 水射流冲击冷却由于具有较高的换热能力,广泛应用于机械和化工行业,以实现工件的快速冷却和控制工件的温度变化.在大型轴类工件(工件直径D=1000～3000 mm)喷水冷却装置中,多喷嘴圆孔自由水射流以特定阵列布置冲击至工件表面,相对短暂的沸腾换热结束后,阵列自由水射流即以强制对流方式实现工件冷却.因此,自由水射流冲击换热特性及多喷嘴布置形式对喷水冷却装置的结构设计至关重要.     

Analysis of solar‐driven gasification of biochar trickling through an interconnected porous structure

The efficient transfer of high‐temperature solar heat to the reaction site is crucial for the yield and selectivity of the solar‐driven gasification of biomass. The performance of a gas‐solid trickle‐bed reactor constructed from a high thermal conductivity porous ceramic packing has been investigated. Beech char particles were used as the model feedstock. A two‐dimensional finite‐volume model coupling chemical reaction with conduction, convection, and radiation of heat within the packing was developed and tested against measured temperatures and gasification rates. The sensitivity of the gasification rate and reactor temperatures to variations of the packing's pore diameter, porosity, thermal conductivity, and particle loading was numerically studied. A numerical comparison with a moving bed projected a more uniform temperature distribution and higher gasification rates due to the increased heat transfer via combined radiation and conduction through the trickle bed. © 2014 American Institute of Chemical Engineers AIChE J, 61: 867–879, 2015     

Analytical studies on the temperature distribution in spiral plate heat exchangers: Straightforward design formulae for efficiency and mean temperature difference

The temperature distribution in spiral plate heat exchangers has been calculated numerically to obtain the efficiency and the logarithmic mean temperature difference (LMTD) correction factors F as a function of the number of transfer units N, the number of turns n, and the heat capacity rate ratio C. It has been found that the LMTD correction factors, when plotted against the number of transfer units per turn N/n, fall approximately on a single curve. That curve for balanced countercurrent operation (C = −1) can be very closely represented by the simple formula F = (n/N)tanh(N/n). From that simple analytic representation of our numerical results it was concluded that a simpler physical model might exist to represent the overall behaviour of a spiral plate heat exchanger equally well. In fact, a countercurrent cascade of n cocurrent heat exchangers does result exactly in the above-mentioned formula for the LMTD correction factor. From that model the F-factors for other heat capacity rate ratios C (−1 < C 0) can also be calculated and they are in sufficient agreement with the numerical results.     

SIMULATION OF HYDRATION/DEHYDRATION OF CaO/Ca(OH)2 CHEMICAL HEAT PUMP REACTOR FOR COLD/HOT HEAT GENERATION

ABSTRACT

A chemical heat pump (CHP) utilizes reversible reactions involving significant endothermic and exothermic heats of reaction in order to develop a heat pump effect by storing and releasing energy while transforming it from chemical to thermal energy and vice versa. In this paper, we present a mathematical model and its numerical solution for the heat and mass transport phenomena occurring in the reactant particle bed of the CHP for heat storage and cold/hot heat generation based on the CaO/Ca(OH)₂ reversible hydration/dehydration reaction

Transient conservation equations of mass and energy transport including chemical kinetics are solved numerically subject to appropriate boundary and initial conditions to examine the influence of the mass transfer resistance on the overall performance of this CHP configuration. These results are presented and discussed with the aim of enhancing the CHP performance in next generation reactor designs.