Investigation of fabrication of gas diffusion substrate for proton exchange membrane fuel cells

The gas diffusion substrate (GDS) is essential in the proton exchange membrane fuel cells. Its fabrication techniques affect the performance significantly and are worthy of investigation. In this study, a manufacturing process of the GDS is proposed to understand the formation process of GDS and promote its structure and performance more pertinently. Different states during the preparation process, raw carbon paper, pre-curing, curing, carbonation, and graphitization, are characterized and measured. Experimental and numerical methods are employed to determine the relationships between microstructure, transport, and mechanical performance variation with the fabricating processes. The results show that its porosity, average pore size, and effective diffusivity decrease first and increase after curing. These parameters after graphitization are lower than that of the carbon paper (CP). The electrical resistivity increases dramatically while pre-curing and decreases gradually after curing, carbonation, and graphitization, and it is much reduced after graphitization. Moreover, mechanical measurement results show that both the picks of tensile strength and flexural modulus occur after curing. Its tensile strength shows little change after graphitization compared to the initial paper's. In contrast, the flexural modulus is improved significantly.     

One-step to prepare high-performance gas diffusion layer (GDL) with three different functional layers for proton exchange membrane fuel cells (PEMFCs)

Gas diffusion layer (GDL) is one of the most important components of fuel cells. In order to improve the fuel cell performance, GDL has developed from single layer to dual layers, and then to multiple layers. However, dual or multi layers in GDL are usually prepared by layer-by-layer methods, which cost too much time, energy, and resources. In this work, we successfully developed a facile one-step method to prepare a GDL with three functional layers by utilizing the different sedimentation rates and filtration rates of short carbon fiber (CF) and carbon nanotube (CNT). The treatment temperature for this GDL is much lower than that of traditional method. The thickness of the GDL can be effectively controlled from as thin as 50 μm to more than 200 μm by simply adjusting the content of CF. The GDL with high flexibility is suitable to develop high performance flexible electronics. The fuel cell with the GDL has the maximum power density 1021 mW cm^?2, which shows 19% improvement comparing to the conventional one. Therefore, this work breaks the traditional concept that GDL for fuel cells only can be prepared by very complex and high-cost procedure.     

Study on ice-melting performance of gradient gas diffusion layer in proton exchange membrane fuel cell

In this study, a three-dimensional model was established using the lattice Boltzmann method (LBM) to study the internal ice melting process of the gas diffusion layer (GDL) of the proton exchange membrane fuel cell (PEMFC). The single-point second-order curved boundary condition was adopted. The effects of GDL carbon fiber number, growth slope of the number of carbon fibers and carbon fiber diameter on ice melting were studied. The results were revealed that the temperature in the middle and lower part of the gradient distribution GDL is significantly higher than that of the no-gradient GDL. With the increase of the growth slope of the number of carbon fiber, the temperature and melting rate gradually increase, and the position of the solid-liquid interface gradually decreases. The decrease in the number of carbon fibers has a similar effect as the increase in the growth slope of the number of carbon fibers. In addition, as the diameter of the carbon fiber increases, the position of the solid-liquid interface gradually decreases first and then increases.     

Numerical simulation of a PEM fuel cell: Effect of tortuosity parameters on the construction of polarization curves

Fuel cells (PEMFCs) are considered a clean alternative for the production of electricity. To improve their efficiency, it is necessary to understand the transport phenomena that determine their performance. This work analyzes the effective diffusion coefficients in the gas diffusion and catalyst layers (GDL and CL, respectively) based on theories such as effective medium approximation and percolation theory of a 3D multiphase nonisothermal model of a single channel PEMFC. We calculate polarization curves with different effective diffusion models and tortuosity factors m, n using OpenFOAM. The best model that approaches to experimental data is the Tomadakis-Sotirchos model with n = 4.0. The exponent m that represents the tortuosity due to the geometry has a high impact on the polarization curve construction compared to the exponent n. High values of diffusibility do not mean high values of current densities showing that there are other phenomena involved, such as water flooding.     

Microstructure,mechanical, tribological,and oxidizing properties of AlCrSiN/AlCrVN/AlCrNbN multilayer coatings with different modulated thicknesses

Multilayer structure design is one of the most promising methods for improving the comprehensive performance of AlCrN-based hard coatings applied to cutting tools. In this study, four types of AlCrSiN/AlCrVN/AlCrNbN multilayer coatings, with different modulated thicknesses, were deposited to investigate their microstructure, mechanical, tribological, and oxidizing properties. All multilayer coatings exhibited grain growth along the crystallographic plane of (200) with a NaCl-type face-centered cubic (FCC) structure. The results show that, as the modulation thickness decreases from ～35 nm to ～10 nm, (1) the grain refinement effect is increasingly evident; (2) all multilayer coatings show a hardness of >30 GPa and an elastic modulus of >300 GPa. Both the ability to resist elastic strain to failure and the plastic deformation of multilayer coatings increase. In addition, their resistance to cracking reduces; (3) the wear rates of these multilayer coatings reduce successively from 1.78 × 10^?16 m³ N^?1 m^?1 to 7.7 × 10^?17 m³ N^?1 m^?1. This is attributed to an increase in self-lubricating VO_x and a decrease in adhesives from the counterparts; (4) the best high-temperature oxidation resistance was obtained for the multilayer coating with a modulated thickness of ～15 nm.     

α钛富氧层增厚动力学模型

为了实现在工业化生产中对α钛富氧层厚度预测和控制,通过实验研究α钛富氧层在高温空气环境中的形成及增厚过程,讨论热处理温度和时间的影响作用,建立高温(750~850℃)空气环境下关于温度、时间的富氧层增厚动力学模型。结果表明:当恒温热处理温度为750~850℃时,α钛富氧层厚度x与保温时间t0.5呈正比例关系,且升高热处理温度可显著提高富氧层增厚速度。在此温度范围内,氧原子的扩散激活能约为203473 J/mol,计算曲线与实验数据吻合性较好。结合文献中已有的扩散系数方程和实验测得的富氧层厚度数据,推导得到5个富氧层增厚动力学方程,其中3个方程的计算曲线与实验数据吻合性较好,可为实际生产中预估富氧层厚度提供理论支持。     

钨中空位及其团簇的能量学和动力学性质参数

在总结前人钨中空位及其团簇的能量学和动力学行为的研究成果基础上，采用第一性原理方法系统计算了钨中空位及其团簇的结合能和扩散能垒。研究发现，交换关联泛函PW91和PBE较PBEsol、AM05和LDA更适合用于计算钨空位的能量学性质。基于第一性原理计算结果对文献中单空位形成能、双空位作用性质等争议性问题进行了讨论，并对钨经验势进行了评估。研究结果表明，钨中孤立单空位间总是相互排斥，而空位团簇（V_n>3）对单空位具有很强的吸引作用，其结合能随着所含空位个数增多呈现波动性增大的趋势。空位团簇稳定结构可通过最小化Wigner-Seitz表面积来确定，其结合能与V_n与V_n-1之间的Wigner-Seitz面积之差呈正比。     

Simulations of NBI fast ion loss in the presence of toroidal field ripple on EAST

NBI fast ion losses in the presence of the toroidal field ripple on EAST have been investigated by using the orbit code GYCAVA and the NBI code TGCO. The ripple effect was included in the upgraded version of the GYCAVA code. It is found that loss regions of NBI fast ions are mainly on the low field side near the edge in the presence of ripple. For co-current NBIs, the synergy effect of ripple and Coulomb collision on fast ion losses is dominant, and fast trapped ions located on the low field side are easily lost. The ripple well loss and the ripple stochastic loss of fast ions have been identified from the heat loads of co-current NBI fast ions. The ripple stochastic loss and the collisioninduced loss are much larger than the ripple well loss. Heat loads of lost fast ions are mainly localized on the right side of the radio frequency wave antennas from the inside view toward the first wall. For counter-current NBIs, the first orbit loss due to the magnetic drift is the dominant loss channel. In addition, fast ion loss fraction with ripple and collision for each NBI linearly increases with the effective charge number, which is related to the pitch angle scattering effect.     

Kinetic analysis of crystallization of zeolite beta synthesized by direct heating

To quantitatively investigate the initial crystallization of zeolite beta synthesized by direct heating, the extent of the reaction was precisely evaluated by X-ray diffraction measurements and Rietveld structural refinement, and a kinetic analysis of crystallization was performed using the Avrami-Erofe'ev equation. The activation energy for crystallization was lower than that for hydrothermal synthesis. Reaction and synthesis time curves revealed that the initial zeolite beta crystallization consisted of three stages. The first was an induction period with nucleation by the generation of building units and the formation of an initial coordinated structure. The second stage was crystal growth by a diffusion-controlled reaction, and the third stage involved slowing down of crystallization by the limitation of dehydrocondensation. These stages could be analyzed by calculation of the rate constant and Avrami exponent for each stage.     

Sizing of multi-stage Op Amps by combining design equations with the gm/ID method

Analog integrated circuit design has as integral parts both analytical reasoning and numerical validation in the process from topology construction to sizing. Given a circuit topology, different circuit sizing results can be obtained from different processes of sizing inference. Sizing methods by simulation-based numerical searching have been a continuously studied subject. However, almost all approaches in this category require an overwhelming number of circuit simulations to arrive at an optimized sizing result. On the other hand, many published manual sizing methods by using the conventional device equations also require repeated SPICE simulations to correct the equation-based sizing results. This paper proposes a systematic gm/ID-based initial sizing method specifically customized for designing multiple-stage operational amplifiers (Op Amps). A main feature of the proposal is to use circuit-level design equations as constraints on the gm/ID table lookup method to substantially reduce the uncertainty in the sizing calculations. As a result, a significant amount of SPICE based correction work can be reduced to complete an initial sizing. The proposed sizing procedure includes a few regular sizing rules customized to the configuration of multi-stage Op Amps. We validate the proposed sizing method by application to several multi-stage Op Amp examples with a capacitive load or Miller compensation. Simulations have justified that the produced initial sizing results can achieve most of the prespecified design targets.