source_term_definition_channel Subroutine

public subroutine source_term_definition_channel(me, sim)

Arguments

Type IntentOptional Attributes Name
type(channel_t), intent(inout) :: me
type(simulation_t), intent(in) :: sim

Calls

proc~~source_term_definition_channel~~CallsGraph proc~source_term_definition_channel source_term_definition_channel proc~signal_v1d signal_t%signal_v1d proc~source_term_definition_channel->proc~signal_v1d proc~sourceterms_notfriction_channel SourceTerms_NOTfriction_channel proc~source_term_definition_channel->proc~sourceterms_notfriction_channel proc~he_prop he_prop proc~sourceterms_notfriction_channel->proc~he_prop proc~jacobian_rot jacobian_roT proc~sourceterms_notfriction_channel->proc~jacobian_rot proc~nusselt_correlation nusselt_t%nusselt_correlation proc~sourceterms_notfriction_channel->proc~nusselt_correlation proc~nusselt_correlation_der_pra nusselt_t%nusselt_correlation_der_Pra proc~sourceterms_notfriction_channel->proc~nusselt_correlation_der_pra proc~nusselt_correlation_der_re nusselt_t%nusselt_correlation_der_Re proc~sourceterms_notfriction_channel->proc~nusselt_correlation_der_re proc~nusselt_correlation_der_t nusselt_t%nusselt_correlation_der_T proc~sourceterms_notfriction_channel->proc~nusselt_correlation_der_t proc~state_roe state_roE proc~sourceterms_notfriction_channel->proc~state_roe proc~state_roe_withr state_roE_withR proc~sourceterms_notfriction_channel->proc~state_roe_withr proc~he_prop->proc~jacobian_rot proc~set_error set_error proc~he_prop->proc~set_error proc~eos_terms eos_terms proc~jacobian_rot->proc~eos_terms proc~state_roe->proc~jacobian_rot proc~t_roe T_roE proc~state_roe->proc~t_roe proc~state_roe_withr->proc~jacobian_rot proc~t_rop T_roP proc~state_roe_withr->proc~t_rop proc~eos_e_terms eos_e_terms proc~eos_terms->proc~eos_e_terms proc~fill_f_terms fill_f_terms proc~eos_terms->proc~fill_f_terms proc~brent brent proc~t_roe->proc~brent proc~fill_e_tpart fill_e_Tpart proc~t_roe->proc~fill_e_tpart proc~fill_ff_terms fill_ff_terms proc~t_roe->proc~fill_ff_terms proc~t_rop->proc~brent proc~t_rop->proc~fill_f_terms proc~fill_g_dreg_terms fill_g_Dreg_terms proc~t_rop->proc~fill_g_dreg_terms proc~brent->proc~set_error f f proc~brent->f proc~zero zero proc~brent->proc~zero proc~eos_e_terms->proc~fill_e_tpart proc~eos_e_terms->proc~fill_ff_terms proc~fill_e_tpart->proc~fill_g_dreg_terms proc~zero->f

Called by

proc~~source_term_definition_channel~~CalledByGraph proc~source_term_definition_channel source_term_definition_channel proc~full_physics_definition_channel full_physics_definition_channel proc~full_physics_definition_channel->proc~source_term_definition_channel proc~main_loop main_loop proc~main_loop->proc~source_term_definition_channel proc~main_loop->proc~full_physics_definition_channel program~reims_p reims_p program~reims_p->proc~main_loop

Source Code

subroutine source_term_definition_channel(me,sim)
    type(channel_t), intent(inout) :: me
    type(simulation_t), intent(in) :: sim

    real(dp) :: QorT(me%HeProp%NbCells)
  
        integer :: i
    
        if(sim%explicit) then
          me%big%bmx=0.0_dp; me%big%bvx=0.0_dp
          
          do i=1,me%HeProp%NbCells
            me%big%bmx(1:Nb_VarC,1:Nb_VarC,i)=id_4x4(:,:)            
          enddo
        endif

        if(me%HeProp%ExtHeating=='flux' .or. me%HeProp%ExtHeating=='temp') QorT=me%HeProp%Temp_or_Q_Wall%v1d()

        do i=1,me%HeProp%NbCells
          ! Adding of source term contributions
            call SourceTerms_NOTfriction_channel(me,i,QorT(i))
      
            me%big%bmx(1:Nb_VarC,Con_Qdm,i)=me%big%bmx(1:Nb_VarC,Con_Qdm,i)-sim%dt*(me%srcPP%DerSu(:,i))
            me%big%bmx(1:Nb_VarC,Con_Ene,i)=me%big%bmx(1:Nb_VarC,Con_Ene,i)-sim%dt*(me%srcPP%DerSe(:))
            me%big%bmx(1:Nb_VarC,Con_R,i)=me%big%bmx(1:Nb_VarC,Con_R,i)-sim%dt*(me%srcPP%DerSr(:,i))
      
            me%big%bvx(Con_Qdm,i)=me%big%bvx(Con_Qdm,i)+sim%dt*(me%srcPP%Su(i))
            me%big%bvx(Con_Ene,i)=me%big%bvx(Con_Ene,i)+sim%dt*(me%srcPP%Se)
            me%big%bvx(Con_R,i)=me%big%bvx(Con_R,i)+sim%dt*(me%srcPP%Sr(i))    
        enddo
        
end subroutine source_term_definition_channel