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 55 records · Page 3

Assessment of Combustion Residual Leachate Volume, Composition, and Treatment Costs

Combustion residuals and the resulting leachate from storage sites represent a large volume of wastewater in the United States (U.S.) that has not been quantified. Here this work estimates the constituents present, volume of wastewater, and costs of treatment for both combustion residual landfill leachate and the leachate from surface impoundment closures. Combustion residual landfill leachate produced from contact with bituminous coal combustion byproducts is generally predicted to be higher in lithium and manganese, whereas landfill leachate produced from contact with subbituminous coal combustion byproducts is generally predicted to be higher in mercury and vanadium. The annual volume of a single landfill with combustion residual leachate can reach more than 800,000 cubic meters. This leachate represents an annual volume of 26.8–42.8 million cubic meters nationally. Closing surface impoundments can yield between 830 and 1040 cubic meters of leachate nationally for a three-year closure period. Costs as low as $1.5/m 3 or as high as $95/m 3 are observed. Treatment trains will need to remove 72% of total suspended solids (TSS), 87% of arsenic, and 64% of mercury from landfill leachate. When applied to impoundments, these treatment trains would need to remove 97% of arsenic.

01 COAL, LIGNITE, AND PEAT↗

Accelerated Beta Radiation Aging of Interlayer Titanium Nitride in Gallium Nitride Contacts

Research conducted under the Laboratory Directed Research and Development Program at Pacific Northwest National Laboratory, a multiprogram national laboratory operated by Battelle for the U.S. Department of Energy. An NNSS author co-authored this paper, and portions of this work were supported by the Mission Support and Test Services, LLC, under Contract No. DE-NA0003624 with the U.S. Department of Energy, National Nuclear Security Administration, NA-10 USDOE NA Office of Defense Programs (NA-10). PNNL-SA-163618 Currently, there is no good way to determine the influence of radiation on the aging of betavoltaic electrical contacts. This work tested a method to accelerate the aging of the contacts inside a betavoltaic by enhancing energy deposition within the interfacial region of interest. An electrical contact of gold-titanium on gallium nitride was aged by exposure to tritium betas and electrons in a beamline. The interface stoichiometry was compared to the electrical performance of the contact. The method to age the betavoltaic component could help predict the lifetime performance of the internal electrical contacts and provide assurance of deployment reliability.

36 MATERIALS SCIENCE↗

Lanthanum Strontium Cobaltite-Alumina Fiber Composite Contact Material for Solid Oxide Fuel Cells

Short alumina fibers were added to LSCo contact material to enhance the contact strength and improve thermal cycle stability. Alumina fiber from 2.5 v% to 20 v% were mixed with LSCo and sintered at elevated temperatures. Sintering study showed retardation in densification of LSCo matrix with increasing alumina content. Sintered bulk strength exhibited similar behavior where strength greatly reduced at high alumina v%. Coefficient of thermal expansion showed consistent results with prediction from rule of mixture. Contact strength using anode supported YSZ bilayer coupons improved greatly with the presence of small amount (2.5-5 v%) alumina fibers and no substantial reduction after 10 deep thermal cycles. Candidate LSCo-Al2O3 contact materials was validated in a generic stack test fixture for 1000h at 800oC followed by 10 thermal cycles. Results of impedance analysis showed substantial improvement in thermal cycle stability as compared to LSCo parent materials.

SOFC cathode, strength, LSCo, alumina↗

Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers

Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.

Alkaline↗

Impact of feed variability on cesium removal with multiple actual waste samples from the Hanford site

Here, the Tank Side Cesium Removal (TSCR) system, under development by Washington River Protection Solutions (WRPS), will send initial low-activity Hanford waste tank supernate feeds to the Hanford Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Facility. In addition to entrained solids removal from the supernate, the primary goal of TSCR is to remove cesium-137 ( 137 Cs) by ion exchange, allowing contact handling of the liquid effluent product at WTP. Crystalline silicotitanate (CST) ion exchange media, manufactured by Honeywell UOP (product IONSIV™ R9140-B), was selected as the ion exchange media at TSCR. This lot of material was found to have superior performance in comparison to historic batches of CST. Ion exchange column and batch contact testing with supernate from Hanford tanks AP-105, AP-107 and AW-102 was performed to assess the impact of feed variability on system performance. These tests demonstrated that batch contact measurements provide a reasonable prediction of column capacities with some deviation in performance attributed to column dynamics. In addition, the variability in CST capacity for cesium in the various actual waste samples does not track with historical understanding of the competing cations, suggesting that tank waste samples contain other components that may significantly impact cesium loading.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Diffusion energy barrier of Au on Bi 2 Se 3 : theory and experiment

The stability and diffusion of ultra-thin thermally deposited Au films on Bi 2 Se 3 was studied using scanning tunneling microscopy and density functional theory calculations. The Au/Bi 2 Se 3 interface is of interest as gold is predicted to provide excellent electrical contact while maintaining the spin-polarized characteristics of the electronic states in Bi 2 Se 3 that make the material attractive for spintronic applications. When deposited at room temperature, Au 10 covers the surface with tightly packed islands of nanometer scale dimension. The surface morphology is stable up to 400K. At this annealing temperature, Au atoms have sufficient energy to diffuse across the surface and aggregate into larger nanostructures. At 550K, the Bi 2 Se 3 surface is only sparsely covered, and the Au has formed clusters with length scales 5-10 times larger than the original islands formed at room temperature. Comparison of the experiment and first principle calculation lead to the conclusion that the diffusion energy barrier for Au on Bi 2 Se 3 is as high as 0.47 eV, 15 which is much larger than diffusion barriers on other van der Waals materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A new method for measuring refractory corrosion of ceramics in glass

Abstract Nuclear waste glass vitrification furnaces are lined with refractory ceramic blocks to contain the molten glass. The refractory liner is susceptible to corrosion and has a finite service lifetime. For this reason, predicting the refractory corrosion in contact with molten glass is integral to estimating melter service lifetime. Standardized laboratory tests varying time and temperature are commonly performed to estimate refractory material loss as a function of glass composition. These data are time and resource‐intensive to collect and are susceptible to considerable measurement error. In order to accelerate glass formulation and design for nuclear waste vitrification, methods are needed to increase laboratory‐scale throughput while maintaining data quality. In this work, a method to remove the residual glass from a corroded coupon using hydrofluoric acid is presented that accelerates the throughput of sample analysis while simultaneously facilitating more accurate measurements.

Amoroso, Jake W. [Savannah River National Laborato↗

Leaf gas exchange and fitted parameters, two sites in Panama, 2022

Photosynthetic CO2 response curves (ACi curves), light response curves (AQ curves), dark adapted dark respiration, conductance curves, and survey measurements for leaves measured in the Parque Natural Metropolitano (PNM) and Gamboa, Panama, from January to April 2022 are presented. Measurements were made on leaves from 46 different tree, shrub and liana species, from sunlit canopy and understory locations. Leaf area index from 8 vertical profiles at PNM are also included. The aim of this measurement campaign was three-fold: (i) to improve our understanding of the vertical variation in leaf-level water use efficiency; (ii) to develop an understanding of the sensitivity of stomata to changes in environmental conditions, especially light, humidity, and temperature; (iii) to improve models which can predict leaf traits from leaf contact spectral measurements. All the gas exchange data and metadata are presented in .csv files and complete instrument output are included in .zip folders. Data and metadata meet the ESS-DIVE leaf-level gas exchange reporting format requirements. The protocol details are provided as pdf documents. In addition to gas exchange data reported here these samples were also used for measurement of leaf optical properties, carbon and nitrogen content, and leaf mass per unit leaf area (LMA). These data can be cross linked using the unique sample ID and are provided in separate data packages.

54 ENVIRONMENTAL SCIENCES↗

Leaf reflectance, leaf transmittance and sky irradiance, at four sites in Panama, 2022

Paired leaf reflectance and transmittance data measured in Panama, at the Parque Natural Metropolitano (PNM), on Barro Colorado Island (BCI), at the San Lorenzo canopy crane site (SLZ), and in Gamboa, Panama from January to April 2022. Measurements were made on leaves from 57 different tree, shrub and liana species, from sunlit canopy and understory locations. Bark reflectance was also measured on the trunks of common tree species at the PNM and SLZ sites. Also included are irradiance measurements and leaf area index (LAI) from eight vertical profiles at PNM. The aim of these data was to refine our understanding of vertical trait and environmental variation in tropical forests, and to improve models which can predict leaf traits from leaf contact spectral measurements. All spectral data and metadata are presented in .csv files and complete instrument output are included in .zip folders. Protocol details are provided as pdf documents. In addition to spectral data described here, some of these samples were also used for measurement of leaf gas exchange, carbon and nitrogen content, and leaf mass per unit leaf area (LMA) these data can be cross linked using the unique sample ID and are provided in separate related data packages.

54 ENVIRONMENTAL SCIENCES↗

Foliar isotopic and elemental biochemistry and leaf mass per area (LMA), Parque Natural Metropolitano, Panama, 2022

Foliar N and C isotope composition, bulk elemental N and C composition, and leaf mass per area (LMA) for leaves sampled from the Parque Natural Metropolitano (PNM), Panama, from January to April 2022 are presented. Measurements were made on leaves from 17 different tree, shrub and liana species, from sunlit canopy and understory locations on 8 vertical profiles. Full sample details and leaf area index (LAI) from the vertical profiles at PNM are also included in the data package. The aim of this measurement campaign was two-fold: to improve our understanding of the vertical variation in leaf-level water use efficiency, and to improve models which can predict leaf traits from leaf contact spectral measurements. Biochemistry data and sample metadata are presented in .csv files. The original isotopic data report, which includes details about standards and data accuracy, is provided in .xlsx format. Data and metadata meet the ESS-DIVE reporting format requirements for file level metadata (FLMD) and comma separated values (csv). The protocol details are provided as pdf documents. In addition to foliar biochemistry data reported here these samples were also used for measurement of leaf gas exchange, and leaf optical properties. These data can be linked using the unique sample ID and are provided in separate data packages (NGT0192 and NGT0193).

54 ENVIRONMENTAL SCIENCES↗

Magnetic braking saturates: evidence from the orbital period distribution of low-mass detached eclipsing binaries from ZTF

ABSTRACT We constrain the orbital period (Porb) distribution of low-mass detached main-sequence eclipsing binaries (EBs) with light-curves from the Zwicky Transient Facility (ZTF), which provides a well-understood selection function and sensitivity to faint stars. At short periods (Porb ≲ 2 d), binaries are predicted to evolve significantly due to magnetic braking (MB), which shrinks orbits and ultimately brings detached binaries into contact. The period distribution is thus a sensitive probe of MB. We find that the intrinsic period distribution of low-mass (0.1 ≲ M1/M⊙ < 0.9) binaries is basically flat (${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^0$) from Porb = 10 d down to the contact limit. This is strongly inconsistent with predictions of classical MB models based on the Skumanich relation, which are widely used in binary evolution calculations and predict ${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^{7/3}$ at short periods. The observed distributions are best reproduced by models in which the magnetic field saturates at short periods with a MB torque that scales roughly as $\dot{J}\propto P_{\rm orb}^{-1}$, as opposed to $\dot{J} \propto P_{\rm orb}^{-3}$ in the standard Skumanich law. We also find no significant difference between the period distributions of binaries containing fully and partially convective stars. Our results confirm that a saturated MB law, which was previously found to describe the spin-down of rapidly rotating isolated M dwarfs, also operates in tidally locked binaries. We advocate using saturated MB models in binary evolution calculations. Our work supports previous suggestions that MB in cataclysmic variables (CVs) is much weaker than assumed in the standard evolutionary model, unless mass transfer leads to significant additional angular momentum loss in CVs.

El-Badry, Kareem (ORCID:0000000268711752)↗

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗

Predictive numerical modeling of plasma-induced surface roughness and wettability evolution in LM-PAEK/CF tape

Plasma surface modification effectively enhances adhesion in thermoplastic composites, yet its impacts on high-performance polymers like low-melting polyaryletherketone (LM-PAEK) remain inadequately quantified. Here, this study integrates experimental analysis and numerical modeling to characterize surface roughness and wettability changes in LM-PAEK/carbon fiber composites treated with atmospheric plasma. Atomic Force Microscopy quantified surface topography (n = 10 per condition for contact angles), while static contact-angle assessments measured wettability. Roughness rapidly increased from ∼0.2 nm to 1.6 nm, and contact angle reduced from ∼90° to 24°, both stabilizing after 25–30 s of exposure. A semi-empirical, physics-informed framework was calibrated to these data, coupling surface chemistry via the Owens–Wendt decomposition with topography via the Wenzel roughness factor, and evaluated using out-of-sample (cross-validated) tests, while static contact angle assessments measured wettability. Numerical predictions matched experimental results closely (RMSE <5%, R 2 > 0.95). Incorporating material-specific parameters, the calibrated model supports plasma-treatment optimization and provides quantitative guidance for improving interfacial adhesion in thermoplastic composite manufacturing.

Atomic Force Microscopy↗

In-Silico Design of Next Generation Cellulose-Derived Packaging Materials (CRADA Final Report)

Developing sustainable solutions for single-use packaging is an important objective to combat the environmental crisis of plastics pollution. Most embodiments of cellulose-based packaging materials, including CellophaneTM, are completely biodegradable in both terrestrial and marine environments. However, petroleum-derived alternatives offer some performance advantages for metrics such as moisture barriers and mechanical properties. This project leverages molecular dynamics simulation to investigate how molecular modifications to cellulose-based polymer assemblies impact their material properties. An important performance criterion for the modified materials was to retain biodegradability; thus, modifications by naturally occurring, biodegradable additives were the focus of this study. Specifically, we developed models with xylan and lignin of varying monomeric compositions into the cellulose matrix. The mechanical properties were investigated by performing stress-strain simulations, and the water barrier and hydrophobicity were investigated by simulating the water contact angle. Our findings indicate that the incorporation of xylan into the cellulose matrix tends to increase the mechanical properties with an optimal loading of ~27 wt%. We also predict that orienting the nanoscale directionality of the xylan chains such that they are perpendicular to the cellulose fibrils will dramatically increase mechanical strength. In contrast, the incorporation of lignin tends to weaken the composite at all loadings investigated. Simulations of water contact angle predicted that coating polymers on the surface of the cellulose assembly creates a more hydrophobic surface than incorporating them throughout the matrix. Of the coatings investigated, lignin resulted in the most hydrophobic surface, followed by pectin and keratin, which both imparted modest increases in hydrophobicity. Future experimental work done by Futamura will focus on designing material prototypes to capitalize on the predictions of performance enhancement obtained from molecular modeling. While substantial progress was made by the simulations performed in this project, there still exists a vast parameter space that we were unable to investigate, including branching, functional group decoration, and degree of polymerization of polymer additives. However, the methods developed in this initial investigation will facilitate more rapid evaluation of the impact of molecular characteristics on the performance of biopolymer composite materials and thereby accelerate future materials discovery efforts in this area.

36 MATERIALS SCIENCE↗

Cesium Ion Exchange Testing Using Crystalline Silicotitanate with Hanford Tank Waste 241-AP-107

At the time of this testing, the Low-Activity Waste Pretreatment System (LAWPS) was to provide for the initial production of immobilized low-activity waste by feeding Hanford tank supernate from tank farms to the Hanford Tank Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Facility for immobilization. Washington River Protection Solutions requested that Hanford tank waste collected from tank 241-AP-107 (hereafter called AP-107) be processed using conceived pretreatment steps (suspended solids removal by filtration, Cs removal by ion exchange) then vitrified. A small-scale test platform to demonstrate the solids filtration, Cs removal, and LAW vitrification was constructed and installed at Pacific Northwest National Laboratory. Bench-scale ion exchange testing with approximately 9 L of AP-107 supernate was conducted using crystalline silicotitanate (CST) ion exchange media. The IONSIV R9140-B CST was provided by Honeywell UOP, LLC in 2018 (Batch 2081000057). The ion exchange media was first tested with simulant and was previously described. This report describes the Cs ion exchange batch contact and column test results with the AP-107 tank waste. Batch contact testing helps to evaluate CST performance on tank waste supernate prior to processing it in the ion exchange columns. Batch contacts were performed with the waste at four Cs concentrations at a phase ratio of 200 (liquid volume to exchanger mass) with AP-107. The distribution coefficient (K d ) at the equilibrium condition of 8.57 µg Cs/mL (AP-107 feed condition) was determined to be 669 mL AP-107/g CST. With a CST bed density of 1.00 g/mL, this K d corresponded to a predicted 50% Cs breakthrough of 669 bed volumes (BVs). The Cs load capacity at the equilibrium feed condition was determined to be 7.5 mg Cs/g dry CST. The column testing was prototypic to the intended LAWPS operations in a lead-lag column format, although on a small-scale basis with 10-mL CST beds. The feed was processed downflow through the lead column and then through the lag column at ~2.2 BV/h. Loading continued until the lag column reached the WTP waste acceptance criteria (WAC) for receiving supernatant waste for vitrification (a function of the Na and 137 Cs concentrations). For AP-107, the WAC is 0.114% of the influent 137 Cs concentration; this required a Cs decontamination factor of 876. The Cs effluent from the lag column reached the WAC after processing ~410 BVs. To keep the subsequent product effluent below the WAC, a replacement lag column was prepared, the lead column was removed from service (after processing a total of 471 BVs), the lag column was put into the lead column position, and the replacement lag column was installed. Feed processing continued and after another ~290 BVs the Cs effluent from the lag column again exceeded the WAC. In both cases, the lead columns only reached 25% Cs breakthrough before removal. Although 50% Cs breakthrough was not reached, this value was estimated and averaged based on extrapolation of the loading curves (640 BVs) and agreed within 4% of the predicted 50% Cs breakthrough from batch contact test results (669 BVs). Table ES.1 summarizes the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Design of an Out-Of-Pile Experimental Facility to Demonstrate the Feasibility of In Situ Thermal Conductivity Measurements of Nuclear Fuels Under Irradiation

There is substantial merit in quantifying nuclear fuel performance under irradiation. At Oak Ridge National Laboratory (ORNL), the MiniFuel irradiation platform has become the primary test vehicle for conducting separate-effects fuel performance irradiation experiments. The MiniFuel experiment is a passively controlled capsule design deployed in the High Flux Isotope Reactor (HFIR) through which fuel performance data is collected post-irradiation. Separate effects fuels irradiation capabilities are being expanded at ORNL by developing instrumented capsule designs that aim to capture fuel performance phenomena in-situ. One such capsule will specifically target fuel specimen thermal conductivity changes as a function of fuel burnup. Due to the complexity of making this measurement on nuclear fuel in-pile, this paper describes the necessary out-of-pile testing conducted on the thermal conductivity capsule (TCC) design. The measurement is ascertained via a thermopile system with heat transferred unidirectionally through a surrogate fuel specimen sandwiched between two conductive materials. The capsules investigated in this study are representative of the in-pile design, with the primary departure from irradiation conditions being the distribution of heat generation within the capsule. In the out-of-pile experiment, an external heater was used to drive heat through the conductive slug materials and into the specimen. This paper expounds the design of the out-of-pile experimental system and the thermal conductivity measurement technique. Predictive models used to determine the sensitivity of the measurement to variables governing thermal contact conductance between the specimen and slug materials and to predict experimental results are also described. Data from the out-of-pile experiment will be used to validate the readiness of the design for insertion into HFIR for irradiation.

Parker, Trevor [ORNL]↗

The Role of Water Volume Fraction on Water Adsorption in Anion Exchange Membranes

Water absorption plays a key role in the performance of polymeric anion exchange membranes. It influences important properties such as ionic conductivity and mechanical strength and alters their performance as solid electrolytes in hydrogen electrochemical devices for energy conversion. However, computational approaches that address the relationship between the polymer design and the absorption process are scarce. In this work, we introduce a simple thermodynamic model to predict the water absorption isotherms of polyelectrolyte membranes in contact with a water vapor reservoir that incorporates the specific chemical design of the polymers. The model accurately predicts the water content and macrostructural properties of polynorbornene membranes as a function of the water activity and successfully captures the effect of various polymer design parameters. The energy of pairwise attractive interactions predicted by our model provides a means to interpret the absorption process at the molecular level. The model also reveals the most significant favorable and unfavorable contributions to the free energy and indicates that their balance is solely governed by the water volume fraction, regardless of the polymer design. This universal behavior leads to important implications in the search for better ion exchange membranes.

Zaldivar, Gervasio [New York University, Brooklyn,↗

Sensitivity of a Subject-specific Ankle Sprain Simulation to Extrinsic Versus Intrinsic Biomechanical Factors

Ankle sprains are the most common musculoskeletal injury in sport and military activity, despite existing prophylactic strategies. The purpose of this report was to develop a probabilistic simulation of lateral ankle sprains during single-limb drop landing, towards accelerating innovation in ankle sprain prevention. A deterministic, subject-specific musculoskeletal model was extended with automation and probabilistic distributions on sprain-related biomechanical factors. Probabilistic simulations were generated using traditional Monte Carlo techniques and the advanced mean value method, a more computationally-efficient approach. Predicted distributions of peak ankle joint rotations, velocities, and moments borne by supporting passive structures agreed favorably with the deterministic model and with reports of real sprain biomechanics. Parameter sensitivities identified that predictions were most strongly influenced by drop height, subtalar joint posture at contact, invertor/evertor co-activation, and passive ankle stiffness. The advanced mean value method predicted confidence bounds comparable to a 1000-trial Monte Carlo simulation, and required only 14 model evaluations and 4-min processing time. The extended probabilistic simulation may be useful to virtually test new prophylactic strategies for ankle sprains, and is made available for open-source use (https://simtk.org/projects/sprain-sim).

60 APPLIED LIFE SCIENCES↗