Namespace for "global" problem parameters. More...

Enumerations
enum	Cases { Spherical_cap_in_cylinder_pinned , All_pinned , Barrel_shape , T_junction_with_nonzero_contact_angle }
	Enumeration for the possible cases. More...

Functions
double	get_exact_kappa ()
	Exact kappa. More...

void	spine_base_function (const Vector< double > &x, Vector< double > &spine_B, Vector< Vector< double > > &dspine_B)
	Spine basis: The position vector to the basis of the spine as a function of the two coordinates x_1 and x_2, and its derivatives w.r.t. to these coordinates. dspine_B[i][j] = d spine_B[j] / dx_i Spines start in the (x_1,x_2) plane at (x_1,x_2). More...

void	spine_function (const Vector< double > &x, Vector< double > &spine, Vector< Vector< double > > &dspine)
	Spine: The spine vector field as a function of the two coordinates x_1 and x_2, and its derivatives w.r.t. to these coordinates: dspine[i][j] = d spine[j] / dx_i. More...

void	setup_dependent_parameters_and_sanity_check ()
	Setup dependent parameters and perform sanity check. More...

std::complex< double >	I (0.0, 1.0)
	Imaginary unit. More...

void	get_exact_u (const Vector< double > &x, Vector< double > &u)
	Exact solution for scattered field (vector returns real and impaginary parts). More...

void	prescribed_incoming_flux (const Vector< double > &x, complex< double > &flux)
	Flux (normal derivative) on the unit disk for a planar incoming wave. More...

Variables
double	Controlled_height = 0.0
	Height control value. More...

double	Alpha_min = MathematicalConstants::Pi/2.0*1.5
	Min. spine angle against horizontal plane. More...

double	Alpha_max = MathematicalConstants::Pi/2.0*0.5
	Max. spine angle against horizontal plane. More...

bool	Use_spines = true
	Use spines (true) or not (false) More...

bool	Use_height_control = true
	Use height control (true) or not (false)? More...

int	Case = All_pinned
	What case are we considering: Choose one from the enumeration Cases. More...

double	Gamma = MathematicalConstants::Pi/4.0
	Contact angle and its cos (dependent parameter – is reassigned) More...

double	Cos_gamma =cos(Gamma)
	Cos of contact angle. More...

Data *	Kappa_pt = 0
	Pointer to Data object that stores the prescribed curvature. More...

double	Kappa_initial = 0.0
	Initial value for kappa. More...

int	Step_sign = 1
	Increase or decrease the value of the control parameters? More...

unsigned	Nsteps = 5
	Number of steps. More...

double	Kappa_increment = -0.05
	Increment for prescribed curvature. More...

double	Controlled_height_increment = 0.1
	Increment for height control. More...

unsigned	Control_element = 0
	Number of element in bulk mesh at which height control is applied. Initialise to 0 – will be overwritte in setup_dependent_parameters_and_sanity_check() More...

double	L_x = 1.0
	Length and width of the domain. More...

double	L_y = 1.0
	Width of domain. More...

unsigned	N_x = 8
	Number of elements in the mesh. More...

unsigned	N_y = 8

double	Beta_min = MathematicalConstants::Pi/2.0
	Min. second spine angle against horizontal plane. More...

double	Beta_max = MathematicalConstants::Pi/2.0
	Max. second pine angle against horizontal plane. More...

bool	Rotate_spines_in_both_directions = true
	Should the spines rotate in the x and y directions (true)? More...

double	K_squared =10.0
	Square of the wavenumber. More...

unsigned	N_fourier =10
	Number of terms used in the computation of the exact solution. More...

bool	DtN_BC =false
	Flag to choose the Dirichlet to Neumann BC or ABC BC. More...

unsigned	ABC_order =3
	Flag to choose wich order to use. More...

double	Outer_radius =1.5
	Radius of outer boundary (must be a circle!) More...

Detailed Description

Namespace for "global" problem parameters.

////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////

Namespace for the Helmholtz problem parameters

Enumeration Type Documentation

◆ Cases

enum GlobalParameters::Cases

Enumeration for the possible cases.

Enumerator
Spherical_cap_in_cylinder_pinned
All_pinned
Barrel_shape
T_junction_with_nonzero_contact_angle

Definition at line 46 of file common_young_laplace_stuff.h.

Function Documentation

◆ get_exact_kappa()

double GlobalParameters::get_exact_kappa ( )

Exact kappa.

Definition at line 54 of file barrel.cc.

References Controlled_height.

Referenced by YoungLaplaceProblem< ELEMENT >::doc_solution().

◆ get_exact_u()

void GlobalParameters::get_exact_u	(	const Vector< double > &	x,
		Vector< double > &	u
	)

Exact solution for scattered field (vector returns real and impaginary parts).

Definition at line 86 of file scattering.cc.

References I(), K_squared, and N_fourier.

Referenced by ScatteringProblem< ELEMENT >::doc_solution().

◆ I()

std::complex< double > GlobalParameters::I	(	0.	0,
		1.	0
	)

Imaginary unit.

Referenced by get_exact_u(), and prescribed_incoming_flux().

◆ prescribed_incoming_flux()

void GlobalParameters::prescribed_incoming_flux	(	const Vector< double > &	x,
		complex< double > &	flux
	)

Flux (normal derivative) on the unit disk for a planar incoming wave.

Definition at line 156 of file scattering.cc.

References I(), K_squared, and N_fourier.

Referenced by ScatteringProblem< ELEMENT >::set_prescribed_incoming_flux_pt().

◆ setup_dependent_parameters_and_sanity_check()

void GlobalParameters::setup_dependent_parameters_and_sanity_check ( )

Setup dependent parameters and perform sanity check.

Definition at line 126 of file common_young_laplace_stuff.h.

References All_pinned, Alpha_max, Alpha_min, Barrel_shape, Case, Control_element, Cos_gamma, Gamma, Kappa_initial, L_x, L_y, N_x, N_y, Rotate_spines_in_both_directions, Spherical_cap_in_cylinder_pinned, and T_junction_with_nonzero_contact_angle.

◆ spine_base_function()

void GlobalParameters::spine_base_function	(	const Vector< double > &	x,
		Vector< double > &	spine_B,
		Vector< Vector< double > > &	dspine_B
	)

Spine basis: The position vector to the basis of the spine as a function of the two coordinates x_1 and x_2, and its derivatives w.r.t. to these coordinates. dspine_B[i][j] = d spine_B[j] / dx_i Spines start in the (x_1,x_2) plane at (x_1,x_2).

Definition at line 72 of file barrel.cc.

Referenced by RefineableYoungLaplaceProblem< ELEMENT >::RefineableYoungLaplaceProblem(), and YoungLaplaceProblem< ELEMENT >::YoungLaplaceProblem().

◆ spine_function()

void GlobalParameters::spine_function	(	const Vector< double > &	x,
		Vector< double > &	spine,
		Vector< Vector< double > > &	dspine
	)

Spine: The spine vector field as a function of the two coordinates x_1 and x_2, and its derivatives w.r.t. to these coordinates: dspine[i][j] = d spine[j] / dx_i.

Spines (and derivatives) are independent of x[0] and rotate in the x[1]-direction

Spines are dependent of x[0] AND x[1] and rotate in both directions

Spines (and derivatives) are independent of x[0] and rotate in the x[1]-direction

Spines are dependent of x[0] AND x[1] and rotate in both directions

Spines (and derivatives) are independent of x[0] and rotate in the x[1]-direction

Definition at line 104 of file barrel.cc.

References Alpha_max, and Alpha_min.

Referenced by RefineableYoungLaplaceProblem< ELEMENT >::RefineableYoungLaplaceProblem(), and YoungLaplaceProblem< ELEMENT >::YoungLaplaceProblem().

Variable Documentation

◆ ABC_order

unsigned GlobalParameters::ABC_order =3

Flag to choose wich order to use.

Definition at line 76 of file scattering.cc.

Referenced by main(), and ScatteringProblem< ELEMENT >::setup_outer_boundary().

◆ Alpha_max

double GlobalParameters::Alpha_max = MathematicalConstants::Pi/2.0*0.5

Max. spine angle against horizontal plane.

Max. first spine angle against horizontal plane.

Definition at line 99 of file barrel.cc.

Referenced by setup_dependent_parameters_and_sanity_check(), and spine_function().

◆ Alpha_min

double GlobalParameters::Alpha_min = MathematicalConstants::Pi/2.0*1.5

Min. spine angle against horizontal plane.

Min. first spine angle against horizontal plane.

Definition at line 96 of file barrel.cc.

Referenced by setup_dependent_parameters_and_sanity_check(), and spine_function().

◆ Beta_max

double GlobalParameters::Beta_max = MathematicalConstants::Pi/2.0

Max. second pine angle against horizontal plane.

Definition at line 116 of file common_young_laplace_stuff.h.

◆ Beta_min

double GlobalParameters::Beta_min = MathematicalConstants::Pi/2.0

Min. second spine angle against horizontal plane.

Definition at line 113 of file common_young_laplace_stuff.h.

◆ Case

int GlobalParameters::Case = All_pinned

What case are we considering: Choose one from the enumeration Cases.

Definition at line 53 of file common_young_laplace_stuff.h.

Referenced by RefineableYoungLaplaceProblem< ELEMENT >::actions_after_adapt(), main(), and setup_dependent_parameters_and_sanity_check().

◆ Control_element

unsigned GlobalParameters::Control_element = 0

Number of element in bulk mesh at which height control is applied. Initialise to 0 – will be overwritte in setup_dependent_parameters_and_sanity_check()

Definition at line 90 of file common_young_laplace_stuff.h.

Referenced by setup_dependent_parameters_and_sanity_check().