Sturm-Liouville Boundary Value Problem for a Sea-Breeze Flow

We present recent results in study of a mathematical model of the Sea-Breeze flow, arising from a general model of the ‘morning glory’ phenomena. Based on analysis of the Dirichlet spectrum of a corresponding Sturm-Liouville problem and application of the Fredholm alternative, we establish conditions of existence/uniqueness of solutions to the given problem.


Introduction
Sea breeze is a atmospheric flow that develops due to a strong temperature contrast between the land and sea surfaces [17]. This flow is caused by the heating of the boundary layer over land, that results in the movement of low-level air from the sea to land (Sea-Breeze) with a return flow called the return current (see Fig. 1).
A typical Sea-Breeze flow is about 300 -1000m thick (with the strongest wind speed at 50-200m hight), whereas the thickness of the return current is about twice that big, leaving for the total vertical dimension of the Sea-Breeze circulation about 1 to 3 km [17]. The mass of air transported by the Sea-Breeze and by the return current are almost the same.
The Sea-Breeze flow strongly depends on the season, latitude where it blows and time of the day. For example, in tropical and subtropical coastal regions it is a regular phenomenon throughout the year, whereas at higher latitudes it is observed during the spring and summer periods (for more details about the Sea-Breeze flows we refer the reader to [7,15-17]).
In this paper we analyze a mathematical model of the Sea-Breeze flow in the Gulf of Carpentaria region (North Australia), that was recently derived from the equations characterizing the "morning glory" phenomena (for more information we refer to [6,8,11]); see Fig. 2.
A remarkable feature of the Sea-Breeze flow in this particular region, what also makes it so interesting to study, is its inland penetration for a distance of about 350km, that is significantly larger than the average distance, typical for these flows (see [3,15]). This is caused by several factors, under which are the intensity of heating, relative dryness of the land surface (except the summer monsoon periods) and a small value of the Coriolis parameter at these latitudes [7].
The mathematical model of the Sea-Breeze flow in the North Australian region is written in the form of the Dirichet boundary value problem, solvability of which we are analysing in this paper. Based on the Sturm-Liouville (S-L) theory [2,12,13] and the Fredholm alternative [5,14] we make conclusions about existence and uniqueness of solutions to the original problem.
The presented here results are new and serve as a natural extension to the analyses in [8]. It helps us to better understand the Sea-Breeze flows in a broader setting, than studied before. Moreover, these results

Mathematical Model of the Breeze-like Flow
In [8] Constantin, Johnson derived a non-dimensional model of the breeze-like flow, that in terms of a horizontal velocity V 0 (z, Φ) reads as: Here • z corresponds to the thickness of the flow; • ρ 0 (z) is the density function; • S, α, C and σ = 2(sin 2 α+C cos 2 α) (1−C) sin α cos α , are the characteristics of the flow in a specific region; • R e ≈ 10 5 is the Reynolds number; • m(z) is the viscosity function; • Φ is a parameter, corresponding to the direction of the flow propagation and • K(z, Φ) is the forcing term in the model. Together with the equation (1) one can consider the physically relevant boundary conditions: V 0 (0, Φ) = 0 -corresponding to the no-slip condition at the surface of the Earth, V 0 (z 0 , Φ) = 0 -standing for the level where the temperature inversion occurs.
In the Gulf of Carpentaria region, that is under our focus in this paper, one can determine the following physically relevant parameter values [8]: Then the differential equation (1) is simplified to with where functions ρ 0 (z) and K(z, Φ) are assumed to be continuous (with Φ being a parameter), while m(z) is continuously differentiable. Additionally, for every fixed Φ one can write the physically relevant constraints of the flow in the form of the homogeneous Dirichlet boundary conditions: In the following section we provide the solvability analysis to the problem (2), (4) using the spectral theory tools.

Sturm-Liouville Problem for a Sea-Breeze Flow Model
Let us associate with the BVP (2), (4) the correspondent S-L problem of the form: wherem(s) does not change its sign.
In [8] the authors showed that if is a self-adjoint unbounded linear operator acting in L 2 (0, 1) with domain the Sobolev space , then S has a discrete spectrum with simple eigenvalues λ 1 < λ 2 < . . . < λ n < . . . accumulating at infinity. Moreover, it was proved that λ 1 > β and thus, for β > 0 there are no zero eigenvalues and the BVP (2), (4) has a unique solution for all functions k 0 (s, Φ).
However, in some settings β admits a negative sign and thus, there exists such N that λ N = 0, which due to the Fredholm alternative leads to certain restrictions on the right hand-side of (2). This case was studied in [8] under assumption that the functionm(s) in the differential equation (2) is constant:m(s) ≡ m 0 . This leads to a sequence of eigenvalues λ k = β + m 0 π 2 k 2 with the correspondent eigenfunctions f n (s) = √ 2 sin(kπs), for k ≥ 1. One can easily see that for a particular value of β the sequence of eigenvalues λ k admits a zero value. This case occurs, for example, in the Calgary region where the sea-breeze-like flow is landlocked.
In this section we extend the spectral analysis presented in [8] by incorporating additional nonlinear profiles ofm(s) in the differential equation (2).

Equation with Constant Coefficients
Assume that function m(z) in the relation (3) is given by where as + b > 0, ∀s ∈ (0, 1), a = 0, b = 0 and a + b = 0. Thenm(s) = (as + b) 2 and the differential equation (2) can be written as Consider now the correspondent Sturm-Liouville BVP and let us find the Dirichlet spectrum for this problem. Note, that the differential equation ( 7) is the Legendre equation admitting its reduction to the equation with constant coefficients.
By introducing an Ansatz x = ln |as + b| we derive that Then the equation (7) can be written as From the characteristic equation corresponding to (8) we obtain that Below we provide a detailed analysis of the Sturm-Liouville problem (7), (4) for different values of λ. Case I. Assume, λ = β. Then μ 1,2 = 0, −1 and or, in terms of the s variable, Substitution of (10) into the boundary conditions (4) leads to or V 0 (s) = c 1 (as where μ 1,2 are defined by (9). The only possibility for V 0 (s) in (11) to satisfy the Dirichlet boundary constraints (4) is for c 1 = c 2 = 0. Thus, V 0 (s) ≡ 0 and there are no eigenvalues for a 2 − 4(λ − β) > 0.
From the last expression we get: and thus, , we conclude that the Dirichlet spectrum is defined by The eigenfunctions in this case are given by: From (16) follows that for all β ≥ 0 all eigenvalues λ n > 0 and accumulate at +∞. Thus, zero is not an eigenvalue of the S-L BVP (7), (4), and by the Fredholm alternative the original BVP (6), (4) has a unique solution [5]. However, if β < 0, then there might exist such N , for which λ N = 0. Then, the BVP (6), (4) has a unique solution iff 0n (s).

Legendre Equation
Consider the Sturm-Liouville problem (5), (4) and let us introduce an Ansatz: Then (5), (4) is written as where we assume that functionm(x)( is given by a relation: In this case the model equation (19) has the form and is coupled with the boundary constraints (20). Based on the analysis of the S-L BVPs for Legendre equations (see [2]), it is easy to check that Additionally, if β > 0 then all λ n > 0 and thus, due to the Fredholm alternative, the original BVP (2), (4) under substitutions (18), (21) has a unique solution. On the other hand, if β ≤ 0 (in particular for β = −4n(n + 1)), then there exists such N that λ N = 0, and we cannot guarantee uniqueness in this case. For a solution to be unique, function k 0 (z, Φ) should satisfy the orthogonality condition: < k 0 (z, Φ), P n (z) > 0, ∀n ∈ N 0 .
Note, that this analysis can be extended towards more complex forms of the viscosity functionm(s) that leads to S-L problems, for which their discrete spectrum can be explicitly found.