Applications of Non-Classical Equations and Their Approaches to the Solution of Some of Classes Equations Arise in the Kelvin-Helmholtz Mechanism and Instability ()
1. Introduction
Many authors (e.g., [1] [2] [3] [4] [5], and the references given therein) investigated mixed type equations and equations of Keldys type (e.g., [6] [7] [8] [9], and the references given therein). Note that a parabolic equation and a mixed equation with changing time direction also have physical applications. The boundary value problems with such sewing conditions appear when modeling, for example, the process of interaction between two reciprocal flows with mutual permeating or when designing certain heat exchangers. Frankly speaking, forward-backward equations (equations of changing time direction) arise in supersonic dynamics, boundary layer theory, and plasma. Therefore, the boundary value problems for equations of mixed hyperbolic-elliptic type with changing time direction (or forward-backward equations) (e.g., [10] [11] and the references therein) present important objects for all investigators. Non-classical models are described in particular by equations of mixed types (for example, the Tricomi equation), degenerate equations (for example, the Keldysh equation), Sobolev type [12] equations (also the Barenblatt-Zsolt-Kachina equation), equations of mixed type with changing time direction, and forward-backward equations. At the same time, we include that the boundary value problems for equations with degeneracy belong to “non-classical” problems of mathematical physics. One of the main difficulties that arise in the theory of degenerate elliptic equations is related to the influence of lower (in the sense of the theory of regular elliptic operators) terms of the equation on the formulation of boundary value problems and their coercive solvability. In this case, degenerating ordinary differential equations (DODE) of second order have physical means as boundary layer as the same physical means with which such sewing conditions appear when modeling, for example, the process of interaction between two reciprocal flows with mutual permeating. These degenerating ordinary differential equations also do not fit into classes of ordinary differential equations (ODE), and may be called non-classical ODE. It means that, in an analogical one-dimensional case, the non-classical equations having the same boundary layer as for ODE are said to be non-classical ordinary differential equations (NODE). This definition of a non-classical ODE can be justified below the presented physical and analytical descriptions at the first entry in this work. The investigators may carry out citation non-classical ordinary differential equations based on the results of the considered problems, as clearly justified in this paper and introduction of this paper subdivided the subsection, in order to show research background (in subsection 1.1), situation of the former researches analyze and the problems existed of works will be summarized (subsection 1.2) and problems we would like to solve in the paper will be listed (subsection 1.3).
1.1. Non-Classical Ordinary Differential Equations Second-Order and Applications Fields
In the mathematical analytical description process of investigation, first of all we need to give historical information for our research into matters. If the degeneration ordinary differential equations (DODEs) second order has physical meaning (for example, boundary layer), analogical to the one-dimensional case of the non-classical equation of mathematical physics is said to be non-classical ODE second order. This definition of non-classical ordinary differential equations of second order entered for the first time by author M.A. Nurmammadov, is proved in the following applications. The fundamental results for degenerate elliptic equations are due to M. V. Keldysh [6]. The results obtained by him were then developed and generalized by O.A Oleinik [13]. Generalized solutions of degenerate second-order elliptic equations were studied in the works of S. G. Mikhlin [14] and M. I. Vishik [15]. For this paper in order we use the main ideas in the derivation of the Friedrichs inequality with minor ones are transferred to the study of the positivity of operators, and obtain for NODEs the Fridric’s-Pouncare inequality, in order to applying for our consideration boundary value problems of ordinary differential equations.
In 1985, the author [4] first considers a new equation and four well-posed boundary value problems, which leads to the equation presented in the form: as it is shown in the work [7] - [17], the theory of boundary value problems for degenerate equations and equations of mixed type, it is a well-known fact that the well-posedness and the class of its correctness essentially depend on the coefficient of the first order derivative (younger member) of equations. Therefore, non-classical ordinary differential equations (NODEs) with second order or for degenerating ODEs also satisfy these facts. In the case of well-posedness, the coefficient of the first order derivative (younger member) of the ODE with second order determines the class of correctness. Accordance to the work from degenerating nonlinear (of power type) mixed type equations considered by author formulated new nonlinear and non-classical ordinary differential equation (NODE) in this paper with physical means (for linear and nonlinear cases is shown in Figure 1 and Figure 2) is expressed. Since the applications of the author of the works are first applied to the problems of astrophysics, planetary and solar systems in space, and a number of phenomena may not be sufficiently studied in works close to this one, there are a large number of open applications that can be an interesting direction for researchers on a global scale, such as new approaches indicated by the scientific field after the appearance of this article. In the Planetary boundary layer, considering object investigations, this work can include the following field applications. Within the lowest portion of the planetary boundary, a semi-empirical log wind profile is commonly used to describe the vertical distribution of horizontal mean wind speeds. The simplified equation that describes it could be from ordinary differential equations. Due to the limitation of observation instruments and the theory of mean values, the levels (z) should be chosen where there is enough difference between the measurement readings. If one has more than two readings, the measurements can be fit to the above equation to determine the shear velocity. The stellar atmosphere is the outer region of the volume of the star, lying above the stellar core, the radiation zone, and the convection zone. The stellar atmosphere is divided into several distinct regions: an unsolved problem in stellar astrophysics is how the corona can be heated to such high temperatures. The answer lies in the magnetic field, but the exact mechanism remains unclear. Another side is also the Kelvin-Helmholtz instability (KHI), which is the instability of a fluid that occurs when there is a velocity shear rate in one continuous fluid or a velocity difference at an interface between two fluids. Therefore, accounting for the above information, we say that Kelvin-Helmholtz instabilities are visible in the atmospheres of planets and moons, for example, in the formation of the Earth or the Red Spot on Jupiter, and in the atmospheres of the Sun and other stars (e.g., [18] - [24] and therein). Shear flows in plasmas are attracting a lot of interest because of their suppression effect on turbulence in magnetically confined fusion devices. The current naïve suggestion regarding this phenomenon holds that the stretching of modes in a shear flow brings about length scale reduction leading to the suppression of fluctuations. This argument for stability, however, ignores the fact that the available free energy associated with a shear flow may be a potent source for the destabilization of some other class of fluctuations (e.g., [18] [19] [20]). The Kelvin-Helmholtz (KH) instability, for instance, is a well-known example of an instability that feeds on the ambient flow-energy [23]. Linear wave patterns in Jupiter’s clouds with wavelengths strongly clustered around 300 km are commonly observed in the planet’s equatorial atmosphere (see [24]). In the work [25] proposed that the preferred wavelength is related to the thickness of an unstable shear layer within the clouds (see also [26]). If numerically analyze the linear stability of wavelike disturbances that have nonzero horizontal phase speeds in Jupiter’s atmosphere and find that, if the static stability in the shear layer is very low (but still nonnegative), a deep vertical shear layer like the one measured by the Galileo probe can generate the instabilities [27] The fastest growing waves grow exponentially within an hour, and their wavelengths match the observations. Close to zero values of static stability that permit the growth of instabilities are within the range of values measured by the Galileo probe in a hot spot and the model probes Jupiter’s equatorial atmosphere below the cloud deck and suggests that thick regions of wind shear and low static stability exist outside hot spots [28]. Note that all corresponding figures are relative to investigations of KHM and instabilities can be founded from the list of astrophysical and astronomical processes that arise from KHI problems (from Wikipedia). About KHI and ODE simultaneously, our ideas have the following forms.
Note that above unmentioned authors works arisen ordinary differential equations has been contained degenerated cases, but they do not obtained exactly solutions for its problems, because at the degenerated points to decomposing to the series impossible, since at these points the derivatives does not exists. In order to establish this situation all above considered problems which is analytical mathematical description of KHI beings to the so-called non-classical ordinary differential equations and its corresponding boundary conditions problems needed to prove solvability. After finding class of solvability we find exactly generalized solution, because generalized derivatives any order not dependent from existence of first order. It may be only applying non classical approaches methods. Therefore our method is new and equations also we used also belongs non-classical ordinary differential equations and having physical mean boundary layer, at the same time do not fit classes theory of ODEs. Therefore, our investigation basic objects is generalized theory than until exists theory of ODEs and its applications of physical problems. For this reason we introduce the following mathematical justification of this theory in subsection 1.2 and physical mean illustrated in subsection 1.3.
1.2. History of Non-Classical Equation which is Arises in KHI Problems
Equation of M. V. Keldysh (with strong degenerating at point
, elliptic equation) [8]:
(1)
Equation of M. V. Keldysh in model form (weak degenerating at point
, elliptic equation):
(2)
Equation of F. Tricomi in model form:
(3)
In 1985 year, at first, obtained by author M.A. Nurmammadov (see [4] [29] in page 2 noted the citation of this situation) new equation and four new well-posed boundary value problems, which model of M. V. Keldysh type(degenerating elliptic, degenerating hyperbolic, equations of mixed type, equation of mixed type changing time of directions) (having perpendicular degenerating two lines(two sonic lines) (having physical means of aerodynamic, gas dynamic, fluid dynamics and so on) [4]:
(4)
(or
)
(5)
(where l, a, b are numbers). Nonlinear case:
where,
. (6)
(7)
(8)
(9)
(
)
(where L, a, b are numbers).
Note that the linear ordinary differential equations of (7) for case
, and (8) for
, are degenerates which is corresponds to non-classical equations of second order (4), (9) (for equation of (9) account into the spaces variables
), respectively. At the same time that the nonlinear ordinary differential equations of (5) in case of
,
, is degenerate which is corresponds to non-classical equations of second order (5). We will deal with equations of (6), (7), (8). While, in the world many author considered: i) case
, ii)-case
(10)
Or in terminology Equation (4):
,
, but
, is number (11)
The mechanism was originally proposed by Kelvin and Helmholtz in the late nineteenth century to explain the source of energy of the Sun. By the mid-nineteenth century, conservation of energy conservation had been accepted, and one consequence of this law of physics is that the Sun must have some energy source to continue to shine. Because nuclear reactions were unknown, the main candidate for the source of solar energy was gravitational contraction problems almost equations being the following degenerating ordinary differential equations when
, also is
,
,
, at the points:
(12)
equation degenerate, which can be considered as a one dimensional alternative of a second order elliptic, mixed types equations degenerating (above written as (1), (4), (5), (6), (7), (8)) on the boundary of the domain (interval, plane suffuse, and so on). Finally, we note that in mathematically senses, this presented investigation equations are new ordinary differential equations and its physical means at first which is no considered with others authors and considered our ordinary differential equations include:
i) Presented difficult solvable of NODE is degenerating coefficient of first and second order derivatives;
ii) Nonlinearity of power type case—so called quasi-linear degenerating coefficient of first and second order derivatives;
iii) In given intervals coefficients of degenerating cases of NODEs are arbitrarily changes;
iv) Correctness of well-posed problems depends from coefficients higher and younger terms of derivatives;
v) Presented NODEs contains all possible cases classes of ordinary differential equations which arises KHI;
vi) Generalized solution at the same time is exact solution and common case can be used for all problems of KHI which being to ODEs;
vii) The main results of our investigations consist of finding generalized solutions to DODE or NODE and proving its solvability (i.e. existence and uniqueness of problems), after showing that our results are applicable to all problems of KHI which is being to the ODE;
viii) In spite of difficult problems that arise when the above aforementioned equation has singularities and degenerating cases simultaneously, at points and in this case, our results also remain valid, in the case of linear wave analysis of Kelvin-Helmholtz instability (KHI) at a tangential discontinuity (TD) of ideal magneto-hydro-dynamic (MHD) plasma, it can be attributed to the presented class, and in this case, as far as we know, solutions for eigen modes of instability KH in MHD plasma that satisfy suitable nonhomogeneous boundary conditions. Our results are establishing continuity at points of degenerating and singularity. Above listed i) - viii) are belongs to our objects investigations.
1.3. Illustration Hydrodynamic, Aerodynamic and Gas-Dynamic Physical Meaningful
In case of main non-classical Equations ((4), (6), (9)) (which is at first time obtained by author), also at first time in complex form transition of subsonic, sonic and supersonic for hydrodynamic, aerodynamic, gas-dynamic meaningful which for each equations using with the number Mache its corresponding physical process was presented (see, Figure 1 is for Equations (4), (5), (9), but Figure 2 for Equations (6), (12)) (also take nonhomogeneous B.V.C.):
Figure 1. Depicts physical process transitions in gas dynamics, aerodynamics, and hydrodynamics for Equations (4), (5), and (9).
Figure 2. Depicts physical process transitions in gas dynamics, aerodynamics, and hydrodynamics for Equations (6) and (12).
2. Main Results: Solvability of Boundary Value of Problems for General Degenerates Ordinary Equations which Is Arise in KHM and Instabilities
2.1. Statements Well-Posed Problem and Auxiliary Lemmas
Boundary value problem. To find solution of Equation (12), such that satisfy the following boundary conditions
(13)
By the symbol
, we denote a class of twice continuously differentiable functions in the closed domain
satisfying the boundary conditions (13), by
, Sobolev’s space (see [10] [11]) with weighted spaces obtained by the class
, which is closed by the norm:
(14)
Lemma 2.1. Let there exists
and
constants, such that (i)
,
; (ii)
,
, (iii)
.
Then there exists constants
are such that for any
the inequalities are holds true:
(15)
As well as known that, ODEs of second order which is obtain from the elliptical types equation second order as one dimensional analogical sense, then in this case is true the inequality of Frederic’s-Poincare, In order to apply non-classical approaches (method functional analyses) it is important prove it. Following as methods of work (e.g, [8] [9] [10] [11]) integrating by parts and account into boundary conditions, using inequality of Cauchy’s we can prove inequalities (15).
Lemma 2.2. (For ODE, Frederic’s-Poincare,). If the function
then the following inequality
is holds true.
Proof. It is sufficient to take
continuously differential function vanishing on
, and to take
, to be the interval
. Then the boundary conditions on
allow us to write
. By the Schwarz inequality we have
(16)
where m is a positive constant which depends on interval
, squaring both sides we obtain
. (17)
If we integrate both sides with repeat to x between –L and L, we multiply the constant on the right by a new constant which also depends on interval, as m(L), then we have
(18)
This completes the proof of Lemma 1.2, which has been proved for the first time for a non-classical ordinary differential Equation (12) (or degenerating ODE second order).
Thus, let us to introduce the following notations:
,
, and preliminary functions
,
(19)
Let
is the closed classes functions of
with the norm:
. Denote the new Hilbert space, which is completed by the norm
. Now, introducing the function (averaging)
(20)
where
(21)
and
chosen from the relations:
and 、
. Since,
,
for
and
;
,
,
; (22)
then for any
which is weighted space,
,
, additionally, may satisfy:
,
. For any function
setting the new functions:
,
, additionally scalar products in the space
the following:
(23)
Definition 2.1. The function
is said to be generalized solution of the problem (12), (13) if for any
and
, the following identity holds:
(24)
Lemma 2.3. If the function
is generalized solution of the problem (12), (13), then there exist the following finite limits:
, hence also is
, for
,
, hence also is
, for
. Proof. From the equality of (24) we have
(25)
Hence, for
we get:
(26)
Since, for
the following integrals are converges:
. By virtue of condition, when
the
and hence, we get
when
,
. Also, by virtue of condition, when
the
and hence, we get
when
. Indeed in cases of when,
,
the integrals
and
can be estimates. Therefore, analogically, in cases
,
the integrals are converges
, and
. (27)
From the lemma 2.1, 2.2, 2.3 it follows that for weakly degenerating case in equation of (12) the term of
is small term and
, keeps meaningful. But for strong degenerating case consider the following notations:
Let
,
and
,
. (28)
Describe the common solution in case of Equation (11) in the form
where
are arbitrarily constants. Then we can write special solution of nonhomogeneous Equation (11) in the form
(29)
Remark 2.1. If
,
is number, then
(30)
Lemma 2.4. Let fulfilled the Lemma 2.1 and Lemma 2.2 Then, for any
(31)
Proof. Let’s decompose (22) in the following
(32)
In equality of (24) the expression for the first bracket tends to zero and the second bracket, after estimation and applying the rule of L’Hopital, can be asserted by Lemma 2.4.
2.2. The Theorems of Existence and Uniqueness of Generalized Solution of Problem (12), (13)
Theorem 2.1. Let be satisfied Lemma 2.1, 2.2, 2.3 2.4 and i)
, ii)
, iii)
,
are satisfied. Then for any function from the class
there exists the generalized solution of boundary value
problem of (12), (13) from the space
.
Proof. Using definition of generalized solution
and
integrate the equality (24):
(33)
and account into results of lemma 2.1, 2.2, 2.3, 2.4, and with aid a prior estimation and boundary condition (13) obtain that generalized solution is
.
Theorem 2.2. Let be satisfied conditions of Theorem 2. Then, for any function of
. The generalized solution of problem (12), (13) from the class
is unique.
Proof. Assume that there are two distinct generalized solutions to problem (12). (13) from the class
. Take the function
using inequality of (14) we can write the following inequality
, hence we get
, and it means that
, i.e.
.
Remark 2.2. Smoothness the generalized solution of problem (12), (13) from the space
can be prove as it shown by author’s in the works (e.g., [13] [14] [15] [16]).
3. Some Class Equation Solutions Arise in the KHM and Instability
3.1. Some of Solution of Class Equations which is Presented by Its Coefficients
Let’s take for possible cases class equations for general non-classical ordinary equations (12), (in particular case of (11)) by selection the coefficients
which is corresponding arises in the different problems of KHM and instabilities
1)
or the number (
)
(34)
2)
,
or the number (
)
(35)
3)
, and others singularities of
,
(36)
4)
or the number (
)
(37)
5)
or the number (
)
(38)
I case:
,
, at points
,
and (dependence of problems sometimes
,
, in case of
).
Generalized solution of (8) is
(39)
where denoted by the
(40)
Hence, by the boundary conditions of (13) we have
(41)
(42)
II case. when
is number then
(43)
where denoted
,
,
,
(44)
III case when
(45)
Hence, by the boundary conditions of (13) we have
(46)
(47)
3.2. Applicable Cases (About 50 Problems Selected We Can Not Able Include Its)
Example 3.1. Comparison between barotropic and Kelvin-Helmholtz instability: the differential equation determining the horizontal structure of a perturbation in barotropic rotational flow has been stated in class. In the form which will be most evocative to us, it is:
,
,
at the point
is degenerate. In the analysis of barotropic instability, this is labeled the Rayleigh equation, and its similarity to the Taylor-Goldstein equation determining vertical structure of irrational flows is striking. In point of fact, they are actually the same equation with different antecedents.
Example 3.2.
, (48)
,
,
where
denotes by the Brunt-Vaisala frequency. The eigenvalue parameter of the problem is c. If the imaginary part of the wave speed wave speed C is positive, then the flow is unstable, and the small perturbation introduced to the system is amplified in time (the coefficients of this equation put in corresponding formulas of solution directly problem is solvable).If top and bottom be rigid walls, then from linearized covering equations of continuity case boundary condition is
, and when
then equation is no degenerating:
, and may be solved as problem of Sturm-Lovell (in simple form). Note that the argument is unchanged, if the top and bottom are at
. The Equation (48) constitute eigenvalue problem where
is the eigenvalue, if
instability occurs. But the generalized solution of (48) can be expressed by means of formulas (39)-(42), if substitute the corresponding coefficients and boundary conditions.
Example 3.3.
case corresponding of equation
(49)
In case of
,
,
,
,
are identically equations and at the points
is degenerating. Moreover, the coefficient
, at point
is also degenerates. In the case of finding the solution of a super linear ordinary differential equation with non-degenerating applications of form
with boundary conditions,
may be used as a theoretical sense of the work (see [30] and therein). Finally, our result for equation of (49) admissible to the Equations (6) and (12) are applicable and at the same time can be found generalized solution, if the coefficients of this equation we substitute in corresponding above solutions formulas and method of prove account into works (see [10] [11]), can asset the Equation (49) with corresponding boundary conditions is solvable.
Example 3.4. A KH instability on the planet Saturn, formed at the interaction of two bands of the planet; atmosphere (role of this example to illustrate boundary layer, see Figure 3).
Example 3.5. The Kelvin-Helmholtz can be seen in the bands of Jupiter (this example illustrates boundary layer, see Figure 4).
Example 3.6. The Kelvin-Helmholtz instability is commonly found in Earth’s
atmosphere and ocean (this example indicate possibility of KHI, see list of KHI in Wikipedia).
Example 3.7. This example illustrates that Kelvin-Helmholtz instability is possible for ocean mixing (see Wikipedia’s list of KHI).
4. Conclusions
In this work we obtain the following results:
Therefore, in the above information, we include that Kelvin-Helmholtz instabilities are visible in the atmospheres of planets and moons, for example, in the formation of the Earth or the Red Spot on Jupiter, and in the atmospheres of the Sun and other stars. For this reason, we investigate degenerate ordinary differential equations, which arise with very difficultly solvable problems of the Kelvin-Helmholtz mechanism and instabilities. In Figure 1 and Figure 2, we give illustrations of hydrodynamic, aerodynamic, and gas-dynamic physical meanings. The considered ODE of (12) is new (because it is degenerating) and corresponds in sense as physical aspects (see Figure 2) to a non-classical equation of mathematical physics (see the physical illustration in Figure 1) and is a one-dimensional analog of (4), as shown in the section introduction. Therefore it is called as non-classical ordinary differential equations (NODE). In Section 3, we proved the theorem of existence and uniqueness of common cases problems (8) and (13) in weighted spaces while also establishing generalized solutions. It was also demonstrated in this section that the results considered problems are applicable for such KHI families’ problems that arise in a variety of fields, including processes in astrophysics, the Sun, galaxies, stars, and planets. Regardless of the difficulties that arise when the above-mentioned equation has singularities and degenerates cases at the same time, at certain points, and in this case, our findings are also applicable to the linear wave analysis of Kelvin-Helmholtz instability (KHI) at a tangential discontinuity (TD) of ideal magneto-hydrodynamics (MHD) plasma, and the obtained results are validity. Note that additionally, the continuities of solutions are at the points singularity and degenerates cases established in the weighted spaces.