Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary Layer Stability”

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

Kinetics and Thermodynamics of Sr Permeation in CeO 2 -Based Barrier Layers for Solid-Oxide Electrolyzer Cells

Solid-oxide electrolyzer cells (SOECs) convert steam to hydrogen efficiently at high temperatures. However, during operation, the diffusion of cations or impurities through the cells due to electrode degradation can cause unwanted secondary phases to form, which may degrade device performance. Here, in this study, we use atomistic and mesoscale simulations coupled with experimental analysis to study the diffusion of Sr through the Gd-doped CeO 2 (GDC) barrier layer used to protect the yttria-stabilized zirconia (YSZ) electrolyte in SOECs. From our atomistic calculations, we find Sr diffusion to be negligibly slow in bulk GDC; however, surface diffusion is much more favorable. Subsequent mesoscale simulations show that Sr diffusion is activated when the porosity of GDC exceeds ∼10% and significantly exceeds diffusion in bulk and grain boundary regions. We also find that SrO-based species can accumulate at GDC surfaces; however, SrO aggregation and coarsening will be limited by the large lattice mismatch between GDC and SrO. Energy-dispersive X-ray spectroscopy (EDS) and electron diffraction confirm that Sr can accumulate within GDC pores and form disperse Sr-containing secondary phases. Altogether, Sr diffusion in dense GDC is unlikely to give rise to thick SrO layers, which would severely limit device performance. The formation of Sr-containing secondary phases can largely be avoided by restricting the porosity of the GDC layer as much as possible.

36 MATERIALS SCIENCE↗

Pit Stability Predictions of Additively Manufactured SS316 Surfaces Using Finite Element Analysis

Stainless steels are susceptible to localized forms of corrosion attack, such as pitting. The size and lifetime of a nucleated pit can vary, depending on a critical potential or current density criterion, which determines if the pit repassivates or continues growing. This work uses finite element method (FEM) modeling to compare the critical pit radii predicted by thermodynamic and kinetic repassivation criteria. Experimental electrochemical boundary conditions are used to capture the active pit kinetics. Geometric and environmental parameters, such as the pit shape and size (analogous to additively manufactured lack-of-fusion pores), solution concentration, and water layer thickness were considered to assess their impact on the pit repassivation criterion. The critical pit radius (the transition point from stable growth to repassivation) predicted for a hemispherical pit was larger when using the repassivation potential (E rp ) criteria, as opposed to the current density criteria (pit stability product). Including both the pit stability product and E rp into its calculations, the analytical maximum pit model predicted a critical radius two times more conservative than the FEA approach, under the conditions studied herein. The complex pits representing lack-of-fusion pores were shown to have minimal impact on the critical radius in atmospheric conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nanometric modulations of the magnetic structure of the element Nd

The rare earth neodymium arguably exhibits the most complex magnetic ordering and series of magnetic phase transitions of the elements. Here we report the results of small-angle neutron scattering (SANS) measurements as a function of temperature and applied magnetic field to study magnetic correlations on nanometer length scales in Nd. We support these results with additional neutron diffraction, thermodynamic and synchrotron x-ray measurements. The SANS measurements reveal the presence of modulation vectors characterizing the ordered spin configuration which exhibit changes in magnitude and direction that are phase dependent. Between 5.9 and 7.6 K an additional modulation vector is observed with magnitude Q = 0.12 Å -1 . This scattering appears to originate from order of the Nd layers which contain a center of inversion. Tracking this modulation vector as a function of magnetic field indicates a phase boundary at ≈1 T. The modulation vectors observed in this study establish the presence of nanometer length scale spin textures which are likely stabilized by frustrated Ruderman–Kittel–Kasuya–Yosida (RKKY) interactions.

36 MATERIALS SCIENCE↗

Aerosols consistently suppress the convective boundary layer development

Aerosols with different vertical distribution and various optical properties induce diverse heating rates and thereby affecting convective boundary layer (CBL) development. Our results showed consistent CBL-suppression of aerosols during daytime with numerical experiments, in which aerosols were specified at different heights with synthesized single scattering albedo from 64 studies and asymmetry factor from 20 studies globally. Absorbing aerosols concentrated below but close to the CBL top had the strongest suppression effect on CBL development relative to that concentrated near surface or above CBL. Aerosol cooling effect by attenuating incident solar radiation and surface heat flux exceeded its warming effect by reheating the atmosphere layer with absorbed shortwave radiation, and eventually declined net heating rate, which inhibited CBL development, lowered mixed-layer potential temperature and stabilized atmospheric stratification. Stove effect of absorbing aerosols (CBL enhancement) under a zero background aerosol extinction coefficient is negligible for dominant dome effect (CBL suppression) which consistently suppresses CBL development regardless of aerosol vertical height and background aerosol extinction coefficient. Our study also highlighted the importance of specifying background aerosol extinction coefficient in numerical experiments for accurate assessment of aerosol radiative forcing and CBL-aerosol interactions.

54 ENVIRONMENTAL SCIENCES↗

Self-consistent theory of the Darrieus–Landau and Rayleigh–Taylor instabilities with self-generated magnetic fields

The Rayleigh–Taylor (RT) and Darrieus–Landau (DL) instabilities are studied in an inertial confinement fusion context within the framework of small critical-to-shell density ratio D R and weak acceleration regime, i.e., large Froude number Fr. The quasi-isobaric analysis in Sanz et al. is completed with the inclusion of non-isobaric and self-generated magnetic-field effects. The analysis is restricted to perturbation wavelengths k~1 larger than the conduction length scale at the ablation front, yet its validity ranges from wavelengths shorter and larger than the conduction layer width (distance between the ablation front and the critical surface). The use of a sharp boundary model leads to a single analytical expression of the dispersion relation encompassing both instabilities. The two new effects come into play by modifying the perturbed mass and momentum fluxes at the ablation front. The momentum flux (perturbed pressure at the spike) is the predominant stabilizing mechanism in the RT instability (overpressure) and the driving mechanism in the DL instability (underpressure). The non-isobaric effects notably modify the scaling laws in the DL limit, leading to an underpressure scaling as ~k –11/15 rather than ~k –2/5 obtained in the quasi-isobaric model. The magnetic fields are generated due to misalignment between pressure and density gradients (Biermann battery effect). They affect the hydrodynamics by bending the heat flux lines. Within the framework of this paper, they enhance ablation, resulting in a stabilizing effect that peaks for perturbation wavelengths comparable to the conduction layer width. Here, the combination of parameters D R Fr 2/3 defines the region of predominance of each instability in the dispersion relation. It is proven that the DL region falls outside of the parameter range in inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterization of wind speed and directional shear at the AWAKEN field campaign site

The American wake experiment (AWAKEN) is taking place in northern Oklahoma, USA, close to the Atmospheric Radiation Measurement Southern Great Plains (ARM SGP) atmospheric observatory. The planning for the deployment of the instruments in this observational field campaign required an assessment of the wind characteristics of the site. This paper analyzes long-term data collected by instruments at the ARM SGP observatory to characterize the winds near the AWAKEN site. The analysis shows that this site experiences high wind shear and veer events with a large number of nocturnal low-level jets. A total of 7086 low-level jet wind profiles over 6 years are examined and found to be dominant from the south and southeast. Significant nocturnal wind veer is observed, which causes southerly wind near the surface to become westerly wind aloft. By identifying a strong relationship between atmospheric stability and wind shear, the wind shear at the site is predicted using the Monin–Obukhov similarity theory (MOST) and validated with the observational data collected by a scanning Doppler lidar. The results show that wind speed at a height of 91 m, a proxy hub height for wind turbines in this area, can be predicted from data collected at a height of 10 m with a bias of −0.35 and 0.65 m s−1 in unstable and stable atmospheric boundary layers, respectively. The bias of the predicted wind speed is mostly in the region of low wind speed, and wind speed above 5 m s−1 at a height of 91 m can be predicted with a bias of less than 0.2 m s−1, and the limitations of the MOST in predicting winds during the stably stratified boundary layer is well-observed.

17 WIND ENERGY↗

A Twenty-Year Analysis of Winds in California for Offshore Wind Energy Production Using WRF v4.1.2

Offshore wind resource characterization in the United States relies heavily on simulated winds from numerical weather prediction (NWP) models, given the lack of hub-height observations offshore. One such NWP data set used extensively by U.S. stakeholders is the Wind Integration National Dataset (WIND) Toolkit, a 7-year time-series data set produced in 2013 by the National Renewable Energy Laboratory. In this study, we present an update to that data set for offshore California that leverages recent advancements in NWP modeling capabilities and extends the period of record to a full 20 years. The data set predicts a significantly larger wind resource (0.25–1.75 m s-1 stronger), including in three Call Areas that the Bureau of Ocean Energy Management is considering for commercial activity. We conduct a set of yearlong simulations to study factors that contribute to this increase in the modeled wind resource. The largest impact arises from a change in the planetary boundary layer parameterization from the Yonsei University scheme to the Mellor-Yamada-Nakanishi-Niino scheme and their diverging wind profiles under stable stratification. Additionally, we conduct a refined wind resource assessment at the three Call Areas, characterizing distributions of wind speed, shear, veer, stability, frequency of wind droughts, and power production. We find that, depending on the attribute, the new data set can show substantial disagreement with the WIND Toolkit, thereby driving important changes in predicted power.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

Perovskite superlattices with efficient carrier dynamics

Compared with their three-dimensional (3D) counterparts, low-dimensional metal halide perovskites (2D and quasi-2D; B 2 A n-1 M n X 3n+1 , such as B = R-NH 3 + , A = HC(NH 2 ) 2 + , Cs + ; M = Pb 2+ , Sn 2+ ; X = Cl - , Br - , I - ) with periodic inorganic-organic structures have shown promising stability and hysteresis-free electrical performance. However, their unique multiple-quantum-well structure limits the device efficiencies because of the grain boundaries and randomly oriented quantum wells in polycrystals. In single crystals, the carrier transport through the thickness direction is hindered by the layered insulating organic spacers. Furthermore, the strong quantum confinement from the organic spacers limits the generation and transport of free carriers. Also, lead-free metal halide perovskites have been developed but their device performance is limited by their low crystallinity and structural instability 11 . Here, in this paper, we report a low-dimensional metal halide perovskite BA 2 MA n-1 Sn n I 3n+1 (BA, butylammonium; MA, methylammonium; n = 1, 3, 5) superlattice by chemical epitaxy. The inorganic slabs are aligned vertical to the substrate and interconnected in a criss-cross 2D network parallel to the substrate, leading to efficient carrier transport in three dimensions. A lattice-mismatched substrate compresses the organic spacers, which weakens the quantum confinement. The performance of a superlattice solar cell has been certified under the quasi-steady state, showing a stable 12.36% photoelectric conversion efficiency. Moreover, an intraband exciton relaxation process may have yielded an unusually high open-circuit voltage (V OC ).

36 MATERIALS SCIENCE↗

In situ energetics modulation enables high-efficiency and stable inverted perovskite solar cells

In contrast to conventional (n-i-p) perovskite solar cells (PSCs), inverted (p-i-n) PSCs offer enhanced stability and integrability with tandem solar cell architectures, which have garnered increasing interest. However, p-i-n cells suffer from energy level misalignment with transport layers, imbalanced transport of photo-generated electrons and holes, and significant defects with the perovskite films. Here we introduce tris(2,4,6-trimethyl-3-(pyridin-3-yl)phenyl)borane (3TPYMB), a nonionic n-type molecule that, through hydrogen bonding and Lewis acid-base reactions with perovskite surfaces or grain boundaries, enables in situ modulation of perovskite energetics, effectively mitigating the key challenges of p-i-n PSCs. The p-i-n PSCs incorporating 3TPYMB achieve a certified quasi-steady-state power conversion efficiency of 24.55 +- 0.33%, with a reverse scan efficiency of 25.58%. Finally, they also exhibit exceptional stability, with unencapsulated devices retaining 97.8% of their initial efficiency after 1,800 h of continuous operation at maximum power point under N 2 atmosphere, 1 sun illumination and 60 °C conditions.

14 SOLAR ENERGY↗

Characterization of oceanic post-cold frontal clouds and their model representation. Final report

Low-level clouds are ubiquitous over the oceans, and their impact on the Earth’ radiative budget is important. However, general circulation models (GCMs) that are used to represent our climate and its evolution still experience problems to represent their correct amount and radiative impact. In the extratropical latitudes (~30-60 N/S), cloud properties are strongly modulated by the occurrence of anti-cyclones (high pressure systems) and cyclones (low pressure systems). While deep clouds dominate in cyclones near the center and along the cold and warm fronts, low-level clouds populate the colder regions of the cyclones, and more specifically the region in the wake of the cold fronts, the post-cold frontal region (PCF). The GCMs mentioned above were found to underestimate the amount of these PCF clouds, and this caused errors in their representation of the amount of solar radiation at the surface of the southern hemisphere oceans, causing an excess in absorption and errors in long-term predictions. To better understand which processes are ill-represented in the GCMs to cause such issues, we used ground-based observations from the Eastern North Atlantic (ENA) ARM site and a high resolution regional model (called WRF) to examine the properties of PCF clouds. With the WRF model, we simulated the passage of a cold front at the ENA site using theoretically distinct representations of convection and the physics of the boundary layer. While these physical representations had little impact on the timing, structure and circulation of the cold front, we found a high sensitivity to how the large scale information was fed into the model. Errors rapidly developed if the incoming flow information was injected too far from the site, but 1000 km was found to be an optimal distance between the outer boundary and the ENA site for a realistic cold front passage. With this model, we then conducted a series of experiment where convection and boundary layer schemes were changed, and revealed that 1) the clouds were more sensitive to the convection than boundary layer schemes and 2) the interaction between the two has a large impact. Furthermore, changing the representation of the physics in the model had greater impacts on PCF clouds than changing the initial conditions. As large-scale models increase their spatial resolution, convection will eventually be resolved but it remains that their representation of the physics of the boundary layer will matter for the correct representation of shallow (low-level) clouds. Observations were used to find which large-scale information had a clear relationship with PCF cloud and precipitation properties, in view of helping with GCM parameterizations and evaluation. Both lower troposphere stability and surface forcing were found to be important for PCF clouds that experience more vigorous dynamics than their tropical or subtropical counterparts. For cloud base and top heights, a simple measure of the change in potential temperature between the surface and the 800 hPa was found to be a good predictor. This relationship was found to hold as well when using data from campaigns recently conducted in the southern ocean south of Tasmania. A recent GCM was found to succeed in representing similar relations. For cloud depth, water content, and different characteristics of precipitation both the temperature contrast between the surface and near-air and the near-surface wind speed had a strong relationship. To summarize, through observational constraints and regional model experiments, the project has provided guidance to evaluate GCMs and tools to improve their representation of low-level clouds in the extratropics. It has also highlighted differences in behavior between these clouds and those typically found in the subtropics, thus opening new avenues of research for shallow clouds in general.

54 ENVIRONMENTAL SCIENCES↗

Enhancing the Reaction Kinetics and Stability of Co‐Free Li‐Rich Cathode Materials via a Multifunctional Strategy

Co-free Li-rich layered oxides (CFLLOs) with anionic redox activity are among the most promising cathode materials for high-energy-density and low-cost lithium-ion batteries (LIBs). However, irreversible oxygen release often causes severe structural deterioration, electrolyte decomposition, and the formation of unstable cathode-electrolyte interface (CEI) film with high impedance. Additionally, the elimination of cobalt elements further deteriorates the reaction kinetics, leading to reduced capacity and poor rate performance. Here, in this study, a multifunctional strategy is proposed, incorporating Li 2 MnO 3 phase content regulation, micro-nano structure design, and heteroatom substitution. The increased content of Li 2 MnO 3 phase enhances the capacity through oxygen redox. The smaller nanoscale primary particles induce greater tensile strain and introduce more grain boundaries, thereby improving the reaction kinetics and reactivity, while the larger micron-sized secondary particles help to reduce interfacial side reactions. Furthermore, Na⁺ doping modulates the local coordination environment of oxygen, stabilizing both the anion framework and the crystal structure. As a result, the designed cathode exhibits enhanced rate performance, delivering a capacity of 158 mAh g⁻¹ at 5.0 C and improved cyclic stability, with a high capacity retention of 99% after 400 cycles at 1.0 C. This multifunctional strategy holds great promise for advancing the practical application of CFLLOs in next-generation LIBs.

Co-free Li-rich layered oxide↗

First-principles investigation of phase stability in layered Na x ⁢CrO 2

Layered oxide intercalation compounds continue to attract interest as electrode materials for Na-ion batteries. However, many of these materials undergo complex phase transitions during cycling that influence battery performance but are still not completely understood. We have conducted a first-principles study of layered Na x ⁢CrO 2 (0 ≤ x ≤ 1) to assess phase stability between various Na-vacancy ordered phases in the O3 and P3 host structures. We predict that many of the low-energy phases belong to families of Na orderings that follow specific patterns. At high x, we identify families of vacancy row orderings in O3, which may also couple to magnetic orderings of the Cr spins. We predict similar orderings at intermediate x in P3 that contain antiphase boundaries between regions of the x = 1/2 ordering. In both cases, the average spacing between rows/boundaries is set by the overall composition. At x = 0, we find a strong preference for charge disproportionation and migration of Cr to tetrahedral sites in the intercalation layers. Finally, we obtain generally good agreement with experimental observations and rationalize key discrepancies.

25 ENERGY STORAGE↗

Processing of SOFC Anodes for Enhanced Intermediate Temperature Catalytic Activity at High Fuel Utilization (Final Report)

The overall objective was to infiltrate anodes with nanoparticle catalysts so that anode catalytic performance can be improved as anode operation temperature is reduced from 800°C to intermediate temperatures of 700°C and 600°C and as anode fuel utilization is increased. The addition of nickel nanoparticle catalysts into the anode improves performance by providing additional triple phase boundaries (TPBs), increasing the density of electrochemical reaction sites. To achieve this goal, both liquid phase and vapor phase methods for depositing nickel nanoparticles in the anode active layer, the region of the anode near the electrolyte where the electrochemical reaction occurs. It is imperative that nanoparticles are stable against coarsening during long-term operation, so the mechanisms and kinetics of nanoparticle stability during cell operation will be investigated. Finally, priority is given to processes that are scalable and easily transferable to industry, leveraging existing Ni-YSZ cermet anode technology. Utilizing nickel nanoparticles at different temperatures has different considerations that will be addressed. At the upper end of the intermediate temperature range (800°C), it is important that performance improvement is maintained over long periods by mitigating the instability of nanoparticles. At the middle of the temperature range (700°C), it is important that performance improvement is maintained at high fuel utilizations. At the lower end of the temperature range (600°C), where the performance of cells is poor due to increased resistances, it is important to improve anode performance by significantly increasing the electrochemical reaction site density, while nanoparticle durability is less of an issue. Finally, the objective was to explore mixed ionic and electronic (MIEC) nanoparticle catalysts like GDC to see if 2-phase boundaries can be created to increase the electrochemical reactions at the anode.

20 FOSSIL-FUELED POWER PLANTS↗

Improving Thermal Stability of Perovskite Solar Cells by Suppressing Ion Migration Using Copolymer Grain Encapsulation

Thermal stability of organic–inorganic hybrid perovskites (OIHPs) remains as one of the critical challenges against the stable operation of perovskite solar cells (PSCs) in direct sunlight with elevated temperatures. Here, we show that the addition of a polystyrene-co-polyacrylonitrile (SAN) copolymer can significantly enhance thermal stability of OIHPs and improve the stability of the corresponding PSCs by suppressing the migration of organic cations in OHIP. The methylammonium lead iodide (MAPI) with SAN incorporated within the perovskite layer featured a superior thermal stability compared to pure MAPI without SAN, only displaying an average of 5–15% decrease in PCE even after continuous thermal aging for 24 h at 100 °C. The secondary ion mass spectrometry revealed that the thermal degradation of the pure MAPI was largely associated with MA+ out-migration. Conducting atomic force microscopy analysis further indicated that the incorporated SAN led to a suppression of ionic currents present at the grain boundaries of the perovskite film, which was understood by high immiscibility between SAN and MA+ components as confirmed by the experimentally estimated Flory–Huggins parameter between them. This study newly identifies a potential of using polymer grain encapsulation for enhancing thermal stability of OIHPs and their solar cell performance by suppressing the out-diffusion of cationic organic components.

36 MATERIALS SCIENCE↗

Implications of Shear and Thermal Stratification on Wind Turbine Tip-Vortex Stability

The interaction between wind turbines in a wind farm through their wakes is a phenomenon that has been studied for decades and is still relevant today. Turbines clustered together in arrays will often operate in the wake of other upstream turbines which may lead to significant power losses and fatigue loads. For modern large-scale wind turbines, the mean shear velocity profile and thermal stratification are major components of the atmospheric boundary layer so it is important to understand their impact on near-wake development. Additionally, veer is present due to the rotation of the Earth. The impact of shear, thermal stratification and veer on the stable wake length of turbines with a dynamic control strategy is studied numerically in this work using a suite of highly resolved large-eddy simulations. Instantaneous flow fields are extracted from the simulations and used to conduct proper orthogonal decomposition (POD) and compute the mean kinetic energy fluxes by different POD modes to better understand the tip-vortex instability mechanisms. Our findings show that the dynamic pitch control scheme is able to shorten the stable wake length to about 1.5R in uniform flow. Shear can significantly affect the break up of wind turbine tip-vortices as well as the shape and stable length of the wake, whereas thermal stratification seems to only have limited contribution to the spatial development of the near-wake field. Veer causes the wake boundary to skew but has a limited impact on the wake length.

length thermal stratification↗

Novel architectures for stabilization of Mn-rich cathodes: a high-valent approach to interfaces

Lithium- and manganese-rich (LMR) layered oxides continue to generate significant interest as promising, earth-abundant cathode materials for next-generation Li-ion batteries. In spite of their attractive capacity and cost advantages, a few long-standing challenges still hamper their widespread adoption, with manganese dissolution being one of the most persistent and vexing issues. In the present study, we explore the incorporation of Sb5+ as a high-valent cation and exploit its ability to form unique lithium-rich surface and grain-boundary structures that can integrate directly with the LMR lattice. When synthesized under appropriate conditions, Sb5+ orders strongly with Li+ to form localized Li+–Sb5+ motifs, which play a key role in restructuring the surface and grain-boundary regions. These restructured regions act as protective, stabilizing entities that substantially suppress electrolyte-driven side reactions, reduce impedance growth, limit manganese dissolution, and help retain cyclable lithium during long-term electrochemical cycling. Further improvements of the electrochemical performance of Sb-treated LMR were achieved using a well-known additive to mitigate Mn dissolution and highlight the synergistic effects of combined strategies. Overall, this work showcases how high-valent elements such as Sb5+ can help tailor the surface and intergranular regions and work in synergy with other modifiers (e.g. electrolyte additives), enhancing the cycle life and practical viability of LMR cathodes for use in graphite-based full cells.

Mallick, Subhadip [Argonne National Laboratory (AN↗

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SEI Formation and Lithium-Ion Electrodeposition Dynamics in Lithium Metal Batteries via First-Principles Kinetic Monte Carlo Modeling

The stabilization and enhanced performance of lithium metal batteries (LMBs) depend on the formation and evolution of the Solid Electrolyte Interphase (SEI) layer as a critical component for regulating the Li metal electrodeposition processes. This study employs a first-principles kinetic Monte Carlo (kMC) model to simulate the SEI formation and Li + electrodeposition processes on a lithium metal anode, integrating both the electrochemical electrolyte reduction reactions and the diffusion events giving place to the SEI aggregation processes during battery charge and discharge processes. The model replicates the competitive interactions between organic and inorganic SEI components, emphasizing the influence of the cycling regime. Results indicate that grain boundaries within the SEI facilitate faster lithium-ion transport compared to crystalline regions, crucial for improving the performance and stability of LMBs. The findings underscore the importance of dynamic SEI modeling for further development of next-generation high-energy-density batteries.

25 ENERGY STORAGE↗