Start-up Operation Research Articles

Transient analysis of an electric vehicle air-conditioning system using CO2 for start-up and cabin pull-down operations

Abstract The energy consumption during the start-up and cabin pull-down operation of an air-conditioning system (A/C) for electric cars was studied in this work. An entire system model consisting of a CO2 vapor compression system (VCS), air handling unit (AHU), and the cabin was established utilizing an acausal, object-oriented, and equation-based modeling platform. A previously validated model was adopted for the VCS. AHU and cabin models were established by paying attention to the latent heat load, humidity, and dehumidification. The models were validated against the data from three different studies. A set of transient simulations was carried out for the start-up and cabin pull-down operation periods. It is found that the initial temperature of the cabin significantly influences the energy consumption and pull-down time. When the cabin target temperature decreases from 25 °C to 21 °C, the compressor consumes 58% more energy, and the pull-down time increases by 23%. Increasing the initial temperature 35 °C to 55 °C leads to 150% longer pull-down time, and 289% more energy consumption. The initial RH of the cabin air has a negligible impact on the pull-down time, but affects energy consumption. Increasing the RH from 33% to 66% results in 9% energy consumption.

Improved Carbon Corrosion and Platinum Dissolution Durability in Automotive Fuel Cell Startup and Shutdown Operation

Three protocols of accelerated startup and shutdown (SU/SD) test were investigated: startup and shutdown with air supply and soak to both anode and cathode (air-SU/SD), hydrogen protected startup and shutdown (H2-SU/SD), and hydrogen protected startup and shutdown with a load (H2-SU/SD with a load). The performance losses, electrochemical surface area (ECSA) reduction, and catalyst layer degradation were characterized and compared for these SU/SD protocols. Air-SU/SD protocol showed much more severe performance loss and catalyst layer degradation than hydrogen protected ones, which confirmed the benefits of hydrogen protection. The temperature effect on air-SU/SD was significant in a broad range from 20 °C to 70 °C, with low temperature greatly reducing the degradation. The mechanism of H2 protection and load drawn in alleviating carbon corrosion was explained based on reactions and charge conservation during SU/SD. This paper provides comprehensive test data and failure analysis to quantify the benefits of H2 protection and load drawn and to facilitate future enhancement of system strategies on SU/SD durability.

One-step start-up and subsequent operation of CANON process in a fixed-bed reactor by inoculating mixture of partial nitrification and Anammox sludge

The feasibility of one-step start-up of CANON process in a FBR by inoculating mixture of partial nitrification and Anammox sludge as well as its subsequent operation performances were investigated in the present study. The FBR was operated for around 3 months. The CANON process was quickly started up within 21 days. The max total nitrogen (TN) removal rate reached 183.61 g m−3 d−1 with the TN removal efficiency of 91.81% on day 95. The CANON process exhibited a good capability for resistance to loading shock and restoration from the unstable state. The mature CANON biofilms displayed a morphology of aggregates and had porous and microporous structure. The structural characteristics of the biofilms were conducive to improve the transferring of substrates and products. AOB and Anammox bacteria absolutely predominated in the mature biofilms and furthermore established a balanced interaction relationship. The microbial community structure contributed to the relatively stable operation performances.

Experimental and simulation study on preheating pulverized coal and air flow by a heat exchanger using thermal oil

Utility boilers consume a lot of fuel oil during the startup and low-load operation, which leads to the rise of power generation costs and environmental pollution problems. A finned shell-and-tube heat exchanger is proposed to preheat pulverized coal and air flow using thermal oil. Pilot-scale heat exchange tests were carried out to study the effects of the oil temperature, air temperature, air velocity, oil flow rate, coal-air ratio on the preheating performance and resistance properties of the system. Meanwhile, Computational fluid dynamic (CFD) numerical simulations were conducted. The tests verified the technical feasibility of the heat exchanger using thermal oil to preheat the primary air. The proposed preheating process is conducive to environmental protection and energy-saving. The test results show that when the coal-to-air ratio is 0.15, the 190 ℃ of thermal oil can preheat the pulverized coal and air flow from 65.2 ℃ to above 115 ℃. This preheating process will improve the ignition performance of the pulverized coal and air flow, which is favorable for startup and stable operation of the pulverized coal boiler, particularly when the boiler load is low.

Refrigerant Charge Prediction of Vapor Compression Air Conditioner Based on Start-Up Characteristics

Refrigerant charge faults, which occur frequently, increase the energy loss and may fatally damage the system. Refrigerant leakage is difficult to detect and diagnose until the fault has reached a severe degree. Various techniques have been developed to predict the refrigerant charge amount based on steady-state operation; however, steady-state experiments used to develop prediction models for the refrigerant charge amount are expensive and time-consuming. In this study, a prediction model was established with dynamic experimental data to overcome these deficiencies. The dynamic models for the condensation temperature, degree of subcooling, compressor discharge temperature, and power consumption were developed with a regression support vector machine (r-SVM) model and start-up experimental data. The dynamic models for the condensation temperature and degree of subcooling can predict the distinct start-up characteristics depending on the refrigerant charge amount. Moreover, the estimated root mean square error (RMSE) of the condensation temperature and degree of subcooling of the test data are 0.53 and 0.84 °C, respectively. The refrigerant charge is one of the predictors that defines the dynamic characteristics. The refrigerant charge can be estimated by minimizing the RMSE of the predicted values of the dynamic models and experimental data. When the dynamic characteristics of the two predictor variables, “condensation temperature” and “degree of subcooling” are used together, the average prediction error of the test data is 2.54%. The proposed method, which uses the dynamic model during start-up operation, is an effective technique for predicting the refrigerant charge amount.

Development of a feedback control system for a differential biofilter degrading toluene contaminated air

In this study, a proportional – integral feedback control system was implemented on a lab-scale differential biofilter to control the gas phase toluene concentration in the soil bed through online manipulation of the inlet toluene concentration. The feedback control system was based on a cascade controller that manipulated the setpoint of an air bath diffusion system to manipulate the inlet toluene concentration. The controller performed well for toluene concentrations in the reactor of 10 – 300 ppm for both setpoint changes and disturbance rejections; however, the system was nonlinear requiring different tuning parameters at different outlet concentrations. Feedback control of the toluene concentration in the differential reactor was used to explore the impact of concentration on start-up and long-term biofilter operation in a rigorous fashion. Starting at an reactor concentration of 20 ppm and then increasing to 65 ppm increased the toluene removal rate (33 ± 1.6 g m−3h−1) compared to starting the reactor at an outlet concentration of 81 ppm before settling at 65 ppm (42 ± 0.9 g m− 3h−1). The toluene removal rate increased with increasing outlet toluene concentration and then eventually decreased when reaching the inhibitory toluene concentration (ranged from 80 to 250 ppm).

1-cycle RMS를 이용한 독립계통용 디젤발전기의 자동전압제어장치 성능 개선에 관한 연구

Since the power system in an island area is generally small-scale and not connected to the main land, it should be able to operate independently. And the power is mainly supplied by diesel generators. Under these circumstances, the frequency variation occurs easily when a variable load with a large capacity is operated during the diesel generator start-up or stand-alone operation. So, AVR(Automatic Voltage Regulator) unit generally used in the main land, which is based on 60 Hz fundamental frequency DFT(Discrete Fourier Transform) algorithm, cannot control the voltage reliably and immediately due to the error of voltage effective value caused by the system frequency fluctuation. In this study, the performance of the AVR unit with 1-cycle RMS(Root Mean Square) algorithm, which enables stable voltage control despite frequency variation, was developed and demonstrated through field application tests.

Basic concept of JA DEMO fuel cycle

In the JA DEMO, the deuterium-tritium (DT) fuel will be supplied by a combination of gas puff, pellet injection and neutral beam injection. Since the combustion rate of fuel in the JA DEMO is approximately 1.7 %, it is necessary to establish the closed deuterium-tritium fuel cycle in the facility for fuel recycling. The fuel cycle for DEMO inevitably differs from that for ITER. Compared to the ITER fuel cycle, which requires high flexibility for fuel supply based on the requirements of plasma experiments, the flexibility required for the fuel cycle decreases in the DEMO reactor. The flow in the fuel cycle becomes steadier. Therefore, the fuel cycle of the DEMO shifts to the continuous processing system. In addition, considering the tritium breeding in the blanket, the fuel cycle must be established so that tritium in the DEMO facility can be balanced. The process components of the fuel cycle that process a large amount of tritium-containing gas are provided with multiple barriers that confine tritium, and the detritiation system processes the confined tritium. The confinement of tritium in the DEMO facility according to nuclear safety regulations is an important safety issue. As a result, the fuel cycle becomes a complex chemical plant. Assuming safety requirements, it is necessary to consider the reduction of the tritium inventory as much as possible at the design stage of the DEMO fuel cycle. The tritium inventory in the fuel cycle is closely related to the requirements to pellet manufacturing and hydrogen isotope separation, and the reduction of the requirements to the hydrogen isotope separation is the major issue in the fuel cycle. For this purpose, it is necessary to minimize the D/T separation based on the concept of direct fuel cycle with the D/T mixture. For the possible shortage in the preparation of the initial loading tritium for the DEMO reactor, it is necessary to design the fuel cycle through dynamic numerical simulation so that the fuel cycle can be started with the minimum initial loading tritium. In the system design of the fuel cycle, it is necessary to support various operation modes from the startup of DT operation to steady operation. The requirements to design the JA DEMO fuel cycle are being investigated in the JA DEMO design activity.

Characterization of commercial vehicles’ start-up operations from in-use data

Diesel engines produce disproportionate levels of emissions when the engine and after-treatment systems are operating at low temperatures. This situation arises most commonly when the vehicle is first started after overnight. To quantify emissions attributable to vehicle starts, a sizable collection of on-road commercial vehicle operating data is analyzed to identify start-up events and inspect associated emissions. Data was obtained from the National Renewable Energy Laboratory’s (NREL’s) Fleet DNA and from the Center for Environmental Research & Technology (CE-CERT). Included are 500 + diesel vehicles with more than 42,000 recorded days, drawn from 25 vocational categories across the United States. Analysis shows that vehicle behavior, in terms of engine cold- and warm-operation, starts per day, soak time, and warm-up duration, differs significantly between vehicle vocations. Also, weighting factors for cold- and hot-starts currently used in the U.S. Environmental Protection Agency’s Federal Test Procedure (FTP) for heavy-duty emissions certification accurately represent real-world operations. Although the FTP includes a comparable fraction of cold operation, the hot fraction is much shorter than real-world operation due to limited test duration. The investigation also revealed that real-world engines operate for a significant amount of time when the engine coolant is in the “hot-stabilized” region, but the selective catalytic reduction (SCR) temperature is below its effective operating temperature of 200 °C. Of the vehicles under investigation, almost 20% of their operational time is within this condition. Thus, novel approaches to raise and maintain SCR temperature are highly required to further reduce engine emissions.

Combustion and emission characteristics of a reprocessed used lubricating oil as a renewable fuel for boiler cold start-up operation

This paper investigated the combustion and emission characteristics of a reprocessed used lubricating oil (RULO) with a pilot combustion test furnace to assess the possibility for using it as a fuel for boiler cold start-up operation. Results were compared with those of a by-product fuel oil (BFO) presently used in the power plant. Fuel analysis and testing showed that both oils had similar fuel properties, but RULO represented less reactivity and larger sprayed particle diameter than BFO. The cold start-up combustion tests revealed that RULO was ignited later and furnace exit temperatures was higher in RULO combustion, compared to those in BFO combustion. NOx emission in RULO combustion was about 10 ppm higher due to the higher fuel-N, and particulate emissions, thought to be agglomerated unburned carbonaceous particles, were severer in RULO combustion test than BFO combustion test. The relative comparisons between two fuel oils would provide a basic information for substituting an existing fuel oil to a renewable fuel oil during the start-up operation of the coal-fired power plant.

Unsteady flow simulation of industrial gas turbine compressor under start-up operation

Unsteady flow simulation of industrial gas turbine compressor under start-up operation

Розробка раціональних технологічних параметрів обігріву коксової батареї з трамбуванням вугільної шихти.

DEVELOPMENT OF RATIONAL TECHNOLOGICAL PARAMETERS FOR HEATING OF A COKE BATTERY WITH A STAMPED CHARGING OF COALS © V.I. Goncharov, I.I. Sikan, Ya.I. Dyachuk, N.V. Mukina (Coke production of PJSC «ArcelorMittal Kryvyi Rih», 50095, Dnipropetrovsk region, Kryvyi Rih, Kryvorizhstal str., 1, Ukraine), I.V. Shulga, PhD in technical sciences (State Enterprise "Ukrainian State Research Coal Chemical Institute (UHIN)", 61023, Kharkov, Vesnina st., 7, Ukraine) The article is devoted to the introduction at the PJSC «ArcelorMittal Kryvyi Rih» coke plant of a promising technology for compaction of a coal charge, which allows obtaining blast furnace coke of high mechanical and "hot" strength from charges characterized by an increased content of lean components of various stages of metamorphism. Under the present circumstances, this is of great importance for increasing the economic indicators of the production of blast furnace coke and the operation of blast furnaces. The start-up of coke oven battery No. 6, designed for the application of this technology, necessitated the determination of rational technological parameters of heating: the temperature level in the control verticals and the location of the burners for supplying coke oven gas in the vertical along the heating wall. The article presents the calculations, the results of which confirm the expediency of rationalizing the actual starting (design) distribution of burners on coke oven battery No. 6. First of all, attention is drawn to the insufficient diameters of the holes in the burners for the extreme and pre-extreme verticals. This, according to the authors, caused underheating of the head zones and resin spills, which were observed during start-up operations on the battery. The replacement of burners at the edges of the walls by large diameters, carried out by the personnel of the coke shop of the coke plant PJSC “ArcelorMittal Kryvyi Rih” and LLC “Koksokhimstanciya”, combined with the full opening of the recirculation windows in these zones, improved the heating of the charge and significantly reduced the intensity of spills. The experience gained was taken into account when locating the control devices at coke oven battery No. 5, which was put into operation later. The article also formulates recommendations for adjusting the temperature regime of coking, depending on the grade and component composition, as well as the technological properties of coal charges. Keywords: coke oven battery, stumping of coal charge, temperature in control verticals, location of burners along the length of the heating wall. Corresponding author Ya.I. Dyachuk, е-mail: Yaroslav.Dyachuk@arcelormittal.com

Stress Calculation and Analysis of a 600MW Supercritical Boiler’s Steam-water Separator

The steam-water separator of a 600MW boiler has been chosen as the research object. By theoretical analysis’s method and finite element analysis, the stress and secure research has been carried out in the steam-water separator as the thick-walled pressure vessel. For its start-up and normal operation of the working process and features and considering the internal pressure load and thermal stress loading conditions, the initiation of separation stress and security analysis have been carried out about the steam-water separator. The primary stress at the continuum of the cylinder structure can be calculated by the formula, and the calculated results are basically in agreement with the simulated graphic results. And the research work provides a theoretical basis for on-line stress monitoring and life assessment.

Mathematical modelling of start-up transients at clustered propulsion system with POGO-suppressors for CYCLON-4M launch vehicle

Liquid-propellant rocket propulsion systems of the first stages of launch vehicles of medium, heavy, and super-heavy class usually include POGO-suppressors, which are one of the most widely used methods to eliminate launch vehicle longitudinal structural vibrations (POGO phenomena). However, until now, the theoretical studies and analysis of the effect of the POGO-suppressors’ installation in the feedlines of main liquid rocket engines on transient processes in systems during rocket engine starting have not been carried out due to the complexity of such analysis and the lack, first of all, reliable nonlinear models of cavitation phenomena in rocket engine pumps. A mathematical model for the start-up of a clustered rocket propulsion of the Cyclone-4M launch vehicle has been developed that takes into account the low-frequency dynamics of the POGO-suppressors and the asynchronous start-up timeline sequences of the rocket engines. The first stage of the launch vehicle propulsion system includes four RD-870 rocket engines. A nonlinear mathematical model of low-frequency dynamic processes of the POGO-suppressor with bellows separation of liquid and gaseous media is presented. A significant effect of cavitation in the pumps of engines and the POGO-suppressor installation to the LOX feedline on the propulsion system dynamic gains is shown. Based on the developed mathematical model of the clustered rocket propulsion start-up, the studies of the Cyclone-4M main engines’ start-up transients were carried out. The asynchronous start-up timeline sequences of the rocket engine and the places of installation of the POGO-suppressors in the LOX feedline branches to the RD-870 rocket engine – near the general feedline collector as standard placement or directly at the entrance to the engines – were investigated. The analysis of start-up transients in the oxidizer feed system of the considered propulsion (the time dependences of the flowrate and pressure at the engine inlet) showed the following. Firstly, while the synchronous start-up of the engines, the installation of the POGO-suppressors near the feedline collector makes it possible to eliminate all engine inlet overpressures that exist in the rocket propulsion system in case of the absence of the POGO-suppressors. Secondly, the RD-870 engine asynchronous start-up operation affects negatively the time dependences of the propellant flowrate and pressure at the engine inlet if the POGO-suppressors are located near the feedline collector. So, in the propulsion system’s start-up timeline interval 0.95 s - 1.35 s, for some computational variants of the initial moments of the engine operation start, an abnormally large drop in the LOX flow rate and the overpressures at the engine inlet is observed. The asynchronous start-up of the RD-870 engines with the installation of the POGO-suppressors at the engine inlet does not significantly change the start-up transients compared to the synchronous starting of the engines. Thirdly, thus, it is shown that the installation of the POGO-suppressors both at the engine inlet and at the RD-870 branches near the collector has a significant positive effect on the quality of start-up transient processes for the main engines of the 1st stage of the Cyclone-4M launch vehicle. Placing the POGO-suppressors at the engine inlets is not standard and is considered without reference to the propulsion system layout. Nevertheless, the POGO-suppressors installed at the inlet to the engines are an effective means of preventing overshoots and dips in the parameters of the liquid-propellant rocket engine, including the conditions of asynchronous starting of the liquid rocket engines in the clustered propulsion system. The results obtained can be used in mathematical modeling of the start-up of the first stage propulsion system either for multistage sustainer rockets used in parallel with booster rockets or for the clustered multi-engine rocket propulsion system containing POGO-suppressors.

Design Adaptive controller and learning controller to control the pacemaker: تصميم متحكم تكيفي ومتحكم متعلم للتحكم بناظم الخطى القلبي

Pacemakers play a vital role in the management of electrical disorders of the heart, as having a controller that presents a signal that generates a correct heart rate without delay or overshoot is important, to provide better quality and long-term patient's life. Many advances in pacemaker control technology have been achieved, making today's pacemakers capable of correcting many electrical heart abnormalities. In this paper, a fuzzy model reference learning control was designed to control the pacemaker, and compared with a controller adopts the adaptive mechanism instead of the learning mechanism, in order to confirm the importance of learning in the controller. The controller automatically controls the heart rate (Heart Rate: HR) to prevent it from falling under a reference value and makes it track preset signal accurately. The learning controller achieved better performance than the adaptive controller and the learning controller achieved when comparing with the similar intelligent ones that proposed to control pacemaker the required response with a small settling time 0.9s, and minimized the overshoot at start-up operation where overshoot was 0.2%

Human error analysis in furnace start-up operation using HEART under intuitionistic fuzzy environment

To address human error in system reliability, Human Reliability Analysis (HRA) is an essential issue. Human Error Assessment and Reduction Technique (HEART) as a rather straightforward technique for HRA has successfully been used in many areas to predict human error probability (HEP). However, knowledge acquisition of experts during assessed proportion of affect (APOA) calculation is subjected to vagueness and ambiguity. To overcome this challenge, in this paper Intuitionistic Fuzzy (IF) set due to their advantage to represent more fuzzy information than a classical fuzzy set adopted through APOA calculation. To demonstrate this hybrid approach short for, IF-HEART, the furnace start-up operation is handled, since analysis shows that most of explosions and losses occur during furnace start-ups operation. Further, a sensitivity analysis is carried out to approve the proposed integrated approach. In addition to its academic contribution, the results of the paper enable to improve the overall safety level of a furnace by taking into account potential human error.

Reliability analysis for a multi-stack solid oxide fuel cell system subject to operation condition-dependent degradation

PurposeA reliability analysis of a solid oxide fuel cell (SOFC) system is presented for applications with strict constant power supply requirements, such as data centers. The purpose is to demonstrate the effect when moving from a module-level to a system-level in terms of reliability, also considering effects during start-up and degradation.Design/methodology/approachIn-house experimental data on a system-level are used to capture the behavior during start-up and normal operation, including drifts of the operation point due to degradation. The system is assumed to allow replacement of stacks during operation, but a minimum number of stacks in operation is needed to avoid complete shutdown. Experimental data are used in conjunction with a physics-based performance model to construct the failure probability function. A dynamic program then solves the optimization problem in terms of time and replacement requirements to minimize the total negative deviation from a given target reliability.FindingsResults show that multi-stack SOFC systems face challenges which are only revealed on a system- and not on a module-level. The main finding is that the reliability of multi-stack SOFC systems is not sufficient to serve as sole power source for critical applications such as data center.Practical implicationsThe principal methodology may be applicable to other modular systems which include multiple critical components (of the same kind). These systems comprise other electrochemical systems such as further fuel cell types.Originality/valueThe novelty of this work is the combination of mathematical modeling to solve a real-world problem, rather than assuming idealized input which lead to more benign system conditions. Furthermore, the necessity to use a mathematical model, which captures sufficient physics of the SOFC system as well as stochasticity elements of its environment, is of critical importance. Some simplifications are, however, necessary because the use of a detailed model directly in the dynamic program would have led to a combinatorial explosion of the numerical solution space.

High frequency signal injection method for sensorless permanent magnet synchronous motor drives

The objective of this project is to design a high frequency signal injection method for sensorless control of permanent magnet synchronous motor (PMSM) drives. Generally, the PMSM drives control requires the appearance of speed and positon sensor to measure the motor speed hence to feedback the information for variable speed drives operation. The usage of the sensor will increase the size, cost, extra hardwire and feedback devices. Therefore, there is motivation to eliminate this type of sensor by injecting high frequency signal and utilizing the electrical parameter from the motor so that the speed and positon of rotor can be estimated. The proposed position and speed sensorless control method using high frequency signal injection together with all the power electronic circuit are modelled using Simulink. PMSM sensorless driveis simulated and the results are analyzed in terms of speed, torque and stator current response without load disturbance but under the specification of varying speed, forward to reverse operation, reverse to forward operation and step change in reference speed. The results show that the signal injection method performs well during start-up and low speed operation.

Economic analysis of a small combined cycle gas turbine power plant in the Italian electricity market

Despite the availability of new alternative energy sources, growing worldwide energy demand and emissions targets lead power plants to work optimizing performances. In this new scenario in which renewable energies are increasingly taking the field, it is also important to produce energy at a low price. Moreover, the variability of the energy market price complicates the analysis. Comparison between the produced energy cost and the market price is necessary to get a return on investment. In this context, the paper investigates the implications of using a gas turbine in an energy system estimating the plant layout and the number of working hours that guarantees a better profit. The analysis focuses on the study of the start-up and shutdown operation mode to find the optimal solution strategy in the Italian electricity market.

Pore Network Modeling of Microporous Layers for Polymer Electrolyte Fuel Cells

Liquid water accumulated in the catalyst layer (CL) and gas diffusion layer (GDL) of a polymer electrolyte fuel cell (PEFC, PEMFC) results in performance deterioration due to inhibition of oxygen transport to the cathode CL (Flooding). To enhance the drainage in the GDL, the application of a microporous layer (MPL) is effective to the CL side of the GDL substrate [1, 2]. To elucidate the effects of two-phase flow on the oxygen transport in an MPL for its optimized designs, we model three dimensional porous structures of in-house MPLs with pore network model (PNM) [3, 4] for convective air permeation and oxygen diffusion. Pore diameter distribution derived by a focused ion beam scanning electron microscope (FIB-SEM) is employed in the PNM. Air permeation measurements combined with oxygen diffusion measurements with changing wetting liquid saturation (Galwick, Porous Materials Inc., USA) by using gas chromatography validate the model.AcknowledgmentsThe authors are indebted to graduate students, Messrs. Yuhang Liu, Kentaro Harano, Dingfeng Chen for help in the modeling and experiment. Thanks are offered to Professor Kohei Ito of Kyushu University for valuable discussions. The authors acknowledge Ms. Chie Uryu of International Research Center for Hydrogen Energy, Kyushu University for operating the FIB-SEM.References Nakajima H, Konomi T and Kitahara T. 2007. Direct water balance analysis on a polymer electrolyte fuel cell (PEFC): Effects of hydrophobic treatment and micro-porous layer addition to the gas diffusion layer of a PEFC on its performance during a simulated start-up operation. Journal of Power Sources 171: 457-463. Kitahara T, Konomi T and Nakajima H. 2010. Microporous layer coated gas diffusion layers for enhanced performance of polymer electrolyte fuel cells. Journal of Power Sources 195: 2202-2211. Gostick J, Aghighi M, Hinebaugh J, Tranter T, Hoeh MA, Day H, Spellacy B, Sharqawy MH, Bazylak A, Burns A, Lehnert W and Putz A. 2016. OpenPNM: A pore network modeling package. Computing in Science & Engineering 18: 60-74. http://openpnm.org/. Accessed on September 14, 2020