TKE SGS model

src/.depends

@@ -150,7 +150,9 @@ io/fluid_output.lo : io/fluid_output.f90 field/field.lo registries/registry.lo s
 io/fld_file_output.lo : io/fld_file_output.f90 io/output.lo device/device.lo config/neko_config.lo field/field_list.lo config/num_types.lo 
 io/chkp_output.lo : io/chkp_output.f90 config/num_types.lo io/output.lo common/checkpoint.lo 
 io/fluid_stats_output.lo : io/fluid_stats_output.f90 math/matrix.lo io/output.lo device/device.lo io/fld_file_data.lo sem/map_2d.lo sem/map_1d.lo config/num_types.lo config/neko_config.lo fluid/fluid_stats.lo 
+io/fluid_sgs_stats_output.lo : io/fluid_sgs_stats_output.f90 math/matrix.lo io/output.lo device/device.lo io/fld_file_data.lo sem/map_2d.lo sem/map_1d.lo config/num_types.lo config/neko_config.lo les/fluid_sgs_stats.lo 
 io/scalar_stats_output.lo : io/scalar_stats_output.f90 math/matrix.lo io/output.lo device/device.lo io/fld_file_data.lo sem/map_2d.lo sem/map_1d.lo config/num_types.lo config/neko_config.lo scalar/scalar_stats.lo 
+io/scalar_sgs_stats_output.lo : io/scalar_sgs_stats_output.f90 math/matrix.lo io/output.lo device/device.lo io/fld_file_data.lo sem/map_2d.lo sem/map_1d.lo config/num_types.lo config/neko_config.lo les/scalar_sgs_stats.lo 
 io/data_streamer.lo : io/data_streamer.F90 comm/comm.lo common/utils.lo sem/coef.lo field/field.lo config/num_types.lo 
 io/output_controller.lo : io/output_controller.f90 common/time_based_controller.lo config/num_types.lo common/profiler.lo common/utils.lo common/log.lo common/time_state.lo comm/comm.lo io/fld_file.lo io/output.lo 
 global_interpolation/global_interpolation.lo : global_interpolation/global_interpolation.f90 common/structs.lo math/math.lo math/matrix.lo math/vector.lo global_interpolation/global_interpolation_comm.lo comm/comm.lo global_interpolation/pe_finder.lo global_interpolation/el_finder.lo global_interpolation/legendre_rst_finder.lo global_interpolation/cartesian_pe_finder.lo global_interpolation/cartesian_el_finder.lo global_interpolation/aabb_el_finder.lo global_interpolation/aabb_pe_finder.lo device/device.lo sem/local_interpolation.lo common/utils.lo common/log.lo sem/dofmap.lo adt/stack.lo sem/space.lo config/neko_config.lo config/num_types.lo 
@@ -220,7 +222,9 @@ common/rhs_maker_fctry.lo : common/rhs_maker_fctry.f90 config/neko_config.lo com
 simulation_components/probes.lo : simulation_components/probes.F90 device/device.lo config/neko_config.lo comm/comm.lo case.lo io/csv_file.lo io/file.lo mesh/point_zone_registry.lo mesh/point_zone.lo math/tensor.lo global_interpolation/global_interpolation.lo common/json_utils.lo sem/dofmap.lo registries/registry.lo simulation_components/simulation_component.lo common/time_state.lo field/field_list.lo common/utils.lo common/log.lo math/matrix.lo config/num_types.lo 
 simulation_components/force_torque.lo : simulation_components/force_torque.f90 device/device.lo comm/comm.lo math/bcknd/device/device_math.lo math/math.lo config/neko_config.lo common/log.lo qoi/drag_torque.lo bc/dirichlet.lo math/vector.lo sem/coef.lo common/json_utils.lo case.lo math/operators.lo field/field.lo common/time_state.lo registries/scratch_registry.lo registries/registry.lo simulation_components/simulation_component.lo common/time_based_controller.lo config/num_types.lo 
 simulation_components/fluid_stats_simcomp.lo : simulation_components/fluid_stats_simcomp.f90 comm/comm.lo common/json_utils.lo common/log.lo common/utils.lo sem/coef.lo case.lo io/fluid_stats_output.lo fluid/fluid_stats.lo field/field.lo common/time_state.lo registries/registry.lo simulation_components/simulation_component.lo config/num_types.lo 
+simulation_components/fluid_sgs_stats_simcomp.lo : simulation_components/fluid_sgs_stats_simcomp.f90 comm/comm.lo common/json_utils.lo common/log.lo common/utils.lo sem/coef.lo case.lo io/fluid_sgs_stats_output.lo les/fluid_sgs_stats.lo field/field.lo common/time_state.lo registries/registry.lo simulation_components/simulation_component.lo config/num_types.lo 
 simulation_components/scalar_stats_simcomp.lo : simulation_components/scalar_stats_simcomp.f90 comm/comm.lo common/json_utils.lo common/log.lo common/utils.lo sem/coef.lo case.lo io/scalar_stats_output.lo scalar/scalar_stats.lo field/field.lo common/time_state.lo registries/registry.lo simulation_components/simulation_component.lo config/num_types.lo 
+simulation_components/scalar_sgs_stats_simcomp.lo : simulation_components/scalar_sgs_stats_simcomp.f90 comm/comm.lo common/json_utils.lo common/log.lo common/utils.lo sem/coef.lo case.lo io/scalar_sgs_stats_output.lo les/scalar_sgs_stats.lo field/field.lo common/time_state.lo registries/registry.lo simulation_components/simulation_component.lo config/num_types.lo 
 simulation_components/user_stats.lo : simulation_components/user_stats.f90 common/time_based_controller.lo common/time_state.lo sem/coef.lo field/mean_field.lo common/json_utils.lo io/mean_field_output.lo case.lo field/field.lo registries/registry.lo simulation_components/simulation_component.lo config/num_types.lo 
 simulation_components/field_writer.lo : simulation_components/field_writer.f90 common/time_based_controller.lo common/json_utils.lo io/fld_file_output.lo case.lo registries/registry.lo common/time_state.lo simulation_components/simulation_component.lo config/num_types.lo 
 common/bcknd/cpu/rhs_maker_cpu.lo : common/bcknd/cpu/rhs_maker_cpu.f90 registries/scratch_registry.lo config/num_types.lo field/field.lo field/field_series.lo common/rhs_maker.lo 
@@ -334,7 +338,7 @@ simulation_components/weak_gradient_simcomp.lo : simulation_components/weak_grad
 simulation_components/gradient_simcomp.lo : simulation_components/gradient_simcomp.f90 common/time_based_controller.lo simulation_components/field_writer.lo common/json_utils.lo case.lo common/time_state.lo math/operators.lo field/field.lo registries/registry.lo simulation_components/simulation_component.lo config/num_types.lo 
 simulation_components/derivative_simcomp.lo : simulation_components/derivative_simcomp.f90 common/time_based_controller.lo common/utils.lo simulation_components/field_writer.lo common/json_utils.lo case.lo math/operators.lo field/field.lo registries/registry.lo common/time_state.lo simulation_components/simulation_component.lo config/num_types.lo 
 simulation_components/les_simcomp.lo : simulation_components/les_simcomp.f90 common/utils.lo simulation_components/field_writer.lo common/json_utils.lo les/les_model.lo common/time_state.lo case.lo simulation_components/simulation_component.lo config/num_types.lo 
-simulation_components/simulation_component_fctry.lo : simulation_components/simulation_component_fctry.f90 simulation_components/spectral_error.lo simulation_components/derivative_simcomp.lo simulation_components/divergence_simcomp.lo simulation_components/gradient_simcomp.lo simulation_components/weak_gradient_simcomp.lo simulation_components/curl_simcomp.lo simulation_components/field_writer.lo common/utils.lo simulation_components/les_simcomp.lo simulation_components/probes.lo simulation_components/lambda2.lo simulation_components/user_stats.lo simulation_components/scalar_stats_simcomp.lo simulation_components/fluid_stats_simcomp.lo simulation_components/force_torque.lo simulation_components/simulation_component.lo 
+simulation_components/simulation_component_fctry.lo : simulation_components/simulation_component_fctry.f90 simulation_components/spectral_error.lo simulation_components/derivative_simcomp.lo simulation_components/divergence_simcomp.lo simulation_components/gradient_simcomp.lo simulation_components/weak_gradient_simcomp.lo simulation_components/curl_simcomp.lo simulation_components/field_writer.lo common/utils.lo simulation_components/les_simcomp.lo simulation_components/probes.lo simulation_components/lambda2.lo simulation_components/user_stats.lo simulation_components/scalar_sgs_stats_simcomp.lo simulation_components/scalar_stats_simcomp.lo simulation_components/fluid_sgs_stats_simcomp.lo simulation_components/fluid_stats_simcomp.lo simulation_components/force_torque.lo simulation_components/simulation_component.lo 
 source_terms/source_term.lo : source_terms/source_term.f90 common/time_state.lo field/field_list.lo sem/coef.lo config/num_types.lo 
 source_terms/user_source_term.lo : source_terms/user_source_term.f90 common/time_state.lo common/user_intf.lo sem/dofmap.lo math/field_math.lo math/bcknd/device/device_math.lo device/device.lo sem/coef.lo field/field_list.lo source_terms/source_term.lo common/utils.lo config/num_types.lo config/neko_config.lo 
 source_terms/coriolis_source_term.lo : source_terms/coriolis_source_term.f90 common/time_state.lo registries/registry.lo field/field.lo source_terms/bcknd/device/coriolis_source_term_device.lo source_terms/bcknd/cpu/coriolis_source_term_cpu.lo common/utils.lo config/neko_config.lo sem/coef.lo source_terms/source_term.lo common/json_utils.lo field/field_list.lo config/num_types.lo 
@@ -366,7 +370,9 @@ filter/bcknd/device/opencl/opencl_mappings.lo : filter/bcknd/device/opencl/openc
 source_terms/gradient_jump_penalty.lo : source_terms/gradient_jump_penalty.f90 math/operators.lo common/time_state.lo registries/registry.lo field/field_list.lo source_terms/source_term.lo source_terms/bcknd/device/device_gradient_jump_penalty.lo math/bcknd/device/device_math.lo device/device.lo gs/gather_scatter.lo sem/space.lo gs/gs_ops.lo mesh/quad.lo mesh/hex.lo mesh/element.lo sem/coef.lo config/neko_config.lo sem/dofmap.lo field/field.lo mesh/point.lo math/math.lo common/json_utils.lo common/utils.lo config/num_types.lo 
 source_terms/bcknd/device/device_gradient_jump_penalty.lo : source_terms/bcknd/device/device_gradient_jump_penalty.F90 device/device.lo common/utils.lo config/num_types.lo 
 les/les_model.lo : les/les_model.f90 comm/comm.lo common/utils.lo math/bcknd/device/device_math.lo math/math.lo device/device.lo config/neko_config.lo gs/gs_ops.lo sem/coef.lo sem/dofmap.lo registries/registry.lo field/field_series.lo field/field.lo common/rhs_maker.lo time_schemes/time_scheme_controller.lo fluid/fluid_pnpn.lo fluid/fluid_scheme_base.lo config/num_types.lo 
-les/les_model_fctry.lo : les/les_model_fctry.f90 common/utils.lo les/wale.lo fluid/fluid_scheme_base.lo les/sigma.lo les/dynamic_smagorinsky.lo les/smagorinsky.lo les/vreman.lo les/les_model.lo 
+les/les_model_fctry.lo : les/les_model_fctry.f90 common/utils.lo les/TKE_SGS.lo les/wale.lo fluid/fluid_scheme_base.lo les/sigma.lo les/dynamic_smagorinsky.lo les/smagorinsky.lo les/vreman.lo les/les_model.lo 
+les/fluid_sgs_stats.lo : les/fluid_sgs_stats.f90 registries/registry.lo common/stats_quant.lo field/field_list.lo field/field.lo sem/coef.lo math/operators.lo math/field_math.lo config/num_types.lo field/mean_field.lo 
+les/scalar_sgs_stats.lo : les/scalar_sgs_stats.f90 registries/registry.lo common/stats_quant.lo field/field_list.lo field/field.lo sem/coef.lo math/operators.lo math/field_math.lo config/num_types.lo field/mean_field.lo 
 les/vreman.lo : les/vreman.f90 common/log.lo registries/registry.lo les/bcknd/device/vreman_device.lo les/bcknd/cpu/vreman_cpu.lo config/neko_config.lo common/json_utils.lo fluid/fluid_scheme_base.lo field/field.lo les/les_model.lo config/num_types.lo 
 les/bcknd/cpu/vreman_cpu.lo : les/bcknd/cpu/vreman_cpu.f90 gs/gs_ops.lo sem/coef.lo math/operators.lo field/field.lo registries/registry.lo registries/scratch_registry.lo math/math.lo field/field_list.lo config/num_types.lo 
 les/bcknd/device/vreman_device.lo : les/bcknd/device/vreman_device.f90 les/bcknd/device/device_vreman_nut.lo math/bcknd/device/device_math.lo gs/gs_ops.lo sem/coef.lo math/operators.lo field/field.lo registries/registry.lo registries/scratch_registry.lo math/math.lo config/num_types.lo 
@@ -387,6 +393,10 @@ les/wale.lo : les/wale.f90 common/log.lo registries/registry.lo les/bcknd/device
 les/bcknd/cpu/wale_cpu.lo : les/bcknd/cpu/wale_cpu.f90 math/math.lo gs/gs_ops.lo sem/coef.lo math/operators.lo field/field.lo registries/registry.lo registries/scratch_registry.lo field/field_list.lo config/num_types.lo 
 les/bcknd/device/wale_device.lo : les/bcknd/device/wale_device.f90 les/bcknd/device/device_wale_nut.lo math/bcknd/device/device_math.lo gs/gs_ops.lo sem/coef.lo math/operators.lo field/field.lo registries/registry.lo registries/scratch_registry.lo math/math.lo config/num_types.lo 
 les/bcknd/device/device_wale_nut.lo : les/bcknd/device/device_wale_nut.F90 comm/comm.lo common/utils.lo config/num_types.lo 
+les/TKE_SGS.lo : les/TKE_SGS.f90 common/log.lo registries/registry.lo les/bcknd/device/TKE_SGS_device.lo les/bcknd/cpu/TKE_SGS_cpu.lo config/neko_config.lo common/json_utils.lo les/les_model.lo fluid/fluid_scheme_base.lo field/field.lo config/num_types.lo common/utils.lo 
+les/bcknd/cpu/TKE_SGS_cpu.lo : les/bcknd/cpu/TKE_SGS_cpu.f90 math/math.lo gs/gs_ops.lo sem/coef.lo math/operators.lo field/field.lo registries/registry.lo registries/scratch_registry.lo field/field_list.lo config/num_types.lo common/utils.lo 
+les/bcknd/device/TKE_SGS_device.lo : les/bcknd/device/TKE_SGS_device.f90 les/bcknd/device/device_TKE_SGS_nut.lo math/bcknd/device/device_math.lo gs/gs_ops.lo sem/coef.lo math/operators.lo field/field.lo registries/registry.lo registries/scratch_registry.lo math/math.lo common/utils.lo config/num_types.lo 
+les/bcknd/device/device_TKE_SGS_nut.lo : les/bcknd/device/device_TKE_SGS_nut.F90 comm/comm.lo common/utils.lo config/num_types.lo 
 wall_models/wall_model.lo : wall_models/wall_model.f90 common/json_utils.lo wall_models/bcknd/device/wall_model_device.lo io/file.lo common/log.lo comm/comm.lo math/math.lo common/utils.lo math/vector.lo device/device.lo config/neko_config.lo sem/coef.lo sem/dofmap.lo registries/registry.lo field/field.lo config/num_types.lo 
 wall_models/bcknd/device/wall_model_device.lo : wall_models/bcknd/device/wall_model_device.F90 common/utils.lo config/num_types.lo 
 wall_models/rough_log_law.lo : wall_models/rough_log_law.f90 registries/scratch_registry.lo wall_models/bcknd/cpu/rough_log_law_cpu.lo wall_models/bcknd/device/rough_log_law_device.lo common/json_utils.lo registries/registry.lo wall_models/wall_model.lo config/neko_config.lo sem/coef.lo config/num_types.lo field/field.lo 

src/.depends_device

@@ -45,6 +45,7 @@ les/bcknd/device/hip/smagorinsky_nut.lo : les/bcknd/device/hip/smagorinsky_nut.h
 les/bcknd/device/hip/vreman_nut.lo : les/bcknd/device/hip/vreman_nut.hip les/bcknd/device/hip/vreman_nut_kernel.h
 les/bcknd/device/hip/sigma_nut.lo : les/bcknd/device/hip/sigma_nut.hip les/bcknd/device/hip/sigma_nut_kernel.h
 les/bcknd/device/hip/wale_nut.lo : les/bcknd/device/hip/wale_nut.hip les/bcknd/device/hip/wale_nut_kernel.h
+les/bcknd/device/hip/TKE_SGS_nut.lo : les/bcknd/device/hip/TKE_SGS_nut.hip les/bcknd/device/hip/TKE_SGS_nut_kernel.h
 multigrid/bcknd/device/hip/amg_cheby.lo : multigrid/bcknd/device/hip/amg_cheby.hip
 krylov/bcknd/device/hip/cheby_aux.lo : krylov/bcknd/device/hip/cheby_aux.hip
 comm/bcknd/device/cuda/nvshmem.lo : comm/bcknd/device/cuda/nvshmem.cu comm/comm.h
@@ -97,6 +98,7 @@ les/bcknd/device/cuda/smagorinsky_nut.lo : les/bcknd/device/cuda/smagorinsky_nut
 les/bcknd/device/cuda/vreman_nut.lo : les/bcknd/device/cuda/vreman_nut.cu les/bcknd/device/cuda/vreman_nut_kernel.h
 les/bcknd/device/cuda/sigma_nut.lo : les/bcknd/device/cuda/sigma_nut.cu les/bcknd/device/cuda/sigma_nut_kernel.h
 les/bcknd/device/cuda/wale_nut.lo : les/bcknd/device/cuda/wale_nut.cu les/bcknd/device/cuda/wale_nut_kernel.h
+les/bcknd/device/cuda/TKE_SGS_nut.lo : les/bcknd/device/cuda/TKE_SGS_nut.cu les/bcknd/device/cuda/TKE_SGS_nut_kernel.h
 device/opencl/check.lo : device/opencl/check.c device/opencl/check.h
 device/opencl/jit.lo : device/opencl/jit.c device/opencl/jit.h
 math/bcknd/device/opencl/math.lo : math/bcknd/device/opencl/math.c math/bcknd/device/opencl/math_kernel.cl.h

src/Makefile.am

@@ -153,7 +153,9 @@ neko_fortran_SOURCES = \
 	io/fld_file_output.f90\
 	io/chkp_output.f90\
 	io/fluid_stats_output.f90\
+	io/fluid_sgs_stats_output.f90\
 	io/scalar_stats_output.f90\
+	io/scalar_sgs_stats_output.f90\
 	io/data_streamer.F90\
 	io/output_controller.f90\
 	global_interpolation/global_interpolation.f90\
@@ -223,7 +225,9 @@ neko_fortran_SOURCES = \
 	simulation_components/probes.F90\
 	simulation_components/force_torque.f90\
 	simulation_components/fluid_stats_simcomp.f90\
+	simulation_components/fluid_sgs_stats_simcomp.f90\
 	simulation_components/scalar_stats_simcomp.f90\
+	simulation_components/scalar_sgs_stats_simcomp.f90\
 	simulation_components/user_stats.f90\
 	simulation_components/field_writer.f90\
 	common/bcknd/cpu/rhs_maker_cpu.f90\
@@ -370,6 +374,8 @@ neko_fortran_SOURCES = \
 	source_terms/bcknd/device/device_gradient_jump_penalty.F90\
 	les/les_model.f90\
 	les/les_model_fctry.f90\
+  les/fluid_sgs_stats.f90\
+  les/scalar_sgs_stats.f90\
 	les/vreman.f90\
 	les/bcknd/cpu/vreman_cpu.f90\
 	les/bcknd/device/vreman_device.f90\
@@ -390,6 +396,10 @@ neko_fortran_SOURCES = \
 	les/bcknd/cpu/wale_cpu.f90\
 	les/bcknd/device/wale_device.f90\
 	les/bcknd/device/device_wale_nut.F90\
+  les/TKE_SGS.f90\
+	les/bcknd/cpu/TKE_SGS_cpu.f90\
+	les/bcknd/device/TKE_SGS_device.f90\
+	les/bcknd/device/device_TKE_SGS_nut.F90\
 	wall_models/wall_model.f90\
 	wall_models/bcknd/device/wall_model_device.F90\
 	wall_models/rough_log_law.f90\
@@ -471,6 +481,7 @@ libneko_la_SOURCES += \
 	les/bcknd/device/hip/smagorinsky_nut.hip\
 	les/bcknd/device/hip/dynamic_smagorinsky_nut.hip\
 	les/bcknd/device/hip/wale_nut.hip\
+	les/bcknd/device/hip/TKE_SGS_nut.hip\
 	wall_models/bcknd/device/hip/wall_model.hip\
 	wall_models/bcknd/device/hip/spalding.hip\
 	wall_models/bcknd/device/hip/rough_log_law.hip\
@@ -531,6 +542,7 @@ libneko_la_SOURCES += \
 	les/bcknd/device/cuda/smagorinsky_nut.cu\
 	les/bcknd/device/cuda/dynamic_smagorinsky_nut.cu\
 	les/bcknd/device/cuda/wale_nut.cu\
+	les/bcknd/device/cuda/TKE_SGS_nut.cu\
 	comm/bcknd/device/cuda/comm_nvshmem.cu\
 	wall_models/bcknd/device/cuda/wall_model.cu\
 	wall_models/bcknd/device/cuda/spalding.cu\
@@ -717,6 +729,7 @@ EXTRA_DIST = \
 	les/bcknd/device/cuda/smagorinsky_nut_kernel.h\
 	les/bcknd/device/cuda/dynamic_smagorinsky_nut_kernel.h\
 	les/bcknd/device/cuda/wale_nut_kernel.h\
+	les/bcknd/device/cuda/TKE_SGS_nut_kernel.h\
 	wall_models/bcknd/device/cuda/wall_model_kernel.h\
 	wall_models/bcknd/device/cuda/spalding_kernel.h\
 	wall_models/bcknd/device/cuda/rough_log_law_kernel.h\
@@ -772,6 +785,7 @@ EXTRA_DIST = \
 	les/bcknd/device/hip/smagorinsky_nut_kernel.h\
 	les/bcknd/device/hip/dynamic_smagorinsky_nut_kernel.h\
 	les/bcknd/device/hip/wale_nut_kernel.h\
+	les/bcknd/device/hip/TKE_SGS_nut_kernel.h\
 	wall_models/bcknd/device/hip/wall_model_kernel.h\
 	wall_models/bcknd/device/hip/spalding_kernel.h\
 	wall_models/bcknd/device/hip/rough_log_law_kernel.h\