### Introduction

This tutorial aims to illustrate the application of the adjoint-based approach for finding invariant solutions of the Navier-Stokes equations that was developed in Ref. [1]. The accompanying MATLAB codes may be used with proper citation to Ref. [1].

### Overview of the method

Consider a set of partial differential equations (PDEs)

where is a (nonlinear) differential operator acting on the function . We assume that the differential operator only contains partial derivatives with respect to the spatial variable . We are interested in finding (relative) equilibrium solutions of PDE (1). These equilibria are typically unstable.

The goal is to derive a second PDE

whose equilibria are stable and contain those of PDE (1) and further the solutions of PDE (2) converge asymptotically to its equilibria for almost all initial conditions . Note that the parameter is a fictitious time and does not in general coincide with time .

It is shown in Ref. [2] in the context of nonlinear waves and in Ref. [1] for the incompressible Navier-Stokes equation that the choice accomplishes the above-mentioned goals. That is the solutions of

converge asymptotically to the equilibria of Eq. (1). We refer to the above equation as the *adjoint PDE*.

Here the operator is the Gateaux derivative (or linearization) of and is its adjoint with respect to some appropriate inner product.

### Example

Consider the Kuramoto-Sivashinsky equation

with periodic boundary conditions.

The domain size plays the role of the “Reynolds number” here. It is known that for a large enough domain size , the invariant solutions of Kuramoto-Sivashinsky equation become unstable and the dynamics becomes chaotic. Ref. [3] presents a thorough study of the state space of the Kuramoto-Sivashinsky equation for . For the rest of this tutorial we set . The movie below shows a solution of the Kuramoto-Sivashinsky equation starting from a random phase initial condition.

For the Kuramoto-Sivashinsky equation, the nonlinear differential operator is given by

whose Gateaux derivative reads

The adjoint of with respect to the inner product is

One can readily check that

Therefore, the adjoint descent PDE reads

We numerically solve this adjoint PDE with initial conditions In the following, we set .

Ref. [3] shows that the Kuramoto-Sivashinsky equation has three non-trivial equilibrium solutions for . These equilibria are labeled as E1, E2 and E3 in Figure 3 of Ref. [3]. We will use the same labeling here and find these equilibrium solutions by numerically integrating the adjoint PDE.

The movie below shows equilibrium E1 obtained from the initial condition . The red dashed curve marks the initial condition and the black curve shows the evolution of the solution of the adjoint PDE. The panel on the right shows the evolution of the residue as the fictitious time increases.

**Note: **The solution of the adjoint PDE converges very quickly to the vicinity of the equilibrium E1 such that, after approximately the first 20 fictitious time units, it is visually indistinguishable from the equilibrium itself.

The movie below shows equilibrium E2 obtained from the initial condition . The red dashed curve marks the initial condition and the black curve shows the evolution of the solution of the adjoint PDE. The panel on the right shows the evolution of the residue as the fictitious time increases.

The movie below shows equilibrium E3 obtained from the initial condition . The red dashed curve marks the initial condition and the black curve shows the evolution of the solution of the adjoint PDE. The panel on the right shows the evolution of the residue as the fictitious time increases.

We found that for initial conditions with , the solutions of the adjoint PDE converge to the trivial solution . The movie below shows an example with , i.e.,

### MATLAB code

This MATLAB code may be used with proper citation to Ref. [1].

### References

- Farazmand, M., An adjoint-based approach for finding invariant solutions of Navier-Stokes equations, J. Fluid Mech., vol. 795, pp. 278-312, 2016
- Yang, J. & Lakoba, T. I., Universally-convergent squared-operator iteration methods for solitary waves in general nonlinear wave equations. Studies in Applied Mathematics 118 (2), 153–197, 2007
- Cvitanovic, P., Ruslan L. D., and Evangelos S., On the state space geometry of the Kuramoto-Sivashinsky flow in a periodic domain.,
*SIAM Journal on Applied Dynamical Systems*9.1, 1-33, 2010