Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cell behavior”

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,621 records · Page 90

Impact of supercritical carbon dioxide on the frictional strength of and the transport through thin cracks in shale

Understanding the mechanical and transport behavior of thin (i.e. small aperture) cracks slipping under supercritical carbon dioxide (sc-CO 2 ) conditions is essential to evaluate the integrity of sealing formations with buoyant sc-CO 2 below and the success of waterless fracturing. The two major items of interest in this work are frictional strength and permeability change of the crack. We used a triaxial cell that permits in situ visualization to conduct and monitor slippage along the faces of narrow cracks subjected to triaxial stresses. Such cracks are analogs to small geological faults. We tested carbonate-rich, 1-inch diameter Wolfcamp shale samples that are saw cut 30° to vertical to create a thin crack. Friction coefficients ranged from about 0.6 to 0.8 consistent with expectations for brittle rocks. The sc-CO 2 generally did not alter friction coefficient over the time scale of experiments. From a transport perspective, saturating cracks with sc-CO 2 substantially decreased permeability of the crack by 26%–52%, while slip resulted in a variety of permeability responses. Overall, the combined impact of sc-CO 2 saturation and slip reduced fault permeability for all tests. Our observations support the notion that the sealing capacity of some caprocks improves when saturated with sc-CO 2 and that some slip of small fractures is not necessarily detrimental to caprock integrity.

58 GEOSCIENCES↗

Analysis of tumor-immune functional responses in a mathematical model of neoantigen cancer vaccines

Cancer neoantigen vaccines have emerged as a promising approach to stimulating the immune system to fight cancer. We propose a simple model including key elements of cancer-immune interactions and conduct a phase plane analysis to understand the immunological mechanisms of cancer neoantigen vaccines. Analytical results are obtained for two widely used functional forms that represent the killing rate of tumor cells by immune cells: the law of mass action (LMA) and the dePillis-Radunskaya Law (LPR). Using the LMA, our results reveal that a slowly growing tumor can escape the immune surveillance and that there is a unique periodic solution. The LPR offers richer dynamics, in which tumor elimination and uncontrolled tumor growth are both present. We show that tumor elimination requires sufficient number of initial activated T cells in relationship to the malignant cells, which lends support to using the neoantigen cancer vaccine as an adjuvant therapy after the primary tumor is surgically removed or treated using radiotherapy. We also derive a sufficient condition for uncontrolled tumor growth under the assumption of the LPR. Here, the juxtaposition of analyses with these two different choices for the killing rate function highlights their importance on model behavior and biological implications, by which we hope to spur further theoretical and experimental work to understand mechanisms underlying different functional forms for the killing rate.

60 APPLIED LIFE SCIENCES↗

Carbon nanotube (CNT) reinforced 316L stainless steel composites made by laser powder bed fusion: Microstructure and wear response

Carbon nanotube (CNT) reinforced 316L stainless steel (SS) composites were fabricated by laser powder bed fusion (LPBF) additive manufacturing (AM). This study focuses on microstructure evolution, hardness, and wear behavior. The wear behavior for AM 316L SS and 1 wt% CNT/316L SS was studied by a dry sliding test under normal loads of 10, 15, and 25 N. The addition of CNT to 316L SS led to the change of solidification mode from cellular/columnar to dendritic growth, producing grain refinement. Nanoscale dendrites were observed. Although TEM characterization confirmed the existence of CNT in cell interior, most CNT were damaged, resulting in carbon material segregation along cellular/dendritic boundaries. The wear rate was found dramatically decreased due to the improved hardness. Compared to the adhesive and abrasive wear in AM 316L SS, the wear mechanism of CNT/316L composite was adhesion and oxidative wear. Furthermore, enhanced cellular/dendritic boundary strengthening by carbon segregation was the main strengthening mechanism of CNT/316L composite.

36 MATERIALS SCIENCE↗

Emerging Potassium Metal Anodes: Perspectives on Control of the Electrochemical Interfaces

Potassium metal acts as the anode in emerging potassium metal batteries (KMBs). It further serves as the counter-electrode for potassium ion battery (KIB) half-cells, with its reliable performance being critical for assessing the working electrode material. This first-of-its-kind critical review focuses on the dual challenge of controlling the potassium metal-substrate and the potassium metal-electrolyte interface so as to prevent dendrites. The discussion begins with a comparison of the physical and chemical properties of K metal anodes versus the much oft studied Li and Na metal anodes. Based on established descriptions for root causes of dendrites, the problem should be less severe for K than for Li or Na, while in fact the opposite is observed. The key reason that the K metal surface rapidly becomes dendritic in common electrolytes is its unstable solid electrolyte interphase (SEI). An unstable SEI layer is defined as being non self-passivating. No SEI is perfectly stable during cycling, and all SEI structures are heterogenous both vertically and horizontally relative to the electrolyte interface. The difference between a "stable" and an "unstable" SEI may be viewed as the relative degree to which during cycling it thickens and becomes further heterogeneous. The unstable SEI on K metal leads to a number of interrelated problems, such as low cycling Coulombic efficiency (CE), a severe impedance rise, large overpotentials, and possibly electrical shorting; all of which have been reported to occur as early as in the first ten plating/stripping cycles. Many of the traditional "interface fixes" employed for Li and Na metal anodes, such as various artificial SEIs, surface membranes, barrier layers, the secondary separators, etc. have not been attempted or optimized for the case of K. This is an important area for further exploration, with an understanding that success may come harder than with Li due to K-based SEI reactivity with both ether and ester solvents. The second critical problem with K metal anodes is that they do not thermally or electrochemically wet a standard (untreated) Cu foil current collector. Published experimental and modeling research directly highlights the weak bonding between the K atoms and a Cu surface. Existing surface treatment approaches that achieve improved K wetting are discussed, along with the general design rules for future studies. Also discussed are geometry-based methods to tune nucleation, as well dual approaches where nucleation and SEI structure are tuned through complementary schemes to achieve extended half-cell and full battery stability. We hypothesize that K metal never achieves a planar wetting morphology even at cycle one making the dendrites "baked-in". We propose that classical thin films growth models, Frank van der Merwe (F-M), Volmer Weber (V-W), and Stranski-Krastanov (S-K) can be employed to describe early stage plating behavior. It is demonstrated that island-like V-W growth is the applicable description for the natural plating behavior of K on pristine Cu. Moving forward, there are three inter-related thrusts to be pursued: First, K salt - based electrolyte formulations have to mature and become further tailored to handle the increased reactivity of a metal rather than an ion anode. Second, the K-based SEI structure needs to be further understood and ultimately tuned to be less reactive. Third, the energetics of the K metal - current collector interface must be controlled to promote planar wetting/dewetting throughout cycling.

25 ENERGY STORAGE↗

Importance of Site Diversity and Connectivity in Electrochemical CO Reduction on Cu

Electrochemical CO 2 reduction on Cu is a promising approach to produce value-added chemicals using renewable feedstocks, yet various Cu preparations have led to differences in activity and selectivity toward single and multicarbon products. Here, we find, surprisingly, that the effective catalytic activity toward ethylene improves when there is a larger fraction of less active sites acting as reservoirs of *CO on the surface of Cu nanoparticle electrocatalysts. In an adaptation of chemical transient kinetics to electrocatalysis, we measure the dynamic response of a gas diffusion electrode (GDE) cell when the feed gas is abruptly switched between Ar (inert) and CO. When switching from Ar to CO, CO reduction (COR) begins promptly, but when switching from CO to Ar, COR can be maintained for several seconds (delay time) despite the absence of the CO reactant in the gas phase. A three-site microkinetic model captures the observed dynamic behavior and shows that Cu catalysts exhibiting delay times have a less active *CO reservoir that exhibits fast diffusion to active sites. The observed delay times and the estimated *CO reservoir sizes are affected by catalyst preparation, applied potential, and microenvironment (electrolyte cation identity, electrolyte pH, and CO partial pressure). Notably, we estimate that the *CO reservoir surface coverage can be as high as 88 ± 7% on oxide-derived Cu (OD-Cu) at high overpotentials (–1.52 V vs SHE) and this increases in reservoir coverage coincide with increased turnover frequencies to ethylene. We also estimate that *CO can travel substantial distances (up to 10s of nm) prior to desorption or reaction. It appears that active C–C coupling sites by themselves do not control selectivity to C 2+ products in electrochemical COR; the supply of CO to those sites is also a crucial factor. More generally, the overall activity of Cu electrocatalysts cannot be approximated from linear combinations of individual site activities. Future designs must consider the diversity of the catalyst network and account for intersite transportation pathways.

30 DIRECT ENERGY CONVERSION↗

Polysulfide entrapment and retardation in gel electrolyte Li–S batteries: experiments and modeling

Gelling the electrolyte in lithium–sulfur (Li–S) batteries provides advantages in terms of safety and mitigating polysulfide shuttling leading to reduced side reactions and dendrite growth on the anode side. In the current study, organic and ceramic crosslinkers have been used to improve the entrapment of polysulfides and prevent them from shuttling to the anode side. However, the performance of gel electrolytes is still far from being on par with liquid electrolytes in terms of raw capacity. The discharge curves exhibit an extended first plateau for gel electrolyte batteries compared to liquid electrolytes, showing that polysulfides are better trapped in the gel system. Yet, in the conversion of soluble polysulfides into insoluble polysulfides, the second plateau was delayed. Qualitative and quantitative analysis of the polysulfide diffusion in the liquid electrolyte and gel electrolyte systems were also conducted to predict the relative diffusivity ratio between the liquid electrolyte and gel electrolyte for simulation. To probe what happens at the end of the discharge, where the mass transport in the gel electrolyte appears to be slow, a modified discharge procedure was used where the cell is first paused after discharge and then discharged a second time before recharging normally. Coin cell tests for two gel electrolytes in addition to liquid electrolyte show up to 7% capacity recovery after the pause, which indicates that the diffusion of species during the pause is responsible for the failure recovery. Additionally, the recovery is higher for gel electrolytes compared to liquid, for higher C-rates, and for longer pause times. To better understand this behavior, a Li–S numerical model with mass transport limited reactions was used to examine the polysulfide diffusion coefficients expected in different electrolyte systems. The model is able to reproduce the trends seen in experiments and yields a higher recovery for smaller Li 2 S 4 diffusion coefficients, suggesting that insufficient Li 2 S 4 mass transport is responsible for failure at the end of discharge. By identifying the effect that gel electrolytes have on the state of the battery during discharge using both experimental and simulation methods, this work is a step towards leveraging the advantages of gel electrolytes and mitigating their weaknesses.

25 ENERGY STORAGE↗

Temporal deposition of copper and zinc in the sediments of metal removal constructed wetlands

The objective of this study was to explore the effects of time, seasons, and total carbon (TC) on Copper (Cu) and Zinc (Zn) deposition in the surface sediments. This study was performed at the H-02 constructed wetland on the Savannah River Site (Aiken, SC, USA). Covering both warm (April-September) and cool (October-March) seasons, several sediment cores were collected twice a year from the H-02 constructed wetland cells from 2007 to 2013. Total concentrations of Cu and Zn were measured in the sediments. Concentrations of Cu and Zn (mean ± standard deviation) in the surface sediments over 7 years of operation increased from 6.0 ± 2.8 and 14.6 ± 4.5 mg kg -1 to 139.6 ± 87.7 and 279.3 ± 202.9 mg kg -1 dry weight, respectively. The linear regression model explained the behavior and the variability of Cu deposition in the sediments. On the other hand, using the generalized least squares extension with the linear regression model allowed for unequal variance and thus produced a model that explained the variance properly, and as a result, was more successful in explaining the pattern of Zn deposition. Total carbon significantly affected both Cu ( p = 0.047) and Zn ( p < 0.001). Time effect on Cu deposition was statistically significant ( p = 0.013), whereas Zn was significantly affected by the season ( p = 0.009).

54 ENVIRONMENTAL SCIENCES↗

Micro-mechanics of micro-composites

The Structural Dynamics branch at NASA LaRC is working on developing an active passive mount system for vibration control. Toward this end a system utilizing piezoelectric actuators is currently being utilized. There are limitations to the current system related to space applications under which it is desired to eliminate deformations in the actuators associated with thermal effects. In addition, a material that is readily formable into complex shapes and whose mechanical properties can be optimized with regards to vibration control would be highly desirable. Microcomposite material are currently under study to service these needs. Microcomposite materials are essentially materials in which particles on the scale of microns are bound together with a polyimide (LaRC Si) that has been developed at LaRC. In particular a micro-composite consisting of LaRC Si binder and piezoelectric ceramic particles shows promise in satisfying the needs of the active passive mount project. The LaRC/ Si microcomposite has a unique combination of piezoelectric properties combined with a near zero coefficient of thermal expansion and easy machinability. The goal of this ASEE project is to develop techniques to analytically determine important material properties necessary to characterize the dynamic properties of actuators and mounts made from the LaRC Si / ceramic microcomposite. In particular, a generalized method of cells micromechanics originally developed at NASA Lewis is employed to analyze the microstructural geometry of the microcomposites and predict the overall mechanical properties of the material. A testing program has been established to evaluate and refine the GMC approach to these materials. In addition, a theory of mixtures analysis is being developed that utilizes the GMC micromechanics information to analyze complex behavior of the microcomposite material which has a near zero CTE.

Donovan, Richard P.↗

Membrane peptides and their role in protobiological evolution

How simple membrane peptides performed such essential protocellular functions as transport of ions and organic matter across membranes separating the interior of the cell from the environment, capture and utilization of energy, and transduction of environmental signals, is a key question in protobiological evolution. On the basis of detailed, molecular-level computer simulations we explain how these peptides fold at water-membrane interfaces, insert into membranes, self-assemble into higher-order structures and acquire functions. We have investigated the interfacial behavior and folding of several peptides built of leucine and glutamine residues and have demonstrated that many of them tend to adopt ordered structures. Further, we have studied the insertion of an alpha-helical peptide containing leucine (L) and serine (S) of the form (LSLLLSL)3 into a model membrane. The transmembrane state is metastable, and approximately 15 kcal mol(-1) is required to insert the peptide into the membrane. Investigations of dimers formed by (LSLLLSL)3 and glycophorin A demonstrate how the favorable free energy of helix association can offset the unfavorable free energy of insertion, leading to self-assembly of peptide helices in the membrane. An example of a self-assembled structure is the tetrameric transmembrane pore of the influenza virus M2 protein, which is an efficient and selective voltage-gated proton channel. Our simulations explain the gating mechanism and provide guidelines how to re-engineer the channel to act as a simple proton pump. In general, emergence of integral membrane proteins appears to be quite feasible and may be easier to envision than the emergence of water-soluble proteins.

Models, Chemical↗

Insights into the Enhanced Reversibility of Graphite Anode Upon Fast Charging Through Li Reservoir

Increasing the charging rate and reducing the charging time for Li-ion batteries are crucial to realize the mainstream of electric vehicles. However, it is formidable to avoid the Li plating on graphite anode upon fast charging. Despite the tremendous progress in Li detection techniques, the fundamental mechanism of Li plating and its chemical/electrochemical responses upon cycling still remains elusive. Herein, we present a comprehensive electrochemical method to investigate the fast charging behavior of graphite electrode. A detailed analysis is directed toward understanding the changes in phase, composition, and morphology of the fast-charged graphite. By applying a resting process, we scrutinize the further reactions of the plated Li, which readily transforms into irreversible (dead) Li. We further develop a modified graphite electrode with a thin Ag coating as the Li reservoir. The plated Li can be "absorbed" by the Ag layer to form the Li-Ag solid solution that suppresses the formation of dead Li and provides structural stability, thus promoting the further lithiation of graphite and enhancing the reversibility. Here this work not only provides additional insights into the fast charging behavior of graphite electrode but also demonstrates a potential strategy to improve the fast charging performance of graphite anode.

25 ENERGY STORAGE↗

In Situ Inference for Earth System Predictability

An understanding of future evolution in precipitation extremes is critical to numerous DOE mission questions. Extreme events are by nature short time-scale events that are difficult to diagnose in available model data. Accurate modeling of extreme events necessarily requires high spatial resolution at the storm scale locally. However, the environment in which storms grow is dependent on global, remote, processes. These complex spatiotemporal relationships are impossible to diagnose at resolutions required to accurately model storms responsible for extreme precipitation. At exascale, climate simulations will produce results at fine enough resolution to investigate these relationships. However, the resulting data from these simulations will be far too large to save for post-simulation analysis. We advocate for fitting statistical models inside the simulations as they run, a context known as in situ, which will facilitate scientific investigations using the full fine-scale data stream. Figure 1 shows an example of the type of model we could consider, a Bayesian hierarchical spatial regression model. Precipitation extremes at each grid cell are modeled using extreme value distributions. Since extremes are rare, fitting models to individual grid cells can result in high variance and poor estimates. Instead, the model can be made more robust by smoothing the parameters of the extreme value model across space. Additionally, the parameters themselves can be functionally linked to other variables elsewhere in the simulation. Thus, we can use the fine-scale data to build more robust models for extremes that link extreme behavior to other climate patterns.

54 ENVIRONMENTAL SCIENCES↗

In Situ Inference for Earth System Predictability

Focal Area: Focal Area 3: Insight gleaned from complex simulated data using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI. Science Challenge: An understanding of future evolution in precipitation extremes is critical to numerous DOE mission questions. Extreme events are by nature short time-scale events that are difficult to diagnose in available model data. Accurate modeling of extreme events necessarily requires high spatial resolution at the storm scale locally. However, the environment in which storms grow is dependent on global, remote, processes. These complex spatiotemporal relationships are impossible to diagnose at resolutions required to accurately model storms responsible for extreme precipitation. At exascale, climate simulations will produce results at fine enough resolution to investigate these relationships. However, the resulting data from these simulations will be far too large to save for post-simulation analysis. We advocate for fitting statistical models inside the simulations as they run, a context known as in situ, which will facilitate scientific investigations using the full fine-scale data stream. Figure 1 shows an example of the type of model we could consider, a Bayesian hierarchical spatial regression model. Precipitation extremes at each grid cell are modeled using extreme value distributions. Since extremes are rare, fitting models to individual grid cells can result in high variance and poor estimates. Instead, the model can be made more robust by smoothing the parameters of the extreme value model across space. Additionally, the parameters themselves can be functionally linked to other variables elsewhere in the simulation. Thus, we can use the fine-scale data to build more robust models for extremes that link extreme behavior to other climate patterns.

54 ENVIRONMENTAL SCIENCES↗

Development and test fuel cell powered on-site integrated total energy systems. Phase 3: Full-scale power plant development

The on-site system application analysis is summarized. Preparations were completed for the first test of a full-sized single cell. Emphasis of the methanol fuel processor development program shifted toward the use of commercial shell-and-tube heat exchangers. An improved method for predicting the carbon-monoxide tolerance of anode catalysts is described. Other stack support areas reported include improved ABA bipolar plate bonding technology, improved electrical measurement techniques for specification-testing of stack components, and anodic corrosion behavior of carbon materials.

Kaufman, A.↗

On Io, all that flickers is not cold

The 4.8 micron flux from the Galilean satellite Io was shown to fluctuate with an amplitude of approximately 10 percent on time scales of every 28 seconds, 40 minutes, night-to-night, and perhaps year-to-year. Such behavior was found to be the result of random fluctuations for which the mean square fluctuation varies inversely with frequency for constant bandwidth measurement. The theory developed for thermionic emission from barium oxide cathodes in vacuum tubes might be applicable to this situation. If so, the fluctuations in the flux from Io's volcanoes may be caused by diffusion of hot convective cells onto the surface of Io. Long term fluctuations may furnish a means by which the Io volcanism can shut down and conserve energy. Thus the discrepancy that exists between measurements of the current heat flow from Io and calculations of tidal dissipation may be resolved. Tests for rapid flickering at 10 microns showed no fluctuations greater than one percent. This agrees with the prediction of the flow model theory in which the 10 microns volcanic thermal emission arises from cooling of old flows.

Sinton, W. M.↗

Studies of Model Immiscible Systems

The objectives are to use model transparent monotectics to obtain fundamental information applicable to two-phase systems in general, to apply this understanding to materials of interest in the Microgravity Science and Applications program, and to interpret results of flight experimental involving monotectic alloys. A number of model immiscible systems are in use to study various aspects of two-phase behavior within the miscibility gap and during solidification. Particle growth, coalescence and particle motions are under investigation using a holographic microscopy system. The system is capable of working with particle densities up to 10 to the 7th power particles/cubic centimeters through a 100 micron depth and can resolve particles of the order of 2 to 3 micron in diameter throughout the entire cell volume. Particle size, distribution changes with respect to time and temperature are observable from sequential holograms. Initial experiments using diethylene glycol/ethyl salicylate (DEG/ES) have demonstrated the usefulness of the technique. The thermal system controls temperature to at least plus or minus 0.001 K over the course of an experiment. A time-lapse film, made from holograms, of a succinonitrile/water solution shows particle size and number distribution changes with time under isothermal conditions. The observations are consistent with Ostwald ripening theory.

Frazier, D. O.↗

Laboratory and theoretical models of planetary-scale instabilities and waves

Research work is proceeding in theoretical, numerical, and experimental geophysical fluid dynamics leading up to a reflight of the GFFC (Geophysical Fluid Flow Cell Experiment) on USML-2. The work is intended not only to generate ideas for future space experiments, but to provide fundamental results concerned with nonlinear and chaotic properties of thermal convection and baroclinic waves in terrestrial and planetary atmospheres. The major efforts are focussed on thermal convection in a rapidly rotating annulus relevant to Jovian atmospheric dynamics, and on the chaotic behavior of baroclinic waves relevant to the Earth's atmosphere. The approach, in preparation for USML-2, is primarily theoretical and numerical. Mechanistic process models are solved numerically in order to identify physical mechanisms that may be observed in the GFFC, and which are important in real geophysical applications. The results from numerical simulations of geophysical fluid flow (subject to rotation and stratification) are compared with previous GFFC experiments on Spacelab-3 and with existing and proposed terrestrial laboratory experiments of various types. Pattern recognition algorithms have been employed to generate low-dimensional descriptions of the highly nonlinear and turbulent numerical simulations. Such empirically truncated descriptions provide for simplified but robust physical interpretations of the dynamics, as well as yielding highly efficient computations of these chaotic flows.

Hart, John E.↗

Intercalation of Lithium in Pitch-Based Graphitized Carbon Fibers Chemically Modified by Fluorine: Soft Carbon With or Without an Oxide Surface

The effects of carbon structure and surface oxygen on the carbon's performance as the anode in lithium-ion battery were studied. Two carbon materials were used for the electrochemical tests: soft carbon made from defluorination of graphite fluoride, and the carbon precursor from which the graphite fluoride was made. In this research the precursor was graphitized carbon fiber P-100. It was first fluorinated to form CF(0.68), then defluorinated slowly at 350 to 450 C in bromoform, and finally heated in 1000 C nitrogen before exposed to room temperature air, producing disordered soft carbon having basic surface oxides. This process caused very little carbon loss. The electrochemical test involved cycles of lithium intercalation and deintercalation using C/saturated LiI-50/50 (vol %) EC and DMC/Li half cell. The cycling test had four major results. (1) The presence of a basic oxide surface may prevent solvent from entering the carbon structure and therefore prolong the carbon's cycle life for lithium intercalation-deintercalation. (2) The disordered soft carbon can store lithium through two different mechanisms. One of them is lithium intercalation. which gives the disordered carbon an electrochemical behavior similar to its more ordered graphitic precursor. The other is unknown in its chemistry, but is responsible for the high-N,oltage portion (less than 0.3V) of the charge-discharge curve. (3) Under certain conditions, the disordered carbon can store more lithium than its precursor. (4) These sample and its precursor can intercalate at 200 mA/g. and deintercalate at a rate of 2000 mA/g without significant capacity loss.

Hung, Ching-Chen↗

On the Role of Friction and Particle Size Distribution in Granular Packings

Packing of particles in a disordered arrangement has tremendous significance in both condensed matter physics and engineering applications. The last three decades have seen remarkable progress in our understanding of the physics of granular packings that has been largely facilitated by a rapid growth in the power of modern computers. Although granular packings are ubiquitous in diverse natural settings, from clogging of powders in hoppers to the crowding of living cells, a significant motivation for modeling granular packings has emerged from a proposal that the transition from a fluid-like state to a solid-like state of a granular material upon increasing volume fraction, called jamming, is intimately related to the origins of glass transition in thermal systems. As such, a majority of modeling efforts have focused on the jamming behavior of an idealized granular material: frictionless, monodisperse sphere. While such studies have illuminated the rich physics of jamming, granular materials in nature and engineering practice are rarely frictionless or monodisperse. The analogous research on the packing of these ‘real-world’ granular materials is still not fully developed. Besides requiring the exploration of a huge parameter space, three key computational considerations have inhibited their modeling: (i) traditional computational methods are not adept at simulating mechanically-stable packings of frictional particles near the jamming transition; (ii) standard algorithms of contact detection in discrete element methods are impractical to simulate granular packings with a wide distribution of particle sizes; (iii) a lack of well-established contact mechanics models of friction that can accurately reproduce experimental data. This chapter will review the latest computational advances to simulate the jamming of size-dispersed frictional particles, and describe the rich microstructural diversity that emerges in their packings.

granular↗