Travelling Waves: Interplay of Low to High Reynolds Number and Tan-Cot Function Method to Solve Burger’s Equations

We study the nonlinear parabolic equations for travelling wave solutions of Burger’s equations. The purpose of the present work is to study various types of Burger’s equations describing waves and those are based on nonlinear equations. We focus on to describe the analytic solution in the special pattern of travelling wave solutions using tan-cot function method. We discuss about inviscid and viscous version of Burger’s equation for fluid flow and investigate the effects of internal friction of a fluid via Reynolds number. By changing the velocity amplitude, the nature of flows with shock wave and disturbance are observed. For numerical solutions, the Crank-Nicolson scheme is introduced to establish the wave solutions.


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 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 where d is the viscosity of the fluid and g is a nonlinear function of u.
Introducing the function ( ) g u , the Burger equation can be written in the form of here u is the fluid density, d is the viscosity, the spatial variable x and time is t.In absence of viscosity, 0 d = , and the Equation (2) becomes

( )
G s 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 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 1 a = and 0 k = 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.

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] where ( ) , , u t y x is a travelling wave solution of (5).Now using the travelling wave transformation, Journal of Applied Mathematics and Physics ( ) ( ) where s mx y nt = + + .Here m and n are real constant and transformation gives us following changes By integrating Formula ( 7) with respect to s, a non-linear partial differential Equation ( 5) transforms to a nonlinear ordinary differential equation such that 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: e G e G e s s + + + = where ( ) to be determined and the integral constant 0 e .
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 where , , α µ δ are parameters and to be determined and µ is known as a wave number.
Remark 1.To get a function

( )
G s 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 Similarly, for cot function, we have 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.

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.

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) Integrating Equation ( 12) once with respect to s, we have Substituting ( )

G s and d d
G s from (10), the Equation ( 13) implies 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) which is the travelling wave solution of one dimensional Burger's equation.
and on integration with respect to s, we obtain ( ) where A is a constant.For equilibrium points ( )( ) Hence the roots are G G > then the travelling wave solutions are possible if ( )( ) integration in (17), we have After solving this we get After that multiplying and dividing by This solution is similar to shock wave profile because of joining the asymptotic states of G 2 and G 1 with the boundary conditions

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 where 0 u and l be the length of velocity amplitude and the calculation area (domain).Hence the non-dimensional form of Equation ( 2) is as a Reynolds number and ignore "-" we can write the above equation as follows: Here the limit conditions are periodic and initial condition at 0 t = are written as ( ) . 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 Journal of Applied Mathematics and Physics 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 where and expanding the term , we have ( ) Similarly for other terms and using the notations t K ∆ = and x h ∆ = , finally Equation ( 22) becomes ( ) ( ) 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 ( 1 2 j + ); and 3) Truncation error of order is O t O x ∆ + ∆ , 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 0.1 h K = = 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 0.1 Re = and 0.25 Re = are the values of the Reynolds number and the maximum velocity amplitude is 1.0.We get similar behavior of flow if we increase the velocity amplitude for 0.01 d = which is large at time 20 t = .Shock arises and the disturbance is removed.Journal of Applied Mathematics and Physics 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 0.0003 d = 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.

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 and c is the wave speed.
After using the transformation Equation ( 8), we obtain the ordinary differential equation from (4) such that Integrating Equation ( 24) with respect to s, we have

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 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.

Figure 1
Figure 1 depicted the physical structure of travelling wave which counts the wave number and distributions.If we use the transformation ( ) ( ) , u t x G s = ,

Figure 1 .
Figure 1.Travelling waves of Burger's equation for 1 d c = = while (left) 0.5 µ = , (middle) 1.0 µ = , and (right) 8.0 µ = at Initial value 0s is not essential and the derivative of ( )Gs 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 0 d = .So the wave solution is

Figure 2
Figure 2 demonstrates the shock waves which appears for viscosity, 10 d =

4 d
= 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 ( 10 &100 Re = ) 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 ( 0.005 d = ) than the previous two illustrations.The viscosity makes the difference of velocity decrease till there is no break.The

Figure 2 .
Figure 2. The profile of the solution of Burger's equation for small Reynolds number (left) 0.1 Re = and (right) 0.25 Re = at time 20 t = over the domain.

Figure 3 .
Figure 3. Solutions profile of Burger's equation for Reynolds number (left) 10 Re = and (right) 100 Re = at time 20 t = over the domain.

Figure 4 .Figure 4 (
Figure 4.The profile of the solution of Burger's equation for large Reynolds number (left) 200 Re = and (right) 200 Re = with time scale 0.01 K = at time 20 t = over the domain.

Figure 5 .
Figure 5.The profile of the solution of Burger's equation for very low viscosity of the fluid (left) 0.001 d = and (right) 0.0001 d = at time 20 t = over the domain.

(
and the coefficients of each pair of the tan functions, Substituting the values of c and α in (10), the exact solution of the Equation (4) is of the form 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 Figure6due to the effect of wave number parameter µ and for the fixed constants 1seen that the travelling waves are visible with the increasing values of µ .Similarly, Figure7shows 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.