26 #ifndef OOMPH_COMMON_YOUNG_LAPLACE_STUFF_DOC
27 #define OOMPH_COMMON_YOUNG_LAPLACE_STUFF_DOC
60 double Gamma = MathematicalConstants::Pi/4.0;
107 double Alpha_min = MathematicalConstants::Pi/2.0;
110 double Alpha_max = MathematicalConstants::Pi/2.0;
133 if ((
N_x%2!=0)||(
N_y%2!=0))
135 cout <<
"n_x n_y should even" << endl;
152 Alpha_max=MathematicalConstants::Pi/2.0*0.5;
158 Alpha_min=MathematicalConstants::Pi/2.0*1.5;
159 Alpha_max=MathematicalConstants::Pi/2.0*0.5;
165 Alpha_min=MathematicalConstants::Pi/2.0*1.5;
166 Alpha_max=MathematicalConstants::Pi/2.0*0.5;
172 Alpha_min=MathematicalConstants::Pi/2.0*1.5;
173 Alpha_max=MathematicalConstants::Pi/2.0*0.5;
178 std::cout <<
"Never get here: Case = " <<
Case << std::endl;
196 Vector<double>& spine_B,
197 Vector< Vector<double> >& dspine_B)
204 dspine_B[0][0] = 1.0 ;
205 dspine_B[1][0] = 0.0 ;
206 dspine_B[0][1] = 0.0 ;
207 dspine_B[1][1] = 1.0 ;
208 dspine_B[0][2] = 0.0 ;
209 dspine_B[1][2] = 0.0 ;
218 Vector<double>& spine,
219 Vector< Vector<double> >& dspine)
223 Vector<double> x(2,0.0);
290 std::cout <<
"No exact solution for this case..." << std::endl;
Namespace for "global" problem parameters.
double Alpha_max
Max. spine angle against horizontal plane.
double L_x
Length and width of the domain.
double Controlled_height
Height control value.
unsigned Control_element
Number of element in bulk mesh at which height control is applied. Initialise to 0 – will be overwrit...
double get_exact_kappa()
Exact kappa.
bool Use_height_control
Use height control (true) or not (false)?
double Controlled_height_increment
Increment for height control.
bool Rotate_spines_in_both_directions
Should the spines rotate in the x and y directions (true)?
double Kappa_increment
Increment for prescribed curvature.
bool Use_spines
Use spines (true) or not (false)
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...
double Beta_max
Max. second pine angle against horizontal plane.
unsigned Nsteps
Number of steps.
unsigned N_x
Number of elements in the mesh.
Data * Kappa_pt
Pointer to Data object that stores the prescribed curvature.
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 a...
double L_y
Width of domain.
double Gamma
Contact angle and its cos (dependent parameter – is reassigned)
Cases
Enumeration for the possible cases.
@ T_junction_with_nonzero_contact_angle
@ Spherical_cap_in_cylinder_pinned
double Kappa_initial
Initial value for kappa.
double Alpha_min
Min. spine angle against horizontal plane.
double Cos_gamma
Cos of contact angle.
int Case
What case are we considering: Choose one from the enumeration Cases.
double Beta_min
Min. second spine angle against horizontal plane.
int Step_sign
Increase or decrease the value of the control parameters?
void setup_dependent_parameters_and_sanity_check()
Setup dependent parameters and perform sanity check.