Powders under compression form mesostructures of particle agglomerations in response to both inter- and intra-particle forces. The ability to computationally predict the resulting mesostructures with reasonable accuracy requires models that capture the distributions associated with particle size and shape, contact forces, and mechanical response during deformation and fracture. The following report presents experimental data obtained for the purpose of validating emerging mesostructures simulated by discrete element method and peridynamic approaches. A custom compression apparatus, suitable for integration with our micro-computed tomography (micro-CT) system, was used to collect 3-D scans of a bulk powder at discrete steps of increasing compression. Details of the apparatus and the microcrystalline cellulose particles, with a nearly spherical shape and mean particle size, are presented. Comparative simulations were performed with an initial arrangement of particles and particle shapes directly extracted from the validation experiment. The experimental volumetric reconstruction was segmented to extract the relative positions and shapes of individual particles in the ensemble, including internal voids in the case of the microcrystalline cellulose particles. These computationally determined particles were then compressed within the computational domain and the evolving mesostructures compared directly to those in the validation experiment. The ability of the computational models to simulate the experimental mesostructures and particle behavior at increasing compression is discussed.
The peridynamic theory of solid mechanics is applied to modeling the deformation and fracture of micrometer-sized particles made of organic crystalline material. A new peridynamic material model is proposed to reproduce the elastic–plastic response, creep, and fracture that are observed in experiments. The model is implemented in a three-dimensional, meshless Lagrangian simulation code. In the small deformation, elastic regime, the model agrees well with classical Hertzian contact analysis for a sphere compressed between rigid plates. Under higher load, material and geometrical nonlinearity is predicted, leading to fracture. The material parameters for the energetic material CL-20 are evaluated from nanoindentation test data on the cyclic compression and failure of micrometer-sized grains.
This Laboratory Directed Research and Development project developed and applied closely coupled experimental and computational tools to investigate powder compaction across multiple length scales. The primary motivation for this work is to provide connections between powder feedstock characteristics, processing conditions, and powder pellet properties in the context of powder-based energetic components manufacturing. We have focused our efforts on multicrystalline cellulose, a molecular crystalline surrogate material that is mechanically similar to several energetic materials of interest, but provides several advantages for fundamental investigations. We report extensive experimental characterization ranging in length scale from nanometers to macroscopic, bulk behavior. Experiments included nanoindentation of well-controlled, micron-scale pillar geometries milled into the surface of individual particles, single-particle crushing experiments, in-situ optical and computed tomography imaging of the compaction of multiple particles in different geometries, and bulk powder compaction. In order to capture the large plastic deformation and fracture of particles in computational models, we have advanced two distinct meshfree Lagrangian simulation techniques: 1.) bonded particle methods, which extend existing discrete element method capabilities in the Sandia-developed , open-source LAMMPS code to capture particle deformation and fracture and 2.) extensions of peridynamics for application to mesoscale powder compaction, including a novel material model that includes plasticity and creep. We have demonstrated both methods for simulations of single-particle crushing as well as mesoscale multi-particle compaction, with favorable comparisons to experimental data. We have used small-scale, mechanical characterization data to inform material models, and in-situ imaging of mesoscale particle structures to provide initial conditions for simulations. Both mesostructure porosity characteristics and overall stress-strain behavior were found to be in good agreement between simulations and experiments. We have thus demonstrated a novel multi-scale, closely coupled experimental and computational approach to the study of powder compaction. This enables a wide range of possible investigations into feedstock-process-structure relationships in powder-based materials, with immediate applications in energetic component manufacturing, as well as other particle-based components and processes.
Microstructures and corresponding properties of compacted powders ultimately depend on the mechanical response of individual particles. In principle, computational simulations can predict the results of powder compaction processes, but the selection of appropriate models for both particle–particle interactions and particle deformations across all relevant length scales remain nontrivial tasks, especially in material systems lacking detailed mechanical property information. The work presented here addresses these issues by conducting uniaxial compressions in situ inside of a scanning electron microscope to characterize the mechanical response of individual micron-sized particles of a molecular crystal, hexanitrohexaazaisowurtzitane (CL-20). This experimental approach enabled the collection of quantitative force and displacement data alongside simultaneous imaging to capture morphology changes. The results reveal information about elastic deformation, yield, plastic deformation, creep, and fracture phenomena. Accordingly, this work demonstrates a generalizable approach for assessing the mechanical response of individual micron-sized molecular crystal particles and utilizing those responses in particle-level models. Graphic abstract: [Figure not available: see fulltext.].
Ignition and material response properties of aluminumized HMX heterogeneous explosive mixtures were explored in a series of planar impact experiments performed over multiple years. This work expands on previous work studying material response to impact in single-component HMX granular materials. The addition of nanometric aluminum is shown to affect the ignition sensitivity and growth to reaction from impact. The gas gun test results are presented here varying parameters of particle size, shock strength, and aluminum mass fraction.
Two key mechanical processes exist in the formation of powder compacts. These include the complex kinematics of particle rearrangement as the powder is densified and particle deformation leading to mechanical failure and fragmentation. Experiments measuring the time varying forces across a densifying powder bed have been performed in powders of microcrystalline cellulose with mean particle sizes between 0.4 and 1.2 mm. In these experiments, diagnostics measured the applied and transmitted loads and the bulk powder density. Any insight into the particle behavior must be inferred from deviations in the smoothly increasing stress-density compaction relationship. By incorporating a window in the compaction die body, simultaneous images of particle rearrangement and fracture at the confining window are captured. The images are post-processed in MATLAB® to track individual particle motion during compression. Complimentary discrete element method (DEM) simulations are presented and compared to experiment. The comparison provides insight into applying DEM methods for simulating large or permanent particle deformation and suggests areas for future study.
Particle characteristics can drastically influence the process-structure-property-performance aspects of granular materials in compression. We aim to computationally simulate the mechanical processes of stress redistribution in compacts including the kinematics of particle rearrangement during densification and particle deformation leading to fragmentation. Confined compression experiments are conducted with three sets of commercial microcrystalline cellulose particles nearly spherical in shape with different mean particle size. Experimentally measured compression curves from tall powder columns are fitted with the Kenkre et al. (J. of American Chemical Society, Vol. 79, No. 12) model. This model provides a basis to derive several common two-parameter literature models and as a framework to incorporate statistical representations of critical particle behaviors. We focus on the low-stress compression data and the model comparisons typically not discussed in the literature. Additional single particle compressions report fracture strength with particle size for comparison to the apparent particle strength extracted from bulk compression data.
A novel experimental methodology is presented to study the deviatoric response of powders in shock regimes. The powders are confined to a cylindrical wedge volume, and a projectile-driven shock wave with a sinusoidally varying front propagates through the powder. The perturbed shock wave exhibits a damping behavior due to irreversible processes of viscosity and strength (deviatoric) of the powder with propagation through increasing powder thicknesses. The inclined surface of the wedge is polished and coated to establish a diffuse surface suitable for reflecting incident laser light into a high-speed camera imaging at 5 MHz. Images of the contrast loss upon shock wave arrival at the observation surface are post-processed for qualitative and quantitative information. New data of shock damping behavior with parameters of perturbation wavelength and initial shock strength are presented for powders of copper, tantalum, and tungsten carbide as well as their mixtures. We present the first full-field images showing additional spatial disturbances on the perturbed shock front that appear dependent on particle material and morphology.
Surrogate laboratory experiments were conducted to understand how cavities surrounding an explosive event in a bounded material influence the transfer of energy to the surface of the material. Exploding bridgewires were detonated in a cavity that had been created in an optically accessible polymer material. The shock propagation was investigated through five paths, three through spherical cavities of varying size filled with air, one without a cavity and one in air. Shock propagation was imaged using the schlieren technique. The path effects were determined by measuring energy at the polymer surface using a pressure sensor. All cavities were successful at reducing the peak pressure nearly 100% from the case without a cavity, a reduction close to air. The degree of reduction was not uniform over all frequencies. Overall, our findings suggest that cavities influence both the amount and frequency content of the energy transfer from source to surface.
A multi-point microwave interferometer (MPMI) concept is presented for non-invasively monitoring the internal transit of a shock, detonation, or reaction front in energetic media. The concept utilizes an electro-optic (EO) crystal to impart a timevarying phase lag onto a laser with a microwave signal. Polarization optics convert this phase lag into an amplitude modulation. A heterodyne interferometer compares the modulated laser beam to a constant reference. This enables the detection of changes in the modulating microwave frequency generated by the motion of the measurement surface. The design is scalable and makes use of the established construction and analysis methods employed in photonic Doppler velocimetry (PDV). The technical challenges associated with the concept are the frequency stability of the lasers, the amount of light return after EO modulation, and the frequency uncertainty of fast Fourier transform (FFT) methods.
The combustion of molten metals is an important area of study with applications ranging from solid aluminized rocket propellants to fireworks displays. This work uses digital in-line holography (DIH) to experimentally quantify the three-dimensional position, size, and velocity of aluminum particles during combustion of ammonium perchlorate (AP) based solid-rocket propellants. In addition, spatially resolved particle temperatures are simultaneously measured using two-color imaging pyrometry. To allow for fast characterization of the properties of tens of thousands of particles, automated data processing routines are proposed. Using these methods, statistics from aluminum particles with diameters ranging from 15 to 900 µm are collected at an ambient pressure of 83 kPa. In the first set of DIH experiments, increasing initial propellant temperature is shown to enhance the agglomeration of nascent aluminum at the burning surface, resulting in ejection of large molten aluminum particles into the exhaust plume. The resulting particle number and volume distributions are quantified. In the second set of simultaneous DIH and pyrometry experiments, particle size and velocity relationships as well as temperature statistics are explored. The average measured temperatures are found to be 2640 ± 282 K, which compares well with previous estimates of the range of particle and gas-phase temperatures. The novel methods proposed here represent new capabilities for simultaneous quantification of the joint size, velocity, and temperature statistics during the combustion of molten metal particles. The proposed techniques are expected to be useful for detailed performance assessment of metalized solid-rocket propellants.