Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable threshold”

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 109 records · Page 6

Resistive switching properties of monolayer h-BN atomristors with different electrodes

Resistive switching properties based on molecular beam epitaxy-grown monolayer hexagonal boron nitride (h-BN) atomristors are studied by using metal insulator metal configurations with different electrode materials. Au/monolayer h-BN/Ni devices demonstrate a forming-free bipolar resistive switching (BRS) behavior, a good endurance with up to 97 cycles at a high compliance current of 100 mA, an average on/off ratio of 10 3 , and a low set/reset voltage variability. Metal/monolayer h-BN/graphite/Co devices exhibit self-compliant current BRS characteristics. Both metal/h-BN/Ni and metal/h-BN/graphite/Co devices show the coexistence of BRS, unipolar resistive switching (URS), and nonvolatile threshold switching (TH) modes. The formation of conductive filaments is attributed to the diffusion and trapping of metal ions on the defect sites driven by the electric field, while the rupture is driven by the electric field in BRS and by Joule heating in URS and TH modes.

74 ATOMIC AND MOLECULAR PHYSICS↗

Insights on seasonal solifluction processes in warm permafrost Arctic landscape using a dense monitoring approach across adjacent hillslopes

Solifluction processes in the Arctic are highly complex, introducing uncertainties in estimating current and future soil carbon storage and fluxes, and assessment of hillslope and infrastructure stability. This study aims to enhance our understanding of triggers and drivers of soil movement of permafrost-affected hillslopes in the Arctic. To achieve this, we established an extensive soil deformation and temperature sensor network, covering 48 locations across multiple hillslopes within a 1 km² watershed on the Seward Peninsula, AK. We report depth-resolved measurements down to 1.8 m depth for May to September 2022, a period conducive to soil movement due to deepening thaw layers and frequent rain events. Over this period, surface movements of up to 334 mm were recorded. In general, these movements occur close to the thawing front, and are initiated as thawing reaches depths of 0.4 to 0.75 m. The largest movements were observed at the top of the south-east facing slope, where soil temperatures are cold (mean annual soil temperatures averaging -1.13°C) and slopes are steeper than 15°. Our analysis highlights three primary factors influencing movements: slope angle, soil thermal conditions, and thaw depth. The latter two significantly impact the generation of pore water pressures at the thaw–freeze interface. Specifically, soil thermal conditions govern the liquid water content, while thaw depth influences both the height of the water column and, consequently, the pressure at the thawing front. These factors affect soil properties, such as cohesion and internal friction angle, which are crucial determinants of slope stability. This underscores the significance of a precise understanding of subsurface thermal conditions, including spatial and temporal variability in soil temperature and thaw depth, when assessing and predicting slope instabilities. Based on our observations, we developed a Factor of Safety proxy that consistently falls below the triggering threshold for all probes exhibiting displacements exceeding 50 mm. This study offers novel insights into patterns and triggers of hillslope movements in the Arctic and provides a venue to evaluate their impact on soil redistribution.

54 ENVIRONMENTAL SCIENCES↗

Using Maximum Circular Velocity in Halo Occupation Distribution Models to Predict Galaxy Clustering

Abstract Standard halo occupation distribution (HOD) models predict how galaxies occupy halos based on a single property, the virial mass, M vir . To incorporate galaxy assembly bias, we consider alternative halo mass proxies to use in the HOD. A promising replacement is maximum circular velocity, V max , which measures the halo’s potential well, and successfully predicts galaxy clustering in subhalo abundance matching models. We fit galaxy clustering measurements from SDSS DR7 with two modified HODs besides the standard M vir -based model: one where V max is used as the halo property that determines galaxy occupation, and the other adopting an adjustable combination of M vir and V max . Except for the reduced volume M r < − 20 threshold sample fit with V max alone, neither change improves the model performance in predicting the projected two-point correlation function, w p . We conclude that switching the primary variable to V max does not significantly improve the HOD model’s ability to fit galaxy clustering.

Mezini, Lorena↗

Techno-Economic Evaluation of Electrified Vehicle Options in Drayage Fleets

The electrification of drayage fleets offers potential economic and operational benefits, but the financial viability of electrified vehicles remains sensitive to battery cost, energy price, and fleet usage patterns. While total cost of ownership (TCO) is a useful benchmark, fleet operators and investors are equally concerned with investment performance metrics such as payback period (PB) and Internal Rate of Return (IRR), which better reflect financial risks and investment return timelines. This study develops a unified techno-economic framework that jointly evaluates TCO, PB, and IRR to determine when electrified trucks become cost-effective alternatives to diesel trucks. Building on a previously developed cost modeling tool and using real-world telematics data from a Class 8 drayage fleet at the Port of Savannah, the analysis incorporates projected battery cost trajectories, electricity and diesel price trends, vehicle efficiency improvements, and multiple battery capacities. Parameter ranges reflect widely cited projections and observed drayage-duty-cycle variability. A surrogate-modeling method approximates economic performance across thousands of battery cost–electricity price combinations, enabling high-resolution identification of conditions that achieve TCO parity, acceptable PB thresholds, and target IRR levels. Additionally, the study estimates the evolving share of the fleet that can feasibly electrify over time under multiple economic metrics. This integrated framework offers a novel, data-driven approach to inform risk-aware decision-making for fleet electrification and supports investment planning under evolving cost and operational conditions.

Sun, Ruixiao [ORNL] (ORCID:0000000341768676)↗

A Kinetic Model-Driven Techno-Economic Analysis of Plastic Pyrolysis: Linking Process Dynamics to Economic Viability

This study employs a kinetic model integrated into Aspen Plus to predict pyrolysis product distribution under various conditions. A techno-economic assessment calculated the minimum selling price (MSP) of pyrolysis oil under different operating conditions for the baseline capacity of 100 kta, and across eight processing capacities ranging from 30 to 150 kta. The lowest MSP under the baseline capacity is estimated at $\$$420/ton, which is 33% lower than the 2023 average US crude oil price ($\$$74.6/bbl, equivalent to $\$$634/ton based on the density of pyrolysis oil). Under Monte Carlo simulation, accounting for variability in key economic and technical parameters, the mean MSP is estimated at $\$$1137/ton. The economic viability depends on feedstock price remaining below $\$$320/ton, defining the break-even feedstock price threshold. Sensitivity analysis further identifies capital investment and transportation cost as key economic drivers. Capacities beyond 90 kta show limited economies of scale benefits. Reducing product storage time cuts capital costs by 7% but raises operational risk. Uncertainty analysis suggests the economic feasibility of pyrolysis oil is unlikely to compete with crude oil without policy incentives.

petrochemicals↗

Warming of the lower Columbia River , 1853 to 2018

Abstract Water temperature is a critical ecological indicator; however, few studies have statistically modeled century‐scale trends in riverine or estuarine water temperature, or their cause. Here, we recover, digitize, and analyze archival temperature measurements from the 1850s onward to investigate how and why water temperatures in the lower Columbia River are changing. To infill data gaps and explore changes, we develop regression models of daily historical Columbia River water temperature using time‐lagged river flow and air temperature as the independent variables. Models were developed for three time periods (mid‐19th, mid‐20th, and early 21st century), using archival and modern measurements (1854–1876; 1938–present). Daily and monthly averaged root‐mean‐square errors overall are 0.89°C and 0.77°C, respectively for the 1938–2018 period. Results suggest that annual averaged water temperature increased by 2.2°C ± 0.2°C since the 1850s, a rate of 1.3°C ± 0.1°C/century. Increased water temperatures are seasonally dependent. An increase of approximately 2.0°C ± 0.2°C/century occurs in the July–Dec time‐frame, while springtime trends are statistically insignificant. Rising temperatures change the probability of exceeding ecologically important thresholds; since the 1850s, the number of days with water temperatures over 20°C increased from ~5 to 60 per year, while the number below 2°C decreased from ~10 to 0 days/per year. Overall, the modern system is warmer, but exhibits less temperature variability. The reservoir system reduces sensitivity to short‐term atmospheric forcing. Statistical experiments within our modeling framework suggest that increased water temperature is driven by warming air temperatures (~29%), altered river flow (~14%), and water resources management (~57%).

54 ENVIRONMENTAL SCIENCES↗

Thermal Model for Horizontal In-Water and In-Air Mini-Plate Experiments in the Advanced Test Reactor

Mini-Plate (MP) irradiation tests are a fueled experiment designed for irradiation in multiple test locations in the Advanced Test Reactor (ATR) at the Idaho National Laboratory (INL). MP experiments are considered a non-instrumented drop-in test where small aluminum-clad fuel plate samples are cooled directly by the ATR Primary Coolant System (PCS) water. These fueled experiments contain aluminum-clad fuel mini plates consisting of monolithic U-10Mo. Previous thermal analyses have been documented concerning oxide growth[1] and thermal safety margin[2] for the MP-2 experiment. A thermal analysis has been performed on the MP-2 experiment using a new technique for calculating the necessary time after the experiment has been irradiated and reactor shutdown necessary to not exceed the 95% lower bound blister threshold temperature for the fuel in a horizontal in air configuration. A horizontal in water model has also been analyzed and will be discussed. After the experiment is irradiated in the ATR, it is necessary to calculate the fuel meat temperatures at 12 hours after shutdown. This is done to assure that at the 12 hour mark, if this experiment is removed from the reactor and accidentally dropped onto the working tray and is in the most severe geometrical configuration of horizontal (minimum heat transfer) that fission products will not be released. The blister threshold temperature is considered this safety limit. The horizontal in air analysis is also performed considering that the experiment is horizontal in a drained canal and filled with air and natural convection heat transfer with air on the outside. The initial mini-plate testing campaign known as MP-1 [3] has been irradiated and its purpose was to form the base fuel specimens fabricated by various processes so that a single process can be down-selected for qualification. MP-2 will then serve to populate a comprehensive fuel performance data set for the down-selected fabrication process. The MP-2 experiment will be irradiated in several locations of the ATR as shown in Figure 1. The fuel plate thickness and geometry is identical to the MP-1 experiment [3], with the only variable being the fuel meat thickness. This model of the B-12 position is treated as encompassing the I-21, I-24, and B-9 positions since it has more fuel and higher total heat than the other positions. Since the blister threshold temperature is a function of fission density, a series of fission densities is calculated for the B-12 position to cover the other three positions. The B 12 position specifically uses fuel meat that is 0.016 in thick by 3.25 in long and 0.75 in tall. The total plate thickness is 0.050 in. Aluminum cladding forms the barrier between the PCS and the fuel meat. The purpose of this paper is to demonstrate this new technique developed to calculate the necessary decay time after shutdown for a series of fission densities using the ABAQUS [4] finite element and heat transfer code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Dynamics of streamflow permanence in a headwater network: Insights from catchment-scale model simulations

The hillslope and channel dynamics that govern streamflow permanence in headwater systems have important implications for ecosystem functioning and downstream water quality. Recent advancements in process-based, semi-distributed hydrologic models that build upon empirical studies of streamflow permanence in well-monitored headwater catchments show promise for characterizing the dynamics of streamflow permanence in headwater systems. However, few process-based models consider the continuum of hillslope-stream network connectivity as a control on streamflow permanence in headwater systems. The objective of this study was to expand a process-based, catchment-scale hydrologic model to better understand the spatiotemporal dynamics of headwater streamflow permanence and to identify controls of streamflow expansion and contraction in a headwater network. Further, we aimed to develop an approach that enhanced the fidelity of model simulations, yet required little additional data, with the intent that the model might be later transferred to catchments with limited long-term and spatially explicit measurements. This approach facilitated network-scale estimates of the controls of streamflow expansion and contraction, albeit with higher degrees of uncertainty in individual reaches due to data constraints. Our model simulated that streamflow permanence was highly dynamic in first-order reaches with steep slopes and variable contributing areas. The simulated stream network length ranged from nearly 98±2% of the geomorphic channel extent during wet periods to nearly 50±10% during dry periods. The model identified a discharge threshold of approximately 1 mm d –1 , above which the rate of streamflow expansion decreases by nearly an order of magnitude, indicating a lack of sensitivity of streamflow expansion to hydrologic forcing during high-flow periods. Overall, we demonstrate that process-based, catchment-scale models offer important insights on the controls of streamflow permanence, despite uncertainties and limitations of the model. We encourage researchers to increase data collection efforts and develop benchmarks to better evaluate such models.

54 ENVIRONMENTAL SCIENCES↗

Superscaling variable and neutrino energy reconstruction from theoretical predictions to experimental limitations

We introduce the novel approach of using the superscaling variable as an observable and an analysis tool in the context of charged current neutrino-nucleus interactions. We study the relation between the superscaling variable and the removal energy, in addition to other fundamental parameters of the neutrino-nucleus interaction models. In the second half of the paper, we discuss the experimental viability of this measurement following a study of neutrino energy and missing momentum reconstruction. We show that the superscaling variable is measurable in neutrino interaction experiments provided that the proton is detected in the final state. We discuss the resolution of this measurement, and the limitation imposed by the proton’s detection threshold. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Nuclear Thermal Energy Storage Configurations for Industrial Combined Heat and Power Supply: Conceptual Study and Engineering Designs

The industries examined in this report primarily rely on moderate-temperature heat provided by gas- or coal-fired boilers and combined heat and power (CHP) plants, delivered through standard process steam systems. High-temperature energy demands are often industry-specific and typically exceed the capabilities of high-temperature gas-cooled reactors (HTGRs). While it is technically feasible to replace process steam from fossil-based heat sources with nuclear energy, certain industries, such as methanol production and pulp and paper, face technoeconomic challenges in integrating nuclear energy without major changes or a technological shift. This is mainly due to the limited external energy demand remaining after the use of internal byproducts, waste heat recovery, and simple efficiency improvements. Achieving full decarbonization of these processes with nuclear energy would require significant technological advancements, involving experimental technology and substantial investments, making widespread adoption in existing industrial plants unlikely in the near term. This study reviews TES options in the context of enabling a flexible CHP supply while maintaining a steady nuclear heat input. Heat storage systems that interface between the reactor primary fluid and the CHP system offer superior performance and flexibility. Specifically, steam extraction downstream of the reheater with a two-tank molten-salt TES appears as the best solution regarding thermodynamic system benefits and system drawbacks. Using selected system configurations, a conceptual design of an industrial energy park was developed for industries with varying energy demands, such as steel production plants utilizing electric arc furnaces (EAFs) and chemical plants, as well as for those with constant energy demands, like petroleum refineries. This design highlights the capabilities of TES and explores its potential business cases. The study also conceptually develops the potential for integrating additional energy sources with nuclear systems through the implementation of TES. The potential of the HTGR-TES-CHP system was also evaluated considering key uncertainties such as industrial demand profiles, external grid access availability, and eligible tax credit levels, using the Holistic Energy Resource Optimization Network. Sensitivity of net present value to these uncertainties was analyzed to determine the optimal number of nuclear reactors (and CHP systems) and the suitable TES capacity. The results were interpreted from a decision-maker’s perspective, focusing on three key areas: deployment strategy (oversized units vs. undersized units with TES support), industrial process characteristics (thermal-intensive single profiles vs. electricity-intensive combined profiles), and operational goals (maximizing profits vs. minimizing natural gas (NG) consumption or external grid dependence). The optimization results indicate that the HTGR-TES-CHP system significantly reduces reliance on NG boilers for individual industrial processes by 9-60% (in NG capacity factor), with an average reduction of 38%, compared to standalone NG boiler operation case (Business As Usual [BAU]). For combined industrial processes, the reduction ranges from 37-77%, with an average of 60%. Additionally, the system greatly reduces dependence on external grids. In meeting industrial electrical demands, a 33-100% self-sufficient internal electricity supply is achieved for single industrial process, with an average of 74%, compared to the BAU scenario, where 100% of electricity is imported. For combined processes, 35-100% of internal electricity demands are met by the reactor, with an average of 73%. At last, the relative NG price levels at which the proposed HTGR-TES-CHP system can cost-effectively enter the market currently dominated by existing NG boilers were estimated. For a moderate HTGR CAPEX level ($\$$2500/kWth, $\$$6329/kWe), the analysis suggests that NG prices must be 2.5 to 7 times higher than HTGR variable operating and maintenance costs for single industrial process, and 5.5 to 9.5 times higher for a combined process scenario. Tax credit modeling shows that the Investment Tax Credit significantly reduces the price threshold needed to break even, making the system competitive with NG boilers in certain cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Automated identification of local contamination in remote atmospheric composition time series

Abstract. Atmospheric observations in remote locations offer a possibility of exploring trace gas and particle concentrations in pristine environments. However, data from remote areas are often contaminated by pollution from local sources. Detecting this contamination is thus a central and frequently encountered issue. Consequently, many different methods exist today to identify local contamination in atmospheric composition measurement time series, but no single method has been widely accepted. In this study, we present a new method to identify primary pollution in remote atmospheric datasets, e.g., from ship campaigns or stations with a low background signal compared to the contaminated signal. The pollution detection algorithm (PDA) identifies and flags periods of polluted data in five steps. The first and most important step identifies polluted periods based on the derivative (time derivative) of a concentration over time. If this derivative exceeds a given threshold, data are flagged as polluted. Further pollution identification steps are a simple concentration threshold filter, a neighboring points filter (optional), a median, and a sparse data filter (optional). The PDA only relies on the target dataset itself and is independent of ancillary datasets such as meteorological variables. All parameters of each step are adjustable so that the PDA can be “tuned” to be more or less stringent (e.g., flag more or fewer data points as contaminated). The PDA was developed and tested with a particle number concentration dataset collected during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic. Using strict settings, we identified 62 % of the data as influenced by local contamination. Using a second independent particle number concentration dataset also collected during MOSAiC, we evaluated the performance of the PDA against the same dataset cleaned by visual inspection. The two methods agreed in 94 % of the cases. Additionally, the PDA was successfully applied to a trace gas dataset (CO2), also collected during MOSAiC, and to another particle number concentration dataset, collected at the high-altitude background station Jungfraujoch, Switzerland. Thus, the PDA proves to be a useful and flexible tool to identify periods affected by local contamination in atmospheric composition datasets without the need for ancillary measurements. It is best applied to data representing primary pollution. The user-friendly and open-access code enables reproducible application to a wide suite of different datasets. It is available at https://doi.org/10.5281/zenodo.5761101 (Beck et al., 2021).

54 ENVIRONMENTAL SCIENCES↗

Machine-learning identification of the variability of mean velocity and turbulence intensity for wakes generated by onshore wind turbines: Cluster analysis of wind LiDAR measurements

Light detection and ranging (LiDAR) measurements of isolated wakes generated by wind turbines installed at an onshore wind farm are leveraged to characterize the variability of the wake mean velocity and turbulence intensity during typical operations, which encompass a breadth of atmospheric stability regimes and rotor thrust coefficients. The LiDAR measurements are clustered through the k-means algorithm, which enables identifying the most representative realizations of wind turbine wakes while avoiding the imposition of thresholds for the various wind and turbine parameters. Considering the large number of LiDAR samples collected to probe the wake velocity field, the dimensionality of the experimental dataset is reduced by projecting the LiDAR data on an intelligently truncated basis obtained with the proper orthogonal decomposition (POD). The coefficients of only five physics-informed POD modes are then injected in the k-means algorithm for clustering the LiDAR dataset. The analysis of the clustered LiDAR data and the associated supervisory control and data acquisition and meteorological data enables the study of the variability of the wake velocity deficit, wake extent, and wake-added turbulence intensity for different thrust coefficients of the turbine rotor and regimes of atmospheric stability. Furthermore, the cluster analysis of the LiDAR data allows for the identification of systematic off-design operations with a certain yaw misalignment of the turbine rotor with the mean wind direction.

17 WIND ENERGY↗

Variable voltage electron microscopy: Toward atom-by-atom fabrication in 2D materials

The precisely focused electron beam (e-beam) in scanning transmission electron microscopy has been found to be a versatile tool for patterning matter on the atomic level and enabling atom-by-atom fabrication, however, beam-induced phenomena are known to be extremely sensitive to the e-beam energy and are closely correlated with the proximity of the knock-on threshold. Here we provide a method to control the energy transferred to the sample while facilitating much faster changes in accelerating voltages by maintaining a nearly constant temperature in the electron microscope lenses. We use this method to demonstrate in-situ nano-milling of a graphene film followed rapidly by “gentler” imaging at lower energy. Additionally, the insertion and controlled movement of silicon dopants in graphene is demonstrated by employing a different e-beam energy during each process. We believe the incorporation of variable e-beam energy will broaden the potential for atomic scale e-beam fabrication.

36 MATERIALS SCIENCE↗

Knowledge spillovers between PV installers can reduce the cost of installing solar PV

We analyze pricing in PV systems using data from 2008-2014 to identify the effects of knowledge spillovers in reducing the installed cost of PV. This paper estimates the size of the effects of learning by doing and knowledge spillovers using multiple formulations of spillover related variables. We found knowledge spillovers between firms within a county to be a significant and substantial factor in reducing the costs of PV. However, these spillovers reduce costs only for firms over a certain size threshold, and no cost-reducing spillovers were found for smaller installers. Geographic spillovers, within a firm from one county to another were also significant although not as large as the local between-firm effects. We ran 43 specifications on the data and generally found these main results to be robust, although not in every specification. One implication of these results is that policies that subsidize demand for PV are leading to the creation of new knowledge that would not exist at lower levels of demand. Furthermore, that the spillovers become stronger at higher levels of experience suggests that these subsidies would need to be substantial, particularly in new markets with many small firms.

14 SOLAR ENERGY↗

Shock wave formation in radiative plasmas

The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thick and thin regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thick and thin solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive testcase has been designed to benchmark radiative hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics-based Induced Earthquake Forecasting: Process Understanding, and Hazards Mitigation

Disposal of saltwater co-produced with oil and gas is linked to elevated seismicity in the Central and Midwest US. There is a concern that these events may lead to widespread damage and an overall increase in seismicity. Thus an improved understanding of the spatially and temporally variable deformation and stress field associated with fluid injection operation is critically necessary for evaluating time-varying seismic hazards. Despite the improvements in seismic monitoring capacity and the resulting decrease in the magnitude detection threshold, estimates of induced earthquake probability remain elusive due to insufficient models incapable of accounting for the complex physics governing the process of induced seismicity. The proposed research effort will comprehensively analyze, integrate, and interpret geodetic, injection, and seismic data in the vicinity of the injection sites in Oklahoma to resolve the 4-dimensional distribution of pore pressure and stress in the shallow crust. This project, in particular, is focused on exploring the statistical relation between injection operation and increased earthquake hazard. The amplitude of and the extent to which pore pressure changes are determined by some factors, in particular, the hydrogeological properties of the rocks, such as diffusivity. Thus the available deformation data will be used to constrain hydrogeological properties of the medium, to accurately resolve the evolution of crustal stresses due to fluid injection. Having the time-varying models of stress changes, a statistical framework will be implemented to estimate the time-dependent probability of large earthquakes on the nearby fault systems. These data and models help to improve seismic hazard estimates and aid in constructing operational-induced earthquake forecast models. This information can also be integrated into the updated U.S. National Seismic Hazard Map, which local communities and authorities use in their earthquake risk estimates and mitigation efforts.

58 GEOSCIENCES↗

Understanding Microstructure Variability in Vapor-Deposited Energetic Materials by Using Phase-Field Methods

Critical components, such as detonators, in Sandia's stockpile contain heterogeneous materials whose performance and reliability depend on accurate, predictive models of coupled, complex phenomena to predict their synthesis, processing, and operation. Ongoing research in energetic materials has shown that microstructural properties, such as density, pore-size, morphology, and specific surface area are correlated to their initiation threshold and detonation behavior. However, experiments to study these specific characteristics of energetic materials are challenging and time consuming. Therefore, in this work, we turn to mesoscale modeling methods that may be capable of reproducing some observed phenomena to refine and predict outcomes beforehand. Even so, we have no physics-based modeling capability to predict how the microstructure of an energetic material will evolve over near- and long-term time scales. Thus, the goal of this work is to (i) identify any knowledge gaps in how the underlying microstructure forms and evolves during the synthesis process, and (ii) develop and test a mesoscale phase-field model for vapor deposition to capture critical mechanisms of microstructure formation, evolution, and variability in vapor-deposited energetic materials, such as processing conditions, material properties, and substrate interactions. Based on this work, the phase-field method is shown to be a valuable tool for developing the necessary models containing coupled, complex phenomena to investigate and understand the synthesis and processing of energetic materials.

36 MATERIALS SCIENCE↗

A systematic strategy for estimating hERG block potency and its implications in a new cardiac safety paradigm

Introduction: hERG block potency is widely used to calculate a drug's safety margin against its torsadogenic potential. Previous studies are confounded by use of different patch clamp electrophysiology protocols and a lack of statistical quantification of experimental variability. Since the new cardiac safety paradigm being discussed by the International Council for Harmonisation promotes a tighter integration of nonclinical and clinical data for torsadogenic risk assessment, a more systematic approach to estimate the hERG block potency and safety margin is needed. Methods: A cross-industry study was performed to collect hERG data on 28 drugs with known torsadogenic risk using a standardized experimental protocol. A Bayesian hierarchical modeling (BHM) approach was used to assess the hERG block potency of these drugs by quantifying both the inter-site and intra-site variability. A modeling and simulation study was also done to evaluate protocol-dependent changes in hERG potency estimates. Results: A systematic approach to estimate hERG block potency is established. The impact of choosing a safety margin threshold on torsadogenic risk evaluation is explored based on the posterior distributions of hERG potency estimated by this method. The modeling and simulation results suggest any potency estimate is specific to the protocol used. Discussion: This methodology can estimate hERG block potency specific to a given voltage protocol. The relationship between safety margin thresholds and torsadogenic risk predictivity suggests the threshold should be tailored to each specific context of use, and safety margin evaluation may need to be integrated with other information to form a more comprehensive risk assessment.

60 APPLIED LIFE SCIENCES↗