Recent progresses in the development of tubular segmented-in-series solid oxide fuel cells:Experimental and numerical study

Solid oxide fuel cells(SOFCs)have attracted a great deal of interest because they have the highest efficiency without using any noble metal as catalysts among all the fuel cell technologies.However,traditional SOFCs suffer from having a higher volume,current leakage,complex connections,and difficulty in gas sealing.To solve these problems,Rolls-Royce has fabricated a simple design by stacking cells in series on an insulating porous support,resulting in the tubular segmented-in-series solid oxide fuel cells(SIS-SOFCs),which achieved higher output voltage.This work systematically reviews recent advances in the structures,preparation methods,perform-ances,and stability of tubular SIS-SOFCs in experimental and numerical studies.Finally,the challenges and future development of tubular SIS-SOFCs are also discussed.The findings of this work can help guide the direction and inspire innovation of future development in this field.

NUMERICAL STUDIES FOR GAS SOLID TWO PHASE STEADY MIXED CONVECTION PROBLEMS WITH PHASE CHANGE

A mathematical model for describing gas solid two phase steady mixed convection with phase change has been developed and numerical calculation methods presented.A melting liquid droplet failing a counter gas currenl expe- riences three processes,cooling of liquid droplet,solidification and cooling of the solid particle.The turbulent model used for Rayleigh number greater than 10~6 is a two equation(k—ε)model of turbulence.For phase change,an improved enthalpy method with varied time step is proposed.The gas particle two phase flow is described by using Eulerian-Lagrangian approach.Modified SIMPLE algorithm and Runge-Kutta method are used in interative calcu- lation.As an example of calculation,the flow in a special 2-dimensional axi-symmetrical prilling tower of diameter 20 m and height 50 m has been performed.Buoyancy effect is important for moving droplet with phase change. The model to be developed and analysis of results obtained in this paper are useful for engineering design in indus- try.

Evaluation of Water Losses by Evaporation in the Nakanbe Basin

A numerical approach to heat and mass transfer in a large water reservoir is presented. This water reservoir is likened to a parallelepiped reservoir whose vertical and lower walls are adiabatic and impermeable. The equations that govern natural convection in water are solved by the finite volume method and Thomas’salgorithm. The adequacy between the velocity and pressure fields is ensured by the SIMPLE algorithm. We are going to evaluate the water losses by evaporation from three dams in the Nakanbé basin in Burkina Faso for a period of thirty years, that is to say from January 1, 1991, to March 15, 2020. The three dams have a rate of evaporation greater than 40% of the volume of water stored. Indeed the rate of evaporation in each dam increases with the water filling rate in the reservoir: we have observed the following results for each dam in the Nakanbé basin;for the date of 02/27/1988 to 03/13/2020., the Loumbila dam received a total volume of stored water of 22.02 Mm3 and 10.57 Mm3 as the total volume of water evaporated at the same date. At the Ouaga dam (2 + 3), it stored a water volume of 4.06 Mm3 and evaporated 2.03 Mm3 of its storage volume from 01/01/1988 to 05/07/2016. Finally, with regard to the Bagré dam, it stored 745.16 Mm3 of water and 365.13 Mm3 as the volume of water evaporated from 01/01/1993 to 03/31/2020.

Tapered Suction Caissons:A Numerical Investigation into Their Pull-out Performance

Suction caissons have been widely used as anchors and foundations for floating and fixed offshore platforms. The pull-out performance of conventional suction caissons （with upright walls） has been investigated by a number of researchers. However, no attention has been paid to tapered suction caissons. This paper deals with the performance of tapered suction caissons under vertical pull-out loads. A numerical approach is used for this purpose. The numerical model is first verified against test data available for common upright caissons. The verified model is then used to study the pullout performance of tapered suction caissons. It is noticed that the pull-out capacities exhibited by tapered suction caissons are in general considerably higher than those from their corresponding traditional upright caissons. To obtain an insight into this superior behaviour, effects from certain soil/caisson/drainage parameters on the pull-out capacity of tapered suction caissons are studied. Soil cohesion is noticed to have a linear improving effect on the capacity of both upright and ta- pered suction caissons. The soil internal friction angle is noticed to have an exponential increasing effect on the pull-out capacity. With a constant caisson diameter, an increase in the aspect ratio is seen to particularly influence the pull-out capacity. With a constant caisson length, an increase in the aspect ratio is discovered to result in non-linear decrease in the pull-out capacity. Under undrained conditions, tapered models generally show less sensitivity to above mentioned soil/caisson parameters as compared with those under drained conditions.

A numerical study on the impact of tidal waves on the storm surge in the north of Liaodong Bay

A storm surge is an abnormal sharp rise or fall in the seawater level produced by the strong wind and low pressure field of an approaching storm system.A storm tide is a water level rise or fall caused by the combined effect of the storm surge and an astronomical tide.The storm surge depends on many factors,such as the tracks of typhoon movement,the intensity of typhoon,the topography of sea area,the amplitude of tidal wave,the period during which the storm surge couples with the tidal wave.When coupling with different parts of a tidal wave,the storm surges caused by a typhoon vary widely.The variation of the storm surges is studied.An once-in-a-century storm surge was caused by Typhoon 7203 at Huludao Port in the north of the Liaodong Bay from July 26th to 27th,1972.The maximum storm surge is about 1.90 m.The wind field and pressure field used in numerical simulations in the research were derived from the historical data of the Typhoon 7203 from July 23rd to 28th,1972.DHI Mike21 is used as the software tools.The whole Bohai Sea is defined as the computational domain.The numerical simulation models are forced with sea levels at water boundaries,that is the tide along the Bohai Straits from July 18th to 29th（2012）.The tide wave and the storm tides caused by the wind field and pressure field mentioned above are calculated in the numerical simulations.The coupling processes of storm surges and tidal waves are simulated in the following way.The first simulation start date and time are 00：00 July 18th,2012; the second simulation start date and time are 03：00 July 18th,2012.There is a three-hour lag between the start date and time of the simulation and that of the former one,the last simulation start date and time are 00：00 July 25th,2012.All the simulations have a same duration of 5 days,which is same as the time length of typhoon data.With the first day and the second day simulation output,which is affected by the initial field,being ignored,only the 3rd to 5th day simulation results are used to study the rules of the storm surges in the north of the Liaodong Bay.In total,57 cases are calculated and analyzed,including the coupling effects between the storm surge and a tidal wave during different tidal durations and on different tidal levels.Based on the results of the 57 numerical examples,the following conclusions are obtained：For the same location,the maximum storm surges are determined by the primary vibration（the storm tide keeps rising quickly） duration and tidal duration.If the primary vibration duration is a part of the flood tidal duration,the maximum storm surge is lower（1.01,1.05 and 1.37 m at the Huludao Port,the Daling Estuary and the Liaohe Estuary respectively）.If the primary vibration duration is a part of the ebb tidal duration,the maximum storm surge is higher（1.92,2.05 and 2.80 m at the Huludao Port,the Daling Estuary and the Liaohe Estuary respectively）.In the mean time,the sea level restrains the growth of storm surges.The hour of the highest storm tide has a margin of error of plus or minus 80 min,comparing the high water hour of the astronomical tide,in the north of the Liaodong Bay.

Study on Numerical Simulation of Mold Filling and Heat Transfer in Die Casting Process

A 3-D mathematical model considering turbulence phenomena has been established based on a computational fluid dynamics technique, so called 3-D SOLA-VOF (Solution Algorithm-Volume of Fluid), to simulate the fluid flow of mold filling process of die casting. In addition, the mathematical model for simulating the heat transfer in die casting process has also been established. The computation program has been developed by the authors with the finite difference method (FDM) recently. As verification, the mold filling process of a S-shaped die casting has been simulated and the simulation results coincide with that of the benchmark test. Finally, as a practical application, the gating design of a motorcycle component was modified by the mold filling simulation and the dies design of another motorcycle component was optimized by the heat transfer simulation. All the optimized designs were verified by the production practice.

Numerical simulation on the mechanism of the normal impact of two droplets onto a thin film

The normal impingement process of two water droplets upon a thin film on the solid surface was numerically investigated. The numerical treatment was based on the finite volume solution of the Navier-Stokes （N-S） equations combined with the volume of fluid method （VOF）. Physically reasonable results for the process of two droplets impacting on the thin film were obtained. The effects of the droplet velocity, the fihn thickness and the spacing between the two droplets on the splash and spread process of the impact were examined.

Numerical study on three-dimensional flow field of continuously rotating detonation in a toroidal chamber

Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.

Numerical Study of A Round Buoyant Jet Under the Effect of JONSWAP Random Waves

This paper presents a numerical study on the hydrodynamic behaviours of a round buoyant jet under the effect of JONSWAP random waves. A three-dimensional large eddy simulation （LES） model is developed to simulate the buoyant jet in a stagnant ambient and JONSWAP random waves. By comparison of velocity and concentration fields, it is found that the buoyant jet exhibits faster decay of centerline velocity, wider lateral spreading and larger initial dilution under the wave effect, indicating that wave dynamics improves the jet entrainment and mixing in the near field, and subsequently mitigate the jet impacts in the far field. The effect of buoyancy force on the jet behaviours in the random waves is also numerically investigated. The results show that the wave effect on the jet entrainment and mixing is considerably weakened under the existence of buoyancy force, resulting in a slower decay rate of centerline velocity and a narrower jet width for the jet with initial buoyancy.

A Numerical Study on the Effect of an Extratropical Cyclone on the Evolution of a Midlatitude Front

The extratropical transition （ET） of tropical cyclone （TC） Haima （2004） was simulated to understand the impact of TC on midlatitude frontal systems. Two experiments were conducted using the Advanced Research version of the Weather Research and Forecast （WRF） model. In the control run （CTL）, a vortex was extracted from the 24-hour pre-run output and then inserted into the National Centers for Environmental Prediction （NCEP） global final （FNL） analysis as an initial condition, while TC circulation was removed from the initial conditions in the sensitivity run （NOTC）. Comparisons of the experiments demonstrate that the midlatitude front has a wider meridional extent in the NOTC run than that in the CTL run. Furthermore, the CTL run produces convection suppression to the southern side of the front due to strong cold advection related to the TC circulation. The easterly flow north of the TC not only decelerates the eastward displacement of the front and contracts its zonal scale but also transports more moisture westward and lifts the air along equivalent potential temperature surfaces ahead of the front. As a result, the ascending motion and diabatic heating are enhanced in the northeastern edge of the front, and the anticyclonic outflow in the upper-level is intensified. The increased pressure gradient and divergent ftow aloft strengthen the upper-level jet and distort the trough axis in a northwest-southeast orientation. The thermal contrast between the two systems and the dynamic contribution related to the TC circulation can facilitate scalar and rotational frontogenesis to modulate the frontal structure.

A Study on Numerical Simulation of Core-Shooting Process

In order to evaluate the main factors influencing the core-shooting process and to optimize the design of core boxes, the fluid-particle (air-sand) model has been built based on the two-phase flow theory. The fluid phase, air, and the particulate phase, sand granules, have been treated as a continuum. By using this model, it is possible to simulate the flow and compaction behavior of sand particles during the core-shooting process. To benchmark the calculated results, the shooting process has been recorded with a digital high speed camera, and the inlet condition of sand particles has also been achieved by using the camera. The preliminary results have showed that the calculation is in agreement with the testing results.

Numerical Study of Fluid Dynamics and Heat Transfer Induced by Plasma Discharges

A numerical investigation is conducted to explore the evolution of a plasma discharge and its interaction with the fluid flow based on a self-consistent fluid model which couples the discharge dynamics with the fluid dynamics.The effects of the applied voltage on the distribution of velocity and temperature in initially static air are parainetrically studied.Furthermore,the spatial structure of plasma discharge and the resulting force contours in streamwise and normal directions are discussed in detail.The result shows that the plasma actuator produces a net force that should always be directed away from the exposed electrode,which results in an ionic wind pushing particles into a jet downstream of the actuator.When the energy added by the plasma is taken into account,the ambient air temperature is increased slightly around the electrode,but the velocity is almost not affected.Therefore it is unlikely that the induced flow is buoyancy driven.For the operating voltages considered in this paper,the maximum induced velocity is found to follow a power law,i.e.,it is proportional to the applied voltage to the 3.5 power.This promises an efficient application in the flow control with plasma actuators.

A Numerical Study on Forecasting the Henan Extraordinarily Heavy Rainfall Event in August 1975

This study is essentially an experiment on the control experiment in the August 1975 catastrophe which was the heaviest rainfall in China's Mainland with a maximum 24-h rainfall of 1060.3 mm, and it significantly demonstrates that the limited area model can still skillfully give reasonable results even only the conventional data are available. For such a heavy rainfall event, a grid length of 90 km is too large while 45 km seems acceptable. Under these two grid sizes, the cumulus parameterization scheme is evidently superior to the explicit scheme since it restricts instabilities such as CISK to limited extent. The high resolution scheme for the boundary treatment does not improve forecasts significantly.The experiments also revealed some interesting phenomena such as the forecast rainfall being too small while affecting synoptic system so deep as compared with observations. Another example is the severe deformation of synoptic systems both in initial conditions and forecast fields in the presence of complicated topography. Besides, the fixed boundary condition utilized in the experiments along with current domain coverage set some limitations to the model performances.

Numerical and Experimental Study of the 3D Effect on Connecting Arm of Vertical Axis Tidal Current Turbine

Vertical axis tidal current turbine is a promising device to extract energy from ocean current. One of the important components of the turbine is the connecting arm, which can bring about a significant effect on the pressure distribution along the span of the turbine blade, herein we call it 3D effect. However, so far the effect is rarely reported in the research, moreover, in numerical simulation. In the present study, a 3D numerical model of the turbine with the connecting arm was developed by using FLUENT software compiling the UDF（User Defined Function） command. The simulation results show that the pressure distribution along the span of blade with the connecting arm model is significantly different from those without the connecting arm. To facilitate the validation of numerical model, the laboratory experiment has been carried out by using three different types of NACA aerofoil connecting arm and circle section connecting arm. And results show that the turbine with NACA0012 connecting arm has the best start-up performance which is 0.346 m/s and the peak point of power conversion coefficient is around 0.33. A further study has been performed and a conclusion is drawn that the aerofoil and thickness of connecting arm are the most important factors on the power conversion coefficient of the vertical axis tidal current turbine.

Natural Convection in a Vertical Open-Ended Channel： Comparison between Experimental and Numerical Results

The present study deals with natural convection flow in a vertical open-ended channel with wall constant heat flux. The experimental and numerical investigations are both conducted using water as the working fluid. The numerical code is developed using finite differences scheme to solve the Navier-Stokes equations under the Boussinesq assumption. Concerning the experimental apparatus, it consists of two heated walls immersed in water. Temperature and velocity measurements are provided for different modified Rayleigh numbers based on the walls spacing b Rab = 1.67 x 10 6,3.6 x 10 6,8.97 x10 6,1.69 x 10 7,4.29 x 10 7. The numerical code is first validated with a numerical benchmark. Then, comparison between experimental and numerical results is performed. The code provides a satisfactory prediction of main quantities compared to the experimental results but only for the lowest Rayleigh numbers. For higher modified Rayleigh numbers, the flow becomes three-dimensional and turbulent. Therefore, 2D numerical simulations fail to predict flow and heat transfer for this range of modified Rayleigh number.

Numerical Study of Natural Convection in a Two-Dimensional Enclosure with a Sinusoidal Boundary Thermal Condition Utilizing Nanofluid

Nanofluids are considered to offer important advantages over conventional heat transfer fluids. A model is developed to analyze the behavior of nanofluids taking into account the solid fraction χ. The Navier-Stokes equations are solved numerically with Alternating Direct Implicit method (ADI method) for various Grashof numbers 104 and 105;we have an excellent agreement between our numerical code and previously published works. Copper-Water nanofluid is used with Pr = 6.2 and solid volume fraction χ is varied as 0%;5%;10%;15% and 20%. The problem considered is a two-dimensional heat transfer in a square cavity. The vertical walls are differentially heated, the left is maintained at hot con- dition (sinusoidal) when the right one is cold. The horizontal walls are assumed to be insulated, non conducting and impermeable to mass transfer. The nanofluid in the enclosure is Newtonian, incompressible and laminar. The nanopar- ticles are assumed to have a uniform shape and size. Moreover, it is assumed that both the fluid phase and nanoparticles are in thermal equilibrium state and they flow at the same velocity. The thermophysical properties of the nanofluid are assumed to be constant except for the density variation in the buoyancy force, which is based on the Boussinesq approximation. Different correlations are proposed for predicting heat transfer for uniform and sinusoidal boundary thermal conditions.

Numerical Study on the Drag Coefficients of Sphere and Double Spheres

Winds on the earth are commonly strong enough to erode transport and deposit sediment. The modes of sand transport by the wind are greatly different from those by water flow. On the other hand wind-blown sands are of a material circulation process of the earth surface. They affect wind-sand transport flux and sand ejection of a flux, the damage of grains formed cannot be neglected in engineering. Because of the complexity of windblown sand flux system, the understanding of its basic mechanics is not yet clear. The key forces in sand salutation mainly includes： the valid gravity, air drag force ＇Magnus force＇ Saffman force ＇Basset force＇ additional quality force and scatter force among grains. The most important force in sand salutation is the air drag force. Computation of the single sphere drag coefficient and double spheres drag coefficient is presented for the distance between two spheres being smaller than twelve times of the sphere diameter and the spheres being at different angular positions. The flow interference of two spheres was investigated for the distance s = 0.08 d to 12d, angular position 0 = 0 to 360 and Reynolds number 15≤Re≤1000.

A study on the numerical prediction method for the vertical thermal structure in the Bohai Sea and the Huanghai Sea-I.One-dimensional numerical prediction model

In this paper, on the basis of the heat conduction equation without consideration of the advection and turbulence effects, one-dimensional model for describing surface sea temperature ( T1), bottom sea temperature ( Tt ) and the thickness of the upper homogeneous layer ( h ) is developed in terms of the dimensionless temperature θT and depth η and self-simulation function θT - f(η) of vertical temperature profile by means of historical temperature data.The results of trial prediction with our one-dimensional model on T, Th, h , the thickness and gradient of thermocline are satisfactory to some extent.

Numerical study on residual current and its impact on mass transport in the Hangzhou Bay and the Changjiang Estuary Ⅱ. Residual current and its impact on mass transport in winter

Some observational characteristics of residual current and mass transport in the Hangzhou Bay and the Changjiang Estuary in winter are analyzed. The residual current and its impact on mass transport are simulated with a 3 - D joint model for the Hangzhou Bay and the Changjiang Estuary, in which the impacts of river flux, wind, baroclinic pressure gradient (BPG), background current in the East China Sea and tide (including M2, S2, K1 and O1) are taken into account. Based on there studies, further simulations are made to analyze the dynamical mechanisms of the observational characteristics.

Numerical Study on the Energy Extraction Performance of Coupled Tandem Flapping Hydrofoils

Tidal current energy is a promising renewable energy source for future electricity supply.The flapping hydrofoil is regarded as a useful tool to extract the tidal current energy in shallow water.A concept of coupled tandem flapping hydrofoils under semi-activated mode was proposed in the present study.A two-dimensional numerical model,based on the computational fluid dynamics software ANSYS-Fluent,was established to investigate the power extraction performance of the coupled tandem flapping hydrofoils.The effects of the reduced frequency,pitching amplitude,moment of inertia,damping coefficient,and longitudinal distance between hydrofoils were studied.The vortices,pressure distribution,and kinetic characteristics of hydrofoils under various conditions were analyzed to reveal the interaction between the shedding vortex and hydrofoils.The energy extraction mechanism and hydrodynamic performance were analyzed.The positive interactions for energy harvesting were identified for improvements of the further performance.The peak values of efficiency and power coefficient were achieved at 0.69 and 2.13,respectively.