Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Simulation”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

SOLEDGE3X full vessel plasma boundary simulations of ITER non-active phase plasmas

The onset of detachment in the ITER machine is analyzed in this work through the help of 2D-axisymmetric boundary plasma simulations with the SOLEDGE3X-EIRENE code, which features a numerical domain for the plasma solver extending up to the first wall. The plasma boundary is computed in scenarios from the first non-active phase of ITER, in pure H and at 20 MW. This set of simulations is used in two aspects: first, to study the plasma detachment in the divertor, and second, the plasma conditions, fluxes, and beryllium erosion at the first wall. Here, the code results are also compared to those obtained with the well-established SOLPS-ITER code, which includes a plasma numerical domain only covering the main SOL. Results show an increase in the SOL width λ q with increasing density, and a detailed analysis is carried out, for the first time, on each of the different plasma-neutral interactions in the code’s physics model in EIRENE. The gross beryllium erosion rates of first wall panels are estimated from 2D simulations, with the aim of assessing their sensitivity to two parameters: the divertor density regime, and the presence of density shoulders in the far-SOL formed by enhanced perpendicular transport at this location. The erosion contributions from neutrals and ions are considered in each case, and the charge-exchange atoms fluxes and energy distributions are provided, highlighting the two atom populations (cold and charge-exchange).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Observing system simulation experiments related to space-borne LIDAR wind profiling. I - Forecast impacts of highly idealized observing systems. II - Sensitivity to atmospheric and instrumental influences

It is pointed out that extended-range numerical weather prediction requires complete global and accurate three-dimensional fields of temperature, moisture, and winds as a function of pressure throughout the troposphere, lower stratosphere, and at the earth's surface. Recent developments in the use of incoherent optical radar (LIDAR) systems indicate the feasibility of obtaining pressure as a function of height with presumably very high vertical resolution and accuracy. In evaluating the feasibility of an employment of spaceborne LIDAR systems as an aid for obtaining data required for weather prediction, simulation studies could provide important information. The present investigation is concerned with the conduction of such simulation experiments. It is found that the results of a composite system of temperature and pressure show significant improvements over those of temperature only. The effects of simulating several sources of LIDAR wind error for two different representations of nature are also examined.

Dlouhy, R.↗
NASA NTRS1999

Climate Change and the Brewer-Dobson Circulation

Using GCM (General Circulation Model) simulations of the troposphere-stratosphere circulation in different climate regimes, the sensitivity of the Brewer-Dobson Circulation to climate change will be discussed. Included will be simulations of the stratosphere with increasing/doubled CO2 as done with different versions of the GISS (Goddard Institute for Space Studies) Global Climate Middle Atmosphere Model, with and without ozone response. To put the results in a broader perspective, stratospheric circulation changes from several paleoclimate simulations will also be presented, including the Last Glacial Maximum and the Paleocene, each with several different (plausible) latitudinal sea surface temperature gradients. Hence the climate changes investigated will include strong variations in baroclinic and topographic tropospheric planetary wave and gravity wave forcing, as well as global mean temperature and atmospheric CO2 levels. Results will be analyzed in terms of wave driving via E-P flux convergences and gravity wave effects, and will be shown to be strongly dependent on wave propagation characteristics.

Rind, David↗
NASA NTRS1999

Climate Change and the Brewer-Dobson Circulation

Using GCM simulations of the troposphere-stratosphere circulation in different climate regimes, the sensitivity of the Brewer-Dobson Circulation to climate change will be discussed. Included will be simulations of the stratosphere with increasing/doubled C02 as done with different versions of the GISS Global Climate Middle Atmosphere Model, with and without ozone response. To put the results in a broader perspective, stratospheric circulation changes from several paleoclimate simulations will also be presented, including the Last Glacial Maximum and the Paleocene, each with several different (plausible) latitudinal sea surface temperature gradients. Hence the climate changes investigated will include strong variations in baroclinic and topographic tropospheric planetary wave and gravity wave forcing, as well as global mean temperature and atmospheric C02 levels. Results will be analyzed in terms of wave driving via E-P flux convergences and gravity wave effects, and will be shown to be strongly dependent on wave propagation characteristics.

Rind, David↗
NASA NTRS1992

Flight simulator fidelity assessment in a rotorcraft lateral translation maneuver

A model-based methodology for assessing flight simulator fidelity in closed-loop fashion is exercised in analyzing a rotorcraft low-altitude maneuver for which flight test and simulation results were available. The addition of a handling qualities sensitivity function to a previously developed model-based assessment criteria allows an analytical comparison of both performance and handling qualities between simulation and flight test. Model predictions regarding the existence of simulator fidelity problems are corroborated by experiment. The modeling approach is used to assess analytically the effects of modifying simulator characteristics on simulator fidelity.

Hess, R. A.↗
NASA NTRS1995

A mesoscale modeling study of the atmospheric circulation of high southern latitudes

The meteorology of high southern latitudes during winter is simulated using a cloud-free version of The Pennsylvania State University-National Center for Atmospheric Research Mesoscale Model version 4 (MM4) with a 100-km horizontal resolution. Comparisons between idealized simulations of Antarctic with MM4 and with the mesoscale model of Parish and Waight reveal that both models produce similarly realistic velocity fields in the boundary layer. The latter model tends to produce slightly faster drainage winds over East Antarctica. The intensity of the katabatic winds produced by MM4 is sensitive to parameterizations of boundary layer fluxes. Two simulations performed with MM4 using analyses from the European Center for Medium-Range Weather Forecasts (ECMWF) for June 1988 as initial and boundary conditions. A simulation of the period from 000 UTC 2 June to 0000 UTC 8 June produces realistic synoptic phenomena including ridge development over East Antarctica, frontogenesis over the Amundsen Sea, and a katabatic surge over the Ross Ice Shelf. The simulated time-averaged fields for June 1988, particularly that of a 500-hPa height, are in good agreement with time-averaged fields analyzed by the ECMWF. The results of the simulations provide detailed features of the Antarctic winter boundary layer along the steeply sloping terrain. Highest boundary layer wind speeds averaged over the month-long simulation are approximately 20 m/s. The lack of latent heating in the simulations apparently results in some bias in the results. In particular, the cloud-free version of MM4 underpredicts the intensity of lows in the sea level pressure field.

Hines, K. M.↗
NASA NTRS2015

An Assessment of Antarctic Circumpolar Current and Southern Ocean Meridional Overturning Circulation During 1958-2007 in a Suite of Interannual CORE-2 Simulations

In the framework of the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II), we present an analysis of the representation of the Antarctic Circumpolar Current (ACC) and Southern Ocean Meridional Overturning Circulation (MOC) in a suite of seventeen global ocean-sea ice models. We focus on the mean, variability and trends of both the ACC and MOC over the 19582007 period, and discuss their relationship with the surface forcing. We aim to quantify the degree of eddy saturation and eddy compensation in the models participating in CORE-II, and compare our results with available observations, previous fine-resolution numerical studies and theoretical constraints. Most models show weak ACC transport sensitivity to changes in forcing during the past five decades, and they can be considered to be in an eddy saturated regime. Larger contrasts arise when considering MOC trends, with a majority of models exhibiting significant strengthening of the MOC during the late 20th and early 21st century. Only a few models show a relatively small sensitivity to forcing changes, responding with an intensified eddy-induced circulation that provides some degree of eddy compensation, while still showing considerable decadal trends. Both ACC and MOC interannual variability are largely controlled by the Southern Annular Mode (SAM). Based on these results, models are clustered into two groups. Models with constant or two-dimensional (horizontal) specification of the eddy-induced advection coefficient show larger ocean interior decadal trends, larger ACC transport decadal trends and no eddy compensation in the MOC. Eddy-permitting models or models with a three-dimensional time varying show smaller changes in isopycnal slopes and associated ACC trends, and partial eddy compensation. As previously argued, a constant in time or space is responsible for a poor representation of mesoscale eddy effects and cannot properly simulate the sensitivity of the ACC and MOC to changing surface forcing. Evidence is given for a larger sensitivity of the MOC as compared to the ACC transport, even when approaching eddy saturation. Future process studies designed for disentangling the role of momentum and buoyancy forcing in driving the ACC and MOC are proposed.

Southern Ocean dynamics↗
DOE OSTI2024

Projecting Changes in the Frequency and Magnitude of Ozone Pollution Events Under Uncertain Climate Sensitivity

Abstract Climate change is projected to worsen ozone pollution over many populated regions, with larger impacts at higher concentrations. More intense and frequent ozone episodes risk setbacks to human health and environmental policy achievements. However, assessing these changes is complicated by uncertain climate sensitivity, closely related to climate model response, and internal variability in simulations projecting climate's influence on air quality. Here, leveraging a global modeling framework that one‐way couples a human activity model, an Earth system model of intermediate complexity, and an atmospheric chemistry model, we investigate the role of climate sensitivity in climate‐induced changes to high ozone pollution episodes in the United States using multiple greenhouse gas emissions scenarios, representations of climate sensitivity, and initial condition members. We bias correct and evaluate historical model simulations, identifying modeled and observed O 3 episodes using extreme value theory, and extend the approach to projections of mid‐ and end‐century climate impacts. Results show that the influence of climate sensitivity can be as significant as that of greenhouse gas emissions scenario absent precursor emissions changes. Climate change is projected to increase the magnitude of the highest annually occurring O 3 concentrations by over 2.3 ppb on average across the U.S. at mid‐century under a high climate sensitivity and moderate emissions scenario, but the increase is limited to less than 0.3 ppb under lower climate sensitivity. Further, we show that areas in the U.S. currently meeting air quality standards risk being pushed into non‐compliance due to a climate‐induced increase in frequency of high ozone days.

Environmental Sciences & Ecology↗
DOE OSTI2022

Comparison of Electron Capture Rates in the N = 50 Region using 1D Simulations of Core-collapse Supernovae

Abstract Recent studies have highlighted the sensitivity of core-collapse supernovae (CCSNe) models to electron-capture (EC) rates on neutron-rich nuclei near the N = 50 closed-shell region. In this work, we perform a large suite of one-dimensional CCSN simulations for 200 stellar progenitors using recently updated EC rates in this region. For comparison, we repeat the simulations using two previous implementations of EC rates: a microphysical library with parametrized N = 50 rates (LMP), and an older independent-particle approximation (IPA). We follow the simulations through shock revival up to several seconds post-bounce, and show that the EC rates produce a consistent imprint on CCSN properties, often surpassing the role of the progenitor itself. Notable impacts include the timescale of core collapse, the electron fraction and mass of the inner core at bounce, the accretion rate through the shock, the success or failure of revival, and the properties of the central compact remnant. We also compare the observable neutrino signal of the neutronization burst in a DUNE-like detector, and find consistent impacts on the counts and mean energies. Overall, the updated rates result in properties that are intermediate between LMP and IPA, and yet slightly more favorable to explosion than both.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Periodicity significance testing with null-signal templates: reassessment of PTF’s SMBH binary candidates

Periodograms are widely employed for identifying periodicity in time series data, yet they often struggle to accurately quantify the statistical significance of detected periodic signals when the data complexity precludes reliable simulations. We develop a data-driven approach to address this challenge by introducing a null-signal template (NST). The NST is created by carefully randomizing the period of each cycle in the periodogram template, rendering it non-periodic. It has the same frequentist properties as a periodic signal template, and we show with simulations that the distribution of false positives is the same as with the original periodic template, regardless of the underlying data. Thus, performing a periodicity search with the NST acts as an effective simulation of the null (no-signal) hypothesis, without having to simulate the noise properties of the data. We apply the NST method to the supermassive black hole binaries (SMBHB) search in the Palomar Transient Factory (PTF), where Charisi et al. had previously proposed 33 high signal-to-noise candidates utilizing simulations to quantify their significance. Our approach reveals that these simulations do not capture the complexity of the real data. There are no statistically significant periodic signal detections above the non-periodic background. To improve the search sensitivity, we introduce a Gaussian quadrature based algorithm for the Bayes Factor with correlated noise as a test statistic. We show with simulations that this improves sensitivity to true signals by more than an order of magnitude. However, the Bayes Factor approach also results in no statistically significant detections in the PTF data.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2018

The PMIP4 Contribution to CMIP6 - Part 1: Overview and Over-Arching Analysis Plan

This paper is the first of a series of four GMD (Geoscientific Model Development) papers on the PMIP4-CMIP6 (Paleoclimate Modelling Intercomparison Project - Phase 4 -- Coupled Model Intercomparison Project - Phase 6) experiments. Part 2 (Otto-Bliesner et al., 2017) gives details about the two PMIP4-CMIP6 interglacial experiments, Part 3 (Jungclaus et al., 2017) about the last millennium experiment, and Part 4 (Kageyama et al., 2017) about the Last Glacial Maximum experiment. The mid-Pliocene Warm Period experiment is part of the Pliocene Model Intercomparison Project (PlioMIP) - Phase 2, detailed in Haywood et al. (2016). The goal of the Paleoclimate Modelling Intercomparison Project (PMIP) is to understand the response of the climate system to different climate forcings for documented climatic states very different from the present and historical climates. Through comparison with observations of the environmental impact of these climate changes, or with climate reconstructions based on physical, chemical, or biological records, PMIP also addresses the issue of how well state-of-the-art numerical models simulate climate change. Climate models are usually developed using the present and historical climates as references, but climate projections show that future climates will lie well outside these conditions. Palaeoclimates very different from these reference states therefore provide stringent tests for state-of-the-art models and a way to assess whether their sensitivity to forcings is compatible with palaeoclimatic evidence. Simulations of five different periods have been designed to address the objectives of the sixth phase of the Coupled Model Intercomparison Project (CMIP6): the millennium prior to the industrial epoch (CMIP6 name: past1000); the mid-Holocene, 6000 years ago (midHolocene); the Last Glacial Maximum, 21,000 years ago (lgm); the Last Interglacial, 127,000 years ago (lig127k); and the mid-Pliocene Warm Period, 3.2 million years ago (midPliocene-eoi400). These climatic periods are well documented by palaeoclimatic and palaeoenvironmental records, with climate and environmental changes relevant for the study and projection of future climate changes. This paper describes the motivation for the choice of these periods and the design of the numerical experiments and database requests, with a focus on their novel features compared to the experiments performed in previous phases of PMIP and CMIP. It also outlines the analysis plan that takes advantage of the comparisons of the results across periods and across CMIP6 in collaboration with other MIPs.

climate↗
DOE OSTI2020

Critical nuclei at hetero-phase interfaces

Two-step nucleation, in which a metastable intermediate phase acts as a precursor for nucleating a thermodynamically stable phase, has been widely observed in many materials systems and solid-state reactions. Among the advantages of two-step nucleation is that the stable phase may nucleate heterogeneously at the hetero-phase interface between the original and the precursory phases. Although heterogeneous nucleation (HN) theories for homo-phase grain boundaries and inert surfaces are well established, our understanding of HN at reactive hetero-phase interfaces remains incomplete. This deficiency stems from the discontinuity of the chemical potential driving force across the hetero-phase interface, which profoundly affects the fundamental properties of the nucleus in a way that is not properly accounted for in existing models. In this work, we incorporate these effects to extend the classical nucleation theory to HN at hetero-phase interfaces. Our extended model demonstrates that the nucleus shape along the minimum energy path is strongly size-dependent, and this additional degree of freedom can result in the reduction of the critical nucleus volume and associated activation energy barrier by orders of magnitude relative to conventional predictions. The simulation results are used to construct a sensitivity map in the parameter space of interfacial energy and bulk driving force ratios, which quantifies the difference in nucleation barriers predicted by different models.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge

In this work, we present the results of two-dimensional modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge. A multiscale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC) based on PASSKEy discharge modeling package and compressible multi-component reactive flow solver ASURF+ is developed and validated. This model is applied to simulate the impact of non-equilibrium plasma excitation and electrode geometry and heat loss on the dynamics of the discharge from streamer to spark and ignition kernel development in a H 2 /air mixture with a pair of cylindrical electrodes. The results show that the plasmagenerated species (N 2 (A), N 2 (B), N 2 (a'), N 2 (C), O( 1 D), O and H) in the spark and afterglow significantly accelerate the ignition kernel development. The increase of discharge voltage at the same total discharge energy promotes the non-equilibrium active species production. It is found that the production of electronically excited species at higher reduced electric field strength is more efficient in enhancing ignition in comparison to the vibrational excitation and heating. Moreover, the 2D simulation clearly reveals that the electric field and active species distribution are highly non-uniform. The streamers are initiated at the sharp outer edges of the negative and positive electrodes by a strong electric field while the electric field is much weaker at the centerline of the electrodes. Furthermore, the simulations reveal that the ignition enhancement is sensitive to the variation of electrode shape, diameter, and gap size due to the changes of electric field distribution and location of streamer formation. A cylindrical electrode produces a larger discharge volume and ignition kernel than the parabolic and spherical electrodes, when the discharge is localized near the axis of the gap. It is found that there is a non-monotonic dependence of ignition kernel size on the electrode diameter and inter-electrode distance. The increase of electrode diameter and gap size above the optimal conditions leads to the reduction of ignition kernel volume, due to the decrease of active species concentration and gas temperature. At a larger electrode surface area and electrode diameter as well as smaller electrode gap size, the heat loss to electrode plays a greater role in reducing the ignition kernel size and slowing ignition kernel development. This work provides insights and guidance to understand the kinetic enhancement of non-equilibrium plasma and the effects of electrode geometries on ignition for the optimization ignitors in advanced engines.

42 ENGINEERING↗
DOE OSTI2020

Designing a drone delivery network with automated battery swapping machines

Drones are projected to alter last-mile delivery, but their short travel range is a concern. In this study, we propose a drone delivery network design using automated battery swapping machines (ABSMs) to extend ranges. The design minimizes the long-term delivery costs, including ABSM investment, drone ownership, and cost of the delivery time, and locates ABSMs to serve a set of customers. We build a mixed-integer nonlinear program that captures the nonlinear waiting time of drones at ABSMs. To solve the problem, we create an exact solution algorithm that finds the globally optimal solution using a derivative-supported cutting-plane method. To validate the applicability of our program, we conduct a case study on the Chicago Metropolitan area using cost data from leading ABSM manufacturer and geographical data from the planning and operations language for agent-based regional integrated simulation (more commonly known as POLARIS). A sensitivity analysis identifies that ABSM service times and costs are the key parameters impacting the long-term adoption of drone delivery.

25 ENERGY STORAGE↗
DOE OSTI2022

An electrochemical mesoscale tool for modeling the corrosion of structural alloys by molten salt

Understanding the impact of microstructure on corrosion rates can aid the development of corrosion-resistant alloys for molten salt reactors. Here in this work, we develop an electrochemical phase-field model for capturing the microstructure-dependent corrosion of structural alloys by molten salts. As a demonstration problem, we apply this model to capture the selective depletion of Cr from Ni-Cr grain boundaries during corrosion in molten FLiBe salt. We perform sensitivity analysis and model verification on 1D simulations to confirm that the model predicts diffusion-limited kinetics. The model is validated using 1D, 2D, and 3D simulations against experimental data for Ni-5Cr and Ni-20Cr corrosion in molten FLiBe. The 1D simulations predict the corrosion behavior with reasonable accuracy when using an effective diffusion coefficient that accurately represents the grain boundary diffusion. 2D simulations that represent the grain structure underpredict the corrosion. 3D simulations that represent the grain structure predict the corrosion with reasonable accuracy. The corrosion rate predicted by the 3D simulations is proportional to the average grain size at the alloy/salt interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Snapshot multispectral imaging using a diffractive optical network

Abstract Multispectral imaging has been used for numerous applications in e.g., environmental monitoring, aerospace, defense, and biomedicine. Here, we present a diffractive optical network-based multispectral imaging system trained using deep learning to create a virtual spectral filter array at the output image field-of-view. This diffractive multispectral imager performs spatially-coherent imaging over a large spectrum, and at the same time, routes a pre-determined set of spectral channels onto an array of pixels at the output plane, converting a monochrome focal-plane array or image sensor into a multispectral imaging device without any spectral filters or image recovery algorithms. Furthermore, the spectral responsivity of this diffractive multispectral imager is not sensitive to input polarization states. Through numerical simulations, we present different diffractive network designs that achieve snapshot multispectral imaging with 4, 9 and 16 unique spectral bands within the visible spectrum, based on passive spatially-structured diffractive surfaces, with a compact design that axially spans ~72 λ m , where λ m is the mean wavelength of the spectral band of interest. Moreover, we experimentally demonstrate a diffractive multispectral imager based on a 3D-printed diffractive network that creates at its output image plane a spatially repeating virtual spectral filter array with 2 × 2 = 4 unique bands at terahertz spectrum. Due to their compact form factor and computation-free, power-efficient and polarization-insensitive forward operation, diffractive multispectral imagers can be transformative for various imaging and sensing applications and be used at different parts of the electromagnetic spectrum where high-density and wide-area multispectral pixel arrays are not widely available.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Measurement-induced entanglement and teleportation on a noisy quantum processor

Measurement has a special role in quantum theory: by collapsing the wavefunction, it can enable phenomena such as teleportation and thereby alter the ‘arrow of time’ that constrains unitary evolution. When integrated in many-body dynamics, measurements can lead to emergent patterns of quantum information in space–time that go beyond the established paradigms for characterizing phases, either in or out of equilibrium. For present-day noisy intermediate-scale quantum (NISQ) processors, the experimental realization of such physics can be problematic because of hardware limitations and the stochastic nature of quantum measurement. Here we address these experimental challenges and study measurement-induced quantum information phases on up to 70 superconducting qubits. By leveraging the interchangeability of space and time, we use a duality mapping to avoid mid-circuit measurement and access different manifestations of the underlying phases, from entanglement scaling to measurement-induced teleportation. We obtain finite-sized signatures of a phase transition with a decoding protocol that correlates the experimental measurement with classical simulation data. The phases display remarkably different sensitivity to noise, and we use this disparity to turn an inherent hardware limitation into a useful diagnostic. Our work demonstrates an approach to realizing measurement-induced physics at scales that are at the limits of current NISQ processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The charge reduction rate for multiply charged polymer ions via ion–ion recombination at atmospheric pressure

The charge reduction of multiply charged macromolecular ions via recombination with small ions in the gas phase is commonly employed to modulate the charge on macromolecules prior to mass spectrometric and mobility analyses. We employ a recently developed continuum-Molecular Dynamics (MD) calculation approach to determine the recombination rate coefficient of multiply charged (1 to 7 excess positive charged) polyethylene glycol ions (mass of 4600 Da) with smaller singly charged anions, modeled as NO 2 – ions. The continuum-MD approach accounts explicitly for the influence of the background gas on the recombination process, accounts explicitly for ion translational, vibrational, and rotational motion, and enables recombination rate coefficient calculation in nitrogen near atmospheric pressure, wherein neither low pressure nor high pressure recombination theories are strictly applicable. Continuum-MD simulations yield recombination rate coefficients near 3.9 × 10 –14 m 3 s –1 for singly charged ions, increasing to 3.0 × 10 –11 m 3 s –1 for the +7 ion. Pre-existing collision rate coefficient expressions for rigid ions are found to be within a factor of 2–5 of calculations for all charge states, but their use requires knowledge of an appropriate collision distance, which is not well-defined for flexible polymer ions. Continuum-MD-inferred rate coefficients are incorporated into a model of charge reduction, and the charge state distribution versus anion concentration determined with it is compared to charge reduction measurements made via atmospheric pressure differential mobility analysis. Furthermore, good agreement is observed between simulations and experiments; although results are extremely sensitive to the recombination rate coefficients, experimental results are bound by models utilizing rates within a factor of 2 (0.5–2.0×) of the continuum-MD rates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗