Effect of bipolar-plates design on corrosion,mass and heat transfer in proton-exchange membrane fuel cells and water electrolyzers:A review

Attaining a decarbonized and sustainable energy system,which is the core solution to global energy issues,is accessible through the development of hydrogen energy.Proton-exchange membrane water electrolyzers(PEMWEs)are promising devices for hydrogen production,given their high efficiency,rapid responsiveness,and compactness.Bipolar plates account for a relatively high percentage of the total cost and weight compared with other components of PEMWEs.Thus,optimization of their design may accelerate the promotion of PEMWEs.This paper reviews the advances in materials and flow-field design for bipolar plates.First,the working conditions of proton-exchange membrane fuel cells(PEMFCs)and PEMWEs are compared,including reaction direction,operating temperature,pressure,input/output,and potential.Then,the current research status of bipolar-plate substrates and surface coatings is summarized,and some typical channel-rib flow fields and porous flow fields are presented.Furthermore,the effects of materials on mass and heat transfer and the possibility of reducing corrosion by improving the flow field structure are explored.Finally,this review discusses the potential directions of the development of bipolar-plate design,including material fabrication,flow-field geometry optimization using threedimensional printing,and surface-coating composition optimization based on computational materials science.

Turbine Passage Secondary Flow Dynamics and Endwall Heat Transfer Under Different Inflow Turbulence

This study presents endwall hydrodynamics and heat transfer in a linear turbine cascade at Re 5×105 at low and high intensities of turbulence.Results are numerically predicted using the standard SST model and Reθ-γtransition model as well as using the high-resolution LES separately.The major secondary flow components,comprising the horseshoe,corner,and passage vortices are recognized and the impact on heat or mass transfer is investigated.The complicated behavior of turbine passage secondary flow generation and establishment are impacted by the perspective of boundary layer attributes and inflow turbulence.The passage vortex concerning the latest big leading-edge vane is generated by the enlargement of the circulation developed at the first instance adjacent to the pressure side becomes powerful and mixes with other vortex systems during its migration towards the suction side.The study conclusions reveal that substantial enhancements are attained on the endwall surface,for the entire spanwise blade extension on the pressure surface,and in the highly 3-D region close to the endwall on the suction surface.The forecasted suction surface thermal exchange depicts great conformity with the measurement values and precisely reproduces the enhanced thermal exchange owing to the development and lateral distribution of the secondary flows along the midspan of the blade passage downstream.The impacts of the different secondary flow structures on the endwall thermal exchange are described in depth.

Numerical modeling and parametric sensitivity analysis of heat transfer and two-phase oil and water flow characteristics in horizontal and inclined flowlines using OpenFOAM

Estimating the oil-water temperatures in flowlines is challenging especially in deepwater and ultra-deepwater offshore applications where issues of flow assurance and dramatic heat transfer are likely to occur due to the temperature difference between the fluids and the surroundings. Heat transfer analysis is very important for the prediction and prevention of deposits in oil and water flowlines, which could impede the flow and give rise to huge financial losses. Therefore, a 3D mathematical model of oil-water Newtonian flow under non-isothermal conditions is established to explore the complex mechanisms of the two-phase oil-water transportation and heat transfer in different flowline inclinations. In this work, a non-isothermal two-phase flow model is first modified and then implemented in the InterFoam solver by introducing the energy equation using OpenFOAM® code. The Low Reynolds Number (LRN) k-ε turbulence model is utilized to resolve the turbulence phenomena within the oil and water mixtures. The flow patterns and the local heat transfer coefficients (HTC) for two-phase oil-water flow at different flowlines inclinations (0°, +4°, +7°) are validated by the experimental literature results and the relative errors are also compared. Global sensitivity analysis is then conducted to determine the effect of the different parameters on the performance of the produced two-phase hydrocarbon systems for effective subsea fluid transportation. Thereafter, HTC and flow patterns for oil-water flows at downward inclinations of 4°, and 7° can be predicted by the models. The velocity distribution, pressure gradient, liquid holdup, and temperature variation at the flowline cross-sections are simulated and analyzed in detail. Consequently, the numerical model can be generally applied to compute the global properties of the fluid and other operating parameters that are beneficial in the management of two-phase oil-water transportation.

Analysis of periodic pulsating nanofluid flow and heat transfer through a parallel-plate channel in the presence of magnetic field

In this paper,we focus on the two-dimensional pulsating nanofluid flow through a parallel-plate channel in the presence of a magnetic field.The pulsating flow is produced by an applied pressure gradient that fluctuates with a small amplitude.A kind of proper transformation is used so that the governing equations describing the momentum and thermal energy are reduced to a set of non-dimensional equations.The analytical expressions of the pulsating velocity,temperature,and Nusselt number of nanofluids are obtained by the perturbation technique.In the present study,the effects of the Cu-H2O and Al_(2)O_(3)-H2O nanofluids on the flow and heat transfer in pulsating flow are compared and analyzed.The results show that the convective heat transfer effect of Cu-H2O nanofluids is better than that of Al_(2)O_(3)-H2O nanofluids.Also,the effects of the Hartmann number and pulsation amplitude on the velocity,temperature,and Nusselt number are examined and discussed in detail.The present work indicates that increasing the Hartmann number and pulsation amplitude can enhance the heat transfer of the pulsating flow.In addition,selecting an optimal pulsation frequency can maximize the convective heat transfer of the pulsating flow.Therefore,improved understanding of these fundamental mechanisms is conducive to the optimal design of thermal systems.

Effects of tube cross-sectional shapes on flow pattern, liquid film and heat transfer of n-pentane across tube bundles

The heat transfer of hydrocarbon refrigerant across tube bundles have been widely used in refrigeration.Three-dimensional simulation model using volume of fluid(VOF) was presented to study the effects of tube shapes on flow pattern, film thickness and heat transfer of n-pentane across tube bundles, including circle, ellipse-shaped, egg-shaped and cam-shaped tube bundles. Simulation results agree well with experimental data in the literature. The liquid film thickness of sheet flow and heat transfer for different tube shapes were obtained numerically. The flow pattern transition occurs lower vapor quality for ellipse-shaped tube than other tube shapes. For sheet flow, the liquid film on circle tube and ellipseshaped tube is symmetrically distributed along the circumferential direction. However, the liquid film on egg-shaped tube at circumferential angles(θ) = 15°–60° is thicker than θ = 135°–165°. The liquid film on cam tube at θ = 15°–60° is slightly thinner than θ = 135°–165°. The liquid film thickness varies from thinner to thicker for ellipse-shaped, cam-shaped, egg-shape and circle within θ = 15°–60°. The effect of tube shape is insignificant on thin liquid film thickness. Ellipse-shaped tube has largest heat transfer coefficient for sheet flow. In practical engineering, the tube shape could be designed as ellipse to promote heat transfer.

Research on shell-side heat and mass transfer with multi-component in LNG spiral-wound heat exchanger under sloshing conditions

The spiral-wound heat exchanger(SWHE) is the primary low-temperature heat exchanger for large-scale LNG plants due to its high-pressure resistance, compact structure, and high heat exchange efficiency. This paper studied the shell-side heat and mass transfer characteristics of vapor-liquid two-phase mixed refrigerants in an SWHE by combining a multi-component model in FLUENT software with a customized multicomponent mass transfer model. Besides, the mathematical model under the sloshing condition was obtained through mathematical derivation, and the corresponding UDF code was loaded into FLUENT as the momentum source term. The results under the sloshing conditions were compared with the relevant parameters under the steady-state condition. The shell-side heat and mass transfer characteristics of the SWHE were investigated by adjusting the component ratio and other working conditions. It was found that the sloshing conditions enhance the heat transfer performance and sometimes have insignificant effects. The sloshing condition is beneficial to reduce the flow resistance. The comprehensive performance of multi-component refrigerants has been improved and the improvement is more significant under sloshing conditions, considering both the heat transfer and pressure drop.These results will provide theoretical support for the research and design of multi-component heat and mass transfer enhancement of LNG SWHE under ocean sloshing conditions.

A new model of flow over stretching(shrinking)and porous sheet with its numerical solutions

The viscous fluid flow and heat transfer over a stretching(shrinking)and porous sheets of nonuniform thickness are investigated in this paper.The modeled problem is presented by utilizing the stretching(shrinking)and porous velocities and variable thickness of the sheet and they are combined in a relation.Consequently,the new problem reproduces the different available forms of flow motion and heat transfer maintained over a stretching(shrinking)and porous sheet of variable thickness in one go.As a result,the governing equations are embedded in several parameters which can be transformed into classical cases of stretched(shrunk)flows over porous sheets.A set of general,unusual and new variables is formed to simplify the governing partial differential equations and boundary conditions.The final equations are compared with the classical models to get the validity of the current simulations and they are exactly matched with each other for different choices of parameters of the current problem when their values are properly adjusted and manipulated.Moreover,we have recovered the classical results for special and appropriate values of the parameters(δ_(1),δ_(2),δ_(3),c,and B).The individual and combined effects of all inputs from the boundary are seen on flow and heat transfer properties with the help of a numerical method and the results are compared with classical solutions in special cases.It is noteworthy that the problem describes and enhances the behavior of all field quantities in view of the governing parameters.Numerical result shows that the dual solutions can be found for different possible values of the shrinking parameter.A stability analysis is accomplished and apprehended in order to establish a criterion for the determinations of linearly stable and physically compatible solutions.The significant features and diversity of the modeled equations are scrutinized by recovering the previous problems of fluid flow and heat transfer from a uniformly heated sheet of variable(uniform)thickness with variable(uniform)stretching/shrinking and injection/suction velocities.

Stability Scrutinization of Agrawal Axisymmetric Flow of Nanofluid through a Permeable Moving Disk Due to Renewable Solar Radiation with Smoluchowski Temperature and Maxwell Velocity Slip Boundary Conditions

The utilization of solar energy is essential to all living things since the beginning of time.In addition to being a constant source of energy,solar energy(SE)can also be used to generate heat and electricity.Recent technology enables to convert the solar energy into electricity by using thermal solar heat.Solar energy is perhaps the most easily accessible and plentiful source of sustainable energy.Copper-based nanofluid has been considered as a method to improve solar collector performance by absorbing incoming solar energy directly.The goal of this research is to explore theoretically the Agrawal axisymmetric flow induced by Cu-water nanofluid over a moving permeable disk caused by solar energy.Moreover,the impacts of Maxwell velocity and Smoluchowski temperature slip are incorporated to discuss the fine points of nanofluid flow and characteristics of heat transfer.The primary partial differential equations are transformed to similarity equations by employing similarity variables and then utilizing bvp4c to resolve the set of equations numerically.The current numerical approach can produce double solutions by providing suitable initial guesses.In addition,the results revealed that the impact of solar collector efficiency enhances significantly due to nanoparticle volume fraction.The suction parameter delays the boundary layer separation.Moreover,stability analysis is performed and is found that the upper solution is stable and physically trustworthy while the lower one is unstable.

Fluid Flow and Heat Transfer Characteristic of Outer and Inner Half Coil Jackets

The physical models of the outer and inner half coil jackets were simplified to two types of coiled ducts.The mathematic models of incompressible fluid at the condition of laminar flow and heat transfer in the two types of jackets for cooling process reactor were set up and solved by the semi-implicit method for pressure linked equa-tions consistent （SIMPLEC） algorithm based on a control volume method.The flow and temperature fields were given and the effects of Dean and Prandtl numbers on flow and heat transfer were studied.The results show that flow in the inner half coil jacket is found to exhibit transition of secondary flow pattern from two vortices to four vortices when the Dean number increases,but that in the outer half coil jacket is not found.The critical Dean num-ber is about 96.The inner half coil jacket has stronger heat transfer ability than the outer half coil jacket and this superiority is more evident with larger Prandtl number.However,as the Dean number is greater than 105,the flow resistance enhances more severely in the inner jacket than the outer jacket.For both jackets,the centers of the heated wall are the poorest for heat transfer.

NUMERICAL ANALYSIS OF THREE-DIMENSIONAL FLUID FLOW AND HEAT TRANSFER IN TIG WELD POOL WITH FULL PENETRATION

A model is established to analyze three-dimensional fluid flow and heat transfer in TICweld pools with full penetration.It considers the deformation of the molten pool surfaceat the condition of full penetrated workpieees,takes the are pressure as the drivingforce of the pool surface deformation,and determines the surface configuration of weldpool based on the dynamic balance of arc pressure,pool gravity and surface tension atdeformed weld pool surface. The SIMPLER algorithm is used to calculate the fluid flowfield and temperature distribution in TIG weld pools of stainless steel workpieces.TIGwelding experiments are made to verify the validity of the model.It shows the calculatedresults by the model are in good agreement with experimental measurements. professor,Dept of Welding Engineering,Harbin Institute of Technology,Harbin 150006,China

Numerical investigation of laminar heat transfer and nanofluid flow between two porous horizontal concentric cylinders

In this study, the laminar heat transfer and nanofluid flow between two porous horizontal concentric cylinders was investigated. The problem is investigated in two different geometries and the Re=10, 25, 50, 75, 100 and volume fraction 0, 0.2%, 0.5%, 2% and 5% that related to copper nanoparticles, and porous medium porosity of 0.5 and 0.9. Compared to the first geometry, the convective coefficient in the second geometry increases by 8.3%, 7% and 5.5% at Reynolds numbers of 100, 75 and 50, respectively. Comparison of the outlet temperatures for two heat fluxes of 300 and 1200 W/m^2 indicates a 2.5% temperature growth by a fourfold increase in the heat fluxes. Also, the higher Nusselt number is associated with the second geometry occurring at porosities of 0.9 and 0.5, respectively. In both geometries, the Nusselt number values at the porosity of 0.9 are higher, which is due to the increased nanofluid convection at higher porosities. The velocity of the nanofluid experiences a two-fold increase at the outlet compared to its inlet velocity in the first geometry and for both porosities. Similarly, a three-fold increase was achieved in the second geometry and for both porosities.

Modelling Study to Compare the Flow and Heat Transfer Characteristics of Low-Power Hydrogen,Nitrogen and Argon Arc-Heated Thrusters

A modelling study is performed to compare the plasma flow and heat transfer characteristics of low-power arc-heated thrusters （arcjets） for three different propellants： hydrogen, nitrogen and argon. The all-speed SIMPLE algorithm is employed to solve the governing equations, which take into account the effects of compressibility, Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. The temperature, velocity and Mach number distributions calculated within the thruster nozzle obtained with different propellant gases are compared for the same thruster structure, dimensions, inlet-gas stagnant pressure and arc currents. The temperature distributions in the solid region of the anode-nozzle wall are also given. It is found that the flow and energy conversion processes in the thruster nozzle show many similar features for all three propellants. For example, the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appearing near the cathode tip; the flow transition from the subsonic to supersonic regime occurs within the constrictor region; the highest axial velocity appears inside the nozzle; and most of the input propellant flows towards the thruster exit through the cooler gas region near the anode-nozzle wall. However, since the properties of hydrogen, nitrogen and argon, especially their molecular weights, specific enthMpies and thermal conductivities, are different, there are appreciable differences in arcjet performance. For example, compared to the other two propellants, the hydrogen arcjet thruster shows a higher plasma temperature in the arc region, and higher axial velocity but lower temperature at the thruster exit. Correspondingly, the hydrogen arcjet thruster has the highest specific impulse and arc voltage for the same inlet stagnant pressure and arc current. The predictions of the modelling are compared favourably with available experimental results.

A New Experimental Rig of Testing Flow Boiling Heat Transfer of Refrigerantand Lubricant Mixture

This paper proposed a new experimental rig of testing flow boiling heat transfer of refrigerant and lubricant oil mixture. The quantity of oil in the test section can be controlled and regulated conveniently and accurately by connecting separate lubricant oil circuit with test section in parallel. It was built up by retrofitting a multiple air-conditioner and installing three oil-separators in serials at the compressor outlet. And so the lubricant oil in the discharged refrigerant gas of compressor can be removed completely.The refrigerant flow rate through test section can be bypassed by the by-path circuit of indoor unit.This experimental rig has advantages such as on-line and continuous oil injection, short time of obtaining stability, flexible operation, simple control, which lead to high efficiency in the research of flow boiling heat transfer of refrigerant and lubricant oil mixture.

Effect of viscous dissipation and heat source on flow and heat transfer of dusty fluid over unsteady stretching sheet

This paper investigates the problem of hydrodynamic boundary layer flow and heat transfer of a dusty fluid over an unsteady stretching surface. The study considers the effects of frictional heating （viscous dissipation） and internal heat generation or ab- sorption. The basic equations governing the flow and heat transfer are reduced to a set of non-linear ordinary differential equations by applying suitable similarity transformations. The transformed equations are numerically solved by the Runge-Kutta-Fehlberg-45 order method. An analysis is carried out for two different cases of heating processes, namely, variable wall temperature （VWT） and variable heat flux （VHF）. The effects of various physical parameters such as the magnetic parameter, the fluid-particle interaction pa- rameter, the unsteady parameter, the Prandtl number, the Eckert number, the number density of dust particles, and the heat source/sink parameter on velocity and temperature profiles are shown in several plots. The effects of the wall temperature gradient function and the wall temperature function are tabulated and discussed.

MHD flow and heat transfer of micropolar fluid between two porous disks

A numerical study is carried out for the axisymmetric steady laminar incompressible flow of an electrically conducting micropolar fluid between two infinite parallel porous disks with the constant uniform injection through the surface of the disks. The fluid is subjected to an external transverse magnetic field. The governing nonlinear equations of motion are transformed into a dimensionless form through yon Karman＇s similarity transformation. An algorithm based on a finite difference scheme is used to solve the reduced coupled ordinary differential equations under associated boundary conditions. The effects of the Reynolds number, the magnetic parameter, the micropolar parameter, and the Prandtl number on the flow velocity and temperature distributions are discussed. The results agree well with those of the previously published work for special cases. The investigation predicts that the heat transfer rate at the surfaces of the disks increases with the increases in the Reynolds number, the magnetic parameter, and the Prandtl number. The shear stresses decrease with the increase in the injection while increase with the increase in the applied magnetic field. The shear stress factor is lower for micropolar fluids than for Newtonian fluids, which may be beneficial in the flow and thermal control in the polymeric processing.

Large-eddy Simulation of Fluid Flow and Heat Transfer in a Mixing Tee Junction

The temperature fluctuation caused by thermal striping phenomena of hot and cold fluids mixing results in cyclical thermal stress fatigue failure of the pipe wall. Mean temperature difference between hot and cold fluids was often used as thermal load in previous analysis of thermal fatigue failure, thereby the influences of the amplitude and frequency of temperature fluctuation on thermal fatigue failure were neglected. Based on the mechanism of flow and heat transfer which induces thermal fatigue, the turbulent mixing of hot and cold water in a tee junction is simulated with FLUENT platform by using the Large-eddy simulation（LES） turbulent flow model with the sub-grid scale（SGS） model of Smagorinsky-Lilly（SL） to capture the amplitude and frequency of temperature fluctuation. In a simulation case, hot water with temperature of 343.48 K and velocity of 0.15 m/s enters the horizontal main duct with the side length of 100 mm, while cold water with temperature of 296.78 K and velocity of 0.3 m/s enters the vertical branch duct with the side length of 50 mm. The numerical results show that the mean and fluctuating temperatures are in good agreement with the previous experimental data, which describes numerical simulation with high reliability and accuracy; the power spectrum density（PSD） on top wall is higher than that on bottom wall（as the frequency less than 1 Hz）, while the PSD on bottom wall is relatively higher than that on top wall （as the frequency of 1-10Hz）. The temperature fluctuations in full mixing region of the tee junction can be accurately captured by LES and can provide the theoretical basis for the thermal stress and thermal fatigue analyses.

Heat flow pattern,base of methane hydrates stability zones and BSRs in Shenhu Area,northern South China Sea

Using the collected 433 heat flow values, we estimated the bases of methane hydrate stability zone （BHSZ）, in northern South China Sea （NSCS）. Through comparing BHSZs with the depths of bottom simulating reflectors （BSRs）, in Shenhu Area （SA）, we found that there are big differences between them. In the north of SA, where the water depth is shallow, many slumps developed and the sedimentation rate is high, it appears great negative difference （as large as -192%）. However, to the southeast of SA, where the water depth is deeper, sedimentation rate is relatively low and uplift basement topography exists, it changes to positive difference （as large as ＋45%）. The differences change so great, which haven＇t been observed in other places of the world. After considering the errors from the process of heat flow measurement, the BSR depth, the relationship of thermal conductivity with the sediments depth, and the fluid flow activities, we conclude that the difference should be not caused by these errors. Such big disagreement may be due to the misunderstanding of BSR. The deviant ＂BSRs＂ could represent the paleo-BSRs or just gas-bearing sediment layers, such as unconformities or the specific strata where have different permeability, which are not hydraterelated BSRs.

Laminar Flow and Heat Transfer Characteristics in Jackets of Triangular Flow Channels

Laminar flow and heat transfer characteristics of jacketed vessel with triangular flow channels were numerically studied under hydrodynamically and thermally fully developed conditions. Constant heat flux at theheated wall was assumed. The numerical program code interms of vorticity, stream function, axial velocity com ponent and energy equations was written based on a finite volume method. Based on the numerical results, the flow and temperature field were given, and the effects of Dean and Prandtl numbers on flow and heat transfer were ex amined, and the correlations of flow resistance and mean Nusselt number were developed for the jacket. The results show that the structure of secondary flow is steady two vortices in the investigated range of dimensionless curvatureratio and Reynolds number. Two peaks of local Nusselt number increase significantly with Prandtl and Dean num ber increasing, but the local Nusselt numbers near two ends and at the center of the heated wall increase only slightly. The center and two ends of heated wall are the poor positions for heat transfer in the jacket. Compared with the outer half coil jacket at the same area of heated wall, curvature radius, Reynolds number and Prandtl number, e jacket of triangular flow chmnel has lower flow resistance and less mean Nusselt number.

Dual solutions of MHD stagnation point flow and heat transfer over a stretching/shrinking sheet with generalized slip condition

An analysis was made to study the steady momentum and heat transfer characteristics of a viscous electrically conducting fluid near a stagnation point due to a stretching/shrinking sheet in the presence of a transverse magnetic field and generalized slip condition. Two flow problems corresponding to the planar and axisymmetric stretching/shrinking sheet were considered. By means of similarity transformations, the obtained resultant nonlinear ordinary differential equations were solved numerically using a shooting method for dual solutions of velocity and temperature profiles. Some important physical features of the flow and heat transfer in terms of the fluid velocity, the temperature distribution, the skin friction coefficient and the local Nusselt number for various values of the controlling governing parameters like velocity slip parameter, critical shear rate, magnetic field, ratio of stretching/shrinking rate to external flow rate and Prandtl number were analyzed and discussed. An increase of the critical shear rate decreases the fluid velocity whereas the local Nusselt number increases. The comparison of the present numerical results with the existing literature in a limiting case is given and found to be in an excellent agreement.

Second-order slip MHD flow and heat transfer of nanofluids with thermal radiation and chemical reaction

The effects of the second-order velocity slip and temperature jump boundary conditions on the magnetohydrodynamic （MHD） flow and heat transfer in the presence of nanoparticle fractions are investigated. In the modeling of the water-based nanofluids containing Cu and A1203, the effects of the Brownian motion, thermophoresis, and thermal radiation are considered. The governing boundary layer equations are transformed into a system of nonlinear differential equations, and the analytical approximations of the solutions axe derived by the homotopy analysis method （HAM）. The reliability and efficiency of the HAM solutions are verified by the residual errors and the numerical results in the literature. Moreover, the effects of the physical factors on the flow and heat transfer are discussed graphically.