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.     

Effect of wind condition on unintended hydrogen release in a hydrogen refueling station

Ambient condition, especially the wind condition, is an important factor to determine the behavior of hydrogen diffusion during hydrogen release. However, only few studies aim at the quantitative study of the hydrogen diffusion in a wind-exist condition. And very little researches aiming at the variable wind condition have been done. In this paper, the hydrogen diffusion in different wind condition which including the constant wind velocity and the variable wind velocity is investigated numerically. When considering the variable wind velocity, the UDF (user defined function) is compiled. Characteristics of the FGC (flammable gas cloud) and the HMF (hydrogen mass fraction) are analyzed in different wind condition and comparisons are made with the no-wind condition. Results indicate that the constant wind velocity and the variable wind velocity have totally different effect for the determination of hydrogen diffusion. Comparisons between the constant wind velocity and the variable wind velocity indicate that the variable wind velocity may cause a more dangerous situation since there has a larger FGC volume. More importantly, the wind condition has a non-negligible effect when considering the HMF along the radial direction. As the wind velocity increases, the distribution of the HMF along the radial direction is not Gaussian anymore when the distance between the release hole and the observation line exceeds to a critical value. This work can be a supplement of the research on the hydrogen release and diffusion and a valuable reference for the researchers.     

Effects of upholstery materials on the burning behavior of real‐scale residential upholstered furniture mock‐ups

Fire spread and growth on real‐scale four cushion mock‐ups of residential upholstered furniture (RUF) were investigated with the goal of identifying whether changes in five classes of materials (barrier, flexible polyurethane foam, polyester fiber wrap, upholstery fabric, and sewing thread), referred to as factors, resulted in statistically significant changes in burning behavior. A fractional factorial experimental design plus practical considerations yielded a test matrix with 20 material combinations. Experiments were repeated a minimum of two times. Measurements included fire spread rates derived from video recordings and heat release rates (HRRs). A total of 13 experimental parameters (3 based on the videos and 10 on the HRR results), referred to as responses, characterized the measurements. Statistical analyses based on Main Effects Plots (main effects) and Block Plots (main effects and factor interactions) were used. The results showed that three of the factors resulted in statistically significant effects on varying numbers of the 13 responses. The Barrier and Fabric factors had the strongest main effects with roughly comparable magnitudes. Foam was statistically significant for fewer of the responses and its overall strength was weaker than for Barrier and Fabric. No statistically significant main effects were identified for Wrap or Thread. Multiple two‐term interactions between factors were identified as being statistically significant. The Barrier*Fabric interaction resulted in the highest number of and strongest statistically significant effects. The existence of two‐term interactions means that it will be necessary to consider their effects in approaches designed to predict the burning behavior of RUF.     

Hall current effects on MHD flow of chemically reacting fluid through a porous medium with heat source

We considered the magnetohydrodynamic (MHD) free convective flow of an incompressible electrically conducting viscous fluid past an infinite vertical permeable porous plate with a uniform transverse magnetic field, heat source and chemical reaction in a rotating frame taking Hall current effects into account. The momentum equations for the fluid flow during absorbent medium are controlled by the Brinkman model. Through the undisturbed state, both the plate and fluid are in a rigid body rotation by the uniform angular velocity perpendicular to an infinite vertical plate. The perpendicular surface is subject to the homogeneous invariable suction at a right angle to it and the heat on the surface varies about a non-zero unvarying average whereas the warmth of complimentary flow is invariable. The systematic solutions of the velocity, temperature, and concentration distributions are acquired systematically by utilizing the perturbation method. The velocity expressions consist of steady-state and fluctuating situations. It is revealed that the steady part of the velocity field has a three-layer characteristic while the oscillatory part of the fluid field exhibits a multi-layer characteristic. The influence of various governing flow parameters on the velocity, temperature, and concentration are analyzed graphically. We also discuss computational results for the skin friction, Nusselt number, and Sherwood number in the tabular forms.     

Estimates of daily oxygen consumption,carbon dioxide and methane emissions,and heat production for beef and dairy cattle using spot gas sampling

A simulation study was conducted to examine accuracy of estimating daily O₂ consumption, CO₂ and CH₄ emissions, and heat production (HP) using a spot sampling technique and to determine optimal spot sampling frequency (FQ). Data were obtained from 3 experiments where daily O₂ consumption, emissions of CO₂ and CH₄, and HP were measured using indirect calorimetry (respiration chamber or headbox system). Experiment 1 used 8 beef heifers (ad libitum feeding; gaseous exchanges measured every 30 min over 3 d in respiration chambers); Experiment 2 used 56 lactating Holstein-Friesian cows (restricted feeding; gaseous exchanges measured every 12 min over 3 d in respiration chambers); Experiment 3 used 12 lactating Jersey cows (ad libitum feeding; gaseous exchanges measured every hour for 1 d using headbox style chambers). Within experiment, averages of all measurements (FQALL) and averages of measurements selected at time points with 12, 8, 6, or 4 spot sampling FQ (i.e., sampling every 2, 3, 4, and 6 h in a 24-h cycle, respectively; FQ12, FQ8, FQ6, and FQ4, respectively) were compared. Within study a mixed model was used to compare gaseous exchanges and HP among FQALL, FQ12, FQ8, FQ6, and FQ4, and an interaction of dietary treatment by FQ was examined. A regression model was used to evaluate accuracy of spot sampling within study [i.e., FQALL (observed) vs. FQ12, FQ8, FQ6, or FQ4 (estimated)]. No interaction of diet by FQ was observed for any variables except for CH₄ production in experiment 1. No FQ effect was observed for gaseous exchanges and HP except in experiment 2 where CO₂ production was less (5,411 vs. 5,563 L/d) for FQ4 compared with FQALL, FQ12, and FQ8. A regression analysis between FQALL and each FQ within study showed that slopes and intercepts became farther from 1 and 0, respectively, for almost all variables as FQ decreased. Most variables for FQ12 and FQ8 had root mean square prediction error (RMSPE) less than 10% of the mean and concordance correlation coefficient (CCC) greater than 0.80, and RMSPE increased and CCC decreased as FQ decreased. When a regression analysis was conducted with combined data from the 3 experiments (mixed model with study as a random effect), results agreed with those from the analysis for the individual studies. Prediction errors increased and CCC decreased as FQ decreased. Generally, all the estimates from FQ12, FQ8, FQ6, and FQ4 had RMSPE less than 10% of the means and CCC greater than 0.90 except for FQ6 and FQ4 for O₂ consumption and CH₄ production. In conclusion, the spot sampling simulation with 3 indirect calorimetry experiments indicated that FQ of at least 8 samples (every 3 h in a 24-h cycle) was required to estimate daily O₂ consumption, CO₂ and CH₄ production, and HP and to detect changes in those in response to dietary treatments. This sampling FQ may be considered when using techniques that measure spot gas exchanges such as the GreenFeed and face mask systems.     

2.3–21 GHz broadband and high linearity distributed low noise amplifier

This paper focuses on the design of a 2.3–21 GHz Distributed Low Noise Amplifier (LNA) with low noise figure (NF), high gain (S₂₁), and high linearity (IIP3) for broadband applications. This distributed amplifier (DA) includes S/C/X/Ku/K-band, which makes it very suitable for heterodyne receivers. The proposed DA uses a 0.18 μm GaAs pHEMT process (OMMIC ED02AH) in cascade architecture with lines adaptation and equalization of phase velocity techniques, to absorb their parasitic capacitances into the gate and drain transmission lines in order to achieve wide bandwidth and to enhance gain and linearity. The proposed broadband DA achieved an excellent gain in the flatness of 13.5 ± 0.2 dB, a low noise figure of 3.44 ± 1.12 dB, and a small group delay variation of ±19.721 ps over the range of 2.3–21 GHz. The input and output reflection coefficients S₁₁ and S₂₂ are less than −10 dB. The input compression point (P_1dB) and input third-order intercept point (IIP3) are −1.5 dBm and 11.5 dBm, respectively at 13 GHz. The dissipated power is 282 mW and the core layout size is 2.2 × 0.8 mm².     

Numerical study on laminar burning velocity of ammonia flame with methanol addition

The combustion characteristics of ammonia/methanol mixtures were investigated numerically in this study. Methanol has a dramatic promotive effect on the laminar burning velocity (LBV) of ammonia. Three mechanisms from literature and another four self-developed mechanisms constructed in this study were evaluated using the measured laminar burning velocities of ammonia/methanol mixtures from Wang et al. (Combust.Flame. 2021). Generally, none of the selected mechanisms can precisely predict the measured laminar burning velocities at all conditions. Aiming to develop a simplified and reliable mechanism for ammonia/methanol mixtures, the constructed mechanism utilized NUI Galway mechanism (Combust.Flame. 2016) as methanol sub-mechanism and the Otomo mechanism (Int. J. Hydrogen. Energy. 2018) as ammonia sub-mechanism was optimized and reduced. The reduced mechanism entitled ‘DNO-NH3’, can accurately reproduce the measured laminar burning velocities of ammonia/methanol mixtures under all conditions. A reaction path analysis of the ammonia/methanol mixtures based on the DNO-NH3 mechanism shows that methanol is not directly involved in ammonia oxidation, instead, the produced methyl radicals from methanol oxidization contribute to the dehydrogenation of ammonia. Besides, NO_x emission analysis demonstrates that 60% methanol addition results in the highest NO_x emissions. The most important reactions dominating the NO_x consumption and production are identified in this study.     

Effect of heat transfer through the release pipe on simulations of cryogenic hydrogen jet fires and hazard distances

Jet flames originated by cryo-compressed ignited hydrogen releases can cause life-threatening conditions in their surroundings. Validated models are needed to accurately predict thermal hazards from a jet fire. Numerical simulations of cryogenic hydrogen flow in the release pipe are performed to assess the effect of heat transfer through the pipe walls on jet parameters. Notional nozzle exit diameter is calculated based on the simulated real nozzle parameters and used in CFD simulations as a boundary condition to model jet fires. The CFD model was previously validated against experiments with vertical cryogenic hydrogen jet fires with release pressures up to 0.5 MPa (abs), release diameter 1.25 mm and temperatures as low as 50 K. This study validates the CFD model in a wider domain of experimental release conditions - horizontal cryogenic jets at exhaust pipe temperature 80 K, pressure up to 2 MPa ab and release diameters up to 4 mm. Simulation results are compared against such experimentally measured parameters as hydrogen mass flow rate, flame length and radiative heat flux at different locations from the jet fire. The CFD model reproduces experiments with reasonable for engineering applications accuracy. Jet fire hazard distances established using three different criteria - temperature, thermal radiation and thermal dose - are compared and discussed based on CFD simulation results.     

Ultrasonic tomographic velocimeter for visualization of axial flow fields in pipes

An ultrasonic tomographic velocimeter to provide quantitative images of axial flow fields in pipes is developed and presented in this work. To detect the flow in various directions and positions, a novel transducer configuration strategy is proposed. All-in-one transducers are mounted in two sectional planes of the pipe. In each plane, N transducers are equally spaced along the circumference. Overlapped propagation paths are introduced by the configuration strategy, and the influence of the vortex flow can be eliminated theoretically by averaging the line velocities of the overlapped paths. To achieve a fast detection speed, the projection data is collected via an electrical scan in a fan-beam mode. After rearrangement and interpolation of the projection data, the parallel beam filtered back projection (FBP) algorithm is implemented to reconstruct the axial flow field. Numerical simulations with the theoretical velocity profiles were performed. The compensation method for the vortex flow is proved to be effective and necessary, and the number of transducers required for reconstruction of common flow profiles was estimated. Accordingly, an ultrasonic tomographic velocimeter consisting of 2×12 transducers was fabricated. Experiments were conducted in the straight pipe and downstream of a single bend pipe and compared with the computational fluid dynamics (CFD) simulation results. As demonstrated, the ultrasonic tomographic velocimeter was capable of visualizing both symmetric and asymmetric axial flow fields with high reliability.     

Inverse identification of constitutive parameters from heat source fields: A local approach applied to hyperelasticity

In this paper, a new inverse identification method of constitutive parameters is developed from full kinematic and thermal field measurements. It consists in reconstructing the heat source field from two different approaches by using the heat diffusion equation. The first one requires the temperature field measurement and the value of the thermophysical parameters. The second one is based on the kinematic field measurement and the choice of a thermo-hyperelastic model that contains the parameters to be identified. The identification is carried out at the local scale, ie, at any point of the heat source field, without using the boundary conditions. In the present work, the method is applied to the challenging case of hyperelasticity from a heterogeneous test. Due to large deformations undergone by the rubber specimen tested, a motion compensation technique is developed to plot the kinematic and the thermal fields at the same points before reconstructing the heterogeneous heat source field. In the present case, the constitutive parameter of the Neo-Hookean model has been identified, and its distribution has been characterized with respect to the strain state at the surface of a cross-shaped specimen.