user_module.f90

Go to the documentation of this file.

 
module topology_optimization_user_module
  use case, only: case_t
  use field, only: field_t
  use json_file_module, only: json_file
  use json_utils, only: json_get
  use num_types, only: rp
  use design_module, only: topopt_permeability_force
  use coefs, only: coef_t
  use initial_conditions, only: scalar_z_split_ic
 
  implicit none
 
  private
  public :: neko_user_init
 
contains
 
  subroutine neko_user_init(neko_case)
    type(case_t), intent(inout) :: neko_case
 
    ! Set the properties for the fluid
    neko_case%usr%material_properties => set_material_properties
    neko_case%usr%scalar_user_ic => scalar_z_split_ic
 
  end subroutine neko_user_init
 
  subroutine set_material_properties(t, tstep, rho, mu, cp, lambda, params)
    real(kind=rp), intent(in) :: t
    integer, intent(in) :: tstep
    real(kind=rp), intent(inout) :: rho, mu, cp, lambda
    type(json_file), intent(inout) :: params
 
    real(kind=rp) :: re, pe
 
    call json_get(params, 'case.fluid.Re', re)
    call json_get(params, 'case.scalar.Pe', pe)
 
    mu = 1.0_rp / re
    lambda = 1.0_rp / pe
    rho = 1.0_rp
    cp = 1.0_rp
  end subroutine set_material_properties
 
 
 
 
end module topology_optimization_user_module