sub_navier_stokes.f90

Go to the documentation of this file.

MODULE update_navier_stokes
PUBLIC :: navier_stokes_decouple
PRIVATE

CONTAINS

  SUBROUTINE  navier_stokes_decouple(comm_one_d_ns,time, vv_3_LA, pp_1_LA, &
       list_mode, pp_mesh, vv_mesh, incpn_m1, incpn, pn_m1, pn, un_m1, un, &
       vvz_per, pp_per, hn_p2, bn_p2, density_m2, density_m1, density, visco_dyn, tempn, &
       level_set_m1, level_set, visc_entro_level)

    USE def_type_mesh
    USE periodic
    USE input_data
    USE subroutine_ns_with_u
    USE subroutine_ns_with_m
    USE my_util
    IMPLICIT NONE
    TYPE(mesh_type),                INTENT(IN)     :: vv_mesh, pp_mesh
    INTEGER,      DIMENSION(:),     INTENT(IN)     :: list_mode
    REAL(KIND=8), DIMENSION(vv_mesh%np,6,SIZE(list_mode)), INTENT(INOUT)  :: un_m1, un
    REAL(KIND=8), DIMENSION(vv_mesh%np,2,SIZE(list_mode)), INTENT(IN)     :: density_m2, density_m1, density
    REAL(KIND=8), DIMENSION(vv_mesh%np,2,SIZE(list_mode)), INTENT(IN)     :: visco_dyn
    REAL(KIND=8), DIMENSION(:,:,:),                        INTENT(IN)     :: tempn
    REAL(KIND=8), DIMENSION(pp_mesh%np,2,SIZE(list_mode)), INTENT(INOUT)  :: pn_m1, pn, incpn_m1, incpn
    REAL(KIND=8), DIMENSION(:,:,:),                        INTENT(IN)     :: hn_p2
    REAL(KIND=8), DIMENSION(:,:,:),                        INTENT(IN)     :: bn_p2
    REAL(KIND=8), DIMENSION(:,:,:,:),                      INTENT(IN)     :: level_set_m1, level_set
    REAL(KIND=8),                   INTENT(IN)     :: time  
    REAL(KIND=8), DIMENSION(:,:),   INTENT(OUT)    :: visc_entro_level
    TYPE(periodic_type),            INTENT(IN)     :: vvz_per, pp_per
    TYPE(petsc_csr_la)                             :: vv_3_la, pp_1_la
#include "petsc/finclude/petsc.h"
    mpi_comm, DIMENSION(:), POINTER  :: comm_one_d_ns

    IF (inputs%if_navier_stokes_with_u) THEN
       CALL bdf2_ns_stress_bc_with_u(comm_one_d_ns,time,vv_3_la, pp_1_la, vvz_per, pp_per, &
            inputs%dt, inputs%Re, list_mode, pp_mesh, vv_mesh, incpn_m1, incpn, &
            pn_m1, pn, un_m1, un, hn_p2, bn_p2, opt_tempn=tempn)
    ELSE
       CALL three_level_ns_tensor_sym_with_m(comm_one_d_ns,time,vv_3_la, pp_1_la, &
            inputs%dt, inputs%Re, list_mode, pp_mesh, vv_mesh, incpn_m1, incpn, &
            pn_m1, pn, un_m1, un, hn_p2, bn_p2, density_m2, density_m1, density, visco_dyn, tempn, &
            level_set_m1, level_set, visc_entro_level)
    END IF
  END SUBROUTINE navier_stokes_decouple

END MODULE update_navier_stokes