News center
Premium-quality after-sales assistance

Hydrodynamic Kelvin

Dec 11, 2023

Scientific Reports volume 13, Article number: 2686 (2023) Cite this article

463 Accesses

Metrics details

Kelvin–Helmholtz instability on metallic surface is relevant to intense oblique impact in many physical processes such as explosive welding, Inertial Confinement Fusion and planetary impact events. Evolution of instability results in the formation of wavy morphology leading to material bonding or even mixing. However, mostly due to lack method to describe the dynamic behavior, instability mechanism controlled by elastoplastic properties of metal remains elusive. Here, we introduce a theory to reveal the evolution characteristics aroused by tangential velocity. Our simulations find that the unstable metallic surfaces exhibit amplitude growth and tangential motion by overcoming the depression of yield strength to generate wavy morphology. For diverse loading velocities, corrugated surfaces and material properties, an instability boundary distinguishes all unstable evolutions. Our analytical method with scale-independent variables reproducing numerical findings reveals plentiful characteristics of instability in strength materials. For designed loading velocities and material in oblique impact experiment in laboratory, the property of corrugated surfaces becomes an important factor to determine instability evolution.

Kelvin–Helmholtz instability(KHI)1,2 due to shear at metallic surface remains scarcely understood, which especially deserve to interpret as metal suffering intense oblique impact in High Velocity Impact Welding (HVIW)3,4,5, Inertial Confinement Fusion (ICF)6,7, planetary impact events8,9,10, etc. The wavy structures aroused by tangential velocity jump at the instance of surface collision with angles indicate material bonding or even potential mixing5,8. Although KHI between fluids are studied extensively11,12, characteristics of KHI evolution associated with depression effects of elastic–plastic (EP) properties of metal13 merit thoroughly digging.

Detection of KHI on metallic surface is a serious challenge due to practical difficulties of sustainment of high-speed shear flow in experimental facilities14. The characteristics of wavy morphologies are usually discussed with the help of high velocity oblique impact experiment whose results can only be imaged at the end of experiments which do not reveal the evolution processes3,4,15,16, not mentioning another issue of recovering sample without severe fragmentations under high velocity loading17. Although oblique impact processes can be exhibited by computer simulations, besides acquiring adequately fine mesh distributions, the accuracy of calculations is largely determined by different arithmetic of capturing material interface15,18,19,20. For KHI on metals, it is surprising that relevant simulations have not been shown for now but only theoretical analysis with traditional normal mode method which merely presents growth rate and involves impossibility of analytical treatments due to nonlinear governing equations and nonlinear constitutive relations of metal5,18. As a result, we especially lack the descriptions of evolution characteristics of metallic perturbed surface under the operation of tangential velocity discontinuity.

For the purpose of investigating surface behavior of KHI on solid, we have proposed a theoretical analysis with a potential flow method to describe the growth rate and amplitude evolution by analytical formulas21. The properties of resisting shear deformation of solid material influence the instability evolution of surface as flushed by tangential flow. The amplitude growth is prevented by EP properties of solid to become behavior of oscillating around. Although the depression effect of EP properties has been detected in the amplitude evolution, it is interesting that the growth rate is the same as KHI for different ideal fluids, i.e. \(k\sqrt{{\rho }_{1}{\rho }_{2}{u}_{0}^{2}}/({\rho }_{1}{+\rho }_{2})\), which is always positive to indicate continuous growth of amplitude. The traditional method to estimate whether the surface is stable or unstable by growth rate18,22 seems invalid for solid. Besides, the relation between EP transition and instability evolution can not also be exhibited by growth rate and amplitude. In present work, we attempt to illuminate a method to estimate whether instability develop named instability boundary and to explain the effect of EP transition on instability by EP division.

Here, we consider the instability for the configuration of an ideal fluid with constant tangential velocity u0 gliding over a quiescent perfectly EP solid (Fig. 1). For simplicity our discussion is restricted in two-dimensional plane with y-axis perpendicular to flow direction x. The small perturbation can be represented by η(x,t) = ξ(t)eikx where ξ(t) is amplitude with initial value of ξ(0) = ξ0 and k = 2π/λ is wave number for wavelength λ. The metal and fluid may commonly have different densities of ρ1 and ρ2. In this system, u0 is the inducement of instability, yet the surface may be stable under the suppression of constant shear modulus G1 and constant yield stress Y before and after plastic deformation. Then several dimensionless variables which characterize a KHI system are defined: AT = (ρ1 - ρ2) / (ρ1 + ρ2) is Atwood number; M02 = ρ1u02/G1 is Mach number; z = ξ(t)/ξ0 is growth factor; τ = tku0, \(\widehat{\lambda }\hspace{0.17em}\)= 2πξ0/λ and \(\widehat{Y}\hspace{0.17em}\)= ρ1u02/Y are respectively dimensionless time, wavelength and yield strength.

The configuration of ideal fluid flowing over a perturbed surface of perfectly EP metal. The solid with density ρ1 is quiescent in the 2D Cartesian coordinate system and the fluid with density ρ2 has a constant tangential velocity u0 in x direction. The initial perturbation has a cosinusoidal form ξ0coskx with periodic wavelength λ and amplitude ξ0. Both thicknesses h1 and h2 of materials are large enough to make sure approaching infinite medium, i.e. kh1 >> 1 and kh2 >> 1.

We first execute extensive numerical calculations to reveal the time evolutions of metallic perturbations. The simulations perform the generation of various wavy morphologies displaying continuous growth of amplitude and tangential motion of wave crest which implies system instability. Based on the characteristics of instability evolution, an instability boundary separating all stable and unstable cases is achieved, and EP division splitting all cases without and with plastic behaviors is also obtained. Furthermore, a theoretical method is carried out to disclose analytical expressions of instability boundary and EP division which are evaluated by simulation findings. The identical results show that the seemingly complex instability can be quantitatively described by the defined scale-independent variables to elucidate the characteristics of instability evolution by the suppression of shear modulus and yield strength from a straightforward perspective. That makes the theory a potential versatile method to estimate similar KHI system on any other strength medium whether the surface is unstable to form wavy patterns after high intense oblique impact.

The KHI configuration in Fig. 1 is simulated by finite element method23 (see "Methods" section). At initial moment, water with constant tangential velocity closely contacts to perturbed copper surface. Scales of variables are chosen by consulting oblique impact experiments24,25.

After abundant calculations, we performed four typical time evolutions of metallic surfaces including stable and unstable cases. Due to the depression of shear modulus G1 and strength Y, the stable surface maintains similarity to the initial perturbation as time passing (Fig. 2a), and the amplitude presents oscillating behavior around a small range in vertical direction (Fig. 2e). Under the operation of velocity in x direction, the unstable evolutions perform y direction growth accompanying tangential direction motion such as slight tangential shifting (Fig. 2b), visible tangential movement (Fig. 2c), or surface rolling up (Fig. 2d). The growth factors of unstable wavy patterns all exhibit the trends of increasing with different rates (Fig. 2e). By the way, the wavy morphologies of unstable surfaces in our simulations at time 0.15 μs and 0.2 μs (Figs. 2b–d) are similar to those observed in oblique impact experiments, including different tangential movement behavior and curling structure24,25,26,27,28.

Four cases of surface temporal evolutions. ( a–d) Maps of surface temporal morphologies of quiescent Cu plate surface (blue) flowed by ideal fluid H2O (grey) with tangential velocity u0, as calculated by finite element methods23. Each Cu plate has fixed density ρ1 = 8.9 kg/m3, fixed shear modulus 39.39 GPa and the same initial perturbation with 250 μm wavelength and 10 μm amplitude. The four cases are obtained by varying tangential velocity, fluid densities and yield strength of Cu for case a u0 = 1.0 mm/μs, ρ2 = 2.0 kg/m3, Y = 500 MPa, case (b) u0 = 1.0 mm/μs, ρ2 = 3.0 kg/m3, Y = 500 MPa, case (c) u0 = 2.0 mm/μs, ρ2 = 1.0 kg/m3, Y = 500 MPa and case (d) u0 = 2.0 mm/μs, ρ2 = 1.0 kg/m3, Y = 100 MPa. Surface morphologies of four cases at time 0.05 μs, 0.1 μs, 0.15 μs, 0.2 μs are given respectively. The red arrow in each surface at time 0.2 μs is the schematic of surface motion direction. (e) The growth factor z of each case is extracted from simulation.

Due to obvious disparity of amplitude development, we find a boundary to partition stable and unstable evolutions for all variable combinations (Fig. 3). Different dimensionless variables AT, M0, \(\widehat{\lambda }\) and \(\widehat{Y}\) in simulations are obtained by changing fluid density, tangential velocity, initial amplitude, wavelength, shear modulus and yield strength respectively (see Supplementary). According to stable oscillation and continuous increase of growth factors, the boundary is achieved by fixing \(\widehat{\lambda }\) and varying \(\widehat{Y}\) with enough point densities to approach the position of instability margin. Then, adopting the same procedure to more \(\widehat{\lambda }\) values, a line can be drawn to split the domains. As shown in Fig. 3 with different AT and M0, the areas containing all variable combinations below and above instability boundary respectively indicate stable and unstable surface. Similarly, EP division partitioning amplitude motion with and without plastic behavior is also achieved. The area below EP division implies elastic deformation, as plastic behavior above the division line. We also observe that EP division is beneath instability boundary hinting that continuous growth of perturbation must experience plastic transformation to overcome the resistance of deforming from the initial perturbation.

Instability boundary and EP division by simulations and theory. Three combinations of AT and M0 are calculated i.e. (a) AT = 0.7980 and M0 = 0.4754, (b) AT = 0.4958 and M0 = 0.3803, (c) AT = 0.0 and M0 = 0.2377. The axis 2πξ0/λ and ρ1u02/Y in the figures are \(\hat{\lambda }\) and \(\hat{Y}\). Solid point and hollow point in each figure represents simulation results of instability boundary and EP division respectively. Solid line and dash line means results by theory for instability boundary and EP division. The results with identical AT and M0 by numerical simulation and analytical theory are plotted in one figure for comparison.

Instability analysis starts from the governing equations of continuity and momentum with the assumptions of incompressible and irrotational flow. Potential flow method is adopted to represent the velocity field which should be continuous in the normal direction of material interface with the perturbation of η(x,t) = ξ(t)eikx. Then, the motion equation for describing evolution of interface amplitude ξ(t) is achieved with the condition of force equilibrium in the normal direction of material interface. We give a specific case of the process of establishing the instability analysis in "Methods" section. With perfectly EP properties of solid and Cauchy stress of viscous fluid, the motion equation of interface between EP solid and viscous fluid is obtained. The instability of solid is focused in the present study to ignore the viscosity of fluid. The motion equation of amplitude is changed to be a dimensionless form (Eq. (21) in "Methods" section)

where zp is growth factor at when EP transition takes place and

The symbol Λ and Χ contain influences of shear modulus and yield strength represented by dimensionless variables M0 and \(\hat{Y}\). In Eq. (1) the first branch controls amplitude motion with elastic behavior before growth factor z arrives at zp beyond which plastic deformation occurs to depict the amplitude behavior by the second branch. Starting from Eq. (1) with mathematical deriving (see "Methods" section), we disclose the analytical formulas of instability boundary (Eq. (32) in "Methods" section)

and EP division (Eq. (38) in "Methods" section)

The specific formulations of xp ,τp and τe (see Eqs. (30), (28) and (36) in "Methods" section) are relative to M1, M2 and Λ, therefore instability boundary and EP division with the relation of \(\widehat{Y}\)= f (\(\widehat{\lambda }\)) are determined by AT and M0. Lines calculated by Eqs. (3) and (4) with same dimensionless variables as simulations are also plotted in Fig. 3 which show identical results.

Furthermore, we perform more properties about instability boundary, EP division and amplitude evolution (solutions of Eq. (1) in "Methods" section) by our theory to understand the instability behavior.

Instability boundary and EP division divide \(\widehat{Y}\) vs \(\hat{\lambda }\) field into three parts. Some points in each part are picked to perform amplitude time evolutions (Fig. 4). Figure 4a plots instability boundary and EP division for AT = 0.5 and M0 = 0.4. The area below EP division means amplitude motion only with elastic behavior whose growth factor is also determined by AT and M0 (Eq. (39a) in "Methods" section), therefore all parameter combinations in the area below EP division possesses the same time evolution which is controlled by G1 leading to vibrate in a small range (Fig. 4b).The area between instability boundary and EP division implies the amplitude is stable with plastic deformation. The growth factor initially elastically oscillate to exceed zp to plastic stage, then is suppressed by Y to a maximum value and oscillates around (Eq. (39b) in "Methods" section). For the regions above instability boundary, the growth factor also elastically vibrates to plastic stage, but the effect of Y is evidently weak as shown continuous increase of amplitude in Fig. 4b (Eq. (39c) in "Methods" section). Except for the features illustrated above, some others are also detected, such as growth factor exhibit more stable motion as \(\widehat{Y}\) or \(\widehat{\lambda }\) lessen, and identical \(\widehat{Y}\widehat{\lambda }\) specifies the same evolution in plastic stage (Eqs. (39b) and (39c) in "Methods" section). The corresponding values of \(\widehat{\lambda }\) and \(\widehat{Y}\), EP state and surface state are summarized in Tables 1 and 2 for Fig. 4. Besides, we also give instability boundary, EP division and amplitude evolution for another AT = 0.9 and M0 = 0.2 (Fig. 4c,d). The regularity of surface development are similar to those with AT = 0.5 and M0 = 0.4 except that the range of oscillation is obliviously narrower.

Two groups of instability boundary, EP division and growth factor by theory. (a,b) The first group of instability boundary, EP division and growth factor for arbitrarily picked points with AT = 0.5, M0 = 0.4. (c,d) The second group with AT = 0.9, M0 = 0.2. The solid points meaning picked variable combinations locate above instability boundary, and the hollow points are between instability boundary and EP division. \(\hat{Y}\) and \(\hat{\lambda }\) values of each point are marked in bracket after point. Some of them have the same \(\hat{\lambda }\), some have the same \(\hat{Y}\). For the points with the same shape and color, their \(\hat{Y}\hat{\lambda }\) is identical. The color of the curve of growth factor z in (b) and (d) figures are completely the same as point color. The line labeled as "Elastic area" is the growth factor for the region beneath EP division, and other lines are labeled with corresponding \(\hat{Y}\) and \(\hat{\lambda }\) values in bracket.

Due to AT and M0 deciding instability boundary and EP division, we also illustrate the influences. By decreasing AT gradually for M0 = 0.2 (Fig. 5a), instability boundary and EP division are simultaneously moving close to the coordinate axis which show the area of amplitude growth and plastic motion are enlarged. For a fixed combination of \(\hat{Y}\) and \(\hat{\lambda }\) which is designated in Fig. 5a, as varying AT the point locates at different areas. For AT = 0.9, the point is in the elastic area and the corresponding amplitude elastically oscillates (black line in Fig. 5b). As AT = 0.7, the point is between instability boundary and EP division and the growth is controlled by yield stress (green line in Fig. 5b). After AT = 0.3, the point is in the instability area to exhibit a continuous increasing amplitude (blue line in Fig. 5b). Besides, we also find that as AT decrease the area between instability boundary and EP division diminishes which indicates that the surface is unstable once plastic transformation occurs for small AT. Then, the features of instability boundary and EP division affected by M0 are plotted in Fig. 5c,d. The influences of M0 are not that evident especially on EP division and as increasing M0 instability boundary presents appreciable movement in the direction of approaching to axis. For a point between instability boundary of M0 = 0.6 and 0.8, the growth factor must undergo plastic deformation, and the amplitude is stable for M0 = 0.6 (black line in Fig. 5d) and unstable for M0 = 0.8 (blue line in Fig. 5d). It seems that the state of amplitude evolution of variable combinations near the instability boundary is sensitive to the change of M0.

Influences of AT and M0 on instability boundary and EP division. (a) For a fixed M0 = 0.2, four groups of instability boundary (solid line) and EP division (dash line) with AT = 0.9, 0.7, 0.3, − 0.5 are plotted. A point between instability boundary and EP division of AT = 0.7 is designated with red star. (b) The growth factor of red star in figure (a) as AT varying. (c) For a fixed AT = 0.9, five groups of instability boundary (solid line) and EP division (dash line) with M0 = 0.2, 0.4, 0.6, 0.8, 1.0 are shown. Also a red star between M0 = 0.6 and 0.8 is designated. (d) The growth factor of red star in figure (c) as M0 varying.

The prominent characteristic of solid is the ability to resist shear deformation. Solid body deforms when it is subjected to external forces. If the deformation is reversible, i.e. the deformation vanishes instantaneously as soon as external forces are removed, this type of deformation is elastic. If the deformation is permanent, i.e. solid yields, solid deforms with plastic behavior. The resistance to shear deformation obviously determines the instability of solid surface which are shown in the above results.

Those results show that a whole picture of KHI in solid are becoming clear through instability boundary, EP division and evolution of growth factor which are denoted by dimensionless variables AT, M0, \(\widehat{\lambda }\) and \(\widehat{Y}\). The instability boundary and EP division are built in the plane of \(\widehat{\lambda }\) and \(\widehat{Y}\) which are divided into three parts. The relative position of instability boundary and EP division is that instability boundary locates above EP division in \(\widehat{\lambda }\) and \(\widehat{Y}\) plane which indicates physical insights of the characteristics of instability evolution. The first part in \(\widehat{\lambda }\) and \(\widehat{Y}\) plane is the region below EP division in which solid does not reach the yield point and the surface growth is restrained by shear modulus. The amplitude shows vibration around a small range due to resisting shear deformation in elastic stage. The morphology of surface may be reversible to original state if the inducement of instability disappears. The second part in \(\widehat{\lambda }\) and \(\widehat{Y}\) plane is the region above EP division and beneath instability boundary in which yield occurs and the surface growth is controlled by strength. In the second part, elasticity of solid can not suppress amplitude growth and deformation becomes large to plastic stage. Yet, the continuous growth of amplitude is restrained by yield strength after EP transition. Therefore, the characteristics of amplitude evolution contain three processes including elastic vibration to plastic deformation, growth suppressed by strength and vibration around. The morphology of surface may be permanent if the inducement of instability disappears. Those two parts constitute all cases of stable surfaces. The third part in \(\widehat{\lambda }\) and \(\widehat{Y}\) plane is the region above instability boundary in which surfaces are unstable with plastic deformation. In the third part, the inducement of instability drives the surface deforming from elastic stage to plastic stage, and strength can not prevent increase of amplitude. The surface evolution exhibits large deformation to form wavy morphology. Although, instability boundary and EP division are affected by AT and M0 with different degrees, the above physical regularity of instability evolution does not change. According to the analysis, it is found that EP transition is the necessary condition to amplitude growth to form wavy morphology, which imply that the minimal energy required to engender effective bonding or mixing between materials after oblique impact must exceed the energy leading to plastic deformation.

For KHI system we discussed above, the simulations demonstrate that the continuous growth of amplitude of metallic surface in early stage by overcoming the depression effects of EP properties is a signature of dynamic instability to form wavy morphology, and there is a boundary to distinguish all stable and unstable variable combinations for diverse velocities, perturbation geometries (initial amplitude and wavelength) and material properties (densities, shear modulus and yield strength).

Instability boundary and EP division predicted by our analytical formulas are intriguingly consistent with those obtained by simulations, and the model identifies dimensionless variables AT, M0, \(\widehat{\lambda }\) and \(\widehat{Y}\) as characteristic parameters to describe surface instability. The analytical results also quantitatively illustrate how the evolution processes are affected by those dimensionless variables from plentiful perspectives. Since our mathematical formulas do not include scale-independent terms, the model appears to be able to predict evolutions for broader KHI systems not only cases we discussed here.

Comparing the late-time morphology by our simulations with that after oblique impact experiment24,25,26,27,28, the characteristics of metallic wavy pattern are similar including amplitude growth and tangential motion which reproves that indeed the formation mechanism of wavy interface is the evolution of KHI5,25. The configuration in our study resembles the situation of a strength material collided with angle by a metal undergoing thermal softening18. Tangential velocity is decomposed from impact velocity according to impact angle between impactor and target. By extracting perturbation (e.g. surface roughness)24,25 and material properties of impactor and target, dimensionless variables of the configuration are easily calculated. Smaller AT and larger M0 mean an impactor with larger density and velocity for the same target. Smaller \(\widehat{\lambda }\) and \(\widehat{Y}\) corresponds to a higher strength target with more gentle perturbation. We can predict the surface evolution state through locating the position between calculated (\(\widehat{\lambda }\),\(\widehat{Y}\)) and instability boundary above which means wavy morphology formation. Thus, the presented theoretical method offers a potential versatile tool over a broad range of scales to help estimate the evolution of metallic surface after oblique impact to explain collisional bonding or even mixing29,30. For a more comprehensive evolution between two EP metals, the characteristics require further investigations.

The common assumption of material constitutive properties recognizing scaling strength linearly with the shear modulus31,32 and the constitutive behavior at high pressures33,34 indicate the connection between shear modulus and yield strength. For a designed condition of dynamic loading (tangential velocity) and samples (densities, shear modulus and yield strength), AT , \(\widehat{Y}\), M0 except \(\widehat{\lambda }\) (initial amplitude and wavelength) can be roughly estimated, namely the instability state of surface are determined by the properties of perturbation including origin and dimension which are another complex problems referring to machining techniques, material microstructures etc. The wavy morphology with plastic deformation after oblique impact takes as evidence of bonding between dissimilar materials29,30, therefore perturbation becomes a crucial factor of effective bonding.

Numerical simulation taken in this paper is a 2D Lagrange finite element method which has been adopted to simulate dynamic impact and surface instability frequently18,24. The advantage of Lagrange method of capturing material interface is utilized to obtain a clear surface of solid. The surfaces of materials are initially contact and defined by a sliding_only contact which is a two-surface method. The details of framework and reliability of the Lagrange method is performed in previous study21.

The solid and fluid are both simulated with EOS of the Mie-Grüneisen state with a coefficient γ = ρ0γ0/ρ where γ0 is a parameter characteristic and ρ0 is the initial density. The relation between shock velocity vs and particle velocity vp is vs = c0 + svp where c0 is the bulk sound velocity and s is a characteristic constants. For copper, ρ0 = 8.9 g/cm3, γ0 = 2.02, c0 = 3.94 cm/μs, and s = 1.49 are used, while ρ0 = 1.0 g/cm3, γ0 = 0.4934, c0 = 1.48 cm/μs, and s = 2.56 are used for water35,36. In order to be consistent with the theoretical model, a perfectly elastic and rigid plastic model is adopted to characterize solid with a constant shear modulus G1 and a constant yield strength Y.

The amplitude 2ξ0 is the distance from wave crest to trough in y direction. The length of the interface in x direction contains twenty wavelengths to diminish the effect of lateral extents. Square meshes with side length 2.5 μm are distributed at initial moment. The solid plate is quiescent and the tangential velocity of the fluid plate is set to be u0 at the initial time.

For a two-dimensional configuration, the amplitude motion equation is derived starting from the governing equations of continuity and momentum without conservative and nonconservative forces acting on the interface

where i = 1 and i = 2 represents two materials respectively, ui characterizes the irrotational perturbed velocity of material i, and pi is pressure. Potential flow theory is adopted to establish the instability analysis. For irrotational flow at initial time and considering the tangential velocities, we have

where ϕi and Φi satisfies Laplace equation. The mathematical formula of pressure can be obtained by integrating Eq. (6) from y = 0 to the instantaneous interface y = η (x, t) in y direction with C1 = C221

With kinematic conditions in the normal direction

and potential function

we determine the coefficient Ai(t) in terms of ξ(t) by considering the 2D coordinates system fixed on solid

The force equilibrium is also acquired at the interface in normal direction

where Fy(j) represent force per unitary surface acting on interface. For system of EP solid and viscous fluid, the upper formula can be expressed as

where S1,yy(ep) represents the vertical component of the deviatoric stress for EP solid and S2,yy(v) is the vertical component of the deviatoric part of the Cauchy stress tensor σij = -pδij + Sij for fluid.

The solid in elastic stage is assumed to be Hookean solid having linear constitutive relation37

where D1,ij is the strain rate tensor. The deviatoric stress tensors of fluid have the form of

where μ2 is dynamic viscosity. The vertical component of the deviatoric stress for elastic solid and viscous fluid are obtained with Eqs. (7), (11), (12), (15a), (15b), (16a) and (16b)

Then, the motion of amplitude for elastic equation can be described by Eq. (18) with substituting Eqs. (17a) and (17b) into Eq. (14)

The perturbation amplitude ξp at when solid yields is determined by the condition of EP transition, i.e. the effective stress \(\tilde{\sigma }=\sqrt{3{S}_{1,ij}{S}_{1,ij}/2}\) arrives at Y. Then, the mathematical expression for plastic solid is derived with Eqs. (15a) and (15b) and with the fact that the plastic deformation only occurs at a small layer away form surface such as y ~ k-1. The integrated motion equation of amplitude for EP solid and viscous fluid is

In order to focus on the instability on solid surface, viscosity effect of fluid is not discussed temporarily. The amplitude motion of KHI between EP solid and ideal fluid is described by taking μ2 = 0 in Eq. (19)

whose dimensionless form is

with the definition of dimensionless variables AT = (ρ1-ρ2) / (ρ1 + ρ2), M02 = ρ1u02/G1, z = ξ(t)/ξ0, τ = tku0, \(\widehat{\lambda }\)= 2πξ0/λ and \(\widehat{Y}\)= ρ1u02/Y. zp is growth factor at when EP transition takes place and

\(\widehat{\lambda }\) represents the characteristic of interface and \(\widehat{Y}\) denotes the inducement and resistance on instability.

The condition of stability is that the amplitude must have a maximum value at a certain time τ = τm which implies

The deriving are started from Eq. (21) with initial conditions z(0) = 1 and ż(0) = 0. The transformations

are introduced into Eq. (21) to achieve

and the initial conditions and continuous conditions are

where τp is the time when the solid transients from elasticity to plasticity. Integrating Eq. (25a) with Eq. (26a) and evaluating Eqs. (26b) and (26c) at time τp when the marginal stable state transits to plastic regime, it has

Then, performing integration of Eq. (25a) twice with Eq. (26a) and (26b),we get

By taking a first integration on Eq. (25b) with Eq. (26b) and (26c) and evaluating x2(τm) = ẋ2(τm) = 0 due to Eq. (23), we have

Then, by carrying the first derivatives of Eqs. (24a) and (24b) and evaluating them at τ = τp, and combining Eqs. (27) and (29), we obtain xp

Evaluating Eq. (26b) into Eq. (24a) and combining zp in Eq. (22), it has

With the definition of Χ in Eq. (22), the instability boundary is achieved

EP transition happens when the maximum amplitude zme = z(τe) of pure elasticity becomes equal to the amplitude zp for the occurrence of the plastic flow. From Eq. (24a), it has

where τe is the time when EP transition takes place. Since ż (τe) = 0, the derivative of Eq. (24a) at τ = τe is

and evaluating the first integration of Eq. (25a) at τ = τe , it has

In order to obtain the transition time τe, integrating Eq. (25a) twice and evaluating at τ = τe with Eq. (35), we get

Therefore, combining Eq. (33) with zme = zp, it has

and using the definition of Χ again in Eq. (22), the expression of EP division is found

Growth factor are the solutions of Eq. (21) and the solving processes are similar to previous work21. Here, we list the dimensionless forms of the solutions. Two stable ones are

for the purely elastic case below zp and

with EP transition. Two unstable solutions are, respectively,


The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.

Lord, K. Hydrokinetic solutions and observations. Philos. Mag. 42, 362–377 (1871).

Article Google Scholar

von Helmholtz, H. On discontinuous movements of fluid. Philos. Mag. 36, 337 (1868).

Article MATH Google Scholar

Mousavi, A. A. & Al-Hassani, S. T. S. Numerical and experimental studies of the mechanism of the wavy interface formations in explosive/impact welding. J. Mech. Phys. Solids 53(11), 2501–2528 (2005).

Article ADS Google Scholar

Ben-Artzy, A., Stern, A., Frage, N., Shribman, V. & Sadot, O. Wave formation mechanism in magnetic pulse welding. Int. J. Impact Eng. 37(4), 397–404 (2010).

Article Google Scholar

Nassiri, A., Kinsey, B. & Chini, G. Shear instability of plastically-deforming metals in high-velocity impact welding. J. Mech. Phys. Solids 95, 351–373 (2016).

Article ADS MathSciNet CAS Google Scholar

Wang, L. F., Xue, C., Ye, W. H. & Li, Y. J. Destabilizing effect of density gradient on the Kelvin–Helmholtz instability. Phys. Plasmas 16, 112104 (2009).

Article ADS Google Scholar

Flippo, K. A. et al. Late-time mixing and turbulent behavior in high-energy-density shear experiments at high Atwood numbers. Phys. Plasmas 25, 056315 (2018).

Article ADS Google Scholar

Ravid, A., Citron, R. I. & Jeanloz, R. Hydrodynamic instability at impact interfaces and planetary implications. Nat. Commun. 12(1), 1–5 (2021).

Article Google Scholar

Raman, K. S. et al. Three-dimensional modeling and analysis of a high energy density Kelvin-Helmholtz experiment. Phys. Plasmas 19, 092112 (2012).

Article ADS Google Scholar

Zhou, Y. Rayleigh–Taylor and Richtmyer–Meshkov instability induced flow, turbulence and mixing I and II. Phys. Rep. 723–725, 1–160 (2017).

ADS MathSciNet MATH Google Scholar

Funada, T. & Joseph, D. D. Viscous potential flow analysis of Kelvin–Helmholtz instability in a channel. J. Fluid Mech. 445, 263–283 (2001).

Article ADS MathSciNet MATH Google Scholar

Kim, H., Padrino, J. C. & Joseph, D. D. Viscous effects on Kelvin–Helmholtz instability in a channel. J. Fluid Mech. 680, 398–416 (2011).

Article ADS MathSciNet MATH Google Scholar

Mikhailov, A. L. Hydrodynamic instabilities in solid media-from the object of investigation to the investigation tool. Phys. Mesomech. 10(5–6), 265–274 (2007).

Article Google Scholar

Capelli, D. et al. Development of indirectly driven shock tube targets for counter-propagating shear-driven Kelvin–Helmholtz experiments on the national ignition facility. Fusion Sci. Technol. 70, 316–323 (2016).

Article ADS Google Scholar

Nassiri, A., Chini, G., Vivek, A., Daehn, G. & Kinsey, B. Arbitrary Lagrangian–Eulerian finite element simulation and experimental investigation of wavy interfacial morphology during high velocity impact welding. Mater. Des. 88, 245–358 (2015).

Article Google Scholar

Chen, S. H. et al. Interfacial characteristics of Ti/Al joint by vaporizing foil actuator welding. J. Mater. Process. Tech. 263, 73–81 (2019).

Article CAS Google Scholar

Drennov, O. B., Davydov, A. I., Mikhailov, A. L. & Raevskii, V. A. Shear instability at the ‘"explosion product–metal"’ interface for sliding detonation of an explosive charge. Int. J. Impact Eng. 32, 155–160 (2005).

Article Google Scholar

Drennov, O. B., Mikhailov, A. L., Nizovtsev, P. N. & Raevskii, V. A. Perturbation evolution at a metal-metal interface subjected to an oblique shock wave: supersonic velocity of the point of contact. Tech. Phys. 48(8), 1001–1008 (2003).

Article CAS Google Scholar

Raoelison, R. N., Sapanathan, T., Padayodi, E., Buiron, N. & Rachik, M. Interfacial kinematics and governing mechanisms under the influence of high strain rate impact conditions: Numerical computations of experimental observations. J. Mech. Phys. Solids 96, 147–161 (2016).

Article ADS Google Scholar

Gleason, G., Sunny, S., Sadeh, S., Yu, H. L. & Malik, A. Eulerian modeling of plasma-pressure driven laser impact weld processes. Procedia Manuf. 48, 204–214 (2020).

Article Google Scholar

Wang, X., Hu, X. M., Wang, S. T., Pan, H. & Yin, J. W. Linear analysis of Atwood number effects on shear instability in the elastic–plastic solids. Sci. Rep. 11, 18049 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Gou, J. N., Sun, Y. B. & Zeng, R. H. An analytical model for the Kelvin–Helmholtz instability in solids. Int. J. Mech. Sci. 222, 107252 (2022).

Article Google Scholar

Yin, J.-W. Study on the Growth Regularity of Richtmyer-Meshkov Flow in Solid Medias with Strength (Beijing Institute of Technology, 2018).

Google Scholar

Drennov, O. B., Mikhailov, A. L., Nizovtsev, P. N. & Raevskii, V. A. Instability of an interface between steel layers acted upon by an oblique shock wave. Int. J. Impact Eng. 32, 161–172 (2005).

Article Google Scholar

Drennov, O. B. Effect of an oblique shock wave on the interface between metals. J. Appl. Mech. Tech. Phys. 56(3), 377–380 (2015).

Article ADS Google Scholar

Manikandan, P., Hokamoto, K., Deribas, A. A., Raghukandan, K. & Tomoshige, R. Explosive welding of titanium/stainless steel by controlling energetic conditions. Mater. Trans. 47(8), 2049–2055 (2006).

Article CAS Google Scholar

Vivek, A., Hansen, S. R., Liu, B. C. & Daehn, G. S. Vaporizing foil actuator: a tool for collision welding. J. Mater. Process. Technol. 213(12), 2304–2311 (2013).

Article CAS Google Scholar

Mendes, R., Ribeiro, J. B. & Loureiro, A. Effect of explosive characteristics on the explosive welding of stainless steel to carbon steel in cylindrical configuration. Mater. Des. 51, 182–192 (2013).

Article CAS Google Scholar

Vaidyanathan, P. V. & Ramanathan, A. R. Design for quality explosive welding. J. Mater. Process. Technol. 32(1), 439–448 (1992).

Article Google Scholar

Balakrishna, H.K., Venkatesh, V.C.& Philip, P.K. Influence of Collision Parameters on the Morphology of Interface in Aluminum-Steel Explosion Welds. In Shock Waves and High-strain-rate Phenomena in Metals, Springer US, 975–988 (Plenum Press, New York 1981).

Steinberg, D. J., Cochran, S. G. & Guinan, M. W. A constitutive model for metals applicable at high-strain rate. J. Appl. Phys. 51, 1498–1504 (1980).

Article ADS CAS Google Scholar

McClintock, F. A. & Argon, A. S. Mechanical Behavior of Materials (Addison-Wesley, 1966).

Google Scholar

Brown, J. L. et al. Experimental evaluation of shear modulus scaling of dynamic strength at extreme pressures. J. Appl. Phys. 128, 045901 (2020).

Article ADS CAS Google Scholar

Brown, J. L., Davis, J. P. & Seagle, C. T. Multi megabar dynamic strength measurements of Ta, Au, Pt, and Ir. J. Dyn. Behav. Mat. 7, 196–206 (2021).

Article Google Scholar

López Ortega, A., Hill, D. J., Pullin, D. I. & Meiron, D. I. Linearized Richtmyer–Meshkov flow analysis for impulsively accelerated incompressible solids. Phys. Rev. E 81, 066305 (2010).

Article ADS Google Scholar

Liu, M. B., Liu, G. R., Lam, K. Y. & Zong, Z. Smoothed particle hydrodynamics for numerical simulation of underwater explosion. Comput. Mech. 30, 106–118 (2003).

Article CAS MATH Google Scholar

Landau, L. D. & Lifshits, E. M. Theory of Elasticity 3rd edn. (Pergamon, 1986).

Google Scholar

Download references

The research was supported by the National Natural Science Foundation of China (No. 11902039).

Institute of Applied Physics and Computational Mathematics, 100094, Beijing, People's Republic of China

Xi Wang, Xiao-Mian Hu, Sheng-Tao Wang, Hao Pan & Jian-Wei Yin

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

X. W. and H. P. wrote the manuscript and worked on the methodology; X. W., X. M. H. and H. P. supervised and conceived the idea; S. T. W. in revised manuscript and arranged the funds; X. W. and J. W. Y. did the software work and simulations.

Correspondence to Hao Pan or Jian-Wei Yin.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit

Reprints and Permissions

Wang, X., Hu, XM., Wang, ST. et al. Hydrodynamic Kelvin–Helmholtz instability on metallic surface. Sci Rep 13, 2686 (2023).

Download citation

Received: 20 June 2022

Accepted: 10 February 2023

Published: 15 February 2023


Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

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.