DBC+Implementation+in+OPTIMA2D

This page describes the implementation of dissipation-based continuation (DBC) in OPTIMA2D.

=DBC Modifications To the PROBE Flow Solver=

We wish to solve the non-linear problem

math R ( \widehat{Q}) = 0 math

To find a solution Q^ that satisfies the above equation, Newton's method is used:

math \mathcal{A}^{(n)}\Delta \widehat{Q}^{(n)} = -R^{(n)} math

where

math \mathcal {A}^{(n)} = \left( \frac{\partial{R}}{\partial{\widehat{Q}}} \right)^{(n)} math

The Newton update is then calculated as

math \widehat{Q}^{(n+1)} = \widehat{Q}^{(n)} + \Delta \widehat{Q}^{(n)} math

Right Hand Side Modifications:
The right hand side (RHS) of the linear system arising at each Newton-Krylov iteration is -R, where R is the residual of the discretized governing equations. For example, the residual of the discretized Euler equations is given by

math R ( \widehat{Q}) = \delta_{\xi} \widehat{E} + \delta_{\eta} \widehat{F} - \nabla_{\xi}D_{\xi} - \nabla_{\eta}D_{\eta} math

Let H be the modified residual given by

math H ( \widehat{Q}) = R ( \widehat{Q}) + \lambda \hspace{0.03in} D ( \widehat{Q}) math

The following subroutines are affected:

BCRHS VBCRHS RHSX RHSY COEF24X COEF24Y FILTERX FILTERY PROBE

Input Parameters for Controlling DBC Performance:
DISSCON LAMDISSMAX DCTOL LAMKILL