Parametric Dirac Delta to Simplify the Solution of Linear and Nonlinear Problems with an Impulsive Forcing Function

The Laplace transform is a very useful tool for the solution of problems involving an impulsive excitation, usually represented by the Dirac delta, but it does not work in nonlinear problems. In contrast with this, the parametric representation of the Dirac delta presented here works both in linear and nonlinear problems. Furthermore, the parametric representation converts the differential equation of a problem with an impulsive excitation into two equations: the first equation referring to the impulse instant (absent in the conventional solution) and the second equation referring to post-impulse time. The impulse instant equation contains fewer terms than the original equation and the impulse is represented by a constant, just as in the Laplace transform, the post-impulse equation is homogeneous. Thus, the solution of the parametric equations is considerably simpler than the solution of the original equation. The parametric solution, involving the equations of both the dependent and independent variables in terms of the parameter, is readily reconverted into the usual equation in terms of the dependent and independent variables only. This parametric representation may be taught at an earlier stage because the principle on which it is based is easily visualized geometrically and because it is only necessary to have a knowledge of elementary calculus to understand it and use it.


Introduction
The purpose of this paper is to present a more rigorous derivation of the parametric Dirac delta than that which was presented in [1], and also to illustrate its application in partial differential equations and nonlinear problems.
According to distribution theory, the Dirac delta is the result of differentiating the Heaviside unit step.The particular parameterization presented in [1] permits this differentiation to be carried out by means of elementary calculus and the resulting pair of parametric equations are exact and closed.
The delta equations have the same function values as those specified in the definition; the area involved has a unit value; they comply with the fundamental property and yield the correct Laplace and Fourier transforms [1].In the solution of differential equations, they are handled exclusively by calculus and algebra, both at an elementary level.The parameterized representation can be read-ily visualized geometrically.These two features should make these parametric equations particularly convenient as a useful research tool, and also, for the purpose of teaching the Dirac delta concept at an early stage in undergraduate school.

Basic Concepts
For the sake of simplicity and readability, in this paper, the Dirac delta will be derived and applied considering it to represent a time concentration, and not a space concentration, and that the point of concentration occurs at time equal to zero, i.e., in the form often referred to as the "Impulse Function".

The Cauchy Limiting Coefficient
Cauchy proposed a limiting coefficient and represented it by the following equation [7]: This coefficient he used to delimit the interval of validity of a function.
It is worth pointing out that both in the third definition of the Heaviside unit step, Equation (1), and in the Cauchy limiting coefficient, Equation ( 2), the value of the jump point is undefined; as a matter of fact, their graphical representation is identical, Figure 1.However, they differ in that the derivative of the Cauchy coefficient is zero, it is not the Dirac delta, this is easily confirmed by Mathematica, Maple and by the TI 92 calculator.

Unit Step with a Riser
Consider a variant of the Heaviside unit step which, unlike it, the jump point is filled with a vertical straight line.We will call this the unit step with a riser, H R , Figure 2.

Derivation of the Parametric Dirac Delta
Consider the approximation of the unit step with a near vertical riser, H Ra , shown in Figure 3.
It is clear that: From Figure 3, the equation for H Ra is easily established: Notice that, in this paper, λ is used as a switch, i.e., to switch on functions at the beginning of their interval of validity and to switch them off at the end of their interval of validity.In this manner, various different functions are linked together into a single composite function.

Parametric Representation
Now: Carrying out the slash division in a piecewise fashion yields: tion of H Ra , with the length of the curve w as the parameter.
From Figures 4 and 5 we obtain the equations: or simply and which means: (13) These two functions would be continuous were it not for the fact that they are undetermined at the points and 0 w  , however, since their left limit is the same as their right limit at those points, they will be treated as if they were continuous because this "••• is generally inconsequential in applications", [6], see also [4,5,7].This is true also for most of the functions contained in the equations in the rest of this paper.It is significant that this conversion of a discontinuous function into a "continuous" function by means of a parameterization has been used also to eliminate the Gibbs phenomenon [8].but this is the Dirac delta, But also, d d Ra a H t may be expressed as: taking the limit yields: in accordance with Equation ( 14): From this point on, all the derivatives have been verified by using Mathematica.
17) Differentiating Equations ( 5) and (6): Taking the limit of as a t 0   of Equation (6) yields: Thus the pair of Equations ( 17) and ( 18) is the parametric representation of the Dirac delta.However, an abbreviated version of these equations for the application referring to the impulse function will be established in Section 3.4.

Illuminating Plots
Some very illuminating plots result if we invert the order of the operations on parametric Equations ( 5) and ( 6), i.e., we first carry them to the limit as 0 and (d), and we differentiate them afterwards, Figures 6(e) and (f).
Step.Carrying them to the limit yields: see Figures 6(c) and (d).

Displaced Point Plots
If we now plot  

Plot of a True Single Point
It is interesting to compare Figures 7(a  that the size of the displaced point remains the same, since it really refers to a single point.
Considering that there is no negative time, the plots of Figure 6 are conveniently substituted by those of Figure 9.
From Figure 9 we obtain the following abbreviated equations: Thus the remarkably simple pair of parametric Equation (28) represent what is often called the "impulse function".These are the equations to be used in the solution of problems.

Example 1
Consider a one dimensional rod subject to an impulsive heat source with initial temperature of 0˚C along its full length and with the ends kept at 0˚C throughout the whole process.The governing equation is: (here q has units of energy/volume) Subject to the boundary conditions: and to the initial condition:  mass density.Following the method of separation of variables [9]: Substituting Equations ( 25) and (33) into Equation (30) yields: Introducing the parameter w into Equation (34): Open Access JAMP   Substituting Equations ( 26) and (27) into Equation ( 36), yields what we will call the control equation:

Impulse Instant
During the impulse instant, designated as interval i, , , Equation ( 23), accordingly Equation (37) becomes or according to Equation (33): Comparing Equation (39) with the differential Equation (30) it stands out that, during the impulse instant, the term referring to conduction is absent.This is as it should be, since there is no time for conduction to take place, in perfect agreement with physical reality.Also the Dirac delta, as such, is absent and the forcing function is simply a constant.
Dividing Equation (38) by where is a separation constant, then Substituting Equations (42) into Equation (33): At the "beginning" of the impulse instant, w = 0, and according to the initial condition, Equation (33), Substituting this into Equation (43) requires that: and thus Equation (43) reduces to: Equation ( 45) governs during the impulse instant, thus: the change of temperature during the impulse instant is , substituting these values in Equation (47) yields: since m remains constant, we have the principle or in words: "The heat impulse is equal to the change in sensible heat".This is somewhat similar to the mechanical impulse and change in momentum principle.

Post Impulse Time
At post impulse time, designated as interval p, w ≥ 1, , Equation (23), and thus Equation (37) becomes: Physical considerations require that: To make things clear, it is convenient to emphasize that: But in this interval, , substituting this and making use of Equation (33), Equation ( 50) is equivalent to the following homogeneous equation: and in view of the second of Equations( 46) and Equation (51), the initial temperature of post-impulse time is: and, of course, the boundary conditions are the same as Open Access JAMP the original ones, Equation (31), thus: The problem made up of Equations ( 53)-( 55) has the well known conventional solution, see for instance [9]:

The Complete Solution
Using Equations ( 45) and (56) the complete parametric solution is obtained thus: Notice that the first term of Equation ( 57) is the impulse instant solution and the second term is the post impulse solution.
Notice that since t is a function of w, both Equations (57) are functions only of the parameter, w.See Figure 10.
It is interesting to compare the plot of the parametric solution, Equation (57), Figure 10( which is not in agreement with the specified condition, Equation (32).Also in Figure 11 the pseudo initial condition is plagued by the spurious oscillations due to the Gibbs phenomenon.

Example 2
Consider a mass-spring non-linear system subjected to an initial impulse where the spring force is given by: as in the Duffing equation.    with initial conditions: I is the magnitude of the impulse.
Converting the second order Equation (59) into two first order equations by means of the state variables: Introducing the parameter w into Equation (61): Open Access   multiplying by d d substituting Equations ( 26) and (27) into Equation (63) yields the control equations: (64)

Impulse Instant
During the impulse instant: , Equation (23), and thus Equations(64) become: (The subscript i refers to the impulse instant).
Notice that the parameterization resulted in the simplification of the forcing function which is now simply the constant I, i.e., the Dirac delta, as such, is absent from the previous equations.Also absent is the term referring to the spring force, thus representing faithfully the physical reality, i.e., during the impulse instant, there is no time for the displacement of the mass, and thus the spring does not participate in this process.
Integrating Equations (65): from the initial conditions the integration constants are evaluated, thus: and, therefore, the governing equations for the impulse instant are: and so at the "end" of the impulse instant, i.e., at 1, w  0 t  , Equation (23): The 2 nd of Equations (67) states that there is no motion during the impulse instant.Now This is the principle of impulse and momentum and it was applied automatically by virtue of the parameterization.

Post Impulse
It is convenient to reduce the pair of Equations ( 72) into the single equation: Resorting to Equations (68) we establish the initial conditions of post-impulse time: the phase-plane solution of Equation ( 73) is , , In terms of these variables Equation (75) becomes In the first of Equations (78) the first term refers to the impulse instant and the second term refers to the post impulse time.
See Figure 12.Equation (77) leads to: the solution of Equation ( 79) is the hypergeometric function see Figure 13(a).

Conclusions
The parametric representation of the Dirac delta as proposed in [1] has been reviewed with the purpose of clarifying the concept.
The parametric Dirac delta, contained in the forcing function of a differential equation referring to an impulsive process, has an operator action that splits this equation into two: the first referring to the impulse instant and the second referring to post-impulse time.In the impulse instant, the operator action converts the delta forcing function into a constant; and furthermore, it cancels the terms referring to the phenomena that cannot take place instantaneously.
Also, the operator action makes the post-impulse equation homogeneous.All of this makes the process of obtaining the solution considerably easier.As it has been mentioned before, in the case of Example 1 referring to the metal rod heated impulsively, the term associated with heat conduction disappears from the impulse instant equation.This is reconciled with the physical reality, i.e., there is no time for this process to take place.
In the case of the second order mass-nonlinear spring problem of Example 2, the term representing the spring force disappears from the impulse instant equation.This illustrates the faithfulness of the mathematical model, since in the impulse instant there is no time for the displacement of the mass to take place.
In Examples 1 and 2, the parametric representation permits the separation of variables even though the original differential equations are non-homogeneous.
In the parametric solution, the real initial conditions are used and the processes that take place during the impulse instant are fully represented as such.In contrast to this, in the solution of problems involving impulsive processes in some textbooks, the forcing function is sometimes replaced by equivalent initial conditions [10].
During the impulse instant: time, t, does not flow, but pseudo-time (parametric time), w, does flow and this is what makes possible the establishment of the equations that represent the processes that take place instantaneously.
In the parametric solution, the impulse instant becomes a finite interval in terms of the parameter.Thus there is no need to deal with infinitesimals.
This parametric representation may be taught at an earlier stage because the principle on which it is based is easily visualized geometrically and because it is only necessary to have a knowledge of elementary calculus to understand it and use it.
By virtue of the parametric representation of the Dirac delta, the principle of impulse and momentum was applied automatically in Example 2, also in Example 1 a similar result was obtained: the heat impulse is equal to the change in sensible heat.This suggests that the parametric representation of the Dirac delta may turn out to be a valuable research tool.

Figures 4
Figures 4 and 5 together are the parametric representa-

Figure 1 .
Figure 1.Graphical representation of both the third definition of the Heaviside unit step, Equation (1), and the Cauchy limiting coefficient, Equation (2).Note: Matlab was used for this plot because, unlike other software, it does not leave a trace at the jump point.

Figure 2 .
Figure 2. Unit step with a riser.

Figure 3 .
Figure 3. Approximate unit step with a riser.

Figure 4 .
Figure 4. Unit step with a near vertical riser H Ra as a function of the parameter w.

Figure 5 .
Figure 5. Approximate time t a versus the parameter w.

Figure 6 (
a) illustrates the usual unit step with a gap at the jump point.

Figure 6 (
b) shows the fundamental idea of a unit step with a riser, the parameter w is the length along this function.Figures6(c) and (d) represent the parameterized unit step with a riser, these are the plots of the resulting virtually continuous functions.It is significant that the single point, , of Figure 6(b) has been expanded into the finite interval, , in Figures 6 t, and not versus w, the finite intervals, from 0 to 1, of Figures6(e) and (f) become the single points 0 t  of Figures7(a) and (b).As can be seen in Figure7(c), even though, the scale of the ordinate is very compressed (up to 10 × 10 30 ), no displaced point appears since it is

Figure 6 .
Figure 6.(a) Plot of either the Heaviside unit step, H(t − 0), third version, Equation (1), or Cauchy's limiting coefficient, λ(t − 0), Equation (2).(b) Parametric plot of the unit step with a riser, Equation (19) versus Equation (20).(c) Plot of Equation (20).(d) Plot of Equation (19).(e) Plot of Equation (22).(f) Plot of Equation (21).Note: Matlab was used for these plots because it makes a clear distinction between the step with a riser and the step without a riser and so does the TI 92 graphics calculator.A plot of the riserless step, Figure 6(a), in another software or another graphics calculator may result in a trace at the jump point making it indistinguishable from the step with a riser.

Figure 7 .
Figure 7. (a) Plot of Equation (21) versus Equation (20).(b) Plot of Equation (22) versus Equation (20).The plots (a) and (b) may be called displaced point plots.The circles help to locate the gaps which are very narrow and the points which are very faint.(c) Plot of Equation (17) versus Equation (18).located at infinity.
) and 8, both are plots of Equation (21) versus Equation (20), but the plotting increment of Figure8is much greater than that of Figure7(a), and consequently the gap of Figure8is also much greater than that of Figure7(a).Notice however,

Figure 8 .
Figure 8. Direct plot of Equation (21) versus Equation (20), with a very large plotting increment.Notice that the gap is much larger than in plot 7(a), however the displaced point remains the same size and just as faint since it is truly a single point.
a), with the plot of the conventional solution Equation (56), Figure 11.The parametric solution clearly shows the initial condition to be zero, i.e., in accordance with the specified condition, Equation (32), and, furthermore it shows the initial, instantaneous process of the change of temperature, i.e.: the plot of the conventional solution shows the initial temperature to be

Figure 10 .
Figure 10.Plots of the parametric solution, Equations (57).(a) "Front view" plot.Notice that the instantaneous change of the value of the temperature from its initial value, 0, to 1 (in non-dimensional terms) is clearly shown by the front "wall" which is not present in the plot of the conventional solution, Figure 11.(b) "Rear view" plot.

Figure 11 .
Figure 11.Plot of the conventional solution, Equation (56).The value of the initial temperature was given as zero, but in this plot it is shown as one (dimensionless) and the instantaneous change of temperature process that takes place during the impulse instant (the "wall" of Figure 10(a)) is not shown.The oscillations in the value of the pseudo initial temperature along the full length of the rod are spurious, they are due to the Gibbs phenomenon.Compare this with the plot of the parametric solution, Figure 10.
phase-plane solution may be expressed as:

Figure 12 .
Figure 12.Plot of the phase-plane solution.