38 use json_module,
only : json_file
57 real(kind=
rp) :: kappa = 0.41_rp
59 real(kind=
rp) :: b = 0.0_rp
61 real(kind=
rp) :: z0 = 0.0_rp
71 procedure, pass(this) :: init_from_components => &
90 character(len=*),
intent(in) :: scheme_name
91 type(
coef_t),
intent(in) :: coef
92 integer,
intent(in) :: msk(:)
93 integer,
intent(in) :: facet(:)
94 integer,
intent(in) :: h_index
95 type(json_file),
intent(inout) :: json
96 real(kind=
rp) :: kappa, b, z0
102 call this%init_from_components(scheme_name, coef, msk, facet, h_index, &
111 type(
coef_t),
intent(in) :: coef
112 type(json_file),
intent(inout) :: json
114 call this%partial_init_base(coef, json)
126 integer,
intent(in) :: msk(:)
127 integer,
intent(in) :: facet(:)
129 call this%finalize_base(msk, facet)
143 facet, h_index, kappa, B, z0)
145 character(len=*),
intent(in) :: scheme_name
146 type(
coef_t),
intent(in) :: coef
147 integer,
intent(in) :: msk(:)
148 integer,
intent(in) :: facet(:)
149 integer,
intent(in) :: h_index
150 real(kind=
rp),
intent(in) :: kappa
151 real(kind=
rp),
intent(in) :: b
152 real(kind=
rp),
intent(in) :: z0
154 call this%init_base(scheme_name, coef, msk, facet, h_index)
166 call this%free_base()
175 real(kind=
rp),
intent(in) :: t
176 integer,
intent(in) :: tstep
181 real(kind=
rp) :: ui, vi, wi, magu, utau, normu
189 this%ind_s_d, this%ind_t_d, this%ind_e_d, &
190 this%n_x%x_d, this%n_y%x_d, this%n_z%x_d, &
191 this%h%x_d, this%tau_x%x_d, this%tau_y%x_d, &
192 this%tau_z%x_d, this%n_nodes, u%Xh%lx, this%kappa, &
193 this%B, this%z0, tstep)
196 this%ind_t, this%ind_e, this%n_x%x, this%n_y%x, this%n_z%x, &
197 this%h%x, this%tau_x%x, this%tau_y%x, this%tau_z%x, &
198 this%n_nodes, u%Xh%lx, u%msh%nelv, this%kappa, &
199 this%B, this%z0, tstep)
__global__ void rough_log_law_compute(const T *__restrict__ u_d, const T *__restrict__ v_d, const T *__restrict__ w_d, const int *__restrict__ ind_r_d, const int *__restrict__ ind_s_d, const int *__restrict__ ind_t_d, const int *__restrict__ ind_e_d, const T *__restrict__ n_x_d, const T *__restrict__ n_y_d, const T *__restrict__ n_z_d, const T *__restrict__ h_d, T *__restrict__ tau_x_d, T *__restrict__ tau_y_d, T *__restrict__ tau_z_d, const int n_nodes, const int lx, const T kappa, const T B, const T z0)
Retrieves a parameter by name or assigns a provided default value. In the latter case also adds the m...
Retrieves a parameter by name or throws an error.
Defines a registry for storing solution fields.
type(field_registry_t), target, public neko_field_registry
Global field registry.
Utilities for retrieving parameters from the case files.
integer, parameter neko_bcknd_device
integer, parameter, public rp
Global precision used in computations.
Implements the CPU kernel for the rough_log_law_t type.
subroutine, public rough_log_law_compute_cpu(u, v, w, ind_r, ind_s, ind_t, ind_e, n_x, n_y, n_z, h, tau_x, tau_y, tau_z, n_nodes, lx, nelv, kappa, b, z0, tstep)
Compute the wall shear stress on CPU using the rough log-law model.
Implements the device kernel for the rough_log_law_t type.
subroutine, public rough_log_law_compute_device(u_d, v_d, w_d, ind_r_d, ind_s_d, ind_t_d, ind_e_d, n_x_d, n_y_d, n_z_d, h_d, tau_x_d, tau_y_d, tau_z_d, n_nodes, lx, kappa, b, z0, tstep)
Compute the wall shear stress on device using the rough log-law model.
Implements rough_log_law_t.
subroutine rough_log_law_init_from_components(this, scheme_name, coef, msk, facet, h_index, kappa, b, z0)
Constructor from components.
subroutine rough_log_law_finalize(this, msk, facet)
Finalize the construction using the mask and facet arrays of the bc.
subroutine rough_log_law_init(this, scheme_name, coef, msk, facet, h_index, json)
Constructor from JSON.
subroutine rough_log_law_partial_init(this, coef, json)
Constructor from JSON.
subroutine rough_log_law_free(this)
Destructor for the rough_log_law_t (base) class.
Defines a registry for storing and requesting temporary fields This can be used when you have a funct...
type(scratch_registry_t), target, public neko_scratch_registry
Global scratch registry.
Coefficients defined on a given (mesh, ) tuple. Arrays use indices (i,j,k,e): element e,...
Wall model based on the log-law for a rough wall. The formula defining the law is ....
Base abstract type for wall-stress models for wall-modelled LES.