Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum scale modeling”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Climate Scenarios for the NASA / USAID SERVIR Project: Challenges for Multiple Planning Horizons

SERVIR, an acronym meaning "to serve" in Spanish, is a joint venture between NASA and the U.S. Agency for International Development (USAID) which provides satellite-based Earth observation data, modeling, and science applications to help developing nations in Central America, East Africa and the Himalayas improve environmental decision making. Anticipating climate variability / climate change impacts has now become an important component of the SERVIR efforts to build capacity in these regions. Uncertainty in hydrometeorological components of climate variations and exposure to extreme events across scales from weather to climate are of particular concern. We report here on work to construct scenarios or outlooks that are being developed as input drivers for decision support systems (DSSs) in a variety of settings. These DSSs are being developed jointly by a broad array NASA Applied Science Team (AST) Investigations and user communities in the three SERVIR Hub Regions, Central America, East Africa and the Himalayas. Issues span hydrologic / water resources modeling, agricultural productivity, and forest carbon reserves. The scenarios needed for these efforts encompass seasonal forecasts, interannual outlooks, and likely decadal / multi-decadal trends. Providing these scenarios across the different AST efforts enables some level of integration in considering regional responses to climate events. We will discuss a number of challenges in developing this continuum of scenarios including the identification and "mining" of predictability, addressing multiple continental regions, issues of downscaling global model integrations to regional / local applications (i.e. hydrologic and crop modeling). We compare / contrast the role of the U.S. National Multi- Model Experiment initiative in seasonal forecasts and the CMIP-5 climate model experiments in supporting these efforts. Examples of these scenarios, their use, and an assessment of their utility as well as limitations will be presented.

Robertson, Franklin R.↗
NASA NTRS1987

Emission-line variability and the broad-line region of quasi-stellar objects. I - Time scales and photon densities

Spectroscopic data for five low-redshift QSOs taken over several years are presented as evidence for emission-line variability in QSOs. The broad hydrogen emission lines in the QSOs 0026 + 129, 1202 + 281, 1612 + 261, and 2141 + 175 are found to vary in step with their underlying continua. Another QSO, 2135 - 147, shows an increase in the intensities of its hydrogen lines five years after the increase in its continuum level. The inferred upper limit to the size of the broad emission line region is smaller than 3 lt-yr, probably as small as 0.4 lt-yr for 0026 + 129, if the line-emitting region is spherically distributed. Such a small size implies a very high density of ionizing photons, of the order of 10 to the 10th/cu cm. A very high nucleon density must then be assumed in order to obtain a value of the ionization parameter consistent with current photoionization models.

Zheng, Wei↗
NASA NTRS1994

The morphology of the north Jovian ultraviolet aurora observed with the Hubble Space Telescope

A series of six images covering a complete rotation of the north polar region of Jupiter were obtained in February 1993 with the Faint Object Camera on board the Hubble Space Telescope (HST). These images provide the first global picture of the morphology of the Jovian ultraviolet aurora observed from Earth orbit. The camera passband was centered near 153 nm, a region dominated by the H2 Lyman bands and continuum. The successive exposures, taken approximately 90 min apart, are used to construct a polar view of the auroral zone. It is found that the auroral emissions do not exactly follow the footprint of a constant L-shell although the size of the oval and its location agree best with the footprints of the approximately equal to 30 R(sub J) field line in the GSFC O6 model of the Jovian magnetic field. The displacement between the observed auroral zone and the theoretical oval may indicate a possible distortion of the Jovian magnetic field lines near the surface. A comparison with two images at the same wavelength obtained 8 months earlier shows that the main morphological features are persistent, in spite of changes in the detailed emission distribution. Small scale features with characteristic sizes of approximately 1000 km are observed along the auroral oval. The change of morphology observed as a function of the System 3 longitude appears as a persistent characteristic of the morphology of the north polar aurora.

Gerard, Jean-Claude↗
NASA NTRSDate not supplied

Design of a Subscale, Inert Gas Test for Plume-Surface Interactions in a Reduced Pressure Environment

Rocket plume-surface interaction is a multi-phase, complex discipline characterized by plume flow physics, erosion physics, and ejecta dynamics. To provide critical data for the validation of predictive modeling capabilities, an inert-gas, subscale ground test has been designed, with one test phase focused on cratering and ejecta and another test phase focused on plume structure and behavior. Using a supersonic, heated N2 jet impinging onto six different mechanical regolith simulants, visual data on the temporal evolution of cratering and ejecta are collected in a vacuum chamber at ambient pressures ranging from 600 Pa down to approximately 7 Pa. Using the same N2 jet impinging onto an instrumented plate, NO planar laser-induced fluorescence data are obtained at lower ambient pressures for the validation of both continuum and rarefied gas-phase physics. This paper discusses the design of this test and data products, with emphasis on the cratering and ejecta experiments, along with the physical scaling, challenges, and relevance of these data to the maturation of the plume-surface interaction discipline.

Ashley M Korzun↗
DOE OSTI2023

Challenges in studying water fluxes within the soil-plant-atmosphere continuum: A tracer-based perspective on pathways to progress

Tracing and quantifying water fluxes in the hydrological cycle is crucial for understanding the current state of ecohydrological systems and their vulnerability to environmental change. Especially the interface between ecosystems and the atmosphere that is strongly mediated by plants is important to meaningfully describe ecohydrological system functioning. Many of the dynamic interactions generated by water fluxes between soil, plant and the atmosphere are not well understood, which is partly due to a lack of interdisciplinary research. This opinion paper reflects the outcome of a discussion among hydrologists, plant ecophysiologists and soil scientists on open questions and new opportunities for collaborative research on the topic "water fluxes in the soil-plant-atmosphere continuum" especially focusing on environmental and artificial tracers. We emphasize the need for a multi-scale experimental approach, where a hypothesis is tested at multiple spatial scales and under diverse environmental conditions to better describe the small-scale processes (i.e., causes) that lead to large-scale patterns of ecosystem functioning (i.e., consequences). Novel in-situ, high-frequency measurement techniques offer the opportunity to sample data at a high spatial and temporal resolution needed to understand the underlying processes. Here we advocate for a combination of long-term natural abundance measurements and event-based approaches. Multiple environmental and artificial tracers, such as stable isotopes, and a suite of experimental and analytical approaches should be combined to complement information gained by different methods. Virtual experiments using process-based models should be used to inform sampling campaigns and field experiments, e.g., to improve experimental designs and to simulate experimental outcomes. On the other hand, experimental data are a pre-requisite to improve our currently incomplete models. Interdisciplinary collaboration will help to overcome research gaps that overlap across different earth system science fields and help to generate a more holistic view of water fluxes between soil, plant and atmosphere in diverse ecosystems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

FE K EMISSION AND ABSORPTION FEATURES IN THE XMM-EPIC SPECTRUM OF THE SEYFERT GALAXY IC 4329A

We present a re-analysis of the XMM-Newton long-look of the X-ray bright Seyfert galaxy IC 4329a. The Fe K bandpass is dominated by two peaks, consistent with emission from neutral or near-neutral Fe Ka and KP. A relativistic diskline model whereby both peaks are the result of one doubly-peaked diskline profile is found to be a poor description of the data. Models using two relativistic disklines are found to describe the emission profile well. A low-inclination, moderately-relativistic dual-diskline model is possible if the contribution from narrow components, due to distant material, is small or absent. A high-inclination, moderately relativistic profile for each peak is possible if there are roughly equal contributions from both the broad and narrow components. Upper limits on Fe XXV and Fe XXVI emission and absorption at the systemic velocity of IC 4329a are obtained. We also present the results of RXTE monitoring of this source obtained so far; the combined XMM-Newton and RXTE data sets allow us to explore the time-resolved spectral behavior of this source on time scales ranging from hours to 2 years. We find no strong evidence for variability of the Fe Ka emission line on any time scale probed, likely due to the minimal level of continuum variability. We detect a narrow absorption line, at a energy of 7.68 keV in the rest frame of the source; its significance has been confirmed using Monte Carlo simulations. This feature is most likely due to absorption from Fe XXVI blueshifted to approximately 0.1c relative to the systemic velocity, making IC 4329a the lowest-redshift AGN known with a high-velocity, highly-ionized outflow component. As is often the case with similar outflows seen in high-luminosity quasars, the estimated mass outflow rate is larger than the inflow accretion rate, signaling that the outflow represents a substantial portion of the total energy budget of the AGN. The outflow could arise from a radiatively-driven disk wind, or it may be in the form of a discrete, transient blob of ejected material.

Markowitz, A.↗
NASA NTRS2010

Detection of the 158 Micrometers[CII] Transition at z=1.3: Evidence for a Galaxy-Wide Starburst

We report the detection of 158 micron [C II] fine-structure line emission from MIPS J142824.0+352619, a hyperluminous ( L(sub IR) approximates 10(exp 13) L (sub solar)) starburst galaxy at z = 1.3. The line is bright, corresponding to a fraction L(sub [Cu II] L(sub Fir) approximates 2 x 10(exp -3) of the far-IR (FIR) continuum. The [C II], CO, and FIR continuum emission may be modeled as arising from photodissociation regions (PDRs) that have a characteristic gas density of n approximates 10(exp 4.2) /cm(exp 3) , and that are illuminated by a far-UV radiation field approximately 10(exp 3.2) times more intense than the local interstellar radiation field. The mass in these PDRs accounts for approximately half of the molecular gas mass in this galaxy. The L(sub [CII])/L(sub FIR) ratio is higher than observed in local ultralummous infrared galaxies or in the few high-redshift QSOs detected in [C II], but the L(sub [CII])/L(sub FIR) and L(sub CO)/L(sub FIR) ratios are similar to the values seen in nearby starburst galaxies. This suggests that MIPS J142824.0+352619 is a scaled-up version of a starburst nucleus, with the burst extended over several kiloparsecs.

Hailey-Dunsheath, S.↗
NASA NTRS2018

TPSAS-NF1676L-30451-DND

The fundamental laws of the physics of fluid flow in the continuum regime are embedded in a somewhat "simple" set of partial differential equations describing the transport of mass, momentum, energy, and/orentropy. These basic equations can be manipulated to construct other important physical quantities such as vorticity that help further elucidate flow physics. Additionally, these equations can be nondimensionalized to identify important nondimensional groups governing compressible, reacting fluid flows, most prominently Reynolds, Mach, Froude, Prandtl, Lewis, Schmidt, Damkohler, and Zeldovich numbers. In this lecture, we present these governing equations and discuss the significance of the different terms and nondimensional groups. Furthermore, since these equations cannot be solved directly on current supercomputers for turbulent flows at Reynolds numbers of engineering interest, a scale reduction by averaging or filtering is introduced, which leads to a "closure" problem. The closure is then provided by introducing some common models. The current introduction is by no means complete, but it and the references herein may serve as a good starting point for learning about the most common approaches to turbulent reacting flows.

Tomasz (Tom) G Drozda↗
DOE OSTI2022

$O(N)$ ab initio calculation scheme for large-scale moiré structures

Here we present a two-step method specifically tailored for band structure calculation of the small-angle moiré-pattern materials which contain tens of thousands of atoms in a unit cell. In the first step, the self-consistent field calculation for the ground state is performed with the O(N) Krylov subspace method implemented in openmx. Second, the crystal momentum-dependent Bloch Hamiltonian and overlap matrix are constructed from the results obtained in the first step and only a small number of eigenvalues near the Fermi energy are solved with shift-invert and Lanczos techniques. By systematically tuning two key parameters, the cutoff radius for electron hopping interaction and the dimension of the Krylov subspace, we obtained the band structures for both rigid and corrugated twisted bilayer graphene structures down to the first magic angle (θ = 1.08°) with high enough accuracy at affordable costs. The band structures are in good agreement with those from tight-binding models, continuum models, plane-wave pseudopotential based ab initio calculations, and experimental observations. This method is also shown to be efficient in twisted double-bilayer graphene and bilayer WSe 2 . We think this two-step method can play a crucial role in other twisted two-dimensional materials, especially those with much more complex band structure and where the effective model is hard to construct.

36 MATERIALS SCIENCE↗
NASA NTRS2022

Recommended Direct Simulation Monte Carlo Collision Model Parameters for Reacting Methane Flows

While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].

M R Gosma↗
NASA NTRS1991

Accretion dynamics and polarized X-ray emission of magnetized neutron stars

The basic ideas of accretion onto magnetized neutron stars are outlined. These are applied to a simple model of the structure of the plasma mound sitting at the magnetic poles of such a star, in which upward diffusion of photons is balanced by their downward advection. This steady flow model of the plasma's dynamical state is used to compute the emission of polarized X-raysfrom the optically thick, birefringent medium. The linear polarization of the continuum radiation emerging from the quasi-static mound is found to be as much as 40 percent at some rotation phases, but is insensitive to the geometry of the accretion flow. The role of the accretion shock, whose detailed polarimetric and spectral characteristics have yet to be calculated, is emphasized as the final determinant of the properties of the emerging X-rays. Some results describing the fully time dependent dynamics of the flow are also presented. In particular, steady flow onto a neutron star is shown to exhibit formation of 'photon bubbles', regions of greatly reduced plasma density filled with radiation which form and rise on millisecond time scale. The possible role of these complex structures in the flow for the formation of the emergent spectrum is briefly outlined.

Arons, Jonathan↗
NASA NTRS2011

Continuum Damage Mechanics Models for the Analysis of Progressive Failure in Open-Hole Tension Laminates

The performance of a state-of-the-art continuum damage mechanics model for interlaminar damage, coupled with a cohesive zone model for delamination is examined for failure prediction of quasi-isotropic open-hole tension laminates. Limitations of continuum representations of intra-ply damage and the effect of mesh orientation on the analysis predictions are discussed. It is shown that accurate prediction of matrix crack paths and stress redistribution after cracking requires a mesh aligned with the fiber orientation. Based on these results, an aligned mesh is proposed for analysis of the open-hole tension specimens consisting of different meshes within the individual plies, such that the element edges are aligned with the ply fiber direction. The modeling approach is assessed by comparison of analysis predictions to experimental data for specimen configurations in which failure is dominated by complex interactions between matrix cracks and delaminations. It is shown that the different failure mechanisms observed in the tests are well predicted. In addition, the modeling approach is demonstrated to predict proper trends in the effect of scaling on strength and failure mechanisms of quasi-isotropic open-hole tension laminates.

Song, Kyonchan↗
DOE OSTI2025

Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling

Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.

Dynamic simulation↗
DOE OSTI2025

A review of thermo-hydro-mechanical modeling of coupled processes in fractured rock: From continuum to discontinuum perspective

Coupled thermo-hydro-mechanical (THM) processes in fractured rock are playing a crucial role in geoscience and geoengineering applications. Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling. This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions. By comparing various models, one can better assess the uncertainties in predictions, particularly those related to the conceptual models. The review explores how these methodologies have significantly enhanced the fundamental understanding of how fractures respond to fluid injection and production, and improved predictive capabilities pertaining to coupled processes within fractured systems. It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience. The selection and calibration of model parameters should be based on typical ranges and applied to the specific conditions of applications. The challenges arising from inherent heterogeneity and uncertainties, nonlinear THM coupled processes, scale dependence, and computational limitations in representing field scale fractures are discussed. Realizing potential advances on computational capacity calls for methodical conceptualization, mathematical modeling, selection of numerical solution strategies, implementation, and calibration to foster simulation outcomes that intricately reflect the nuanced complexities of geological phenomena. Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical field of study.

Coupling scheme↗
DOE OSTI2021

Plant-microbe interactions: from genes to ecosystems using Populus as a model system

Plant-microbe symbioses span a continuum from pathogenic to mutualistic with functional consequences for both organisms in the symbiosis. In order to increase sustainable food and fuel production in the future, it is imperative that we harness these symbioses. The tree genus Populus is an excellent model system for studies examining plant-microbe interactions due to the wealth of genomic information available and the molecular tools that have been developed to manipulate Populus-microbe symbioses. In this review, we highlight how Populus can serve as a model system to explore plant-microbe interactions. Specifically, highlighting research linking Populus-microbe interactions from the gene to the ecosystem level. We explore why Populus is an excellent model for perennial plant systems, the molecular underpinnings of Populus-microbe interactions, how host genetics influence microbial community composition, and how microbial communities vary at fine spatial scales and between Populus species. Further, we explore how the patterns of the microbiome may affect ecosystem level functions in managed and natural ecosystems. Understanding and manipulating these interactions in Populus has the potential to improve plant health and impact ecosystem sustainability and processes as Populus trees function as foundational species in many natural ecosystems and are also deployed in managed ecosystems for various agroforestry applications.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2012

Constraining Aerosol Optical Models Using Ground-Based, Collocated Particle Size and Mass Measurements in Variable Air Mass Regimes During the 7-SEAS/Dongsha Experiment

During the spring of 2010, NASA Goddard's COMMIT ground-based mobile laboratory was stationed on Dongsha Island off the southwest coast of Taiwan, in preparation for the upcoming 2012 7-SEAS field campaign. The measurement period offered a unique opportunity for conducting detailed investigations of the optical properties of aerosols associated with different air mass regimes including background maritime and those contaminated by anthropogenic air pollution and mineral dust. What appears to be the first time for this region, a shortwave optical closure experiment for both scattering and absorption was attempted over a 12-day period during which aerosols exhibited the most change. Constraints to the optical model included combined SMPS and APS number concentration data for a continuum of fine and coarse-mode particle sizes up to PM2.5. We also take advantage of an IMPROVE chemical sampler to help constrain aerosol composition and mass partitioning of key elemental species including sea-salt, particulate organic matter, soil, non sea-salt sulphate, nitrate, and elemental carbon. Our results demonstrate that the observed aerosol scattering and absorption for these diverse air masses are reasonably captured by the model, where peak aerosol events and transitions between key aerosols types are evident. Signatures of heavy polluted aerosol composed mostly of ammonium and non sea-salt sulphate mixed with some dust with transitions to background sea-salt conditions are apparent in the absorption data, which is particularly reassuring owing to the large variability in the imaginary component of the refractive indices. Extinctive features at significantly smaller time scales than the one-day sample period of IMPROVE are more difficult to reproduce, as this requires further knowledge concerning the source apportionment of major chemical components in the model. Consistency between the measured and modeled optical parameters serves as an important link for advancing remote sensing and climate research studies in dynamic aerosol-rich environments like Dongsha.

Bell, Shaun W.↗
DOE OSTI2020

Breakdown of electroneutrality in nanopores

Ion transport in extremely narrow nanochannels has gained increasing interest in recent years due to unique physical properties at the nanoscale and the technological advances that allow us to study them. It is tempting to approach this confined regime with the theoretical tools and knowledge developed for membranes and microfluidic devices, and naively apply continuum models, such as the Poisson-Nernst-Planck and Navier-Stokes equations. However, it turns out that some of the most basic principles we take for granted in larger systems, such as the complete screening of surface charge by counter-ions, can break down under extreme confinement. Additionally, we show that in a truly one-dimensional system of ions interacting with three-dimensional electrostatic interactions, the screening length is exponentially large, and can easily exceed the macroscopic length of a nanotube. Without screening, electroneutrality breaks down within the nanotube, with fundamental consequences for ion transport and electrokinetic phenomena. In this work, we build a general theoretical framework for electroneutrality breakdown in nanopores, focusing on the most interesting case of a one-dimensional nanotube, and show how it provides an elegant interpretation for the peculiar scaling observed in experimental measurements of ionic conductance in carbon nanotubes.

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

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗