Basilisk CFD
Adaptive Cartesian mesh PDE framework
Loading...
Searching...
No Matches
NASG.h File Reference

Go to the source code of this file.

Macros

#define PIGAMMA
 

Functions

double sound_speed (Point point)
 These functions are provided by the Equation Of State.
 
double average_pressure (Point point)
 
double bulk_compressibility (Point point)
 
double internal_energy (Point point, double fc)
 
double average_temperature (Point point, double p)
 These functions are provided by the Equation Of State.
 
double thermal_expansion (Point point)
 

Variables

double gamma1 = 1.4 [0]
 
double gamma2 = 1.4 [0]
 
double PI1 = 0.
 
double PI2 = 0.
 
double b1 = 0.
 
double b2 = 0.
 
double q1 = 0.
 
double q2 = 0.
 
double cv1 = 0.
 
double cv2 = 0.
 

Macro Definition Documentation

◆ PIGAMMA

#define PIGAMMA
Value:
double invgammaavg = (fc - frho1[]*b1)/(gamma1 - 1.) + \
(1. - fc - frho2[]*b2)/(gamma2 - 1.), \
PIGAMMAavg = (fc - frho1[]*b1)*PI1*gamma1/(gamma1 - 1.) + frho1[]*q1 + \
(1. - fc - frho2[]*b2)*PI2*gamma2/(gamma2 - 1.) + frho2[]*q2
double gamma2
Definition NASG.h:18
double b2
Definition NASG.h:19
double b1
Definition NASG.h:19
double gamma1
Definition NASG.h:18
double PI1
Definition NASG.h:18
double q2
Definition NASG.h:20
double PI2
Definition NASG.h:18
double q1
Definition NASG.h:20
int x
Definition common.h:76
scalar frho2[]
Definition two-phase.h:57
scalar frho1[]
Definition two-phase.h:57

Average pressure

Definition at line 57 of file NASG.h.

Function Documentation

◆ average_pressure()

double average_pressure ( Point  point)

Definition at line 62 of file NASG.h.

References clamp(), dimension, f, fE1, fE2, frho1, frho2, PIGAMMA, q, sq(), vector::x, and x.

Referenced by event_properties().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ average_temperature()

double average_temperature ( Point  point,
double  p 
)

These functions are provided by the Equation Of State.

Average temperature

Definition at line 99 of file NASG.h.

References b1, b2, clamp(), cp1, cp2, cv1, cv2, f, frho1, frho2, p, PI1, PI2, and x.

Referenced by event_acceleration().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ bulk_compressibility()

double bulk_compressibility ( Point  point)

Bulk compressibility of the mixture

i.e. \(\rho c^2\).

Definition at line 78 of file NASG.h.

References b1, b2, clamp(), f, frho1, frho2, gamma1, gamma2, p, PI1, PI2, and x.

Referenced by event_properties().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ internal_energy()

double internal_energy ( Point  point,
double  fc 
)

Internal energy

Definition at line 90 of file NASG.h.

References p, PIGAMMA, and x.

Referenced by fE_refine().

Here is the caller graph for this function:

◆ sound_speed()

double sound_speed ( Point  point)

These functions are provided by the Equation Of State.

Sound speed

In mixture cells, this function returns the maximum between the speeds in both phases.

Definition at line 29 of file NASG.h.

References clamp(), dimension, f, fE1, fE2, frho1, frho2, gamma1, gamma2, max, p, PI1, PI2, q, sq(), vector::x, and x.

Referenced by event_stability().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ thermal_expansion()

double thermal_expansion ( Point  point)

Thermal expansion coefficient

Definition at line 111 of file NASG.h.

References b2, clamp(), cv2, f, gamma2, PI2, ps, Ts, and x.

Referenced by event_acceleration().

Here is the call graph for this function:
Here is the caller graph for this function:

Variable Documentation

◆ b1

◆ b2

◆ cv1

double cv1 = 0.

Definition at line 21 of file NASG.h.

Referenced by average_temperature().

◆ cv2

double cv2 = 0.

Definition at line 21 of file NASG.h.

Referenced by average_temperature(), and thermal_expansion().

◆ gamma1

double gamma1 = 1.4 [0]

The Noble-Abel Stiffened-Gas (NASG) Equation Of State

This EOS is typically used in combination with the two-phase compressible solver with thermal effects.

The general form of the NASG EOS (Le Métayer & Saurel, 2016) is

\[ \rho_i e_i = \frac{p_i + \Gamma_i \Pi_i}{\Gamma_i - 1}(1 - \rho_i b_i) + \rho_i q_i \]

with \(\rho_i\), \(e_i\) and \(p_i\) the densities, internal energies and pressures of each phase.

These are the coefficients of the NASG EOS for each phase.

Definition at line 18 of file NASG.h.

Referenced by bulk_compressibility(), and sound_speed().

◆ gamma2

double gamma2 = 1.4 [0]

Definition at line 18 of file NASG.h.

Referenced by bulk_compressibility(), sound_speed(), and thermal_expansion().

◆ PI1

double PI1 = 0.

Definition at line 18 of file NASG.h.

Referenced by average_temperature(), bulk_compressibility(), and sound_speed().

◆ PI2

double PI2 = 0.

Definition at line 18 of file NASG.h.

Referenced by average_temperature(), bulk_compressibility(), sound_speed(), and thermal_expansion().

◆ q1

scalar q1 = 0.

Definition at line 20 of file NASG.h.

Referenced by event_vof().

◆ q2

scalar q2 = 0.

Definition at line 20 of file NASG.h.

Referenced by event_vof().