Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coarse-grained 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Correction to: Investigating the rheology of fluidized and non-fluidized gas-particle beds: implications for the dynamics of geophysical flows and substrate entrainment

We report natural geophysical mass flows are among the most complex granular systems and their dynamics are often modified by the presence of an interstitial fluid. Prediction of their runout requires the development of models estimating the solid stresses in these hazardous currents wherein excess pore-fluid pressure can develop. We use discrete element modelling (DEM-CFD) with a Coarse-Graining post-processing technique (CG) to investigate the rheology of unsteady gas-particle fluidized to non-fluidized granular beds placed on horizontal and inclined planes. Similar to fluidized beds immersed in viscous fluids, the effective friction coefficient of air-fluidized beds can be defined as a function of the classic μ(I)-rheology and the non-dimensional fluid or solid pressure to explain the failure and dynamics of granular flows with excess pore pressure on inclines. However, dilation imposed by fluid drag and particle collisions in gas-particle fluidized beds can drastically change its effective frictional properties. In contrast with the common assumption in water-particle flows that granular temperature is negligible, in our gas-particle simulations, the contribution of the velocity fluctuations to the stress tensor is significant. Hence, the shear stress is found to be non-zero even when the flow is fully fluidized in the inertial regime. These results suggest the need to better understand velocity fluctuations to predict the effective viscosity of sheared fluidized mixtures and are relevant for many applications. Notably, a unified approach is useful for many geophysical flows that encompass a range of fluidization conditions in a single flow such as pyroclastic density currents and snow avalanches.

58 GEOSCIENCES↗
DOE OSTI2022

An approximate inertial manifold (AIM) based closure for turbulent flows

A closure model for turbulent flows is developed based on a dynamical system theory. An appropriately discretized formulation of the governing equations is considered for this process. The key ingredient is an approximation of the system’s attractor, where all the trajectories in phase space are confined. This approximate inertial manifold based approach provides a path to track trajectories of the system in a lower-dimensional subspace. Unlike conventional coarse-graining approaches, the turbulent field is decomposed into resolved and unresolved dynamics using the properties of the governing equations. The novelty of the approach relies on the reconstruction of the unresolved field constrained by the governing equations. A posteriori tests for homogeneous isotropic turbulence and the Kuramoto–Sivashinsky equation show promising results for considerable dimension reduction with strong convergence properties. The proposed model outperforms the dynamic Smagorinsky model, and the computational overhead is competitive with similar approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Scale interactions and anisotropy in Rayleigh–Taylor turbulence

We study energy scale transfer in Rayleigh–Taylor (RT) flows by coarse graining in physical space without Fourier transforms, allowing scale analysis along the vertical direction. Two processes are responsible for kinetic energy flux across scales: baropycnal work Λ, due to large-scale pressure gradients acting on small scales of density and velocity; and deformation work Π, due to multiscale velocity. Our coarse-graining analysis shows how these fluxes exhibit self-similar evolution that is quadratic-in-time, similar to the RT mixing layer. We find that Λ is a conduit for potential energy, transferring energy non-locally from the largest scales to smaller scales in the inertial range where Π takes over. In three dimensions, Π continues a persistent cascade to smaller scales, whereas in two dimensions Π rechannels the energy back to larger scales despite the lack of vorticity conservation in two-dimensional (2-D) variable density flows. This gives rise to a positive feedback loop in 2-D RT (absent in three dimensions) in which mixing layer growth and the associated potential energy release are enhanced relative to 3-D RT, explaining the oft-observed larger α values in 2-D simulations. Despite higher bulk kinetic energy levels in two dimensions, small inertial scales are weaker than in three dimensions. Moreover, the net upscale cascade in two dimensions tends to isotropize the large-scale flow, in stark contrast to three dimensions. Furthermore, our findings indicate the absence of net upscale energy transfer in three-dimensional RT as is often claimed; growth of large-scale bubbles and spikes is not due to ‘mergers’ but solely due to baropycnal work Λ.

42 ENGINEERING↗
DOE OSTI2023

Intergranular Hotspots: A Molecular Dynamics Study on the Influence of Compressive and Shear Work

Numerous crystal- and microstructural-level mechanisms are at play in the formation of hotspots, which are known to govern high explosives initiation behavior. Most of these mechanisms, including pore collapse, interfacial friction, and shear banding, involve both compressive and shear work done within the material and have thus far remained difficult to separate. Here, we assess hotspots formed at shocked crystal–crystal interfaces using quasi-1D molecular dynamics simulations that isolate effects due to compression and shear. Two high explosive materials are considered (TATB and PETN) that exhibit distinctly different levels of molecular conformational flexibility and crystal packing anisotropy. Temperature and intramolecular strain energy localization in the hotspot are assessed through parametric variation of the crystal orientation and two velocity components that respectively modulate compression and shear work. The resulting hotspots are found to be highly localized to a region within 5–20 nm of the crystal–crystal interface. Compressive work plays a considerably larger role in localizing temperature and intramolecular strain energy for both materials and all crystal orientations considered. Shear induces a moderate increase in energy localization relative to unsheared cases only for relatively weak compressive shock pressures of approximately 10 GPa. These results help isolate and rank the relative importance of hotspot generation mechanisms and are anticipated to guide the treatment of crystal–crystal interfaces in coarse-grained models of polycrystalline high-explosive materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗
DOE OSTI2023

Shadow Molecular Dynamics and Atomic Cluster Expansions for Flexible Charge Models

Here, a shadow molecular dynamics scheme for flexible charge models is presented where the shadow Born–Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, is modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born–Oppenheimer molecular dynamics (BOMD). XL-BOMD provides stable dynamics while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO 2 ) and a molecular system of liquid water. The combined ACE+XL-QEq molecular dynamics simulations are stable over a wide range of temperatures both for the oxide and for the molecular systems and provide a precise sampling of the Born–Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Radiation breakage of DNA: a model based on random-walk chromatin structure

Monte Carlo computer software, called DNAbreak, has recently been developed to analyze observed non-random clustering of DNA double strand breaks in chromatin after exposure to densely ionizing radiation. The software models coarse-grained configurations of chromatin and radiation tracks, small-scale details being suppressed in order to obtain statistical results for larger scales, up to the size of a whole chromosome. We here give an analytic counterpart of the numerical model, useful for benchmarks, for elucidating the numerical results, for analyzing the assumptions of a more general but less mechanistic "randomly-located-clusters" formalism, and, potentially, for speeding up the calculations. The equations characterize multi-track DNA fragment-size distributions in terms of one-track action; an important step in extrapolating high-dose laboratory results to the much lower doses of main interest in environmental or occupational risk estimation. The approach can utilize the experimental information on DNA fragment-size distributions to draw inferences about large-scale chromatin geometry during cell-cycle interphase.

Chromatin/radiation effects↗
DOE OSTI2021

Multiscale Concurrent Atomistic-Continuum (CAC) modeling of multicomponent alloys

We report strengthening in complex multicomponent systems such as solid solution alloys is controlled primarily by the dynamic interactions between dislocation lines and heterogeneously distributed solute species. Modeling of extended defect length scales in such multicomponent systems becomes prohibitively expensive, motivating the development of reduced order approaches. This work explores the application of the Concurrent Atomistic-Continuum (CAC) method to model dislocation mobility in random alloys at extended length scales. By employing recently developed average-atom interatomic potentials, the average “bulk” material response in coarse-grained regions interacts with true random solute species in the atomistic-scale domain. We demonstrate that spurious stresses in domain resolution transition regions are eliminated entirely due to the CAC formulation. Simultaneously, the key details of local stress fluctuation due to randomness in the dislocation core region are captured, and fluctuating stress smoothly decays to the long-range dislocation stress field response. Dislocation mobility calculations, for line lengths over 400 nm, are computed as a function of alloy composition in the model FeNiCr system and compared to full molecular dynamics (MD). The results capture the composition-dependent trends, while reducing degrees of freedom by nearly 40%. This approach can be readily extended to any system described by an EAM potential and facilitates the study of large-scale defect dynamics in complex solute environments to support computational alloy design.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Optical and TEM study of shock metamorphism from the Sedan test site

Thus far, detailed petrologic studies of shock metamorphism have been performed on samples recovered from laboratory experiments and on a few natural impactites. The loading history of these samples is quite different: In particular, laboratory experiments spend only a short time (less than 1 microsec) at peak pressure, whereas natural impactites may have stress pulses from 0.1 - 1 ms. On the other hand, laboratory experiments have known stress histories; natural impactites do not. Natural samples are also subjected to thousands or millions of years of postshock annealing and/or weathering. A useful intermediate case is that of nuclear detonation. Stress pulses for these events can reach 0.1 ms or higher, and samples are obtained in pristine condition. All three types of loading produce stresses of hundreds of kilobars. Samples studied were taken from the Sedan nulcear test site, and consist of a coarse-grained granodiorite containing quartz, K-feldspar, cordierite, and hornblende. Samples were studied optically in this section, then were thinned with an ion mill and studied by transmission electron microscopy (TEM). Optically, quartz and K-feldspar displayed numerous sets of planar deformation features (PDF's) identical to the nondecorated PDF's seen in laboratory samples and many natural impactites. TEM study showed that the PDF's in quartz and feldspar corresponded to densely packed wide transformation lamellae identical to those described in laboratory studies. The transformation lamellae in both minerals were amorphous, with no sign of high-pressure phases. In the case of K-feldspar only, narrow sublamellae extended outward from some wide lamellae. Quartz, which was more abundant and studied more extensively, contained no shock-induced dislocations. Some planar features were also seen in cordierite, but could not be identified due to rapid beam damage. No shock defects were seen in hornblende in TEM. The shock-induced defects present at the Sedan site are very similar to those seen in shock recovery experiments, and also to those present at certain natural events (e.g., Meteor Crater). This suggests that shock deformation in quartz is not strongly dependent on shock pulse duration, and that laboratory recovery experiments are useful simulations of natural impact events.

Gratz, A. J.↗
DOE OSTI2022

A Coarse-Grained Model of Affinity Maturation Indicates the Importance of B-Cell Receptor Avidity in Epitope Subdominance

The elicitation of broadly neutralizing antibodies (bnAbs) is a major goal in the design of vaccines against rapidly-mutating viruses. In the case of influenza, many bnAbs that target conserved epitopes on the stem of the hemagglutinin protein (HA) have been discovered. However, these antibodies are rare, are not boosted well upon reinfection, and often have low neutralization potency, compared to strain-specific antibodies directed to the HA head. Different hypotheses have been proposed to explain this phenomenon. We use a coarse-grained computational model of the germinal center reaction to investigate how B-cell receptor binding valency affects the growth and affinity maturation of competing B-cells. We find that receptors that are unable to bind antigen bivalently, and also those that do not bind antigen cooperatively, have significantly slower rates of growth, memory B-cell production, and, under certain conditions, rates of affinity maturation. The corresponding B-cells are predicted to be outcompeted by B-cells that bind bivalently and cooperatively. We use the model to explore strategies for a universal influenza vaccine, e.g., how to boost the concentrations of the slower growing cross-reactive antibodies directed to the stem. The results suggest that, upon natural reinfections subsequent to vaccination, the protectiveness of such vaccines would erode, possibly requiring regular boosts. Collectively, our results strongly support the importance of bivalent antibody binding in immunodominance, and suggest guidelines for developing a universal influenza vaccine.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Coarse-Grained Density Functional Theory Predictions via Deep Kernel Learning

Scalable electronic predictions are critical for soft materials design. Recently, the Electronic Coarse-Graining (ECG) method was introduced to renormalize all-atom quantum chemical (QC) predictions to coarse-grained (CG) resolutions using deep neural networks (DNNs). While DNNs can learn complex representations that prove challenging for kernel-based methods, they are susceptible to overfitting and the overconfidence of uncertainty estimations. Here, we develop ECG within a GPU-accelerated Deep Kernel Learning (DKL) framework to enable CG QC predictions using range-separated hybrid density functional theory (DFT), obtaining a 107 speedup relative to naive all-atom QC. By treating the predicted electronic properties as random Gaussian Processes, DKL incorporates CG mapping degeneracy by learning the distribution of electronic energies as a function of CG configuration. DKL-ECG accurately reproduces molecular orbital energies from range-separated DFT while facilitating efficient training via active learning using the uncertainties provided by DKL. Further, we show that while active learning algorithms enable efficient sampling of a more diverse configurational space relative to random sampling, all explored query methods exhibit comparable performance for the examined system. We attribute this result to the significant overlap of the feature space and output property distributions across multiple temperatures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Shadow molecular dynamics and atomic cluster expansions for flexible charge models

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Effect of microstructure on fatigue crack propagation in additive manufactured nickel-based superalloy Haynes 282: an experiment and crystal plasticity study

Haynes 282 is a γ' precipitation-strengthened nickel-based superalloy known for its exceptional high-temperature creep resistance and excellent fabricability. Recent advancements in powder-bed fusion-based additive manufacturing (PBF-AM) have enabled the fabrication of Haynes 282 with accurate, site-specific control of the grain orientation and morphology at the microscale. This ability opens up new avenues for microstructure design, to improve the material’s fatigue crack resistance and service life. This paper investigates the fatigue crack growth behavior of hybrid microstructure Haynes 282 fabricated via PBF-AM. Previous experiments revealed a higher crack propagation rate in the coarse columnar-grained microstructural regions when compared against fine-grained areas. Here, a strain gradient crystal plasticity model was adapted to study the fracture-related mechanical fields at the crack tip in the two microstructures. The simulation results showed a consistent influence of grain structure and texture on crack propagation, as was seen in the experiment. The model analysis revealed higher crack propagation driving force along crack direction in the coarse-grained sharply textured microstructure and higher driving force for crack kinking in fine-grained more diffusely textured microstructure, which is ascribed to the combined effect of yield stress, hardening rate, texture and grain morphology. The presented modeling approach will facilitate the development of the AM-based accurate microstructure design by deepening the fundamental study in AM-specific microstructure-properties relations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Grain boundary softening from stress assisted helium cavity coalescence in ultrafine-grained tungsten

The formation of helium cavities in coarse-grained materials produces hardening proportional to the number density and size of the cavities and due to the interaction of dislocations with intragranular helium defects. In nanostructured metals containing a high density of interfacial sinks, preferential cavity formation in the grain boundaries instead produces softening that is often attributed to enhanced interfacial plasticity. Here, employing two grades of ultrafine-grained tungsten, we explore this effect using targeted implantation studies to map cavity evolution as a function of the irradiation conditions and quantify its impact on the mechanical response through nanoindentation. Softening is reported at implantation temperatures above the threshold for preferential grain boundary cavity formation but at a sufficiently low fluence prior to the growth of intragranular cavities. Collective changes in the mean cavity size, density, and morphology beneath a residual impression on an implanted surface indicate that cavity coalescence accompanied the reduction in hardness. Complementary atomistic simulations demonstrate that, in tungsten grain structures exhibiting softening, grain boundary bubble coalescence is driven by stress concentrations that further act to localize strain in the grain boundaries through cooperative deformation processes involving local atomic shuffling and sliding, dislocation emission, and even the nucleation of unstable twinning events.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A classical equation that accounts for observations of non-Arrhenius and cryogenic grain boundary migration

Observations of microstructural coarsening at cryogenic temperatures, as well as numerous simulations of grain boundary motion that show faster migration at low temperature than at high temperature, have been troubling because they do not follow the expected Arrhenius behavior. This work demonstrates that classical equations, that are not simplified, account for all these oddities and demonstrate that non-Arrhenius behavior can emerge from thermally activated processes. According to this classical model, this occurs when the intrinsic barrier energies of the processes become small, allowing activation at cryogenic temperatures. Additional thermal energy then allows the low energy process to proceed in reverse, so increasing temperature only serves to frustrate the forward motion. This classical form is shown to reconcile and describe a variety of diverse grain boundary migration observations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Grain boundary effects in high-temperature liquid-metal dealloying: a multi-phase field study

Abstract A multi-phase field model is employed to study the microstructural evolution of an alloy undergoing liquid dealloying, specifically considering the role of grain boundaries. A semi-implicit time-stepping algorithm using spectral methods is implemented, which enables simulating large 2D and 3D domains over long time scales while still maintaining a realistic interfacial thickness. Simulations reveal a mechanism of coupled grain–boundary migration to maintain equilibrium contact angles with the topologically complex solid–liquid interface, which locally accelerates diffusion-coupled growth of a liquid channel into the precursor. This mechanism asymmetrically disrupts the ligament connectivity of the dealloyed structure in qualitative agreement with published experimental observations. The grain boundary migration-assisted corrosion channels form even for precursors with small amounts of the dissolving alloy species, below the parting limit . The activation of this grain boundary dealloying mechanism depends strongly on grain boundary mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Arbitrarily accurate, nonparametric coarse graining with Markov renewal processes and the Mori–Zwanzig formulation

Stochastic dynamics, such as molecular dynamics, are important in many scientific applications. However, summarizing and analyzing the results of such simulations is often challenging due to the high dimension in which simulations are carried out and, consequently, due to the very large amount of data that are typically generated. Coarse graining is a popular technique for addressing this problem by providing compact and expressive representations. Coarse graining, however, potentially comes at the cost of accuracy, as dynamical information is, in general, lost when projecting the problem in a lower-dimensional space. This article shows how to eliminate coarse-graining error using two key ideas. First, we represent coarse-grained dynamics as a Markov renewal process. Second, we outline a data-driven, non-parametric Mori–Zwanzig approach for computing jump times of the renewal process. Numerical tests on a small protein illustrate the method.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Lunar resources: Oxygen from rocks and soil

The first set of hydrogen reduction experiments to use actual lunar material was recently completed. The sample, 70035, is a coarse-grained vesicular basalt containing 18.46 wt. percent FeO and 12.97 wt. percent TiO2. The mineralogy includes pyroxene, ilmenite, plagioclase, and minor olivine. The sample was crushed to a grain size of less than 500 microns. The crushed basalt was reduced with hydrogen in seven tests at temperatures of 900-1050 C and pressures of 1-10 atm for 30-60 minutes. A capacitance probe, measuring the dew point of the gas stream, was used to follow reaction progress. Experiments were also conducted using a terrestrial basalt similar to some lunar mare samples. Minnesota Lunar Simulant (MLS-1) contains 13.29 wt. percent FeO, 2.96 wt. percent Fe2O3, and 6.56 wt. percent TiO2. The major minerals include plagioclase, pyroxene, olivine, ilmenite, and magnetite. The rock was ground and seived, and experiments were run on the less than 74- and 500-1168-micron fractions. Experiments were also conducted on less than 74-micron powders of olivine, pyroxene, synthetic ilmenite, and TiO2. The terrestrial rock and mineral samples were reduced with flowing hydrogen at 1100 C in a microbalance furnace, with reaction progress monitored by weight loss. Experiments were run at atmospheric pressure for durations of 3-4 hr. Solid samples from both sets of experiments were analyzed by Mossbauer spectroscopy, petrographic microscopy, scanning electron microscopy, tunneling electron microscopy, and x-ray diffraction. Apollo 17 soil 78221 was examined for evidence of natural reduction in the lunar environment. This sample was chosen based on its high maturity level (I sub s/FeO = 93.0). The FeO content is 11.68 wt. percent and the TiO2 content is 3.84 wt. percent. A polished thin section of the 90-150 micron size fraction was analyzed by petrographic microscopy and scanning electron microscopy.

Allen, C. C.↗