Next Article in Journal
Using Parameter Elimination to Solve Discrete Linear Chebyshev Approximation Problems
Next Article in Special Issue
Stabilization and the Design of Switching Laws of a Class of Switched Singularly Perturbed Systems via the Composite Control
Previous Article in Journal
Exponential and Hypoexponential Distributions: Some Characterizations
Previous Article in Special Issue
A Ginzburg–Landau Type Energy with Weight and with Convex Potential Near Zero
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Homogenization of a 2D Tidal Dynamics Equation

by
Giuseppe Cardone
1,*,
Aurelien Fouetio
2 and
Jean Louis Woukeng
3
1
Department of Engineering, Università del Sannio, Corso Garibaldi, 107, 84100 Benevento, Italy
2
Department of Mathematics, Higher Teacher Training College, University of Ngaoundere, P.O. Box 652 Bertoua, Cameroon
3
Department of Mathematics and Computer Science, University of Dschang, P.O. Box 67 Dschang, Cameroon
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(12), 2209; https://doi.org/10.3390/math8122209
Submission received: 19 November 2020 / Revised: 9 December 2020 / Accepted: 10 December 2020 / Published: 12 December 2020
(This article belongs to the Special Issue Asymptotic Analysis and Homogenization of PDEs)

Abstract

:
This work deals with the homogenization of two dimensions’ tidal equations. We study the asymptotic behavior of the sequence of the solutions using the sigma-convergence method. We establish the convergence of the sequence of solutions towards the solution of an equivalent problem of the same type.
MSC:
35B40; 46J10

1. Introduction

Ocean tides have been investigated by many authors, starting from [1,2]. The last decades have been marked by rapid progress in both theoretical and experimental studies of ocean tidal phenomena. Today, experimental and theoretical knowledge of ocean tides are used in order to address important problems in Oceanography, Atmospheric Sciences, Geophysics, as well as in Electronics and Telecommunications. It is important to point out that Laplace [3] was the first author to make the first major input in the theoretical formulation for tides of water on a rotating globe. Indeed, he formulated a system of partial differential equations that relate the horizontal flow to the surface height of the ocean. The existence and uniqueness of solutions of the deterministic tidal equation while using the classical compactness method have been proven in [2,4]. In our work, we consider the deterministic counterpart of a model of tidal dynamics that was studied by Manna et al. [5] and originally proposed by Marchuk and Kagan [2], where they considered the model of tidal dynamics derived by taking the shallow water model on a rotating sphere, with the latter being a slight generalization of the one considered earlier by Laplace.
Our objective is to carry out the homogenization of (2)–(5) under a suitable structural assumption on the coefficients of the operator that is involved in (2). These assumptions cover several physical behaviors, such as the periodicity, the almost periodicity, and much more. In order to achieve our goal, we shall use the concept of sigma-convergence [6], which is roughly a formulation of the well-known two-scale convergence method [7,8,9,10,11,12,13] in the context of algebras with mean value [6,14,15,16]. Therefore, our study falls within the framework of homogenization beyond the periodic setting, but including the periodic study as a special case.
The outline of the paper is as follows. The statement of the model problem, together with the derivation of appropriate uniform estimates, are the objectives of Section 2. Section 3 deals with the fundamentals of the sigma-convergence method. The homogenization process is performed in Section 4, while, in Section 5, we provide some applications of the main homogenization result.

2. Setting of the Problem and Uniform Estimates

2.1. Statement of the Problem

The tidal dynamics system that was developed by Manna et al. [5] for suitably normalized velocity u and tide height z reads as
u t + A ( u ) + B ( u ) + g z = f in Q = Ω × ( 0 , T ) z t + div ( h u ) = 0 in Q u = 0 on Ω × ( 0 , T ) u ( x , 0 ) = u 0 ( x ) and z ( x , 0 ) = z 0 ( x ) in Ω ,
where Ω is an open bounded subset, where A and B are defined by
A = α η η α Δ and B u = γ u + ω 0 u + ω 0 ,
α and η (the Coriolis parameter) being positive constants, ω 0 a given function and γ ( x ) = r / h ( x ) .
In our work, we ignore the Coriolis parameter ( η = 0 ), so that A ( u ) = α Δ u . However, instead of the Laplace operator, we rather consider a general linear elliptic operator of order 2 in divergence form, leading to the investigation of the limiting behavior (when 0 < ε 0 ) of the generalized sequence ( u ε , z ε ) ε of solution to the system (2)–(5), below
u ε t div A 0 ε u ε + B u ε + g z ε = f in Q
z ε t + div h u ε = 0 in Q
u ε = 0 on Ω × 0 , T
u ε x , 0 = u 0 x and z ε x , 0 = z 0 x in Ω ,
where Ω is a Lipschitz bounded domain of R 2 and T is a positive real number. Here, u ε and z ε stand, respectively, for the total transport 2-D vector (the vertical integral of the velocity) and the deviation of the free surface with respect to the ocean bottom. We have chosen the Dirichlet boundary condition in order to simplify the presentation. It is worth noticing that other boundary conditions can be considered, such as the Robin one: the difficulties may be only of technical types, with the method being the same.
In (2)–(5), ∇ (resp. div) is the gradient (resp. divergence) operator in Ω and the functions A 0 ε , h, u 0 , z 0 , and B are constrained, as follows:
(A1)
The function A 0 ε is defined by A 0 ε x = A 0 ( x , x / ε ) ( x Ω ), where A 0 C ( Ω ¯ , L ( R y 2 ) ) 2 × 2 is a symmetric matrix with
A 0 ( x , y ) ξ · ξ α ξ 2 for all ξ R 2 , x Ω ¯ and a . e . y R 2 ,
where α > 0 is a given constant independent of x, y and ξ .
(A2)
The operator B is defined on L 4 ( Ω ) 2 by B v = γ v + ω 0 v + ω 0 ( v L 4 Ω 2 ), where ω 0 L 2 ( 0 , T ; H 0 1 ( Ω ) 2 ) is a given function and γ ( x ) = r / h ( x ) (for a fixed real number r > 0 ), h being a continuously differentiable function satisfying
min x Ω h ( x ) = β > 0 , max x Ω h ( x ) = μ and max x Ω h ( x ) M ,
where M > 0 is a given constant, which is equal to zero at a constant ocean depth. The functions u 0 , z 0 and f are such that u 0 L 2 Ω 2 , z 0 L 2 Ω , f L 2 ( 0 , T ; H 1 ( Ω ) 2 ) , and g is the gravitational constant.
(A3)
We assume further that, for all x Ω ¯ , the matrix-function A 0 x , · has its entries in B A 2 ( R 2 ) , where A is an algebra with mean value on in R 2 , while B A 2 ( R 2 ) stands for the generalized Besicovitch space that is associated to A .
Remark 1.
The operator B continuously sends L 4 Ω 2 into L 2 Ω 2 with the following properties (see ([5], Lemma 3.3)): for u , v L 4 Ω 2 , we have
( B u B v , u v ) 0 ;
B u L 2 Ω 2 γ u L 4 Ω 2 ;
B u B v L 2 Ω 2 γ u L 4 Ω 2 + v L 4 Ω 2 u v L 4 ( Ω ) 2 .
The assumption (A3), which depends on the algebra with a mean value A , is crucial in the homogenization process. It shows how the microstructures are distributed in the medium Ω and, therefore, allows for us to pass to the limit.
Before dealing with the well-posedness of (2)–(5), we first need to define the concept of solutions that we will deal with.
Definition 1.
Let u 0 L 2 Ω 2 , z 0 L 2 Ω , f L 2 0 , T ; H 1 ( Ω ) 2 , ω 0 L 2 0 , T ; H 0 1 ( Ω ) 2 and 0 < T < . The couple u ε , z ε ε > 0 is a weak solution to the problem (2)–(5) if
u ε C 0 , T ; L 2 Ω 2 L 2 0 , T ; H 0 1 ( Ω ) 2 ; u ε t L 2 0 , T ; H 1 ( Ω ) 2 ; z ε L 0 , T ; L 2 Ω , z ε t L 2 0 , T ; L 2 Ω ;
and for all φ L 2 0 , T ; H 0 1 ( Ω ) 2 and ψ L 2 0 , T ; L 2 Ω , we have
0 T u ε t , φ d t + Q A 0 ε u ε · φ d x d t + Q B ( u ε ) φ d x d t + Q g z ε φ d x d t = 0 T ( f ( t ) , φ ( t ) ) d t
and
0 T z ε t , ψ d t + 0 T div ( h u ε ) , ψ d t = 0 .
In the above definition, ( · , · ) stands for the duality pairings between any Hilbert space X and its topological dual X . We also recall that the operator div A 0 ε u ε acts in a diagonal way, which is, for v = ( v 1 , v 2 ) H 0 1 ( Ω ) 2 , we have
div A 0 ε u ε , v = Ω A 0 ε u ε · v d x i = 1 2 Ω A 0 ε u ε i · v i d x
where u ε = ( u ε i ) 1 i 2 . This being so, the following existence and uniqueness result holds.
Theorem 1.
Assume that (A1)–(A2) are satisfied. Subsequently, there exists (for each ε > 0 ) a unique weak solution u ε , z ε to the problem (2)–(5) in the sense of Definition 1.
Proof. 
We note that, in the problem stated in [5], if we replace the Laplace operator by div ( A 0 ε u ε ) and we neglect therein the Coriolis parameter, then the proof follows exactly the lines of that of ([5], Propositions 3.6 and 3.7). □

2.2. A Priori Estimates

The following result will be useful in deriving the uniform estimates for u ε , z ε
Lemma 1
([5], Lemma 3.1). For any real-valued smooth function φ compactly supported in R 2 , we have
φ L 4 Ω 4 2 φ L 2 Ω 2 φ L 2 Ω 2 .
The following lemma provides us with the a priori estimates.
Lemma 2.
Under assumptions (A1)–(A2), the weak solution u ε , z ε of problem (2)–(5) in the sense that Definition 1 satisfies the following estimates
sup 0 t T u ε t L 2 Ω 2 C ;
0 T u ε t H 0 1 ( Ω ) 2 2 d t C ;
u ε t L 2 0 , T ; H 1 ( Ω ) 2 C ;
sup 0 t T z ε t L 2 Ω C ;
z ε t L 2 0 , T ; L 2 Ω C
where the positive constant C is independent of ε.
Proof. 
We first deal with Equation (2). In the variational form of (2), we choose the test function u ε ( t ) associated to (4) to obtain
1 2 d d t u ε t L 2 ( Ω ) 2 + A 0 ε u ε t , u ε t + B u ε t , u ε t + g z ε t , u ε t = f t , u ε t .
By the divergence theorem, we have
g z ε t , u ε t = g z ε t , div u ε t .
Applying Young’s inequality in the form
a b δ 2 a 2 + 1 2 δ b 2
to (19) (with δ = 2 g α ), we obtain
g z ε t , u ε t = g z ε t , div u ε t g 2 2 g α z ε t L 2 Ω 2 + α 2 g div u ε t L 2 Ω 2 g 2 2 g α z ε t L 2 Ω 2 + α 2 g u ε t H 0 1 Ω 2 2 .
In (7) if we take u = u ε and v = 0 to get B u ε t γ ω 0 2 , u ε t 0 , which yields
B u ε t , u ε t = B u ε t γ ω 0 2 , u ε t + γ ω 0 t 2 , u ε t γ ω 0 t 2 , u ε t r β ω 0 t L 4 Ω 2 2 u ε t L 2 Ω 2 r 2 β ω 0 t L 4 Ω 2 4 + u ε t L 2 Ω 2 2 .
Using again (20), but this time with δ = 1 , we get
f ( t ) , u ε t 1 2 f t L 2 Ω 2 2 + u ε t L 2 Ω 2 2 .
Putting together (6), (21)–(23), we derive, from (18), the following
d d t u ε t L 2 ( Ω ) 2 2 + 2 α u ε t H 0 1 Ω 2 2 f t L 2 Ω 2 + r β ω 0 t L 4 Ω 2 4 + u ε t L 2 Ω 2 2 + g 2 g α z ε t L 2 Ω 2 + α 2 g u ε t H 0 1 Ω 2 2 + u ε t L 2 Ω 2 2 = 1 + r β u ε t L 2 Ω 2 2 + 2 g 2 α z ε t L 2 Ω 2 + r β ω 0 t L 4 Ω 2 4 + α 2 u ε t H 0 1 Ω 2 2 + f t L 2 Ω 2 2 .
Integrating (24) with respect to t, we obtain
u ε t L 2 ( Ω ) 2 2 + 2 α 0 t u ε s H 0 1 Ω 2 2 d s 1 + r β 0 t u ε s L 2 Ω 2 2 d s + 2 g 2 α 0 t z ε s L 2 Ω 2 d s + r β 0 t ω 0 s L 4 Ω 2 4 d s + α 2 0 t u ε s H 0 1 Ω 2 2 d s + 0 t f s L 2 Ω 2 2 d s + u 0 t L 2 ( Ω ) 2 .
Next. dealing with (3), which we multiply by z ε ( t ) and then integrate the resulting equality over Ω , we get
1 2 d d t z ε t L 2 ( Ω ) 2 + div h u ε t , z ε t = 0 .
However,
div h u ε t , z ε t = h div u ε t , z ε t + u ε t · h , z ε t h div u ε t , z ε t + u ε t · h , z ε t h u ε t H 0 1 Ω 2 z ε t L 2 ( Ω ) + M u ε t L 2 Ω 2 z ε t L 2 ( Ω ) μ 2 α 2 μ u ε t H 0 1 Ω 2 2 + 2 μ α z ε t L 2 ( Ω ) 2 + M 2 u ε t L 2 Ω 2 2 + z ε t L 2 ( Ω ) 2 .
Taking (27) into account and integrating (26) in t gives
z ε t L 2 ( Ω ) 2 M 0 t u ε s L 2 Ω 2 2 d s + 2 μ 2 α + M 0 t z ε s L 2 Ω 2 d s + α 2 0 t u ε s H 0 1 Ω 2 2 d s + z 0 L 2 ( Ω ) .
Summing up inequalities (25) and (28), gives readily
u ε t L 2 ( Ω ) 2 2 + z ε t L 2 ( Ω ) 2 + α 0 t u ε s H 0 1 Ω 2 2 d s λ 1 0 t u ε s L 2 ( Ω ) 2 2 + z ε s L 2 ( Ω ) 2 d s + r β 0 t ω 0 s L 4 Ω 2 4 d s + λ 2 ,
where
λ 1 = max 1 + M + r β , 2 μ 2 α + M + 2 g 2 α
and
λ 2 = 0 T f s L 2 ( Ω ) 2 2 d s + u 0 L 2 ( Ω ) 2 2 + z 0 L 2 ( Ω ) 2 .
Now, appealing to inequality (12) (in Lemma 1) and owing to the fact that ω 0 L 2 ( 0 , T ; H 0 1 ( Ω ) 2 ) , we have
ω 0 ( s ) L 4 ( Ω ) 4 C ω 0 ( s ) L 4 ( Ω ) 2 ω 0 ( s ) L 2 ( Ω ) 2 C ω 0 ( s ) H 0 1 ( Ω ) 4 for a . e . s ( 0 , T ) ,
so that
ω 0 ( s ) L 2 ( 0 , T ; L 4 ( Ω ) 2 ) C ω 0 ( s ) L 2 ( 0 , T ; H 0 1 ( Ω ) 2 ) C .
We are, therefore, led to
u ε ( t ) L 2 ( Ω ) 2 2 + z ε ( t ) L 2 ( Ω ) 2 + α 0 t u ε s H 0 1 Ω 2 2 d s C + λ 1 0 t ( u ε ( s ) L 2 ( Ω ) 2 2 + z ε ( s ) L 2 ( Ω ) 2 ) d s .
Applying the Gronwall inequality leads to
sup 0 t T u ε t L 2 Ω 2 C , sup 0 t T z ε t L 2 Ω C , 0 T u ε t H 0 1 Ω 2 2 d t C .
From (10), we obtain, for all φ L 2 ( 0 , T ; H 0 1 Ω 2 ) ,
u ε t , φ C u ε L 2 ( 0 , T ; H 0 1 Ω 2 ) φ L 2 ( 0 , T ; H 0 1 Ω 2 ) + B u ε L 2 ( Q ) 2 φ L 2 ( Q ) 2 + C z ε L 2 ( Q ) φ L 2 ( 0 , T ; H 0 1 Ω 2 ) + f L 2 ( 0 , T ; H 1 Ω 2 ) φ L 2 ( 0 , T ; H 0 1 Ω 2 ) .
Next, while using the embedding H 0 1 ( Ω ) L 4 ( Ω ) , we have
B ( u ε ) L 2 ( Q ) C u ε L 2 ( 0 , T ; L 4 Ω 2 ) C u ε L 2 ( 0 , T ; H 0 1 Ω 2 ) .
Therefore, we infer, from (29), that
u ε t , φ C φ L 2 ( 0 , T ; H 0 1 ( Ω ) 2 ) ,
from which
u ε t L 2 ( 0 , T ; H 1 ( Ω ) 2 ) C .
We follow the same way of reasoning to see that
z ε t L 2 ( 0 , T ; L 2 ( Ω ) ) C .
This concludes the proof. □

3. Fundamentals of the Sigma-Convergence Method

Here, we are concerned with the main features of the sigma-convergence method, which we define. The reader may find more details in [14,15].
We first recall that by an algebra with mean value A on R d is meant any Banach algebra consisting of bounded uniformly contiuous functions on R d , satisfying
(i)
A contains the constants;
(ii)
u ( · + a ) A for any u A and a R d ; and,
(iii)
for any u A , the limit M ( u ) = lim R BR u(y)dy exists and is called the mean value of u.
In (iii), above, BR is the integral mean over the open ball BR centered at 0 and of radius R: BR = |BR|−1BR.
For obvious purposes, we define the generalized Besicovitch space B A p ( R d ) ( 1 p < ) associated to a given algebra with mean value A , as the completion with respect to the seminorm · p defined on A by u p = ( M ( u p ) ) 1 / p ( u A ). It is worth noticing that · p is well defined, since u p A for any p > 0 and u A . We may also define the Banach counterpart B A p ( R d ) of B A p ( R d ) by cutting with the kernel of the seminorm · p : B A p ( R d ) = B A p ( R d ) / N with N = { u B A p ( R d ) : u p = 0 } .
In the current work, we assume that all of the algebras with a mean value are ergodic (see [6,17] for the definition). We also need a further space, say B # A 1 , p ( R d ) , which is defined, as follows:
B # A 1 , p ( R d ) = { u W l o c 1 , p ( R d ) : y u B A p ( R d ) d and M ( y u ) = 0 } .
We identify two elements of B # A 1 , p ( R d ) by their gradients, which is, u = v in B # A 1 , p ( R d ) if y ( u v ) p = 0 . Equipped with the gradient norm u # , p = y u p , B # A 1 , p ( R d ) is a Banach space ([18], Theorem 3.12).
We are now able to define the concept of sigma-convergence.
Definition 2.
A sequence u ε ε > 0 L p ( Q ) ) 1 p < is said to:
(i) 
weakly Σ-converge in L p ( Q ) to u 0 L p ( Q ; B A p ( R d ) ) as if ε 0 , we have
Q u ε x , t f x , t , x ε d x d t Q M ( u 0 ( x , t , · ) f ( x , t , · ) ) d x d t
for every f L p ( Q ; A ) , 1 p + 1 p = 1 . We express this by writing u ε u 0 in L p ( Q ) -weak Σ;
(ii) 
strongly Σ-converge in L p ( Q ) to u 0 L p ( Q ; B A p ( R d ) ) if (30) holds and further
u ε L p Q u 0 L p Q ; B A p R d .
We express this by writing u ε u 0 in L p ( Q ) -strong Σ.
Remark 2.
(1) We can prove that the weak Σ-convergence in L p Q implies the weak convergence in L p Q . (2) The convergence (30) still holds true for f C ( Q ¯ ; B A p , R d ) , where B A p , R d = B A p R d L R d , 1 p + 1 p = 1 .
The following results are the main properties of sigma-convergence and they can be found in [6,14,17]. Before we can state them, we need to define what we call a fundamental sequence. By a fundamental sequence, we term any ordinary sequence ( ε n ) n 1 (denoted here, below, by E) of real numbers satisfying 0 < ε n 1 and ε n 0 when n .
(SC) 1
For 1 < p < , any sequence that is bounded in L p Q possesses a weakly Σ -convergent subsequence.
(SC) 2
Let ( u ε ) ε E L p ( 0 , T ; W 0 1 , p ( Ω ) ) ( 1 < p < ) be a bounded sequence in L p ( 0 , T ; W 0 1 , p ( Ω ) ) . Afterwards, up to a subsequence E from E, there exists a couple u 0 , u 1 with u 0 L p ( 0 , T ; W 0 1 , p ( Ω ) ) and u 1 L p ( Q ; B # A 1 , p ( R d ) ) , such that, as E ε 0 ,
u ε u 0 in L p Q - weak Σ
and
u ε x i u 0 x i + u 1 y i in L p ( Q ) - weak Σ , 1 i d .
(SC) 3
Let 1 < p , q < and r 1 be such that 1 r = 1 p + 1 q 1 . Assume that u ε ε > 0 L q ( Q ) is weakly Σ -convergent in L q ( Q ) to some u 0 L q ( Q ; B A q ( R d ) ) and v ε ε > 0 L p ( Q ) is strongly Σ -convergent in L p ( Q ) to some v 0 L p ( Q ; B A p ( R d ) ) . Subsequently, the sequence u ε v ε ε > 0 is weakly Σ -convergent in L r ( Q ) to u 0 v 0 .

4. Homogenization Result

4.1. Passage to the Limit

First, we set
V = u L 2 ( 0 , T ; H 0 1 ( Ω ) 2 ) : u = u t L 2 ( 0 , T ; H 1 ( Ω ) 2 ) ; H = H 1 ( 0 , T ; L 2 ( Ω ) ) .
The spaces V and H are Hilbert spaces with obvious norms. Moreover, the imbedding V L 2 ( 0 , T ; L 2 ( Ω ) 2 ) is compact.
Now, in view of a priori estimates in Lemma 2, the sequences ( u ε ) ε and ( z ε ) ε are bounded in V and in H , respectively. Thus, given a fundamental sequence E, there exist a subsequence E of E and a couple ( u 0 , z 0 ) V × H , such that, as E ε 0 ,
u ε u 0 in V - weak ;
u ε u 0 in L 2 ( 0 , T ; L 2 ( Ω ) 2 ) - strong ;
z ε z 0 in H - weak ;
Using the estimates (13)–(17), it follows that there exist a subsequence of E (not relabeled) and a function u 1 L 2 ( Q ; B # A 1 , 2 ( R 2 ) 2 ) , such that, as E ε 0 ,
u ε x i u 0 x i + u 1 y i in L 2 ( Q ) 2 - weak Σ , i = 1 , 2 .
It follows that ( u 0 , u 1 ) F 0 1 = V × L 2 ( Q ; B # A 1 , 2 ( R 2 ) 2 ) .
Now, for an element v = v 0 , v 1 F 0 1 , we set
D v = v 0 + y v 1 = ( D i v ) 1 i 2 where D i v = v 0 x i + v 1 y i , i = 1 , 2
with v 0 x i + v 1 y i = v 0 j x i + v 1 j y i 1 j 2 . The smooth counterpart of F 0 1 is defined by F 0 = C 0 ( Q ) 2 C 0 ( Q ; ( A / R ) 2 ) .
Proposition 1.
Let u = u 0 , u 1 F 0 1 and z 0 H . Afterwards, u and z 0 solve the following variational problem:
Q u 0 φ 0 t d x d t + Q M A 0 D u · D φ d x d t + Q B u 0 φ 0 d x d t + Q g z 0 φ 0 d x d t = 0 T f ( t ) , φ 0 ( t ) d t
Q z 0 ψ 0 t d x d t Q h u 0 · ψ 0 d x d t = 0
for all φ = φ 0 , φ 1 F 0 and ψ 0 C 0 ( Q ) .
Proof. 
Let φ = φ 0 , φ 1 and ψ 0 be, as above, and define
φ ε ( x , t ) = φ 0 ( x , t ) + φ 1 x , t , x ε for ( x , t ) Q .
Taking φ ε , ψ 0 as a test function in the variational form of (2)–(5), we obtain
Q u ε φ ε t d x d t + Q A 0 ε u ε · φ ε d x d t + Q B ( u ε ) φ ε d x d t + Q g z ε φ ε d x d t = 0 T f ( t ) , φ ε ( t ) d t
and
Q z ε ψ 0 t d x d t Q h u ε · ψ 0 d x d t = 0 .
While using the identities
φ ε t = φ 0 t + ε φ 1 t ε and φ ε = φ 0 + ( y φ 1 ) ε + ε ( φ 1 ) ε ,
we infer that, as ε 0 ,
φ ε t φ 0 t in L 2 0 , T ; H 1 ( Ω ) 2 - weak
φ ε φ 0 + y φ 1 in L 2 ( Q ) 2 × 2 - strong Σ
φ ε φ 0 in L 2 ( Q ) 2 - strong .
Let us consider each of the Equations (39) and (40) separately. We first consider (39) and, using the convergence results (34) and (41), we obtain
Q u ε φ ε t d x d t Q u 0 φ 0 t d x d t .
Because A 0 C ( Q ¯ ; B A 2 , R 2 2 × 2 ) , we use it as test function together with property (SC) 3 (recall that we have (36) and (42)) to obtain
Q A 0 ε u ε · φ ε d x d t Q M A 0 D u · D φ d x d t .
Let us show that
Q B ( u ε ) φ ε d x d t Q B ( u 0 ) φ 0 d x d t .
First, we have, from (34), that, up to a subsequence of E not relabeled, u ε u 0 a.e. in Q. Hence, from the continuity of B , we entail
B ( u ε ) B ( u 0 ) a . e . in Q .
we infer, from the boundedness of the sequence ( B ( u ε ) ) ε > 0 , that B ( u ε ) B ( u 0 ) in L 2 ( Q ) 2 -weak. Putting this together with (43), we obtain (46). We also easily obtain
0 T f ( t ) , φ ε ( t ) d t 0 T ( f ( t ) , φ 0 ( t ) ) d t .
Next, the convergence results (35) and (43) yield
Q g z ε φ ε d x d t Q g z 0 φ 0 d x d t .
As for Equation (40), we use the weak convergence (35) that is associated to (43) to obtain
Q z ε ψ 0 t d x d t Q z 0 ψ 0 t d x d t .
Concerning the second term in (40), we infer, from (34), that
Q h u ε · ψ 0 d x d t Q h u 0 · ψ 0 d x d t ,
thereby completing the proof of the proposition. □

4.2. Homogenized Problem

Here, we intend to derive the problem whose the couple u 0 , z 0 is solution. In order to achieve this, we first uncouple Equation (37), which is equivalent to the system consisting of (49) and (50), below:
Q u 0 φ 0 t d x d t + Q M A 0 D u · φ 0 d x d t + Q B u 0 φ 0 d x d t + Q g z 0 · φ 0 d x d t = 0 T f ( t ) , φ 0 ( t ) d t ;
Q M A 0 D u · y φ 1 d x d t = 0 .
Choosing in (50)
φ 1 ( x , t , y ) = θ ( x , t ) v ( y ) where θ C 0 Q , v ( A ) 2 ,
we obtain
M A 0 D u · v = 0 for all v ( A ) 2 .
Let us deal with (52). To this end, fix ξ R 2 × 2 and consider the corrector problem:
Find ß ξ C ( Ω ¯ ; B # A 1 , 2 ( R 2 ) 2 ) such that: div y [ A 0 ( x , · ) ( ξ + y ß ( ξ ) ) ] = 0 in R 2 .
Subsequently, in view of the properties of the matrix A 0 ( x , · ) , we infer, from [19,20], that (53) possesses a unique solution in C ( Ω ¯ ; B # A 1 , 2 ( R 2 ) 2 ) . Coming back to (53) and taking there ξ = u 0 ( x , t ) , testing the resulting equation with φ 1 as in (51), we get, by the uniqueness of the solution of (53), that u 1 ( x , t , y ) = ß ( u 0 ( x , t ) ( x , y ) . This shows that ß ( u 0 ) belongs to L 2 ( 0 , T ; C ( Ω ¯ ; B # A 1 , 2 ( R 2 ) 2 ) ) . Clearly, if χ j is the solution of (53) corresponding to ξ = ξ j = ( δ i j δ k ) 1 i , k 2 (that is all of the entries of ξ are zero, except the entry occupying the jth row and the th column, which is equal to 1), then
u 1 = j , = 1 2 u 0 x j χ j where u 0 = ( u 0 ) 1 2 .
We recall again that χ j depends on x, as it is the case for A 0 . In the variational form of (49), we insert the value of u 1 that was obtained in (54) to obtain the equation
u 0 t div ( A ^ 0 ( x ) u 0 ) + B u 0 + g z 0 = f in Q .
where A ^ 0 ( x ) = ( a ^ i j k ( x ) ) 1 i , j , k , 2 , a ^ i j k ( x ) = a hom ( χ j + P j , χ i k + P i k ) with P j = y j e ( e the th vector of the canonical basis of R 2 ) and
a hom ( u , v ) = i , j , k = 1 2 M a i j u k y j v k y i where A 0 = ( a i j ) 1 i , j 2 .
Additionally, Equation (38) is equivalent to
z 0 t + div ( h u 0 ) = 0 in Q .
Finally, putting together the Equations (55) and (56) associated to the boundary and initial conditions, we are led to the homogenized problem, viz.
u 0 t div ( A ^ 0 ( x ) u 0 ) + B u 0 + g z 0 = f in Q z 0 t + div ( h u 0 ) = 0 in Q u 0 = 0 on Ω × 0 , T u 0 ( x , 0 ) = u 0 ( x ) , z 0 ( x , 0 ) = z 0 ( x ) in Ω .
It can be easily shown that the matrix A ^ 0 of homogenized coefficients has entries in C ( Ω ¯ ) , and it is uniformly elliptic, so that, under the conditions (A1)–(A2), the problem (57) possesses a unique solution u 0 , z 0 with u 0 L 2 0 , T ; H 0 1 ( Ω ) 2 and z 0 L 2 0 , T ; L 2 ( Ω ) . Because the solution of (57) is unique, we infer that the whole sequence u ε , z ε converges in a suitable space towards u 0 , z 0 , as stated in the following result, which is the main result of the work.
Theorem 2.
Assume that (A1) to (A3) hold. For any ε > 0 , let ( u ε , z ε ) be the unique solution of problem (2) to (5). Subsequently, u ε u 0 in L 2 ( Q ) 2 -strong and z ε z 0 in L 2 ( Q ) -weak, where ( u 0 , z 0 ) is the unique solution of problem (57).
Proof. 
The proof is a consequence of the previous steps. □

5. Some Concrete Applications of Theorem 2

The homogenization of problem has been made possible under the fundamental assumption (A3). Some of the physical situations that lead to (A3) are listed below.
Problem 1 (Periodic Homogenization).
The homogenization of (2)–(5) holds under the periodicity assumption that the matrix-function A 0 x , · is periodic with period 1 in each coordinate, for any x Ω ¯ . In that case, we have A = C p e r ( Y ) , where Y = ( 0 , 1 ) 2 and C p e r ( Y ) is the algebra of continuous Y-periodic functions defined in R 2 . It is easy to see that B A 2 ( R 2 ) = L p e r 2 ( Y ) { u L l o c 2 ( R 2 ) : u is Y-periodic}, and the mean value expresses as M ( u ) = Y u ( y ) d y . Hence, the homogenized matrix is defined by A ^ 0 ( x ) = ( a ^ i j k ( x ) ) 1 i , j , k , 2 , a ^ i j k ( x ) = a hom ( χ j + P j , χ i k + P i k ) with P j = y j e ( e the th vector of the canonical basis of R 2 ) and
a hom ( u , v ) = i , j , k = 1 2 Y a i j u k y j v k y i d y w h e r e A 0 = ( a i j ) 1 i , j 2 .
where, here, χ j is the solution of the cell problem
χ j ( x , · ) H # 1 ( Y ) 2 : div y ( A 0 ( x , · ) ( ξ j + y χ j ( x , · ) ) ) = 0 i n Y
with H # 1 ( Y ) = { v H p e r 1 ( Y ) : Y v d y = 0 } , H p e r 1 ( Y ) = { v L p e r 2 ( Y ) : y v L p e r 2 ( Y ) 2 } and ξ j = ( δ i j δ k ) 1 i , k 2 .
Problem 2 (Almost periodic Homogenization).
We may consider the homogenization problem for (2)–(5) under the assumption that the coefficients of the matrix A 0 x , · are Besicovitch almost periodic functions [21]. In that case, hypothesis (A3) holds true, with A = AP ( R 2 ) , where AP R 2 is the algebra of Bohr almost periodic functions on R 2 [22]. The mean value of a function u AP R 2 is the unique constant that belongs to the close convex hull of the family of the translates u · + a a R 2 .
Problem 3 (Weakly almost periodic Homogenization).
We may solve the homogenization problem for (2)–(5) under the assumption: the function A 0 x , · is weakly almost periodic, which is, the matrix A 0 x , · has its entries in the algebra with mean value A = WAP R 2 (where WAP R 2 is the algebra of continuous weakly almost periodic functions on R 2 ; see, e.g., [23]).

Author Contributions

G.C., A.F. and J.L.W. contributed in writing the original draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The third author is very grateful to the Alexander von Humboldt Foundation for the Equipment’s Grant he received.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lamb, H. Hydrodynamics; Dover Publications: New York, NY, USA, 1932. [Google Scholar]
  2. Marchuk, G.I.; Kagan, B.A. Ocean Tides, Mathematical Models and Numerical Experiments; Pergamon Press: Elmsford, NY, USA, 1984. [Google Scholar]
  3. Laplace, P.S. Recherches sur quelques points du système du monde. Mem. Acad. Roy. Sci. Paris 1775, 88, 75–182. [Google Scholar]
  4. Ipatova, V.M. Solvability of a tide dynamics model in adjacent seas. Russ. J. Numer. Anal. Math. Mod. 2005, 20, 67–79. [Google Scholar] [CrossRef]
  5. Manna, U.; Menaldi, J.L.; Sritharan, S.S. Stochastic Analysis of Tidal Dynamics Equation. In Infinite Dimensional Stochastic Analysis: In Honor of Hui-Hsiung Kuo; World Scientific Publishers: Hackensack, NJ, USA, 2008. [Google Scholar]
  6. Sango, M.; Svanstedt, N.; Woukeng, J.L. Generalized Besicovitch spaces and application to deterministic homogenization. Nonlinear Anal. TMA 2011, 74, 351–379. [Google Scholar] [CrossRef] [Green Version]
  7. Bunoiu, R.; Cardone, G.; Kengne, R.; Woukeng, J.L. Homogenization of 2D Cahn-Hilliard-Navier-Stokes system. J. Elliptic Parabolic Equat. 2020, 6, 377–408. [Google Scholar] [CrossRef]
  8. Cardone, G.; Woukeng, J.L. Corrector problem in the deterministic homogenization of non- linear elliptic equations. Appl. Anal. 2019, 98, 118–135. [Google Scholar] [CrossRef]
  9. Cardone, G.; Perugia, C.; Timofte, C. Homogenization results for a coupled system of reaction-diffusion equations. Nonlinear Anal. 2019, 188, 236–264. [Google Scholar] [CrossRef] [Green Version]
  10. Cardone, G.; Panasenko, G.; Sirakov, Y. Asymptotic analysis and numerical modeling of mass transport in tubular structures. Math. Models Meth. Appl. Sc. 2010, 20, 1–25. [Google Scholar] [CrossRef] [Green Version]
  11. Borisov, D.; Cardone, G. Complete asymptotic expansions for the eigenvalues of the Dirichlet Laplacian in thin three-dimensional rods. Control Optim. Calc. Var. 2011, 17, 887–908. [Google Scholar] [CrossRef]
  12. Cardone, G.; Durante, T.; Nazarov, S.A. Water-waves modes trapped in a canal by a body with the rough surface. Z. Angew. Math. Mech. 2010, 90, 983–1004. [Google Scholar] [CrossRef] [Green Version]
  13. Cardone, G.; Minutolo, V.; Nazarov, S.A. Gaps in the essential spectrum of periodic elastic waveguides. Z. Angew. Math. Mech. 2009, 89, 729–741. [Google Scholar] [CrossRef]
  14. Woukeng, J.L. Homogenization in algebra with mean value. Banach J. Math. Anal. 2015, 9, 142–182. [Google Scholar] [CrossRef] [Green Version]
  15. Woukeng, J.L. Introverted algebras with mean value and applications. Nonlin. Anal. TMA 2014, 99, 190–215. [Google Scholar] [CrossRef] [Green Version]
  16. Zhikov, V.V.; Krivenko, E.V. Homogenization of singularly perturbed elliptic operators. Math. Notes Acad. Sc. USSR 1983, 33, 294–300. [Google Scholar] [CrossRef]
  17. Nguetseng, G.; Sango, M.; Woukeng, J.L. Reiterated ergodic algebras and applications. Commun. Math. Phys 2010, 300, 835–876. [Google Scholar] [CrossRef]
  18. Casado Diaz, J.; Gayte, I. The two-scale convergence method applied to generalized Besicovitch spaces. Proc. R. Soc. Lond. A 2002, 458, 2925–2946. [Google Scholar] [CrossRef]
  19. Jäger, W.; Tambue, A.; Woukeng, J.L. Approximation of homogenized coefficients in deterministic homogenization and convergence rates in the asymptotic almost periodic setting. arXiv 2019, arXiv:1906.11501. [Google Scholar]
  20. Shen, Z. Convergence rates and Hölder estimates in almost-periodic homogenization of elliptic systems. Anal. PDE 2015, 8, 1565–1601. [Google Scholar] [CrossRef]
  21. Besicovitch, A.V. Almost Periodic Functions; Dover Publications: Cambridge, UK, 1954. [Google Scholar]
  22. Bohr, H. Almost Periodic Functions; Chelsea: New York, NY, USA, 1947. [Google Scholar]
  23. Eberlein, W.F. Abstract ergodic theorems and weak almost periodic functions. Trans. Amer. Math. Soc. 1949, 67, 217–240. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cardone, G.; Fouetio, A.; Woukeng, J.L. Homogenization of a 2D Tidal Dynamics Equation. Mathematics 2020, 8, 2209. https://doi.org/10.3390/math8122209

AMA Style

Cardone G, Fouetio A, Woukeng JL. Homogenization of a 2D Tidal Dynamics Equation. Mathematics. 2020; 8(12):2209. https://doi.org/10.3390/math8122209

Chicago/Turabian Style

Cardone, Giuseppe, Aurelien Fouetio, and Jean Louis Woukeng. 2020. "Homogenization of a 2D Tidal Dynamics Equation" Mathematics 8, no. 12: 2209. https://doi.org/10.3390/math8122209

APA Style

Cardone, G., Fouetio, A., & Woukeng, J. L. (2020). Homogenization of a 2D Tidal Dynamics Equation. Mathematics, 8(12), 2209. https://doi.org/10.3390/math8122209

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop