Additional fields

Compute additional fields (normal, surface tension, smooth volume fraction, etc.) More...

Functions
subroutine	mod_multiphase_compute_additional_fields::compute_phases_additional_fields ()
	Compute the phases' additional fields.
subroutine	mod_multiphase_compute_additional_fields::compute_phase_additional_fields (phase, representation, boundary_condition)
	Compute the phase's additional fields if needed.
subroutine	mod_multiphase_compute_curvature::compute_curvature (phase, representation)
	Compute the phase curvature field.
subroutine	mod_multiphase_height_function_curvature::compute_curvature_height_function (phase)
	Compute the phase curvature field with a height function method.
subroutine	mod_multiphase_height_function_curvature_2d::compute_curvature_height_function_hybrid_2d (curvature, vf, nn, nh, order)
	Compute the phase curvature field with a hybrid height function method in 2d.
subroutine	mod_multiphase_height_function_curvature_3d::compute_curvature_height_function_hybrid_3d (curvature, vf, nn, nh)
	Compute the phase curvature field with a hybrid height function method in 3d.
subroutine	mod_multiphase_compute_neutral_phase::compute_neutral_phase (phases)
	Compute the neutral phase from other volume fractions.
double precision function	mod_multiphase_compute_phase_volume::compute_phase_volume (volume_fraction)
	Compute the volume of the phase.
subroutine	mod_smooth_volume_fraction::smooth_volume_fraction_inverse_distance (phase, boundary_condition)
	Smooth volume fraction thanks to an inverse distance weighting method.
subroutine	mod_smooth_volume_fraction::smooth_volume_fraction_cubic_b_spline_kernel (phase, boundary_condition)
	Smooth volume fraction thanks to a cubic B-spline kernel convolution.
subroutine	mod_smooth_volume_fraction::compute_smoothing_coef_2d (smoothing_coef, extension, h, i, j)
	Compute smoothing coefficient with cubic_b_spline_kernel (2D)
subroutine	mod_smooth_volume_fraction::compute_smoothing_coef_3d (smoothing_coef, extension, h, i, j, k)
	Compute smoothing coefficient with cubic_b_spline_kernel (3D)
double precision function	mod_smooth_volume_fraction::cubic_b_spline_kernel (x)
	Cubic B-spline kernel.
subroutine	mod_compute_surface_tension_height_function::compute_surface_tension_height_function (phase, sigma, surface_tension_term)
	Compute the surface tension term following the Height function method.

Detailed Description

Compute additional fields (normal, surface tension, smooth volume fraction, etc.)

Function Documentation

compute_curvature()

subroutine mod_multiphase_compute_curvature::compute_curvature	(	type(t_phase_geometry), intent(inout)	phase,
		integer, intent(in)	representation )

Compute the phase curvature field.

Parameters

[in,out]	phase	the phase for which to compute the curvature
[in]	representation	see enum_phase_representation

compute_curvature_height_function()

subroutine mod_multiphase_height_function_curvature::compute_curvature_height_function

(

type(t_phase_geometry), intent(inout), target

phase

)

Compute the phase curvature field with a height function method.

Parameters

[in,out]

phase

the phase for which to compute the curvature

compute_curvature_height_function_hybrid_2d()

subroutine mod_multiphase_height_function_curvature_2d::compute_curvature_height_function_hybrid_2d	(	double precision, dimension(:,:,:), intent(inout)	curvature,
		double precision, dimension(:,:,:), intent(in)	vf,
		integer, intent(in)	nn,
		integer, intent(in)	nh,
		integer, intent(in)	order )

Compute the phase curvature field with a hybrid height function method in 2d.

Parameters

[in,out]	curvature	the curvature computed
[in]	vf	the volume fraction for which to compute the curvature
[in]	nn	the number of stencils
[in]	nh	the height of each stencil
[in]	order	the derivative order

compute_curvature_height_function_hybrid_3d()

subroutine mod_multiphase_height_function_curvature_3d::compute_curvature_height_function_hybrid_3d	(	double precision, dimension(:,:,:), intent(inout)	curvature,
		double precision, dimension(:,:,:), intent(in)	vf,
		integer, intent(in)	nn,
		integer, intent(in)	nh )

Compute the phase curvature field with a hybrid height function method in 3d.

Parameters

[in,out]	curvature	the curvature computed
[in]	vf	the volume fraction for which to compute the curvature
[in]	nn	the number of stencils
[in]	nh	the height of each stencil

compute_neutral_phase()

subroutine mod_multiphase_compute_neutral_phase::compute_neutral_phase

(

type(t_phase_geometry), dimension(:), intent(inout)

phases

)

Compute the neutral phase from other volume fractions.

FIXME remove this

Parameters

phases

all the phases, n=1 being the neutral phase

compute_phase_additional_fields()

subroutine mod_multiphase_compute_additional_fields::compute_phase_additional_fields	(	type(t_phase_geometry), intent(inout)	phase,
		integer, intent(in)	representation,
		type(t_boundary_condition), intent(in)	boundary_condition )

Compute the phase's additional fields if needed.

Parameters

[in,out]	phase	the phase for which to compute the normal
[in]	representation	see enum_phase_representation
[in]	boundary_condition	boundary condition

compute_phase_volume()

double precision function mod_multiphase_compute_phase_volume::compute_phase_volume

(

double precision, dimension(:,:,:), intent(in)

volume_fraction

)

Compute the volume of the phase.

Parameters

volume_fraction

volume fraction

compute_surface_tension_height_function()

subroutine mod_compute_surface_tension_height_function::compute_surface_tension_height_function	(	type(t_phase_geometry), intent(inout)	phase,
		double precision, intent(in)	sigma,
		type(t_face_field), intent(inout)	surface_tension_term )

Compute the surface tension term following the Height function method.

Parameters

[in,out]	phase	the phase for which to compute the normal
[in]	sigma	the surface tension coefficient
[in]	surface_tension_term	the surface tension term added into RHS

cubic_b_spline_kernel()

double precision function mod_smooth_volume_fraction::cubic_b_spline_kernel

(

double precision, intent(in)

x

)

Cubic B-spline kernel.

remark: coefficient \(40/(7 \pi h^2)\) is not taken into account because of the normalization step

smooth_volume_fraction_cubic_b_spline_kernel()

subroutine mod_smooth_volume_fraction::smooth_volume_fraction_cubic_b_spline_kernel	(	type(t_phase_geometry), intent(inout)	phase,
		type(t_boundary_condition), intent(in)	boundary_condition )

Smooth volume fraction thanks to a cubic B-spline kernel convolution.

remark: to be improved for irregular grids