scholarly journals Fractal-Like Flow-Fields with Minimum Entropy Production for Polymer Electrolyte Membrane Fuel Cells

Entropy ◽  
2020 ◽  
Vol 22 (2) ◽  
pp. 176 ◽  
Author(s):  
Natalya Kizilova ◽  
Marco Sauermoser ◽  
Signe Kjelstrup ◽  
Bruno G. Pollet

The fractal-type flow-fields for fuel cell (FC) applications are promising, due to their ability to deliver uniformly, with a Peclet number Pe~1, the reactant gases to the catalytic layer. We review fractal designs that have been developed and studied in experimental prototypes and with CFD computations on 1D and 3D flow models for planar, circular, cylindrical and conical FCs. It is shown, that the FC efficiency could be increased by design optimization of the fractal system. The total entropy production (TEP) due to viscous flow was the objective function, and a constant total volume (TV) of the channels was used as constraint in the design optimization. Analytical solutions were used for the TEP, for rectangular channels and a simplified 1D circular tube. Case studies were done varying the equivalent hydraulic diameter (Dh), cross-sectional area (DΣ) and hydraulic resistance (DZ). The analytical expressions allowed us to obtain exact solutions to the optimization problem (TEP→min, TV=const). It was shown that the optimal design corresponds to a non-uniform width and length scaling of consecutive channels that classifies the flow field as a quasi-fractal. The depths of the channels were set equal for manufacturing reasons. Recursive formulae for optimal non-uniform width scaling were obtained for 1D circular Dh -, DΣ -, and DZ -based tubes (Cases 1-3). Appropriate scaling of the fractal system providing uniform entropy production along all the channels have also been computed for Dh -, DΣ -, and DZ -based 1D models (Cases 4-6). As a reference case, Murray’s law was used for circular (Case 7) and rectangular (Case 8) channels. It was shown, that Dh-based models always resulted in smaller cross-sectional areas and, thus, overestimated the hydraulic resistance and TEP. The DΣ -based models gave smaller resistances compared to the original rectangular channels and, therefore, underestimated the TEP. The DZ -based models fitted best to the 3D CFD data. All optimal geometries exhibited larger TEP, but smaller TV than those from Murray’s scaling (reference Cases 7,8). Higher TV with Murray’s scaling leads to lower contact area between the flow-field plate with other FC layers and, therefore, to larger electric resistivity or ohmic losses. We conclude that the most appropriate design can be found from multi-criteria optimization, resulting in a Pareto-frontier on the dependencies of TEP vs TV computed for all studied geometries. The proposed approach helps us to determine a restricted number of geometries for more detailed 3D computations and further experimental validations on prototypes.


2020 ◽  
Vol 22 (13) ◽  
pp. 6993-7003 ◽  
Author(s):  
Marco Sauermoser ◽  
Signe Kjelstrup ◽  
Natalya Kizilova ◽  
Bruno G. Pollet ◽  
Eirik G. Flekkøy

We show how we can improve bio-inspired flow field patterns for use in PEMFCs by deviating from Murray's law.



Author(s):  
Yi Han ◽  
Feng Liu ◽  
Xin Ran

In the production process of large-diameter seamless steel pipes, the blank heating quality before roll piercing has an important effect on whether subsequently conforming piping is produced. Obtaining accurate pipe blank heating temperature fields is the basis for establishing and optimizing a seamless pipe heating schedule. In this paper, the thermal process in a regenerative heating furnace was studied using fluent software, and the distribution laws of the flow field in the furnace and of the temperature field around the pipe blanks were obtained and verified experimentally. The heating furnace for pipe blanks was analyzed from multiple perspectives, including overall flow field, flow fields at different cross sections, and overall temperature field. It was found that the changeover process of the regenerative heating furnace caused the temperature in the upper part of the furnace to fluctuate. Under the pipe blanks, the gas flow was relatively thin, and the flow velocity was relatively low, facilitating the formation of a viscous turbulent layer and thereby inhibiting heat exchange around the pipe blanks. The mutual interference between the gas flow from burners and the return gas from the furnace tail flue led to different flow velocity directions at different positions, and such interference was relatively evident in the middle part of the furnace. A temperature “layering” phenomenon occurred between the upper and lower parts of the pipe blanks. The study in this paper has some significant usefulness for in-depth exploration of the characteristics of regenerative heating furnaces for steel pipes.



2004 ◽  
Vol 50 (170) ◽  
pp. 342-352 ◽  
Author(s):  
Perry Bartelt ◽  
Othmar Buser

AbstractAn essential problem in snow science is to predict the changing form of ice grains within a snow layer. Present theories are based on the idea that form changes are driven by mass diffusion induced by temperature gradients within the snow cover. This leads to the well-established theory of isothermal- and temperature-gradient metamorphism. Although diffusion theory treats mass transfer, it does not treat the influence of this mass transfer on the form — the curvature radius of the grains and bonds — directly. Empirical relations, based on observations, are additionally required to predict flat or rounded surfaces. In the following, we postulate that metamorphism, the change of ice surface curvature and size, is a process of thermodynamic optimization in which entropy production is minimized. That is, there exists an optimal surface curvature of the ice grains for a given thermodynamic state at which entropy production is stationary. This state is defined by differences in ice and air temperature and vapor pressure across the interfacial boundary layer. The optimal form corresponds to the state of least wasted work, the state of minimum entropy production. We show that temperature gradients produce a thermal non-equilibrium between the ice and air such that, depending on the temperature, flat surfaces are required to mimimize entropy production. When the temperatures of the ice and air are equal, larger curvature radii are found at low temperatures than at high temperatures. Thus, what is known as isothermal metamorphism corresponds to minimum entropy production at equilibrium temperatures, and so-called temperature-gradient metamorphism corresponds to minimum entropy production at none-quilibrium temperatures. The theory is in good agreement with general observations of crystal form development in dry seasonal alpine snow.



2013 ◽  
Vol 291-294 ◽  
pp. 1981-1984
Author(s):  
Zhang Xia Guo ◽  
Yu Tian Pan ◽  
Yong Cun Wang ◽  
Hai Yan Zhang

Gunpowder was released in an instant when the pill fly out of the shell during the firing, and then formed a complicated flow fields about the muzzle when the gas expanded sharply. Using the 2 d axisymmetric Navier-Stokes equation combined with single equation turbulent model to conduct the numerical simulation of the process of gunpowder gass evacuating out of the shell without muzzle regardless of the pill’s movement. The numerical simulation result was identical with the experimental. Then simulated the evacuating process of gunpowder gass of an artillery with muzzle brake. The result showed complicated wave structure of the flow fields with the muzzle brake and analysed the influence of muzzle brake to the gass flow field distribution.



1979 ◽  
Vol 21 (1) ◽  
pp. 1-6 ◽  
Author(s):  
D. Adler ◽  
Y. Levy

A laser-Doppler technique is successfully applied to measure the flow field inside a closed, backswept impeller, through a rotating window. Results show that, in contrast to the flow in many radial-exit impellers, the flow in the backswept impeller is stable and attached. Further, comparison with an open impeller demonstrates the fundamental difference in the flow fields near the shroud.



2006 ◽  
Vol 4 (3) ◽  
pp. 365-368 ◽  
Author(s):  
Guo-Bin Jung ◽  
Ay Su ◽  
Cheng-Hsin Tu ◽  
Fang-Bor Weng ◽  
Shih-Hung Chan

The flow-field design of direct methanol fuel cells (DMFCs) is an important subject about DMFC performance. Flow fields play an important role in the ability to transport fuel and drive out the products (H2O,CO2). In general, most fuel cells utilize the same structure of flow field for both anode and cathode. The popular flow fields used for DMFCs are parallel and grid designs. Nevertheless, the characteristics of reactants and products are entirely different in anode and cathode of DMFCs. Therefore, the influences of flow fields design on cell performance were investigated based on the same logic with respect to the catalyst used for cathode and anode nonsymmetrically. To get a better and more stable performance of DMFCs, three flow fields (parallel, grid, and serpentine) utilized with different combinations were studied in this research. As a consequence, by using parallel flow field in the anode side and serpentine flow-field in the cathode, the highest power output was obtained.



2009 ◽  
Vol 132 (1) ◽  
Author(s):  
Alessandro Armellini ◽  
Filippo Coletti ◽  
Tony Arts ◽  
Christophe Scholtes

The present contribution addresses the aerothermal, experimental, and computational studies of a trapezoidal cross-sectional model simulating a trailing edge cooling cavity with one rib-roughened wall. The flow is fed through tilted slots on one side wall and exits through straight slots on the opposite side wall. The flow field aerodynamics is investigated in Part I of the paper. The reference Reynolds number is defined at the entrance of the test section and set at 67,500 for all the experiments. A qualitative flow model is deduced from surface-streamline flow visualizations. Two-dimensional particle image velocimetry measurements are performed in several planes around midspan of the channel and recombined to visualize and quantify three-dimensional flow features. The crossing-jets issued from the tilted slots are characterized and the jet-rib interaction is analyzed. Attention is drawn to the motion of the flow deflected by the rib-roughened wall and impinging on the opposite smooth wall. The experimental results are compared with the numerical predictions obtained from the finite volume Reynolds-averaged Navier–Stokes solver, CEDRE.



Author(s):  
Ayman Al-Sukhon ◽  
Mostafa SA ElSayed

In this paper, a novel multiscale and multi-stage structural design optimization procedure is developed for the weight minimization of hopper cars. The procedure is tested under various loading conditions according to guidelines established by regulatory bodies, as well as a novel load case that considers fluid-structure interaction by means of explicit finite elements employing Smoothed Particle Hydrodynamics. The first stage in the design procedure involves topology optimization whereby optimal beam locations are determined within the design space of the hopper car wall structure. This is followed by cross-sectional sizing of the frame to concentrate mass in critical regions of the hopper car. In the second stage, hexagonal honeycomb sandwich panels are considered in lower load regions, and are optimized by means of Multiscale Design Optimization (MSDO). The MSDO drew upon the Kreisselmeier–Steinhausser equations to calculate a penalized cost function for the mass and compliance of a hopper car Finite Element Model (FEM) at the mesoscale. For each iteration in the MSDO, the FEM was updated with homogenized sandwich composite properties according to four design variables of interest at the microscale. A cost penalty is summed with the base cost by comparing results of the FEM with the imposed constraints. Efficacy of the novel design methodology is compared according to a baseline design employing conventional materials. By invoking the proposed methodology in a case study, it is demonstrated that a mass savings as high as 16.36% can be yielded for a single hopper car, which translates into a reduction in greenhouse gas emissions of 13.09% per car based on available literature.



2021 ◽  
Author(s):  
Rishabh Prakash Sharma ◽  
Max P. Cooper ◽  
Anthony J.C. Ladd ◽  
Piotr Szymczak

<p>Dissolution of porous rocks by reactive fluids is a highly nonlinear process resulting in a variety of dissolution patterns, the character of which depends on physical conditions such as flow rate and reactivity of the fluid. Long, finger-like dissolution channels, “wormholes”, are the main subject of interest in the literature, however, the underlying dynamics of their growth remains unclear. </p><p>While analyzing the tomography data on wormhole growth.  one open question is to define the exact position of the tip of the wormhole. Near the tip the wormhole gradually thins out and the proper resolution of its features is hindered by the finite spatial resolution of the tomographs. In particular, we often observe in the near-tip region several disconnected regions of porosity growth, which - as we hypothesized - are connected by the dissolution channels at subpixel scale. In this study, we show how these features can be better resolved by using numerically calculated flow fields in the reconstructed pore-space. </p><p>We used 70 micrometers, 16-bit grayscale X-ray computed microtomography (XCMT) time series scans of limestone cores, 14mm in diameter and 25mm in length. Scans were performed during the entire dissolution experiment with an interval of 8 minutes. These scans were further processed using a 3-phase segmentation proposed by Luquot et al.[1], in which grayscale voxels are converted to macro-porosity, micro-porosity and grain phases from their grayscale values. The macro-porous phase is assigned a porosity of 1, while the grain phase is assigned 0. Micro-porous regions are assigned an intermediate value determined by linear interpolation between pore and grain threshold using grayscale values. An OpenFOAM based, Darcy-Brinkman solver, porousFoam, is then used to calculate the flow field in this extracted porosity field. </p><p>Porosity contours reconstructed from the tomographs show some disconnected porosity growth near the tip region which later become part of the wormhole in subsequent scans. We have used a novel approach by including the micro-porosity phase in pore-space to calculate the flow-fields in the near-tip region. The calculated flow fields clearly show an extended region of focused flow in front of the wormhole tip, which is a manifestation of the presence of a wormhole at the subpixel scale. These results show that micro-porosity plays an important role in dissolution and 3-phase segmentation combined with the flow field calculations is able to capture the sub-resolved dissolution channels. </p><p> </p><p> [1] Luquot, L., Rodriguez, O., and Gouze, P.: Experimental characterization of porosity structure and transport property changes in limestone undergoing different dissolution regimes, Transport Porous Med., 101, 507–532, 2014</p>



Author(s):  
I. Janajreh ◽  
C. Ghenai

Large scale wind turbines and wind farms continue to evolve mounting 94.1GW of the electrical grid capacity in 2007 and expected to reach 160.0GW in 2010 according to World Wind Energy Association. They commence to play a vital role in the quest for renewable and sustainable energy. They are impressive structures of human responsiveness to, and awareness of, the depleting fossil fuel resources. Early generation wind turbines (windmills) were used as kinetic energy transformers and today generate 1/5 of the Denmark’s electricity and planned to double the current German grid capacity by reaching 12.5% by year 2010. Wind energy is plentiful (72 TW is estimated to be commercially viable) and clean while their intensive capital costs and maintenance fees still bar their widespread deployment in the developing world. Additionally, there are technological challenges in the rotor operating characteristics, fatigue load, and noise in meeting reliability and safety standards. Newer inventions, e.g., downstream wind turbines and flapping rotor blades, are sought to absorb a larger portion of the cost attributable to unrestrained lower cost yaw mechanisms, reduction in the moving parts, and noise reduction thereby reducing maintenance. In this work, numerical analysis of the downstream wind turbine blade is conducted. In particular, the interaction between the tower and the rotor passage is investigated. Circular cross sectional tower and aerofoil shapes are considered in a staggered configuration and under cross-stream motion. The resulting blade static pressure and aerodynamic forces are investigated at different incident wind angles and wind speeds. Comparison of the flow field results against the conventional upstream wind turbine is also conducted. The wind flow is considered to be transient, incompressible, viscous Navier-Stokes and turbulent. The k-ε model is utilized as the turbulence closure. The passage of the rotor blade is governed by ALE and is represented numerically as a sliding mesh against the upstream fixed tower domain. Both the blade and tower cross sections are padded with a boundary layer mesh to accurately capture the viscous forces while several levels of refinement were implemented throughout the domain to assess and avoid the mesh dependence.



Sign in / Sign up

Export Citation Format

Share Document