Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compositional simulation”

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 127 records · Page 7

Validation of Material Models For Automotive Carbon Fiber Composite Structures Via Physical And Crash Testing (VMM Composites Project)

As automotive fuel economy requirements increase, the push for reducing overall vehicle weight will likely include the consideration of materials that have not previously been part of mainstream vehicle design and manufacturing, including carbon fiber composites. Vehicle manufacturers currently rely on computer-aided engineering (CAE) methods as part of the design and development process, so going forward, the ability to accurately and predictably model carbon fiber composites will be necessary. If composites are to be used for structural components, this need applies to both, crash and quasi-static modeling. This final report covers the results of a five-year, $6.89M, 50% cost-shared research project between Department of Energy (DOE) and the US Advanced Materials Partnership (USAMP) under Cooperative Agreement DE-EE-0005661 known as “Validation of Material Models for Automotive Carbon Fiber Composite Structures Via Physical and Crash Testing (VMM).” The objective of the VMM Composites Project was to validate and assess the ability of physics-based material models to predict crash performance of automotive primary load-carrying carbon fiber composite structures. Simulation material models that were evaluated included micro-mechanics based meso-scale models developed by the University of Michigan (UM) and micro-plane models by Northwestern University (NWU) under previous collaborations with the DOE and Automotive Composites Consortium/USAMP, as well as five commercial crash codes: LS-DYNA, RADIOSS, VPS/PAM-CRASH, Abaqus, and GENOA-MCQ. CAE predictions obtained from seven organizations were compared with experimental results from quasi-static testing and dynamic crash testing of a thermoset carbon fiber composite front-bumper and crush-can (FBCC) system gathered under multiple loading conditions. This FBCC design was developed to demonstrate progressive crush, virtual simulation, tooling, fabrication, assembly, non-destructive evaluation and crash testing advances in order to assess the correlation of the predicted results to the physical tests. The FBCC was developed to meet a goal of 30-35% mass reduction while aiming for equivalent energy absorption as a steel component for which baseline experimental results were obtained from testing in the same crash modes. The project also evaluated crash performance of thermoplastic composite structures fabricated from commercial prepreg materials and low cost carbon fiber sourced from Oak Ridge National Laboratory. The VMM Project determined that no set of predictions from a CAE supplier were found to be universally accurate among all the six crash modes evaluated. In general, crash modes that were most dependent on the properties of the prepreg were more accurate than those that were dependent on the behavior of the joints. The project found that current CAE modeling methods or best practices for carbon fiber composites have not achieved standardization, and accuracy of CAE is highly reliant on the experience of its users. Coupon tests alone are not sufficient to develop an accurate material model, but it is necessary to bridge the gap between the coupon data and performance of the actual structure with a series of subcomponent level tests. Much of the unreliability of the predictions can be attributed to shortcomings in our ability to mathematically link the effects of manufacturing and material variability into the material models. This is a subject of ongoing research in the industry. The final report is organized by key technical tasks to describe how the validation project developed, modeled and compared crash data obtained on the composite FBCC to the multiple sets of CAE predictions. Highlights of the report include a discussion of the quantitative comparison between predictions and experimental data, as well as an in-depth discussion of remaining technological gaps that exist in the industry, which are intended to spur innovations and improvements in CAE technology.

36 MATERIALS SCIENCE↗

Material and component developments for the DEMO divertor using fibre reinforcement and additive manufacturing

Within the research along the European Fusion Roadmap, water-cooled divertor PFCs are foreseen in the design of a first fusion demonstration power plant (DEMO) in order to provide reliable heat removal capability. In the frame of this concerted attempt, the Max Planck Institute for Plasma Physics is concentrating on the development and testing of composite materials based on tungsten (W, preferred armour material) and copper (Cu, preferred heat sink material). W fibres (W f ) as monofilaments and yarns as reinforcement play a central role in these investigations due to their extraordinary properties concerning ductility already at room temperature and high tensile strength. Recent investigations on the impact of radiation damage suggest that the fibres retain their ductility upon irradiation. W reinforced with W fibres (W f /W) allows to overcome the intrinsic brittleness of W. Quantitative mechanical fracture tests of W f /W confirm the basic mechanisms of fibre reinforcement and the increased resistance to mechanical fatigue. The good wettability of W with liquid Cu and the absence of any metallurgical solubility make up an ideal material pairing for composite production. W fibre-reinforced Cu (W f /Cu) cooling tubes provide a rather high thermal conductivity (> 250 W mK -1 ) and at least twice the strength of CuCrZr in hoop direction in the temperature range up to at least 500 °C. Very recent neutron irradiation experiments confirm the sustainment of ductility of the W f /Cu composite. Numerical simulations suggest that thermal stresses in W-Cu PFCs could be strongly reduced by tailoring the local W and Cu volume fraction. This 'freely' distributed material composition can be achieved by means of additively manufactured W skeletons consecutively infiltrated by Cu. Investigations with W preforms produced by Laser Beam Powder Bed Fusion and infiltrated by Cu demonstrate the feasibility of this approach while testing of specifically prepared specimen is ongoing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Integrated Approach to Predicting Ash Deposition and Heat Transfer in Coal-Fired Boilers

The overall goal of this project is to develop via measurements and simulations an advanced online technology to predict, monitor and manage fireside ash deposition in a coal-fired boiler allowing for more efficient operations under a range of load conditions and fuel property variability. With this in place fuel sorting and blending can be done upstream and operations can be optimized to compensate for load and fuel properties. In support of this objective, three experimental campaigns were undertaken during the course of the project to measure ash deposition rates within the boiler at different fuel flow rates and its ash composition variability. Simulations of the experimental conditions representing actual geometry, operational scenarios in terms of air flow rates, coal flow rates as well as coal compositions, heating values, and particle size distributions were also carried out. Deposition rates were predicted using a unique particle kinetic energy and viscosity based ash deposition methodology whose validity was ascertained by comparing against deposition rate measurements for widely varying operating conditions and ash compositions in a lab-scale furnace. With a unique end-to-end combustion modeling methodology established and different simulation scenarios carried out, the results from our computational fluid dynamic (CFD) simulations in conjunction with the plant data summarized in this report were used to refine Microbeam Technology Incorporated’s MTI CSPI-CT Tool to predict and monitor fire-side ash deposition under a range of load conditions and fuel property variability in real time.

01 COAL, LIGNITE, AND PEAT↗

Delamination imaging in composites using cross-correlation method by non-contact air-coupled Lamb waves

Abstract This paper presents a cross-correlation damage detection technique using damage scattered waves of all directions cross-correlate with the incident waves for delamination imaging in carbon fiber reinforced polymer (CFRP) composite plate. The presented imaging can not only detect the delamination but also precisely quantify the damage with locations and sizes. First, a CFRP composite sample with interlaminar delamination by inserting the Teflon layer is designed and manufactured in house using the hot-press machine. And a three-dimensional model of composite sample is simulated using the finite element method. Next, the cross-correlation imaging algorithm is introduced in detail, and the cross-correlation imaging proof-of-concept study is conducted with the simulated Lamb waves in the composite sample. Finally, a fully non-contact air-coupled transducer and scanning laser Doppler vibrometer system with a single-mode Lamb wave method is established to actuate and sense the interrogating Lamb waves in the structure. The imaging method is experimentally implemented for the one delamination and two delaminations imaging and quantification.

Instruments & Instrumentation↗

High-Temperature Interactions Between Vanadium-Titanium Magnetite Carbon Composite Hot Briquettes and Pellets Under Simulated Blast Furnace Conditions

The high-temperature interactions between vanadium-titanium magnetite carbon composite hot briquettes (VTM-CCBs) and pellets were systematically investigated under simulated blast furnace conditions with respect to the reduction behavior, softening–melting–dripping characteristics, gas permeability, and Ti(C, N) precipitation mechanisms. The results showed that VTM-CCB charging can promote the reduction of the pellet in the packed bed and decrease the compressive strength of the pellet after reduction. The compressive strength of the VTM-CCB after reduction decreased with increasing temperature when the FC/O ratio (the ratio of the fixed carbon mol(C) in coal to the reducible oxygen mol(O) in iron oxides) was higher than 1.0. With an FC/O ratio lower than 1.0, the compressive strength of the VTM-CCB initially decreased and then increased. The FC/O ratio has a significant influence on the softening–melting interaction mechanism between the VTM-CCB and the pellet. With an FC/O ratio of 0.8, the bonding layer at the interface between the pellet and the VTM-CCB (consisting of molten fayalitic slag) can promote the softening process, thereby decreasing the softening start and end temperatures. By increasing the FC/O ratio to 1.4, a dense metallic iron shell with relatively high strength formed at the interface and restricted the collapse of the packed bed, thereby increasing the softening start and end temperatures and ensuring the transport of the reduction gas through the packed bed. The melting point of the primary slag phase increased with increasing FC/O ratio due to a decrease in the FeO content, which resulted in an increase in the melting start temperature from 1273 °C to 1294 °C (1546 K to 1567 K). The gas permeability in the cohesive zone increased with an increasing FC/O ratio of the VTM-CCB due to a combination of the skeletal role performed by the residual VTM-CCB and the decrease in the liquid slag proportion. In addition, as the FC/O ratio increased to 1.4, unconsumed carbon promoted the precipitation of Ti(C, N) at the slag–carbon and slag–metal interfaces, which resulted in a substantial increase in the dripping temperature and deterioration of the dripping behavior of the packed bed. Therefore, to suppress the precipitation of Ti(C, N) and improve the dripping behavior of the packed bed, the FC/O ratio of the charged VTM-CCB should be controlled within an appropriate range.

36 MATERIALS SCIENCE↗

Hybrid physics-based and machine learning tools for materials assessment

In this work, we develop novel physics-based and machine learning computational techniques to predict fundamental properties of metallic systems that affect radiation damage behavior in reactor structural materials. Oftentimes, atomistic predictions of engineering alloys simplify chemical compositions to a single element to reduce computational cost and complexity, introducing large sources of uncertainty and potentially missing important behavior. In addition, engineering alloys such as SS 316 are particularly challenging to simulate with first principles methods because of the additional degrees of freedom introduced by magnetic natures of the constituent elements, and very little data of this type exists within the literature. These novel methods aim to improve the qualitative prediction of radiation damage in engineering alloys by more accurately simulating their compositions, both by accelerating the computations and by developing novel analyses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

CGILS Phase 2 LES intercomparison of response of subtropical marine low cloud regimes to CO 2 quadrupling and a CMIP3 composite forcing change: Large eddy simulation of cloud feedbacks

We extended Phase 1 of the CGILS large-eddy simulation (LES) intercomparison in order to understand if subtropical marine boundary-layer clouds respond to idealized climate perturbations consistently in six LES models. Here the responses to quadrupled carbon dioxide (“fast adjustment”) and to a composite climate perturbation representative of CMIP3 multimodel mean 2×CO 2 near-equilibrium conditions are analyzed. As in Phase 1, the LES is run to equilibrium using specified steady summertime forcings representative of three locations in the Northeast Pacific Ocean in shallow well-mixed stratocumulus, decoupled stratocumulus, and shallow cumulus cloud regimes. Our results are generally consistent with a single-LES study of Bretherton et al. (2013) on which this intercomparison was based. Both quadrupled CO 2 and the composite climate perturbation result in less cloud and a shallower boundary layer for all models in well-mixed stratocumulus and for all but a single LES in decoupled stratocumulus and shallow cumulus, corroborating similar findings from global climate models (GCMs). For both perturbations, the amount of cloud reduction varies across the models, but there is less intermodel scatter than in GCMs. Furthermore, the cloud radiative effect changes are much larger in the stratocumulus-capped regimes than in the shallow cumulus regime, for which precipitation buffering may damp the cloud response. In the decoupled stratocumulus and cumulus regimes, both the CO 2 increase and CMIP3 perturbations reduce boundary-layer decoupling, due to the shallowing of inversion height.

54 ENVIRONMENTAL SCIENCES↗

Multimaterial Aerosol Jet Printed Magnetic Nanocomposites for Microwave Circuits

Here this article demonstrates the process and characterization of additively manufactured polymer-matrix magnetic nanocomposites. These nanocomposites are composed of nickel-zinc ferrite (Ni 0.5 Zn 0.5 Fe 2 O 4 ) nanoparticles in a polyimide matrix. These materials are printed using multimaterial aerosol jet printing where aerosols of nickel-zinc ferrite and polyimide inks are mixed in place to fabricate structures with nonunity relative permeability. To characterize these materials, we use energy-dispersive X-ray spectroscopy to estimate the content of ferrite in the composite, fitted simulation models, and ring resonator measurements to estimate the permeability. In fully aerosol jet printed 1-mm 2 inductors, we demonstrate an increase in inductance density of 40% in a 4.5-turn microstrip spiral inductor using this printed composite. These fully printed inductors demonstrate an inductance density of 4.2-nH/mm 2 . EDS analysis shows films with 18 wt.% Ni 0.5 Zn 0.5 Fe 2 O 4 . Finally, we show microstrip transmission lines on printed composite compared to the same on polyimide subjected to a dc magnetic bias, displaying a magnetic resonance behavior.

42 ENGINEERING↗

Insulation Activation and Contamination in Low-Power Experiments

Among the diverse advanced nuclear reactor concepts being developed are designs that expect to operate at much higher temperatures than conventional light-water reactors, and therefore require alternative materials for components capable of withstanding these extreme temperatures. Insulation is one such component where alternative high-temperature industrial composites are being considered. Low-power experiments can provide an opportunity for advancing our understanding of material behavior under irradiation; however, these experiments are not without occupational hazards. The focus of this paper is the potential for activation of insulating composites when exposed to neutron fluence during reactor experiments, and to discuss the subsequent contamination of the testing area. This paper will also discuss practical measures for preventing the creation and inhalation of activated particulate matter (i.e., dust), such as engineering and administrative controls and personal protective equipment. This paper evaluates the suitability of several insulation materials for use in reactor experiments based on their activation from neutron irradiation. The materials considered are Pyrogel XTE by Aspen Aerogel, Cerablanket and Kaowool by Morgan Advanced Materials, Thermo-12 Gold by Johns Manville, and Maxsil CF6-2000 by McAllister Mills. Although these materials are not yet widely used in the nuclear industry, they are common in other industries due to their thermal insulating properties. These materials produce dust while handling that can pose health and safety hazards, which could be further exacerbated by activation. The analysis presented in this paper utilizes a deterministic method developed to predict the activation source term of the provided insulation materials. The method calculates the atomic density of each element within the composites and simulates their exposure to a specified neutron flux density representative of low power experimental conditions. During the irradiation phase and throughout the extensive cooling period, neutron absorptions and nuclear decays are tracked for relevant isotopes. The residual activity within each composite is then assessed to identify the material with the lowest potential radiological hazard to workers. Dust will be created during installation and removal of insulating materials and may also be generated during the experiment. Contamination of surrounding surfaces due to dispersion is highly likely unless preventative measures are taken. The results of this analysis are used to estimate an internal dose from inhalation of activated dust, which is significantly lower than the allowable whole-body dose. However, the risk of internal exposure may increase if the dust becomes a transport vector for nuclear material. The goals of this paper are to enhance the understanding of the potential activation and contamination hazards from the use of these insulating materials in low-power nuclear experiments and to foster safer work environments through informed decision-making and strategic planning. This work ultimately contributes to a broader understanding of how new non-nuclear materials for advanced reactors impact the radiological safety of co-located workers, helping to ensure the responsible management of radioactive materials during installation and decommissioning activities.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

Evaluation of WRF simulation of deep convection in the US Southern Great Plains

The Southern Great Plains (SGP) exhibits a relatively high frequency of periods with extremely high rainfall rates (RR) and hail. Here, seven months of 2017 are simulated using the Weather Research and Forecasting (WRF) model applied at convection permitting resolution with the Mibrandt-Yau microphysics scheme. Simulation fidelity is evaluated, particularly during intense convective events, using data from ASOS stations, dual-polarization RADAR, gridded data sets and observations at the DoE Atmospheric Radiation Measurement site. The spatial gradients and temporal variability of precipitation and the cumulative density functions for both RR and wind speeds exhibit fidelity. Odds ratios >1 indicate WRF is also skillful in simulating high composite reflectivity (cREF, used as a measure of widespread convection) and RR > 5 mmhr –1 over the domain. Detailed analyses of the ten days with highest spatial coverage of cREF >30 dBZ show spatially similar reflectivity fields and high RR in both RADAR data and WRF simulations. However, during periods of high reflectivity, WRF exhibits a positive bias in terms of very high RR (> 25 mmhr –1 ) and hail occurrence, and during the summer and transition months, maximum hail size is underestimated. For some renewable energy applications fidelity is required with respect to the joint probabilities of wind speed and RR and/or hail. While partial fidelity is achieved for the marginal probabilities, performance during events of critical importance to these energy applications is currently not sufficient. Further research into optimal WRF configurations in support of potential damage quantification for these applications is warranted.

54 ENVIRONMENTAL SCIENCES↗

Nepheline structural and chemical dependence on melt composition

Nepheline crystallizes upon slow-cooling in some melts concentrated in Na 2 O and Al 2 O 3 , which can result in a residual glass phase of low chemical durability. Nepheline can incorporate many components often found in high-level waste radioactive borosilicate glass, including glass network ions (e.g., Si, Al, Fe), alkali metals (e.g., Cs, K, Na, and possibly Li), alkaline-earth metals (e.g., Ba, Sr, Ca, Mg), and transition metals (e.g., Mn, and possibly Cr, Zn, Ni). When crystallized from melts of different compositions, nepheline chemistry varies as a function of starting glass composition. Five simulated high level nuclear waste borosilicate glasses shown to crystallize large fractions of nepheline on slow cooling, were selected for study. These melts constituted a range of Al 2 O 3 , B 2 O 3 , CaO, Na 2 O, K 2 O, Fe 2 O 3 , and SiO 2 compositions. Compositional analyses of nepheline crystals in glass by electron probe micro-analysis (EPMA) indicate that boron is unlikely to be present in any significant concentration, if at all, in nepheline. Here, several models are also presented for calculating the fraction of vacancies in the nepheline structure.

36 MATERIALS SCIENCE↗

Review of strongly-coupled composite dark matter models and lattice simulations

We review models of new physics in which dark matter arises as a composite bound state from a confining strongly-coupled non-Abelian gauge theory. We discuss several qualitatively distinct classes of composite candidates, including dark mesons, dark baryons, and dark glueballs. Here, we highlight some of the promising strategies for direct detection, especially through dark moments, using the symmetries and properties of the composite description to identify the operators that dominate the interactions of dark matter with matter, as well as dark matter self-interactions. We briefly discuss the implications of these theories at colliders, especially the (potentially novel) phenomenology of dark mesons in various regimes of the models. Throughout the review, we highlight the use of lattice calculations in the study of these strongly-coupled theories, to obtain precise quantitative predictions and new insights into the dynamics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The origin of passivity in aluminum-manganese solid solutions

In this work, the effects of manganese on the aqueous corrosion of aluminum-manganese alloys were investigated by experiments and atomistic simulations. Electrochemical measurements, x-ray photoelectron spectroscopy, and atom probe tomography analysis indicate that manganese addition enhanced the corrosion resistance of aluminum without participating in the surface oxidation. Furthermore, the selective dissolution of manganese was believed to increase the free volume at the metal/oxide interface to facilitate the formation of a denser, thinner oxide layer with closer to stoichiometry composition. Atomistic simulations showed that the oxide layer compactness increased, and defect density decreased with increasing free volume content in Al, resulting in enhanced barrier characteristics.

36 MATERIALS SCIENCE↗

Design of an additively manufactured functionally graded material of 316 stainless steel and Ti-6Al-4V with Ni-20Cr, Cr, and V intermediate compositions

This study presents a method for designing a computationally informed gradient pathway to fabricate a functionally graded material (FGM) with terminal alloys of 316 stainless steel (SS316) and Ti-6Al-4V via directed energy deposition additive manufacturing with powder feedstock. In this work, the grading is accomplished through the introduction of intermediate elements and alloys (Ni-20Cr, Cr, and V) to avoid the brittle Fe-Ti intermetallic phases that form in the direct liquid phase joining of Ti-alloys and stainless steels. Using a combination of equilibrium calculations and Scheil-Gulliver simulations, a compositional pathway was designed to avoid deleterious phases. FGM samples were fabricated and experimentally characterized to determine the viability of the pathway. A change in phases from fcc to bcc was predicted to occur within the Ni-20Cr/Cr gradient region, and this was validated through experimental characterization. No detrimental phases (intermetallic, Laves, or σ phases) formed along the gradient path, demonstrating a successful computationally-informed design and fabrication of an FGM from SS316 to Ti-6Al-4V.

36 MATERIALS SCIENCE↗

Electronic structure simulations in the cloud computing environment

The transformative impact of modern computational paradigms and technologies, such as high-performance computing, quantum computing, and cloud computing, has opened up profound new opportunities for scientific simulations. Scalable computational chemistry is one beneficiary of this technological progress. The main focus of this paper is on the performance of various quantum chemical formulations, ranging from low-order methods to high-accuracy approaches, implemented in different computational chemistry packages, such as NWChem, NWChemEx, SPEC, ExaChem, and FLOSIC codes on the Azure Quantum Element (AQE) Microsoft cloud services. We pay particular attention to the intricate workflows for performing composite chemistry simulations, associated data curation, and mechanisms for accuracy assessment, as defined by the enabling cloud Computational Chemistry as a Service (CCaaS). Our focus also extends to Arrows' automated workflow for high throughput simulations. Finally, we provide a perspective on the role of cloud computing in supporting the mission of leadership computational facilities (LCFs).

computational chemistry, electronic structure, Clo↗

Accelerated Modeling of Lithium Diffusion in Solid State Electrolytes using Artificial Neural Networks

Abstract Previous efforts to understand structure‐function relationships in high ionic conductivity materials for solid state batteries have predominantly relied on density functional theory (DFT‐) based ab initio molecular dynamics (MD). Such simulations, however, are computationally demanding and cannot be reasonably applied to large systems containing more than a hundred atoms. Here, an artificial neural network (ANN) is trained to accelerate the calculation of high accuracy atomic forces and energies used during such MD simulations. After carefully training a robust ANN for four and five element systems, nearly identical lithium ion diffusivities are obtained for Li 10 GeP 2 S 12 (LGPS) when benchmarking the ANN‐MD results with DFT‐MD. Applying the ANN‐MD approach, the effect of chlorine doping on the lithium diffusivity is calculated in an LGPS‐like structure and it is found that a dopant concentration of 1.3% maximizes ionic conductivity. The optimal concentration balances the competing consequences of effective atomic radii and dielectric constants on lithium diffusion and agrees with the experimental composition. Performing simulations at the resolution necessary to model experimentally relevant and optimal concentrations would be infeasible with traditional DFT‐MD. Systems that require a large number of simulated atoms can be studied more efficiently while maintaining high accuracy with the proposed ANN‐MD framework.

Rao, Karun K.↗

LOCAL SIMULATIONS OF INSTABILITIES DRIVEN BY COMPOSITION GRADIENTS IN THE ICM

The distribution of helium in the intracluster medium (ICM) permeating galaxy clusters is not well constrained due to the very high plasma temperature. Therefore, the plasma is often assumed to be homogeneous. A nonuniform helium distribution can, however, lead to biases when measuring key cluster parameters. This has motivated one-dimensional models that evolve the ICM composition assuming that the effects of magnetic fields can be parameterized or ignored. Such models for nonisothermal clusters show that helium can sediment in the cluster core, leading to a peak in concentration offset from the cluster center. The resulting profiles have recently been shown to be linearly unstable when the weakly collisional character of the magnetized plasma is considered. In this paper, we present a modified version of the MHD code Athena, which makes it possible to evolve a weakly collisional plasma subject to a gravitational field and stratified in both temperature and composition. We thoroughly test our implementation and confirm excellent agreement against several analytical results. In order to isolate the effects of composition, in this initial study we focus our attention on isothermal plasmas. We show that plasma instabilities, feeding off gradients in composition, can induce turbulent mixing and saturate by rearranging magnetic field lines and alleviating the composition gradient. Composition profiles that increase with radius lead to instabilities that saturate by driving the average magnetic field inclination to roughly 45°. We speculate that this effect may alleviate the core insulation observed in homogeneous settings, with potential consequences for the associated cooling flow problem.

79 ASTRONOMY AND ASTROPHYSICS↗

The atomistic origin of the extraordinary oxygen reduction activity of Pt 3 Ni 7 fuel cell catalysts

Recently Debe et al. reported that Pt 3 Ni 7 leads to extraordinary Oxygen Reduction Reaction (ORR) activity. However, several reports show that hardly any Ni remains in the layers of the catalysts close to the surface (“Pt-skin effect”). This paradox that Ni is essential to the high catalytic activity with the peak ORR activity at Pt 3 Ni 7 while little or no Ni remains close to the surface is explained here using large-scale first-principles-based simulations. We make the radical assumption that processing Pt–Ni catalysts under ORR conditions would leach out all Ni accessible to the solvent. To simulate this process we use the ReaxFF reactive force field, starting with random alloy particles ranging from 50% Ni to 90% Ni and containing up to ~300 000 atoms, deleting the Ni atoms, and equilibrating the resulting structures. We find that the Pt 3 Ni 7 case and a final particle radius around 7.5 nm lead to internal voids in communication with the exterior, doubling the external surface footprint, in fair agreement with experiment. Then we examine the surface character of these nanoporous systems and find that a prominent feature in the surface of the de-alloyed particles is a rhombic structure involving 4 surface atoms which is crystalline-like but under-coordinated. Using density-functional theory, we calculate the energy barriers of ORR steps on Pt nanoporous catalysts, focusing on the O ad -hydration reaction (O ad + H 2 O ad → OH ad + OH ad ) but including the barriers of O 2 dissociation (O 2ad → O ad + O ad ) and water formation (OH ad + H ad → H 2 O ad ). We find that the reaction barrier for the O ad -hydration rate-determining-step is reduced significantly on the de-alloyed surface sites compared to Pt(111). Moreover we find that these active sites are prevalent on the surface of particles de-alloyed from a Pt–Ni 30 : 70 initial composition. These simulations explain the peak in surface reactivity at Pt 3 Ni 7 , and provide a rational guide to use for further optimization of improved catalytic and nanoporous materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗