MFC
Exascale flow solver
Loading...
Searching...
No Matches
m_conduction Module Reference

Fourier heat conduction, div(k grad T), as a face-centered source flux on the energy equation. The term has no cross-derivatives, so the direction-split face difference below is exact for it. More...

Functions/Subroutines

subroutine, public s_compute_conduction_source_flux (idir, q_prim_qp, q_t_sf, flux_src_vf, irx, iry, irz)
 Accumulate -k*dT/dx_idir into the energy source flux at each idir-normal face.
subroutine, public s_compute_conduction_axis_source (q_prim_vf, q_t_sf, tau_re_vf, ix, iy, iz)
 Cell-centered -k*dT/dr for the axis cell of a cylindrical grid. The generic geometric source in m_rhs uses face values of flux_src(E) in this role; the axis cell has no face pair, so it reads tau_Re_vf(E) instead. Accumulates: call after s_compute_viscous_stress_cylindrical_boundary, which zeroes tau_Re_vf(mombeg:E); without viscosity this routine does that zeroing itself.

Variables

type(int_bounds_infoisc1
type(int_bounds_infoisc2
type(int_bounds_infoisc3
integer, dimension(3) offsets_c

Detailed Description

Fourier heat conduction, div(k grad T), as a face-centered source flux on the energy equation. The term has no cross-derivatives, so the direction-split face difference below is exact for it.

Function/Subroutine Documentation

◆ s_compute_conduction_axis_source()

subroutine, public m_conduction::s_compute_conduction_axis_source ( type(scalar_field), dimension(sys_size), intent(in) q_prim_vf,
type(scalar_field), intent(in) q_t_sf,
type(scalar_field), dimension(1:sys_size), intent(inout) tau_re_vf,
type(int_bounds_info), intent(in) ix,
type(int_bounds_info), intent(in) iy,
type(int_bounds_info), intent(in) iz )

Cell-centered -k*dT/dr for the axis cell of a cylindrical grid. The generic geometric source in m_rhs uses face values of flux_src(E) in this role; the axis cell has no face pair, so it reads tau_Re_vf(E) instead. Accumulates: call after s_compute_viscous_stress_cylindrical_boundary, which zeroes tau_Re_vf(mombeg:E); without viscosity this routine does that zeroing itself.

Definition at line 493 of file m_conduction.fpp.f90.

◆ s_compute_conduction_source_flux()

subroutine, public m_conduction::s_compute_conduction_source_flux ( integer, intent(in) idir,
type(scalar_field), dimension(sys_size), intent(in) q_prim_qp,
type(scalar_field), intent(in) q_t_sf,
type(scalar_field), dimension(sys_size), intent(inout) flux_src_vf,
type(int_bounds_info), intent(in) irx,
type(int_bounds_info), intent(in) iry,
type(int_bounds_info), intent(in) irz )

Accumulate -k*dT/dx_idir into the energy source flux at each idir-normal face.

Definition at line 383 of file m_conduction.fpp.f90.

Variable Documentation

◆ isc1

type(int_bounds_info) m_conduction::isc1

Definition at line 355 of file m_conduction.fpp.f90.

◆ isc2

type(int_bounds_info) m_conduction::isc2

Definition at line 355 of file m_conduction.fpp.f90.

◆ isc3

type(int_bounds_info) m_conduction::isc3

Definition at line 355 of file m_conduction.fpp.f90.

◆ offsets_c

integer, dimension(3) m_conduction::offsets_c

Definition at line 367 of file m_conduction.fpp.f90.