Travelling Waves: Interplay of Low to High Reynolds Number and Tan-Cot Function Method to Solve Burger’s Equations ()
1. Introduction
The wave propagation is one of the important pillar of both linear and nonlinear partial differential equations. A wave is prominently observable which is transported from one segment of the medium to another segment with a recognizable speed of propagation. The mathematical term of wave is a function of the form
where c is a constant known as wave speed and greater than zero; u is a wave function depends on two variables x and t. Here t represents the time, the initial term translated to the right by ct spatial units.
Many important sectors of sciences such as chemical, physical, economical and biological are introduced by Burger’s type nonlinear partial differential equations. It is also noted that in the area of nonlinear modeling of fluid flow and population dynamics the equations pattern like Burger’s is rigorously used to present various numerous scientific outcome. Exact solutions to nonlinear differential equations play an important role in physical science since they can provide much more physical information and insight which leads to further applications [1] [2] [3] [4] [5] . Nonlinear wave propagation phenomena of dispersal, dissipation, dispersion, and reaction-diffusion-convection are significant in wave equations. Some of the applications are drawn in the following series:
• Fluid mechanics (water waves, aerodynamics);
• Acoustics (sound waves in air and liquids);
• Elasticity (stress waves, earthquakes);
• Electromagnetic theory (optics, electromagnetic waves);
• Biology (epizootic waves);
• Chemistry (combustion and detonation waves).
To solve the nonlinear equations, there are many analytical and numerical methods in the literature and among them are tan-cot scheme [1] [2] , implicit finite difference scheme [3] , Hopf-Cole transformation [4] , hyperbolic conservation laws [5] [6] [7] [8] , sine-cosine method [9] [10] and tanh-sech method [11] .
To find the solutions of different patterns of Burgers’ equations, variant numerical methods are introduced. Anwar and Ali used tanh and tan-cot schemes for exact complex solutions of some different types of nonlinear partial differential equations [1] [2] whereas Soliman [3] used implicit finite difference scheme. The system of two-dimensional equations is solved by the method called Hopf-Cole transformation [4] . Fletcher used this transformation to give an analytic solution for the system of required equations [4] . In viscous fluid, the approximate theory of flow is used through a shock wave traveling. The persistence of viscous term supports to stop the wave-breaking, and smooth out shock discontinuities and finally the prediction gained a well-behaved and smooth solution. Solutions and characteristics of the Burger’s equation with non-zero viscosity, considering shock speed with weak solution are interpreted in [4] - [9] . Several hyperbolic conservation laws with various methods such as Riemann, conservative methods for nonlinear problems which are generalized for the advection equation are discussed in [11] - [16] .
In this paper, we discuss the travelling wave solutions of different types of Burger’s equations analytically using tan-cot function method. For numerical illustration, we employ finite difference method which is based on Crank-Nicolson scheme to solve one dimensional Burger’s equation. Initially, we list out the equations for further analysis. Nonlinear differential equations such as Burger’s equation have various types [1] [3] [4] [5] and two well known of them are
1) Travelling wave solution for Burger’s equation with viscosity, and
2) Travelling wave solution of the KdV-Burger’s equation.
A simple introduction of these two different forms of Burger’s equations is given below:
Consider the Burger equation of the form
(1)
where d is the viscosity of the fluid and g is a nonlinear function of u. Introducing the function
, the Burger equation can be written in the form of
(2)
here u is the fluid density, d is the viscosity, the spatial variable x and time is t. In absence of viscosity,
, and the Equation (2) becomes
(3)
which is a non-viscid Burger’s equation and the travelling wave solution does not exist. The inviscid (non-viscid) Burgers’ equation is a quasilinear wave conservation equation. The analytic traveling wave solution of this equation is an implicit relation that provided characteristics and do not intersect. If the characteristics do intersect, a classical solution to this equation does not exist and leads to the formation of a shock wave. For smooth initial data, the solution satisfied the condition
is always satisfied for sufficiently small time t, where
is the solution. Otherwise, the solution develops a singularity (discontinuity) and arise the shock wave.
Let us now consider the KdV-Burger’s equation of the form
(4)
which is also known as Korteweg-de Vries Burger’s equation and here a and k are arbitrary constants and d is the viscosity. If
and
then the Equation (4) turn to the Burger’s equation.
In the following section, we will discuss the tan-cot function method to get the exact solution of Burger’s equation.
2. Tan-Cot Function Method
For travelling wave solution of Burger’s equation, we used the tan-cot function method as introduced in [1] [2] . The key idea of this method is to discuss in the following steps. Consider the nonlinear partial differential equation in the following form [15] [16]
(5)
where
is a travelling wave solution of (5). Now using the travelling wave transformation,
(6)
where
. Here m and n are real constant and transformation gives us following changes
(7)
By integrating Formula (7) with respect to s, a non-linear partial differential Equation (5) transforms to a nonlinear ordinary differential equation such that
(8)
As an instant example, integrate the KdV Burger’s equation, (4) with respect to dummy variable s instead of x produces the following differential equation:
where
to be determined and the integral constant
.
The derivative terms of Equation (8) can be reduced by using integration while the integration constants are negligible and we will obtain the solutions of many nonlinear equations represented by the form
(9)
where
are parameters and to be determined and
is known as a wave number.
Remark 1. To get a function
from the ordinary differential Equation (8), the key principle is to integrate this formula as long as all terms contain derivatives and setting the constant of integration to be zero.
First we can write their solutions and derivatives using tan function such that
(10)
Similarly, for cot function, we have
(11)
Replacing the Equations (10) or (11) in (8) and after balancing the terms of either tan or cot functions, we can solve the resulting system of algebraic equations by using computerized program. This method is defined as tan-cot function method.
Next, we will explore the tan-cot function method to find the pattern of exact solutions.
3. Exploration of Burger’s Equations with Respect to Travelling Wave Solution
In this section, we are interested to solve the Burger’s type equations as defined in the earlier section for finding the travelling wave solutions.
3.1. Travelling Wave Solution for Burger’s Equation with Viscocity
Let us recall the equation in (2) and the wave transformation is
where wave variable
and c is the wave speed. At first, we consider the wave variable,
and after using the transformation Equation (8), which yields the following ordinary differential equation from (2)
(12)
Integrating Equation (12) once with respect to s, we have
(13)
Substituting
and
from (10), the Equation (13) implies
(14)
Equating the exponents and the coefficients of each pair of the tan functions, we find that
and we obtain the following relation
Substituting the values of
and
in (10), the solution of the Equation (2) is of the form
(15)
which is the travelling wave solution of one dimensional Burger’s equation.
Figure 1 depicted the physical structure of travelling wave which counts the wave number and distributions. If we use the transformation
, where
then Equation (2) reduces to
(16)
and on integration with respect to s, we obtain
Figure 1. Travelling waves of Burger’s equation for
while (left)
, (middle)
, and (right)
at time
over the domain
.
(17)
where A is a constant. For equilibrium points
Hence the roots are
Let us consider
We observe that
and obviously
then the travelling wave solutions are possible if
for
. Now taking integration in (17), we have
After solving this we get
where
Now we can write
as follows
(18)
After that multiplying and dividing by
, we obtain
Initial value
is not essential and the derivative of
is closer to zero because of smaller m and larger d. Moreover, we obtain inviscid Burgers equation which cannot give us continuous travelling waves for
. So the wave solution is
(19)
This solution is similar to shock wave profile because of joining the asymptotic states of G2 and G1 with the boundary conditions
3.2. Numerical Solutions for Low to High Reynolds Number
For numerical schematics, we choose finite-difference method based on Crank-Nicolson implicit time differencing [14] [15] [16] . At first we consider the Equation (2) and have to find the dimensionless form of this equation. For obtaining the dimensionless equation of (2), let
(20)
where
and l be the length of velocity amplitude and the calculation area (domain). Hence the non-dimensional form of Equation (2) is
If we consider
as a Reynolds number and ignore “-” we can write the above equation as follows:
(21)
Here the limit conditions are periodic and initial condition at
are written as
with the dimensionless variables
. We can observe that the high-speed fluid catches up with the slow-moving one so that to create a velocity break. This observation is known as shock. If we declare the disturbance as
, the slope would have decreased. As we say at first for numerical schematics of this dimensionless Equation (21), we use Crank-Nicolson method for simulation investigation. We can elaborate these partial terms of (21) and obtain
(22)
where
and expanding the term
, we have
Similarly for other terms and using the notations
and
, finally Equation (22) becomes
(23)
Remark 2. Due to some limitations over Explicit method, mainly regarding convergence and stability, the another scheme of Finite-difference method is Crank-Nicolson method. The most common features of this method are as follows:
1) Implicit method, unconditionally convergent and stable;
2) The method derived by introducing a fictitious time level at (
); and
3) Truncation error of order is
, and hence less computation cost.
We can nicely observe the simulation part of this calculation using MATLAB program and some graphical representation. It is noted that everywhere in the numerical simulations, we consider
unless notified the different values.
In both left and right diagrams of Figure 2, let us consider the very small Reynolds number compared to later one and the shock creates and disturbance is erased as soon as it becomes visible. Here
and
are the values of the Reynolds number and the maximum velocity amplitude is 1.0.
Figure 2 demonstrates the shock waves which appears for viscosity,
in the left figure while
for the rest one, respectively, where time scales are equal in both figures.
The similar scenarios’ are observed in Figure 3 for large Reynolds number (
) compared to Figure 2 although the amplitude structures change its directions and density.
The illustration in Figure 4 (left) establishes that the shock arises and the viscosity is too small (
) than the previous two illustrations. The viscosity makes the difference of velocity decrease till there is no break. The
Figure 2. The profile of the solution of Burger’s equation for small Reynolds number (left)
and (right)
at time
over the domain.
Figure 3. Solutions profile of Burger’s equation for Reynolds number (left)
and (right)
at time
over the domain.
Figure 4. The profile of the solution of Burger’s equation for large Reynolds number (left)
and (right)
with time scale
at time
over the domain.
disturbance is increased highly for viscosity which was too small in earlier diagrams, for example in Figure 2 and Figure 3. It creates oscillations around the shock.
If we increase velocity amplitude as
, we can see the following graphical scenario depicted in Figure 4 (right). The graphical illustration of Figure 4 (right) is similar to Figure 3 (left, right) while
, the time scale is small, even the Reynolds number is equal and
in both figures of 4. We get similar behavior of flow if we increase the velocity amplitude for
which is large at time
. Shock arises and the disturbance is removed.
Finally in Figure 5, the viscosity is smaller than the previous diagrams for
. The shock tends to appear and the disturbance is increasing for decreasing of viscosity. More specifically, in Figure 5 (right), we take the value of viscosity
which is too small than Figure 5 (left) and notice that the following schematics time scale is still similar.
Investigating those results, we can say that bigger values of viscosity erase the disturbance in the flow. That is, the disturbance has no more issue, if the viscosity is huge. On the other hand, the smaller the viscosity is, the shock appears and the faster disturbance has more effect. So we have concluded that shock arises for small values of viscosity and produces disturbance as oscillations and also it creates sooner when velocity amplitude is bigger.
4. Travelling Wave Solution of the KdV-Burger’s Equation
We can similarly able to find the travelling wave solutions of KdV-Burger’s equation which is described in this portion shortly using only tan-cot function method [1] [2] . Consider the Equation (4), where a, b and d are arbitrary constants. In order to solve (4) by tan-cot function method, let the wave transformation is
where
and c is the wave speed.
After using the transformation Equation (8), we obtain the ordinary differential equation from (4) such that
(24)
Integrating Equation (24) with respect to s, we have
(25)
Then substitute
and
from (10), the Equation (25) translated to
Figure 5. The profile of the solution of Burger’s equation for very low viscosity of the fluid (left)
and (right)
at time
over the domain.
(26)
Equating the exponents and the coefficients of each pair of the tan functions, we find that
For
, we obtain the following relations from Equation (26)
(27)
and it is easy to find the following solutions
Substituting the values of c and
in (10), the exact solution of the Equation (4) is of the form
(28)
which is the travelling wave solution of KdV-Burgers equation.
It is remarked that the next two figures (Figure 6 & Figure 7) show the significance of different parametric values on exact travelling wave solutions.
The travelling wave solutions pattern are displayed in Figure 6 due to the effect of wave number parameter
and for the fixed constants
. It is seen that the travelling waves are visible with the increasing values of
.
Similarly, Figure 7 shows the effect of constant d when all the parameters are unchanged. It is remarked that wave length is too small and number of wave is easily countable as long as d is very small, see 7 (left) and simultaneously the wave number and length are increasing for larger values of d.
Figure 6. Travelling waves of KdV’s equation for
while (left)
, (middle)
and (right)
at time
over the domain
.
Figure 7. Travelling waves of KdV’s equation for
while (left)
, (middle)
and (right)
at time
over the domain
.
5. Conclusion
In this study, the Burger’s equation is solved by using tan-cot function method. Numerically we present the effect of viscosity on amplitude that corresponds to the Reynolds number. The Burger’s equation would shock up and tend to break without the presence of viscous terms. Mainly viscous term helps to suppress this breaking effect by countering the nonlinearity and the larger viscosity creates the smaller disturbance. Some exact solutions are presented graphically to observe the travelling waves.
Acknowledgements
We thank the Editor and the anonymous referees for their comments and suggestions to rich the manuscript.