Numerical simulation of water and heat transport in the cathode channel of a PEM fuel cell

Cathode channel of a PEM fuel cell is the critical domain for the transport of water and heat. In this study, a mathematical model of water and heat transport in the cathode channel is established by considering two-phase flow of water and air as well as the phase change between water and vapor. The transport process of the species of air is governed by the convection-diffusion equation. The VOSET (coupled volume-of-fluid and level set method) method is used to track the interface between air and water, and the phase equilibrium method of water and vapor is employed to calculate the mass transfer rate on the two-phase interface. The present model is validated against the results in the literature, then applied to investigate the characteristics of two-phase flow and heat transfer in the cathode channel. The results indicate that in the inlet section, water droplets experience three evolution stages: the growing stage, the coalescence stage and the generation stage of dispersed water drops. However, in the middle and outlet sections of the channel, there are only two stages: the growth of water droplets, and the formation of a water film. The mass transfer rate of phase change in the inlet section of the channel varies over time, exhibiting an initial increase, a decrease followed, and a stabilization finally, with the maximum and stable values of 1.78 × 10^?4 kg/s and 1.52 × 10^?4 kg/s for Part 1, respectively. In the middle and outlet sections, the mass transfer rate increase firstly and then keeps stable gradually. Furthermore, regarding the distribution of the temperature and vapor mass fraction in the channel, near the upper surface of the channel, the temperature and vapor mass fraction first change slightly (x < 0.03 m) and then rapidly decrease with fluctuations (x > 0.03 m). In the middle of the channel, the temperature and vapor mass fraction slowly decrease with fluctuation.     

Cold and wet: Diatoms dominate the phytoplankton community during a year of anomalous weather in a Great Lakes estuary

As sentinels of climate change and other anthropogenic forces, freshwater lakes are experiencing ecosystem disruptions at every level of the food web, beginning with the phytoplankton, a highly responsive group of organisms. Most studies regarding the effects of climate change on phytoplankton focus on a potential scenario in which temperatures continuously increase and droughts intersperse heavy precipitation events. Like much of the conterminous United States in 2019, the Muskegon River watershed (Michigan, USA) experienced record-breaking rainfall accompanied by unusually cool temperatures, affording an opportunity to explore how an alternate potential climate scenario may affect phytoplankton. We conducted biweekly sampling of environmental variables and phytoplankton in Muskegon Lake, a Great Lakes Area of Concern that connects to Lake Michigan. We compared environmental variables in 2019 to the previous eight years using long-term data from the Muskegon Lake Observatory buoy, and annual monitoring excursions provided historical phytoplankton data. Under cold and wet conditions, diatoms were the single dominant division throughout the entire growth season – an unprecedented scenario in Muskegon Lake. In 10 of the 13 biweekly sampling days in 2019, diatoms comprised over 75% of the phytoplankton community in the lake by count, indicating that the spring diatom bloom persisted through the fall. Additionally, phytoplankton seasonal succession and abundance patterns typically seen in this lake were absent. In a world experiencing reduced predictability, increased variability, and regional climate anomalies, studying periods of extreme weather events may offer insight into how natural systems will be affected and respond under future climate scenarios.     

Influence of temperature on acoustic emission source location accuracy in underground structure

An acoustic emission (AE) experiment was carried out to explore the AE location accuracy influenced by temperature. A hollow hemispherical specimen was used to simulate common underground structures. In the process of heating with the flame, the pulse signal of constant frequency was stimulated as an AE source. Then AE signals received by each sensor were collected and used for comparing localization accuracy at different temperatures. Results show that location errors of AE keep the same phenomenon in the early and middle heating stages. In the later stage of heating, location errors of AE increase sharply due to the appearance of cracks. This provides some beneficial suggestions on decreasing location errors of structural cracks caused by temperature and improves the ability of underground structure disaster prevention and control.     

Weighted Superposition Attraction (WSA): A swarm intelligence algorithm for optimization problems – Part 2: Constrained optimization

This paper is the second one of the two papers entitled “Weighted Superposition Attraction (WSA) Algorithm”, which is about the performance evaluation of the WSA algorithm in solving the constrained global optimization problems. For this purpose, the well-known mechanical design optimization problems, design of a tension/compression coil spring, design of a pressure vessel, design of a welded beam and design of a speed reducer, are selected as test problems. Since all these problems were formulated as constrained global optimization problems, WSA algorithm requires a constraint handling method for tackling them. For this purpose we have selected 6 formerly developed constraint handling methods for adapting into WSA algorithm and analyze the effect of the used constraint handling method on the performance of the WSA algorithm. In other words, we have the aim of producing concluding remarks over the performance and robustness of the WSA algorithm through a set of computational study in solving the constrained global optimization problems. Computational study indicates the robustness and the effectiveness of the WSA in terms of obtained results, reached level of convergence and the capability of coping with the problems of premature convergence, trapping in a local optima and stagnation.     

A Global Maximum Error Controller Based Nonlinearity Aware TPWL Method for Reducing Nonlinear Systems

Abstract

Model order reduction is a common practice to reduce large order systems so that their simulation and control become easy. Nonlinearity aware trajectory piecewise linear is a variation of trajectory piecewise linearization technique of order reduction that is used to reduce nonlinear systems. With this scheme, the reduced approximation of the system is generated by weighted sum of the linearized and reduced sub-models obtained at certain linearization points on the system trajectory. This scheme uses dynamically inspired weight assignment that makes the approximation nonlinearity aware. Just as weight assignment, the process of linearization points selection is also important for generating faithful approximations. This article uses a global maximum error controller based linearization points selection scheme according to which a state is chosen as a linearization point if the error between a current reduced model and the full order nonlinear system reaches a maximum value. A combination that not only selects linearization points based on an error controller but also assigns dynamic inspired weights is shown in this article. The proposed scheme generates approximations with higher accuracies. This is demonstrated by applying the proposed method to some benchmark nonlinear circuits including RC ladder network and inverter chain circuit and comparing the results with the conventional schemes.     

毒品滥用流行病模型的稳定性分析

现有毒品滥用流行病模型假设吸毒者康复后对毒品拥有永久"免疫"力，而忽视了其再次成为毒品易感者的可能性。针对这一问题，通过考虑社区治疗和隔离治疗两种措施，分析了毒品滥用人群的演化过程，提出了基于暂时"免疫"力的毒品滥用流行病模型，并计算了模型的基本再生数，讨论了模型平衡点的存在性和稳定性。当基本再生数小于1时，模型存在一个局部渐进稳定的无毒平衡点；当基本再生数大于1时，模型存在唯一的地方病平衡点，并利用几何方法证明了地方病平衡点的全局稳定性；当基本再生数等于1时，如果满足一定条件，模型出现后向分支现象。数值模拟验证了上述所有结果。研究结果表明提高隔离治疗率、改善社区治疗效果和降低接触传染率可以有效抑制毒品滥用的流行。     

Impact assessment of reservoir desiltation measures for downstream riverbed migration in climate change: A case study in northern Taiwan

Typhoon Aere in 2004 induced severe sedimentation and loss of storage capacity of the Shihmen Reservoir in northern Taiwan. The resulting dramatic increase in the turbidity of the water seriously affected the water supply. To effectively maintain the stability of the water supply and maintain the reservoir’s storage capacity, the government of Taiwan began to plan and construct a series of improvement measures, such as a sediment flushing tunnel, the JhongJhuang Bank-Side Reservoir, and the Amuping Desilting Tunnel. However, previous studies only focused on the impact of the sediment flushing tunnel and the Amuping Desilting Tunnel on the downstream riverbed, and did not consider the possibility of increasing sediment discharge after the completion of the JhongJhuang Bank-Side Reservoir. In addition, climate change will cause the intensity of extreme rainfall to increase enormously in the future. That rainfall and extra sediment flushing will challenge the existing flood prevention facilities. Therefore, this study considered that the JhongJhuang Bank-Side Reservoir will increase sediment discharge of the Shihmen Reservoir, and used dynamical downscaling extreme typhoon data of climate change under the RCP 8.5 scenario to explore the flood prevention and riverbed migration of the main channels of the Dahan and Tamsui Rivers in the future. We used the rainfall–runoff model of Hydrologic Modeling System to simulate rainfall and runoff, and used the hydraulic and sediment transport model of CCHE1D to holistically simulate flood events and consequent river scouring and deposition behaviors. Our results showed that the projected peak discharge during the late 21st century (2075 to 2099) will be at least 50% higher than that during the baseline (1979 to 2003) period. In terms of flood prevention, the potential of overbank flooding will increase in the downstream area, and the trend of long-term change in the riverbed will be dominated by degradation (-0.489 ± 0.743 m) in the future. The improvement measures will have a limited impact on riverbed migration (0.011 ± 0.094 m) in the Dahan and Tamsui Rivers. After the operation of the JhongJhuang Bank-Side Reservoir, the Shihmen Reservoir is expected to increase the sediment discharge ratio by 70% during floods, and it will not cause excessive water turbidity that may affect downstream water supply.     

Classic and inverse compositional Gauss-Newton in global DIC

Today, effective implementations of digital image correlation (DIC) are based on iterative algorithms with constant linear operators. A relevant idea of the classic finite element (or, more generally, global) DIC solver consists in replacing the gradient of the deformed state image with that of the reference image, so as to obtain a constant operator. Different arguments (small strains, small deformations, equality of the two gradients close to the solution, etc) have been given in the literature to justify this approximation, but none of them are fully accurate. Indeed, the convergence of the optimization algorithm has to be investigated from its ability to produce descent directions. Through such a study, this paper attempts to explain why this approximation works and what is its domain of validity. Then, an inverse compositional Gauss-Newton implementation of finite element DIC is proposed as a cost-effective and mathematically sound alternative to this approximation.     

On global Hessenberg based methods for solving Sylvester matrix equations

In the first part of this paper, we investigate the use of Hessenberg-based methods for solving the Sylvester matrix equation $AX+XB=C$. To achieve this goal, the Sylvester form of the global generalized Hessenberg process is presented. Using this process, different methods based on a Petrov–Galerkin or on a minimal norm condition are derived. In the second part, we focus on the SGl-CMRH method which is based on the Sylvester form of the Hessenberg process with pivoting strategy combined with a minimal norm condition. In order to accelerate the SGl-CMRH method, a preconditioned framework of this method is also considered. It includes both fixed and flexible variants of the SGl-CMRH method. Moreover, the connection between the flexible preconditioned SGl-CMRH method and the fixed one is studied and some upper bounds for the residual norm are obtained. In particular, application of the obtained theoretical results is investigated for the special case of solving linear systems of equations with several right-hand sides. Finally, some numerical experiments are given in order to evaluate the effectiveness of the proposed methods.