fluid_scheme_compressible.f90

Go to the documentation of this file.

! Copyright (c) 2025, The Neko Authors
! All rights reserved.
!
! Redistribution and use in source and binary forms, with or without
! modification, are permitted provided that the following conditions
! are met:
!
!   * Redistributions of source code must retain the above copyright
!     notice, this list of conditions and the following disclaimer.
!
!   * Redistributions in binary form must reproduce the above
!     copyright notice, this list of conditions and the following
!     disclaimer in the documentation and/or other materials provided
!     with the distribution.
!
!   * Neither the name of the authors nor the names of its
!     contributors may be used to endorse or promote products derived
!     from this software without specific prior written permission.
!
! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
! "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
! LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
! FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
! COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
! INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
! BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
! LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
! CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
! LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
! ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
! POSSIBILITY OF SUCH DAMAGE.
!
module fluid_scheme_compressible
  use dirichlet, only : dirichlet_t
  use field, only : field_t
  use field_math, only : field_cfill, field_col2, field_col3, &
       field_cmult2, field_cmult, field_addcol3, field_add2
  use field_registry, only : neko_field_registry
  use fluid_scheme_base, only : fluid_scheme_base_t
  use json_module, only : json_file
  use logger, only : log_size
  use num_types, only : rp
  use mesh, only : mesh_t
  use scratch_registry, only : scratch_registry_t
  use space, only : space_t, gll
  use user_intf, only : user_t
  use usr_inflow, only : usr_inflow_eval
  use json_utils, only : json_get_or_default
  implicit none
  private
 
  type, public, abstract, extends(fluid_scheme_base_t) :: fluid_scheme_compressible_t
     type(field_t), pointer :: m_x => null() 
     type(field_t), pointer :: m_y => null() 
     type(field_t), pointer :: m_z => null() 
     type(field_t), pointer :: e => null() 
 
     real(kind=rp) :: gamma
 
     type(scratch_registry_t) :: scratch
 
   contains
     procedure, pass(this) :: scheme_init => fluid_scheme_compressible_init
     procedure, pass(this) :: scheme_free => fluid_scheme_compressible_free
 
     procedure, pass(this) :: validate => fluid_scheme_compressible_validate
     procedure, pass(this) :: compute_cfl &
          => fluid_scheme_compressible_compute_cfl
     procedure, pass(this) :: update_material_properties => &
          fluid_scheme_compressible_update_material_properties
 
  end type fluid_scheme_compressible_t
  type, public, abstract, extends(fluid_scheme_base_t) :: fluid_scheme_compressible_t …
 
contains
  subroutine fluid_scheme_compressible_init(this, msh, lx, params, scheme, user)
    class(fluid_scheme_compressible_t), target, intent(inout) :: this
    type(mesh_t), target, intent(inout) :: msh
    integer, intent(in) :: lx
    character(len=*), intent(in) :: scheme
    type(json_file), target, intent(inout) :: params
    type(user_t), target, intent(in) :: user
 
    !
    ! SEM simulation fundamentals
    !
 
    this%msh => msh
 
    if (msh%gdim .eq. 2) then
       call this%Xh%init(gll, lx, lx)
    else
       call this%Xh%init(gll, lx, lx, lx)
    end if
 
    call this%dm_Xh%init(msh, this%Xh)
 
    call this%gs_Xh%init(this%dm_Xh)
 
    call this%c_Xh%init(this%gs_Xh)
 
    ! Local scratch registry
    this%scratch = scratch_registry_t(this%dm_Xh, 10, 2)
 
    ! Case parameters
    this%params => params
 
    ! Assign a name
    call json_get_or_default(params, 'case.fluid.name', this%name, "fluid")
 
    ! Fill mu and rho field with the physical value
    call this%mu_field%init(this%dm_Xh, "mu")
    call this%rho_field%init(this%dm_Xh, "rho")
    call field_cfill(this%mu_field, 0.0_rp, this%mu_field%size())
 
    ! Assign momentum fields
    call neko_field_registry%add_field(this%dm_Xh, "m_x")
    call neko_field_registry%add_field(this%dm_Xh, "m_y")
    call neko_field_registry%add_field(this%dm_Xh, "m_z")
    this%m_x => neko_field_registry%get_field("m_x")
    this%m_y => neko_field_registry%get_field("m_y")
    this%m_z => neko_field_registry%get_field("m_z")
    call this%m_x%init(this%dm_Xh, "m_x")
    call this%m_y%init(this%dm_Xh, "m_y")
    call this%m_z%init(this%dm_Xh, "m_z")
 
    ! Assign energy field
    call neko_field_registry%add_field(this%dm_Xh, "E")
    this%E => neko_field_registry%get_field("E")
    call this%E%init(this%dm_Xh, "E")
 
    ! ! Assign velocity fields
    call neko_field_registry%add_field(this%dm_Xh, "u")
    call neko_field_registry%add_field(this%dm_Xh, "v")
    call neko_field_registry%add_field(this%dm_Xh, "w")
    this%u => neko_field_registry%get_field("u")
    this%v => neko_field_registry%get_field("v")
    this%w => neko_field_registry%get_field("w")
    call this%u%init(this%dm_Xh, "u")
    call this%v%init(this%dm_Xh, "v")
    call this%w%init(this%dm_Xh, "w")
    call neko_field_registry%add_field(this%dm_Xh, 'p')
    this%p => neko_field_registry%get_field('p')
    call this%p%init(this%dm_Xh, "p")
 
    ! !! Initialize time-lag fields
    call this%ulag%init(this%u, 1)
    call this%vlag%init(this%v, 1)
    call this%wlag%init(this%w, 1)
 
    !
    ! Setup right-hand side fields.
    !
    allocate(this%f_x)
    allocate(this%f_y)
    allocate(this%f_z)
    call this%f_x%init(this%dm_Xh, fld_name = "fluid_rhs_x")
    call this%f_y%init(this%dm_Xh, fld_name = "fluid_rhs_y")
    call this%f_z%init(this%dm_Xh, fld_name = "fluid_rhs_z")
 
    ! Compressible parameters
    call json_get_or_default(params, 'case.fluid.gamma', this%gamma, 1.4_rp)
  subroutine fluid_scheme_compressible_init(this, msh, lx, params, scheme, user) …
  end subroutine fluid_scheme_compressible_init
 
  subroutine fluid_scheme_compressible_free(this)
    class(fluid_scheme_compressible_t), intent(inout) :: this
    call this%dm_Xh%free()
    call this%gs_Xh%free()
    call this%c_Xh%free()
    call this%Xh%free()
 
    call this%mu_field%free()
    call this%rho_field%free()
    call this%m_x%free()
    call this%m_y%free()
    call this%m_z%free()
    call this%E%free()
 
    nullify(this%m_x)
    nullify(this%m_y)
    nullify(this%m_z)
 
    nullify(this%E)
 
    nullify(this%u)
    nullify(this%v)
    nullify(this%w)
    nullify(this%p)
 
    call this%ulag%free()
    call this%vlag%free()
    call this%wlag%free()
  subroutine fluid_scheme_compressible_free(this) …
  end subroutine fluid_scheme_compressible_free
 
  subroutine fluid_scheme_compressible_validate(this)
    class(fluid_scheme_compressible_t), target, intent(inout) :: this
    integer :: n
    type(field_t), pointer :: temp
    integer :: temp_indices(1)
 
    n = this%dm_Xh%size()
    call this%scratch%request_field(temp, temp_indices(1))
 
    call field_col3(this%m_x, this%u, this%rho_field)
    call field_col3(this%m_y, this%v, this%rho_field)
    call field_col3(this%m_z, this%w, this%rho_field)
 
    call field_cmult2(this%E, this%p, 1.0_rp/(this%gamma - 1.0_rp), n)
    call field_col3(temp, this%u, this%u, n)
    call field_addcol3(temp, this%v, this%v, n)
    call field_addcol3(temp, this%w, this%w, n)
    call field_col2(temp, this%rho_field, n)
    call field_cmult(temp, 0.5_rp, n)
    call field_add2(this%E, temp, n)
 
    call this%scratch%relinquish_field(temp_indices)
 
  subroutine fluid_scheme_compressible_validate(this) …
  end subroutine fluid_scheme_compressible_validate
 
  function fluid_scheme_compressible_compute_cfl(this, dt) result(c)
    class(fluid_scheme_compressible_t), intent(in) :: this
    real(kind=rp), intent(in) :: dt
    real(kind=rp) :: c
 
    c = 0.1_rp
 
  function fluid_scheme_compressible_compute_cfl(this, dt) result(c) …
  end function fluid_scheme_compressible_compute_cfl
 
  subroutine fluid_scheme_compressible_update_material_properties(this)
    class(fluid_scheme_compressible_t), intent(inout) :: this
  subroutine fluid_scheme_compressible_update_material_properties(this) …
  end subroutine fluid_scheme_compressible_update_material_properties
end module fluid_scheme_compressible