Dark Solitons in Acoustic Transmission Line Metamaterials

We study dark solitons, namely density dips with a phase jump across the density minimum, in a one-dimensional, weakly lossy nonlinear acoustic metamaterial, composed of a waveguide featuring a periodic array of side holes. Relying on the electroacoustic analogy and the transmission line approach, we derive a lattice model which, in the continuum approximation, leads to a nonlinear, dispersive and dissipative wave equation. The latter, using the method of multiple scales, is reduced to a defocusing nonlinear Schrödinger equation, which leads to dark soliton solutions. The dissipative dynamics of these structures is studied via soliton perturbation theory. We investigate the role—and interplay between—nonlinearity, dispersion and dissipation on the soliton formation and dynamics. Our analytical predictions are corroborated by direct numerical simulations.


Introduction
A dark soliton is an envelope soliton that has the form of a density dip with a phase jump across its density minimum.This localized waveform is supported on the top of a stable continuous wave background.Since dark solitons are fundamental nonlinear excitations of a universal model, the defocusing nonlinear Schrödinger (NLS) equation, they have been studied extensively in diverse branches of physics.These include chiefly nonlinear optics [1] and Bose-Einstein condensates [2], but also discrete mechanical systems [3], thin magnetic films [4], complex plasmas [5], water waves [6], and so on.The interest of dark solitons arises from the fact that-since they are composite objects, consisting of a background wave and a soliton-they have a number of interesting properties that may find important applications.Indeed, as compared to the bright solitons (which are governed by the focusing NLS equation), they can be generated by a thresholdless process, are less affected by loss and background noise, and are more robust against various perturbations [1,2].In addition, in optics, dark solitons have potential applications, e.g., in inducing steerable waveguides in optical media, or for ultradense wavelength-division-multiplexing [1].
Recently, important nonlinear effects-including higher harmonic generation [7][8][9], self-demodulation [10,11], and soliton formation [12][13][14][15][16], and others-have attracted much attention in acoustics, and especially in the context of acoustic metamaterials; these are periodic structures-thus featuring dispersion-which are characterized by negative effective parameters, i.e., negative mass density [17] and negative bulk modulus [18].Usually, such metamaterials are acoustic waveguides loaded with resonators, presenting strong dispersion around the resonance frequency and Bragg band [19,20].Dissipation naturally plays a key role in such structures, especially in super-absorbing sound materials [21,22] (see also work in Refs.[23,24]).The combination of dispersion and dissipation [20,23,[25][26][27] has mainly been investigated in the linear regime; nevertheless, the additional effect of nonlinearity, which is naturally introduced at high acoustic levels has not been studied in detail in acoustic metamaterials.This is also the case of studies on acoustic solitons of various types-pulse-like [12,13] and envelope bright ones [12,14].As far as dark solitons are concerned, they were predicted to occur in acoustic waveguides loaded with an array with Helmholtz resonators [12], while the effect of dissipation was neglected.To the best of our knowledge, the dissipation-induced dynamics of dark solitons in acoustic metamaterials, has not been investigated so far.
In this work, we analytically and numerically study envelope dark solitons in a one-dimensional (1D) acoustic metamaterial composed of a waveguide featuring a periodic array of side holes.Based on the transmission line (TL) approach, which is widely used in acoustics [28][29][30][31][32], we model this system by a nonlinear dynamical lattice with losses.In the continuum approximation, the latter leads to a nonlinear, dispersive and dissipative wave equation.Employing the multiple scales perturbation method, we derive an effective defocusing nonlinear Schrödinger (NLS) equation which supports dark soliton solutions.The effect of dissipation is also taken into regard, and is studied analytically by means of perturbation theory for solitons [33].Numerical simulations, at the level of the nonlinear lattice model, are also performed showing a very good agreement with the analytical predictions relying on the NLS picture.The methodology of this work paves the way for relevant studies in double-negative acoustic metamaterials.
The paper is structured as follows.In Section 2, we introduce our setup, derive the 1D nonlinear lattice model, and study the linear properties of the system.In Section 3, we employ the multiple scale perturbation method to obtain the NLS equation and its envelope dark soliton solutions; both the lossless and dissipative cases are studied.Finally, in Section 4, we present our conclusions and discuss some future research directions.

Setup and Model
We consider plane wave propagation in a structure composed of a waveguide of radius r periodically loaded with an array of side holes of radius r H and length l H , and the distance between two consecutive cells is d, as shown in Figure 1a.We choose to work in a frequency range well below the first cut-off frequency of the higher propagating modes in the waveguide, and therefore the problem is considered as 1D.
In order to theoretically analyze the problem, we employ the electro-acoustic analogue modeling based on the transmission line (TL) approach.This allows us to derive a nonlinear, dispersive and dissipative discrete wave equation, which-in the continuum limit-can be treated analytically via the multiple scales perturbation method.The corresponding unit-cell circuit of the TL model is shown in Figure 1b, where the voltage v and the current i correspond to the acoustic pressure p and to the volume velocity u through the waveguide's cross-sectional area, respectively [12,28].
The first part of the unit-cell circuit, corresponding to the waveguide, is modeled by the inductance L ω , the resistance R ω and shunt capacitance C ω .The linear part of the inductance and the capacitance are L ω0 = ρ 0 d/S and C ω0 = Sd/(ρ 0 c 2 0 ), respectively, where S = πr 2 , ρ 0 and c 0 are the cross section of the waveguide, the density and the sound velocity of the fluid in the system, respectively.We consider weakly nonlinear wave propagation, where the celerity c NL ≈ c 0 1 + β 0 p/ρ 0 c 2 0 , with β 0 being the nonlinear parameter, depends on the amplitude.Thus, we consider the capacitance C ω to be nonlinear, depending on the pressure p, while the inductance L ω linear [12].The pressure-dependent capacitance [14] can be expressed as C ω = C ω0 − C ω p n , where C ω = 2β 0 C ω0 /(ρ 0 c 2 0 ).The resistance R ω = Im(kZ c )d, corresponds to the viscothermal losses at the boundaries of the waveguide wall: is the wavenumber and is the characteristic impedance of the waveguide, with γ being specific heat ratio, Pr the Prandtl number, s = ωρ 0 r 2 /η, and η the shear viscosity.
The second part of the unit-cell circuit describes the effects of the holes that are considered to be linear (assuming sufficiently smoothed edges of the holes) [34,35].The side holes could be modeled by a shunt LR circuit composed of the series combination of an inductance L H and a resistance R H .We study the regime where the wavelength λ of the sound wave is much bigger than the geometric characteristic of the side holes, i.e., kr H 1. The corresponding inductance could be given by L H = ρ 0 (l H + ∆l Hi + ∆l Ho )/S H , where ∆l Hi and ∆l Ho are length corrections due to the radiation inside the waveguide and to the outer environment, respectively [36][37][38] (see more details in Appendix A) and S H is the area of the side holes.The resistance R H = ρ 0 c 0 (kr H ) 2 /(2S H ) = ρ 0 ω 2 /(2πc 0 ) corresponds to the radiation losses to the outer environment and the propagation losses due to viscous and thermal effects in the side holes [36].
Here, we should mention that we approximate the frequency-dependent viscothermal losses of the waveguide and the radiation losses due to the side holes by resistance, R ω and R H , respectively, with the corresponding constant value fixed at the characteristic frequency of the dark solitons (see below).
We now use Kirchhoff's voltage and current laws to derive the discrete nonlinear dissipative evolution equation for the pressure p n in the nth cell of the lattice (see details in Appendix B):

Continuum Limit
We now focus on the continuum limit of Equation (3), corresponding to n → ∞ and d → 0 (but with nd being finite); in such a case, the pressure becomes p n (t) → p(x, t), where x = nd is a continuous variable, and i.e., , where subscripts denote partial derivatives.In order to understand more accurately the dynamics of the system, we keep terms up to order O(d 4 ) coming from the periodicity of the side holes [14].This way, we obtain the corresponding partial differential equation (PDE), It is also convenient to express our model in dimensionless form by introducing the normalized variables τ and χ and normalized pressure P, which are defined as follows: τ is time in units of , where ω B = πc 0 /d is the Bragg frequency; χ is space in units of c 0 /ω B and p = P 0 P, where P 0 = ρ 0 c 2 0 and 0 < 1 is a formal dimensionless small parameter.Then, Equation ( 5) is reduced to the following dimensionless form: where, It is interesting to identify various limiting cases of Equation ( 6).First, in the linear limit (β 0 = 0, or p 2 1), in the absence of side holes (m 2 → 0, γ H → 0) and without considering viscothermal losses (γ ω → 0) and higher order spatial derivatives, Equation ( 6) is reduced to the linear wave equation, P ττ − P χχ = 0.In the linear limit, in the presence of side holes, in the long wavelength approximation and without considering viscothermal losses (γ ω → 0), radiation losses (γ H → 0), and higher order spatial derivatives (ζ → 0), Equation ( 6) takes the form of the linear Klein-Gordon equation [39,40], P ττ − P χχ + m 2 P = 0. Finally, in the nonlinear lossless regime, and in the absence of sides holes, without considering higher order spatial derivatives, Equation ( 6) is reduced to the well-known Westervelt equation, P ττ − P χχ − β 0 P 2 ττ = 0, which is a common nonlinear model describing 1D acoustic wave propagation [10].

Linear Limit
We now consider the linear limit of Equation ( 6) in order to obtain the dispersion relation of the system.Assuming plane waves of the form P ∝ exp[i(kχ − ωτ)], we obtain the following complex dispersion relation connecting the wavenumber k and frequency ω, where γ ω and γ H account for the viscothermal losses in the waveguide and the losses (radiation losses and viscothermal losses) due to the side holes, respectively.In the absence of losses, Equation ( 8) is reduced to: which is the familiar dispersion relation of the linear Klein-Gordon model [39,40], with a higher-order spatial derivative term ζk 4 accounting for the influence of the periodicity of the system to the dispersion relation.Although this term appears to lead to instabilities for large values of k, in the long-wavelength limit where k is sufficiently small, such instabilities do not occur.For frequencies between 0 ≤ ω < m, there is a band gap, and, for m < ω < ω B , there is a propagating band, with the dispersion curve ω(k) having the form of hyperbola (see Figure 2).Considering that all quantities in the dispersion relation (Equations ( 8) and ( 9)) are dimensionless, it is also relevant to express them in physical units.The frequency ω ph and wavenumber k ph in physical units are connected with their dimensionless counterparts through ω = ω ph /ω B and k = k ph c/ω B .Then, we rewrite Equation ( 8) with physical units as follows: The real and imaginary parts of Equation (10), shown in Figure 2, are found to be almost identical with the lossless case i.e., γ ω = 0 and γ H = 0 (dashed blue lines).Thus, the medium is considered to be weakly lossy.In particular, below we neglect viscothermal losses in the waveguide (γ H = 0), which are small compared to radiation losses.Furthermore, we assume that the remaining losses are sufficiently small, such that γ H = 2 γ H .
On the other hand, the green stars in Figure 2 show the corresponding results for the lossy dispersion relation obtained by using the transfer matrix method (TMM) [20]: where k and Z c are given by Equations ( 1) and ( 2), respectively, and Z H = iω ph L H + R H is the input impedance of the side holes.For the lossless case, k, Z c and Z H are respectively reduced to k = ω/c 0 , Z c = ρ 0 c 0 /S and Z H = iω ph L H .The corresponding results are shown as thick green continuous lines in Figure 2. We find that the linear dispersion relation resulting from the TLM is in very good agreement with the one obtained by using the TMM, as shown in Figure 2 in the low frequencies regime.
Here, we mention that results shown in Figure 2 (as well as in Figures 3 and 4

Dark Solitons
In this section, we use multiple scales perturbation method to describe the evolution of the pressure using an effective lossy NLS equation.From the latter equation, we can obtain the approximate analytical dark soliton solutions of our lattice system and study numerically their dynamics.
We start our analysis by introducing the slow variables, and express P as an asymptotic series in , Then, substituting Equations ( 12) and ( 13) into Equation ( 6), we obtain a hierarchy of equations at various orders in (see more details in Appendix C).
The leading order O( 0 ), Equation (A8) in Appendix C, possesses a linear plane wave solution of the form where A is an unknown envelop function, θ = k r χ 0 − ωτ 0 with the wave number k r and the frequency ω satisfying the lossless dispersion relation, Equation ( 9), and c.c. denotes a complex conjugate.
Next, at the order O( 1 ), the solvable condition is that the secular part (i.e., the term ∝ exp(iθ)) vanishes, which yields, where is the inverse group velocity.Equation ( 15) is satisfied as long as A depends on the variables χ 1 and τ 1 through the traveling-wave coordinate τ1 At the same order, we could obtain the form of the field p 1 , where B is an unknown function that can be found at a higher-order approximation.Finally, following the same process as above, the nonsecularity condition at O( 2 ) yields a lossy NLS equation for the envelop function A, where the dispersion, nonlinearity and dissipation coefficients are respectively given by, The sign of the product σ ≡ sgn(qk ) determines the nature of the NLS equation [39,40].A focusing NLS equation with σ = +1 supports bright solitons, localized waves with vanishing tails towards infinity.A defocusing NLS equation with σ = −1 supports dark solitons, which are density dips, with a phase jump across the density minimum, on top of a non-vaninishing continuous wave background.Figure 2c shows that the product qk is always negative, i.e., σ = −1 in our case, which means that our system only supports dark solitons.
In the absence of losses (Λ = 0), the analytical dark soliton solution for the envelope function A is of the form, where η 0 is a free parameter setting the amplitude of the dark soliton background, and parameters a 0 = sin φ, b 0 = cos φ are connected by a relation, a 2 0 + b 2 0 = 1, with the effective angle 2φ corresponding to the phase shift across the dark soliton.The parameter b 0 = cos φ characterizes the soliton intensity at the center.A black soliton is a special case of a dark soliton, with π phase shift and zero intensity at the center.When φ = 0, the minimum intensity of a dark soliton does not equal zero, and we call it a gray soliton.

Black Solitons
We start with the black soliton, with φ = 0, which at the NLS level is expressed as: The corresponding approximate black soliton solution of Equation ( 6) is as follows: which is a function of parameters χ and τ.In the original coordinates, space x and time t, the approximate black soliton solution for the pressure p reads where k r /c 0 is the inverse group velocity at the carrier frequency that is independent of the amplitude of the background 2 η 0 P 0 .
To ensure the balance between nonlinearity and dispersion, we should also pay attention to the dispersion length L D and the nonlinearity length L NL providing the length scales over which dispersive or nonlinear effects become important for pulse evolution.These lengths are expressed as follows [14]: and where t 0 and A 0 are the characteristic width and the amplitude of the propagating envelope, respectively.From the black soliton solution, Equation (25), we find t 0 = η 0 q/k r ω B −1 and A 0 = η 0 P 0 .Substituting them in Equations ( 26) and ( 27), we obtain L NL /L D ∼ 1.Now, we study numerically the evolution of the approximate black soliton solution of Equation ( 25) in the lattice.We start by integrating the lossless version of Equation ( 3) with a driver of the form given by Equation ( 25) at x = 0 with = 0.029 (2 ηP 0 = 8816 Pa) (170 dB) and carrier frequency f = 600 Hz.The numerical results are shown in Figure 3a,b.We observe that the black soliton with zero intensity at the center propagates with a constant velocity, amplitude and width as shown in the contour plot (Figure 3a).Good agreement between the simulations (blue line) calculated at t = 1.41 s and the corresponding analytical solution of Equation (25) (green dashed line) is shown in Figure 3b.In order to confirm that the black solitons are due to the balance between dispersion and nonlinearity, we also compare Figure 3a,b to the unbalanced case where the boundary condition has the same width but a smaller amplitude, = 0.0095 (2 ηP 0 = 2888 Pa) (160 dB), shown in Figure 3c,d.In Figure 3d, it is observed that the soliton is not formed and the initial wavepacket spreads.

Black Solitons under Dissipation
We next study the dynamics of the black solitons under the presence of the radiation and viscothermal losses.The effective NLS (18), which includes a linear loss term, can be studied via the direct perturbation theory for dark solitons [33].The loss term does not vanish for τ1 → ∞ and, thus, in this limit (where ∂ 2 A/∂ τ2 1 → 0), the evolution of the black soliton's background A ≡ η(χ 2 ) is determined by the equation: Solutions of Equation ( 28) can be sought in the form η(χ 2 ) = η 0 (χ 2 ) exp[iφ(χ 2 )], where η 0 (χ 2 ) describes the loss-induced change of the background amplitude, and φ(χ 2 ) is the spatically varying phase of the background.Substituting the above ansatz into Equation ( 28), and separating real and imaginary parts, we derive the equation for the amplitude of the background as follows: Integrating Equation ( 29), we find that the background amplitude decays according to the following exponential law: In terms of the original space coordinate, the amplitude of the background of the black soliton decreases exponentially as: where is the dissipation length.Thus, the envelope of the approximate black soliton solution reads During propagation, the amplitude of background η decreases (due to the presence of loss), while the minimum intensity is always zero.The black soliton does not move against the background, and the velocity of the background only depends on the carrier frequency.Thus, the linear loss only affects the black soliton's background amplitude and width.

Dissipation Length Nonlinearity Length and Dispersion Length
We numerically integrate the lattice model, Equation (3), using a driver corresponding to the black soliton shown in Figure 3 with = 0.029 (2 ηP 0 = 8816 Pa) (170 dB) and f = 600 Hz.The numerical result (blue line) shown in Figure 4 is in good agreement with the analytical result (green dashed line), Equation (31).The radiation and viscothermal losses strongly attenuate the background amplitude of the black soliton, according to the value of the dissipation length L l = 9.12 m, which is much smaller than the nonlinearity and dispersion lengths, given by L NL = L D = 21.5 m, illustrated by the light green thick line in Figure 4.

Dissipation Length ≥ Nonlinearity Length and Dispersion Length
Next, we choose the parameters of our setup so that L l ≥ L NL , L D .This is achieved by replacing air to water, and by decreasing the radius of the side holes.We now use the following parameters for the lattice: d = 0.05 m, r = 0.025 m, r H = 0.002 m, l H = 0.002 m, at 25 • C, with nonlinear parameter β 0 = 3.6, velocity c 0 = 1483 m/s, density ρ 0 = 998.3kg/m 3 , specific heat ration γ = C p /C v = 1.01, specific heat (constant pressure) C p = 75.327KJ/mol/K, specific heat (constant volume) C v = 74.53KJ/mol/K, Prandtl number Pr = 7.01 and dynamic viscosity η = 1.002 10 −3 kg/m/s.We use an envelope given by Equation ( 25) with = 0.01 (2 ηP 0 = 4.39 × 10 7 Pa) and f = Hz.In this case, the background amplitude of the dark soliton is weakly attenuated, and the numerical results (blue lines) are in good agreement with the corresponding analytical ones (green dashed lines) in Figure 5, Equation (31).Below, we adopt the parameters corresponding to the water-filled acoustic metamaterial.

Gray Solitons
Apart from the stationary black soliton (characterized by a zero density minimum), there exists the moving gray soliton solution, with φ = 0 (with a nonzero density minimum).The gray soliton solution of Equation ( 6) reads: and is a function of parameters χ and τ.In terms of the original coordinates x and t, the approximate gray soliton solution for the pressure p(x, t) is given by: Here, the velocity v g is given by: and is the sum of the background velocity, v b = c 0 /k r , controlled by the carrier frequency, and the gray soliton's velocity, v s = c 0 / a 0 η 0 qk r , which depends on a 0 = sin φ and η 0 .
We now focus on the numerical study of the evolution of the gray soliton in the lattice.We start with the lossless case, by integrating Equation ( 3) with a driver of the form given by Equation ( 35) at x = 0, with = 0.01 (2 ηP 0 = 4.39 × 10 7 Pa) and f = 1300 Hz.The phase angle in the simulations is chosen as φ = π/4, φ = π/8 and φ = −π/8-see, respectively, Figure 6a-f.We observe that the gray solitons propagate at a constant velocity with a constant amplitude, width and minimum intensity, as shown in the contour plot in Figure 6a,c,e.The numerical spatial profiles calculated at t = 0.35 s (blue lines) are in good agreement with the corresponding analytical ones, Equation ( 35) (green dashed lines), as shown in Figure 6b,d,f.The phase angle φ determines the minimum intensity and the velocity of the gray soliton − → v gray , as well as the group velocity, v g = v b + v s .Upon changing φ (from φ = π/4 to φ = π/8 ), we observe that the group velocity and the minimum intensity decrease (see Figure 6b,d).The change of sign of φ (from φ = π/8 to φ = −π/8 ) does not affect the minimum intensity of the gray soliton, as shown in Figure 6d,f.For φ = −π/8, the gray soliton moves in a direction opposite to that of the background (see Equation ( 36)).Our analytical prediction of Equation (36), illustrated by black dashed lines in Figure 6a,c,e, is in full agreement with the numerical results.

Gray Solitons under Dissipation
According to the perturbation theory for dark solitons [33], the gray soliton of Equation ( 18) evolves as: where a = a 0 exp(−Λχ 2 ), The background amplitude decays exponentially with the same rate as in the black soliton case, Equation (30).The first-order correction term is given by: and σ 0 is an extra phase induced by the perturbation.In the original coordinates, x and t, the approximate form of the gray soliton reads: with the velocity During gray soliton propagation, both the background amplitude and the minimum intensity decrease exponentially due to losses, while the width increases (see Equation ( 41)).
We numerically integrate the lattice model, Equation (3), with a gray soliton driver of the form given by Equation (35) or Equation (41), with = 0.01 (2 ηP 0 = 4.39 × 10 7 Pa) and f = 1300 Hz.The numerical profile found at t = 0.58 s, with φ = π/8, is shown in Figure 7a (blue line), is in a good agreement with the analytical prediction (green dashed line), Equation (41).The numerically obtained gray soliton trajectory [blue circles in Figure 7b], and is in good agreement with the analytical expression given in Equation (42).

Conclusions
In conclusion, we have theoretically and numerically studied envelope dark (black and gray) solitons in 1D nonlinear acoustic metamaterial composed of a waveguide with a periodic array of side holes, featuring viscothermal and radiation losses.Based on the electro-acoustic analogy and transmission line approach, we derived a nonlinear lattice model.We studied the linear dispersion relation of our system, which was found to be in good agreement with the one obtained by using the transfer matrix method.We used a multiple scale perturbation method to analytically treat the problem, and derived an effective NLS model describing the evolution of the pressure.We predicted the existence of dark solitons and studied analytically their evolution in the presence of losses.We investigated the interplay between dissipation, nonlinearity and dispersion, as described by pertinent characteristic length scales, in the case of a water-filled acoustic metamaterial with a small radius of side holes.The numerical results were found to be in good agreement with the analytical predictions.
Our results and methodology pave the way for studies on nonlinear phenomena in double-negative acoustic metamaterials, e.g., waveguides periodically loaded with side holes and clamped elastic plates.Our works will also help in designing new nonlinear acoustic metamaterials supporting bright, gap, black and grey solitons.According to the properties of the solitons, i.e., propagation with no distortion of high amplitude waves, our studies could also pave the way to designing new devices in medical applications or to designing non destructive sensors.

Figure 1 .
Figure 1.(color online) (a) waveguide loaded with an array of side holes; (b) corresponding unit-cell circuit.

Figure 2 .
Figure 2. (color online)(a,b) show the real and imaginary parts of the complex dispersion relation respectively.The blue circles (green stars) show the results from the TL (TMM) approach from Equation(10) [Equation(11)]; the blue dashed (green continuous) lines show the lossless dispersion relation obtained from the TL (TMM) approach from the lossless limit of Equation (10) (Equation (11)); (c) the frequency dependent qk , the product of the nonlinearity and the dispersion coefficients of NLS equations.

Figure 3 .
Figure 3. (color online) (a) a contour plot of black soliton of the form of Equation (25), obtained by numerically integrating the lossless version of Equation (3) (R ω = 0, R H = 0) with = 0.029 (2 ηP 0 = 8816 Pa) (170 dB) and carrier frequency f = 600 Hz; (b) numerical spatial profile of black soliton calculated at t = 1.41 s (blue line).Green dashed lines are the corresponding analytical envelope result of Equation (25).The light green thick line denotes the nonlinear length L NL and dispersion length L D , with L NL = L D = 21.55 m; (c) a numerical contour plot for dispersive effect, obtained by numerically integrating the lossless version of Equation (3) (R ω = 0, R H = 0).Keeping the same width, we only decrease the amplitude of the driver to = 0.0095 (2 ηP 0 = 2888 Pa) (160 dB) such as the nonlinearity can not balance the dispersion; (d) numerical spatial profile of dispersive wave calculated at t = 1.41 s.

Figure 4 .
Figure 4. (color online) Effects of radiation and viscothermal losses on black soliton propagation in an air-filled acoustic metamaterial with d = 0.05 m, r = 0.025 m, r H = 0.005 m, l H = 0.002, calculated at t = 1.41 s.Numerical dissipative black soliton (blue line) is obtained by integrating our lattice model Equation (3) with the driver corresponding to = 0.029 (2 ηP 0 = 8816 Pa) (170 dB) and f = 600 Hz, Equation(25) or Equation(33).The green dashed line is the corresponding analytical envelope result, Equation(33).The light green thick line denotes the nonlinear length L NL and dispersion length L D , with L NL = L D = 21.55 m.

Figure 5 .
Figure 5. (color online) Effects of radiation and viscothermal losses on black soliton propagation in an water-filled acoustic metamaterial with d = 0.05 m, r = 0.025 m, r H = 0.002 m, l H = 0.002 m, calculated at (a) t = 0.56 s; (b) t = 0.58 s and (c) t = 0.6 s, respectively.A numerical dissipative black soliton (blue line) is obtained by integrating the lattice model Equation (3) with the driver corresponding to = 0.01 (2 ηP 0 = 4.39 × 10 7 Pa) and f = 1300 Hz. Green dashed line is the corresponding analytical envelope results, Equation (33).The light green thick line denotes the nonlinear length L NL and dispersion length L D , with L NL = L D = 36.54m.

Figure 7 .
Figure 7. (color online) (a) numerical spatial profile of dissipative gray soliton calculated at t = 0.58 s with φ = π/8 (blue line).Green dashed lines are the corresponding analytical envelope, Equation (41); (b) the corresponding space-time diagram of dissipative gray soliton (blue circles).The slope of the black line shows the analytical velocity of dissipative gray solitons, Equation (42).