Boundary Value Problems for Burgers Equations , through Nonstandard Analysis

In this paper we study inviscid and viscid Burgers equations with initial conditions in the half plane x R T ∈ > , 0 . First we consider the Burgers equations with initial conditions admitting two and three shocks and use the HOPF-COLE transformation to linearize the problems and explicitly solve them. Next we study the Burgers equation and solve the initial value problem for it. We study the asymptotic behavior of solutions and we show that the exact solution of boundary value problem for viscid Burgers equation as viscosity parameter is sufficiently small approach the shock type solution of boundary value problem for inviscid Burgers equation. We discuss both confluence and interacting shocks. In this article a new approach has been developed to find the exact solutions. The results are formulated in classical mathematics and proved with infinitesimal technique of non standard analysis.


Introduction
The nonlinear parabolic partial differential equation was first introduced by J. M. Burgers [1] [2] as the simplest model for fluid flow, this equation combining both nonlinear propagation effects and diffusive effect.If ε is not null, we approach to the Naveir's Stokes equations in one dimension.Burgers equation has a wide variety of applications in the modeling of water in unsaturated soil, dynamics of soil water, statistics of flow problems mixing and turbulent diffusion cosmology and seismology.

S. Bendaas
When ε is null, this equation approaches to the Euler's equations in one dimension who governs the flows of perfect fluids.It's the viscid equation.it has the form 0.
If the viscous term is dropped from the Burgers equation, discontinuities may appear in finite time; even if the initial condition is smooth, they give rise to the phenomenon of shock waves with important application in physics [3].This property makes Burgers equation a proper model for testing numerical algorithms in flows where severe gradients or shocks are anticipated [4]- [6].Discretization methods are well-known techniques for solving Burgers equation.Ascher and McLachlan established many methods as multi-symplectic box sheme.For the boundary value problem, Sinai [7] was interested to the initial condition case: null on R − and Brow- nian on R + .She, Aurell and Frich [8] with numerical calculations particularly examined the initial conditions of Brownian fraction nair to the asymptotic behavior.
A remarkable feature of viscid Burgers equation is that its solutions with initial conditions of the form ( ) ( ) can be explicitly written down.Hopf [9] and Cole [10] independently showed that the Equation (1.1) can be linearized through the transformation 2 .
x u εϕ ϕ = − (1.4) Then Hopf [9] showed that if ( ) Then ( ) , u x t defined by Equation (1.4) solves (1.1) and (1.3).Conversely, if ( ) , u x t is a solution of problem (1.1) and (1.3) then ( ) , x t ϕ defined by Equation (1.4) is a solution of problem (1.5) and (1.6).Solving for ϕ from (1.5) and (1.6) and substituting it into Equation (1.3) we obtained explicit formula for the solution of problem (1.1) and (1.2) namely: and studied the asymptotic behavior of ( ) , u x t .Explicit solutions of the Burgers equation (1.1) in the quarter plane with integrable initial data and piecewise constant boundary data were constructed by [11] using Hopf-Cole transformation [9].He obtained a formula for its weak limit as viscosity parameter goes to 0. Although, there are many results for initial value problem has been studied less.Using maximum principle, this formula for weak limit was extended to general boundary data.ε is a positive parameter small enough.The problem is considered by [12].As the fact that ε multiplies the largest derivative, one is in the presence of a singular perturbation problem.The purpose of Singular Perturbation Theory is to investigate the behavior of solutions of (1.1) as 0.
ε → The aim of the present article is to study solutions of Inviscid and Viscid Burgers equation if the initial condition admits several singular points, i.e. in the case of a finite number of shocks.A simple formulation is given for the asymptotic behavior based on the evaluation of integrals which is a method of the non standard perturbation theory of differential equations proposed by Imm Van Den Berg [13] and improved by Lutz and Callot.
Historically the subject non standard was developed by Robinson, Reeb, Lutz and Goze [14].The nonstandard perturbation theory of differential equations, which is today a well-established tool in asymptotic theory, has its roots in the seventies, when the Reebian school (see [14] [15]) introduced the use of nonstandard analysis into the field of perturbed differential equations.Our goal in this paper is to generalize these techniques on EDP and our general purpose is to describe the asymptotic behavior of solutions in boundary value problem with a small parameter ε and to discuss in particular the cases of confluence and the interacting shocks with new technical infinitesimal of non-standard analysis.We can conclude that the solutions of the problem: (1.1) and (1.3) are infinitely close to the solutions of problems (1.2) and (1.3), as ε is a parameter positive sufficiently small.
In Section 2, we treat the boundary value problem for inviscid Burgers equation, solve it and study it.Section 3 is devoted to useful lemmas for our main results.In Section 4, we study viscid Burgers equation, solve exactly the initial value problems for it, and describe the asymptotic behavior of solutions with a non standard form.Some components, such as multi-leveled equations, graphics, and tables are not prescribed, although the various table text styles are provided.The formatter will need to create these components, incorporating the applicable criteria that follow.

Shock Fitting
We consider the inviscid Burgers equation: 0.
is continuous.This problem not admits the regular solutions but some weak solutions with certain regularity exist.The Burgers equation on the whole line is known to possess traveling wave solutions.The solution of (1.2) and (1.3) may be given in a parametric form and shocks must be fitted in such that: where 1 ξ and 2 ξ are the value of ξ on the two sides of the shock [16].
According to Equation (1.2), the solution at time t is obtained from the initial profile ( ) This is the differential equation for the line cord of shock that checks the condition of entropy such as [16] [17] Corresponding to the two inflection points.

Confluence of Shocks
When a number of shocks are produced, in general it is possible for one of them to overtake the shock ahead.Then they combine and continue as a single shock.This is also described by our shock solution.
Consider the curve given by f Figure 1, then two shocks are formed corresponding to the inflection points p and q with families of equal area chords, typified by 1 2 p p and 1 2 q q .As time goes, the points 1 q and 2 p approach each other until the stage in Figure 2 is reached where a common chord cuts off lobes of equal area for both humps.
At this stage the characteristics corresponding to 2 p′ and 1 q′ are the same and therefore the shocks have just combined into one as shown in Figure 3. Characteristics 2 p′ and 1 q′ combined.All the characteristics be- tween 2 q′ and 1 p′ have now been absorbed by one or other of the shocks.A single shock proceeds using chords 1 2 p q ′′ ′′ .
In the plane ( )

Preliminaries
In this section we present some lemmas that are important to prove our main result.Proposition 3.1.Let ϕ ϕ ϕ = + be the analytic solution for the initial value problem for the heat equation (1.5) and (1.6).If the initial condition is given as in Figure 1, then u given by the Formula (1.4) is a solution for the initial value problem (1.1) and (1.3).It is explicitly given by: Proof: When a shock overtakes another shock, they merge into a single shock of increased strength as described in inviscid solution ( ) 0 ε = , on the f curve in Figure 1; It is possible to give a simple solution of Burg- ers equation that describes this process for arbitrary ε .The solution for a single shock is given in [3] [4] and the corresponding expression for ϕ may be written in the form: In the expression of solution for a single shock given in [3] [4], the parameters b₁, b₂ witch locate the initial position of the shock are taken to be zero.The expressions 1 2 ϕ ϕ are clearly solutions of the heat equation ( ) Corresponding to the initial conditions: ( ) ( )

S. Bendaas
Then the solutions of the heat equation are given as: Using Equations (1.4) and (3.1) we obtain the expression (3.2).And to prove our results, we use the non standard analysis techniques, that we consider the following lemma.
Lemma 3.2.(The Van.Den.Berg lemma [14]): Let h be a standard function, defined and increasing on where δ for v ≃ 0. ( ) ( ) And let ϕ be an intern function de-  fined on ]0,+∞[ such that : ( ) ( ) where a, r are positive standard, m and q are the both positive δ is an infinitesimal.b and s are standard, To give estimation to the solution, given by (3.2), we state the following lemma: Lemma 3.3.Let ε be a positive real small enough.And let ϕ and h be two standard functions such that: h, is a C 2 class function verified the Lemma 3, and admits on the ξ point a unique absolute minimum It is S-continuous on ξ and satisfies the conditions of the Lemma 3 in the both ways.Then δ is an infinitesimal.Proof: To prove this lemma, we use the "Van Den Berg" method, lemmas: (5.6), (5.7) [14].It consists in the following steps 1) Search for the absolute minimum (maximum) of the function under the exponential sign and bring it out.
2) Bring back the minimum (maximum) to the zero.
3) Searching the galaxy as well as the main galaxy where the function in the exponential sign is appreciable.4) Calculate the integral.
As consequence we have the following lemma.Lemma 3.4.Let f the initial condition as Figure 1.Assume: Then for x and t fixed, the functions defined as: has at most two minima i ξ + and i ξ − relative to the variable i η .These two minima satisfy the equations: ( ), Proof: Let f the initial condition given as in Figure 1, two shocks are formed corresponding to the inflection points of ( ) f ξ .for some time we see appear an area of three values for the solution ( ) , u x t .At first the wave breaks on the feature for which: Inside the zone of each shock and for a point ( ) , x t , there are three characteristics corresponding to two minimum framing a maximum.When a shock overtakes another, they merge into a single shock of increased strength as describe for the inviscid solution ( ) , each couple frames a maximum.Let: ( ) , be the minimum for the functions given by the Formula (3.10), is such that: and within (3.10) this condition can be written as: But i ξ − and i ξ + both verify the equations: ( ) The condition of the shock is expressed by (3.12), is the same condition of shock given by (2.2) for invicid Burgers equation.

Confluence of Shocks
Our general purpose now is to show that the exact solution of (1.1) and (1.3) endorse the ideas regarding shocks in Section 2, we want to confirm that as ε is small enough, the solution of (1.1) and (1.3) reduce to solution of (1.2) and (1.3), with discontinuous shocks which satisfy the condition (2.2), and the shocks are located at the positions determined in Section 2. The shocks are formed corresponding to the inflection points of the initial condition ( ) ( ) , who assume the assumptions (H 1 ), (H 2 ) in the lemma (3.4).Then we proved the following result: Theorem 4.1.Under the assumptions: (H 1 ), (H 2 ) in lemma (3.4), the problem (1.1) and (1.3) admits a unique solution for 0 t > given by:

Such a solution is confluence of shocks and for ε sufficiently small, this solution is infinitely close to the solution of the reduced problem given in (2.2).
Proof: 1) From Figure 2 and Figure 3, after some time the two shocks combine and continue into one, and this is the lowest minimum that carries this amount to the case of a single shock.From the proposition (3.1), the problem (1.1) and (1.3) admits a single solution explicitly given by Formula (4.1).Uniqueness is due to the condition of entropy which restricts the set of solutions to one which is stable with a singular perturbation dissipative nature.
2) Let ( ) , x t be a standard point outside the line of the shock.From the solutions of the equation , there is only one denoted i ξ is the absolute minimum of the function given by the expression (3.10).Using the expressions (3.1) and (3.2) of the Proposition 3.1 ( ) , u x t is given as: From the Lemma 3.2 we have where δ > 0 is an infinitesimal.And we will have the following estimate To conclude we have the following corollary.3) verifies at the point ( ) ( ) Proof: Using lemma (3.4), outside the region of each shock.For ( ) ϕ and i i u ϕ are given by (3.6) and (3.7).Then ( ) And it follows that: ϕ function dominates when x is infinitely large positive and there is obtained ϕ function dominates when x is infinitely large negative and there is obtained

Interacting Shocks
In this section, we discuss the interacting shocks case; before going further in this case we need the following proposition and lemma.Now since any i ϕ a solution of the heat equation, we may clearly add further terms in (3.3) and generate more general solution of burgers' equation.Such solution represents interacting shocks.As consequence we have the following Proposition 4.3.Let , be the analytical solution for the problem (1.5) and (1.3).If f is the initial condition admitting three inflection points, then the solution of (1.1) and (1.3) is explicitly given by ( ) Proof: When a shock overtakes another, they merge into a single shock of increased strength as described in inviscid solution ( ) For arbitrary ε , it is possible to provide a simple solution to the Burgers equation that describes this process.The solution for a single shock is given in [4] and the corresponding expression for ϕ is written as: .
In the solution for a single shock given in [4], the parameters 1 2 , b b which locate the initial position of the shock are taken to be zero and ( ) And S. Bendaas Corresponding to the initial conditions: ( ) ( ) Then the solutions of the heat equation are given as: and ( ) ( ) Using (1.4) and (4.4) we obtain the expression (4.3).Then we have the following.Theorem 4.4.For 0 t > , and for 3 2 1 0 u u u > > > , the problems (1.1) and (1.3) have a unique solution given by: ( ) Such a solution is interacting shocks and for ε sufficiently small, it is infinitely close to the solution of the reduced problems (1.2) and (1.3).
Proof. 1) In the interacting shock case we have three shocks, when a shock overtakes another they merge into a single shock of increased strength and the lowest minimum dominating.Then we go back to the single shock case.Using the proposition (4.3), we deduce the uniqueness of solution explicitly given by (3.4).The uniqueness is due to the entropic condition [3] [4], which restricts the set of solutions to one, who is stable with singular perturbation with dissipative nature.
2) Let ( ) , x t be a standard point outside the line of the shock.From the solutions of the equation , there is only one denoted that i ξ is the absolute minimum of the function given by the expres- sion (3.10).( ) From the Lemma 3.2 we have: where δ > 0 is an infinitesimal.And we will have the following estimate from which the following corollary Corollary 4.5.Let ( ) , x t be a standard point outside the line of the shock.So there exists a unique solution denoted i ξ which is the absolute minimum of the function given by (3.10).If ε is reasonably small, we can recognise shock transition between the states 1 2 3 , , u u u by noting in which regions the corresponding ϕ do- minates.At , , , in 0 , u x t u x t f x , , , in 1 0, u x t u x t f x ξ = − < <  ( ) ( ) ( ) where , i i i u ϕ ϕ are given by (4.8) and (4.9); and using the Lemma 3.2 there are equivalent to those ( ) ( ) ( ) where δ is an infinitesimal positive real.
As shown in Figure 5, If ε is reasonably small, we can recognise shock transition between the states , , u u u by noting in which regions the corresponding ϕ dominates.At 0 t = , in: 0 x < , 1 ϕ dominates and we have:  Since, ( ) ( ) the discontinuity line, it is a straight line chord property still holds.The cord on the f curve cuts off lobes of equal area.The shock determination can then be describe entirely on the fixe ( ) f ξ curve by drawing all the chords with the equal area property can be written analytically as between the points 1 ξ .Moreover since areas are preserved under the mapping, the equal area x t the shocks can be represented by Figure 4.

Figure 1 .
Figure 1.Graphic representation of the initial condition f.

Figure 2 .
Figure 2. Graphical representation from the merger of the two shocks.The characteristics 2 p′ and 1 q′ merge.

Figure 3 .
Figure 3. Construction for merging shocks in a final stage.

Figure 4 .
Figure 4.The (x, t) diagram for merging shocks corresponding to Figure 1.

2 ):
There exist a, b, c, d and e in R, with a < b < c < d < e, such that

Figure 2 . 3 .
The characteristics between 2 q′ and 1 p′ are ab- sorbed by one Figure At this stage there are two stationary values that satisfy the equations (3.11) we noted them by i

Corollary 4 . 2 .η
Let ( ) , x t , be a standard point outside each line of shock.Among the solutions of the equation given by (3.10), and further the solution of (1.1) and (1. fixed, outside the shock region and by the use of the lemma (3.2), each function i h ad- mits one absolute minimum on the point i ξ .In Formula (3.4),( )

, 3 ϕ
dominates and we obtain: The symbol "≃" means infinitely close to[8].Thus we have a shock transitions from 1 u to 2 u centred at 0 x = , and one from 2 u to 3 u centered at 1 x = − .For 0 t > , for early times the transition from 1