physics module for 1D,2D and 3D non-isothermal Euler equations More...
Data Types | |
interface | cons2prim |
type | physics_euler |
interface | prim2cons |
interface | setflux |
interface | statevector_euler |
Functions/Subroutines | |
subroutine | initphysics (this, Mesh, config, IO) |
constructor of physics_euler class More... | |
subroutine | new_statevector (this, new_sv, flavour, num) |
allocate and initialize new non-isothermal state vector More... | |
pure subroutine | convert2primitive_all (this, cvar, pvar) |
Converts conservative to primitive variables on the whole mesh. More... | |
pure subroutine | convert2primitive_subset (this, i1, i2, j1, j2, k1, k2, cvar, pvar) |
Converts conservative to primitive variables on a subset of the data. More... | |
pure subroutine | convert2conservative_all (this, pvar, cvar) |
Converts primitive to conservative variables on the whole mesh. More... | |
pure subroutine | convert2conservative_subset (this, i1, i2, j1, j2, k1, k2, pvar, cvar) |
Converts primitive to conservative variables on a subset of the data. More... | |
pure subroutine | calcfluxesx (this, Mesh, nmin, nmax, prim, cons, xfluxes) |
Calculate Fluxes in x-direction. More... | |
pure subroutine | calcfluxesy (this, Mesh, nmin, nmax, prim, cons, yfluxes) |
Calculate Fluxes in y-direction. More... | |
pure subroutine | calcfluxesz (this, Mesh, nmin, nmax, prim, cons, zfluxes) |
Calculate Fluxes in z-direction. More... | |
pure subroutine | calculatecharsystemx (this, Mesh, i1, i2, pvar, lambda, xvar) |
pure subroutine | calculatecharsystemy (this, Mesh, j1, j2, pvar, lambda, xvar) |
pure subroutine | calculatecharsystemz (this, Mesh, k1, k2, pvar, lambda, xvar) |
pure subroutine | calculateboundarydatax (this, Mesh, i1, i2, xvar, pvar) |
extrapolate pvar using characteristic pseudo variables (absorbing boundaries) More... | |
pure subroutine | calculateboundarydatay (this, Mesh, j1, j2, xvar, pvar) |
pure subroutine | calculateboundarydataz (this, Mesh, k1, k2, xvar, pvar) |
pure subroutine | calculateprim2riemannx (this, Mesh, i, pvar, lambda, Rinv) |
Conversion from primitive to riemann invariants for farfield boundaries. More... | |
pure subroutine | calculateprim2riemanny (this, Mesh, j, pvar, lambda, Rinv) |
Conversion from primitive to riemann invariants for farfield boundaries. More... | |
pure subroutine | calculateprim2riemannz (this, Mesh, k, pvar, lambda, Rinv) |
Conversion from primitive to riemann invariants for farfield boundaries. More... | |
pure subroutine | calculateriemann2primx (this, Mesh, i, Rinv, pvar) |
Convert Riemann invariants to primitives for farfield boundaries. More... | |
pure subroutine | calculateriemann2primy (this, Mesh, j, Rinv, pvar) |
Convert Riemann invariants to primitives for farfield boundaries. More... | |
pure subroutine | calculateriemann2primz (this, Mesh, k, Rinv, pvar) |
Convert Riemann invariants to primitives for farfield boundaries. More... | |
pure subroutine | geometricalsources (this, Mesh, pvar, cvar, sterm) |
Calculates geometrical sources. More... | |
subroutine | externalsources (this, accel, pvar, cvar, sterm) |
compute momentum and energy sources given an external force More... | |
subroutine | viscositysources (this, Mesh, pvar, btxx, btxy, btxz, btyy, btyz, btzz, sterm) |
pure subroutine | addbackgroundvelocityx (this, Mesh, w, pvar, cvar) |
Adds a background velocity field for fargo routines. More... | |
pure subroutine | addbackgroundvelocityy (this, Mesh, w, pvar, cvar) |
Adds a background velocity field for fargo routines. More... | |
pure subroutine | addbackgroundvelocityz (this, Mesh, w, pvar, cvar) |
Adds a background velocity field for fargo routines. More... | |
pure subroutine | subtractbackgroundvelocityx (this, Mesh, w, pvar, cvar) |
Substracts a background velocity field for fargo routines. More... | |
pure subroutine | subtractbackgroundvelocityy (this, Mesh, w, pvar, cvar) |
Substracts a background velocity field for fargo routines. More... | |
pure subroutine | subtractbackgroundvelocityz (this, Mesh, w, pvar, cvar) |
Substracts a background velocity field for fargo routines. More... | |
pure subroutine | addfargosourcesx (this, Mesh, w, pvar, cvar, sterm) |
sources terms for fargo advection along x-direction More... | |
pure subroutine | addfargosourcesy (this, Mesh, w, pvar, cvar, sterm) |
sources terms for fargo advection along y-direction More... | |
pure subroutine | addfargosourcesz (this, Mesh, w, pvar, cvar, sterm) |
sources terms for fargo advection along z-direction More... | |
pure subroutine | updatesoundspeed (this, pvar) |
subroutine | finalize (this) |
Destructor of the physics_euler class. More... | |
type(statevector_euler) function | createstatevector (Physics, flavour, num) |
Constructor of statevector_euler. More... | |
subroutine | assignmarray_0 (this, ma) |
assigns one state vector to another state vector More... | |
subroutine | finalize_statevector (this) |
actual destructor of the statevector_eulerisotherm type More... | |
elemental real function | getsoundspeed (gamma, density, pressure) |
elemental subroutine | setroeaverages (gamma, rhoL, rhoR, ul, uR, vL, vR, pL, pR, eL, eR, u, cs) |
elemental subroutine | setwavespeeds (cs, v, minwav, maxwav) |
set minimal and maximal wave speeds More... | |
elemental subroutine | seteigenvalues1d (gamma, rho, v, P, l1, l2, l3) |
set all eigenvalues for 1D transport (used in absorbing boundary conditions) More... | |
elemental subroutine | seteigenvalues2d (gamma, rho, v, P, l1, l2, l3, l4) |
set all eigenvalues for 2D transport (used in absorbing boundary conditions) More... | |
elemental subroutine | seteigenvalues3d (gamma, rho, v, P, l1, l2, l3, l4, l5) |
set all eigenvalues for 3D transport (used in absorbing boundary conditions) More... | |
elemental subroutine | setintermediatestate (rhoL, rhoR, uL, uR, vl, vR, pL, pR, eL, eR, amin, amax, rho, mu, mv, e, a) |
elemental subroutine | setcharvars1d (gamma, rho1, rho2, u1, u2, P1, P2, l1, l3, xvar1, xvar2, xvar3) |
compute characteristic variables for 1D transport More... | |
elemental subroutine | setcharvars2d (gamma, rho1, rho2, u1, u2, v1, v2, P1, P2, l1, l4, xvar1, xvar2, xvar3, xvar4) |
compute characteristic variables for 2D transport More... | |
elemental subroutine | setcharvars3d (gamma, rho1, rho2, u1, u2, v1, v2, w1, w2, P1, P2, l1, l5, xvar1, xvar2, xvar3, xvar4, xvar5) |
compute characteristic variables for 3D transport More... | |
elemental subroutine | setboundarydata1d (delta, gamma, rho1, u1, P1, xvar1, xvar2, xvar3, rho2, u2, P2) |
extrapolate primitive variables using characteristic pseudo pevariables 1D transport More... | |
elemental subroutine | setboundarydata2d (delta, gamma, rho1, u1, v1, P1, xvar1, xvar2, xvar3, xvar4, rho2, u2, v2, P2) |
extrapolate primitive variables using characteristic pseudo pevariables 2D transport More... | |
elemental subroutine | setboundarydata3d (delta, gamma, rho1, u1, v1, w1, P1, xvar1, xvar2, xvar3, xvar4, xvar5, rho2, u2, v2, w2, P2) |
extrapolate primitive variables using characteristic pseudo pevariables 2D transport More... | |
elemental subroutine | prim2riemann1d (gamma, rho, vx, p, l1, l2, l3, Rminus, Rs, Rplus) |
elemental subroutine | prim2riemann2d (gamma, rho, vx, vy, p, l1, l2, l3, l4, Rminus, Rs, Rvt, Rplus) |
elemental subroutine | prim2riemann3d (gamma, rho, vx, vy, vz, p, l1, l2, l3, l4, l5, Rminus, Rs, Rvt, Rwt, Rplus) |
elemental subroutine | riemann2prim1d (gamma, Rminus, Rs, Rplus, rho, vx, p) |
elemental subroutine | riemann2prim2d (gamma, Rminus, Rs, Rvt, Rplus, rho, vx, vy, p) |
elemental subroutine | riemann2prim3d (gamma, Rminus, Rs, Rvt, Rwt, Rplus, rho, vx, vy, vz, p) |
elemental subroutine | setflux1d (rho, u, P, mu, E, f1, f2, f3) |
set mass, 1D momentum and energy flux for transport along the 1st dimension More... | |
elemental subroutine | setflux2d (rho, u, P, mu, mv, E, f1, f2, f3, f4) |
set mass, 2D momentum and energy flux for transport along the 1st dimension More... | |
elemental subroutine | setflux3d (rho, u, P, mu, mv, mw, E, f1, f2, f3, f4, f5) |
set mass, 3D momentum and energy flux for transport along the 1st dimension More... | |
elemental subroutine | cons2prim1d (gamma, rho_in, mu, E, rho_out, u, P) |
Convert from 1D conservative to primitive variables. More... | |
elemental subroutine | cons2prim2d (gamma, rho_in, mu, mv, E, rho_out, u, v, P) |
Convert from 2D conservative to primitive variables. More... | |
elemental subroutine | cons2prim3d (gamma, rho_in, mu, mv, mw, E, rho_out, u, v, w, P) |
Convert from 3D conservative to primitive variables. More... | |
elemental subroutine | prim2cons1d (gamma, rho_in, u, P, rho_out, mu, E) |
Convert from 1D primitive to conservative variables. More... | |
elemental subroutine | prim2cons2d (gamma, rho_in, u, v, P, rho_out, mu, mv, E) |
Convert from 2D primitive to conservative variables. More... | |
elemental subroutine | prim2cons3d (gamma, rho_in, u, v, w, P, rho_out, mu, mv, mw, E) |
Convert from 3D primitive to conservative variables. More... | |
elemental real function | getgeometricalsourcex (cxyx, cxzx, cyxy, czxz, vx, vy, vz, P, my, mz) |
geometrical momentum source terms P is the either isothermal pressure rho*cs**2 or the real pressure. More... | |
elemental real function | getgeometricalsourcey (cxyx, cyxy, cyzy, czyz, vx, vy, vz, P, mx, mz) |
y-momentum geometrical source term More... | |
elemental real function | getgeometricalsourcez (cxzx, cyzy, czxz, czyz, vx, vy, vz, P, mx, my) |
z-momentum geometrical source term More... | |
Variables | |
character(len=32), parameter | problem_name = "Euler" |
Detailed Description
physics module for 1D,2D and 3D non-isothermal Euler equations
Function/Subroutine Documentation
◆ addbackgroundvelocityx()
|
private |
Adds a background velocity field for fargo routines.
Calculates
\begin{eqnarray*} E &=& E' + m_x' w + \frac{1}{2}\varrho w^2 \\ v_x &=& v_x' + w \\ m_x &=& m_x' + \varrho w, \end{eqnarray*}
with \( E, v_y, m_y \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift.
Definition at line 1955 of file physics_euler.f90.
◆ addbackgroundvelocityy()
|
private |
Adds a background velocity field for fargo routines.
Calculates
\begin{eqnarray*} E &=& E' + m_y' w + \frac{1}{2}\varrho w^2 \\ v_y &=& v_y' + w \\ m_y &=& m_y' + \varrho w, \end{eqnarray*}
with \( E, v_y, m_y \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift.
Definition at line 2003 of file physics_euler.f90.
◆ addbackgroundvelocityz()
|
private |
Adds a background velocity field for fargo routines.
Calculates
\begin{eqnarray*} E &=& E' + m_z' w + \frac{1}{2}\varrho w^2 \\ v_z &=& v_z' + w \\ m_z &=& m_z' + \varrho w, \end{eqnarray*}
with \( E, v_z, m_z \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift.
Definition at line 2059 of file physics_euler.f90.
◆ addfargosourcesx()
pure subroutine physics_euler_mod::addfargosourcesx | ( | class(physics_euler), intent(inout) | this, |
class(mesh_base), intent(in) | Mesh, | ||
real, dimension(mesh%jgmin:mesh%jgmax,mesh%kgmin:mesh%kgmax), intent(in) | w, | ||
class(marray_compound), intent(inout) | pvar, | ||
class(marray_compound), intent(inout) | cvar, | ||
class(marray_compound), intent(inout) | sterm | ||
) |
sources terms for fargo advection along x-direction
If the background velocity \(\vec{w}=w\,\hat{e}_\xi\) with \(w\) independent of \(\xi\) and \(t\) is subtracted from the overall velocity of the flow, an additional source term occurs in the \(\xi\)-momentum equation:
\[ S_\mathrm{Fargo} = -\varrho \vec{v} \cdot \nabla \vec{w} \]
where \(\vec{v}\) is the real velocity (including the background velocity \(\vec{w}\) ).
Definition at line 2267 of file physics_euler.f90.

◆ addfargosourcesy()
pure subroutine physics_euler_mod::addfargosourcesy | ( | class(physics_euler), intent(inout) | this, |
class(mesh_base), intent(in) | Mesh, | ||
real, dimension(mesh%igmin:mesh%igmax,mesh%kgmin:mesh%kgmax), intent(in) | w, | ||
class(marray_compound), intent(inout) | pvar, | ||
class(marray_compound), intent(inout) | cvar, | ||
class(marray_compound), intent(inout) | sterm | ||
) |
sources terms for fargo advection along y-direction
If the background velocity \(\vec{w}=w\,\hat{e}_\eta\) with \(w\) independent of \(\eta\) and \(t\) is subtracted from the overall velocity of the flow, an additional source term occurs in the \(\eta\)-momentum equation:
\[ S_\mathrm{Fargo} = -\varrho \vec{v} \cdot \nabla \vec{w} \]
where \(\vec{v}\) is the real velocity (including the background velocity \(\vec{w}\) ).
Definition at line 2368 of file physics_euler.f90.

◆ addfargosourcesz()
pure subroutine physics_euler_mod::addfargosourcesz | ( | class(physics_euler), intent(inout) | this, |
class(mesh_base), intent(in) | Mesh, | ||
real, dimension(mesh%igmin:mesh%igmax,mesh%jgmin:mesh%jgmax), intent(in) | w, | ||
class(marray_compound), intent(inout) | pvar, | ||
class(marray_compound), intent(inout) | cvar, | ||
class(marray_compound), intent(inout) | sterm | ||
) |
sources terms for fargo advection along z-direction
If the background velocity \(\vec{w}=w\,\hat{e}_phi\) with \(w\) independent of \(\phi\) and \(t\) is subtracted from the overall velocity of the flow, an additional source term occurs in the \(\phi\)-momentum equation:
\[ S_\mathrm{Fargo} = -\varrho \vec{v} \cdot \nabla \vec{w} \]
where \(\vec{v}\) is the real velocity (including the background velocity \(\vec{w}\) ).
Definition at line 2468 of file physics_euler.f90.

◆ assignmarray_0()
|
private |
assigns one state vector to another state vector
Definition at line 2632 of file physics_euler.f90.
◆ calcfluxesx()
|
private |
Calculate Fluxes in x-direction.
Definition at line 492 of file physics_euler.f90.
◆ calcfluxesy()
|
private |
Calculate Fluxes in y-direction.
Definition at line 548 of file physics_euler.f90.
◆ calcfluxesz()
|
private |
Calculate Fluxes in z-direction.
Definition at line 604 of file physics_euler.f90.
◆ calculateboundarydatax()
|
private |
extrapolate pvar using characteristic pseudo variables (absorbing boundaries)
Definition at line 1072 of file physics_euler.f90.

◆ calculateboundarydatay()
|
private |
◆ calculateboundarydataz()
|
private |
◆ calculatecharsystemx()
|
private |
◆ calculatecharsystemy()
|
private |
◆ calculatecharsystemz()
|
private |
◆ calculateprim2riemannx()
|
private |
Conversion from primitive to riemann invariants for farfield boundaries.
Definition at line 1266 of file physics_euler.f90.

◆ calculateprim2riemanny()
|
private |
Conversion from primitive to riemann invariants for farfield boundaries.
Definition at line 1360 of file physics_euler.f90.

◆ calculateprim2riemannz()
|
private |
Conversion from primitive to riemann invariants for farfield boundaries.
Definition at line 1456 of file physics_euler.f90.

◆ calculateriemann2primx()
|
private |
Convert Riemann invariants to primitives for farfield boundaries.
Definition at line 1555 of file physics_euler.f90.

◆ calculateriemann2primy()
|
private |
Convert Riemann invariants to primitives for farfield boundaries.
Definition at line 1606 of file physics_euler.f90.

◆ calculateriemann2primz()
|
private |
Convert Riemann invariants to primitives for farfield boundaries.
Definition at line 1657 of file physics_euler.f90.

◆ cons2prim1d()
|
private |
Convert from 1D conservative to primitive variables.
Definition at line 3033 of file physics_euler.f90.
◆ cons2prim2d()
|
private |
Convert from 2D conservative to primitive variables.
Definition at line 3048 of file physics_euler.f90.
◆ cons2prim3d()
|
private |
Convert from 3D conservative to primitive variables.
Definition at line 3064 of file physics_euler.f90.
◆ convert2conservative_all()
|
private |
Converts primitive to conservative variables on the whole mesh.
Definition at line 367 of file physics_euler.f90.
◆ convert2conservative_subset()
|
private |
Converts primitive to conservative variables on a subset of the data.
Definition at line 405 of file physics_euler.f90.
◆ convert2primitive_all()
|
private |
Converts conservative to primitive variables on the whole mesh.
Definition at line 243 of file physics_euler.f90.
◆ convert2primitive_subset()
|
private |
Converts conservative to primitive variables on a subset of the data.
Definition at line 281 of file physics_euler.f90.
◆ createstatevector()
type(statevector_euler) function physics_euler_mod::createstatevector | ( | class(physics_euler), intent(in) | Physics, |
integer, intent(in), optional | flavour, | ||
integer, intent(in), optional | num | ||
) |
Constructor of statevector_euler.
- methods of class statevector_euler
This is not a class member itself, instead its an ordinary module procedure. The function name is overloaded with the class name.
Definition at line 2593 of file physics_euler.f90.

◆ externalsources()
subroutine physics_euler_mod::externalsources | ( | class(physics_euler), intent(in) | this, |
class(marray_base), intent(in) | accel, | ||
class(marray_compound), intent(in) | pvar, | ||
class(marray_compound), intent(in) | cvar, | ||
class(marray_compound), intent(inout) | sterm | ||
) |
compute momentum and energy sources given an external force
Definition at line 1849 of file physics_euler.f90.
◆ finalize()
subroutine physics_euler_mod::finalize | ( | class(physics_euler), intent(inout) | this | ) |
Destructor of the physics_euler class.
Definition at line 2577 of file physics_euler.f90.
◆ finalize_statevector()
|
private |
actual destructor of the statevector_eulerisotherm type
Definition at line 2666 of file physics_euler.f90.

◆ geometricalsources()
|
private |
Calculates geometrical sources.
Definition at line 1706 of file physics_euler.f90.

◆ getgeometricalsourcex()
|
private |
geometrical momentum source terms P is the either isothermal pressure rho*cs**2 or the real pressure.
- Attention
- These elemental functions exist multiple times for performance reasons (inlining). Please keep this in mind for changes. Other modules with this function:
x-momentum geometrical source term
Definition at line 3128 of file physics_euler.f90.


◆ getgeometricalsourcey()
|
private |
y-momentum geometrical source term
Definition at line 3139 of file physics_euler.f90.


◆ getgeometricalsourcez()
|
private |
z-momentum geometrical source term
Definition at line 3150 of file physics_euler.f90.


◆ getsoundspeed()
|
private |
- elemental non-class subroutines / functions
Definition at line 2682 of file physics_euler.f90.

◆ initphysics()
|
private |
constructor of physics_euler class
- Todo:
- remove / improve in future version set array indices for 1st,2nd,3rd non-vanishing velocities this may actually not coincide with the x,y and z-velocities
Definition at line 135 of file physics_euler.f90.
◆ new_statevector()
subroutine physics_euler_mod::new_statevector | ( | class(physics_euler), intent(in) | this, |
class(marray_compound), pointer | new_sv, | ||
integer, intent(in), optional | flavour, | ||
integer, intent(in), optional | num | ||
) |
allocate and initialize new non-isothermal state vector
Definition at line 225 of file physics_euler.f90.
◆ prim2cons1d()
|
private |
Convert from 1D primitive to conservative variables.
Definition at line 3081 of file physics_euler.f90.
◆ prim2cons2d()
|
private |
Convert from 2D primitive to conservative variables.
Definition at line 3093 of file physics_euler.f90.
◆ prim2cons3d()
|
private |
Convert from 3D primitive to conservative variables.
Definition at line 3106 of file physics_euler.f90.
◆ prim2riemann1d()
|
private |
◆ prim2riemann2d()
|
private |
Definition at line 2913 of file physics_euler.f90.


◆ prim2riemann3d()
|
private |
Definition at line 2930 of file physics_euler.f90.


◆ riemann2prim1d()
|
private |
◆ riemann2prim2d()
|
private |
Definition at line 2965 of file physics_euler.f90.


◆ riemann2prim3d()
|
private |
Definition at line 2983 of file physics_euler.f90.


◆ setboundarydata1d()
|
private |
extrapolate primitive variables using characteristic pseudo pevariables 1D transport
Definition at line 2843 of file physics_euler.f90.


◆ setboundarydata2d()
|
private |
extrapolate primitive variables using characteristic pseudo pevariables 2D transport
Definition at line 2865 of file physics_euler.f90.


◆ setboundarydata3d()
|
private |
extrapolate primitive variables using characteristic pseudo pevariables 2D transport
Definition at line 2879 of file physics_euler.f90.


◆ setcharvars1d()
|
private |
compute characteristic variables for 1D transport
Definition at line 2799 of file physics_euler.f90.

◆ setcharvars2d()
|
private |
compute characteristic variables for 2D transport
Definition at line 2818 of file physics_euler.f90.


◆ setcharvars3d()
|
private |
compute characteristic variables for 3D transport
Definition at line 2830 of file physics_euler.f90.


◆ seteigenvalues1d()
|
private |
set all eigenvalues for 1D transport (used in absorbing boundary conditions)
Definition at line 2729 of file physics_euler.f90.


◆ seteigenvalues2d()
|
private |
set all eigenvalues for 2D transport (used in absorbing boundary conditions)
Definition at line 2746 of file physics_euler.f90.


◆ seteigenvalues3d()
|
private |
set all eigenvalues for 3D transport (used in absorbing boundary conditions)
Definition at line 2757 of file physics_euler.f90.


◆ setflux1d()
|
private |
set mass, 1D momentum and energy flux for transport along the 1st dimension
Definition at line 2999 of file physics_euler.f90.

◆ setflux2d()
|
private |
set mass, 2D momentum and energy flux for transport along the 1st dimension
Definition at line 3011 of file physics_euler.f90.

◆ setflux3d()
|
private |
set mass, 3D momentum and energy flux for transport along the 1st dimension
Definition at line 3022 of file physics_euler.f90.
◆ setintermediatestate()
|
private |
Definition at line 2769 of file physics_euler.f90.
◆ setroeaverages()
|
private |
- Todo:
- NOT VERIFIED only for advanced wavespeeds
Definition at line 2693 of file physics_euler.f90.
◆ setwavespeeds()
|
private |
set minimal and maximal wave speeds
Definition at line 2717 of file physics_euler.f90.
◆ subtractbackgroundvelocityx()
|
private |
Substracts a background velocity field for fargo routines.
Calculates
\begin{eqnarray*} E' &=& E - m_x w + \frac{1}{2}\varrho w^2 \\ v_x' &=& v_x - w \\ m_x' &=& m_x - \varrho w, \end{eqnarray*}
with \( E, v_x, m_x \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift. ATTENTION: the "+" before the 3rd term in the energy transformation is correct!
Definition at line 2109 of file physics_euler.f90.
◆ subtractbackgroundvelocityy()
|
private |
Substracts a background velocity field for fargo routines.
Calculates
\begin{eqnarray*} E' &=& E - m_y w + \frac{1}{2}\varrho w^2 \\ v_y' &=& v_y - w \\ m_y' &=& m_y - \varrho w, \end{eqnarray*}
with \( E, v_y, m_y \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift. ATTENTION: the "+" before the 3rd term in the energy transformation is correct!
Definition at line 2159 of file physics_euler.f90.
◆ subtractbackgroundvelocityz()
|
private |
Substracts a background velocity field for fargo routines.
Calculates
\begin{eqnarray*} E' &=& E - m_z w + \frac{1}{2}\varrho w^2 \\ v_z' &=& v_z - w \\ m_z' &=& m_z - \varrho w, \end{eqnarray*}
with \( E, v_z, m_z \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift. ATTENTION: the "+" before the 3rd term in the energy transformation is correct!
Definition at line 2217 of file physics_euler.f90.
◆ updatesoundspeed()
|
private |
◆ viscositysources()
|
private |
Definition at line 1887 of file physics_euler.f90.
Variable Documentation
◆ problem_name
|
private |
Definition at line 52 of file physics_euler.f90.