Search NASA⌕ Search

SEARCH · Search NASA

Results for “Contact prediction”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Minimal entanglement and emergent symmetries in low-energy QCD

Here, we study low-energy scattering of spin-$\frac{1}{2}$ baryons from the perspective of quantum information science, focusing on the correlation between entanglement minimization and the appearance of accidental symmetries. The baryon transforms as an octet under the SU⁡(3) flavor symmetry and its interactions below the pion threshold are described by contact operators in an effective field theory (EFT) of QCD. Despite there being 64 channels in the 2-to-2 scattering, only six independent operators in the EFT are predicted by SU⁡(3). We show that successive entanglement minimization in SU⁡(3)-symmetric channels are correlated with increasingly large emergent symmetries in the EFT. In particular, we identify scattering channels whose entanglement suppression are indicative of emergent SU⁡(6), SO⁡(8), SU⁡(8), and SU⁡(16) symmetries. We also observe the appearance of non-relativistic conformal invariance in channels with unnaturally large scattering lengths. Improved precision from lattice simulations could help determine the degree of entanglement suppression, and consequently the amount of accidental symmetry, in low-energy QCD.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Development of Numerical Model of Metal Foam with PCM for the Estimation of Effective Thermal Conductivity

Global warming due to climate change is a threat to humankind. Nuclear energy is one of the promising solutions to reduce fossil fuel usage. Nuclear energy can handle the base load, compensating for the volatility of renewable energy. If nuclear energy could achieve load following capability, the combination with renewable energy would be more suitable. Thermal energy storage (TES) is one of the options for enabling load following of nuclear reactors. The TES makes it possible to store surplus nuclear thermal energy and release it later as needed. In Idaho National Laboratory (INL), a new concept of latent heat TES integrated with high-temperature heat pipe has been proposed and is under development, which is called Heat pipe-Integrated Thermal Battery (HITB). HITB exchanges thermal energy between the reactor system and TES via heat pipe. The heat transferred to TES medium, made of phase change material (PCM), stores energy as sensible heat and/or latent heat. As PCM typically has poor thermal conductivity, however, various heat transfer enhancement techniques are required to achieve a rapid charging cycle. There are many techniques to enhance the heat transfer ability of TES medium such as disk, fin, and metal foam. Among them, metal foam is an appropriate option to enhance the heat transfer because it maximizes the heat transfer area through metal wicks. Metal foam is a lightweight metal structure that has a high porosity of over 0.9. The typical materials for metal foam are Aluminum, Copper, Nickel, and Silicon Carbide (SiC). Metal foam not only enhances heat transfer via conduction but also increases contact surface area. In the HITB design , the metal foam is being considered as one of the options to enhance the heat transfer of TES medium (PCM) [1]. To predict the enhanced thermal performance of TES, one should properly estimate the effective thermal conductivity of metal foam combined with PCM material or calculate heat transfer in distributed model. There are many experimental works that provides effective thermal conductivity of metal foam with various PCM [2,3]. Also, many theoretical models were developed based on the unit cell model of metal foam [4,5]. With a distributed model, on the other hand, detail heat transfer characteristics between metal foam and PCM material can be analyzed considering the geometry or buoyancy effect. However, due to the complex geometry of metal foam pores, the computational cost for three-dimensional modeling highly increases. Therefore, if metal foam structure can be modeled in simple and repetitive design, the computational cost would decrease Among the various metal foam models [2], lattice model is one of the simple and extendable design. The porosity and pores per inch (PPI) can be characterized by the size and spatial distance of lattice structure. If the three-dimensional metal foam model consists of lattice structure could properly estimate the heat transfer, which is characterized by effective thermal conductivity, it would be a good option to assess the thermal performance of metal foam with PCM. In this study, a three-dimensional numerical model was developed to simulate conductive heat transfer between metal foam and PCM. The three-dimensional lattice structure of square pillars was selected as a basic structure of the metal foam. The calculation result was characterized by the effective thermal conductivity of the whole domain. A sensitivity study was conducted for mesh size, domain size, and PPI to check whether the calculation result gives a converged result or not. Lastly, the effective thermal conductivity from the lattice model was compared with existing experimental data to validate the model result

25 ENERGY STORAGE↗
DOE OSTI2022

Thermo-Mechanical Distortion of Tungsten-Coated Steel During High Heat Flux Testing Using Plasma Arc Lamps

An experimental setup and a test section were designed and fabricated for high heat flux testing (HHFT) of neutron-irradiated specimens using water-wall plasma arc lamps. Because of the radiological considerations and limitations of reactor irradiation, the size of the test articles was limited to disks less than 10 mm in diameter. The specimen was clamped onto an actively cooled block, and clamping allowed the insertion of several thermocouples on the back surface of the specimen through a copper (Cu) block. Five vacuum plasma sprayed tungsten (W)–coated F82H steel specimens were subjected to HHFT. Surface profilometry measurements, which were conducted after HHFT, revealed central bowing of the top W surface. This type of residual distortion occurred for all of the specimens, and the larger the specimens were, the larger was the distortion.In an attempt to understand specimen distortion and address the science questions related to the testing of subsize specimens during HHFT, a simplified thermo-mechanical model was developed. By using a measured temperature in the Cu as an isothermal boundary condition, the model eliminated the need for coupling cooling fluid flow models with stress models, greatly simplifying the analysis. The main variable in the proposed model is hC, i.e., the thermal contact conductance between the F82H and the Cu washer. Inelastic properties, including hardening properties, were considered for F82H steel and Cu. Numerical simulation results demonstrated a buildup of residual deformation during HHFT and a very complex state of stress and deformation during typical heat flux (HF) cycling. Additionally, hoop stress evolution during a high heat flux cycle reveals that F82H at an interface with W would be mainly in compression during HF application and experienced a transition to a tension state during cooldown. Also, specimen distortion evolves during each HF cycle, as the specimen bows downward during HF application and upward during the cooldown period between HF cycles. The final specimen distortion, i.e., upward bowing of the specimen center, was qualitatively predicted for hC values of 4000 to 5000 W/(m 2 ·K). This hC range of values, for which bulging is obtained, is at the lower spectrum of the range of values for hC, consistent with the low thermal contact conductance expected from the unpolished F82H surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effect of plasticity on the experimental characterization of the creep property of a cement paste: Spherical vs. conical microindentation

Microindentation is today a powerful tool for rapidly characterizing the elastic modulus and the logarithmic basic creep of a cement paste. However, a theoretical work proved that a microindentation test with a sharp indenter shape involves stress concentration and delayed plasticity which can cause an overestimation of the estimated creep property. This work aims at experimentally investigating the effect of the shape indenter on the identification of the visco-elastic property of a cement paste by microindentation. First, several microindentation grids were carried out on a cement paste with both a Berkovich indenter and a spherical indenter at different penetration depths or load levels for load control (creep) tests or displacement control (relaxation) tests, respectively. As for a visco-elastic material the duality principle between creep and relaxation holds, a numerical method was developed for predicting the relaxation response from the assumed logarithmic creep compliance. The present results show that only spherical microindentation satisfy the duality principle for which creep and relaxation response can be predicted by a unique logarithmic creep compliance function. Finally, spherical microindentation creep tests at low load as well as spherical microindentation relaxation tests provide a greater indentation modulus M and contact creep modulus C{sub v} as likely less affected by possible plasticity effects.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluating route to impact convergence of the harmonic balance method for piecewise-smooth systems

Here in this work, we investigate the applicability of the harmonic balance method (HBM) to predict periodic solutions of a single degree-of-freedom forced Duffing oscillator with freeplay nonlinearity. By studying the route to impact, which refers to a parametric study as the contact stiffness increases from soft to hard, the convergence behavior of the HBM can be understood in terms of the strength of the non-smooth forcing term. HBM results are compared to time-integration results to facilitate an evaluation of the accuracy of nonlinear periodic responses. An additional contribution of this study is to perform convergence and stability analysis specifically for isolas generated by the non-smooth nonlinearity. Residual error analysis is used to determine the approximate number of harmonics required to get results accurate to a given error tolerance. Hill’s method and Floquet theory are employed to compute the stability of periodic solutions and identify the types of bifurcations in the system.

42 ENGINEERING↗
DOE OSTI2022

Influence of intermolecular interactions on the infrared complex indices of refraction for binary liquid mixtures

This paper investigates the accuracy of deriving the composite optical constants of binary mixtures from only the complex indices of refraction of the neat materials. These optical constants enable the reflectance spectra of the binary mixtures to be modeled for multiple scenarios (e.g., different substrates, thicknesses, volume ratios), which is important for contact and standoff chemical detection. Using volume fractions, each mixture's complex index of refraction was approximated via three different mixing rules. To explore the impact of intermolecular interactions, these predictions are tested by experimental measurements for two representative sets of binary mixtures: (1) tributyl phosphate combined with n-dodecane, a non-polar medium, to represent mixtures which primarily interact via dispersion forces and (2) tributyl phosphate and 1-butanol to represent mixtures with polar functional groups that can also interact via dipole–dipole interactions, including hydrogen bonding. The residuals and the root-mean-square error between the experimental and calculated index values are computed and demonstrate that for miscible liquids in which the average geometry of the cross-interactions can be considered isotropic (e.g., dispersion), the refractive indices of the mixtures can be modeled using composite n and k values derived from volume fractions of the neat liquids. Conversely, in spectral regions where the geometry of the cross-interactions is more restricted and anisotropic (e.g., hydrogen bonding), the calculated n and k values vary from the measured values. The impact of these interactions on the reflectance spectra are then compared by modeling a thin film of the binary mixtures on an aluminum substrate using both the measured and the mathematically computed indices of refraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Continuous snow depth and temperature measurements from dense network of above-ground distributed temperature profiling systems from 2021-09-23 to 2024-08-23, Seward Peninsula, Alaska

The dataset contains temperature measurements from distributed temperature profiling (DTP) systems (Dafflon et al., 2022; Wielandt et al., 2022; Wang et al., 2024a; Fiolleau et al., 2024) deployed vertically above the ground surface at a large number of locations from 2021 to 2024. The research is designed to improve understanding of the local heterogeneity in snow depth and snow thermal insulation dynamics, as well as their interactions in a discontinuous permafrost region (Wang et al., 2025). The DTP systems were deployed at 96 locations in a watershed along the Nome-Teller road at mile marker 27 (T27) and at 54 locations on a hillslope along the Kougarok road at mile marker 64 (K64) in the Seward Peninsula, Alaska. The probe location information is stored in Probe_locations_T27.csv and Probe_locations_K64.csv. Temperature measurements were recorded at 15-minute intervals using high-precision digital sensors (accuracy: ±0.1°C, resolution: 0.0078°C). The temperature probes, either 1.4 m or 1.6 m long, contain sensors spaced every 5 cm or 10 cm along their length. The temperature data are stored in compressed files following the format: DTP_snow_air_temperature_(site)_(start)_(end).zip, where site is either T27 or K64, and start and end represent the time series period. Within each ZIP file, individual CSV files are named by probe ID and contain temperature records at different heights above the ground surface.This dataset also includes derived snow depth time series over three snow seasons, estimated from temperature measurements. Snow depth was estimated by identifying the consecutive sensor pair that exhibited the largest drop in high-frequency temperature fluctuations (detailed in the methods). These data are stored in: Snow_depths_flags_(site)_(start)_(end).csv, which includes snow depth time series and corresponding quality flags (defined in the methods) from different probes. Additionally, the dataset includes derived metrics and supporting measurements at selected locations over two snow seasons, contributing to the manuscript of Wang et al., 2025. These locations were chosen based on the availability of high-quality snow depth time series during both seasons. The additional data include: (1) Air temperature proxies measured from the top sensors on the pole when they were not buried by snow, stored in Air_temperature_proxies_(site)_(start)_(end).csv (2) Ground interface temperature, recorded at 3 cm above the ground, stored in Ground_interface_temperature_(site)_(start)_(end).csv (3) Site characteristics, including vegetation height, elevation, and the topographic position index (TPI) within a 50 m radius, stored in Selected_probe_locations_gps_vegheight_tpi_elevation_(site).csv. These metrics were derived from 1 m resolution summer LiDAR-based digital elevation models and digital surface models from Singhania et al., 2023, DOI:10.5440/1832016. Metadata files include data descriptions (_dd.csv) for tabular data. All included files are listed and described in xxxx_flmd.csv.This dataset is an updated version of a previous archive (Wang et al., 2024b, DOI: 10.15485/2475020), incorporating multiple seasons and improved snow depth estimation. Please note that due to large amount of information present in this dataset, many specificities associated with the acquisition of snow temperature, air temperature proxy and estimation of snow depth, and the future archiving of additional datasets on the soil temperature, thaw depth and soil characteristics at these locations, the author would welcome being contacted by people planning to use this dataset.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Characterization and Damping Control of Mechanical Connections to Improve Performance of Horn Stripline

Magnetic focusing horns are critical components for creating a stable beam of neutrinos for neutrino facilities, such as the Long Baseline Neutrino Facility (LBNF) and the Neutrinos and Main Injector (NuMI) beam lines at Fermilab. The pulsed magnetic horns are powered by high current electricity through long striplines. In addition to requirements for low inductance and voltage standoff, the striplines must survive in a harsh radiation environment for the operational life of the component, specified as the 100 million pulses requirement for LBNF. Each stripline assembly consists of four (eight layers that splits off to 2 pairs of four layers at the horn interface) layers of Al 6101-T6. As an electromechanical system, the stripline layers are bolted together with ceramic isolators for electrical insulation and connections must facilitate passive cooling and mechanical stability. The striplines experience vibrational force in addition to clamping force, repetitive thermal and electro-magnetic loading. The holes of stripline plates for ceramic joints represent one of the weak links for potential failure [1,2]. By characterizing the contact behavior of these joints and optimizing their damping properties through finite element analysis and experimental modal analysis, stripline performance and longevity can be improved. This study not only helps predict the behavior of the striplines but also improves their performance to meet the required operational lifetime.

Liu, Zunping [Fermilab]↗
DOE OSTI2021

Benchmarking the multiconfigurational Hartree method by the exact wavefunction of two harmonically trapped bosons with contact interaction

Highlights: • Comparison of exact and numerical single-particle density matrices. • Benchmark for multiconfigurational Hartree. • Convergence for large number of orbitals. We consider two bosons in a one-dimensional harmonic trap, interacting by a contact potential, and compare the exact solution of this problem to a self-consistent numerical solution by using the multiconfigurational time-dependent Hartree (MCTDH) method. We thereby benchmark the predictions of the MCTDH method with a few-body problem that has an analytical solution for the most commonly experimentally realized interaction potential in ultracold quantum gases. It is found that exact ground state energy and first order correlations are accurately reproduced by MCTDH up to the intermediate dimensionless coupling strengths corresponding to typical background scattering lengths of magnetically trapped ultracold dilute Bose gases. For larger couplings, established for example by (a combination of) Feshbach resonances and optical trapping, the MCTDH approach overestimates the depth of the trap-induced correlation dip of first order correlations in position space, as well as overestimates the fragmentation, defined as the average relative occupation of orbitals other than the energetically lowest one. We anticipate that qualitatively similar features in the correlation function may arise for larger particle numbers, paving the way for a quantitative assessment of the accuracy of MCTDH by experiments with ultracold atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Effect of the numerical discretization scheme in Shock-Driven turbulent mixing simulations

In this work, we evaluate the effects of distinct numerical strategies in simulations of shock-driven turbulent mixing. An air-SF 6 -air gas curtain subjected to Mach 1.2 shock-waves is computed with Implicit Large-Eddy Simulation based on three numerical schemes: directional-split with Harten-Lax-van Leer (HLL) solver; directional-unsplit using HLL-Contact (HLLC) solver; and directional-unsplit with HLLC solver and a Low Mach number Correction (LMC). The results illustrate the importance of the numerical strategy to the accuracy of the predictions. Whereas both split and unsplit schemes result in a similar spatial development of the initial shock-driven instability, only the unsplit schemes can predict the turbulent mixing transition after reshock observed by the laboratory experiments. Such feature increases the mixing rate of the two fluids, this being particularly pronounced when the LMC is active due to i) the reduced flow Mach number, and ii) the larger effective Reynolds number. Since the selected mixing problem is driven by the deposition of vorticity at the fluids’ interface, the resultant flow physics is analyzed by investigating the contribution of distinct inviscid mechanisms to the production of vorticity: baroclinicity, stretching, and dilatation. As expected, it is observed that the production of vorticity is initially dominated by the baroclinicity mechanism. Yet, the relevance of the remainder mechanisms is enhanced after the reshock and may even surpass the baroclinicity term.

42 ENGINEERING↗
DOE OSTI2022

ENIGMA: Ecosystems and Networks Integrated with Genes and Molecular Assemblies. Final Report

ENIGMA research is performed at a very large scale involving collaborations between multiple research groups to address two major driving concepts: For any microbe, it is possible to discover efficiently the genetic determinants of adaptation to life in dynamic environments in contact with other living members of its ecology and reveal hidden functional phenomena that are only possible by using a systems biology approach; and for any environment, it is possible to predict the temporal changes in geochemistry and accompanying system behavior, for some length of time, given biogeochemical and physical inputs. These driving concepts also map directly to grand challenges the ENIGMA team has outlined: “Define the levels of biological organization most relevant to scaling from single cells to ecosystems and global cycles; capture how that organization varies in time and space; and identify critical interactions that dictate rates of carbon, nutrient, and energy transformation in different environments.”

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Synthesis, crystal and electronic structure of the Zintl phase Ba 16 Sb 11 . A case study uncovering greater structural complexity via monoclinic distortion of the tetragonal Ca 16 Sb 11 structure type

The binary Zintl phase Ba 16 Sb 11 has been synthesized and structurally characterized. Detailed studies via single-crystal X-ray diffraction methods indicate that although Ba 16 Sb 11 appears to crystallize in the tetragonal Ca 16 Sb 11 structure type (space group $P\bar{4}2_1$m with a=13.5647(9) Å, c=12.4124(12)Å, Z=2, R 1 = 3.14%; wR 2 = 4.77%), there exists an extensive structural disorder. Some Ba 16 Sb 11 crystals were found to be monoclinic and the structure was solved and refined in space group P2 1 (a=18.3929(12) Å, b=13.5233(8) Å, c=18.3978(12) Å, β=94.6600(10)°; Z=4, R 1 =5.84 %; wR 2 =9.58 %). The latter corresponds to a 2-fold superstructure of the tetragonal one, which provides a disorder-free structural model. In both descriptions, the disordered tetragonal and the ordered monoclinic superstructure, the basic building units that make up the structure of this Ba-rich compound are pairs of face-shared square antiprisms of Ba atoms, which are centered by Sb atoms. The dimerized antiprisms are linked into parallel chains via square prisms of Ba atoms, which are also centered by Sb atoms. The Zintl concept can be applied in a straightforward manner and as result, the structure of Ba 32 Sb 22 (=2×Ba 16 Sb 11 ) can be rationalized as (Ba 2+ ) 32 (Sb 3– ) 20 [Sb 2 ] 4– . Notably, the partitioning of the valence electrons is done taking into an account the homoatomic Sb–Sb contacts (d=3.01 Å), which can be clearly distinguished in the lower symmetry space group. Electronic structure calculations of Ba 16 Sb 11 are in good accordance with the Zintl rationalization and predict a semiconductor with a band gap of 0.77 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Learning-based approach to plasticity in athermal sheared amorphous packings: Improving softness

The plasticity of amorphous solids undergoing shear is characterized by quasi-localized rearrangements of particles. While many models of plasticity exist, the precise relationship between the plastic dynamics and the structure of a particle’s local environment remains an open question. Previously, machine learning was used to identify a structural predictor of rearrangements called “softness.” Although softness has been shown to predict which particles will rearrange with high accuracy, the method can be difficult to implement in experiments where data are limited and the combinations of descriptors it identifies are often difficult to interpret physically. Here, we address both of these weaknesses, presenting two major improvements to the standard softness method. First, we present a natural representation of each particle’s observed mobility, allowing for the use of statistical models that are both simpler and provide greater accuracy in limited datasets. Second, we employ persistent homology as a systematic means of identifying simple, topologically informed, structural quantities that are easy to interpret and measure experimentally. We test our methods on two-dimensional athermal packings of soft spheres under quasi-static shear. We find that the same structural information that predicts small variations in the response is also predictive of where plastic events will localize. We also find that an excellent accuracy is achieved in athermal sheared packings using simply a particle’s species and the number of nearest neighbor contacts.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Excess carrier concentration in silicon devices and wafers: How bulk properties are expected to accelerate light and elevated temperature degradation

Light and elevated temperature induced degradation (LeTID) is accelerated nearly linearly by the presence of excess carriers. It is therefore important to understand how excess carrier concentration (Δn) changes as a function of exposure conditions, materials properties, and sample structure. We simulate Δn as a function of wafer thickness and bulk minority carrier lifetime (τ) in solar cells and wafers using SCAPS and Quokka3. We also derive closed-form analytic expressions. For wafers, there is a near-linear relationship between Δn and τ or thickness, whereas for solar cells, Δn in the bulk may become limited by rear surface recombination. Thus, LeTID may progress more quickly in wafers than in cells, with a stronger dependence on τ. When comparing experiments, observed degradation rates must be corrected between samples or conditions to account for differences in Δn. This study demonstrates three tools to estimate the magnitude of such corrections, which can aid in the quantitative interpretation of LeTID data and performance predictions. Finally, while each tool yields similar results, there are advantages to each approach that must be weighed in terms of simplicity of inputs versus sophistication of treatment. Incomplete specification of back contact characteristics in commercial products is identified as an important contributor to uncertainty in expected LeTID rates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Amorphization and siliconization of silicon carbide as a first wall material

The understanding and prediction of silicon carbide (SiC) material evolution exposed to SOL plasma conditions is of prime interest because SiC represents a promising main chamber wall plasma-facing material for next-step fusion devices (low hydrogenic diffusion, good mechanical and thermal properties under neutron irradiation). Gross and net Si erosion rates from SiC surfaces in contact with a well-diagnosed L-mode plasma in the DIII-D tokamak have been simulated and the surface concentrations of impurities have been tracked as a function of time. Coupled simulation of surface model and impurity transport demonstrates amorphization of crystalline SiC exposed to L-mode plasma due to the accumulation of displacement damages under ion irradiation. This affects the lifetime of SiC plasma facing components. Surface evolution is tightly coupled to impurity transport in the plasma and therefore needs to be integrated with impurity transport simulations to effectively predict Si erosion rates and sub-surface concentrations as a function of time. The simulation workflow couples a semi-analytical surface model to the impurity transport code GITR. The surface model is a homogeneous mixed-material model that tracks physical & chemical sputtering and reflection of impurities. Gross erosion is primarily influenced by the background plasma parameters and redeposition patterns are mainly influenced by the prompt redeposition due to the gyro-orbits of impurity ions. Although crystalline form of SiC is preferable for fusion wall applications because of resistance to neutron irradiation, this work indicates that crystalline SiC will undergo amorphization under D plasma contact with implications of higher sputtering and fuel retention. These results direct us to explore the effects of amorphization on crystalline SiC and further the physics basis of SiC usage as first wall material for fusion environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗
DOE OSTI2020

Quantitative Understanding and Implementation of Screen-printed p+ Poly-Si/Oxide Passivated Contact to Enhance the Efficiency of p-PERC Cells

Investigate the passivation quality of p-TOPCon to evaluate the promise and viability of upgraded p-PERC with p-TOPCon back. Model and simulate the p-PERC and p-TOPCon cells to predict the gain in cell efficiency by implementation of the p-TOPCon Develop a simple process sequence and fabricate p-TOPCon-PERC cell with ≥ 21.0% efficiency

p-TOPCon, P-PERC↗
DOE OSTI2025

Diffusion and Deformation Mechanism Maps in Li Metal from Atomistic Simulations

The rate of Li transport in the Li metal anode is important for the operation of Li metal-solid state batteries. Here, transport rates due to diffusion and creep are predicted in Li using atomistic simulations. First, molecular dynamics is used to estimate the rate of Li diffusion along dislocations and in grain boundary triple junctions. By combining this data with that from a prior study of grain boundary diffusion the dominant mechanisms and rates of self-diffusion in Li polycrystals are predicted as a function of grain size, grain shape, dislocation density, and temperature. Second, the dominant creep mechanisms are predicted and used to estimate critical current densities and void annihilation times. Grain boundary sliding and coble creep are the dominant mechanisms for micron-sized grains. Lastly, a continuum model for interfacial contact loss reveals that high dislocation densities of ∼10 12 /cm 2 enable achieving battery performance targets for Li grain sizes of ∼10 μm.

Batteries↗