11.4 Session file configuration: Linear stability analysis

Stability analyses of incompressible flow involves solving the linearised Navier-Stokes equations

   ′
∂u--
 ∂t + L(U,u′) = -∇p + ν∇2u′,

where L is a linear operator, its adjoint form, or both. The evolution of the linearised Navier-Stokes operator, which evolves a solution from an initial state to a future time t, can be written as

u(t) = A(t)u(0).

The adjoint evolution operator is denoted as A*. This section details the additional configuration options, in addition to the standard configuration options described earlier, relating to performing this task.

11.4.1 Solver Info

11.4.2 Parameters

The following parameters can be specified in the PARAMETERS section of the session file:

11.4.3 Functions

When using the direct solver for stability analysis it is necessary to specify a Forcing function “StabilityCoupledLNS” in the form:

1<FORCING> 
2   <FORCE TYPE="StabilityCoupledLNS"> 
3   </FORCE> 
4</FORCING>

This is required since we need to tell the solver to use the existing field as a forcing function to the direct matrix inverse as part of the Arnoldi iteration.

Note: Examples of the set up of the direct solver stability analysis (and other incompressible Navier-Stokes solvers) can be found in the regression test directory NEKTAR/solvers/IncNavierStokesSolver/Tests. See for example the files PPF_R15000_ModifiedArnoldi_Shift.tst and PPF_R15000_3D.xml noting that some parameters are specified in the .tst files.