Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed model predictive control”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

A coupled vacancy diffusion-dislocation dynamics model for the climb-glide motion of jogged screw dislocations

Here, we develop a novel model to study the climb/glide motion of jogged screw dislocations within the discrete dislocation dynamics (DDD) framework. We present results for the dependence of the climb velocity on the applied stress and on the jog size and distribution statistics. We show that the model predictions are consistent with experimental data in both γ-TiAl and Zircaloy-4. The ranges of the applied stress and jog spacing that determine the dominance of one of three dislocation mechanisms are identified. These are the jog dragging, dipole dragging, and dipole bypass mechanisms, respectively. The overall dislocation motion in the jog dragging regime is composed of glide of screw segments and climb of jogs, controlling the plastic strain and the creep rate, respectively. Based on current simulations and on a detailed examination of the predicted jog heights compared to experiments, we advance the hypothesis that a combination of jog dragging and dipole bypass mechanisms is necessary to reproduce the high creep rate observed in some experiments.

36 MATERIALS SCIENCE↗

Predicting Outcomes of Nanoparticle Attachment by Connecting Atomistic, Interfacial, Particle, and Aggregate Scales

Predicting nanoparticle aggregation and attachment phenomena requires a rigorous understanding of the interplay among crystal structure, particle morphology, surface chemistry, solution conditions, and interparticle forces, yet no comprehensive picture exists. We used an integrated suite of experimental, theoretical, and simulation methods to resolve the effect of solution pH on the aggregation of boehmite nanoplatelets, a case study with important implications for the environmental management of legacy nuclear waste. Real-time observations showed that the particles attach preferentially along the (010) planes at pH 8.5 and the (101) planes at pH 11. To rationalize these results, we established the connection between key physicochemical phenomena across the relevant length scales. Starting from molecular-scale simulations of surface hydroxyl reactivity, we developed an interfacial-scale model of the corresponding electrostatic potentials, with subsequent particle-scale calculations of the resulting driving forces allowing successful prediction of the attachment modes. Finally, we scaled these phenomena to understand the collective structure at the aggregate-scale . Here, our results indicate that facet-specific differences in surface chemistry produce heterogeneous surface charge distributions that are coupled to particle anisotropy and shape-dependent hydrodynamic forces, to play a key role in controlling aggregation behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of macromolecular architecture on necking in polymer extrusion film casting process

Extrusion film casting (EFC) is an important polymer processing technique that is used to produce several thousand tons of polymer films/coatings on an industrial scale. In this research, we are interested in understanding quantitatively how macromolecular chain architecture (for example long chain branching (LCB) or molecular weight distribution (MWD or PDI)) influences the necking and thickness distribution of extrusion cast films. We have used different polymer resins of linear and branched molecular architecture to produce extrusion cast films under controlled experimental conditions. The necking profiles of the films were imaged and the velocity profiles during EFC were monitored using particle tracking velocimetry (PTV) technique. Additionally, the temperature profiles were captured using an IR thermography and thickness profiles were calculated. The experimental results are compared with predictions of one-dimensional flow model of Silagy et al{sup 1} wherein the polymer resin rheology is modeled using molecular constitutive equations such as the Rolie-Poly (RP) and extended Pom Pom (XPP). We demonstrate that the 1-D flow model containing the molecular constitutive equations provides new insights into the role of macromolecular chain architecture on film necking.{sup 1}D. Silagy, Y. Demay, and J-F. Agassant, Polym. Eng. Sci., 36, 2614 (1996)

36 MATERIALS SCIENCE↗

Modeling the effect of neighboring grains on twin growth in HCP polycrystals

In this paper, we study the dependence of neighboring grain orientation on the local stress state around a deformation twin in a hexagonal close packed (HCP) crystal and its effects on the resistance against twin thickening. We use a recently developed, full-field elasto-visco-plastic formulation based on fast Fourier transforms that accounts for the twinning shear transformation imposed by the twin lamella. The study is applied to Mg, Zr and Ti, since these HCP metals tend to deform by activation of different types of slip modes. The analysis shows that the local stress along the twin boundary are strongly controlled by the relative orientation of the easiest deformation modes in the neighboring grain with respect to the twin lamella in the parent grain. A geometric expression that captures this parent-neighbor relationship is proposed and incorporated into a larger scale, mean-field visco-plastic self-consistent model to simulate the role of neighboring grain orientation on twin thickening. We demonstrate that the approach improves the prediction of twin area fraction distribution when compared with experimental observations.

36 MATERIALS SCIENCE↗

SU-C-BRE-02: BED Vs. Local Control: Radiobiological Effect of Tumor Volume in Monte Carlo (MC) Lung SBRT Planning

Purpose: SBRT with hypofractionated dose schemata has emerged a compelling treatment modality for medically inoperable early stage lung cancer patients. It requires more accurate dose calculation and treatment delivery technique. This report presents the relationship between tumor control probability(TCP) and size-adjusted biological effective dose(sBED) of tumor volume for MC lung SBRT patients. Methods: Fifteen patients who were treated with MC-based lung SBRT to 50Gy in 5 fractions to PTVV100%=95% were studied. ITVs were delineated on MIP images of 4DCT-scans. PTVs diameter(ITV+5mm margins) ranged from 2.7–4.9cm (mean 3.7cm). Plans were generated using non-coplanar conformal arcs/beams using iPlan XVMC algorithm (BrainLABiPlan ver.4.1.2) for Novalis-TX with HD-MLCs and 6MVSRS(1000MU/min) mode, following RTOG-0813 dosimetric guidelines. To understand the known uncertainties of conventional heterogeneities-corrected/uncorrected pencil beam (PBhete/ PB-homo) algorithms, dose distributions were re-calculated with PBhete/ PB-homo using same beam configurations, MLCs and monitor units. Biologically effective dose(BED10) was computed using LQ-model with α/β=10Gy for meanPTV and meanITV. BED10-c*L, gave size-adjusted BED(sBED), where c=10Gy/cm and L=PTV diameter in centimeter. The TCP model was adopted from Ohri et al.(IJROBP, 2012): TCP = exp[sBEDTCD50]/ k /(1.0 + exp[sBED-TCD50]/k), where k=31Gy corresponding to TCD50=0Gy; and more realistic MC-based TCP was computed for PTV(V99%). Results: Mean PTV PB-hete TCP value was 6% higher, but, mean PTV PB-homo TCP value was 4% lower compared to mean PTV MC TCP. Mean ITV PB-hete/PB-homo TCP values were comparable (within ±3.0%) to mean ITV MC TCP. The mean PTV(V99%)had BED10=90.9±3.7%(median=92.2%),sBED=54.1±8.2%(median=53.5%) corresponding to mean MC TCP value of 84.8±3.3%(median=84.9%) at 2- year local control. Conclusion: The TCP model which incorporates BED10 and tumor diameter indicates that radiobiological effect of target volume and dose calculation algorithm significantly affects TCP for lung SBRT patients. Dose calculation using MC-based algorithm is more realistic with tissue heterogeneities and is routinely performed in our clinic. Patients will be followed up to determine whether TCP prediction correlate clinical outcomes.

60 APPLIED LIFE SCIENCES↗

Priors on red galaxy stochasticity from hybrid effective field theory

ABSTRACT We investigate the stochastic properties of typical red galaxy samples in a controlled numerical environment. We use halo occupation distribution (HOD) modelling to create mock realizations of three separate bright red galaxy samples consistent with data sets used for clustering and lensing analyses in modern galaxy surveys. Second-order Hybrid Effective Field Theory (HEFT) is used as a field-level forward model to describe the full statistical distribution of these tracer samples, and their stochastic power spectra are directly measured and compared to the Poisson shot-noise prediction. While all of the galaxy samples we consider are hosted within haloes with sub-Poisson stochasticity, we observe that the galaxy samples themselves possess stochasticities that range from sub-Poisson to super-Poisson, in agreement with predictions from the halo model. As an application of our methodology, we place priors on the expected degree of non-Poisson stochasticity in cosmological analyses using such samples. We expect these priors will be useful in reducing the complexity of the full parameter space for future analyses using second-order Lagrangian bias models. More generally, the techniques outlined here present the first application of HEFT methods to characterize models of the galaxy–halo connection at the field level, revealing new connections between once-disparate modelling frameworks.

79 ASTRONOMY AND ASTROPHYSICS↗

Intercomparison of Large-eddy Simulations of Arctic Mixed-phase Clouds: Importance of Ice Size Distribution Assumptions

Large-eddy simulations of mixed-phase Arctic clouds by 11 different models are analyzed with the goal of improving understanding and model representation of processes controlling the evolution of these clouds. In a case based on observations from the Indirect and Semi-Direct Aerosol Campaign (ISDAC), it is found that ice number concentration, Ni, exerts significant influence on the cloud structure. Increasing Ni leads to a substantial reduction in liquid water path (LWP) and potential cloud dissipation, in agreement with earlier studies. By comparing simulations with the same microphysics coupled to different dynamical cores as well as the same dynamics coupled to different microphysics schemes, it is found that the ice water path (IWP) is mainly controlled by ice microphysics, while the inter-model differences in LWP are largely driven by physics and numerics of the dynamical cores. In contrast to previous intercomparisons, all models here use the same ice particle properties (i.e., mass-size, mass-fall speed, and mass-capacitance relationships) and a common radiation parameterization. The constrained setup exposes the importance of ice particle size distributions (PSD) in influencing cloud evolution. A clear separation in LWP and IWP predicted by models with bin and bulk microphysical treatments is documented and attributed primarily to the assumed shape of ice PSD used in bulk schemes. Compared to the bin schemes that explicitly predict the PSD, schemes assuming exponential ice PSD underestimate ice growth by vapor deposition and overestimate mass-weighted fall speed leading to an underprediction of IWP by a factor of two in the considered case.

Arctic↗

Hydrologic applicability of satellite-based precipitation estimates for irrigation water management in the data-scarce region

Reliable precipitation estimates are crucial for planning and managing water resources, monitoring hydrologic extremes, and fulfilling irrigation water requirements. Accurate precipitation estimates are particularly challenging in complex mountain terrains, where monitoring gauges are often sparsely distributed due to their remote locations, and high installation and long-term operation costs. Recent advances in satellite-based precipitation estimates offer promising opportunities to improve our understanding of hydrologic processes and their applications for irrigation water management. Several datasets are available varying considerably in terms of their data sources, quality control methods, estimation procedure, and spatiotemporal resolutions. Choosing the most suitable dataset for a particular application is a complex task. In this study, we (1) evaluate the performance of six satellite-based precipitation estimates (SPEs): i) CHIRPS v2.0, ii) CMORPH v1.0, iii) ERA5, iv) IMERG v6, v) MSWEP v2.8, and vi) PERSIANN-CDR against the gauge precipitation using continuous statistical and categorical indices, (2) integrate SPEs with a calibrated semi-distributed hydrologic model to predict streamflow, and (3) demonstrate practical implications of improved streamflow prediction for irrigation water management in the central Himalayan region, Nepal. Our results illustrate that satellite-based precipitation estimates have competitive performance in capturing a wide range of rainfall characteristics, with demonstrated variability across river basins and time scales. Further, there are no significant discrepancies observed in satellite-based precipitation estimates for estimating irrigation water requirements for the three major crops (maize, wheat, and paddy) during the cropping period across the selected river basins, showing a greater promise for irrigation water management planning and decision making.

54 ENVIRONMENTAL SCIENCES↗

Radio AGN in nearby dwarf galaxies: the important role of AGN in dwarf galaxy evolution

ABSTRACT We combine deep optical and radio data, from the Hyper Suprime-Cam and the Low-Frequency Array (LOFAR), respectively, to study 78 radio active galactic nuclei (AGN) in nearby (z < 0.5) dwarf galaxies. Comparison to a control sample, matched in stellar mass and redshift, indicates that the AGN and controls reside in similar environments, show similar star formation rates (which trace gas availability) and exhibit a comparable incidence of tidal features (which indicate recent interactions). We explore the AGN properties by combining the predicted gas conditions in dwarfs from a cosmological hydrodynamical simulation with a Monte Carlo suite of simulated radio sources, based on a semi-analytical model for radio-galaxy evolution. In the subset of LOFAR-detectable simulated sources, which have a similar distribution of radio luminosities as our observed AGN, the median jet powers, ages, and accretion rates are ∼1035 W, ∼5 Myr, and ∼10−3.4 M⊙ yr−1, respectively. The median mechanical energy output of these sources is ∼100 times larger than the median binding energy expected in dwarf gas reservoirs, making AGN feedback plausible. Since special circumstances (in terms of environment, gas availability, and interactions) are not necessary for the presence of AGN, and the central gas masses are predicted to be an order of magnitude larger than that required to fuel the AGN, AGN triggering in dwarfs is likely to be stochastic and a common phenomenon. Together with the plausibility of energetic feedback, this suggests that AGN could be important drivers of dwarf galaxy evolution, as is the case in massive galaxies.

Davis, F. (ORCID:0000000285224980)↗

Non-growing Season Plant Nutrient Uptake Controls Arctic Tundra Vegetation Composition: Modeling Archive

This Modeling Archive is in support of a NGEE-Arctic publication: Riley et al. (2021) “Non-growing season plant nutrient uptake controls Arctic tundra vegetation composition under future climate” (Riley et al., 2021). The dataset contains “ecosys” model outputs reported in Riley et al. (2021) that examines the role of non-growing season (NGS) nutrient dynamics on 21st century vegetation composition. The study highlighted the importance of these nutrient dynamics on the emergent plant functional type distributions, focusing on their role in shrub expansion. The study showed that ignoring NGS nutrient dynamics led to a dramatic under-prediction of shrub expansion by year 2100, motivating the need to include these processes in Earth System Models. Included are modeled net primary production (NPP) and nitrogen (N) uptake for two scenarios: (1) baseline and (2) no NGS N uptake. The output is provided (1) at 25 km resolution across the North American tundra; (2) for years 2012 through 2100; and (3) PFT-specific. Dataset contains 20 NetCDF files in a zip file and one pdf.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

3D Printed Parahydrophobic Surfaces as Multireaction Platforms

We report Parahydrophobic surfaces (PHSs) composed of arrays of cubic μ-pillars with a double scale of roughness and variable wettability were systematically obtained in one step and a widely accessible stereolithographic Formlabs 3D printer. The wettability control was achieved by combining the geometrical parameters (H = height and P = pitch) and the surface modification with fluoroalkyl silane compounds. Homogeneous distribution of F and Si atoms onto the pillars was observed by XPS and SEM-EDAX. A nano-roughness on the heads of the pillars was achieved without any post-treatment. The smallest P values lead to surfaces with static contact angles (CAs) >150° regardless of the H utilized. Interestingly, the relationship 0.6 ≤ H/P ≤ 2.6 obtained here was in good agreement with the H/P values reported for nano- and submicron pillars. Furthermore, experimental CAs, advancing and receding CAs, were consistent with the theoretical prediction from the Cassie–Baxter model. Structures covered with perfluorodecyltriethoxysilane with high H and short P lead to PHSs. Conversely, structures covered with perfluorodecyltrimethoxysilane exhibited a superhydrophobic behavior. Finally, several aqueous reactions, such as precipitation, coordination complex, and nanoparticle synthesis, were carried out by placing the reactive agents as microdroplets on the parahydrophobic pillars, demonstrating the potential application as chemical multi-reaction array platforms for a large variety of relevant fields in microdroplet manipulation, microfluidics systems, and health monitoring, among others.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗

Efficient Distance-based Global Sensitivity Analysis for Terrestrial Ecosystem Modeling

Sensitivity analysis in terrestrial ecosystem modeling is important for understanding controlling processes, guiding model development, and targeting new observations to reduce parameter and prediction uncertainty. Complex and computationally expensive terrestrial ecosystem models (TEM) limit the number of ensemble simulations, requiring sophisticated and efficient methods to analyze sensitivities of multiple model responses to different types of parameter uncertainties. In this study, we propose a distance-based global sensitivity analysis (DGSA) method. DGSA first classifies model response samples into a small set of discrete classes and then calculates the distance between parameter frequency distributions in different classes to measure the parameter sensitivity. The principle is that, if the parameter distribution is the same in each class, then the model response is insensitive to the parameter, while a large difference in the distributions indicates the parameter is influential to the response. Built on this idea, DGSA can be applied to analyze sensitivity of a single and a group of responses to different kinds of parameter uncertainties including continuous, discrete and even stochastic. Besides the main-effect sensitivity from a single parameter, DGSA can also quantify the sensitivity from parameter interactions. Additionally, DGSA is computationally efficient which can use a small number of model evaluations to obtain an accurate and statistically significant result. We applied DGSA to two TEMs, one having eight parameters and three kinds of model responses, and the other having 47 parameters and a long-period response. We demonstrated that DGSA can be used for sensitivity problems with multiple responses and high-dimensional parameters efficiently.

Lu, Dan↗

On the Low Risk of SSR in Type III Wind Turbines Operating in Grid-Forming Control: Preprint

We have shown in a previous work that the risk of sub-synchronous resonance (SSR) between a wind power plant with Type III wind turbines and series-compensated transmission lines is low when the wind turbines in the plant are operated in grid-forming mode, instead of the standard grid-following mode. The fundamental mechanism behind the improved damping characteristics is explained in this paper by modeling the positive and negative sequence impedances of Type III wind turbines for GFM operation mode. It is discovered that the GFM control naturally acts against the negative resistance of Type III wind turbines that results from an interaction between the proportional gain of the rotor-side converter current controller and negative slip. The developed sequence impedance models and improved damping behavior are verified using PSCAD simulations of a 2.5-MW Type III GFM wind turbine. The modeling predictions are also supported by experimentally measuring the sequence impedance response of a 2.5-MW Type III wind turbine during operation in GFL and GFM modes.

grid-following turbine↗

Optimizing Control for Efficient Load Shifting with Thermal Energy Storage in Existing HVAC Systems

This project developed the integration of a direct-contact heat exchanger (DCHX) based thermal energy storage (TES) system with a chiller–air handling unit (AHU) plant to evaluate its potential for reducing building energy costs. Detailed physical models of the TES unit, building envelope, and HVAC components were developed alongside simplified control-oriented models to support both high-fidelity simulation and real-time optimization. Two control strategies were implemented and compared: a rule-based control (RBC) aligned with utility time-of-use (ToU) rates, and a model predictive control (MPC) framework leveraging forecasts of building load, weather, and internal gains.Simulation results show that the RBC strategy reduced daily electricity costs by around 30% by shifting cooling production from on-peak to off-peak hours. In contrast, the MPC strategy achieved significantly greater performance, reducing daily operating costs by up to 44% and peak-hour costs by more than 60%. Both strategies maintained indoor thermal comfort within acceptable limits, with MPC further improving load distribution and reducing equipment cycling.The outcomes confirm that TES integration, particularly when coordinated with advanced predictive control, can provide substantial cost savings and on-peak demand reduction. . These findings directly support the U.S. Department of Energy’s goals for grid-interactive efficient buildings and demonstrate the potential of TES-enabled HVAC systems for scalable deployment across the commercial building.

99 GENERAL AND MISCELLANEOUS↗

A Learning-based Supervisory Control Architecture for Electric Vehicle Charging System Paired with Energy Storages

A Machine Learning-based predictive model is developed for optimal dispatching energy storage system integrated with Electric Vehicle battery charger. The model is effectively trained using transfer learning algorithm and successfully validated via measurement data to achieve high fidelity and accuracy. The study findings provide possibilities for future development of the presented approach to implement computational heavy predictive algorithms with multi-step predictions in emerging modular bidirectional charger topologies to realize vehicle-to-grid (V2G) technology. The online tuning can be explored using the new measurements obtained during the operation, henceforth, enhancing the controller performance for integrating multiple DERs into distribution power networks that shows huge interest for the utilities and regulatory agencies. Multiple operation scenatios have been evaluated and simulation results are discussed to verify intended performance of the proposed control solution.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Seminar 13 - Advanced Methods for Grid Integration of High-Performance Residential Communities

High-performance residential communities such as net-zero energy (NZE) and NZE-ready communities are emerging across the United States as many cities and states have committed to achieve ambitious energy efficiency and renewable energy goals between now and 2050. The achievement of these goals will transform the nation's communities into the holistic paradigms consisting of pervasive solar photovoltaics (PVs) and other renewable distributed energy resources, energy storage and flexible loads. However, integrating residential communities with high PV penetration may bring many challenges to the grid. For example, during many hours of the year, electric power is flowing out of PV-enabled homes and back into the grid. When this happens simultaneously in many nearby homes, distribution feeder voltage can be raised significantly. Intermittent solar power can also cause voltage flickers, leading to potential damages on both distribution feeder and home appliances. Existing solutions, such as PV curtailment, battery storage, and direct load control, are not able to address all the challenges in a cost-effective way. Better understanding is needed on how to deploy rooftop PV, how to predict building loads, and how to coordinate the homes in a community to enhance grid reliability. This seminar consists of four integral presentations that cover different aspects of a solution to these challenges. First, we will present an automated method to identify the best surfaces to deploy PV and size PV according to the PV penetration goals in a community. Second, we will present a streamlined way to create gray-box building models that are suitable for large-scale building-grid co-simulation. Next, we will present a data-driven method to estimate the end-use energy consumption to better coordinate the behind-the-meter resources to reduce utility cost, provide grid services, and improve renewable energy integration. Lastly, we will present a hierarchical control system that integrates the building models and the estimated non-dispatchable loads into the decision-making process to enhance the reliability and resilience of the residential community. Simulation results from an NZE-ready community in Fort Collins, Colorado will be presented in all four presentations and the guidelines on better grid integration of residential communities will be provided.

14 SOLAR ENERGY↗

Predicted Impacts of Pt and Ionomer Distributions on Low-Pt-Loaded PEMFC Performance

Low-cost, high performance proton exchange membrane fuel cells (PEMFCs) have been difficult to develop due to limited understanding of coupled processes in the cathode catalyst layer (CCL). Low-Pt-loaded PEMFCs suffer losses beyond those predicted solely due to reduced catalyst area. Although consensus links these losses to thin ionomer films in the CCL, a precise mechanistic explanation remains elusive. In this publication, we present a physically based PEMFC model with novel structure-property relationships for thin-film Nafion, validated against PEMFC data with low Pt loading. Results suggest that flooding exacerbates kinetic limitations in low-loaded PEMFCs, shifting the Faradaic current distribution. As current density increases, protons travel further into the CCL, resulting in higher Ohmic overpotentials. We also present a parametric study of CCL design parameters. We find that graded Pt and ionomer loadings reduce Ohmic losses and flooding, but individually do not provide significant improvements. However, a dual-graded CCL (i.e., graded Pt and ionomer) is predicted to significantly improve the maximum power density and limiting current compared to uniformly loaded CCLs. This work highlights the importance of accurate transport parameters for thin-film Nafion and provides a pathway to low-cost PEMFCs via precise control of CCL microstructures.

08 HYDROGEN↗