36 use num_types,
only: rp
37 use field_list,
only: field_list_t
38 use json_module,
only: json_file
39 use json_utils,
only: json_get, json_get_or_default
40 use source_term,
only: source_term_t
41 use coefs,
only: coef_t
42 use neko_config,
only: neko_bcknd_device
43 use time_state,
only: time_state_t
44 use utils,
only: neko_error
45 use field,
only: field_t
46 use field_math,
only: field_subcol3
54 type(field_t),
pointer :: chi => null()
56 type(field_t),
pointer :: u => null()
58 type(field_t),
pointer :: v => null()
60 type(field_t),
pointer :: w => null()
64 procedure, pass(this) :: init => &
65 simple_brinkman_source_term_init_from_json
67 procedure, pass(this) :: init_from_components => &
68 simple_brinkman_source_term_init_from_components
70 procedure, pass(this) :: free => simple_brinkman_source_term_free
72 procedure, pass(this) :: compute_ => simple_brinkman_source_term_compute
82 subroutine simple_brinkman_source_term_init_from_json(this, json, fields, &
85 type(json_file),
intent(inout) :: json
86 type(field_list_t),
intent(in),
target :: fields
87 type(coef_t),
intent(in),
target :: coef
88 character(len=*),
intent(in) :: variable_name
95 end subroutine simple_brinkman_source_term_init_from_json
103 subroutine simple_brinkman_source_term_init_from_components(this, &
104 f_x, f_y, f_z, chi, u, v, w, coef)
106 type(field_t),
pointer,
intent(in) :: f_x, f_y, f_z
107 type(field_list_t) :: fields
109 real(kind=rp) :: start_time
110 real(kind=rp) :: end_time
111 type(field_t),
intent(in),
target :: u, v, w
112 type(field_t),
intent(in),
target :: chi
117 end_time = 100000000.0_rp
124 call fields%assign(1, f_x)
125 call fields%assign(2, f_y)
126 call fields%assign(3, f_z)
128 call this%init_base(fields, coef, start_time, end_time)
137 end subroutine simple_brinkman_source_term_init_from_components
140 subroutine simple_brinkman_source_term_free(this)
143 call this%free_base()
144 end subroutine simple_brinkman_source_term_free
149 subroutine simple_brinkman_source_term_compute(this, time)
151 type(time_state_t),
intent(in) :: time
152 type(field_t),
pointer :: fu, fv, fw
154 fu => this%fields%get_by_index(1)
155 fv => this%fields%get_by_index(2)
156 fw => this%fields%get_by_index(3)
158 call field_subcol3(fu, this%u, this%chi)
159 call field_subcol3(fv, this%v, this%chi)
160 call field_subcol3(fw, this%w, this%chi)
162 end subroutine simple_brinkman_source_term_compute
Implements the simple_brinkman_source_term_t type.
A simple Brinkman source term.