Next Article in Journal
A Generalized Weighted Monte Carlo Calibration Method for Derivative Pricing
Next Article in Special Issue
Automated Generation of EQ-Algebras through Genetic Algorithms
Previous Article in Journal
Method of Constructing a Nonlinear Approximating Scheme of a Complex Signal: Application Pattern Recognition
Previous Article in Special Issue
Real-Time Estimation of R0 for COVID-19 Spread
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multigrid Method for Optimal Control Problem Constrained by Stochastic Stokes Equations with Noise

by
Muhammad Munir Butt
Department of Mathematics and Statistics, College of Sciences, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia
Mathematics 2021, 9(7), 738; https://doi.org/10.3390/math9070738
Submission received: 4 March 2021 / Revised: 20 March 2021 / Accepted: 23 March 2021 / Published: 29 March 2021
(This article belongs to the Special Issue Numerical Analysis and Scientific Computing)

Abstract

:
Optimal control problems governed by stochastic partial differential equations have become an important field in applied mathematics. In this article, we investigate one such important optimization problem, that is, the stochastic Stokes control problem with forcing term perturbed by noise. A multigrid scheme with three-factor coarsening to solve the corresponding discretized control problem is presented. On staggered grids, a three-factor coarsening strategy helps in simplifying the inter-grid transfer operators and reduction in computation (CPU time). For smoothing, a distributive Gauss–Seidel scheme with a line search strategy is employed. To validate the proposed multigrid staggered grid framework, numerical results are presented with white noise at the end.

1. Introduction

Optimal control problems governed by partial differential equations (PDEs) have been studied extensively for the last three decades [1,2]. Particularly, control problem constrined by Stokes and/or Navier–Stokes equations play an important role in many real-world applications such as data assimilation for ocean flows, weather models, and flow problems. Recently, there has been a growing interest in mathematical analysis and simulations of stochastic PDEs (e.g., see [3,4,5,6,7,8]). The motivation of our work is the stochastic optimal control problem with PDEs as constraints (with random inputs) (see [8,9,10,11] and the references therein).
Multigrid methods are the well known efficient schemes to solve the elliptic boundary value problems [12,13,14]. However, to solve the saddle point problems such as the Stokes equations (after discretization with finite differences or finite element), the multigrid methods are more involved [15,16]. The multigrid methods have also been investigated for the flow problems [17,18,19,20,21,22,23,24]. Moreover, there has already been a lot of work on the multigrid solvers for the optimal control problems governed by PDEs in literature (e.g., see [25]).
This paper investigates a multigrid algorithm that solves the (finite difference) discretized optimality system corresponding to the control problem governed by stochastic Stokes equations (with force driven by additive white noise). In particular, we propose a full multigrid method with three-factor coarsening strategy to solve numerically the distributed optimal control problem constrained by stochastic Stokes equations (driven by the white noise perturbed forcing term). The scope of this paper is to develop an efficient multigrid algorithm without the need of any preconditioner [26], which results in fewer computations (CPU time) for the control problems with stochastic models, such as the one in [8]. For this, we take advantage of our previous works on multigrid methods for distributed optimal control problem governed by Stokes equations [27,28].
The objective of the proposed control formulation is to seek a state and a control that minimize the expectation of the distances between the approximate and desired state. In particular, the target state and the control are deterministic. The main difference between the present study and earlier research is the additive white noise in the forcing term that derives the stochastic Stokes equations. Therefore, we need to discretize the noise and modify the whole multigrid scheme for the optimization problem (stochastic model optimization) accordingly, which is explained in this article.
The paper is organized as follows. In next section. the control problem constrained by the stochastic Stokes equations is presented and first-order optimality conditions are discussed. On staggered grids, finite difference approximations are used to discretize the optimality system in Section 3. In Section 4, the proposed full multigrid method with smoothing and intergrid transfer operators (prolongation and restriction) is explained. In Section 5, numerical experiments for the control problem and stochastic Stokes equations with additive white noise are presented. In the end, the conclusions are given.

2. The Control Problem

We consider the following control problem in a convex polygon domain Ω R 2 with boundary Γ = Ω : Given a white noise W ˙ and a target function U , find an optimal stochastic state u H 1 ( Ω ) and a control f L 2 ( Ω ) while the deterministic cost functional is minimized, i.e.,
Min J ( u , f ) : = E 1 2 Ω | u U | 2 d x + α 2 Ω | f | 2 d x
subject to stochastic Stokes equations:
ν u + p = f + σ W ˙ ,   in Ω
· u = 0 ,   in Ω
u = 0 ,   on Γ .
We denote p L 2 ( Ω ) a pressure, f a control, and U L 2 ( Ω ) the desired velocity. Moreover, ν is a viscous constant, σ is a positive continuous function in Ω , α > 0 is the weight of the cost of control, W ˙ = ( W ˙ 1 , W ˙ 2 ) is the white noise, and
E [ W ˙ j ( x ) W ˙ j ( x ) ] = δ ( x x ) , x , x Ω .
We additionally assume that Ω p d x = 0 , that is, p satisfies the zero mean constraint. Furthermore, we denote u = ( u 1 , , u n ) T R n , | u | = u 1 2 , , u n 2 , unless otherwise stated.
We denote the inner product associated with Lebesque space L 2 ( Ω ) by ( · , · ) and
L 0 2 ( Ω ) = ϕ L 2 ( Ω ) : Ω ϕ d x = 0 , H 0 1 ( Ω ) = ψ H 1 ( Ω ) : ψ = 0 on Ω .
The weak solution ( u , p ) H 0 1 ( Ω ) × L 0 2 ( Ω ) of (2)–(4) is the solution of
ν a ( u , w ) + b ( w , p ) = ( f , w ) + ( σ W ˙ , w ) w H 0 1 ( Ω ) b ( u , ψ ) = 0 ψ L 0 2 ( Ω ) ,
where we denote a : H 0 1 ( Ω ) × H 0 1 ( Ω ) R and b : H 0 1 ( Ω ) × L 0 2 ( Ω ) as the bilinear forms, i.e.,
a ( u , w ) = j = 1 2 Ω u j · w j b ( w , p ) = Ω p · w
The existence and uniqueness of the control problem can easily be proven by following the standard arguments (see [1,29]). We define the set of admissible controls
F a d = { ( u , p , f ) H 0 1 ( Ω ) × L 0 2 ( Ω ) × L 2 ( Ω ) : J ^ ( f ) = J ( u , p , f ) < ; ( u , p , f ) satisfies ( 2 ) ( 3 ) }
The optimal solution ( u * , p * , f * ) F is such that if there exists a ε > 0 such that
J ( u * , p * , f * ) J ( u , p , f ) ( u , p , f ) F a d
and
u * u 2 + p * p 1 + f * f 2 ε
Now, we can formulate the optimal control problem that seeks ( u , p , f ) that minimize the cost functional (1) subject to (2)–(3), i.e.,
min J ( u , p , f )
as a constrained optimization problem with (2)–(3) in Hilbert space such that ( u , p , f ) F a d . In the following, we present the first-order optimality conditions for completeness:
Theorem 1.
The optimization problem (5) has a unique solution ( u * , p * , f * ) F a d .
Proof. 
Since F a d is not empty, for example, ( u , p , 0 ) F a d .
Let f n be a minimizing sequence and set the linear dependency of control to state variables as u n = u ( f n ) and p n = p ( u n ) . Moreover,
lim n J ( u n , p n , f n ) = inf J ( u , p , f )
By definition of F a d , the sequence f n is uniformly bounded in L 2 ( Ω ) . From the weak solutions to state system, we conclude that the sequence ( u n , p n ) is uniformly bounded in L 2 ( Ω ) × H 0 1 ( Ω ) almost surely. Thus, we have a subsequence ( u l , p l , f l ) such that
f l f * , p l p * , u l u *
for some ( u * , p * , f * ) L 2 ( Ω ) × H 0 1 ( Ω ) × L 2 ( Ω ) .
Next, by the weak lower semi-continuity of the cost functional J, we have
inf J ( u n , p n , f n ) = lim n J ( u n , p n , f n ) = J ( u * , p * , f * )
Hence, we see that the optimal solution (in the set F a d ) exists. Next, to prove that this optimal solution is unique, we use the fact that the constraints are linear and the cost functional is convex. This completes the proof. □
Next, the solution to the control problem (1)–(4) is characterized by the Lagrangian [1,30]. That is, the first-order optimality conditions for the (stochastic Stokes equations) optimization problem is obtained by taking the Gâteaux (Fréchet) derivative of the following Lagrangian functional
L ( u , f , p , λ , q ) = J ( u , f ) + ( Δ u + p f σ W ˙ , λ ) ( · u , q ) ,
where we denote the Lagrange multipliers λ H 0 1 ( Ω ) , q L 0 2 ( Ω ) . Thus, we have the optimality system given by [8]:
ν Δ u + p = f + σ W ˙ in Ω , · u = 0 in Ω ( state   system ) u = 0 on Γ ,
Δ λ + q = U u in Ω · λ = 0 in Ω ( adjoint   system ) λ = 0 on Γ ,
α f E [ λ ] = 0 in Ω , ( optimality   conditions )
where · u = div u , p = grad p , u = ( u , v ) , λ = ( λ , μ ) , etc.
From the optimality conditions, we have that f = E [ λ ] α . Therefore, we can write the above optimality system into a reduced optimality system as follows, i.e., we find ( u , p , λ , q ) while satisfying the following system:
Δ u + p = E [ λ ] α + W ˙ in Ω , · u = 0 in Ω , u = 0 on Γ , Δ λ + q = U u in Ω , · λ = 0 in Ω , λ = 0 on Γ , α f E [ λ ] = 0 in Ω .

3. Discretization

In the following, we discretize the optimality system using finite differences on staggered grids. For this, we first discretize the white noise W ˙ . Let { R h } be a family of rectangles Ω ¯ , where h is the mesh size. Define
ξ R j = 1 | R | R 1 d W j , j = 1 , 2
for each rectangle R R h , and | R | is the area of R. Here and in what follows, ξ R R j denotes an independent distributed normal random variables family, i.e., with the mean 0 and variance 1. The piecewise constant approximation to W ˙ j ( x ) is defined as follows
W ˙ h j ( x ) = R R h 1 | R | ξ R j χ R ( x ) ,
where χ R is the characteristic function of R. Moreover, W ˙ = ( W ˙ h 1 , W ˙ h 2 ) ( L 2 ( Ω ) ) 2 but W ˙ h is unbounded as h 0 (see Lemma 1.2 [8]).
Next, we consider the approximation problem (2)–(4) with the discrete white noise forcing term:
ν u h + p h = f h + σ W ˙ h ,   in Ω
· u h = 0 ,   in Ω
u h = 0 ,   on Γ .
Then, by Choi and Lee [8], we have the following estimates (for u h and p h , where h is the mesh size):
Proposition 1.
E u h 2 2 + p h 1 2 C h 2
and for errors u u h and p p h :
Proposition 2.
E u u h 2 + p p h 1 2 C | l n h | h 2 .
Next, we consider a sequence of grids { Ω h } h > 0 for a rectangular domain Ω given by
Ω h = { x R 2 : x i = i h x , y j = j h y , i , j Z } Ω ,
Here and in what follows, we take the mesh size ( h x = h y = h ) and choose h in such a way that the boundaries of Ω lie on the grid lines.
Note that, on the staggered grids, the state (respectively, adjoint) variables are placed on cell edges (vertical or horizontal) and on cell centers. We denote these sets of grid points with Ω h s , s { e v , e h , c } . For example, Ω h e v denotes the grid points defined on center of cell edge-vertical. Moreover, the discrete L 2 -scalar product is given by
( u h , v h ) L h 2 ( Ω h s ) = h x h y x Ω h s u h ( x ) v h ( x ) ,
and the associated norm is given by u h L h 2 ( Ω h s ) = ( u h , u h ) L h 2 ( Ω h s ) 1 / 2 . Let L h 2 ( Ω h s ) be the grid of grid functions u h defined on Ω h s . Moreover, we denote U h , V h and P h the spaces for the state (respectively, adjoint) grid functions variables u h , v h , and p h , respectively.
In the following, the discretization of stochastic Stokes equations (state equation) is presented. Note that, on staggered grids, the variables u , v , p are defined on Ω h e v , Ω h e h , and Ω h c , respectively (see Figure 1). Thus, by finite differences, we have
Δ h u j h + j h p h = f j h + σ W ˙ h j , at   centers   of   j - faces j = 1 2 j h u j h = 0 , at   cell   centers
where Δ h denotes the discrete Laplace operator, i.e., the five-point approximation. We approximate the first-order partial derivative by (second-order) central differences. For example, to discretize the pressure variable p at the cell center ( h x / 2 , h y / 2 ) , we have x h and y h as follows:
x h p ( x , y ) = u ( x + h x / 2 , y + h y / 2 ) u ( x h x / 2 , y + h y / 2 ) h x , y h p ( x , y ) = u ( x + h x / 2 , y + h y / 2 ) u ( x + h x / 2 , y h y / 2 ) h y .
Note that, for a point near a boundary, Δ h u j h ( x ) may involve an exterior value, and, to find these exterior values, quadratic extrapolation is used (cf. [17,27]).
Next, to have a more insight about the discretization using finite differences on staggered grid, we proceed as follows: set x i = ( i 1 ) h and y j = ( j 1 ) h . By p i + 1 / 2 , j + 1 / 2 , we mean the discrete counterpart to p ( x i + h / 2 , y j + h / 2 ) . and we consider the grid indices in lexicographic order. In this way, we come up with the state system given by
Δ h u h + p i 1 / 2 , j + 1 / 2 p i 3 / 2 , j + 1 / 2 h x = f i , j + 1 / 2 + σ W h 1 , on Ω h e v
Δ h v h + p i + 1 / 2 , j 1 / 2 p i + 1 / 2 , j 3 / 2 h y = g i + 1 / 2 , j + σ W h 2 , on Ω h e h
u i + 1 , j + 1 / 2 u i , j + 1 / 2 h x + v i + 1 / 2 , j + 1 v i + 1 / 2 , j h y = 0 , on Ω h c
where (15) and (16) are defined on Ω h e v and Ω h e h , respectively. Moreover, the continuity Equation (17) is defined on Ω h c .
The optimality conditions α f λ = 0 and α g μ = 0 illuminate that, for the adjoint equation, the variables u and λ , f are defined on Ω h e v and the variables v, g, and μ are defined on Ω h e h . Moreover, in the continuity equation for state (respectively, adjoint), the variable p and q appears on mesh cell centers Ω h c . Thus, in this way, we obtain the following discretized adjoint equation:
Δ h λ h + q i 1 / 2 , j + 1 / 2 q i 3 / 2 , j + 1 / 2 h x = ( U u ) i , j + 1 / 2 , on Ω h e v
Δ h μ h + q i + 1 / 2 , j 1 / 2 q i + 1 / 2 , j 3 / 2 h y = ( V v ) i + 1 / 2 , j , on Ω h e h
λ i + 1 , j + 1 / 2 λ i , j + 1 / 2 h x + μ i + 1 / 2 , j + 1 μ i + 1 / 2 , j h y = 0 on Ω h c .
Here, we remark that, with our approach, a direct coupling among the state, adjoint, and control variables has e been implemented with the help of optimality conditions (on staggered grids) and thus requires no additional computations (interpolation).
Summarizing, the discretize optimality system consists of Equations (15)–(20) with the optimality conditions given by
α f i , j + 1 / 2 λ i , j + 1 / 2 = 0 , on Ω h e v α g i + 1 / 2 , j μ i + 1 / 2 , j = 0 , on Ω h e h
with boundary conditions given by
u i , j + 1 / 2 = 0 , v i + 1 / 2 , j = 0 , λ i , j + 1 / 2 = 0 , μ i + 1 / 2 , j = 0 .

4. Multigrid

The multigrid framework to solve the optimality system is presented in this section. In the development of multigrid solver, we face some difficulties due to the coupled nature of the state and adjoint equations on the staggered grid structure.
It is a well known fact that multigrid methods involves discretization grids. For example, we obtain such grids by halving the mesh size starting from the coarsest grid [13,14]. In the context of staggered grid multigrid approach and in our previous works [27,28,31], multigrid schemes for the optimal control problems constrained by Stokes equations and/or PDEs, we note that, by tripling the mesh size, starting from the coarsest staggered grid, gives a nested hierarchy of grid that helps in simplifying the intergrid transfer operators. The details are as follows:
In the following, we define a sequence of grids Ω k of mesh size h k = h 1 / 3 ( k 1 ) , k = 1 , , L , where k = L is the finest level and h 1 = 1 / 2 is the mesh sizes of the coarsest grid. We denote the index k to define the functions (operators) on Ω k . Moreover, with three-factor coarsening, setting a variable X I J k 1 on the coarse grid Ω k 1 has the same location as the variable X i j k on the fine grid Ω k . Figure 1 presents an example.
  • u I , J + 1 / 2 k 1 u i j k , for i = 3 I 2 , j = 3 J 1 .
  • v I + 1 / 2 , J k 1 v i j k , for i = 3 I 1 , j = 3 J 2 .
  • p I + 1 / 2 , J + 1 / 2 k 1 p i j k , for i = 3 I 1 , j = 3 J 1 .

4.1. Smoothing Scheme

In this section, we present the smoothing scheme to relax the state (respectively, adjoint) momentum equations and the state (respectively, adjoint) continuity equation by a distributive relaxation [17]. A line search is used to update the control variable.
We denote ( u h , v h , p h , λ h , μ h , q h , f h , g h ) as the discrete counterpart to the numerical approximation to our optimization problem. Then, we update this approximation step by step iteratively, which is explained as follows. First, we update the control variable by the gradient update step (line search) with gradient of the so-called reduced cost functional, i.e, f J ^ ( f h , g h ) = α f h λ h , g J ^ ( f h , g h ) = α g h μ h as follows
f h f h t f J ^ , g h g h t g J ^ ,
where the step size for the control update is t ( 0 , 1 ] .
Then, we consider the sate equation and let ( u h , v h , p h ) be its approximation. We smooth the residuals of the state system (15)–(16) by the classical Gauss–Seidel (pointwise) relaxation. Next, we use distributive relaxation (cf. [17,27]) to relax the variable p h corresponding to the continuity Equation (17).
Next, we relax the adjoint equation as done for the discretized state equation. That is, the residuals of the adjoint system (18)–(19) are smoothed by the classical pointwise Gauss–Seidel scheme at all the interior points. Moreover, the (adjoint) continuity Equation (20) is relaxed by a distributive relaxation (cf. [27,28]).
In the following, we discuss the intergrid transfer operators for the proposed multigrid framework. We use bilinear interpolation to interpolate the (state and adjoint) variables and residual functions. For example, for the variable u, consider the space U k of u k : Ω k e v R , k = 1 , , L . Between the two grids, the fine grid Ω k and the coarser grid Ω k 1 , we define a prolongation operator, i.e., I k 1 k : U k 1 U k . Here, we assume that (for the discretization) this prolongation operator is consistent with the bilinear finite elements on each rectangular partition.
We use a straight injection operator as our restriction operator in the proposed multigrid solver. Note that, in the three-factor coarsening strategy, the location or position of the state (respectively, adjoint) variables on the coarse-grid and the fine-grid are the same (see Figure 1). Therefore, we use straight injection operator I k k 1 : U k U k 1 to transfer the residuals (solution) functions from fine to coarse grids. We also remark here that one can use the other operators such as half weighting or full weighting but this may increase the computations without any gain (efficiency). We employ straight injection because it is the advantage of the three-factor coarsening strategy on staggered grid and therefore a natural choice in this case.

4.2. Algorithms

Here, we present the algorithms for a full multigrid method (FMG) (see [14]). Since we know that the FMG method improves the computational complexity of a full approximation scheme (FAS), we use the full multigrid method to solve the discretized optimality system and hence the optimization problem.
To develop a multigrid algorithm, we consider the optimality system (15)–(22) corresponding to the stochastic Stokes control problem, at the discretization level k. Let the unknown variable be X k = ( u k , v k , p k , λ k , μ k , q k , f k , g k ) . Then, we can write the optimality system given by
A k ( X k ) = F k .
Let X k ( l ) = S k ( X k ( l 1 ) , F k ) be the smoothing result as given in Section 4.1. Then, starting with X k ( 0 ) , we obtain the approximate solution, e.g. X ˜ k = X k ( m 1 ) , after m 1 -times (iterations).
For completeness, we present the full approximation scheme (FAS) in Algorithm 1:
Algorithm 1: FAS ( m 1 , m 2 ) for solving A k ( X k ) = F k .
  • If k = 1 solve A k ( X k ) = F k exactly.
  • Pre-smoothing steps on the fine grid: X k ( l ) = S k ( X k ( l 1 ) , F k ) , l = 1 , , m 1 ;
  • Computation of the residual: r k = F k A k ( X k ( m 1 ) ) ;
  • Restriction of the residual: r k 1 = I k k 1 r k ;
  • Set X k 1 = I k k 1 X k ( m 1 ) ;
  • Set F k 1 = r k 1 + A k 1 ( X k 1 )
  • Call m times MG ( m 1 , m 2 ) to solve A k 1 ( X k 1 ) = F k 1 ;
  • Coarse-grid correction: X k ( m 1 + 1 ) = X k ( m 1 ) + I k 1 k ( X k 1 I k k 1 X k ( m 1 ) ) ;
  • Post-smoothing steps on the fine grid: X k ( l ) = S k ( X k ( l 1 ) , F k ) , l = m 1 + 2 , , m 1 + m 2 + 1 ;
In this algorithm, we can perform m two-grid iterations at each working level. For m = 1 , we have a V-cycle, and, for m = 2 , we have a W-cycle [14].
As is known, with a better initial guess, we can use fewer iterations [14]. Moreover, the FAS is implemented as a subroutine that improves the nonlinearity issues (if required in the problem). Therefore, a full multigrid (FMG) method (see Figure 2) is employed in this article and in Algorithm 2 form it is given as follows:
Algorithm 2: FMG for solving A L ( X L ) = F L .
  • For l = K < L , set initial approximation u l .
  • If l < L , then interpolate to the next finer working level: X ˜ l + 1 = I l l + 1 X l ;
  • Apply FAS to solve A l + 1 ( X l + 1 ) = F l + 1 , starting with X ˜ l + 1 ;
  • Set l : = l + 1 ; If l < L go to step 2; else stop.

5. Numerical Experiments

We considered the following optimization problem constrained by the stochastic Stokes equations to show the numerical results:
Min J ( u , f ) = 1 2 Ω | u U | 2 d x + α 2 Ω | f | 2 d x , subject   to ν u + p = f + σ W ˙ , in Ω · u = 0 , in Ω u = 0 , on Γ .
We set σ = 0.1 and ν = 1 and solved the control problem on various meshes and regularization parameter α > 0 , in the numerical experiments. We took a rectangular domain Ω = [ 0 , 1 ] 2 with uniform mesh size and used MATLAB(2020b) to develop a code for solving the control problem numerically, on a laptop i7 with 1.3 GHz, 16 GB RAM.
For our control problem, we do not need an exact solution. Therefore, we took the desired velocity [8] U = ( U ( x , y ) , V ( x , y ) ) given by
U ( x , y ) = 10 d d y ψ ( x ) ψ ( y ) and V ( x , y ) = 10 d d x ψ ( x ) ψ ( y )
where
ψ ( z ) = ( 1 z ) 2 ( 1 c o s ( 0.8 π z ) ) .
In Figure 3, a vector field of random velocities with the white noise is presented. Moreover, we depict the target velocity in vector form in Figure 3.
For the numerical solution of this control problem, we used the proposed FMG Algorithm 2. A W-cycle with 4 p r e and 4 p o s t smoothing steps was employed, that is, we took m 1 = m 2 = 4 in Algorithm 2). One can also use the V c y c l e s , but we did not find any improvement during our investigation. Furthermore, we chose t = 1 as the step length in the control update step. For the stopping criteria, we used the discrete L 2 -norm of all the variables ( u h , v h , p h , λ h , μ h , q h ) , that is, we stopped the iterations when
max r j h L 2 ( Ω h ) < 10 6 .
As expected, we have a better tracking for the desired (target) velocity field as α 0 (see Table 1). This fact is illustrated in Figure 4, that is, as α 0 the optimal velocities match the desired vector filed (see Table 2, Figure 3).

6. Conclusions

A multigrid scheme to solve the stochastic Stokes control problem with additive white noise is presented. On staggered grids, a three-factor coarsening strategy is employed, which results in nested hierarchy of staggered-grids. To solve the discretized optimality system with finite differences and noise, a line search strategy to update the control and a distributive relaxation scheme to relax the state (respectively, adjoint) variables is employed. Moreover, to discretize the additive white noise, we simply consider the Gaussian noise at the discretized level in our numerical investigation. Numerical results illuminate the efficiency and the ability of the proposed multigrid solver to track the desired target (velocity vector field) while the Stokes equations were driven by force term with additive white noise.
In comparison to the recent work [8] on the numerical solution of stochastic Stokes control problems with additive white noise, we note that our multigrid solver provides better iteration counts due to three-factor coarsening strategy, i.e., fewer levels, less computational cost, and less CPU time.
The proposed multigrid method can be extended to the three-dimensional case, i.e., stochastic Stokes control problem with noise in 3D staggered-grid multigrid framework. Moreover, one can investigate the Navier–Stokes control problems with noise, which is a natural extension of our present work.

Funding

This research was funded by Deanship of Scientific Research, King Fahd University of Petroleum and Minerals (KFUPM), Project No. SR191005.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author thanks the anonymous reviewers for their useful comments and suggestions that helped to improve the quality of the paper.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PDEsPartial differential equations
FASFull approximation scheme
FMGFull multigrid

References

  1. Lions, J.L. Optimal Control of Systems Governed by Partial Differential Equations; Springer: Berlin, Germany, 1971. [Google Scholar]
  2. Tröltzsch, F. Optimal Control of Partial Differential Equations: Theory, Methods and Applications; AMS: Providence, RI, USA, 2010. [Google Scholar]
  3. Babuska, I.; Chatzipantelidis, P. On solving elliptic stochastic partial differential equations. Comput. Methods Appl. Mech. Eng. 2002, 191, 4093–4122. [Google Scholar] [CrossRef]
  4. Babuska, I.; Tempone, R.; Zouraris, G.E. Galerkin finite element approximations of stochastic elliptic partial differential equations. SIAM J. Numer. Anal. 2004, 42, 800–825. [Google Scholar] [CrossRef]
  5. Babuska, I.; Tempone, R.; Zouraris, G.E. Solving elliptic boundary value problems with uncertain coefficients by the finite element method: The stochastic formulation. Comput. Methods Appl. Mech. Eng. 2005, 194, 1251–1294. [Google Scholar] [CrossRef]
  6. Cao, Y.; Chen, Z.; Gunzburger, M.D. Error analysis of finite element approximations of the stochastic Stokes equations. Adv. Comput. Math. 2010, 33, 215–230. [Google Scholar] [CrossRef]
  7. Chen, P.; Quarteroni, A.; Rozza, G. Multilevel and weighted reduced basis method for stochastic optimal control problems constrained by stokes equations. Numer. Math. 2016, 133, 67–102. [Google Scholar] [CrossRef]
  8. Choi, Y.; Lee, H.-C. Error analysis of finite element approximations of the optimal control problem for stochastic Stokes equations with additive white noise. Appl. Numer. Math. 2018, 133, 144–160. [Google Scholar] [CrossRef]
  9. Gunzburger, M.; Lee, H.-C.; Lee, J. Error estimates of stochastic optimal neumann boundary control problems. SIAM J. Numer. Anal. 2011, 49, 1532–1552. [Google Scholar] [CrossRef]
  10. Hou, L.; Lee, J.; Manouzi, H. Finite element approximations of stochastic optimal control problems constrained by stochastic elliptic PDEs. J. Math. Anal. Appl. 2011, 384, 87–103. [Google Scholar] [CrossRef] [Green Version]
  11. Lee, H.-C.; Gunzburger, M. Comparison of approaches for random PDEoptimization problems based on different matching functionals. Comput. Math. Appl. 2017, 73, 1657–1672. [Google Scholar] [CrossRef]
  12. Brandt, A. Multi-level adaptive solutions to boundary-value problems. Math. Comp. 1977, 31, 333–390. [Google Scholar] [CrossRef]
  13. Hackbusch, W. Multi-Grid Methods and Applications; Springer: New York, NY, USA, 1985. [Google Scholar]
  14. Trottenberg, U.; Oosterlee, C.; Schüller, A. Multigrid; Academic Press: London, UK, 2001. [Google Scholar]
  15. Benzi, M.; Golub, G.H.; Liesen, J. Numerical solution of saddle point problems. Acta Numer. 2005, 14, 1–137. [Google Scholar] [CrossRef] [Green Version]
  16. Zulehner, W. Nonstandard norms and robust estimates for saddle point problems. SIAM J. Matrix Anal. Appl. 2011, 32, 536–560. [Google Scholar] [CrossRef] [Green Version]
  17. Brandt, A.; Dinar, N. Multigrid Solutions to Elliptic Flow Problems, Numerical Methods for Partial Differential Equations; Parter, S., Ed.; Academic Press: New York, NY, USA; London, UK, 1979; pp. 53–147. [Google Scholar]
  18. Brandt, A.; Yavneh, I. On Multigrid Solution of High-Reynolds Incompressible Entering Flows. J. Comput. Phys. 1992, 101, 151–164. [Google Scholar] [CrossRef] [Green Version]
  19. Girault, V.; Raviart, P.-A. Finite Element Methods for Navier–Stokes Equations: Theory and Algorithms; Springer Series in Computational Mathematics; Springer: Berlin, Germany, 1986; Volume 5. [Google Scholar]
  20. Kellogg, R.B.; Osborn, J.E. A regularity result for the Stokes problem in a convex polygon. J. Funct. Anal. 1976, 21, 397–431. [Google Scholar] [CrossRef] [Green Version]
  21. Oosterlee, C.W.; Gaspar, F.J. Multigrid methods for the Stokes system. Comput. Sci. Eng. 2006, 8, 34–43. [Google Scholar] [CrossRef]
  22. Takacs, S. A robust all-at-once multigrid method for the Stokes control problem. Numer. Math. 2015, 130, 517–540. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, W.; Chen, L. Multigrid Methods for the Stokes Equations using Distributive Gauss-Seidel Relaxations based on the Least Squares Commutator. J. Sci. Comput. 2013, 56, 409–431. [Google Scholar] [CrossRef]
  24. Wittum, G. Multi-grid methods for Stokes and Navier-Stokes equations. Numer. Math. 1989, 54, 543–563. [Google Scholar] [CrossRef]
  25. Borzì, A.; Schulz, V. Multigrid methods for PDE optimization. SIAM Rev. 2009, 51, 361–395. [Google Scholar] [CrossRef]
  26. Rees, T.; Wathen, A.J. Preconditioning Iterative Methods for the Optimal Control of the Stokes Equations. SIAM J. Sci. Comput. 2011, 33, 2903–2926. [Google Scholar] [CrossRef]
  27. Butt, M.M. A multigrid solver for Stokes control problems. Int. J. Comput. Math. 2017, 94, 2297–2314. [Google Scholar] [CrossRef]
  28. Butt, M.M.; Yuan, Y. A full multigrid method for distributed control problems constrained by Stokes equations. Numer. Math. Theory Methods Appl. 2017, 10, 639–655. [Google Scholar] [CrossRef]
  29. Noack, B.R.; Morzynski, M.; Tadmor, G. Reduced-Order Modelling for Flow Control, 1st ed.; Springer: Wien, Austria, 2011. [Google Scholar]
  30. Ito, K.; Kunisch, K. Lagrange Multiplier Approach to Variational Problems and Applications; SIAM: Philadelphia, PA, USA, 2008. [Google Scholar]
  31. Butt, M.M.; Borzì, A. Formulation and multigrid solution of Cauchy-Riemann optimal control problems. Comput. Vis. Sci. 2011, 14, 79–90. [Google Scholar] [CrossRef]
Figure 1. Three-factor coarsening on staggered grids.
Figure 1. Three-factor coarsening on staggered grids.
Mathematics 09 00738 g001
Figure 2. Full multigrid (FMG) method on three levels. The red circles denotes the initial guess after prolongation ↗ (interpolation) and ↘ shows the restriction.
Figure 2. Full multigrid (FMG) method on three levels. The red circles denotes the initial guess after prolongation ↗ (interpolation) and ↘ shows the restriction.
Mathematics 09 00738 g002
Figure 3. The realization of the white noise W ˙ h (left) and the target velocity field (right) on 18 2 mesh ( h = 1 / 18 ).
Figure 3. The realization of the white noise W ˙ h (left) and the target velocity field (right) on 18 2 mesh ( h = 1 / 18 ).
Mathematics 09 00738 g003
Figure 4. Optimal solution (velocity vector field): (top left) α = 10 1 ; (top right) α = 10 2 ; (bottom left) α = 10 3 ; and (bottom right) α = 10 4 on 18 2 mesh.
Figure 4. Optimal solution (velocity vector field): (top left) α = 10 1 ; (top right) α = 10 2 ; (bottom left) α = 10 3 ; and (bottom right) α = 10 4 on 18 2 mesh.
Mathematics 09 00738 g004
Table 1. Tracking errors on 18 2 mesh.
Table 1. Tracking errors on 18 2 mesh.
α u h U h L 2 ( Ω h ) J ^ ( f h )
10 1 4.3742 × 10−11.9203 × 10−1
10 2 4.2916 × 10−11.9080 × 10−1
10 3 3.0231 × 10−11.2330 × 10−1
10 4 9.7101 × 10 2 3.8977 × 10 2
10 5 2.3830 × 10 2 6.0480 × 10 3
Table 2. Iteration counts and CPU time for stochastic Stokes control problem.
Table 2. Iteration counts and CPU time for stochastic Stokes control problem.
N x × N y 6 2 18 2 54 2
α = 10 1 141009
CPU(sec)(0.43)(0.63)(1.09)
α = 10 2 14610080
CPU(sec)(0.90)(1.87)(4.99)
α = 10 3 1383999741
CPU(sec)(4.04)(12.88)(35.13)
α = 10 4 856686066740
CPU(sec)(20.86)(73.80)(522.20)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Butt, M.M. Multigrid Method for Optimal Control Problem Constrained by Stochastic Stokes Equations with Noise. Mathematics 2021, 9, 738. https://doi.org/10.3390/math9070738

AMA Style

Butt MM. Multigrid Method for Optimal Control Problem Constrained by Stochastic Stokes Equations with Noise. Mathematics. 2021; 9(7):738. https://doi.org/10.3390/math9070738

Chicago/Turabian Style

Butt, Muhammad Munir. 2021. "Multigrid Method for Optimal Control Problem Constrained by Stochastic Stokes Equations with Noise" Mathematics 9, no. 7: 738. https://doi.org/10.3390/math9070738

APA Style

Butt, M. M. (2021). Multigrid Method for Optimal Control Problem Constrained by Stochastic Stokes Equations with Noise. Mathematics, 9(7), 738. https://doi.org/10.3390/math9070738

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