Search NASA⌕ Search

SEARCH · Search NASA

Results for “Experimental Results”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Evaluation of Horizon of Viability Optimization Engine for Sustained Power to Critical Infrastructure: Preprint

In the aftermath of increasingly frequent catastrophic events, a typical scenario is Critical Infrastructure (CI) units being supported by available backup sources with a weak power grid that can be intermittent or absent. Such a scenario is significantly challenging in the sense of reliable supply of power to CI units. In this article, an intelligent optimization scheme termed as Horizon of Viability (HoV) engine is developed to guarantee the viability of a sustained reliable supply of power to the CI units over a time-horizon. The proposed HoV engine generates a cost-optimal portfolio of the locally available generation sources and the loads over a time horizon using a mixed-integer convex programming problem. A Controller hardware-in-the-loop (CHIL) platform is developed to evaluate the control performance of the HoV engine. The experimental results corroborates the efficacy in maintaining the viability of the CI units after a grid interruption event. Further, the proposed HoV optimization scheme performs better compared to existing net-load management schemes in the literature.

disaster resiliency↗
DOE OSTI2020

Impacts of Experimentally Obtained Harmonic Spectrums of Residential Appliances on Distribution Feeder

Owing to the increased use of power electronic based appliances and energy-efficient home equipment such as modern lighting loads, the percentage of nonlinear loads has increased in the distribution system, which can impact system performance, loss, and stability. Thus a comprehensive knowledge of their power quality and harmonic analysis is essential to improve the load models, voltage stability assessment, determination of possible interactions at harmonic frequencies, protection planning, and the effect of system impedance. This paper focuses on harmonic load flow analysis for multiple residential load types to determine the expected impacts on an modeled distribution secondary. For this study, the household appliances include lighting, power electronic, resistive, and motor loads with a nominal supply voltage of 120 V single-phase or 240 V split phase at a nominal frequency of 60 Hz. The harmonic spectrums, as obtained from the experimental evaluation of real loads, are used to inform the harmonic load flow for a detailed secondary network, including the distribution service transformer extracted from a real distribution feeder's model. Additionally, the impact of voltage harmonics on a three-phase test load is presented with experimental results. The harmonic data of the transformer terminal voltage, as obtained from the former load flow study, are scaled to generate the source voltage for the three-phase configuration of the grid simulator to make the test setup very similar to a typical secondary design in the U.S.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗
DOE OSTI2023

Pressure effects on reactivity and extinction of n -dodecane diffusion cool flame

Low-temperature combustion at high pressure has gained increasing interest owing to the growing demand for clean and efficient propulsion and energy systems. In particular, cool flame structure, extinction limits, and fuel reactivities at high pressure are critical in affecting the engine performance at near-limit conditions. Here, this paper investigates the effects of pressure on the cool flame extinction limit, structure, radical index, reactivity, and oxygen concentration dependence by experiments, analysis, and modeling. A large n-alkane, namely n-dodecane (n C 12 H 26 ), is selected to study its diffusion cool flame dynamics and reactivity in a high-pressure counterflow burner up to 10 atm. The experimental results show that higher pressure increases the cool flame extinction strain rates, and that the pressure-weighted extinction strain rate (aP) is proportional to the square of pressure, P 2 . A scaling analysis explains the relationship between the dependence of flame structure, heat release rate, and pressure-weighted cool flame strain rate on pressure. Furthermore, radical indexes at different pressures are measured by isolating the thermal and transport effects from the chemical contribution to diffusion cool flame extinction. The radical index clearly shows that the low-temperature reactivity increases with pressure. In addition, due to the critical role of multiple oxygen addition reactions in low-temperature chemistry, the relationship between the cool flame extinction limit and the oxygen concentration is explored. It is found that the cool flame extinction limits are proportional to the n th power of the oxygen concentration, [O 2 ] n , and increasing pressure leads to stronger extinction limit dependence (larger n) on the oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions, while the scaling analysis explains well the experimental results.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Performance Evaluation of Vertical Federated Machine Learning Against Adversarial Threats on Wide-Area Control System: Preprint

Federated machine learning (FL) is gaining significant popularity to develop cybersecurity solutions in power grids because of its advanced capability to support decentralized data handing at local devices, its privacy preservation, and its low-bandwidth requirement. However, the evolving adversarial machine learning (AML) threats raise significant concerns for the cybersecurity of FL architectures. The FL-based split neural network (SplitNN) achieves high performance through the decentralized training of local neural network models while preserving data privacy across multiple entities. In this paper, we propose a methodology for evaluating the performance of a vertical FLbased anomaly detector against different types of AML attacks, including denial-of-service attacks, adversarial data injection attacks, and replay attacks on the trained local models deployed in the grid network. For a case study, we consider the modified IEEE 13-bus system, and we develop SplitNN-based binary and multiclass classification models to detect, locate, and identify different types of data integrity attacks on the volt-watt control with two pooling layers: maximum pooling and AvgPool. Our experimental results, computed through performance metrics, reveal that the severity of these AML attacks varies with the integrated pooling mechanism, the type of classification model, and the nature of the cyberattack. Further, the AML attacks negatively impacted the prediction time per sample for the pretrained SplitNN during the online testing.

adversarial threats↗
DOE OSTI2022

Low Regeneration Temperature Sorbents for Direct Air Capture of CO 2

Susteon Inc., in partnership with University of Wyoming and SoCalGas, successfully met all major technical objectives to (1) scale up the ionic liquid catalyst for amine-based sorbents for improved desorption and absorption kinetics, (2) evaluate the catalyzed amine-based sorbents for direct CO 2 capture process to determine CO 2 adsorption and desorption rates and energy requirements, and (3) based on the experimental results, develop a conceptual process design to perform a preliminary economic assessment to evaluate the potential for DAC process cost reduction using the catalyzed sorbents. Amine doped solid sorbents are effective for DAC applications and can be regenerated by heat or by a combination of heat, steam, and vacuum. The best sorbent composition identified was polyethyleneimine (PEI) on fumed silica with 200 ppm ionic liquid catalyst. This sorbent formulation was shown to have a CO 2 breakthrough capacity twice that of the non-catalyzed sorbent, in laboratory tests with air at 75% relative humidity (RH). The CO 2 adsorption rate was also 40% higher than that of the non-catalyzed sorbent. This type of sorbents has the attributes required for lowering the overall cost of DAC with high CO 2 capacity and high rate of adsorption. The combination of an industrially utilized amine-based sorbent with a highly active catalyst to form a new class of materials for DAC provides a technically viable pathway for reducing the cost of DAC to <$100/tonne of CO 2 . Laboratory measurements show that the silica/PEI (polyethyleneimine) sorbents with 100 ppm of ionic liquid catalyst have almost 100% higher CO 2 cyclic capacity and 40% higher adsorption rate. Generally, CO 2 desorption occurred at higher temperatures with a rate of desorption 10 times faster than adsorption (which occurred at ambient conditions). Therefore, adsorption rate is a much more important factor in the cost of DAC because it is directly linked to the CAPEX of the total system and the cycle time (i.e., sorbent productivity in ton/day of CO 2 captured per unit volume of the air contactor). An initial process design, coupled with techno-economic analysis, based on optimal experimental results and preliminary resulting from structured sorbent testing, showed a path to lower the DAC cost from the current cost of over $200/tonne CO 2 to less than $100/tonne with a scale-up, mature state of the technology, with projected material and process improvements. These results demonstrate the effectiveness of the catalyst in silica/PEI sorbents in enhancing sorbents’ CO 2 working capacity, in (a) increasing the rate of adsorption and desorption, and (b) in lowering the CAPEX and OPEX of the DAC system employing the ionic liquid catalyzed sorbents.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Oxidation and Associated Pore Structure Modification During Experimental Alteration of Granite

Weathering plays a crucial role in a number of environmental processes, and the microstructure and evolution of multi-scale pore space is a critically important factor in the weathering of geological formations. In igneous rocks the infiltration of meteoric water into initially relatively dry material can cause the rock to begin to disaggregate, increasing porosity and surface area, and allowing further disaggregation and weathering to occur. These processes, in turn, allow biota to colonize the rock, further enhancing the weathering rate. In some igneous rocks this may be driven by oxidation of primary minerals. One such mineral, biotite, has been repeatedly mentioned as a cause of cracking during oxidation. However, the scale-dependence of the processes by which this occurs are poorly understood. We cannot, therefore, accurately extrapolate laboratory reaction rates to the field in predictive numerical models. In order to better understand the effects of oxidation and test the hypothesis that fracture and disaggregation are initiated by swelling of oxidizing biotites, we reacted granite cores in a selenic acid-rich aqueous solution at 200°C for up to 438 days. Elevated temperatures and selenic acid were used to provide relatively fast reaction rates and highly oxidizing conditions in sealed reaction vessels. These experiments were analyzed using a combination of imaging, X-ray diffraction, Mössbauer spectroscopy, and small- and ultra-small angle neutron scattering to interrogate porosity and microfracture formation. The experimental results show little observable biotite swelling, but significant transport and growth of iron oxides and/or clays along grain boundaries throughout the sample. Significant increases in porosity were also observed at the sample rim, likely associated with feldspar alteration. Fractures and transport were observed throughout the core, suggesting that stresses due to crystallization pressures caused by the growing iron phases may be the initiating factors in granite weathering, possibly followed by biotite swelling after sufficient permeability is achieved.

Anovitz, Lawrence M.↗
DOE OSTI2022

Solid-State Transformer and Hybrid Transformer With Integrated Energy Storage in Active Distribution Grids: Technical and Economic Comparison, Dispatch, and Control

Solid-state transformer (SST) and hybrid transformer (HT) are promising alternatives to the line-frequency transformer (LFT) in smart grids. The SST features medium-frequency isolation, full controllability for voltage regulation, reactive power compensation, and the capability of battery energy storage system (BESS) integration with multiport configuration. The HT has a partially-rated converter for fractional controllability and can integrate a small BESS. Fast grid-edge voltage fluctuations from increased solar photovoltaic (PV) and electric vehicle (EV) penetration are difficult to manage for mechanical load tap changers. Hence, along with the trend towards more BESS in the grid, the controllability and the storage integration capability of the SST and HT are of strong interest. However, a review of literature shows existing SST and HT research is mostly at converter level, while system-level assessments are scarce. Assessing technical and economic impacts is critical to understanding the benefits and role of the SST and HT to guide future research, which is presented for the first time in this article. Experimental results from medium-voltage (MV) SST and MV HT prototypes are shown to confirm equipment-level feasibility, where the voltage controllability waveforms of a MV HT prototype are reported for the first time. Comparative simulations are performed on a modified IEEE 34-bus system. Here, a grid-model-less decentralized grid-edge voltage control method and a day-ahead BESS dispatch method are proposed for the SST and HT. The simulations show that the SST and HT with integrated storage can host more PV, achieve peak shaving, mitigate voltage fluctuation and reverse power flow, and support energy arbitrage for operational cost reduction, as compared to the LFT. Moreover, comprehensive analyses of net present value (NPV) and internal rate of return (IRR) are performed under different installed PV capacities, HT’s partial converter ratings, and BESS capacities. Sensitivities to future cost reductions of the PV and BESS are studied. Although the NPV and IRR are currently negative, 60% capital cost reduction or 150% revenue increase will make the SST and HT economically viable in the use case studied.

14 SOLAR ENERGY↗
DOE OSTI2021

SURFACE FINISHING AND ELECTROLESS NICKEL PLATING OF ADDITIVELY MANUFACTURED (AM) METAL COMPONENTS

This study investigates the application of electroless nickel deposition on additively manufactured stainless steel samples. Current additive manufacturing (AM) technologies produce metal components with a rough surface. Rough surfaces generally exhibit fatigue characteristics, increasing the probability of initiating a crack or fracture to the printed part. For this reason, the direct use of as-produced parts in a finished product cannot be actualized, which presents a challenge. Post-processing of the AM parts is therefore required to smoothen the surface. This study analyzes chempolish (CP) and electropolish (EP) surface finishing techniques for post-processing AM stainless steel components CP has a great advantage in creating uniform, smooth surfaces regardless of size or part geometry EP creates an extremely smooth surface, which reduces the surface roughness to the sub-micrometer level. In this study, we also investigate nickel deposition on EP, CP, and as-built AM components using electroless nickel solutions. Electroless nickel plating is a method of alloy treatment designed to increase manufactured component's hardness and surface resistance to the unrelenting environment. The electroless nickel plating process is more straightforward than its counterpart electroplating.. We use low-phosphorus (2-5% P), medium-phosphorus (6-9% P), and high-phosphorus (10-13% P). These Ni deposition experiments were optimized using the L9 Taguchi design of experiments (TDOE), which compromises the prosperous content in the solution, surface finish, plane of the geometry, and bath temperature. The pre-and post-processed surface of the AM parts was characterized by KEYENCE Digital MicroscopeVHX-7000 and Phenom XL Desktop SEM. The experimental results show that electroless nickel deposition produces uniform Ni coating on the additively manufactured components up to 20 μm per hour. Mechanical properties of as-built and Ni coated AM samples were analyzed by applying a standard 10 N scratch test. Nickel coated AM samples were up to two times scratch resistant compared to the as-built samples. This study suggests electroless nickel plating is a robust viable option for surface hardening and finishing AM components for various applications and operating conditions.

Keywords: additive manufacturing, fatigue, chempol↗
DOE OSTI2024

Hypersonic Wake Velocity Measurements Using Acetone Molecular Tagging Velocimetry

A non-intrusive laser diagnostic known as molecular tagging velocimetry was used to find quantitative off-body velocity measurements in the wake of a sphere in the Mach 7 Ludwieg Tube Wind Tunnel located at The University of Texas at San Antonio. Acetone gas seeded in the flow was excited using the 4th harmonic of a pulse-burst Nd:YAG laser. The experimental results were used to validate and compare to simulations using a continuous Galerkin flow solver. Both the experimental and simulated results agreed on negative velocities in the viscous shear layer. Near the centerline of the sphere, average velocities of -200 m/s were observed due to reverse flow in the recirculation region. Outside of the shear layer, velocities of 800 m/s were observed in both the experimental and simulated results. The overall average uncertainty for the strut-mounted and free-flight case was estimated to be ±27 m/s. The flow features of the wake were visualized using schlieren imaging, and the experimental results agreed well with the simulated results.

Shoemake, Lawton↗
DOE OSTI2026

Dopant adsorption as a function of bulk concentration at a near 42º (100) twist grain boundary in SrTiO3

The enrichment of grain boundaries with dopant atoms is of critical importance for the macroscopic physical properties of materials. In thermodynamic equilibrium the Gibbs adsorption isotherm relates grain boundary excess of dopant atoms, their chemical potential in the adjacent bulk, and the respective interface energy. This study has used bicrystals with a near 42º (100) twist grain boundary in SrTiO3 to demonstrate that different kinetic pathways of Fe dopant atom additions converge towards comparable grain boundary configurations. Despite differences in bulk chemical potentials remarkably similar grain boundary excess quantities were observed. The experimental results indicate the general experimental feasibility to establish quantitative relationships between variations of grain boundary energy, interfacial excess, and overall dopant concentration. Improved experimental counting statistics are needed to distinguish solute interface excess as a function of bulk chemical potential.

Hahn, William [University of California, Davis]↗
DOE OSTI2020

Investigation and validation of the dynamic response of an acoustically levitated particle using the lattice Boltzmann method

The stable levitation of an analyte sample in an acoustic levitator is a primary requirement for accurate x-ray characterization of its scientific structure. A rigid particle oscillates in an under-damped manner when introduced into the node of established standing acoustic waves. This investigation has employed the lattice Boltzmann method (LBM), a computational fluid dynamics technique, for the analysis of such rigid particle dynamics in acoustic levitation. The simulation uses the two dimensional and nine velocity (D2Q9) Bhatnagar-Gross-Krook formulation to levitate a rigid 1.6 mm diameter nylon ($\rho$ = 1150 kg/m 3 ) particle in the air at standard pressure and temperature conditions. The presented work is the first reported simulation of realistic acoustic levitator boundary conditions using the LBM. The simulation can capture the particle-fluid interactions that produce dynamic levitation at less than one-period timescale in the ultrasonic frequency regime. An experiment was conducted by levitating a 1.6mm-diameter nylon sphere to estimate the oscillations, and the oscillating frequency was found to be 50 Hz. The dynamic simulation results are consistent with experimental results for particle oscillations within the same order of magnitude, indicating that LBM formulation can be successfully used to study acoustic levitation to understand and mitigate particle jitter. The distortion of the acoustic field due to a levitating particle's presence was also analyzed to demonstrate how the presence of the particle can disrupt adjacent levitating nodes.

42 ENGINEERING↗
DOE OSTI2022

A method to predict texture effect on ion beam channeling analysis of polycrystals and the application to study the mosaic spreading effect in highly oriented pyrolytic graphite

We propose a method to convert the channeling Rutherford backscattering spectrum yield map of a single crystal to a polycrystal through a matrix rotation technique. The rotation matrix is determined by the deviation of the crystal axial direction from the original z axis. The final yield map is created after averaging the rotated yields using the texture function as the weight factor. For highly oriented pyrolytic graphite (HOPG) exhibiting mosaic spread, the method leads to a Gaussian kernel averaging of the map obtained from a single crystal. The yield map of a single crystal is obtained by a simulation of ion trajectories in a potential field described by Moliere screened Coulomb potentials. Yield maps are calculated under various σ values (standard deviations of mosaic spread). The simulated results are compared with experimental results obtained using 1.2 MeV alpha particle. σ is extracted through the best fitting, demonstrating that the method can be used to obtain texture details. The effects of mosaic spread on minimum yield χ min and the half-width at half maximum of angular scans ψ ½ are systematically modeled and compared with previous theoretical equations. The study also shows that previous theoretical equations are valid only at small σ values. The proposed method can be applied to any type of polycrystal and is not limited to HOPG. It provides near-surface mosaic spread and crystallography information with a longitudinal depth resolution of tens of nanometers and is not influenced by grain shapes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Ephemeral antibubbles: Spatiotemporal evolution from direct numerical simulations

Antibubbles, which consist of a shell of a low-density fluid inside a high-density fluid, have several promising applications. We show, via extensive direct numerical simulations (DNSs), in both two and three dimensions, that the spatiotemporal evolution of antibubbles can be described naturally by the coupled Cahn-Hilliard-Navier-Stokes (CHNS) equations for a binary fluid. Our DNSs capture elegantly the gravity-induced thinning and breakup of an antibubble via the time evolution of the Cahn-Hilliard scalar-order-parameter field $\phi$, which varies continuously across interfaces, so we do not have to enforce complicated boundary conditions at the moving antibubble interfaces. To ensure that our results are robust, we supplement our CHNS simulations with sharp-interface volume-of-fluid DNSs. We track the thickness of the antibubble and calculate the dependence of the lifetime of an antibubble on several parameters; we show that our DNS results agree with various experimental results; in particular, the velocity with which the arms of the antibubble retract after breakup scales as σ 1/2 , where σ is the surface tension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Traffic Shaping to Traffic Engineering in Time-Sensitive OT Network

Modern industrial automation systems increasingly depend on network infrastructures for time-critical communication, driving the need for solutions that guarantee timely and reliable data delivery. IEEE 802.1 Time-Sensitive Networking (TSN) holds significant promise for converging Information Technology (IT) and Operational Technology (OT) networks, enabling interoperability and supporting the coexistence of mixed-critical traffic crucial for Industry 4.0 and IIoT. To achieve deterministic communication, TSN employs various traffic shapers such as the Time-Aware Shaper (TAS), Asynchronous Traffic Shaper (ATS), and Credit-Based Shaper (CBS). However, the effective deployment of TSN in industrial automation faces several challenges. These include the non-trivial mapping of diverse industrial traffic types to specific shapers, the complexity of optimizing shaper configurations. We present a model for effective traffic engineering within TSN enabled OT Network. Our experiments also demonstrate how shaping of certain traffic types get affected in absence of precise time synchronization and propose possible solutions based on experiment results. Based on our experimental results we provide recommendations on how traffic type assignments should be done and which traffic shaping mechanisms should be used for a particular traffic type.

Sarker, Taposh Kumer [University of Texas at El Pa↗
DOE OSTI2025

A Modified DBC Substrate Improving Thermal Performance for Confined Space Applications

High-power modules use substrates to house the semiconductor device and for electrical insulation. These substrates are constructed with thermally conductive dielectric material sandwiched between two metals to extract heat from semiconductor chips. Thus, the required cooling performance of a power module is linked to the substrate’s thermal performance and can vary based on the substrate technologies. Here, in this study, five substrate technologies were evaluated for space-restricted applications: direct-bonded copper (DBC), an insulated metal substrate (IMS), a thermally annealed pyrolytic graphite (TPG)-based IMS, DBC-based double-sided cooling, and direct-bonded aluminum (DBA) where the heat sink is directly attached without thermal interface materials (TIMs). The finite element (FE) analysis results suggest that the popular DBC substrate has thermal performance better than that of the other substrates for space-constrained applications. To further improve the thermal performance, a modified DBC substrate was proposed where a copper block was added between the semiconductor and the DBC substrate to achieve heat spreading underneath the chip. The modified DBC performance was then compared with the aforementioned substrates and the results showed significant thermal performance improvement. The results were verified with experimental results where the proposed substrate showed 20% more loss handling capability compared to an identical DBC substrate.

42 ENGINEERING↗
DOE OSTI2024

New Advances in Optical Stochastic Cooling

Recently, Optical Stochastic Cooling (OSC) became the first demonstrated method for ultra-high-bandwidth stochastic cooling. The initial experiments at Fermilab’s IOTA ring explored the essential physics of the method and demonstrated cooling, heating and manipulation of beams and single particles. Having been validated in practice, with continued development, OSC carries the potential for dramatic advances in the state-of-the-art performance and flexibility for beam cooling and control. The ongoing program at Fermilab is now focused on the development of an OSC system that includes high-gain optical amplification, which promises a two-order-of-magnitude increase in the strength of the OSC force. In this talk, we briefly review the results of the initial experimental campaign, describe the status of the conceptual and hardware designs for the amplified OSC system, report initial experimental results of our high-gain amplifier development, and explore near-term operational plans and use cases.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

FEM Modeling and Simulation of 2-D High Specific Heat Nb 3 Sn Wires

In the past few years, new high specific heat Nb 3 Sn wires have gained much focus at FNAL. Indeed, they have proved to be more stable against thermal perturbations with respect to standard wires. Nevertheless, a trade off exists between their thermal efficiency and production feasibility. In this report I describe the thermal and structural models that I have developed by exploiting ANSYS Mechanical APDL ® , which I got acquainted with at the beginning of my training. The aim has been to optimize the location of high specific heat elements in order to obtain an optimal thermal stability, while minimizing the risk of wire breakage during drawing. FEM results are compared with experimental ones made beforehand. Other minimum quench energy (MQE) thermal models have been developed for expected new experimental results. New data may improve the understanding of the embedded physical uncertainties in the model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Additive Manufactured Composite Phase-Change Material for Thermal Energy Storage Applications

Phase-change materials play a critical role in industrial energy storage applications to drive efficiency improvements, thermal energy management, and carbon emissions reductions. Recently, it has been shown that rapid solidification of alloys with metastable immiscibility in the liquid phase has the potential to form unique microstructures in which a low-melting phase is uniformly distributed in a high-melting matrix. This feature can be exploited using additive manufacturing to produce components with complex geometries containing such unique phase-change microstructures. Phase-field simulations utilizing high-performance computing were used to provide a detailed description of the evolution of the active phase during service in terms of their morphology and composition in different polycrystalline matrix grain morphologies that are typically produced during additive manufacturing. Phase field simulations were performed using, MEUMAPPS-SL (Microstructure Evolution Using Massively Parallel Phase-field Simulations – Solid Liquid) code that was developed in-house by the Oak Ridge National Laboratory. The simulations utilized the capabilities of the Kestrel supercomputer at the National Renewable Energy Laboratory. The simulation results were compared with experimental results generated at Siemens Energy, Inc. The results indicate that the kinetics of liquid spreading along grain boundaries is largely determined by the mobility of the triple line along the intersection of the grain boundary liquid and the grain boundary plane.

25 ENERGY STORAGE↗