Structural Analysis of Electrochemical Hydrogen Compressor End-Plates for High-Pressure Applications

: Transportation and portable applications already use hydrogen as fuel, but it is essential to use highly-efﬁcient hydrogen storage methods to increase its usage in the future. The compressed form is the most utilized for transportation applications, but mechanical compressors have low efﬁciency when compressing low quantities of gas to high pressure. The most suitable device for hydrogen compression is the Electrochemical Hydrogen Compressor (EHC). It has the same structure as a Proton Exchange Membrane Fuel Cell (PEM-FC), but it works at very high-pressure (700 bar). The present work analyses the monopolar plate of an Electrochemical Hydrogen Compressor prone to hydrogen embrittlement. Irregular shape variations generate peaks of stress magnitude and triaxiality, further contributing to decreasing metal ductility at the local scale. The calculation of the stress ﬁeld in such components is essential due to the possibility of failure due to the material embrittlement caused by hydrogen. The paper presents a conceptual design of an EHC operating at 700 bar and focuses on the shape and the mechanical stress of the end-plates to have conservative levels of the nominal stress states, which then are taken as the design parameter for providing adequate structural integrity and mechanical reliability to the component. The FEM analysis with Marc software—of MSC Software Corporation—identiﬁed the optimal end-plates conﬁguration in circular plan view and active area. The plate, sized to have a deﬂection no greater than 0.1mm when the EHC works at 700 bar, should have the minimum thickness of 17mm.


Introduction
Today fossil fuels are the primary energy source, especially in the automotive field. Unfortunately, they are rapidly depleting, and their combustion products are environmentally unfriendly [1]. For a long time, improved fuel efficiency has been the goal for hybrid power-train technologies. Only recently, hybrid vehicles were also considered a solution for reducing carbon-emitting, but they still need a certain amount of fossil fuels to work [2]. For Zero-Emission Vehicles (ZEV), the key is to use hydrogen as a fuel in a fuel cell system [3].
Transportation, portable and stationary applications already use hydrogen as fuel [4]. Except for stationary fuel cells that utilize a fixed platform to store the energy, it is important to use extremely efficient hydrogen storage methods from both gravimetric and volumetric points of view [3,5].
Hydrogen has a specific energy of approximately 120 MJ/kg (almost three times more than gasoline) but has the lowest energy density by volume as gaseous hydrogen at room temperature and pressure [1,2]. The compressed form is the most utilized in transportation applications. The typical pressures required to have the same range of autonomy as traditional vehicles are 350 bar and 700 bar [8,9] (hydrogen refueling stations have to be capable of compressing hydrogen up to 1000 bar [10]), with a refueling time of about 5 min [11].
The structural design of the components in contact with gaseous hydrogen must consider the known phenomenon of hydrogen embrittlement. It deteriorates a particular feature of structural metals, such as their ductility. Surface cracks also occur on unstressed components subjected to cathodic charging [12]. Ref. [13] proposed the adoption of "high entropy alloys" (HEA) as a countermeasure for limiting the decrease of ductility as the hydrogen content increases in time, while [14,15] thoroughly investigated the detrimental effect of hydrogen absorption in AISI stainless steels by also quantifying the decrease of fatigue allowed cycles and monotonic true stress-strain curves promoted by different amounts of adsorbed H2.
The compression of a gas in a mechanical compressor is an adiabatic process that needs more than one stage to reduce temperature [10] and has low efficiency when used to compress low quantities of gas to high pressure. In the best case, it requires energy approximately 1.4 times above the ideal isothermal compression [16]. Furthermore, the contamination of the compressed gas is a risk for such application as fuel cells that requires a 99.99% hydrogen purity [17]. For these cases, another type of device, the Electrochemical Hydrogen Compressor (EHC), is more suitable: it works electrochemically, without moving parts, vibrations, and noise, operating at lower hydrogen flux than traditional compressors [18]; furthermore, its process is isothermal in principle [19] and can operate in a single stage, so the EHC is simpler and more efficient than mechanical compressors [20]. While compressing, it also operates the purification of hydrogen from carbon monoxide if it is present after the reforming process [21].
As for FC design, also the EHC needs FEM (Finite Element Method) computer modeling to optimize its components (electrodes, catalysts, membrane, monopolar and bipolar plates, gaskets, and bolts) and to improve the understanding of the mechanisms that govern the cell [22]. Furthermore, for the EHC is necessary to scale up the structural components to higher pressures. In the literature, a few researchers have been dedicated to this purpose. Baker et al. [23] compared different geometrical configurations of gas distribution channels using computational fluid dynamics to determine the most effective. Other studies investigated improving the EHC structure to increase efficiency and reduce hydrogen back diffusion [24].
This work aims to analyze the monopolar plate of an Electrochemical Hydrogen Compressor operating at high pressure (700 bar) to guarantee conservative levels of the nominal stress states also when the stress peaks (due to the irregularities in the shape of the mechanical component) combined to the ductility decrease (due to the hydrogen embrittlement) can increase the possibility of local failure. The conceptual design starts by comparing three different geometric configurations of the end-plate to find the optimal design that guarantees structural integrity and mechanical reliability.

Electrochemical Hydrogen Compressor
The Electrochemical Hydrogen Compressor (EHC) has the same structure as a Proton Exchange Membrane Fuel Cell (PEM-FC). A membrane, positioned between two catalysts containing electrodes, acts as an electrolyte for protons exchange [25]. Hydrogen is oxidized at the anode and reduced at the cathode by a power supply: in this way, hydrogen is transferred from the low-pressure (LP) anodic compartment ( Figure 1, left side) to the high-pressure (HP) cathodic compartment ( Figure 1, right side).

Figure 1.
How an Electrochemical Hydrogen Compressor works: the polymer electrode membrane is positioned between two catalysts containing electrodes and allows the H2 protons to move from the left to the right side. H2 electrons flow in the same direction by a power supply. Then, at the cathode, the protons and electrons form hydrogen again, but at higher pressure.
Electrochemical reactions taking place in the process are [24]: The hydrogen can be pressurized up to several hundred bars with a single cell [26]. Furthermore, the EHC operates hydrogen separation, compression, and purification in only one compact device.
The main components of an EHC stack are: • MEA (Membrane Electrode Assembly), constituted by two (anodic and cathodic) electrodes, two (anodic and cathodic) catalyst layers, and a proton conductive membrane; • bipolar plates, made of stainless steel. They contain the gas flow field; • end-plates, made of stainless steel or carbon fiber composite; • bolts. On the end-plates, they ensure the gas tightening of the system; • gaskets.

Numerical Analysis
This study aims to optimize the design of an EHC single cell. In this case, the endplates are monopolar plates, which means they contain the flow field only on their internal side. The analysis is organized into two steps: 1.
in the first step, different geometries of the monopolar plate are compared to identify the best configuration; 2.
once the best geometric configuration had been evaluated, the second part of the study concerned the design of the plate.
Simulations has been performed by Marc Software of MSC Software Corporation. MSC Marc finite elements approach is especially suitable for involving nonlinear aspects of the numerical analyses, such as contact, here modeled, or for future developments of the work, such as implicit elastoplastic simulations or damage, and cohesive zone method for modeling damage and fracture under limit conditions. The software also allows for easily customizing the mathematical modeling of the material's response via very flexible Fortran user subroutines.
Geometric details, such as the channels of the flow-field, were not taken into account: they were replaced with a groove whose dimensions correspond to the active area (the area in which chemical reactions occur). The depth of the groove was calculated by summing: • the space for the MEA (overall thickness of 0.35 mm, hence 0.175 mm per plate); • the depth the flow-field channels would have had if represented (0.5 mm); • 1.5 mm corresponding to the thickness of a metallic mesh to put over the low-pressure channels to better support the MEA during operation.
The total depth of the groove is 2.175 mm for the low-pressure side and 0.675 mm for the high-pressure side. The FEM analyses were carried on for the low-pressure plate because they yielded a conservative estimate for the high-pressure plate.
The mesh was created manually along with the geometry using brick elements that have three degrees of freedom per node. Due to the high corrosiveness of the compression process, the chosen material was AISI 316 stainless steel (European designation: X5CrNiMo17-12-2/1.4401). It is an austenitic chromium-nickel stainless steel containing molybdenum; it has excellent ductility with high corrosion resistance and mechanical strength at elevated temperatures. Its elastic properties, listed in Table 1, have references in [27-30]. The working pressure of 700 bar acts on the upper side of the plate, inside the active area.

Part 1: Methodology
The diffusion of hydrogen within mechanical components made of stainless steel is known to induce the embrittlement phenomenon, consisting of the reduction of the plastic strains locally attainable before failure [12][13][14][15]. Consequently, the ultimate tensile stress and the fatigue response of stainless steel are negatively affected by pressure-driven diffusion of hydrogen within the alloy [14], while surface fracture can even be induced by cathodic charging with no outer stress applied [12]. In the application at hand, the contact between the surface of the compressor end-plates and the hydrogen flow is rather limited to the small-sized orifices and channels where the hydrogen is allowed to flow, so no large amounts of hydrogen can be expected to diffuse within the volume of the components. At the same time, such small-sized channels, holes, and irregular shape variations act as notches within the stressed mechanical component at hand, thus generating peaks of stress triaxiality, which further contribute to decreasing the ductility of metals at the local scale [31,32]. According to such considerations, the design of mechanical components to be continuously operated in hydrogen environments must be aimed at conservative levels of the nominal stress states. This is the reason why the stress states have been taken as the design parameter to be evaluated and designed by the finite element method for providing adequate structural integrity to the component.
FC monopolar and bipolar plates at the CNR laboratories utilize the configuration "a" represented in Figure 2: it has a squared plan view and squared active area (the area in which chemical reactions occur). This configuration is usually used in Fuel Cell applications because it maximizes the utilization of the active area. An axisymmetric geometry is more suitable for EHC application due to the high pressure on the active area. For this purpose, the comparison between three different geometric configurations (Figures 2-4) was initially carried on.   The plate sizes were calculated starting from the FC monopolar plates used at the CNR laboratories with the following features: • configuration "a"; • square side of plan view: 100 mm; • square side of active area: 50 mm.
The other configurations have a circular plan view corresponding to the circle that circumscribes the square of the first configuration; a characteristic dimension greater than the one of case "a" (circle diameter against square side) was chosen to make the problem worse (larger size means much deflection). It was calculated as: Finally, the circular active area of the last configuration corresponds to the circle that has the same area as the squared active area of the other two configurations: About the sizing along the z-axis, the thickness of the plate is 30 mm for the three configurations. As stated above, the active area is deeper ( Figure 5) than the plate surface.
All the dimensions are summarized in Table 2.
The comparison between geometries does not consider the bolts' holes that ensure the system's tightening: the plate was constrained by 3D pins placed at 12 points on the plate's lower-not loaded-side. The number of 3D pins was chosen considering the perimeter of the plate.

Part 1: Results
Maximum values of deflection, equivalent strain, and Von Mises stress obtained from FEM simulations about the three configurations are reported in Table 3. The results show that configurations "b" and "c" have a significant decrease in deflection compared to the deflection of configuration "a"; in terms of strains and stresses, excluding the constrained points that distort the actual trend, a drastic increase occurs in configuration "b", instead configuration "c" is the best in terms of stress. Based on the results of this comparison, the plate's design concerned the third configuration.

Part 2: Methodology
The goal of this work section was to size the monopolar plate (configuration "c") to have a deflection no greater than 0.1 mm when the EHC works at its operational pressure (700 bar). This design value should ensure wide physical contact between the MEA and the high-pressure monopolar plate.
As the plate material and the number of constrained points cannot be changed, the only way to ensure a deflection lower than 0.1 mm was to reduce the effective length of the plate as much as possible and to find its optimal (minimum) thickness. For this purpose, bolts' holes and gasket's groove that ensures the gas tightening of the system were taken into account.
Twelve high-strength M8 screws (class 10.9) with nuts and washers (with an external diameter of 16 mm) were chosen. The minimum diameter of the screw was calculated as: where: • p is the pressure of 700 bar; • S is 2507 mm 2 active area; • n is the number of bolts; • γ is a safety factor of 3; • σ y is the yield stress of the bolt material.
The gasket for the gas tightening of the EHC was chosen for axial load with internal pressure ( Figure 6). In this condition, O-rings do not allow the use of Back-up Rings to avoid the extrusion of the gasket in the gap between the components to be connected. Instead, some metal rings, used as a deformable gasket in extreme static sealing situations, can increase their sealing capacity thanks to one or more coatings of materials, such as nickel, silver, gold, copper, and PTFE.
The metal ring was chosen to keep the groove width b as small as possible because it increases the plate's effective length ( Figure 6). The ring has a free height (cross-section diameter of the ring) of 2.38 mm and an external diameter inside the range of 25 ÷ 300 mm. As reported in its technical data-sheet, the groove width is b = 3.57 mm, while the groove depth is h = 1.83 + 0.08 mm. The external diameter of the groove was calculated as: where l MEA is the active area diameter, s = 4 mm is the height of the membrane annulus that it was needed to add to the active area, and b is the groove width.
Once known the external diameter of the groove, it will be possible to calculate the external diameter of the metal ring as: where 0.3 mm is a clearance tolerance correction value, and t c is the maximum coating thickness if used. The optimized sizing of the plate are shown in Figure 7. The thickness of the plate has been set to 17 mm because results will show that the deflection is less than 0.09 mm. Once the geometry was defined, constraints were described more accurately than in the previous case to distribute contact stresses on a surface rather than in a point: both the body's bolts and the washers were represented with infinitely rigid surfaces, the former with cylinders, the latter with circles. Constraints were represented as a touching contact between the deformable surface of the plate and the infinitely rigid surfaces of bolts and washers, ensuring that the surfaces were correctly oriented.
Additional nodes on the bottom surface of the plate were constrained to avoid plate rotation during the simulation; in particular, it was locked: • central node displacements in XY plane; • displacement of a y-axis point along the x-axis; • displacement of an x-axis point along the y-axis.

Part 2: Results
Figures 8-11 shows the Von Mises stress, the deflection, the x-axis displacement, and the equivalent strain distributions.    at the center of the plate (about 250 MPa) due to the high pressure applied orthogonally to the surface; 2.
in the contact areas (the automatic full-scale value, not represented in Figure 8, is 695.7 MPa).
In both cases, stress values are higher than AISI 316 yield stress (σ y = 200 MPa, Table 1). Despite this, the plate still has an elastic behavior.
The first case is a hydrostatic compression that does not plasticize the material: the area with the maximum deflection is unloaded in the central axial view of the plate because it is compressed so much above, so much below. This fact avoids the plastic behavior of the material even under stresses that moderately exceed the yield strength without triggering pressure-dependent brittle failure issues.
In the second case, the stress is due to the Hertzian contact between the plate and the rigid surfaces. It must not be compared with the yield stress but with the maximum allowable surface pressure of the material (p all ). As well as the maximum stress σ z max acts orthogonally to the center of the contact ellipse, the other two main stresses also act, preventing lateral deformation. Consequently the equivalent stress due to Hertzian contact pressures (σ eq ) is lower than the σ z max .
As σ z max = 517.2 MPa and, according to [33], σ eq is only 63-65% σ z max : Experimental relations link the maximum allowable surface pressure p all to the material hardness (HB = 215 for AISI 316, see Table 1). One of these is [34]: or [35]: For AISI 316 p all is between 390 MPa and 550 MPa, higher than the equivalent stress σ eq = 325.8 MPa.
The simulation did not consider the second monopolar plate of the EHC, and infinitely rigid surfaces of washers constrained the top surface of the plate. The FEM simulation without top washers showed a 3.13% increase in deflection; hence the actual displacement along the z-axis will be within the range [0%, 3.13%].

Conclusions and Future Work
This work analyzed the monopolar plate of an Electrochemical Hydrogen Compressor. Small-sized channels and holes behave as notches, significantly increasing the stress triaxiality at the local scale. High-stress triaxiality peaks, eventually combined with the hydrogen embrittlement, can significantly increase the probability of failure and crack initiation at a given material point.
Based on such considerations, the Electrochemical Hydrogen Compressor End-Plate must be designed to have conservative levels of the nominal stress states.
The present work aimed to find the optimal shape of the monopolar plate, starting from three different geometric configurations. The results of the FEM analysis with Marc software-of MSC Software Corporation-identified the optimal configuration in circular plan view and active area. The plate in this configuration has been sized to have a deflection no greater than 0.1 mm when the EHC works at its operational pressure (700 bar). The minimum plate thickness is 17 mm. The deflection can increase up to 3.13% as the upper surface of the plate has been constrained by infinitely rigid washers rather than the deformable surface of the other monopolar plate.
Future developments foresee realizing more accurate FEM simulations both on the monopolar plate and on the bipolar plate, involving new aspects, such as implicit elastoplastic simulations or damage both on a global scale (loads acting on the system are not the operating ones but the exceptional ones) and on a local scale (plasticization occurs even with standard loads, where the cracks appear).