PLUTO  4.4-patch2
Functions
rhd_pressure_fix.c File Reference

Inversion scheme for RHD using a pressure fix. More...

#include "pluto.h"

Functions

int RHD_PressureFix (double *uc, double *vc)
 

Detailed Description

Fix pressure to g_smallPressure and find the four velocity by solving the momentum-velocity relation (squared):

\[ m/(Dh) - u = 0 \]

using secant method. The two guesses are:

Note that this equation should always have a solution.

Authors
A. Mignone
Date
May 3, 2020

Function Documentation

◆ RHD_PressureFix()

int RHD_PressureFix ( double *  uc,
double *  vc 
)
Returns
Error codes are:
  • 0 = success
  • 1 = v^2 > 1
  • 2 = too many iterations