|
MFC
Exascale flow solver
|
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_info) | isc1 |
| type(int_bounds_info) | isc2 |
| type(int_bounds_info) | isc3 |
| integer, dimension(3) | offsets_c |
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.
| 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.
| 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.
| type(int_bounds_info) m_conduction::isc1 |
Definition at line 355 of file m_conduction.fpp.f90.
| type(int_bounds_info) m_conduction::isc2 |
Definition at line 355 of file m_conduction.fpp.f90.
| type(int_bounds_info) m_conduction::isc3 |
Definition at line 355 of file m_conduction.fpp.f90.
| integer, dimension(3) m_conduction::offsets_c |
Definition at line 367 of file m_conduction.fpp.f90.