Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion limitation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Strongly Anharmonic Phonons and Their Role in Superionic Diffusion and Ultralow Thermal Conductivity of Cu 7 PSe 6

We report the quest for advanced superionic materials requires understanding their complex atomic dynamics, but detailed studies of the interplay between lattice vibrations and ionic diffusion remain scarce. Here inelastic and quasielastic neutron scattering measurements in the superionic argyrodite Cu 7 PSe 6 are reported, combined with molecular dynamics (MD) based on ab initio and machine-learned potentials, providing critical insights into the atomistic mechanisms underlying fast ion conduction. The results reveal how long-range Cu diffusion is limited by intercluster hopping, controlled by selective anharmonic phonons of the crystalline framework. Further, the Green–Kubo simulations reproduce the ultralow lattice thermal conductivity and identify contributions from mobile ions, phonons, and their cross-correlations. The mode resolved analysis shows that the thermal conductivity is dominated by low-energy acoustic phonon modes of the overall crystal framework. The analysis of mode-resolved spectral functions further show that vibrational modes with significant Cu contributions are strongly damped, corresponding to the breakdown of associated phonon quasiparticles. These results highlight the importance of strongly anharmonic effects in superionic systems, in which the traditional quasiharmonic phonon picture is insufficient, and pave the way toward combining machine-learning accelerated simulations with neutron scattering experiments to rationalize the complex atomic dynamics underlying ionic and thermal transport.

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

Explicit Monotone Stable Super-Time-stepping Methods for Finite Time Singularities

We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$⁡(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

Near-Limit Flamelet Phenomena in Buoyant Low Stretch Diffusion Flames Beneath a Solid Fuel

A unique near-limit low stretch multidimensional stable flamelet phenomena has been observed for the first time which extends the material flammability limit beyond the one-dimensional low stretch flammability limit to lower burning rates and higher relative heat losses than is possible with uniform flame coverage. During low stretch experiments burning the underside of very large radii (greater than or = 75 cm stretch rate less than or = 3/s) cylindrical cast PMMA samples, multidimensional flamelets were observed, in contrast with a one-dimensional flame that was found to blanket the surface for smaller radii samples ( higher stretch rate). Flamelets were observed by decreasing the stretch rate or by increasing the conductive heat loss from the flame. Flamelets are defined as flames that cover only part of the burning sample at any given time, but persist for many minutes. Flamelet phenomena is viewed as the flame's method of enhancing oxygen flow to the flame, through oxygen transport into the edges of the flamelet. Flamelets form as heat losses (surface radiation and solid-phase conduction) become large relative to the weakened heat release of the low stretch flame. While heat loss rates remain fairly constant, the limiting factor in the heat release of the flame is hypothesized to be the oxygen transport to the flame in this low stretch (low convective) environment. Flamelet extinction is frequently caused by encroachment of an adjacent flamelet. Large-scale whole-body flamelet oscillations at 1.2 - 1.95 Hz are noted prior to extinction of a flamelet. This oscillation is believed to be due a repeated process of excess fuel leakage through the dark channels between the flamelets, fuel premixing with slow incoming oxidizer, and subsequent rapid flame spread and retreat of the flamelet through the premixed layer. The oscillation frequency is driven by gas-phase diffusive time scales.

Olson, S. L.↗
DOE OSTI2023

Product selectivity and mass transport in levulinic acid transfer hydrogenation by monolithic MIL-100, MIL-88B and ZIF-8@Pd MOFs

The diffusion processes between adsorbent and adsorbate naturally play a significant role in the efficiency and selectivity of the heterogenous catalytic process. This paper considers the importance of diffusion processes in the transfer hydrogenation reaction of levulinic acid to γ-valerolactone by MIL-88B, MIL-100, and ZIF8@Pd monolithic catalysts. Over a period of five catalytic cycles, it was shown that the Fe-based catalysts can achieve similar conversions to the ZIF-8 supported Pd, with the only current limitation being the lower aqueous stability of these MOFs. Diffusion constants were calculated using the ZLC method, with micropore diffusion limitation found for ZIF-8 and MIL-100 monolithic frameworks at 2.7 x 10 -8 and 4.6 x 10 -8 cm 2 s -1 respectively. This diffusion limitation was also confirmed by IR spectroscopy with an increasing concentration of C—H bands on the MOF substrate post-reaction. Mass transfer coefficients, also calculated by ZLC, revealed increased mass transport for the hydrophobic ZIF-8 framework, which perhaps aids in the γ-valerolactone selectivity over side products that are produced in the absence of catalytic material, as seen for MIL-88B and MIL-100 after multiple uses.

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

Disordered Rocksalts as High‐Energy and Earth‐Abundant Li‐Ion Cathodes

To address the growing demand for energy and support the shift toward transportation electrification and intermittent renewable energy, there is an urgent need for low‐cost, energy‐dense electrical storage. Research on Li‐ion electrode materials has predominantly focused on ordered materials with well‐defined lithium diffusion channels, limiting cathode design to resource‐constrained Ni‐ and Co‐based oxides and lower‐energy polyanion compounds. Recently, disordered rocksalts with lithium excess (DRX) have demonstrated high capacity and energy density when lithium excess and/or local ordering allow statistical percolation of lithium sites through the structure. This cation disorder can be induced by high temperature synthesis or mechanochemical synthesis methods for a broad range of compositions. DRX oxides and oxyfluorides containing Earth‐abundant transition metals have been prepared using various synthesis routes, including solid‐state, molten‐salt, and sol‐gel reactions. This review outlines DRX design principles and explains the effect of synthesis conditions on cation disorder and short‐range cation ordering (SRO), which determines the cycling stability and rate capability. In addition, strategies to enhance Li transport and capacity retention with Mn‐rich DRX possessing partial spinel‐like ordering are discussed. Finally, the review considers the optimization of carbon and electrolyte in DRX materials and addresses key challenges and opportunities for commercializing DRX cathodes.

Li-ion batteries↗
DOE OSTI2025

Thermodynamic Stability and Site‐Specific Distribution of Graphitic and Pyridinic Nitrogen in Graphene Moiré on Ru(0001)

Abstract Graphene‐like materials are of interest for large‐scale hydrogen storage applications due to their lightweight, durable, and scalable properties. Nitrogen‐doping minimizes kinetic limitations in diffusion and recombination on surfaces, however, the role of graphitic nitrogen (GN) and pyridinic nitrogen (PN) is not well understood. Nitrogen‐doped graphene is synthesized on Ru(0001) using chemical vapor deposition (CVD) of pyridine and ion irradiation. Scanning tunneling microscopy (STM), x‐ray photoelectron spectroscopy (XPS), and density functional theory (DFT) are used to identify the structure, location, and thermodynamic stability of nitrogen species within the graphene moiré. CVD of pyridine results in a low nitrogen concentration (<0.1at%), while the post‐growth nitrogen ion irradiation allows us to increase the concentration further. The concentration of GN and PN is controlled by varying the ion dose and annealing temperature. Comparison of measured and simulated STM images of GN and PN yield an excellent agreement, allowing us to confidently establish that GN is preferentially located near the center of the Atop region, while PN is located in the valley region of the graphene moiré. This report explicitly confirms the site assignments and provides a foundation for the site synthesis and analysis of structural and electronic properties that drive the reactivity of N‐doped graphene.

Gedara, Buddhika S. A. [Physical and Computational↗
DOE OSTI2026

Mechanisms of Alkali Ionic Transport in Amorphous Oxyhalides Solid State Conductors

Amorphous oxyhalides have attracted significant attention due to their relatively high ionic conductivity (1 mS cm –1 ), excellent chemical stability, mechanical softness, and facile synthesis routes via standard solid‐state reactions. These materials exhibit an ionic conductivity that is almost independent of the underlying chemistry, in stark contrast to what occurs in crystalline conductors. In this work, we employ machine learning interatomic potentials to construct large‐scale molecular dynamics trajectories encompassing hundreds of nanoseconds to obtain statistically converged transport properties. We find that the amorphous state consists of chain fragments of metal‐anion tetrahedra of various lengths. By analyzing the residence time of alkali cations migrating around tetrahedrally‐coordinated metals, we find that oxygen anions limit alkali diffusion. By computing the full Einstein expression of the ionic conductivity, we demonstrate that the alkali transference number of these materials is strongly influenced by distinct‐particles correlations, while alkali transport is dictated by uncorrelated self‐diffusion. By extending this analysis to chemical compositions AMX 2.5 O 0.75 , spanning different alkaline (A = Li, Na, K), metallic (M = Al, Ga, In), and halogen (X = Cl, Br, I) species, we clarify why the diffusion properties of these materials remain largely insensitive to variations in atomic isovalent chemistry.

amorphous materials↗
DOE OSTI2022

Fast redox switches lead to rapid transformation of goethite in humid tropical soils: A Mössbauer spectroscopy study

Abstract Humid tropical forest soils experience frequent rainfall, which limits oxygen diffusion and creates redox heterogeneity in upland soils. In this study we gauged the effect of short‐term anoxic conditions on Fe mineralogy of relatively Fe‐ and C‐rich surface soils (C/Fe mole ratio ∼5) from a humid tropical forest in Puerto Rico. Soils subjected to 4‐d oxic/anoxic oscillation were characterized by selective chemical extractions, Mössbauer spectroscopy (MBS), and X‐ray diffraction. Chemical extraction data suggested that rapidly switching redox conditions had subtle effects on bulk Fe mineralogy. Mössbauer, on the other hand, indicated that (a) the soil Fe is a mixture of goethites of varying characteristics with minor contributions from ferrihydrite (<5%) and Fe(III)‐organic matter (OM) phases (∼10%), and (b) anoxic conditions rapidly transformed all forms of goethite to relatively stable Fe oxides. Such fast changes in goethite features could be due to rapid depletion and sorption of bio‐released structural Al or sorbed OM onto residual soil components. The rapid temporal changes in MBS parameters and corresponding pore water nominal oxidation state of C values suggest that Fe‐C transformations in these upland tropical soils are complex and intricately coupled. A comprehensive understanding of the fate of Fe, Al, and OM (Fe‐organic moieties) during redox switches and concurrent changes in pore water chemistry is critical for development of robust transport models in humid tropical soils, which are subject to episodic low‐redox events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Synthesis of wafer-scale SWCNT forests with remarkably invariant structural properties in a bulk-diffusion-controlled kinetic regime

Robust synthesis of vertically aligned carbon nanotubes (VACNT) at large scale is required to accelerate deployment of numerous cutting-edge devices to emerging commercial applications. To address this need, this work demonstrates that the structural characteristics of single-walled CNTs (SWCNTs) produced in a growth regime dominated by bulk diffusion of the gaseous carbon precursor are remarkably invariant over a broad range of process conditions (precursor concentration increased up to 30-fold, catalyst substrate area from 1 cm 2 to 180 cm 2 , growth pressure from 20 to 790 mbar, and gas flowrates up to 8-fold). A simple growth kinetics model is derived that accounts for both reactant bulk diffusion and competing byproduct formation in the presence of excess hydrogen, which quantitatively reproduces the experimental data in the entire isothermal parameter space. The model enables critical predictions for process scale-up optimization including a potential 6-fold increase in CNT production rate, a ~90% carbon conversion efficiency in select conditions, and elimination of reaction rate decay observed at high pressures by using a hydrogen-free growth environment. This model-guided opportunity to substantially improve synthesis throughput and efficiency, in conjunction with the remarkably invariant CNT properties, is appealing for reliable VACNT device fabrication at both small and industrial scales.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamics of coupled thermoacoustic modes: Noise and frequency spacing effects

Practical combustion systems consist of multiple thermoacoustic modes that are linearly stable or unstable at a wide range of frequencies. When coupled, the modal interactions are a strong function of the frequency spacing between them – having very significant implications on interpretation of linear stability or instability of modes from experimental data. Frequency spacing between the modes particularly influences the stability and existence of potential limit cycle solutions. In this paper, we extend earlier deterministic studies to include the effects of noise, studying how deterministic dynamics change with the addition of noise and the simultaneous impact of frequency spacing (i.e., closely, or widely spaced). Specifically, this paper characterizes three key effects of noise, which are (1) distribution of limit cycle amplitudes (“diffusion”), (2) shift in the average limit cycle amplitudes (“noise induced drift”), and (3) introduction of new bifurcations that are not present in the deterministic system. Key conclusions from this work are that in highly noisy practical environments, a deterministic/ low noise analysis will not be sufficient to explain even critical qualitative features of the system – such as existence of stable limit cycles or their amplitudes. Furthermore, the nature of noise-induced effects on nonlinear modal interactions can be fundamentally different for degenerate, closely spaced, and widely spaced modes.

42 ENGINEERING↗
DOE OSTI2020

A fast implicit solver for semiconductor models in one space dimension

Several different approaches are proposed for solving fully implicit discretizations of a simplified Boltzmann-Poisson system with a linear relaxation-type collision kernel. This system models the evolution of free electrons in semiconductor devices under a low-density assumption. At each implicit time step, the discretized system is formulated as a fixed-point problem, which can then be solved with a variety of methods. A key algorithmic component in all the approaches considered here is a recently developed sweeping algorithm for Vlasov-Poisson systems. A synthetic acceleration scheme has been implemented to accelerate the convergence of iterative solvers by using the solution to a drift-diffusion equation as a preconditioner. The performance of four iterative solvers and their accelerated variants has been compared on problems modeling semiconductor devices with various electron mean-free-path.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effect of post-weld heat treatment on microstructure and mechanical properties of rotary inertia friction welded dissimilar 422-4140 martensitic steel piston joints

The study investigated how post‐weld heat treatment (PWHT) temperature affects the microstructure and localized deformation/fracture during bend testing of rotary inertia friction welds (RIFW) between AISI 422 stainless steel and AISI 4140 steel. RIFW produced a fully martensitic interface with approximately 550 HV hardness in both the thermo-mechanically affected (TMAZ) and heat‐affected zones (HAZ). Due to differences in temper resistance, the 4140 TMAZ/HAZ softened progressively under PWHT temperatures from 525 °C to 700 °C, while the 422 TMAZ unexpectedly maintained about 550 HV up to 600 °C before significantly softening at temperatures ≥625 °C. This asymmetric softening generated steep hardness gradients across the interface at temperatures ≤600 °C. Furthermore, carbon migration across the interface was minimal up to 600 °C, moderate at 625 °C, and by 700 °C produced a carbide‐rich eutectoid layer in the 422 TMAZ alongside a carbon‐depleted soft ferrite layer in the 4140 TMAZ. Strain during bending was PWHT‐dependent, concentrating on the 4140 side; in as-welded joints, the high hardness led to deformation and crack initiation in the base metals, whereas in PWHT samples, cracking initiated in the softened 4140 TMAZ near the interface. The intermediate PWHT temperature of 625 °C offered the best balance of limited carbon diffusion across the interface, relatively low peak weld hardness and minimized hardness gradients across the interface, more homogenous deformation, and good bend test performance.

4140 low-alloy steel↗
DOE OSTI2020

Oxide structure of air-passivated U-6Nb alloy thin films

U-6Nb is a uranium alloy containing 6 wt% niobium that possesses high corrosion resistance. The structure and composition of the passivating oxide layer formed on air-aged U-6Nb, which gives the material its corrosion resistant properties, was characterized using surface scattering techniques. Stable oxide layers formed on the surface of a set of U-6Nb alloy thin films under ambient conditions were investigated using neutron reflectometry (NR), x-ray reflectometry (XRR) and grazing incidence x-ray diffraction (GIXD). The passivating oxide was composed of approximately 27% U, 5% Nb, and 68% O, primarily consisting of a thin niobium oxide layer (5.5 ± 0.4 nm) separating a thicker UO2 layer (27.1 ± 2.3 nm) from the underlying U-6Nb alloy. A critical density of enriched niobium oxide at the metal-oxide interface is hypothesized to limit oxygen diffusion and confer high corrosion resistance to the alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Color as a tool for quantitative analysis of heterogeneous polymer degradation

Polymers are widely used industry, but degrade due to environmental conditions, such as elevated temperature or radiation. The degradation of polymers is dominated oxidation in oxygenated environments. Such oxidation is commonly referred to as thermal-oxidation or photo-oxidation and in many situations is heterogeneous through thickness of the polymer. While heterogeneous oxidation of polymers is a common phenomenon, the effect on the lifetime of polymers is not well understood as methods to evaluate such oxidation are time consuming and challenging to implement. Here we show a novel method for quantitative evaluation of heterogeneous oxidation in polymers using a fast and accurate color approach. The proposed color approach was observed to improve the efficiency of obtaining heterogenous oxidation profiles in half the time of common approaches, such as indention or Fourier transform infrared spectroscopy. Furthermore, the method was found to produce trends similar to that observed using common approaches and the applicability of the method was verified by mapping oxygen concentration through thickness using energy dispersive X-ray spectroscopy. Lastly, the utility of the proposed method was demonstrated through the generation of an oxidation design space, which is envisioned to assist with designing experiments towards estimation of polymer lifetime. Finally, we anticipate this work to be a starting point for further evaluation of heterogeneous oxidation in polymers using optical properties.

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

Nonequilibrium Behavior in Isoprene Secondary Organic Aerosol

Recent studies have shown that instantaneous gas-particle equilibrium partitioning assumptions fail to predict SOA formation, even at high relative humidity (~85%), and photochemical aging seems to be one driving factor. In this study, we probe the minimum aging time scale required to observe nonequilibrium partitioning of semivolatile organic compounds (SVOCs) between the gas and aerosol phase at ~50% RH. Seed isoprene SOA is generated by photo-oxidation in the presence of effloresced ammonium sulfate seeds at <1 ppbv NO x , aged photochemically or in the dark for 0.3–6 h, and subsequently exposed to fresh isoprene SVOCs. Our results show that the equilibrium partitioning assumption is accurate for fresh isoprene SOA but breaks down after isoprene SOA has been aged for as short as 20 min even in the dark. Modeling results show that a semisolid SOA phase state is necessary to reproduce the observed particle size distribution evolution. The observed nonequilibrium partitioning behavior and inferred semisolid phase state are corroborated by offline mass spectrometric analysis on the bulk aerosol particles showing the formation of organosulfates and oligomers. Here, the unexpected short time scale for the phase transition within isoprene SOA has important implications for the growth of atmospheric ultrafine particles to climate-relevant sizes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Unveiling the Structure and Dynamics of Water Confined in Colloidal Boehmite Suspensions

Thin fluid layers confined between nanoparticles play an important role in several natural and industrial systems, including radioactive wastes stored in tanks at the U.S. Department of Energy’s Hanford site. Multimodal neutron and computational approaches have been integrated to examine the properties of one, two, and four layers of water (H 2 O or D 2 O) on nanoparticulate, hydrous, or deuterated boehmite (γ-AlOOH or γ-AlOOD). Exposure of deuterated boehmite to H 2 O at 90 °C yielded rapid H/D isotopic exchange likely driven by a Grotthus-like proton-hopping mechanism. The in-plane and out-of-plane vibrations of the structural hydroxyls involved in this exchange were observed for both the hydrated and deuterated boehmite. These observations were confirmed by molecular dynamics simulations that also showed that single water molecules on the (010) surface bond to the structure via four bonds: two from hydrogens (deuteriums) on the water to surrounding oxygens and two from the water’s oxygen to surrounding OD/OH. Bulk water-like properties began to appear once four monolayers of water had been added, but steric crowding limited water diffusion rates once two monolayers had been added. The super-Arrhenius temperature dependence observed at four monolayers indicated glass-like behavior in a well-formed hydrogen-bonding network. Such a network is not sufficiently developed, however, when the surface water coverage is less than four layers. In conclusion, the unique nature of these layers can provide critical information for understanding forces between particles in proximity, and resultant effects on suspension rheology.

Anovitz, Lawrence M. [Oak Ridge National Laborator↗
DOE OSTI2026

A Prodrug Strategy to Conditionally Trap Therapeutic Payloads for Improved Tumor Retention

Altered extracellular proteolysis has been exploited to selectively activate therapeutics in diseases such as cancer; however, once activated, extracellular drugs can diffuse away, limiting efficacy. We address this challenge by coupling proteolytic activation with membrane tethering to retain drugs within diseased tissue. To accomplish this, we developed “restricted interaction peptides” (RIPs), a delivery platform that leverages elevated proteolytic activity to activate membrane-interacting peptides, localizing cargos near the site of proteolysis. We demonstrate that RIPs can deliver diverse therapeutic cargos, including cytotoxins and radioisotopes. As proof of concept, we engineered “FRIP,” a RIP designed for cleavage by fibroblast activation protein (FAP), an endoprotease upregulated in solid tumors and fibrosis. Efficient P4–P4’ substrate sequences were identified and incorporated into FRIPs. Cell-based studies showed that, upon activation, the peptide adhered to membranes rapidly internalized and successfully delivered therapeutic cargos. Consistent with this, FRIPs delivering MMAE inhibited proliferation in an FAP-dependent manner. Imaging studies confirmed tumor targeting with minimal uptake in normal tissues. Finally, FRIPs delivering MMAE or Cu-67 exhibited potent antitumor effects. These findings establish membrane tethering as a strategy to enhance drug retention.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Colloidal Clusters and Networks Formed by Oppositely Charged Nanoparticles with Varying Stiffnesses

Colloidal clusters and gels are ubiquitous in science and technology. Particle softness has a strong effect on interparticle interactions; however, our understanding of the role of this factor in the formation of colloidal clusters and gels is only beginning to evolve. Here, we report the results of experimental and simulation studies of the impact of particle softness on the assembly of clusters and networks from mixtures of oppositely charged polymer nanoparticles (NPs). Experiments were performed below or above the polymer glass transition temperature, at which the interaction potential and adhesive forces between the NPs were significantly varied. Hard NPs assembled in fractal clusters that subsequently organized in a kinetically arrested colloidal gel, while soft NPs formed dense precipitating aggregates, due to the NP deformation and the decreased interparticle distance. Importantly, interactions of hard and soft NPs led to the formation of discrete precipitating NP aggregates at a relatively low volume fraction of soft NPs. A phenomenological model was developed for interactions of oppositely charged NPs with varying softnesses. Furthermore, the experimental results were in agreement with molecular dynamics simulations based on the model. This work provides insight on interparticle interactions before, during, and after the formation of hard–hard, hard–soft, and soft–soft contacts and has impact for numerous applications of reversible colloidal gels, including their use as inks for additive manufacturing.

77 NANOSCIENCE AND NANOTECHNOLOGY↗