https://www.nature.com/articles/s41467-022-29345-x Skip to main content Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript. Advertisement Advertisement Nature Communications * View all journals * Search * My Account Login * Explore content * About the journal * Publish with us * Sign up for alerts * RSS feed 1. nature 2. nature communications 3. articles 4. article Stress distribution and surface shock wave of drop impact Download PDF Download PDF * Article * Open Access * Published: 31 March 2022 Stress distribution and surface shock wave of drop impact * Ting-Pi Sun^1, * Franco Alvarez-Novoa^2, * Klebbert Andrade^2, * Pablo Gutierrez ORCID: orcid.org/0000-0002-4673-1656^3, * Leonardo Gordillo ORCID: orcid.org/0000-0002-9516-1346^2 & * Xiang Cheng ORCID: orcid.org/0000-0002-2759-764X^1 Nature Communications volume 13, Article number: 1703 (2022) Cite this article * 8161 Accesses * 326 Altmetric * Metrics details Subjects * Fluid dynamics * Fluids Abstract Drop impact causes severe surface erosion, dictating many important natural, environmental and engineering processes and calling for substantial prevention and preservation efforts. Nevertheless, despite extensive studies on the kinematic features of impacting drops over the last two decades, the dynamic process that leads to the drop-impact erosion is still far from clear. Here, we develop a method of high-speed stress microscopy, which measures the key dynamic properties of drop impact responsible for erosion, i.e., the shear stress and pressure distributions of impacting drops, with unprecedented spatiotemporal resolutions. Our experiments reveal the fast propagation of self-similar noncentral stress maxima underneath impacting drops and quantify the shear force on impacted substrates. Moreover, we examine the deformation of elastic substrates under impact and uncover impact-induced surface shock waves. Our study opens the door for quantitative measurements of the impact stress of liquid drops and sheds light on the origin of low-speed drop-impact erosion. Introduction Laozi, the ancient Chinese philosopher in the fifth century BCE, has long noticed that water, although the softest and weakest material known in his time, is effectual in eroding hard substances^1. Although Laozi used this attribute of water only as a metaphor to extol the virtues of humbleness, flexibility, and persistence, it was a physically nontrivial observation that dripping water drops, with zero shear modulus and readily deformable, can erode hard solid substrates with finite yield stresses (Fig. 1a, b)^2. Beyond its philosophical meaning, erosion by drop impact is relevant for a wide range of natural, environmental, and engineering processes including soil erosion^3,4,5,6,7, preservation of heritage sites^8, wear of wind and steam turbine blades^9, and cleaning and peening of solid materials (e.g. silicon wafers)^10,11,12. While the impact damage on solid substrates caused by high-speed compressible liquid drops with the impact velocity on the order of a few hundred meters per second has been well explored^13,14, our understanding of the impact erosion of low-speed incompressible drops relevant to most natural and engineering processes is still rudimentary. Why is low-speed drop impact erosive? What are the key dynamic features of drop impact that lead to its unexpected ability in erosion? The answer to these questions would provide not only the testimony to the wisdom of the ancient philosopher but also fundamental insights into the early-time dynamics of drop impact. Fig. 1: Erosion by dripping drops. figure 1 a Erosion crater by drop impact on a granular medium made of 90 mm glass beads. The crater is created by the impact of a single water drop of diameter D = 3 mm at impact velocity U = 2.97 m/s. b Erosion crater by drop impact on a plaster slab. The crater is created by 2500 repeated impacts of water drops of D = 3 mm at U = 2.6 m/s. Scale bars: 5 mm. Full size image Due to its ubiquity in nature and industry, drop impact has evoked long-lasting research interests since the early study of Worthington almost one and a half centuries ago^15. Particularly, significant progress has been made over the last two decades in understanding drop impact thanks to the fast advance of high-speed photography techniques. Nevertheless, limited by direct imaging, most current studies focused on the kinematics of impacting drops such as the splashing threshold, the maximum spreading diameter and contact time, and the formation of cushioning air layers underneath drops (see refs. ^16,17 and references therein). Very few experiments have been performed probing the dynamic properties of drop impact that are directly responsible for drop-impact erosion^18. As an important dynamic factor, the impact force of liquid drops has begun to attract attention in recent years^18,19,20,21,22,23,24,25,26. However, it becomes clear from these recent studies that the unusual ability of an impacting drop in erosion cannot stem from its impact force alone, as the maximum impact force induced by a millimetric water drop falling near its terminal velocity is very weak, more than an order of magnitude smaller than that generated by the impact of a solid sphere of similar size, density, and velocity^23. The average impact pressure of the drop is even smaller due to the large contact area formed by the spreading drop^2. Thus, instead of impact force or average impact stresses, the erosion ability of drop impact must originate from the unique spatiotemporal structure of its impact stresses as well as the dynamic response of impacted substrates under such stresses. Nevertheless, quantitative measurements on the stress distributions of drop impact have not been achieved heretofore due to the limitation of existing experimental tools. Recent attempts using an array of miniature force sensors based on microelectromechanical systems (MEMS) provided only the coarse-grained dynamics of impact pressure without a sufficient spatial resolution and failed to measure the shear stress distribution of impacting drops^27,28. Here, we develop an experimental method--high-speed stress microscopy--to measure the impact stress of drop impact on solid elastic substrates. The method integrates the imaging techniques of traction force microscopy^29, laser-sheet microscopy, and high-speed photography, which allows us to map the temporal evolution of the pressure and shear stress distributions underneath millimeter-sized drops in fast impact events with unprecedented spatiotemporal resolutions. Results High-speed stress microscopy We embed low-concentration (0.23% v/v) fluorescent polystyrene particles of diameter 30 mm in a cross-linked polydimethylsiloxane (PDMS) gel as tracers to track the deformation of the gel under impact. The PDMS gel surface is hydrophobic with the water contact angle ~90^[?]. Young's modulus of the gel is fixed at E = 100 kPa in our experiments, although gels with E up to 420 kPa and with hydrophilic surfaces have also been tested (Methods). A thin laser sheet of 30-mm thickness illuminates the gel from a side and excites the fluorescent tracers within the sheet (Fig. 2a). The sheet is finely adjusted to be normal to the impacted surface and to pass through the center of impacting drops. A high-speed camera focusing on the sheet images the motion of tracers at 40,000 frames per second. Fig. 2: High-speed stress microscopy. figure 2 a A schematic showing the principle of high-speed stress microscopy. A drop falls from a height h and impacts onto the surface of a particle-embedded PDMS gel at an impact velocity U. Lower left inset: An image of the PDMS gel embedded with fluorescent particles under the illumination of the laser sheet. b Comparison of the shear stress (top) and the pressure (bottom) induced by the impact of a steel sphere obtained from high-speed stress microscopy (left) and from finite element simulations (right). The diameter and the impact velocity of the steel sphere are D = 3.16 mm and U = 0.49 m/s, respectively. The stresses are measured at time t = 0.25 ms after the instant when the sphere first touches the solid surface. Full size image We track the displacements of the tracers from the high-speed video using digital image correlation (DIC). An interrogation window of 384 mm by 384 mm with 70% overlap is adopted in DIC, which gives a spatial resolution of 115 mm. The temporal resolution is 0.025 ms, set by the frame rate of high-speed photography. To reduce measurement errors, we first average the cross-correlation fields of DIC from five repeated impacts on the same gel at the same impact location. The surface is fully dried between two consecutive experimental runs. The stress measurements are the outcome of a further average over three different averaged displacement fields for impacts on different gels or impacts on the same gels at different impact locations. Thus, one data point represents the average result of total 15 different experimental runs. Stress fields depend on strain fields, which are the derivative of the displacement fields obtained from DIC. A smoothing procedure is necessary in order to reduce the noise of differentiation. We implement the moving least squares (MLS) interpolation method to obtain a continuously differentiable displacement field u(r, z) = [u[ r](r, z), u[z](r, z)] from the discrete displacement field of DIC^30. A third-order polynomial basis is adopted in the interpolation (Supplementary Information (SI) Section 1). When the deformation is small, the strain components in cylindrical coordinates follow: $${\varepsilon }_{rr}=\frac{\partial {u}_{r}}{\partial r},{\ varepsilon }_{zz}=\frac{\partial {u}_{z}}{\partial z},{\varepsilon }_ {\theta \theta }=\frac{{u}_{r}}{r},{\varepsilon }_{rz}=\frac{1}{2}\ left[\frac{\partial {u}_{r}}{\partial z}+\frac{\partial {u}_{z}}{\ partial r}\right],$$ (1) where we take the advantage of the cylindrical symmetry of the drop-impact geometry. By assuming the PDMS gels are isotropic and linear following the generalized Hooke's law at small strains, we calculate the stress fields using the linear stress-strain relation: $${\sigma }_{ij}=\lambda {\varepsilon }_{b}{\delta }_{ij}+2G{\ varepsilon }_{ij},$$ (2) where l = En/[(1 + n)(1 - 2n)] is the Lame coefficient, G = E/[2(1 + n)] is the shear modulus, d[ij] is the Kronecker delta and e[b] [?] e[z z] + e[rr] + e[thth] is the bulk strain. s[rz] gives the shear stress t , whereas s[zz] gives the pressure p. PDMS gels are nearly incompressible with Poisson's ratio n close to 0.5, which result in a large l. But the bulk strain e[b] is close to 0 in this limit. Therefore, the impact pressure cannot be accurately determined from the product of le[b] in Eq. (2). Instead, we adopt a quasi-steady state assumption to calculate the pressure^30, a procedure detailed further in SI Section 2. We have verified the assumption by comparing the inertial force and the elastic force in the impact process (SI Section 2) and by comparing experimental and numerical results on the impact pressure of solid-sphere impact (see below). The shear stress, on the other hand, is not affected by the nearly incompressible condition. The surface stresses and displacements are finally obtained at a location slightly below the original impacted surface (Methods). As a calibration and the basis of comparison, we first measure the pressure and shear stress induced by the impact of a solid steel sphere of diameter D = 3.16 mm at impact velocity U = 0.49 m/s and compare the results with those from finite element simulations (Methods). Experimental measurements agree well with the numerical results, validating the accuracy of high-speed stress microscopy (Fig. 2b). For drop impact, our drops are made of an aqueous solution of sodium iodide (60% w/w), which has a density r = 2.2 g/ml and a viscosity e = 1.12 mPa s^31. The surface tension of the solution s [?] 81.3 mN/m from pendant-drop tensiometry^32, which is slightly larger than that of water. We fix the diameter D and impact velocity U of drops in our experiments. Drops of D = 3.49 mm impact normally on the surface of PDMS gels at U = 2.97 m/s, yielding a Reynolds number Re = rUD/e = 20,360 and a Weber number We = rDU^2/s = 833. Thus, the drop impact is dominated by fluid inertia at early times. We focus on drop impact at early times below, when the shear stress and pressure of impacting drops are high for strong erosion. Positions and times are reported in dimensionless forms using D and D/U as the corresponding length and time scale, respectively. Impact shear stress Surface erosion is the direct consequence of the shear stress of drop impact. Figure 3a, b compares the temporal evolution of the shear stress of solid-sphere impact and liquid-drop impact. Upon the impact, spatially non-uniform shear stresses quickly develop in both cases. However, while the position of the maximum shear stress of solid-sphere impact is stationary near the impact axis at r = 0.095, the maximum shear stress of drop impact propagates radially with the spreading drop. The kymographs of the surface shear stress, t(r, z = 0, t), of the two impact processes are shown in Fig. 3c, d, highlighting further the fast propagation of the maximum shear stress of drop impact. Fig. 3: Shear stress distribution of solid-sphere impact and drop impact. figure 3 a, c The temporal evolution of the shear stress t(r, z, t) and the kymograph of the surface shear stress t(r, z = 0, t) of solid-sphere impact. The diameter and the impact velocity of the steel sphere are D = 3.16 mm and U = 0.49 m/s, respectively. b, d The temporal evolution of t(r, z, t) and the kymograph of t(r, z = 0, t) of drop impact. D = 3.49 mm and U = 2.97 m/s for the liquid drop. The red line in (d) indicates the position of the turning point. t = 0 corresponds to the instant when the impactors first touch the surface of the PDMS gels. Scale bar in (a) is 1 mm. Full size image To understand the origin of the maximum shear stress of drop impact, we correlate the position of the maximum shear stress r[s] with the shape of impacting drops (Fig. 4a). Two kinematic features are analyzed: the tip of the expanding lamella r[lm] and the turning point r[t], where the drop body connects to the root of the lamella (Fig. 4a inset). It should be emphasized that r[t] is not the contact line of the drop. The ejection of the lamella occurs around t [?] We^-2 /3 = 0.0104^33, which is shorter than the temporal resolution of our experiments. While r[lm] moves fastest, r[s] follows closely behind r [t]. Thus, the maximum shear stress arises from the strong velocity gradient near the turning point, where the flow changes rapidly from the downward vertical direction (the -z direction) within the drop body to the horizontal radial direction (the r direction) inside the narrow lamella^34. Quantitatively, r[t](t) follows the well-known square-root scaling \({r}_{t}(t)=\sqrt{6t}/2\approx 1.22\sqrt{t}\) established by many previous experiments^23,33,35,36,37,38,39. In comparison, r[s](t) also shows a square-root scaling with a slightly smaller prefactor \({r}_{s}(t)\approx \sqrt{t}\). The positions of the maximum shear stress and the turning point are independent of the wettability or Young's modulus of PDMS gels (Fig. 4b). Fig. 4: Surface shear stress of drop impact. figure 4 a The position of the lamellar tip r[lm], the turning point r[t], the maximum shear stress r[s] and the maximum pressure r[p] as a function of time t. The dashed line indicates \(r=\sqrt{6t}/2\). Upper inset: Definition of the kinematic features of an impacting drop. b The position of the maximum shear stress of drop impact r[s](t) on substrates of different Young's moduli and wettability. The error bars are the standard deviation of 15 experimental runs. c The rescaled surface shear stress \(\tau /\sqrt{t}\) as a function of the rescaled radial position \(r/\sqrt{t}\). The dashed line is the prediction of Eq. (3) with the modified scaling function f(x). d, e Two components of the surface shear stress: t[1] = G([?]u[z]/[?]r)|[z = 0] and t[2] = G([?]u[r]/[?]z)|[z = 0]. f The total surface shear stress t = t[1] + t[2]. The shear stresses in (c-f) are nondimensionalized by rU^2 and share the same color code for time, as indicated in (c). The diameter and the impact velocity of the liquid drop are D = 3.49 mm and U = 2.97 m/s. Full size image Philippi et al.^40 proposed that the shear stress of incompressible drops on infinitely rigid substrates possesses a self-similar dynamic structure when t - 0^+, $$\tau (r,z=0,t)=2\sqrt{\frac{6}{{\pi }^{3}Re}}\frac{1}{\sqrt{t}}f\ left(\frac{r}{\sqrt{t}}\right)\quad {{{{{{{\rm{for}}}}}}}}\quad r\le {r}_{t}(t),$$ (3) where the scaling function f(x) = x/(3 - 2x^2) dictates a finite-time singularity at the turning point r[t](t). Here, t(r, t) is non-dimensionalized by the inertial pressure rU^2. Note that the much stronger water-hammer pressure rUc associated with the compression wave occurs on the time scale of a few nanoseconds, which is too short to be relevant in our current experiments^7,20,23. Here, c is the speed of sound in liquid. Inspired by the self-similar hypothesis, we plot \(\tau \sqrt{t}\) versus \(r/\sqrt{t}\) of our experimental results (Fig. 4c), which shows a good collapse at a small r away from the singular region. With a modified scaling function f(x) = x/(1 - x^2) to count the different temporal scalings of r[t] and r[s], the collapsed data quantitatively agrees with Eq. ( 3) (the dashed line in Fig. 4c). Thus, our study provides experimental evidence on the propagation of shear stress of drop impact and demonstrates the self-similar structure of shear stress at early times. Despite the general agreement with Eq. (3) at r < r[t], our experiments also reveal the unique features of drop impact on elastic deformable substrates, absent in the theoretical consideration of drop impact on infinitely rigid substrates. The shear stress on the surface of an elastic substrate is given by t = G([?]u[r]/[?]z + [?]u[z]/[?]r ) (Eqs. (1) and (2)), where G is the shear modulus of the substrate and u[r] and u[z] are the radial and vertical displacement of the substrate surface. We find |[?]u[z]/[?]r| > |[?]u[r]/[?]z| (Fig. 4d-f), suggesting the dominant role of the vertical velocity of the impacting drop at the contact surface v[z](r, z = 0) on the shear stress. Note that \({u}_{z}(r,z=0,t)=\int\nolimits_{0}^{t}{v}_{z}(r,z =0,t)dt\). For drop impact on infinitely rigid substrates, v[z](r, z = 0, t) = 0 because of the no-penetration boundary condition, which inevitably gives [?]u[z]/[?]r = 0. Instead, the shear stress of drop impact on infinitely rigid substrates arises from the gradient of the radial velocity, [?]v[r]/[?]z, within the boundary layer near the contact surface. As u[z](r, z = 0) is mainly determined by the pressure distribution on the contact surface at high Re, this finding illustrates the intrinsic coupling between the impact pressure and shear stress of drop impact on elastic substrates. The effect of the finite stiffness of the impacted substrate also manifests in the shear force of impacting drops. By integrating the shear stress over the contact area, we obtain the shear force, \({F}_ {d}(t)=2\pi \int\nolimits_{0}^{{r}_{lm}}\tau (r,z=0,t)rdr\), which quantifies the total erosion strength of drop impact. Although F[d](t ) is independent of the wettability of the impacted surface, it increases with Young's modulus following a scaling \({F}_{d} \sim \ sqrt{E}\) within the range of our experiments (Fig. 5a, b). Thanks to the spreading of the maximum shear stress, drop impact and solid-sphere impact show comparable peak shear forces under similar impact conditions (Fig. 5a). Fig. 5: Shear force of drop impact. figure 5 a The shear force F[d](t) of drop impact and solid-sphere impact on the impacted surface of different Young's moduli and wettability. The leftmost curve is for solid-sphere impact and the rest are for drop impact as indicated in the legend. b The dimensionless shear force F[ d] of drop impact scaled by \({E}^{1/2}{(\rho {U}^{2})}^{1/2}{D}^{2} \), as suggested by the theoretical analysis of drop impact on elastic substrates (SI Section 3). The diameter and the impact velocity of the liquid drop are D = 3.49 mm and U = 2.97 m/s. D = 3.16 mm and U = 0.49 m/s for the steel sphere. Full size image Impact pressure and surface shock wave Although subject to larger experimental errors due to the nearly incompressibility of PDMS, the pressure (i.e. normal stress) distribution underneath impacting drops p(r) can be also measured by high-speed stress microscopy (SI Section 2). Similar to the shear stress, we observe a non-central pressure maximum propagating radially with the spreading drop (Fig. 6b, d). The dynamics are again in sharp contrast to the pressure of solid-sphere impact, where the maximum impact pressure is fixed at the impact axis r = 0 (Fig. 6a, c). The existence of the propagating non-central pressure maximum has been predicted by several theories and simulations of drop impact^38, 40,41,42,43. Nevertheless, to the best of our knowledge, such a counter-intuitive prediction has not been directly verified in experiments heretofore. While our measurements qualitatively confirm the prediction, we find that the maximum pressure falls behind the maximum shear stress (Fig. 4a), a feature unexpected from drop impact on infinitely rigid substrates^40. Fig. 6: Pressure distribution of solid-sphere impact and drop impact. figure 6 a, c The temporal evolution of the pressure p(r, z, t) and the kymograph of the surface pressure p(r, z = 0, t) of solid-sphere impact. The diameter and the impact velocity of the steel sphere are D = 3.16 mm and U = 0.49 m/s. b, d The temporal evolution of p(r, z, t) and the kymograph of p(r, z = 0, t) of drop impact. D = 3.49 mm and U = 2.97 m/s. The red line in d indicates the position of the turning point. The time when the negative pressure emerges t[c] = 0.106 is also indicated. Scale bar in (a) is 1 mm. Note that the pressure scale is kPa for the solid-sphere impact and Pa for the drop impact. Full size image More interestingly, a negative pressure emerges in front of the turning point r[t] at t[c] [?] 0.106 (Fig. 6b, d). In the meantime, we also observe the propagation of surface disturbance on the gel surface away from the stress maxima above t[c] (Fig. 7a). Both suggest the formation of a surface acoustic wave--the classic Rayleigh wave--in the gel. Since the speed of the turning point \({V}_{t}(t)=d {r}_{t}/dt=\sqrt{6}/(4\sqrt{t})\) increases with decreasing t, the stress maxima associated with the turning point spread supersonically at early times (Fig. 7b). Thus, a shock front forms near r[t] on the impacted surface when t < t[c]. The Rayleigh wave finally overtakes the turning point and is released in front of the spreading drop in an explosion-like process above t[c], giving rise to the negative pressure and the propagation of surface disturbance. Based on the above picture, the speed of the surface wave can be estimated as \ ({V}_{t}({t}_{c})=\sqrt{6}/(4\sqrt{{t}_{c}})=1.88\), which quantitatively matches the speed of the Rayleigh wave^44 $${V}_{R}=\frac{1}{M}\left[\sqrt{\frac{1}{2(1+\nu )}}\frac{0.862+1.14 \nu }{1+\nu }\right]=1.89.$$ (4) Here, the Mach number \(M\equiv U\sqrt{{\rho }_{s}/E}=0.292\) with r[ s] = 0.965 g/cm^3 and n = 0.49 is the density and Poisson's ratio of PDMS. Fig. 7: Surface shock wave of drop impact. figure 7 a The kymograph of the radial displacement of gel surface u[r](r, t) (top) and the vertical displacement of gel surface u[z](r, t) (bottom) induced by drop impact. The diameter and the impact velocity of the drop are D = 3.49 mm and U = 2.97 m/s. The displacements are nondimensionalized by D. b Competition between the speed of the turning point V[t] and the speed of the Rayleigh wave V[R]. The supersonic regime before t[c], where V[t] > V[R], is indicated. Symbols are from experiments, where the error bars are the standard deviation of 15 experimental runs. The red dashed line shows \({V}_ {t}=\sqrt{6}/(4\sqrt{t})\). c The numerical solution of u[r](r, t) (top) and u[z](r, t) (bottom) induced by drop impact. The red dashed lines in (a) and (c) indicate the position of the maximum pressure that drives the Rayleigh wave. Note that the maximum pressure propagates with the turning point in the simulation, which is faster than that observed in experiments (Fig. 6d). See SI Section 3 for more details about the numerical solution. Full size image Encouraged by the quantitative agreement between Eq. (3) and experiments, we couple the theoretical impact pressure and shear stress of incompressible drops on infinitely rigid surfaces with the Navier-Lame equation of semi-infinite elastic media (SI Section 3). The dimensional analysis of the governing equation and the boundary conditions suggests that the shear force should scale as \({F}_{d} \ sim {E}^{1/2}{(\rho {U}^{2})}^{1/2}{D}^{2}\), agreeing with our measurements at different E (Fig. 5b). Moreover, the numerical solution of the coupled equations qualitatively reproduces the formation of the shock-induced Rayleigh wave of drop impact, where a sharp surface wave with a well-defined peak emerges at t[c] [?] 0.1 and propagates with V[R] (Fig. 7c). The strong and sharp surface wave is produced by the mechanical resonance occurring when the speed of the stress maxima approaches the speed of the Rayleigh wave near t[c]. Such a resonant phenomenon does not exist for solid-sphere impact with stationary stress maxima. As a result, the surface Rayleigh wave of solid-sphere impact is more diffusive (Supplementary Fig. 2). Discussion Taken together, our high-speed stress microscopy reveals three unique dynamic features of drop impact. These features contribute to the ability of drop impact to erode solid surfaces and result in the distinct nature of drop-impact erosion, which is qualitatively different from that of solid-sphere-impact. (i) The spatiotemporal stress distributions of impacting drops are highly non-uniform. The radially propagating stress maxima simultaneously press and scrub impacted substrates, leading to a large erosion area and a high shear force. (ii) Because of the fast speed of the turning point at short times, a shock wave forms on impacted substrates, which substantially increases the erosion strength. Each impacting drop behaves like a tiny bomb, releasing its kinetic impacting energy explosively. (iii) A sharp shock-induced surface wave finally emerges from the explosion process. The resulting decompression wave weakens the cohesion of surface materials before the arrival of the shear stress maximum. Can drop impact induce surface shock waves on substrates stiffer than the PDMS gels used in our experiments? Although the turning point r[t ](t) exhibits a square-root temporal scaling \(\sqrt{6t}/2\), which suggests a divergent speed as t - 0^+ and therefore supports surface shock waves on infinitely rigid substrates, such a singular behavior is regularized at small times due to the development of compression waves and/or the effect of air cushioning. Nevertheless, even at the shortest time of our experiments of t = 0.021, the square-root scaling still holds well (Fig. 4a). Previous experiments on the kinematics of drop impact have shown the scaling at even smaller times down to t [?] 5 x 10^-4 before the ejection of lamella (Fig. 7 in ref. ^39), where the turning point is coincident with the apparent contact line of impacting drops. The speed of the turning point at this time is \({V}_{t}=\sqrt{6}/(4\sqrt{t})=27.4\). If we set the V[t ] to be the speed of the Rayleigh wave (Eq. (4)), the upper limit of Young's modulus of elastic substrates where we still expect to observe surface shock waves is E [?] 20 MPa. Here, n = 0.49 and r[s] = 1000 kg/m^3 are taken in the estimate. Compression waves develop within impacting drops at very early times t [?] U^2/c^2 [?] 4 x 10^-6 ^ 20. The speed of the turning point at such a small time scale is V[t] = 300, which gives a theoretical upper limit of Young's modulus E [?] 3 GPa, approaching the modulus of sedimentary rocks. The effect of air cushioning is more complicated. By preventing the on-axis contact, air cushioning eliminates both the compression waves and the divergent speed of the turning point as t - 0^+^42,43. Nevertheless, as discussed above, numerous experiments have repeatedly confirmed the square-root temporal scaling of the turning point in the ambient air at small time scales^23,33,35,36,37,38,39, supporting the fast propagation of the turning point in air at short times. More importantly, theories and simulations have both shown that air cushioning does not annihilate the fast propagation of the non-central pressure peak^40,42,43. Particularly, the underlying air layer induces a micron or sub-micron dimple-like deformation bounded by a kink structure at the bottom of an impacting drop^42,45. The pressure reaches the maximum underneath the kink^42,43, which propagates radially outwards at a speed 50 times higher than the impact velocity at early times of t ~ 10^-4.^42,45 Hence, surface shock waves induced by the fast propagation of the pressure maximum should persist in the presence of air cushioning on substrates of Young's modulus at least up to ~70 MPa. It is an open question if surface shock waves can sustain on even stiffer substrates. The radius and the thickness of the air dimple decrease at reduced ambient pressure^46, which eventually leads to the vanishing of air cushioning and the recovery of the singular dynamics of the contact line at short times. Lastly, it is worth noting that shock propagation within impacting drops and along impacted surfaces has been investigated theoretically for compressible drops with the liquid Mach number U/c ~ O(1)^13,14. Nevertheless, the shock process of incompressible drops with U/c [?] 1 that are relevant to most natural and industrial processes has not been discussed heretofore. The application of high-speed stress microscopy in our study demonstrates its great potential to measure the impact stress of low-speed incompressible drops in more diverse situations such as drop impact on patterned substrates, at reduced ambient pressure and with non-Newtonian drops^16,17. Methods Preparation and characterization of PDMS substrates PDMS elastomers were prepared from two-part Sylgard 184 silicone elastomer kits (Dow Corning). To adjust Young's modulus, E, of the cured PDMS gels, we mixed siloxane monomers with crosslinkers at a controlled mass ratio. In most experiments reported in this paper, we used a monomer-to-crosslinker mass ratio of 30:1, which yielded PDMS gels of E = 100 kPa. PDMS gels of E = 250 kPa (mass ratio 20:1) and 420 kPa (mass ratio 18:1) have also been tested to assess the effect of gel stiffness on the impact pressure and shear stress (Figs. 4b and 5). Fluorescent polystyrene (PS) particles of diameter 30 mm (ThermoFisher) were mixed into the two-part PDMS mixtures at a volume fraction of 0.23% before curing. The PDMS-particle mixture was finally vacuumed to remove air bubbles and placed in an oven at 90 ^[?] C overnight for curing. The fully cured gels have a fixed thickness of 6.5 mm and an area of 24 mm by 24 mm, which is much larger than the maximum spreading area of impacting drops. We modified the wettability of the surface of the PDMS gels via plasma modification. Without the modification, the untreated PDMS surfaces are hydrophobic with a contact angle of around 90^[?]. After the modification, the surfaces become hydrophilic with a contact angle of less than 10^[?]. Untreated hydrophobic gels were used in all the experiments reported in this study unless stated otherwise. We measured Young's modulus of the cured gels, E, by surface indentation. Specifically, we used a steel ball of radius R = 2.5 mm as an indenter. Additional weight was also applied on the top of the sphere to indent the gels. The indentation length d at a given indentation force F was measured from the side view. E was then calculated via the Hertzian contact law, E = 9F/(16R^1/2d^3/2)^47, where we assume the deformation of the steel ball is negligible and Poisson's ratio of the gel is n [?] 0.5. We also estimated E independently by matching the experimental stress distributions of solid-sphere impact with those from finite-element simulations. E measured using these two different methods match well with at most 20% difference from different experimental runs. The values also agree quantitatively with the previous studies^48. Poisson's ratio of the gels, n, is more difficult to assess accurately. We used n = 0.49 based on previous studies of the mechanical properties of PDMS gels under small strains^49. Surface stresses and displacements Two approximations were taken to estimate the surface stresses and displacements. First, we defined z = 0 as the horizontal plane through the lowest contact point between impactors and impacted substrates. As the maximum deformation of the PDMS gels under drop impact was only 92 mm comparable to the spatial resolution of DIC at 115 mm, we did not expect the approximation leads to large experimental errors. Second, to avoid any potential boundary-induced artifacts in DIC, we took the stresses and displacements at z = -315 mm below the surface as the surface stresses and displacements. We verified that the stress distributions at shallower heights show quantitatively similar features although noisier. Particularly, the locations of the maximum shear stress and pressure do not vary with z over this range. As a calibration, the procedure yielded a good approximation of the surface stress of solid-sphere impact at a larger maximum deformation of 192 mm (Fig. 2b). Finite element analysis To test the accuracy of high-speed stress microscopy (Fig. 2b), we used the commercial finite-element software ABAQUS to simulate solid-sphere impact, which has a shorter time scale than drop impact under similar impact conditions (Fig. 5a). The impact geometry was axisymmetric. The element shape for the meshing of the impacted surface was quadrilateral with adjustable sizes. Near the impact point, the element size was 0.2 mm. The diameter and the impact velocity of the impacting solid sphere were 3.16 mm and 0.49 m/s, matching the impact condition of the experiments. Since Young's modulus of the steel sphere is much larger than that of the impacted PDMS substrate, the sphere is assumed to be rigid in the simulation. The substrate is isotropic and linearly elastic. While Poisson's ratio of the substrate was fixed at 0.49, Young's modulus of the substrate was chosen to match the outcome of the experiments. Thus, the simulation allowed us to assess Young's modulus of the PDMS gels, independent of the indentation measurements discussed above. The two methods yield quantitatively similar results, agreeing well with the literature value. To numerically solve the coupled equations of the impact stress of impacting drops and the deformation of elastic media (Fig. 7c), we used the partial differential equation toolbox in Matlab. Specifically, we adopted the transient axisymmetric geometry in the structural mechanics analysis of the toolbox. The pressure and shear stress distributions of incompressible drops on infinitely rigid substrates (Eqs. (14) and (15) in SI) were assigned as the boundary loads on an elastic medium, which has Young's modulus and Poisson's ratio matching those of experiments. We then numerically calculated the radial and vertical displacements of the surface of the elastic medium as a function of time (Fig. 7c). The dimension and the grid size of the medium are chosen so that the results are convergent, independent of these parameters. A detailed discussion of the coupled differential equations and their numerical solutions for both drop impact and solid-sphere impact can be found in SI Sec. 3. Data availability The data supporting the main findings of this study are available in the paper and its Supplementary Information. Any additional data can be available from the corresponding author upon request. Code availability The codes that support the findings of this study are available from the corresponding author upon reasonable request. References 1. Laozi & Legge, J. Tao Te Ching (Oxford Univ. Press, 1891). 2. de Jong, R., Zhao, S.-C., Garcia-Gonzalez, D., Verduijn, G. & van der Meer, D. Impact cratering in sand: comparing solid and liquid intruders. Soft Matter 17, 120-125 (2021). ADS PubMed Google Scholar 3. Pimentel, D. et al. Environmental and economic costs of soil erosion and conservation benefits. Science 267, 1117-1123 (1995). ADS CAS PubMed Google Scholar 4. Kinnell, P. I. A. Raindrop-impact-induced erosion processes and prediction: a review. Hydrol. Process. 19, 2815-2844 (2005). ADS Google Scholar 5. Zhao, R., Zhang, Q., Tjugito, H. & Cheng, X. Granular impact cratering by liquid drops: understanding raindrop imprints through an analogy to asteroid strikes. Proc. Natl. Acad. Sci. USA 112, 342-347 (2015). ADS CAS PubMed Google Scholar 6. Zhang, Q., Gao, M., Zhao, R. & Cheng, X. Scaling of liquid-drop impact craters in wet granular media. Phys. Rev. E 92, 042205 (2015). ADS Google Scholar 7. Bako, A. N., Darboux, F., James, F., Josserand, C. & Lucas, C. Pressure and shear stress caused by raindrop impact at the soil surface: Scaling laws depending on the water depth. Earth Surf. Process. Landforms 41, 1199-1210 (2016). ADS Google Scholar 8. Erkal, A., D'Ayala, D. & Sequeira, L. Assessment of wind-driven rain impact, related surface erosion and surface strength reduction of historic building materials. Build. Environ. 57, 336-348 (2012). Google Scholar 9. Zhou, Q. et al. Liquid drop impact on solid surface with application to water drop erosion on turbine blades, part ii: Axisymmetric solution and erosion analysis. Int. J. Mech. Sci. 50 , 1543-1558 (2008). MATH Google Scholar 10. Gamero-Castano, M., Torrents, A., Valdevit, L. & Zheng, J.-G. Pressure-induced amorphization in silicon caused by the impact of electrosprayed nanodroplets. Phys. Rev. Lett. 105, 145701 (2010). ADS PubMed Google Scholar 11. Mitchell, B. R., Klewicki, J. C., Korkolis, Y. P. & Kinsey, B. L. Normal impact force of rayleigh jets. Phys. Rev. Fluids 4, 113603 (2019). ADS Google Scholar 12. Kondo, T. & Ando, K. Simulation of high-speed droplet impact against a dry/wet rigid wall for understanding the mechanism of liquid jet cleaning. Phys. Fluids 31, 013303 (2019). ADS Google Scholar 13. Alder, W. F. In Erosion (ed. Preece, C. M.), pp. 127-184 (Academic Press, 1979). 14. Lesser, M. B. & Field, J. E. The impact of compressible liquids. Annu. Rev. Fluid Mech. 15, 97-122 (1983). ADS Google Scholar 15. Worthington, A. M. On the forms assumed by drops of liquids falling vertically on a horizontal plate. Proc. Royal Soc. Lond. 25, 261-272 (1877). Google Scholar 16. Yarin, A. Drop impact dynamics: splashing, spreading, receding, bouncing.... Annu. Rev. Fluid Mech. 38, 159-192 (2006). ADS MathSciNet MATH Google Scholar 17. Josserand, C. & Thoroddsen, S. Drop impact on a solid surface. Annu. Rev. Fluid Mech. 48, 365-391 (2016). ADS MathSciNet MATH Google Scholar 18. Cheng, X., Sun, T.-P. & Gordillo, L. Drop impact dynamics: impact force and stress distributions. Annu. Rev. Fluid Mech. 54, 57-81 (2022). Google Scholar 19. Grinspan, A. S. & Gnanamoorthy, R. Impact force of low velocity liquid droplets measured using piezoelectric pvdf film. Colloids Surfaces A: Physicochem. Eng. Aspects 356, 162-168 (2010). Google Scholar 20. Soto, D., De Lariviere, A. B., Boutillon, X., Clanet, C. & Quere, D. The force of impacting rain. Soft Matter 10, 4929-4934 (2014). ADS CAS PubMed Google Scholar 21. Li, J., Zhang, B., Guo, P. & Lv, Q. Impact force of a low speed water droplet colliding on a solid surface. J. Appl. Phys. 116, 214903 (2014). ADS Google Scholar 22. Zhang, B., Li, J., Guo, P. & Lv, Q. Experimental studies on the effect of reynolds and weber numbers on the impact forces of low-speed droplets colliding with a solid surface. Exp. Fluids 58 , 125 (2017). Google Scholar 23. Gordillo, L., Sun, T.-P. & Cheng, X. Dynamics of drop impact on solid surfaces: evolution of impact force and self-similar spreading. J. Fluid Mech. 840, 190-214 (2018). ADS CAS Google Scholar 24. Mitchell, B. R., Klewicki, J. C., Korkolis, Y. P. & Kinsey, B. L. The transient force profile of low-speed droplet impact: measurements and model. J. Fluid Mech. 867, 300-322 (2019). ADS CAS Google Scholar 25. Zhang, R., Zhang, B., Lv, Q., Li, J. & Guo, P. Effects of droplet shape on impact force of low-speed droplets colliding with solid surface. Exp. Fluids 60, 64 (2019). Google Scholar 26. Schmid, G. et al. On the measurement and prediction of rainfall noise. Appl. Acoust. 171, 107636 (2021). Google Scholar 27. Thanh-Vinh, N., Matsumoto, K. & Shimoyama, I. Pressure distribution on the contact area during the impact of a droplet on a texture surface. In IEEE 29th International Conference on Micro Electro Mechanical Systems (MEMS), pp. 177-180 (IEEE, 2016). 28. Thanh-Vinh, N. & Shimoyama, I. Maximum pressure caused by droplet impact is dependent on the droplet size. In 20th International Conference on Solid-State Sensors, Actuators and Microsystems and Eurosensors XXXIII, pp. 813-816 (IEEE, 2019). 29. Gerber, J., Lendenmann, T., Eghlidi, H., Schutzius, T. M. & Poulikakos, D. Wetting transitions in droplet drying on soft materials. Nat. Commun. 10, 4776 (2019). ADS PubMed PubMed Central Google Scholar 30. Hall, M. S., Long, R., Hui, C.-Y. & Wu, M. Mapping three-dimensional stress and strain fields within a soft hydrogel using a fluorescence microscope. Biophys. J. 102, 2241-2250 (2012). ADS CAS PubMed PubMed Central Google Scholar 31. Abdulagatov, I. M., Zeinalova, A. B. & Azizov, N. D. Viscosity of aqueous electrolyte solutions at high temperatures and high pressures. viscosity b-coefficient. sodium iodide. J. Chem. Eng. Data 51, 1645-1659 (2006). CAS Google Scholar 32. Chen, H., Li, Z., Wang, F., Wang, Z. & Li, H. Investigation of surface properties for electrolyte solutions: measurement and prediction of surface tension for aqueous concentrated electrolyte solutions. J. Chem. Eng. Data 62, 3783-3792 (2017). CAS Google Scholar 33. Riboux, G. & Gordillo, J. M. Experiments of drops impacting a smooth solid surface: a model of the critical impact speed for drop splashing. Phys. Rev. Lett. 113, 024507 (2014). ADS PubMed Google Scholar 34. Wildeman, S., Visser, C. W., Sun, C. & Lohse, D. On the spreading of impacting drops. J. Fluid Mech. 805, 636-655 (2016). ADS MathSciNet CAS Google Scholar 35. Rioboo, R., Marengo, M. & Tropea, C. Time evolution of liquid drop impact onto solid, dry surfaces. Exp. Fluids 33, 112-124 (2002). Google Scholar 36. Mongruel, A., Daru, V., Feuillebois, F. & Tabakova, S. Early post-impact time dynamics of viscous drops onto a solid dry surface. Phys. Fluids 21, 032101 (2009). ADS MATH Google Scholar 37. Visser, C. W. et al. Dynamics of high-speed micro-drop impact: numerical simulations and experiments at frame-to-frame times below 100 ns. Soft Matter 11, 1708-1722 (2015). ADS CAS PubMed Google Scholar 38. Howland, C. J. et al. It's harder to splash on soft solids. Phys. Rev. Lett. 117, 184502 (2016). ADS PubMed Google Scholar 39. Zhang, J. M., Li, E. Q. & Thoroddsen, S. T. Fine radial jetting during the impact of compound drops. J. Fluid Mech. 883, A46 (2020). ADS Google Scholar 40. Philippi, J., Lagree, P.-Y. & Antkowiak, A. Drop impact on a solid surface: short-time self-similarity. J. Fluid Mech. 795, 96-135 (2016). ADS CAS MATH Google Scholar 41. Smith, F. T., Li, L. & Wu, G. X. Air cushioning with a lubrication/inviscid balance. J. Fluid Mech. 482, 291-318 (2003). ADS MathSciNet MATH Google Scholar 42. Mandre, S., Mani, M. & Brenner, M. P. Precursors to splashing of liquid droplets on a solid surface. Phys. Rev. Lett. 102, 134502 (2009). ADS PubMed Google Scholar 43. Mani, M., Mandre, S. & Brenner, M. P. Events before droplet splashing on a solid surface. J. Fluid Mech. 647, 163-185 (2010). ADS MathSciNet CAS MATH Google Scholar 44. Freund, L. B. Dynamic Fracture Mechanics (Cambridge Univ. Press, 1990). 45. Li, E. Q. & Thoroddsen, S. T. Time-resolved imaging of a compressible air disc under a drop impacting on a solid surface. J. Fluid Mech. 780, 636-648 (2015). ADS Google Scholar 46. Li, E. Q., Langley, K. R., Tian, Y. S., Hicks, P. D. & Thoroddsen, S. T. Double contact during drop impact on a solid under reduced air pressure. Phys. Rev. Lett. 119, 214502 (2017). ADS PubMed Google Scholar 47. Landau, L. D. & Lifshitz, E. M. Theory of Elasticity, 3rd edn (Butterworth-Heinemann, 1986). 48. Park, J. Y. et al. Increased poly(dimethylsiloxane) stiffness improves viability and morphology of mouse fibroblast cells. BioChip J. 4, 230-236 (2010). CAS Google Scholar 49. Johnston, I. D., McCluskey, D. K., Tan, C. K. L. & Tracey, M. C. Mechanical characterization of bulk sylgard 184 for microfluidics and microengineering. J. Micromech. Microeng. 24, 035017 (2014). ADS CAS Google Scholar Download references Acknowledgements We thank Ben Druecke, Tai-yin Chiu, and Grace Lee for help with experiments and data analysis and Michelle Driscoll for fruitful discussions. This research is supported by the US National Science Foundation CBET-2017071 and 2002817 and ACS Petroleum Research Fund 60668-ND9. T.-P.S. acknowledges the partial financial support of the PPG fellowship via UMN IPRIME and the Government Scholarship to Study Abroad from Taiwan. F.A.-N., K.A., L.G., and P.G. acknowledge the financial support of the grants ANID/CONICYT Fondecyt Iniciacion No. 11170700 and 11191106. Author information Authors and Affiliations 1. Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, MN, 55455, USA Ting-Pi Sun & Xiang Cheng 2. Departamento de Fisica, Facultad de Ciencia, Universidad de Santiago de Chile (USACH), Santiago, Chile Franco Alvarez-Novoa, Klebbert Andrade & Leonardo Gordillo 3. Instituto de Ciencias de la Ingenieria, Universidad de O'Higgins, Rancagua, Chile Pablo Gutierrez Authors 1. Ting-Pi Sun View author publications You can also search for this author in PubMed Google Scholar 2. Franco Alvarez-Novoa View author publications You can also search for this author in PubMed Google Scholar 3. Klebbert Andrade View author publications You can also search for this author in PubMed Google Scholar 4. Pablo Gutierrez View author publications You can also search for this author in PubMed Google Scholar 5. Leonardo Gordillo View author publications You can also search for this author in PubMed Google Scholar 6. Xiang Cheng View author publications You can also search for this author in PubMed Google Scholar Contributions T.-P.S. and X.C. designed the research and developed the high-speed stress microscopy. T.-P.S. performed the impact stress measurements. T.-P.S., L.G., and X.C. discussed and analyzed experimental data. L.G. solved the numerical solution of drop impact on elastic media with input from X.C. F.A.-N., K.A., P.G., and L.G. studied drop-impact erosion on plaster slabs. X.C. conceived and supervised the project. T.-P.S. and X.C. cowrote the manuscript. All authors discussed and commented on the manuscript. Corresponding author Correspondence to Xiang Cheng. Ethics declarations Competing interests The authors declare no competing interests. Peer review Peer review information Nature Communications thanks Emile van der Heide and the other, anonymous, reviewers for their contribution to the peer review of this work. Peer reviewer reports are available. Additional information Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary information Supplementary Information Peer Review File Rights and permissions Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. Reprints and Permissions About this article Verify currency and authenticity via CrossMark Cite this article Sun, TP., Alvarez-Novoa, F., Andrade, K. et al. Stress distribution and surface shock wave of drop impact. Nat Commun 13, 1703 (2022). https://doi.org/10.1038/s41467-022-29345-x Download citation * Received: 21 August 2021 * Accepted: 07 March 2022 * Published: 31 March 2022 * DOI: https://doi.org/10.1038/s41467-022-29345-x Share this article Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative Comments By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Download PDF Associated Content Focus Applied physics and mathematics Advertisement Advertisement Explore content * Research articles * Reviews & Analysis * News & Comment * Videos * Collections * Subjects * Follow us on Facebook * Follow us on Twitter * Sign up for alerts * RSS feed About the journal * Aims & Scope * Editors * Journal Information * Open access * Editorial Values Statement * Journal Metrics * Editors' Highlights * Contact * Editorial policies * Top Articles Publish with us * For authors * For Reviewers * Submit manuscript Search Search articles by subject, keyword or author [ ] Show results from [All journals] Search Advanced search Quick links * Explore articles by subject * Find a job * Guide to authors * Editorial policies Nature Communications (Nat Commun) ISSN 2041-1723 (online) nature.com sitemap Nature portfolio * About us * Press releases * Press office * Contact us * * * Discover content * Journals A-Z * Articles by subject * Nano * Protocol Exchange * Nature Index Publishing policies * Nature portfolio policies * Open access Author & Researcher services * Reprints & permissions * Research data * Language editing * Scientific editing * Nature Masterclasses * Nature Research Academies * Research Solutions Libraries & institutions * Librarian service & tools * Librarian portal * Open research * Recommend to library Advertising & partnerships * Advertising * Partnerships & Services * Media kits * Branded content Career development * Nature Careers * Nature Conferences * Nature events Regional websites * Nature Africa * Nature China * Nature India * Nature Italy * Nature Japan * Nature Korea * Nature Middle East Legal & Privacy * Privacy Policy * Use of cookies * Manage cookies/Do not sell my data * Legal notice * Accessibility statement * Terms & Conditions * California Privacy Statement Springer Nature (c) 2022 Springer Nature Limited Close Nature Briefing Sign up for the Nature Briefing newsletter -- what matters in science, free to your inbox daily. Email address [ ] Sign up [ ] I agree my information will be processed in accordance with the Nature and Springer Nature Limited Privacy Policy. Close Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing * *