Research Article | Open Access

Hongwei Yang, Baoshu Yin, Yunlong Shi, Qingbiao Wang, "Forced ILW-Burgers Equation as a Model for Rossby Solitary Waves Generated by Topography in Finite Depth Fluids", *Journal of Applied Mathematics*, vol. 2012, Article ID 491343, 17 pages, 2012. https://doi.org/10.1155/2012/491343

# Forced ILW-Burgers Equation as a Model for Rossby Solitary Waves Generated by Topography in Finite Depth Fluids

**Academic Editor:**Turgut Öziş

#### Abstract

The paper presents an investigation of the generation, evolution of Rossby solitary waves generated by topography in finite depth fluids. The forced ILW- (Intermediate Long Waves-) Burgers equation as a model governing the amplitude of solitary waves is first derived and shown to reduce to the KdV- (Korteweg-de Vries-) Burgers equation in shallow fluids and BO- (Benjamin-Ono-) Burgers equation in deep fluids. By analysis and calculation, the perturbation solution and some conservation relations of the ILW-Burgers equation are obtained. Finally, with the help of pseudospectral method, the numerical solutions of the forced ILW-Burgers equation are given. The results demonstrate that the detuning parameter holds important implications for the generation of the solitary waves. By comparing with the solitary waves governed by ILW-Burgers equation and BO-Burgers equation, we can conclude that the solitary waves generated by topography in finite depth fluids are different from that in deep fluids.

#### 1. Introduction

Solitary waves are finite-amplitude waves of permanent form which own their existence to a balance between nonlinear wave-steepening processes and linear wave dispersion. Typically, they consist of a single isolated wave, whose speed is an increasing function of the amplitude. The theoretical description of nonlinear solitary waves that have emerged in the mathematics and physics has received a great deal of attention in the last decades [1]. In atmospheric and oceanic circulation dynamics, Long [2] and Benney [3] first studied the barotropic Rossby waves employing horizontal shear velocity and obtained the conclusion that the amplitude of nonlinear solitary waves satisfied the KdV equation. The research on baroclinic Rossby waves was carried out by Redekopp in [4], the modified KdV (MKdV) equation as a model governing the amplitude of solitary waves was derived. Furthermore, Yamagata, Grimshaw, Boyd, and so on also discussed the solitary waves by virtue of the KdV equation [5–7]. Yano and Tsujimura classified the KdV-type Rossby solitary waves which were governed by the KdV and MKdV equations, this kind of solitary waves were also called classical solitary waves [8]. With the development of solitary waves theory, people realized that in addition to the classical solitary waves, there were other types of solitary waves. Hiroaki Ono considered the solitary waves in stratified fluids and found that in the case of great of depth, the behavior of long nonlinear waves was governed by an integrodifferential equation of dispersive type (BO equation) instead of the KdV equation [9]. This kind of solitary waves was called algebraic solitary waves. After the KdV theory and BO theory, a more general evolution equation for solitary waves in a finite-depth fluid, which reduces to the KdV equation in the shallow-fluid limit and to the BO equation in the infinitely-deep-fluid limit, was given by Kubota et al. [10]. The equation was called intermediate long waves (ILW) equation. Based on the above researches, many people carried on the research of solitary waves and explained some phenomena that occur in ocean and atmosphere, such as blocking and dipole [11–15].

As everyone knows ocean and atmosphere are driven by external forcing. The motion of ocean and atmosphere must be taken as a forced nonlinear system. Here the forcing factors include topographic forcing and external source forcing and so on. There have been lots of researches on the effect of topographic forcing on Rossby waves [12–18]. On the other hand, the real oceanic and atmospheric motion is dissipative, otherwise the motion would grow explosively because of the constant injecting of the external forcing energy. In particular, when the long-time evolution is considered, the small dissipation effects can become of crucial importance. They can stop self-similar expansion of the shock so that it tends to some stationary wave structure which propagates as a whole with constant velocity. But in many researches dissipation effect is ignored.

In the present paper, we will study the Rossby solitary waves generated by topography in finite depth fluids, especially we will consider the dissipation effect. Here we need emphasis that it seems very few ones on Rossby solitary waves excited by topography with a small dissipation in finite deep fluids by employing the forced ILW-Burgers equation are available up to now. This paper is organized as follows: in Section 2, a forced ILW-Burgers equation will be derived by using a perturbation method from the geostrophic potential vorticity equation with dissipation and topography effect. It will reduce to the KdV-Burgers equation in shallow fluids and to the BO-Burgers equation when the depth . So it is easy to find that the ILW-Burgers equation is an extension of the former researches. This is followed in Section 3 by the generation of the perturbation solution of the ILW-Burgers equation and the discussion of dissipation effect. Section 4 is devoted to a study of the conservation relations associated with the equation and the conservation quantities of Rossby solitary waves. The numerical solutions of the forced ILW-Burgers equation are given for a topographic forcing by using the pseudospectral method in Section 5. The solitary waves generated by topography are simulated, the effect of detuning parameter is analyzed. Especially, we carry out the comparison of Rossby solitary waves in finite deep fluids with that in infinite deep fluids. Finally, some conclusions are given in Section 6.

#### 2. Mathematics Model

The theoretical basis is found in the paper of Pedlosky [19] treating a model on a beta plane. In this model the vorticity equation governing the inviscid, quasi-geostrophic fluid motion with topography and turbulent dissipation, in the nondimensional form, is given by where is the dimensionless stream function; denotes the two-dimensional Laplace operator; , is the northward gradient of the planetary vorticity, and are the characteristic horizontal length and velocity scales; expresses the vorticity dissipation which is caused by the Ekman boundary layer; is a dissipative coefficient and ; is the external source, which is taken to be a function of .

To restrict our consideration to a near-resonant system, we assume where is a small disturbance in the basic flow and reflects the proximity of the system to a resonate state; is a constant, which is regarded as a Rossby waves phase speed; denotes disturbance stream-function. The substitution of (2.2) into (2.1) yields

Here, we divide the region into two parts: the domain and the domain . In the domain , in order to consider the role of nonlinearity, we must assume that shear of zonal flow exists, that is ; while in the domain , the parameter is smaller than that in the domain , we can neglect constant in the domain , meanwhile the zonal flow is assumed to be uniform, that is , here is constant. Furthermore, in the domain the topography and turbulent dissipation is absent and only consider the features of disturbances generated. For simplicity, is assumed to be smooth across .

In the domain , in order to achieve a balance among topography effect, turbulent dissipation and nonlinearity, we take Substituting (2.4) into (2.3), yields In the domain , the governing equations is The boundary conditions are where denotes the disturbance streamfunction in the domain , is that in the domain .

Following, let us consider the asymptotic expansions of (2.5) and (2.6). We introduce the following stretching transformation and the perturbation expansion in (2.5): Substituting (2.8) into (2.5) leads to the following perturbation equations: For the linear solution to be separable, we take the solution of (2.9) in the form then from (2.9) we can obtain The Equation (2.12) is an eigenvalue problem and describes the space structure of the wave along direction. The boundary conditions can be obtained from (2.6). Equation (2.12) will be solved in the latter section.

For (2.10), we have Multiplying the both sides of (2.13) by and integrating it with respect to from 0 to , utilizing the boundary conditions (2.7), we get The left hand of (2.14) can be determined by employing (2.6).

In the domain , we adopt the transformations in the following forms: and the disturbance streamfunction may be expressed by Introducing (2.15) and (2.16) into (2.6), we can get the lowest-order equation for the domain as follows: When taking the integration constant zero, we have According to Zhou [20], by virtue of the boundary condition , we can obtain the solution of (2.18) as following: where denotes the Fourier transformation of . Differentiating (2.19) with respect to , we can obtain Assuming that the solution matches smoothly at , then from (2.19) and (2.20), we get From (2.21), it is easy to find that Based on (2.20) and (2.23), we get where , and we take big enough to satisfy . So it is easy to obtain . Then we obtain from (2.22) and (2.24) Substituting (2.23) and (2.25) into (2.14) leads to where , , , . Equation (2.26) is a forced integrodifferential equation including dissipation and topography terms. In the absence of the topographic forcing and dissipation effect, (2.26) becomes the ILW equation. Here the term denotes the dissipation effect and has the same physical meaning with the term in Burgers equation, so we call (2.26) forced ILW-Burgers equation. As we know that the so-called ILW-Burgers equation is first obtained here. In the absence of forcing and dissipation, (2.26) becomes the normal ILW equation The solitary wave solution of (2.27) is where , , is the propagation speed of the solitary waves, is an amplitude parameter.

Here it should be noted that in some atmospheric and oceanic applications, the depth , (2.26) reduces to the BO-Burgers equation; in the opposite limit, (2.26) reduces to the KdV-Burgers equation. So we can conclude that the conditions describe by the KdV-Burgers equation and the BO-Burgers equation are special cases in the paper.

#### 3. Perturbation Solution of ILW-Burgers Equation

In this section, in order to study the evolutional characters of Rossby solitary waves under the influence of dissipation, we need to seek for the solution of ILW-Burgers equation and the topography effect will be studied in the latter section. In the absence of the topographic forcing, (2.26) is reduced to the ILW-Burgers equation where Next we study the perturbation solution of (3.1). Assume , , let us take a new space coordinate Assuming , then we have Taking two time scales and expanding the solution as follows: we can obtain the following approximate equations: Putting into (3.7) gives Equation (3.8) is a normal ILW equation, based on (2.28), its solution is where . So we can obtain the perturbation solution of the ILW-Burgers equation (3.1) as follows:

In order to seek the expression of , we must consider the equation (3.9). Assuming from (3.9), we get where . Equation (3.13) can be solved only when it satisfies Here satisfies the equation It is easy to find that in the case of , the solution of (3.15) is . By employing the expression of and (3.14), we obtain where is initial amplitude. So the perturbation solution of ILW-Burgers equation is where From Figures 1 and 2, it is obvious to find that the effect of a small dissipation is to cause the amplitude and moving speed of the solitary waves to decrease slowly with time.

In this section, we get the perturbation solution of ILW-Burgers equation and analyze the effect of dissipation on the amplitude and moving speed of solitary waves. While there are many methods to solve the nonlinear partial differential equation, both integrable and nonintegrable equation. In [21], a kind of new solutions of KdV equation called complexitons is presented; in [22, 23], the multiple expfunction method is successfully applied to many nonlinear partial differential equation, the rational function combinations of exponential functions will present good approximates to exact solutions; in [24, 25], a linear superposition principle of exponential traveling waves is analyzed for Hirota bilinear equations, with an aim to construct a specific subclass of N-soliton solutions formed by linear combinations of exponential traveling waves. We wonder if the above-mentioned methods can be applied to obtain the exact solution of ILW-Burgers equation. We will study these problems in the future.

#### 4. Conservation Laws

Conservation laws are a common feature of mathematical physics, where they describe the conservation of fundamental physical quantities. The role played in the science by linear and nonlinear evolution equations, in particular, by conservation laws thereof, is hard to overestimate. It is well known that the KdV equation has an infinite number of conserved quantities, some multicomponent Burgers type equations, which possess infinitely many symmetries but not infinitely many conservation laws [26]. While four conservation quantities of BO equation are found by Ono in [9]. Next we will derive the conservation relations associated with (3.1).

Following the general procedure of Ono [9], we assume that , , vanish as and integrating (3.1) and (3.1) by multiplying with powers of , we find the following two conservation relations: where we employ the property of operator . Equation (4.1) show that and are two time-invariant quantities as the dissipation effect is neglected, that is . By analogy with the KdV equation, and are regarded as the mass and momentum of the solitary waves, respectively. So we conclude that the mass and momentum of the solitary waves are conserved without dissipation. Meanwhile, we can also get that the mass and momentum of the waves decrease exponentially with the increasing of time and the dissipative coefficient in the presence of dissipation effect. Furthermore, the rate of decline of momentum is faster than the rate of mass.

Assume and construct by adding Equation (3.1) to Equation (3.1) and integrating it, by virtue of the relation , after tedious calculation, we find that . Here expresses the energy of the solitary waves. So we can conclude that the energy of the solitary waves is conserved without dissipation.

Defining a quantity related to the phase of the solitary waves we can deduce . While here we are interested in the quantity , which expresses the velocity of the center of gravity for the ensemble of such waves according to [9]. Then by using and , we obtain , which shows that the velocity of the center of gravity is conserved without dissipation.

In fact, beside the above four conservation relations, we can also verify the invariance of the following quantity:

In this section, we obtain five conservation relations of ILW-Burgers equation and draw the conclusion that the dissipation effect causes the mass, the momentum, the energy, and the velocity of the center of gravity to vary. In fact, after the above five conservation laws are given, we can wonder whether there exist other conservation laws? Is there an infinite number of conservation laws like the KdV equation? What deserves to be studied in the future?

#### 5. Numerical Simulation and Discussion

In Section 4, we have obtained the conservation relations of the homogenous ILW-Burgers equation. In this section, we will take into account the influence of topography for the Rossby solitary waves. Rossby solitary waves excited by topography with dissipation will be discussed. In order to resolve the above problems, we need to solve the forced ILW-Burgers equation. There is no analytical solution for this equation, here we will present numerical solutions of (2.26) by using the pseudospectral method [27]. First, let us simply introduce the method.

The pseudospectral method uses a Fourier transform treatment of the space dependence together with a leap-forg scheme in time. For ease of presentation the spatial period is normalized to . This interval is divided into points, then . The function can be transformed to the Fourier space by The inversion formula is These transformations can use fast Fourier Transform algorithm to efficiently perform. With this scheme, can be evaluated as , as , as and so on. Combined with a leap-frog time step, (2.26) would be approximated by The computational cost for (5.3) is six fast Fourier transforms per time step.

Before solving (2.26) by using the pseudospectral method, we need to obtain the coefficients , , and .

Taking , combining (2.12) and the boundary conditions (2.7), we can obtain the following eigenvalue problem: Assume the weak shear flow to be , where is a constant, shows that the flow shear is weak. We may get approximate expressions for and in powers of as follows: Substituting (5.5) into (5.4), we obtain The solutions of in (5.6) Substituting them into (5.6), we have Taking approximately , . Employing (5.7) and (5.8), we can get the approximate solutions for and .

As an example of calculation, we take and the topography profile is taken to be the Gaussian distribution as follows: where , , , .

The time evolutions of the solitary waves generated by the topography in the absence of dissipation are shown in Figure 3 for different detuning parameter . From Figure 3, it is easy to find that the detuning parameter holds important implications for the evolution feature of the solitary waves. In the case of , the topographic forcing generates a steady solitary wave in the forcing region. The amplitude of solitary waves increases slightly with time at the beginning and then remains invariable. A modulated cnoidal wave-train with small amplitude is in the downstream region and there is no wave in the upstream region. Between the solitary wave and the modulated wave-train, there exists a buffer region. With the decreasing of , the modulated cnoidal wave-train approaches to the topographic forcing region. The buffer region between the solitary wave and modulated cnoidal wave-train disappears slowly.

**(a)**

**(b)**

**(c)**

When , the system gets to a resonant state (Figure 3(b)). In this case, a nonsteady solitary wave is formed in the forcing region whose amplitude is larger than that in Figure 3(a). There is also a modulated cnoidal wave-train in the downstream region and there is no wave in the upstream region. The amplitude of modulated cnoidal wave-train becomes larger and the wavelength becomes shorter.

As decreases further and becomes to be negative value (Figure 3(c)), there generates a complex nonsteady wave near the forcing region. The amplitude of wave in the forcing region is larger than that in Figure 3(a) and smaller than that in Figure 3(b). The wavelength of the modulated cnoidal wave-train which exists in the downstream region is the smallest among the Figure 3. Similar with Figure 3(b), there is no buffer region between the wave in the forcing and modulated cnoidal wave-train. Unlike the Figures 3(a) and 3(b), there is also modulated cnoidal wave-train in the upstream region.

In fact, by comparing the Figure 3 in this paper with the Figure 1 in [16], which the amplitude of solitary waves is governed by the BO-Burgers equation, we can find that the solitary waves governed by ILW-Burgers equation and BO-Burgers equation exist obvious difference. For the case in [16], there is still a standing solitary wave in the forcing region, but it becomes a negative one, when , a negative standing solitary wave is formed in the forcing region, which is a stationary solitary wave disturbance. Based on the above discussion, we can conclude that the solitary waves generated by topography in finite depth fluids are different from that in infinite depth fluids.

Next we will discuss the solitary waves generated by topography in the presence of dissipation. Here we only consider the case that the detuning parameter is positive, the other cases are omitted. Comparing Figure 4 with Figure 3(a), we can find that a solitary wave is also generated in the forcing region, but it is nonsteady, its amplitude decreases with time because of the dissipation effect. Meanwhile, the dissipation causes the modulated cnoidal wave-train in the downstream region to be dissipated.

#### 6. Summary

In the paper, a new governing equation (forced ILW-Burgers) is derived to describe the amplitude of Rossby solitary waves generated by topography under the influence of dissipation. Neglecting the topography effect, the perturbation solution of ILW-Burgers is given. From the perturbation solution, we can find that the effect of a small dissipation is to cause the amplitude and moving speed of the solitary waves to decrease slowly with time. Following, five conserved quantities are obtained from the ILW equation. Based on these conserved quantities, we can conclude that the mass, momentum, energy, velocity of the center of gravity of Rossby solitary waves are conserved. Finally, the forced ILW-Burgers equation is solved numerically by using the pseudospectral method. When the dissipation is absent, the numerical results show that the detuning parameter have an important effect on the solitary waves generated by topography. The smaller the parameter is, the amplitude of solitary wave in the forcing region is larger and the solitary wave is less steady. In addition, comparing with the solitary waves governed by ILW-Burgers equation and BO-Burgers equation, we can conclude that the solitary waves generated by topography in deep fluids are different from that in infinite depth fluids. By discussing the solitary waves generated by topography in the presence of dissipation, we can find that the dissipation causes the amplitude of the solitary wave in the forcing region to decrease and the modulated cnoidal wave-train to be dissipated.

#### Acknowledgments

This work was supported by National Natural Science Foundation of China (no. 40976008), Innovation Project of Chinese Academy of Sciences (no. KZCX2-EW-209), National Natural Science Foundation of China (no. 51174128), and Graduate Innovation Foundation from Shandong University of Science and Technology (no. YCA120212).

#### References

- A. C. Scott, F. Y. F. Chu, and D. W. McLaughlin, “The soliton: a new concept in applied science,”
*Procedings of the IEEE*, vol. 61, pp. 1443–1483, 1973. View at: Google Scholar - R. R. Long, “Solitary waves in the westerlies,”
*Journal of the Atmospheric Sciences*, vol. 21, pp. 197–200, 1964. View at: Google Scholar - D. J. Benney, “Long non-linear waves in fluid flows,”
*Journal of Mathematical Physics*, vol. 45, pp. 52–63, 1966. View at: Google Scholar | Zentralblatt MATH - L. G. Redekopp, “On the theory of solitary Rossby waves,”
*Journal of Fluid Mechanics*, vol. 82, no. 4, pp. 725–745, 1977. View at: Publisher Site | Google Scholar | Zentralblatt MATH - R. Grimshaw, “Nonlinear aspects of long shelf waves,”
*Geophysical and Astrophysical Fluid Dynamics*, vol. 8, no. 1, pp. 3–16, 1977. View at: Google Scholar - J. P. Boyd, “Equatorial solitary waves. Part 1: rossby solitons,”
*Journal of Physical Oceanography*, vol. 10, no. 11, pp. 1699–1717, 1980. View at: Google Scholar - T. Yamagata, “On nonlinear planetary waves: a class of solutions missed by the traditional quasi-geostrophic approximation,”
*Journal of the Oceanographical Society of Japan*, vol. 38, no. 4, pp. 236–244, 1982. View at: Publisher Site | Google Scholar - J. I. Yano and Y. N. Tsujimura, “The domain of validity of the KdV-type solitary rossby waves in the shallow water
*β*-plane model,”*Dynamics of Atmospheres and Oceans*, vol. 11, no. 2, pp. 101–129, 1987. View at: Google Scholar - H. Ono, “Algebraic solitary waves in stratified fluids,”
*Journal of the Oceanographical Society of Japan*, vol. 39, no. 4, pp. 1082–1091, 1975. View at: Google Scholar - T. Kubota, D. R. S. Ko, and L. D. Dobbs, “Weakly-nonlinear, long internal gravity waves in stratified fluids of finite depth,”
*Journal of Hydrology*, vol. 12, no. 4, pp. 157–165, 1978. View at: Google Scholar - D. H. Luo, “Algebraic solitary Rossby wave in the atmosphere,”
*Acta Meteorologica Sinica*, vol. 49, pp. 269–277, 1991. View at: Google Scholar - G. Gottwald and R. Grimshaw, “The effect of topography on the dynamics of interacting solitary waves in the context of atmospheric blocking,”
*Journal of the Atmospheric Sciences*, vol. 56, no. 21, pp. 3663–3678, 1999. View at: Google Scholar - D. Luo, “A barotropic envelope Rossby soliton model for block-eddy interaction. I. Effect of topography,”
*Journal of the Atmospheric Sciences*, vol. 62, no. 1, pp. 5–21, 2005. View at: Publisher Site | Google Scholar - O. E. Polukhina and A. A. Kurkin, “Improved theory of nonlinear topographic Rossby waves,”
*Oceanology*, vol. 45, no. 5, pp. 607–616, 2005. View at: Google Scholar - Z. X. Luo, N. E. Davidson, F. Ping, and W. C. Zhou, “Multiple-scale interactions affecting tropical cyclone track changes,”
*Advances in Mechanical Engineering*, vol. 2011, Article ID 782590, 9 pages, 2011. View at: Publisher Site | Google Scholar - L. Meng and K. L. Lv, “Dissipation and algebraic solitary long-waves excited by localized topography,”
*Chinese Journal of Computational Physics*, vol. 19, pp. 159–167, 2002. View at: Google Scholar - L. Yang, C. Da, J. Song, H. Zhang, H. Yang, and Y. Hou, “Rossby waves with linear topography in barotropic fluids,”
*Chinese Journal of Oceanology and Limnology*, vol. 26, no. 3, pp. 334–338, 2008. View at: Publisher Site | Google Scholar - H. W. Yang, B. S. Yin, D. Z. Yang, and Z. H. Xu, “Forced solitary Rossby waves under the influence of slowly varying topography with time,”
*Chinese Physics B*, vol. 20, no. 12, Article ID 120203, 2011. View at: Publisher Site | Google Scholar - J. Pedlosky,
*Geophysical Fluid Dynamics*, Springer, New York, NY, USA, 1979. - Q. F. Zhou, “The second-order solitary waves in stratified fluids with finite depth,”
*Scientia Sinica Series A*, vol. 28, no. 2, pp. 159–171, 1985. View at: Google Scholar | Zentralblatt MATH - W.-X. Ma, “Complexiton solutions of the Korteweg-de Vries equation with self-consistent sources,”
*Chaos, Solitons & Fractals*, vol. 26, no. 5, pp. 1453–1458, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH - W. X. Ma, T. Huang, and Y. Zhang, “A multiple exp-function method for nonlinear differential equations and its application,”
*Physica Scripta*, vol. 82, no. 6, Article ID 065003, 2010. View at: Publisher Site | Google Scholar - W.-X. Ma and Z. Zhu, “Solving the (3 + 1)-dimensional generalized KP and BKP equations by the multiple exp-function algorithm,”
*Applied Mathematics and Computation*, vol. 218, no. 24, pp. 11871–11879, 2012. View at: Publisher Site | Google Scholar - W.-X. Ma and E. Fan, “Linear superposition principle applying to Hirota bilinear equations,”
*Computers & Mathematics with Applications*, vol. 61, no. 4, pp. 950–959, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH - W.-X. Ma, Y. Zhang, Y. Tang, and J. Tu, “Hirota bilinear equations with linear subspaces of solutions,”
*Applied Mathematics and Computation*, vol. 218, no. 13, pp. 7174–7183, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH - W.-X. Ma and Z.-X. Zhou, “Coupled integrable systems associated with a polynomial spectral problem and their Virasoro symmetry algebras,”
*Progress of Theoretical Physics*, vol. 96, no. 2, pp. 449–457, 1996. View at: Publisher Site | Google Scholar - B. Fornberg,
*A Practical Guide to Pseudospectral Methods*, vol. 1, Cambridge University Press, Cambridge, UK, 1996. View at: Publisher Site

#### Copyright

Copyright © 2012 Hongwei Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.