Variables
character, parameter	dflt_char = ' '
	Default string value.

real(wp), parameter	dflt_real = -1e6_wp
	Default real value.

real(wp), parameter	sgm_eps = 1e-16_wp
	Segmentation tolerance.

real(wp), parameter	small_alf = 1e-11_wp
	Small alf tolerance.

real(wp), parameter	pi = 3.141592653589793_wp
	Pi.

real(wp), parameter	verysmall = 1.e-12_wp
	Very small number.

integer, parameter	num_stcls_min = 5
	Minimum # of stencils.

integer, parameter	path_len = 400
	Maximum path length.

integer, parameter	name_len = 50
	Maximum name length.

integer, parameter	dflt_int = -100
	Default integer value.

integer, parameter	fourier_rings = 5
	Fourier filter ring limit.

integer, parameter	num_fluids_max = 10
	Maximum number of fluids in the simulation.

integer, parameter	num_probes_max = 10
	Maximum number of flow probes in the simulation.

integer, parameter	num_patches_max = 10

integer, parameter	pathlen_max = 400

integer, parameter	nnode = 4
	Number of QBMM nodes.

integer, parameter	gp_layers = 3
	Number of ghost point layers for IBM.

real(wp), parameter	capillary_cutoff = 1e-6
	color function gradient magnitude at which to apply the surface tension fluxes

real(wp), parameter	acoustic_spatial_support_width = 2.5_wp
	Spatial support width of acoustic source, used in s_source_spatial.

real(wp), parameter	dflt_vcfl_dt = 100._wp
	value of vcfl_dt when viscosity is off for computing adaptive timestep size

real(wp), parameter	broadband_spectral_level_constant = 20._wp
	The constant to scale the spectral level at the lower frequency bound.

real(wp), parameter	broadband_spectral_level_growth_rate = 10._wp
	The spectral level constant to correct the magnitude at each frqeuency to ensure the source is overall broadband.

real(wp), parameter	dflt_t_guess = 1200._wp

integer, parameter	ifactor_2d = 50
	Multiple factor of the ratio (edge to cell width) for interpolation along edges for 2D models.

integer, parameter	ifactor_3d = 5
	Multiple factor of the ratio (edge to cell width) for interpolation along edges for 3D models.

integer, parameter	ifactor_bary_3d = 20
	Multiple factor of the ratio (triangle area to cell face area) for interpolation on triangle facets for 3D models.

integer, parameter	num_ray = 20
	Default number of rays traced per cell.

real(wp), parameter	ray_tracing_threshold = 0.9_wp
	Threshold above which the cell is marked as the model patch.

real(wp), parameter	threshold_vector_zero = 1e-10
	Threshold to treat the component of a vector to be zero.

real(wp), parameter	threshold_edge_zero = 1e-10
	Threshold to treat two edges to be overlapped.

real(wp), parameter	threshold_bary = 1e-1
	Threshold to interpolate a barycentric facet.

real(wp), parameter	initial_distance_buffer = 1e12_wp
	Initialized levelset distance for the shortest path pair algorithm.

integer, parameter	mapcells = 3
	Number of cells around the bubble where the smoothening function will have effect.

real(wp), parameter	r_uni = 8314._wp

integer, parameter	num_ts_rkck = 6
	Number of time-stages in the RKCK stepper.

real(wp), parameter	rkck_c1 = 0._wp

real(wp), parameter	rkck_c2 = 0.2_wp

real(wp), parameter	rkck_c3 = 0.3_wp

real(wp), parameter	rkck_c4 = 0.6_wp

real(wp), parameter	rkck_c5 = 1._wp

real(wp), parameter	rkck_c6 = 0.875_wp

real(wp), dimension(6), parameter	rkck_coef1 = (/0.2_wp, 0._wp, 0._wp, 0._wp, 0._wp, 0._wp/)

real(wp), dimension(6), parameter	rkck_coef2 = (/3._wp/40._wp, 9._wp/40._wp, 0._wp, 0._wp, 0._wp, 0._wp/)

real(wp), dimension(6), parameter	rkck_coef3 = (/0.3_wp, -0.9_wp, 1.2_wp, 0._wp, 0._wp, 0._wp/)

real(wp), dimension(6), parameter	rkck_coef4 = (/-11._wp/54._wp, 2.5_wp, -70._wp/27._wp, 35._wp/27._wp, 0._wp, 0._wp/)

real(wp), dimension(6), parameter	rkck_coef5 = (/1631._wp/55296._wp, 175._wp/512._wp, 575._wp/13824._wp, 44275._wp/110592._wp, 253._wp/4096._wp, 0._wp/)

real(wp), dimension(6), parameter	rkck_coef6 = (/37._wp/378._wp, 0._wp, 250._wp/621._wp, 125._wp/594._wp, 0._wp, 512._wp/1771._wp/)

real(wp), dimension(6), parameter	rkck_coefe = (/37._wp/378._wp - 2825._wp/27648._wp, 0._wp, 250._wp/621._wp - 18575._wp/48384._wp, 125._wp/594._wp - 13525._wp/55296._wp, -277._wp/14336._wp, 512._wp/1771._wp - 0.25_wp/)

real(wp), parameter	verysmall_dt = 1e-14_wp
	Very small dt, stop stepper.

real(wp), parameter	safety = 0.9_wp
	Next dt will be maximum 0.9*dt if truncation error is above tolerance.

real(wp), parameter	rnddec = 1e8_wp
	Round calculated dt (16th digit) to avoid the inclusion of random decimals.

real(wp), parameter	pshrnk = -0.25_wp
	Factor to reduce dt when truncation error above tolerance.

real(wp), parameter	shrnkdt = 0.5_wp
	Factor to reduce dt due to negative bubble radius.

real(wp), parameter	errcon = 1.89e-4_wp
	Limit to slightly increase dt when truncation error is between ERRCON and 1.

real(wp), parameter	pgrow = -0.2_wp
	Factor to increase dt when truncation error is between ERRCON and 1.

integer, parameter	case_file_error_code = 22

Variable Documentation

◆ acoustic_spatial_support_width

real(wp), parameter m_constants::acoustic_spatial_support_width = 2.5_wp

Spatial support width of acoustic source, used in s_source_spatial.

◆ broadband_spectral_level_constant

real(wp), parameter m_constants::broadband_spectral_level_constant = 20._wp

The constant to scale the spectral level at the lower frequency bound.

◆ broadband_spectral_level_growth_rate

real(wp), parameter m_constants::broadband_spectral_level_growth_rate = 10._wp

The spectral level constant to correct the magnitude at each frqeuency to ensure the source is overall broadband.

◆ capillary_cutoff

real(wp), parameter m_constants::capillary_cutoff = 1e-6

color function gradient magnitude at which to apply the surface tension fluxes

◆ case_file_error_code

integer, parameter m_constants::case_file_error_code = 22

◆ dflt_char

character, parameter m_constants::dflt_char = ' '

Default string value.

◆ dflt_int

integer, parameter m_constants::dflt_int = -100

Default integer value.

◆ dflt_real

real(wp), parameter m_constants::dflt_real = -1e6_wp

Default real value.

◆ dflt_t_guess

real(wp), parameter m_constants::dflt_t_guess = 1200._wp

◆ dflt_vcfl_dt

real(wp), parameter m_constants::dflt_vcfl_dt = 100._wp

value of vcfl_dt when viscosity is off for computing adaptive timestep size

◆ errcon

real(wp), parameter m_constants::errcon = 1.89e-4_wp

Limit to slightly increase dt when truncation error is between ERRCON and 1.

◆ fourier_rings

integer, parameter m_constants::fourier_rings = 5

Fourier filter ring limit.

◆ gp_layers

integer, parameter m_constants::gp_layers = 3

Number of ghost point layers for IBM.

◆ ifactor_2d

integer, parameter m_constants::ifactor_2d = 50

Multiple factor of the ratio (edge to cell width) for interpolation along edges for 2D models.

◆ ifactor_3d

integer, parameter m_constants::ifactor_3d = 5

Multiple factor of the ratio (edge to cell width) for interpolation along edges for 3D models.

◆ ifactor_bary_3d

integer, parameter m_constants::ifactor_bary_3d = 20

Multiple factor of the ratio (triangle area to cell face area) for interpolation on triangle facets for 3D models.

◆ initial_distance_buffer

real(wp), parameter m_constants::initial_distance_buffer = 1e12_wp

Initialized levelset distance for the shortest path pair algorithm.

◆ mapcells

integer, parameter m_constants::mapcells = 3

Number of cells around the bubble where the smoothening function will have effect.

◆ name_len

integer, parameter m_constants::name_len = 50

Maximum name length.

◆ nnode

integer, parameter m_constants::nnode = 4

Number of QBMM nodes.

◆ num_fluids_max

integer, parameter m_constants::num_fluids_max = 10

Maximum number of fluids in the simulation.

◆ num_patches_max

integer, parameter m_constants::num_patches_max = 10

◆ num_probes_max

integer, parameter m_constants::num_probes_max = 10

Maximum number of flow probes in the simulation.

◆ num_ray

integer, parameter m_constants::num_ray = 20

Default number of rays traced per cell.

◆ num_stcls_min

integer, parameter m_constants::num_stcls_min = 5

Minimum # of stencils.

◆ num_ts_rkck

integer, parameter m_constants::num_ts_rkck = 6

Number of time-stages in the RKCK stepper.

◆ path_len

integer, parameter m_constants::path_len = 400

Maximum path length.

◆ pathlen_max

integer, parameter m_constants::pathlen_max = 400

◆ pgrow

real(wp), parameter m_constants::pgrow = -0.2_wp

Factor to increase dt when truncation error is between ERRCON and 1.

◆ pi

real(wp), parameter m_constants::pi = 3.141592653589793_wp

Pi.

◆ pshrnk

real(wp), parameter m_constants::pshrnk = -0.25_wp

Factor to reduce dt when truncation error above tolerance.

◆ r_uni

real(wp), parameter m_constants::r_uni = 8314._wp

◆ ray_tracing_threshold

real(wp), parameter m_constants::ray_tracing_threshold = 0.9_wp

Threshold above which the cell is marked as the model patch.

◆ rkck_c1

real(wp), parameter m_constants::rkck_c1 = 0._wp

◆ rkck_c2

real(wp), parameter m_constants::rkck_c2 = 0.2_wp

◆ rkck_c3

real(wp), parameter m_constants::rkck_c3 = 0.3_wp

◆ rkck_c4

real(wp), parameter m_constants::rkck_c4 = 0.6_wp

◆ rkck_c5

real(wp), parameter m_constants::rkck_c5 = 1._wp

◆ rkck_c6

real(wp), parameter m_constants::rkck_c6 = 0.875_wp

◆ rkck_coef1

real(wp), dimension(6), parameter m_constants::rkck_coef1 = (/0.2_wp, 0._wp, 0._wp, 0._wp, 0._wp, 0._wp/)

◆ rkck_coef2

real(wp), dimension(6), parameter m_constants::rkck_coef2 = (/3._wp/40._wp, 9._wp/40._wp, 0._wp, 0._wp, 0._wp, 0._wp/)

◆ rkck_coef3

real(wp), dimension(6), parameter m_constants::rkck_coef3 = (/0.3_wp, -0.9_wp, 1.2_wp, 0._wp, 0._wp, 0._wp/)

◆ rkck_coef4

real(wp), dimension(6), parameter m_constants::rkck_coef4 = (/-11._wp/54._wp, 2.5_wp, -70._wp/27._wp, 35._wp/27._wp, 0._wp, 0._wp/)

◆ rkck_coef5

real(wp), dimension(6), parameter m_constants::rkck_coef5 = (/1631._wp/55296._wp, 175._wp/512._wp, 575._wp/13824._wp, 44275._wp/110592._wp, 253._wp/4096._wp, 0._wp/)

◆ rkck_coef6

real(wp), dimension(6), parameter m_constants::rkck_coef6 = (/37._wp/378._wp, 0._wp, 250._wp/621._wp, 125._wp/594._wp, 0._wp, 512._wp/1771._wp/)

◆ rkck_coefe

real(wp), dimension(6), parameter m_constants::rkck_coefe = (/37._wp/378._wp - 2825._wp/27648._wp, 0._wp, 250._wp/621._wp - 18575._wp/48384._wp, 125._wp/594._wp - 13525._wp/55296._wp, -277._wp/14336._wp, 512._wp/1771._wp - 0.25_wp/)

◆ rnddec

real(wp), parameter m_constants::rnddec = 1e8_wp

Round calculated dt (16th digit) to avoid the inclusion of random decimals.

◆ safety

real(wp), parameter m_constants::safety = 0.9_wp

Next dt will be maximum 0.9*dt if truncation error is above tolerance.

◆ sgm_eps

real(wp), parameter m_constants::sgm_eps = 1e-16_wp

Segmentation tolerance.

◆ shrnkdt

real(wp), parameter m_constants::shrnkdt = 0.5_wp

Factor to reduce dt due to negative bubble radius.

◆ small_alf

real(wp), parameter m_constants::small_alf = 1e-11_wp

Small alf tolerance.

◆ threshold_bary

real(wp), parameter m_constants::threshold_bary = 1e-1

Threshold to interpolate a barycentric facet.

◆ threshold_edge_zero

real(wp), parameter m_constants::threshold_edge_zero = 1e-10

Threshold to treat two edges to be overlapped.

◆ threshold_vector_zero

real(wp), parameter m_constants::threshold_vector_zero = 1e-10

Threshold to treat the component of a vector to be zero.

◆ verysmall

real(wp), parameter m_constants::verysmall = 1.e-12_wp

Very small number.

◆ verysmall_dt

real(wp), parameter m_constants::verysmall_dt = 1e-14_wp

Very small dt, stop stepper.

Variables