variables_diagnostic_quantities Module Reference

Diagnostic quantity switches. More...

Functions/Subroutines
subroutine	initialize_variables_diagnostic_quantities ()
	Allocate memory for diagnostic quantites.

Variables
logical	is_final_diagnostic_enabled = .true.
	Enable the computation of diagnostic quantities only after the time loop (default: true)
IB surface.
logical	diagnostic_compute_ib_surface = .false.
logical	diagnostic_print_ib_surface = .false.
double precision, dimension(:), allocatable	diagnostic_ib_surface
logical	diagnostic_compute_ib_drag_lift = .false.
logical	diagnostic_print_ib_drag_lift = .false.
double precision, dimension(:,:), allocatable	diagnostic_ib_pressure_drag_lift
double precision, dimension(:,:), allocatable	diagnostic_ib_shear_drag_lift
logical	diagnostic_compute_ib_nusselt_number = .false.
logical	diagnostic_print_ib_nusselt_number = .false.
double precision, dimension(:), allocatable	diagnostic_ib_nusselt_number
Volume of phases.
logical	diagnostic_compute_phases_volume = .false.
logical	diagnostic_print_phases_volume = .false.
double precision, dimension(:), allocatable	diagnostic_phases_volume
Kinetic energy
logical	diagnostic_compute_kinetic_energy = .false.
logical	diagnostic_print_kinetic_energy = .false.
Turbulent kinetic energy
logical	diagnostic_compute_turbulent_kinetic_energy = .false.
logical	diagnostic_print_turbulent_kinetic_energy = .false.
Kinetic energy dissipation rate
logical	diagnostic_compute_kinetic_energy_dissipation_rate = .false.
logical	diagnostic_print_kinetic_energy_dissipation_rate = .false.
Mean kinetic energy.
logical	diagnostic_compute_mean_kinetic_energy = .false.
logical	diagnostic_print_mean_kinetic_energy = .false.
double precision	diagnostic_mean_kinetic_energy = 0d0
Kinetic energy integral
logical	diagnostic_compute_kinetic_energy_integral = .false.
logical	diagnostic_print_kinetic_energy_integral = .false.
double precision	diagnostic_kinetic_energy_integral = 0d0
Mean longitudinal pressure difference.
logical	diagnostic_compute_mean_longitudinal_pressure_difference = .false.
logical	diagnostic_print_mean_longitudinal_pressure_difference = .false.
double precision	diagnostic_mean_longitudinal_pressure_difference = 0d0
Wall shear stress.
logical	diagnostic_compute_wall_shear_stress = .false.
logical	diagnostic_print_wall_shear_stress = .false.
double precision	diagnostic_wall_shear_stress_bottom = 0d0
double precision	diagnostic_wall_shear_stress_top = 0d0
Mean time wall shear stress.
logical	diagnostic_compute_wall_shear_stress_time_averaged = .false.
logical	diagnostic_print_wall_shear_stress_time_averaged = .false.
double precision	diagnostic_wall_shear_stress_time_averaged_bottom = 0d0
double precision	diagnostic_wall_shear_stress_time_averaged_top = 0d0
Wall Nusselt numbers.
logical	diagnostic_compute_nusselt_number = .false.
logical	diagnostic_print_nusselt_number = .false.
logical	diagnostic_print_nusselt_number_time_averaged = .false.
double precision	diagnostic_left_nusselt = 0d0
double precision	diagnostic_right_nusselt = 0d0
double precision	diagnostic_bottom_nusselt = 0d0
double precision	diagnostic_top_nusselt = 0d0
double precision	diagnostic_back_nusselt = 0d0
double precision	diagnostic_front_nusselt = 0d0
double precision	diagnostic_left_nusselt_n = 0d0
double precision	diagnostic_right_nusselt_n = 0d0
double precision	diagnostic_bottom_nusselt_n = 0d0
double precision	diagnostic_top_nusselt_n = 0d0
double precision	diagnostic_back_nusselt_n = 0d0
double precision	diagnostic_front_nusselt_n = 0d0
double precision	diagnostic_ib_wall_nusselt = 0d0
double precision	diagnostic_ib_wall_nusselt_left = 0d0
double precision	diagnostic_ib_wall_nusselt_right = 0d0
double precision	diagnostic_ib_wall_nusselt_bottom = 0d0
double precision	diagnostic_ib_wall_nusselt_top = 0d0
double precision	diagnostic_ib_wall_nusselt_back = 0d0
double precision	diagnostic_ib_wall_nusselt_front = 0d0
double precision	diagnostic_nusselt_delta_temperature = 1d0
Wall Sherwood numbers.
logical, dimension(:), allocatable	diagnostic_compute_sherwood_number
logical, dimension(:), allocatable	diagnostic_print_sherwood_number
double precision, dimension(:), allocatable	diagnostic_left_sherwood
double precision, dimension(:), allocatable	diagnostic_right_sherwood
double precision, dimension(:), allocatable	diagnostic_bottom_sherwood
double precision, dimension(:), allocatable	diagnostic_top_sherwood
double precision, dimension(:), allocatable	diagnostic_back_sherwood
double precision, dimension(:), allocatable	diagnostic_front_sherwood
double precision	diagnostic_sherwood_delta_concentration = 1d0
Wall shear velocity.
logical	diagnostic_compute_wall_shear_velocity = .false.
logical	diagnostic_print_wall_shear_velocity = .false.
double precision	diagnostic_wall_shear_velocity_bottom = 0d0
double precision	diagnostic_wall_shear_velocity_top = 0d0
Mean time wall shear velocity.
logical	diagnostic_compute_wall_shear_velocity_time_averaged = .false.
logical	diagnostic_print_wall_shear_velocity_time_averaged = .false.
double precision	diagnostic_wall_shear_velocity_time_averaged_bottom = 0d0
double precision	diagnostic_wall_shear_velocity_time_averaged_top = 0d0
Wall skin friction coefficients.
logical	diagnostic_compute_wall_skin_friction_coef = .false.
logical	diagnostic_print_wall_skin_friction_coef = .false.
double precision	diagnostic_wall_skin_friction_coef_bottom = 0d0
double precision	diagnostic_wall_skin_friction_coef_top = 0d0
Mean time wall skin friction coefficients.
logical	diagnostic_compute_wall_skin_friction_coef_time_averaged = .false.
logical	diagnostic_print_wall_skin_friction_coef_time_averaged = .false.
double precision	diagnostic_wall_skin_friction_coef_time_averaged_bottom = 0d0
double precision	diagnostic_wall_skin_friction_coef_time_averaged_top = 0d0
logical	diagnostic_compute_y_plus = .false.
logical	diagnostic_compute_u_plus = .false.
logical	diagnostic_compute_k_plus = .false.
logical	diagnostic_compute_local_wall_skin_friction_coef = .false.
logical	diagnostic_compute_local_wall_shear_stress = .false.
logical	diagnostic_compute_local_nusselt = .false.
Domain volume for the energy equation.
logical	diagnostic_compute_domain_volume_energy = .false.
logical	diagnostic_print_domain_volume_energy = .false.
double precision	diagnostic_domain_volume_energy = 0d0
Front position of Stefan problem (melting)
logical	diagnostic_compute_stefan_problem_front_position_melting = .false.
logical	diagnostic_print_stefan_problem_front_position_melting = .false.
double precision	diagnostic_stefan_problem_front_position_melting = 0d0
Front position of Stefan problem (solidification)
logical	diagnostic_compute_stefan_problem_front_position_solidification = .false.
logical	diagnostic_print_stefan_problem_front_position_solidification = .false.
double precision	diagnostic_stefan_problem_front_position_solidification = 0d0
Domain volume for the Navier-Stokes equation.
logical	diagnostic_compute_domain_volume_navier = .false.
logical	diagnostic_print_domain_volume_navier = .false.
double precision	diagnostic_domain_volume_navier_cell = 0d0
double precision	diagnostic_domain_volume_navier_face_u = 0d0
double precision	diagnostic_domain_volume_navier_face_v = 0d0
double precision	diagnostic_domain_volume_navier_face_w = 0d0
Domain volume for the species transport equation.
logical	diagnostic_compute_domain_volume_species_transport = .false.
logical	diagnostic_print_domain_volume_species_transport = .false.
double precision	diagnostic_domain_volume_species_transport = 0d0
Domain volume for the phase advection equation.
logical	diagnostic_compute_domain_volume_phase_advection = .false.
logical	diagnostic_print_domain_volume_phase_advection = .false.
double precision	diagnostic_domain_volume_phase_advection = 0d0
Min species concentration.
logical, dimension(:), allocatable	diagnostic_compute_min_species_concentration
logical, dimension(:), allocatable	diagnostic_print_min_species_concentration
double precision, dimension(:), allocatable	diagnostic_min_species_concentration
Max species concentration.
logical, dimension(:), allocatable	diagnostic_compute_max_species_concentration
logical, dimension(:), allocatable	diagnostic_print_max_species_concentration
double precision, dimension(:), allocatable	diagnostic_max_species_concentration
Mean species concentration.
logical, dimension(:), allocatable	diagnostic_compute_mean_species_concentration
logical, dimension(:), allocatable	diagnostic_print_mean_species_concentration
double precision, dimension(:), allocatable	diagnostic_mean_species_concentration
Min volume fraction.
logical, dimension(:), allocatable	diagnostic_compute_min_volume_fraction
logical, dimension(:), allocatable	diagnostic_print_min_volume_fraction
double precision, dimension(:), allocatable	diagnostic_min_volume_fraction
Max volume fraction.
logical, dimension(:), allocatable	diagnostic_compute_max_volume_fraction
logical, dimension(:), allocatable	diagnostic_print_max_volume_fraction
double precision, dimension(:), allocatable	diagnostic_max_volume_fraction
Mean volume fraction.
logical, dimension(:), allocatable	diagnostic_compute_mean_volume_fraction
logical, dimension(:), allocatable	diagnostic_print_mean_volume_fraction
double precision, dimension(:), allocatable	diagnostic_mean_volume_fraction
Min quantites.
logical, dimension(diagnostic_field_size)	diagnostic_compute_min = .false.
logical, dimension(diagnostic_field_size)	diagnostic_print_min = .false.
double precision, dimension(diagnostic_field_size)	diagnostic_min = 0d0
Max quantites.
logical, dimension(diagnostic_field_size)	diagnostic_compute_max = .false.
logical, dimension(diagnostic_field_size)	diagnostic_print_max = .false.
double precision, dimension(diagnostic_field_size)	diagnostic_max = 0d0
Mean quantites.
logical, dimension(diagnostic_field_size)	diagnostic_compute_mean = .false.
logical, dimension(diagnostic_field_size)	diagnostic_print_mean = .false.
double precision, dimension(diagnostic_field_size)	diagnostic_mean = 0d0
Volume loss
logical	diagnostic_compute_materials_volume_loss = .false.
logical	diagnostic_print_materials_volume_loss = .false.
double precision, dimension(:), allocatable	diagnostic_materials_initial_volume
double precision, dimension(:), allocatable	diagnostic_materials_volume_loss
Mass conservation
logical	diagnostic_compute_mass_loss = .false.
logical	diagnostic_print_mass_loss = .false.
double precision	diagnostic_initial_mass
double precision	diagnostic_mass_loss
Energy conservation
logical	diagnostic_compute_energy_budget = .false.
logical	diagnostic_print_energy_budget = .false.
double precision	diagnostic_internal_energy_domain = 0d0
double precision	diagnostic_energy_boundaries = 0d0
double precision	diagnostic_energy_loss = 0d0
Momentum quantities.
logical	diagnostic_compute_momentum_integral = .false.
logical	diagnostic_print_momentum_integral = .false.
double precision, dimension(:), allocatable	diagnostic_momentum_integral
Reference solution of the density
logical	diagnostic_compute_reference_solution_error_density = .false.
logical	diagnostic_print_reference_solution_error_density = .false.
double precision	diagnostic_reference_solution_error_density_l1 = 0d0
double precision	diagnostic_reference_solution_error_density_l2 = 0d0
double precision	diagnostic_reference_solution_error_density_linf = 0d0
Reference solution of the temperature
logical	diagnostic_compute_reference_solution_error_temperature = .false.
logical	diagnostic_print_reference_solution_error_temperature = .false.
double precision	diagnostic_reference_solution_error_temperature_l1 = 0d0
double precision	diagnostic_reference_solution_error_temperature_l2 = 0d0
double precision	diagnostic_reference_solution_error_temperature_linf = 0d0
Reference solution of the pressure
logical	diagnostic_compute_reference_solution_error_pressure = .false.
logical	diagnostic_print_reference_solution_error_pressure = .false.
double precision	diagnostic_reference_solution_error_pressure_l1 = 0d0
double precision	diagnostic_reference_solution_error_pressure_l2 = 0d0
double precision	diagnostic_reference_solution_error_pressure_linf = 0d0
Reference solution of the velocity
logical	diagnostic_compute_reference_solution_error_velocity = .false.
logical	diagnostic_print_reference_solution_error_velocity = .false.
double precision	diagnostic_reference_solution_error_velocity_l1 = 0d0
double precision	diagnostic_reference_solution_error_velocity_l2 = 0d0
double precision	diagnostic_reference_solution_error_velocity_linf = 0d0
Reference solution of the concentration
logical, dimension(:), allocatable	diagnostic_compute_reference_solution_error_concentration
logical, dimension(:), allocatable	diagnostic_print_reference_solution_error_concentration
double precision, dimension(:), allocatable	diagnostic_reference_solution_error_concentration_l1
double precision, dimension(:), allocatable	diagnostic_reference_solution_error_concentration_l2
double precision, dimension(:), allocatable	diagnostic_reference_solution_error_concentration_linf
Reference solution of the volume fraction
logical, dimension(:), allocatable	diagnostic_compute_reference_solution_error_volume_fraction
logical, dimension(:), allocatable	diagnostic_print_reference_solution_error_volume_fraction
double precision, dimension(:), allocatable	diagnostic_reference_solution_error_volume_fraction_l1
double precision, dimension(:), allocatable	diagnostic_reference_solution_error_volume_fraction_l2
double precision, dimension(:), allocatable	diagnostic_reference_solution_error_volume_fraction_linf
Reference solution of area of symmetric difference.
logical	diagnostic_compute_reference_solution_symmetric_difference_area = .false.
logical	diagnostic_print_reference_solution_symmetric_difference_area = .false.
double precision	diagnostic_reference_solution_symmetric_difference_area = 0d0
Q criterion
logical	diagnostic_compute_q_criterion = .false.
logical	diagnostic_print_q_criterion = .false.
Lagrangian acceleration
logical	diagnostic_compute_lagrangian_acceleration = .false.
logical	diagnostic_print_lagrangian_acceleration = .false.
Strain rate magnitude
logical	diagnostic_compute_strain_rate_magnitude = .false.
logical	diagnostic_print_strain_rate_magnitude = .false.
Vorticity
logical	diagnostic_compute_vorticity = .false.
logical	diagnostic_print_vorticity = .false.
Turbulent kinetic energy dissipation rate
logical	diagnostic_compute_tke_dissipation_rate = .false.
logical	diagnostic_print_tke_dissipation_rate = .false.
Kolmogorov scale
logical	diagnostic_compute_kolmogorov_scale = .false.
logical	diagnostic_print_kolmogorov_scale = .false.
Batchelor scale
logical	diagnostic_compute_batchelor_scale = .false.
logical	diagnostic_print_batchelor_scale = .false.
Turbulent Kolmogorov timescale
logical	diagnostic_compute_kolmogorov_timescale = .false.
logical	diagnostic_print_kolmogorov_timescale = .false.

Detailed Description

Diagnostic quantities are physical tools that returns a scalar value (Nusselt, mean quantities, etc.).

Function/Subroutine Documentation

initialize_variables_diagnostic_quantities()

subroutine variables_diagnostic_quantities::initialize_variables_diagnostic_quantities

Must be called when the number of immersed boundaries, the number of phases and the number of species is known.