The Influence of Electrolyte Flow Hydrodynamics on the Performance of a Microfluidic Dye-Sensitized Solar Cell

The dye-sensitized solar cells microfluidically integrated with a redox flow battery (µDSSC-RFB) belong to a new emerging class of green energy sources with an inherent opportunity for energy storage. The successful engineering of microfluidically linked systems is, however, a challenging subject, as the hydrodynamics of electrolyte flow influences the electron and species transport in the system in several ways. In the article, we have analyzed the microflows hydrodynamics by means of the lattice-Boltzmann method, using the algebraic solution of the Navier-Stokes equation for a duct flow and experimentally by the micro particle image velocimetry method. Several prototypes of µDSSC were prepared and tested under different flow conditions. The efficiency of serpentine µDSSC raised from 2.8% for stationary conditions to 3.1% for electrolyte flow above 20 mL/h, while the fill factor increased about 13% and open-circuit voltage from an initial 0.715 V to 0.745 V. Although the hexagonal or circular configurations are the straightforward extensions of standard photo chambers of solar cells, those configurations are hydrodynamically less predictable and unfavorable due to large velocity gradients. The serpentine channel configuration with silver fingers would allow for the scaling of the µDSSC-RFB systems to the industrial scale without loss of performance. Furthermore, the deterioration of cell performance over time can be inhibited by the periodic sensitizer regeneration, which is the inherent advantage of µDSSC.

Direct numerical simulation of turbulent duct flow with inclined or V-shaped ribs mounted on one wall

In this research, highly disturbed turbulent flow of distinct three-dimensional characteristics in a square duct with inclined or V-shaped ribs mounted on one wall is investigated using direct numerical simulation. The turbulence field is highly sensitive to not only the rib geometry but also the boundary layers developed over the side and top walls. In a cross-stream plane secondary flows appear as large longitudinal vortices in both inclined and V-shaped rib cases due to the confinement of four sidewalls of the square duct. However, owing to the difference in the pattern of cross-stream secondary flow motions, the flow physics is significantly different in these two ribbed duct cases. It is observed that the mean flow structures in the cross-stream directions are asymmetrical in the inclined rib case but symmetrical in the V-shaped rib case, causing substantial differences in the momentum transfer across the spanwise direction. The impacts of rib geometry on near-wall turbulence structures are investigated using vortex identifiers, joint probability density functions between the streamwise and vertical velocity fluctuations, statistical moments of different orders, spatial two-point autocorrelations and velocity spectra. It is found that near the leeward and windward rib faces, the mean inclination angle of turbulence structures in the V-shaped rib case is greater than that of the inclined rib case, which subsequently enhances momentum transport between the ribbed bottom wall and the smooth top wall.

The effect of Reynolds number on near-wall reverse flow in a turbulent duct flow

The influence of the Reynolds number on the statistics of a near-wall reverse flow phenomenon, taking place in a turbulent duct flow, is studied. An increase in the NWRF probability is found in both the core and corner regions of the duct walls for higher Reynolds number. The mechanism of the NWRF formation, described recently by Zaripov et al. [1, 2], is validated for higher Reynolds number flows.

Magnetoconvection in a horizontal duct flow at very high Hartmann and Grashof numbers

Direct numerical simulations and linear stability analysis are carried out to study mixed convection in a horizontal duct with constant-rate heating applied at the bottom and an imposed transverse horizontal magnetic field. A two-dimensional approximation corresponding to the asymptotic limit of a very strong magnetic field effect is validated and applied, together with full three-dimensional analysis, to investigate the flow's behaviour in the previously unexplored range of control parameters corresponding to typical conditions of a liquid metal blanket of a nuclear fusion reactor (Hartmann numbers up to $10^4$ and Grashof numbers up to $10^{10}$ ). It is found that the instability to quasi-two-dimensional rolls parallel to the magnetic field discovered at smaller Hartmann and Grashof numbers in earlier studies also occurs in this parameter range. Transport of the rolls by the mean flow leads to magnetoconvective temperature fluctuations of exceptionally high amplitudes. It is also demonstrated that quasi-two-dimensional structure of flows at very high Hartmann numbers does not guarantee accuracy of the classical two-dimensional approximation. The accuracy deteriorates at the highest Grashof numbers considered in the study.

Estimation of transient boundary heat flux by using hybrid algorithm in a two dimensional laminar duct flow

Conjugate gradient method (CGM) with adjoint problem is a popular optimization algorithm in solving the inverse heat transfer problems. It starts with the initial guess solution of unknown parameters to be estimated which would be updated in an iterative procedure. The initial guess solution is one of the significant factors for the accuracy of estimation. In the current study, the Jaya algorithm has been developed to provide the initial guess solution to CGM. The resultant algorithm is called hybrid algorithm. The test problem considered here for the study is the estimation of transient boundary heat flux for two-dimensional hydrodynamically developed and thermally developing forced convective laminar duct flow. The hybrid algorithm is found to be robust and accurate than CGM.

Pattern Transition on Inertial Focusing of Neutrally Buoyant Particles Suspended in Rectangular Duct Flows

The continuous separation and filtration of particles immersed in fluid flows are important interests in various applications. Although the inertial focusing of particles suspended in a duct flow is promising in microfluidics, predicting the focusing positions depending on the parameters, such as the shape of the duct cross-section and the Reynolds number (Re) has not been achieved owing to the diversity of the inertial-focusing phenomena. In this study, we aimed to elucidate the variation of the inertial focusing depending on Re in rectangular duct flows. We performed a numerical simulation of the lift force exerted on a spherical particle flowing in a rectangular duct and determined the lift-force map within the duct cross-section over a wide range of Re. We estimated the particle trajectories based on the lift map and Stokes drag, and identified the particle-focusing points appeared in the cross-section. For an aspect ratio of the duct cross-section of 2, we found that the blockage ratio changes transition structure of particle focusing. For blockage ratios smaller than 0.3, particles focus near the centres of the long sides of the cross-section at low Re and near the centres of both the long and short sides at relatively higher Re. This transition is expressed as a subcritical pitchfork bifurcation. For blockage ratio larger than 0.3, another focusing pattern appears between these two focusing regimes, where particles are focused on the centres of the long sides and at intermediate positions near the corners. Thus, there are three regimes; the transition between adjacent regimes at lower Re is found to be expressed as a saddle-node bifurcation and the other transition as a supercritical pitchfork bifurcation.