Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex simulations”

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,729 records · Page 96

Solid-state NMR and theoretical studies illuminate lanthanum borohydride C–H borylation catalysts confined within a zeolite

Zeolite-supported single-site lanthanum borohydride catalyzes C–H borylation of hydrocarbons, while the related silica-supported complex is inactive under comparable conditions. The identical composition of support-La(BH 4 ) 2 (THF) 2 sites in the two materials implies that the support bestows underlying structural features onto La that are required for bond activation catalysis, yet the nature of these effects, which could include confinement effects and/or electronic modulation of the site itself, remains to be identified. We used solid-state nuclear magnetic resonance (SSNMR) spectroscopy and molecular dynamics simulations with machine-learning potentials (ML-MD) to analyze the electronic and steric effects imparted by the faujasite support on the precatalyst structure to correlate with catalytic activity. ML-MD simulations show that THF dissociates from La under the influence of confinement, leading to coordinatively unsaturated sites in the zeolite pores. Then, the La complex grafts on Brønsted acid sites (La BAS ) or isolated silanols (La SiO ) or remains physisorbed in the zeolite pores. Catalytic studies comparing compounds supported on faujasite zeolites containing or lacking BAS and/or silanols show that only the former complexes lead to active sites, ruling out confinement as the sole requirement for catalysis. The DFT calculations and ML-MD simulations also reveal that the surface-lanthanum coordination number is two (bidentate) for La BAS , with the metal center forming long, flexible bonds to two oxygen atoms bridging Si and Al, but only one oxygen atom (monodentate) for La SiO . The structure–activity relationship identifies confined, BAS-grafted species as active sites and provides important guidance for the design of enhanced atom-efficient catalysts.

Cui, Jinlei [Ames Laboratory, and Iowa State Univ.↗

Impact of boundary layer simulation on predicting radioactive pollutant dispersion: a case study for HANARO Research Reactor using the WRF-MMIF-CALPUFF modeling system

Wind plays an important role in cases of unexpected radioactive pollutant dispersion, deciding distribution and concentration of the leaked substance. The accurate prediction of wind has been challenging in numerical weather prediction models, especially near the surface because of the complex interaction between turbulent flow and topographic effect. As such, in this study, we investigated the characteristics of atmospheric dispersion of radioactive material (i.e. 137 Cs) according to the simulated boundary layer around the HANARO research nuclear reactor in Korea using the Weather Research and Forecasting (WRF)-Mesoscale Model Interface (MMIF)-California Puff (CALPUFF) model system. We examined the impacts of orographic drag on wind field, stability calculation methods, and planetary boundary layer parameterizations on the dispersion of radioactive material under a radioactive leaking scenario. We found that inclusion of the orographic drag effect in the WRF model improved the wind prediction most significantly over the complex terrain area, leading the model system to estimate the radioactive concentration near the reactor more conservatively. We also emphasized the importance of the stability calculation method and employing the skillful boundary layer parameterization to ensure more accurate low atmospheric conditions, in order to simulate more feasible spatial distribution of the radioactive dispersion in leaking scenarios.

54 ENVIRONMENTAL SCIENCES↗

Plasma Surface Interactions: Predicting the Performance and Impact of Dynamic PFC Surfaces

The objective of this project is to develop, and integrate, high-performance simulation tools capable of predicting plasma-facing component (PFC) operating lifetime and the impact of the evolving surface morphology of tungsten-based PFCs on plasma contamination, including the dynamic recycling of fuel species and tritium retention, in future magnetic fusion devices. Establishing a fundamental physical understanding and developing predictive capabilities of plasma-surface interactions (PSI) requires simultaneously addressing complex and diverse physics occurring over a wide range of length (Angstroms to meters) and time (femtoseconds to years) scales, as well as integrating extensive physical processes across the plasma–surface interface. This requires development of not only detailed physics models and computational strategies at each scale, but also algorithms and methods to couple them effectively in a way that can be robustly validated. Deploying these tools requires the continued development and coupling of leadership-scale computational codes to describe the boundary plasma and the evolving PFC surface, as well as a host of simulations that bridge disparate scales to address complex physical and computational issues at the plasma–surface interface in multi-component materials systems for magnetic fusion energy development beyond ITER.

36 MATERIALS SCIENCE↗

DECA: Discrete Event inspired Cellular Automata for grain structure prediction in additive manufacturing

Microstructure largely dictates macroscopic material properties and is strongly affected by processing. Therefore, the simulation of microstructure evolution in response to thermal fields during processing is of significant interest within the computational materials science community. Additive manufacturing (AM) has emerged as a technique for producing complex geometries and unique microstructures. Yet, complex and rapid thermal cycles in AM pose computational challenges for existing microstructure models. This work proposes a discrete event inspired cellular automata (CA) approach, titled DECA, to accelerate simulation of grain structure evolution in AM. In contrast to conventional time-stepped CA models, this model directly solves the times capture events would take place allowing for stepping in events rather than time (a technique also found in the field of discrete-event simulation). In comparison to purely serial discrete-event models, DECA allows for temporary violation of the causality constraint, but detects and corrects these violations, leading to an emergent phenomenon dubbed causality rippling, in which previously calculated capture events are overwritten. The amount of repeated calculations, defined by the capture ratio, is taken as a measure of computational inefficiency, and the model parameters that affect this ratio are evaluated. The new DECA approach was found to be more computationally efficient than conventional time-stepped CA models while guaranteeing an accurate solution, which can only be achieved in the conventional models for vanishingly small time steps. Finally, opportunities for parallelization and scaling of the new approach are discussed.

36 MATERIALS SCIENCE↗

An atomistic survey of shear coupling in asymmetric tilt grain boundaries and interpretation using the disconnections framework

Grain Boundaries (GB) play an important role in determining the behavior of polycrystalline materials. While the mechanisms of motion and associated shear response for symmetric tilt grain boundaries (STGBs) are well studied, the same is not true for asymmetric tilt grain boundaries (ATGBs) despite their greater prevalence in polycrystals. Here, this study aims to investigate the shear response of a large collection of asymmetric tilt grain boundaries (ATGBs) using molecular dynamics (MD) simulations and interpret the data using a discrete disconnections model that works remarkably well for STGBs. MD simulations of shear-driven ATGBs show that the plastic shear (shear coupling factor) in the region swept by a GB exhibits a complex dependence on the inclination angle, and this dependence changes with the misorientation of the GB. In addition, the shear response was observed to be highly sensitive to the applied shear rate and temperature. Recognizing the spatial and temporal scale limitations of MD simulations, we extended the discrete disconnections mesoscale model of Khateret al. (2012) to calculate the nucleation barriers of disconnection modes and predict the effective shear coupling of an ATGB. We observed that the mesoscale model’s predictions of the shear coupling factor of ATGBs do not agree with those observed in MD simulations. Finally, we examine the hypotheses of our mesoscale model that contribute to disagreements between MD simulations and the mesoscale model and propose improvements to the mesoscale model for future work.

36 MATERIALS SCIENCE↗

Image Distinguishability Analysis Testing Through Principal Components and Its Application to Hot Spot Scale Invariance

Hot spots are spatial regions of intense energy localization that govern initiation of secondary high explosives. Studies that characterize or compare simulated hot spots are frequently either qualitatively descriptive or resort to quantitative distribution functions that neglect stochastic variations and spatial correlations—effects that are also neglected in common comparison tests like the Kolmogorov–Smirnov test. To this end, we develop an image distinguishability analysis (IDA) test based on principal component (PC) analysis that makes pixel-by-pixel comparisons between small, for example, O(<10), image data sets. The IDA test makes comparisons through a generalized distance metric in the PC space and a test statistic that is derived to calculate mathematical equation-values. Here, we derive a statistical distribution and criticality criterion to determine whether images are distinguishable from established baselines. We apply the IDA test on images generated from molecular dynamics simulations of hot spots from pore collapse in TATB to assess scale invariance in the complex patterns of hot spots that form in a representative high explosive crystal. The IDA test shows that TATB hot spot spatial temperature fields and their derived temperature histograms exhibit scale-invariant features over specific intervals of shock orientation, strength, and initial pore diameter. However, the IDA test also shows that qualitatively different conclusions regarding invariance can be reached depending on whether the hot spot is treated as a spatially correlated field as opposed to a distribution function that lacks spatial information.

organic↗

Understanding interactions between biomolecules and two-dimensional nanomaterials using in silico microscopes

We report two-dimensional (2D) nanomaterials such as graphene are increasingly used in research and industry for various biomedical applications. Extensive experimental and theoretical studies have revealed that 2D nanomaterials are promising drug delivery vehicles, yet certain materials exhibit toxicity under biological conditions. So far, it is known that 2D nanomaterials possess strong adsorption propensities for biomolecules. To mitigate potential toxicity and retain favorable physical and chemical properties of 2D nanomaterials, it is necessary to explore the underlying mechanisms of interactions between biomolecules and nanomaterials for the subsequent design of biocompatible 2D nanomaterials for nanomedicine. The purpose of this review is to integrate experimental findings with theoretical observations and facilitate the study of 2D nanomaterial interaction with biomolecules at the molecular level. We discuss the current understanding and progress of 2D nanomaterial interaction with proteins, lipid membranes, and DNA based on molecular dynamics (MD) simulation. In this review, we focus on the 2D graphene nanosheet and briefly discuss other 2D nanomaterials. With the ever-growing computing power, we can image nanoscale processes using MD simulation that are otherwise not observable in experiment. We expect that molecular characterization of the complex behavior between 2D nanomaterials and biomolecules will help fulfill the goal of designing effective 2D nanomaterials as drug delivery platforms.

2D nanomaterials↗

Ten questions concerning agent-based modeling of occupant behavior for energy and environmental performance of buildings

We report the complexity of occupant behavior is one of the major contributors to uncertainty in building performance simulation. Agent-based modeling (ABM), a computational simulation technique, has gained attention in the occupant modeling field due to its capability and flexibility to capture the heterogeneity and dynamics of human behavior and the emergent effects. While multiple efforts in the past decade have demonstrated the usefulness of the ABM approach for simulating occupants and their impacts on building performance, several crucial matters in the ABM research still remain unexplored. This paper presents ten questions that highlight the most important issues regarding ABM research and applications for occupant behavior in the context of building performance simulation. The questions and answers aim to provide insights into current and future ABM research, and more importantly to inspire new significant questions from young researchers in the field. This research is part of the IEA EBC Annex 79 project, occupant-centric building design and operation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Electro-reduction of organics on metal cathodes: A multiscale-modeling study of benzaldehyde on Au (111)

Here, we present a multiscale modeling study on the electro-reduction of benzaldehyde (BZY) on Au(111) under realistic electrochemical conditions. To model the electrochemical cell, we adopt a capacitor model in which complex solvents are confined between a Au cathode and a carbon anode. Classical molecular dynamics simulations reveal that electrode charge density and the presence of alcohol show strong effects on the density, adsorption geometry and dynamics of benzaldehyde on the Au electrode. Under charging conditions, the surface concentration of benzaldehyde on the Au electrode decreases, while the content of other species increases. Finally, we proposed a scheme that correlates the electric current running through the Au/solvent interface with the applied bias. This study provides a molecular level understanding of how solvent composition, in this case water/alcohol content, controls the activity of electrocatalytic hydrogenation

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanical behavior analysis using small-size graphite disc compression testing with digital image correlation methods

The ASTM International Standard Test Method for Tensile Stress-Strain of Carbon and Graphite (ASTM C749) requires dog-bone shaped specimens. Unfortunately, specimen sizes required in ASTM C749 are much too large for irradiation within a test reactor and are incompatible with irradiation capsule volumes, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed as an alternative that provides a means for testing tensile properties using smaller specimens when geometric constraints are present. However, the important stress–strain relation throughout the test loading history that could be recorded using ASTM C749 is lost by using ASTM D8289 because of the difficulties involved with using traditional strain gauges or extensometers on such a small specimen. In this work, a digital image correlation (DIC) method, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 splitting tensile test on eight small samples of graphite grade Mersen 2114 to track their deformation during compressive loading. Initial results showed that the DIC technique can measure the surface displacement/strain on these small (Ø6 × 3 mm) graphite specimens with high accuracy and good repeatability, indicating that the DIC technique, combined with the small-disc splitting tensile test in ASTM D8289, can capture the strain evolution in the test. In conclusion, further comparison of measured strains with analytical and numerical simulation results provides additional understanding of the mechanical behaviour of this nuclear graphite material in complex stress states.

36 MATERIALS SCIENCE↗

Dynamics of polymer segments, polymer chains, and nanoparticles in polymer nanocomposite melts: A review

The addition of nanoparticles (NPs) to a polymer matrix, forming a polymer nanocomposite (PNC), is known to alter the microscopic dynamic processes of both species which leads to unique macroscopic material properties of the PNC. Because the NPs and polymers have overlapping characteristic length, time, and energy scales, the interactions within these materials are complex and the dynamics are interrelated. In this review, we present an overview of experimental, simulation, and theoretical results that probe multi-scale polymer and nanoparticle dynamics in polymer nanocomposites and navigate the dense parameter space presented by these multicomponent systems. Although a variety of PNC systems are mentioned, we focus this discussion on linear thermoplastics filled with hard spherical or cylindrical NPs in the melt state. We begin by introducing PNCs, the dynamic processes within them, and the importance of dynamics for properties and processing. At the smallest length and time scales, we discuss segmental dynamics in PNCs, including the role of polymer attributes, NP attributes, and NP-polymer interactions. Then, we present measurements of collective motions and intermediate (Rouse) dynamics in various PNC materials. At longer length and time scales, we discuss polymer center-of-mass diffusion in PNCs with either spherical or anisotropic NPs. We then explore the dynamics of the NPs in PNCs and polymer melts, including theoretical predictions, simulation results, and experimental observations. Lastly, we note some of the remaining challenges in probing dynamics in PNC materials and fundamentally studying PNCs more generally.

36 MATERIALS SCIENCE↗

Toward accurate measurement of electromagnetic field by retrieving and refining the center position of non-uniform diffraction disks in Lorentz 4D-STEM

Recent advancement in scanning transmission electron microscopy (STEM) allows the use of 4D-STEM, a technique that captures an electron diffraction pattern at each scan point in STEM, to measure electrostatic and magnetic potential and field in materials. However, accurate measurement, separation of the magnetic and electric signals, and removal of artifacts remain challenging, especially in the presence of complex non-uniform diffraction contrast within the disks. In this work, based on dynamic simulations of 4D-STEM patterns built upon superstructures consisting of millions of atoms to account for different sample thickness and edge geometries, we show how the shape and intensity distribution of the central disk are affected by multiple scattering. We propose a robust refinement procedure through iteration of the spin-sensitive peak position of the disk-center in the circular Hough transform filtered images from experimental Lorentz 4D-STEM dataset after minimizing the possible artifacts, such as those due to the change of thickness, dynamic scattering, and scanning process. We verify that caution must be taken as in practice the rigid-disk-shift model used to reconstruct induction maps can easily break down due to disk-protrusion when there exists a nonconstant phase gradient or thickness within the width of the probe. Through quantitative analysis and comparing experiment with calculation the effect of the non-spin-related intensity distribution inside the disk as well as that causes the disk shift due to the intensity-protrusion can be removed, and high-quality magnetic field mapping is possible.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Restricted Rotational Flexibility of the C5α-Methyl-Substituted Carbapenem NA-1-157 Leads to Potent Inhibition of the GES-5 Carbapenemase

Carbapenem antibiotics are used as a last-resort treatment for infections caused by multidrug-resistant bacteria. The wide spread of carbapenemases in Gram-negative bacteria has severely compromised the utility of these drugs and represents a serious public health threat. To combat carbapenemase-mediated resistance, new antimicrobials and inhibitors of these enzymes are urgently needed. Here, we describe the interaction of the atypically C5α-methyl-substituted carbapenem, NA-1-157, with the GES-5 carbapenemase. MICs of this compound against Escherichia coli, Klebsiella pneumoniae, and Acinetobacter baumannii producing the enzyme were reduced 4–16-fold when compared to MICs of the commercial carbapenems, reaching clinically sensitive breakpoints. When NA-1-157 was combined with meropenem, a strong synergistic effect was observed. Kinetic and ESI-LC/MS studies demonstrated that NA-1-157 is a potent inhibitor of GES-5, with a high inactivation efficiency of (2.9 ± 0.9) × 105 M –1 s –1 . Acylation of GES-5 by NA-1-157 was biphasic, with the fast phase completing within seconds, and the slow phase taking several hours and likely proceeding through a reversible tetrahedral intermediate. Deacylation was extremely slow (k 3 = (2.4 ± 0.3) × 10 –7 s –1 ), resulting in a residence time of 48 ± 6 days. MD simulation of the GES-5-meropenem and GES-5-NA-1-157 acyl-enzyme complexes revealed that the C5α-methyl group in NA-1- 157 sterically restricts rotation of the 6α-hydroxyethyl group preventing ingress of the deacylating water into the vicinity of the scissile bond of the acyl-enzyme intermediate. Furthermore, these data demonstrate that NA-1-157 is a potent irreversible inhibitor of the GES-5 carbapenemase.

60 APPLIED LIFE SCIENCES↗

Impact of ion fluxes across thylakoid membranes on photosynthetic electron transport and photoprotection

In photosynthetic thylakoid membranes the proton motive force (pmf) not only drives ATP synthesis, in addition it is central to controlling and regulating energy conversion. As a consequence, dynamic fine-tuning of the two pmf components, electrical (Δψ) and chemical (ΔpH), is an essential element for adjusting photosynthetic light reactions to changing environmental conditions. Good evidence exists that the Δψ/ΔpH partitioning is controlled by thylakoid potassium and chloride ion transporters and channels. However, a detailed mechanistic understanding of how these thylakoid ion transporter/channels control pmf partitioning is lacking. In this work, we combined functional measurements on potassium and chloride ion transporter and channel loss-of-function mutants with extended mathematical simulations of photosynthetic light reactions in thylakoid membranes to obtain detailed kinetic insights into the complex interrelationship between membrane energization and ion fluxes across thylakoid membranes. The data reveal that potassium and chloride fluxes in the thylakoid lumen determined by the K + /H + antiporter KEA3 and the voltage-gated Cl – channel VCCN1/Best1 have distinct kinetic responses that lead to characteristic and light-intensity-dependent Δψ/ΔpH oscillations. These oscillations fine-tune photoprotective mechanisms and electron transport which are particularly important during the first minutes of illumination and under fluctuating light conditions. By employing the predictive power of the model, we unravelled the functional consequences of changes in KEA3 and VCCN1 abundance and regulatory/enzymatic parameters on membrane energization and photoprotection.

59 BASIC BIOLOGICAL SCIENCES↗

Multi-model Hydroclimate Projections for the Alabama-Coosa-Tallapoosa River Basin in the Southeastern United States

Abstract This study uses a high-resolution, process-based modeling framework to assess the impacts of changing climate on water resources for the Alabama-Coosa-Tallapoosa River Basin in the southeastern United States. A 33-member ensemble of hydrologic projections was generated using 3 distributed hydrologic models (Precipitation-Runoff Modeling System, Variable Infiltration Capacity, and Distributed Hydrology Soil Vegetation Model) of different complexity. These hydrologic models were driven by dynamically downscaled and bias-corrected future climate simulations from 11 Coupled Model Intercomparison Project Phase 5 global climate models under Representative Concentration Pathway 8.5 emission scenario, with 40 years each in baseline (1966–2005) and future (2011–2050) periods. The hydroclimate response, in general, projects an increase in mean seasonal precipitation, runoff, and streamflow. The high and low flows are projected to increase and decrease, respectively, in general, suggesting increased likelihood of extreme rainfall events and intensification of the hydrologic cycle. The uncertainty associated with the ensemble hydroclimate response, analyzed through an analysis of variance technique, suggests that the choice of climate model is more critical than the choice of hydrologic model for the studied region. This study provides in-depth insights of hydroclimate response and associated uncertainties to support informed decisions by water resource managers.

54 ENVIRONMENTAL SCIENCES↗

Investigation of local distortion effects on X-ray absorption of ferroelectric perovskites from first principles simulations

Understanding the role of ferroelectric polarization in modulating the electronic and structural properties of crystals is critical for advancing these materials for overcoming various technological and scientific challenges. However, due to difficulties in performing experimental methods with the required resolution, or in interpreting the results of methods therein, the nanoscale morphology and response of these surfaces to external electric fields has not been properly elaborated. Here, in this work, we investigate the effect of ferroelectric polarization and local distortions in a BaTiO 3 perovskite, using two widely used computational approaches which treat the many-body nature of X-ray excitations using different philosophies, namely the many-body, delta-self-consistent-field determinant (mb-ΔSCF) and the Bethe–Salpeter equation (BSE) approaches. We show that in agreement with our experiments, both approaches consistently predict higher excitations of the main peak in the O–K edge for the surface with upward polarization. However, the mb-ΔSCF approach mostly fails to capture the L 2,3 separations at the Ti–L edge, due to the absence of spin–orbit coupling in Kohn–Sham density functional theory (KS-DFT) at the generalized gradient approximation level. On the other hand, and most promising, we show that application of the GW/BSE approach successfully reproduces the experimental XAS, both the relative peak intensities as well as the L 2,3 separations at the Ti–L edges upon ferroelectric switching. Thus simulated XAS is shown to be a powerful method for capturing the nanoscale structure of complex materials, and we underscore the need for many-body perturbation approaches, with explicit consideration of core-hole and multiplet effects, for capturing the essential physics in these systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Dust and gas absorption in the high mass X-ray binary IGR J16318-4848

With an absorption column density on the order of 10 24 cm -2 , IGR J16318-4848 is one of the most extreme cases of a highly obscured high mass X-ray binary. In addition to the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. We aim at a self-consistent description of the continuum absorption, the absorption edges, and the fluorescence lines to constrain the properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescence emission from dust grains. We used XMM-Newton and NuSTAR spectra to first empirically constrain the incident continuum and fluorescence lines. Next we used XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In the third step we modeled the absorption and fluorescence in a dusty olivine absorber and employed both a simple analytical model for the fluorescence line emission and a Monte Carlo simulation of radiative transfer that generates line fluxes, which are very close to the observational data. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is able to produce the observed continuum absorption and edges. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318-4848. In particular, for low optical depth of individual grains, which is the case for typical volume densities and grain size distribution models, the dust will contribute significantly to the fluorescence emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Self-consistent calculations of the electric charge, ion drag force, and the drift velocity of spherical grains using Langevin dynamics and comparisons against canonical experiments

We present trajectory simulation-based modeling to capture the interactions between ions and charged grains in dusty or complex plasmas. Our study is motivated by the need for a self-consistent and experimentally validated approach for accurately calculating the ion drag force and grain charge that determine grain collective behavior in plasmas. We implement Langevin dynamics in a computationally efficient predictor–corrector approach to capture multiscale ion and grain dynamics. Predictions of grain velocity, grain charge, and ion drag force are compared with prior measurements to assess our approach. The comparisons reveal excellent agreement to within ±20% between predicted and measured grain velocities [Yaroshenko et al., Phys. Plasmas 12, 093503 (2005) and Khrapak et al., Europhys. Lett. 97, 35001 (2012)] for 0.64, 1.25 μm grains at ∼20−500 Pa. Comparisons with the measured grain charge [Khrapak et al., Phys. Rev. E 72, 016406 (2005)] under similar conditions reveal agreement to within ∼20% as well. Measurements of the ion drag force [Hirt et al., Phys. Plasmas 11, 5690 (2004); IEEE Trans. Plasma Sci. 32, 582 (2004)] are used to assess the viability of the presented approach to calculate the ion drag force experienced by grains exposed to ion beams of well-defined energy. Excellent agreement between calculations and measurements is obtained for beam energies >10 eV, and the overprediction below 10 eV is attributed to the neglect of charge exchange collisions in our modeling. Along with critical assessments of our approach, suggestions for future experimental design to probe charging of and momentum transfer onto grains that capture the effect of space charge concentration and external fields are outlined.

Physics↗