Neko 1.99.6
A portable framework for high-order spectral element flow simulations
Loading...
Searching...
No Matches
compressible_res_device Module Reference

Data Types

type  compressible_res_device_t
 
interface  compressible_res_part_coef_mult_hip
 
interface  compressible_res_part_rk_sum_hip
 
interface  compressible_res_part_visc_hip
 
interface  inviscid_res_part_e_flux_hip
 
interface  inviscid_res_part_mx_flux_hip
 
interface  inviscid_res_part_my_flux_hip
 
interface  inviscid_res_part_mz_flux_hip
 

Functions/Subroutines

subroutine advance_primitive_variables_device (rho_field, m_x, m_y, m_z, e, p, u, v, w, ax, ax_stress, coef, gs, h, artificial_visc, mu, kappa, bcs_vel, time, rk_scheme, dt)
 
subroutine evaluate_rhs_device (rhs_rho_field, rhs_m_x, rhs_m_y, rhs_m_z, rhs_e, rho_field, m_x, m_y, m_z, e, p, u, v, w, ax, ax_stress, coef, gs, h, artificial_visc, mu, kappa)
 
subroutine add_navier_stokes_flux_device (visc_m_x, visc_m_y, visc_m_z, visc_e, rho_field, p, u, v, w, mu, kappa, ax, ax_stress, coef)
 Add the physical Navier-Stokes flux contribution to the viscous residual.
 

Variables

logical compressible_res_device_add_physical_flux = .false.
 Whether physical Navier-Stokes fluxes are active for the current step.
 
logical compressible_res_device_add_physical_stress = .false.
 Whether physical viscous stress is active for the current step.
 
real(kind=rp), public compressible_res_device_gamma = 1.4_rp
 Module variable to store thermodynamic parameter set by factory.
 

Function/Subroutine Documentation

◆ add_navier_stokes_flux_device()

subroutine compressible_res_device::add_navier_stokes_flux_device ( type(field_t), intent(inout visc_m_x,
type(field_t), intent(inout visc_m_y,
type(field_t), intent(inout visc_m_z,
type(field_t), intent(inout visc_e,
type(field_t), intent(in rho_field,
type(field_t), intent(in p,
type(field_t), intent(in u,
type(field_t), intent(in v,
type(field_t), intent(in w,
type(field_t), intent(in mu,
type(field_t), intent(in kappa,
class(ax_t), intent(inout ax,
class(ax_t), intent(inout ax_stress,
type(coef_t), intent(inout coef 
)
Parameters
visc_m_xViscous residual for x-momentum.
visc_m_yViscous residual for y-momentum.
visc_m_zViscous residual for z-momentum.
visc_EViscous residual for total energy.
rho_fieldDensity field.
pPressure field.
uX-velocity field.
vY-velocity field.
wZ-velocity field.
muDynamic viscosity field.
kappaThermal conductivity field.
AxMatrix-vector product operator for scalar diffusion.
Ax_stressMatrix-vector product operator for full stress.
coefSpatial discretization coefficients.

Definition at line 804 of file compressible_res_device.F90.

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

◆ advance_primitive_variables_device()

subroutine compressible_res_device::advance_primitive_variables_device ( type(field_t), intent(inout rho_field,
type(field_t), intent(inout m_x,
type(field_t), intent(inout m_y,
type(field_t), intent(inout m_z,
type(field_t), intent(inout e,
type(field_t), intent(in p,
type(field_t), intent(in u,
type(field_t), intent(in v,
type(field_t), intent(in w,
class(ax_t), intent(inout ax,
class(ax_t), intent(inout ax_stress,
type(coef_t), intent(inout coef,
type(gs_t), intent(inout gs,
type(field_t), intent(in h,
type(field_t), intent(in artificial_visc,
type(field_t), intent(in mu,
type(field_t), intent(in kappa,
type(bc_list_t), intent(inout bcs_vel,
type(time_state_t), intent(in time,
class(runge_kutta_time_scheme_t), intent(in rk_scheme,
real(kind=rp), intent(in dt 
)

Definition at line 416 of file compressible_res_device.F90.

Here is the call graph for this function:

◆ evaluate_rhs_device()

subroutine compressible_res_device::evaluate_rhs_device ( type(field_t), intent(inout rhs_rho_field,
type(field_t), intent(inout rhs_m_x,
type(field_t), intent(inout rhs_m_y,
type(field_t), intent(inout rhs_m_z,
type(field_t), intent(inout rhs_e,
type(field_t), intent(inout rho_field,
type(field_t), intent(inout m_x,
type(field_t), intent(inout m_y,
type(field_t), intent(inout m_z,
type(field_t), intent(inout e,
type(field_t), intent(in p,
type(field_t), intent(in u,
type(field_t), intent(in v,
type(field_t), intent(in w,
class(ax_t), intent(inout ax,
class(ax_t), intent(inout ax_stress,
type(coef_t), intent(inout coef,
type(gs_t), intent(inout gs,
type(field_t), intent(in h,
type(field_t), intent(in artificial_visc,
type(field_t), intent(in mu,
type(field_t), intent(in kappa 
)

rho = rho - dt * div(m)

m = m - dt * div(rho * u * u^T + p*I)

E = E - dt * div(u * (E + p))

Definition at line 605 of file compressible_res_device.F90.

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

Variable Documentation

◆ compressible_res_device_add_physical_flux

logical compressible_res_device::compressible_res_device_add_physical_flux = .false.

Definition at line 63 of file compressible_res_device.F90.

◆ compressible_res_device_add_physical_stress

logical compressible_res_device::compressible_res_device_add_physical_stress = .false.

Definition at line 65 of file compressible_res_device.F90.

◆ compressible_res_device_gamma

real(kind=rp), public compressible_res_device::compressible_res_device_gamma = 1.4_rp

Definition at line 67 of file compressible_res_device.F90.