Public Member Functions | Private Member Functions | Private Attributes | List of all members
InterfaceProblem< ELEMENT, TIMESTEPPER > Class Template Reference
Inheritance diagram for InterfaceProblem< ELEMENT, TIMESTEPPER >:

Public Member Functions

 InterfaceProblem ()
 Constructor. More...
 
 ~InterfaceProblem ()
 Destructor (empty) More...
 
void set_initial_condition ()
 Set initial conditions. More...
 
void set_boundary_conditions ()
 Set boundary conditions. More...
 
void doc_solution (DocInfo &doc_info)
 Document the solution. More...
 
void unsteady_run (const double &t_max, const double &dt)
 Do unsteady run up to maximum time t_max with given timestep dt. More...
 
 InterfaceProblem (const unsigned &n_x, const unsigned &n_y1, const unsigned &n_y2, const double &l_x, const double &h1, const double &h2)
 Constructor for two fluid interface problem. More...
 
 ~InterfaceProblem ()
 Destructor (empty) More...
 
void set_initial_condition ()
 Set initial conditions. More...
 
void set_boundary_conditions ()
 Set boundary conditions. More...
 
void doc_solution (DocInfo &doc_info)
 Doc the solution. More...
 
void unsteady_run (const double &t_max, const double &dt)
 Do unsteady run up to maximum time t_max with given timestep dt. More...
 

Private Member Functions

void actions_before_newton_solve ()
 No actions required before solve step. More...
 
void actions_after_newton_solve ()
 No actions required after solve step. More...
 
void actions_before_implicit_timestep ()
 Actions before the timestep: For maximum stability, reset the current nodal positions to be the "stress-free" ones. More...
 
void actions_before_adapt ()
 Strip off the interface elements before adapting the bulk mesh. More...
 
void actions_after_adapt ()
 Rebuild the mesh of interface elements after adapting the bulk mesh. More...
 
void create_interface_elements ()
 Create the 1d interface elements. More...
 
void delete_interface_elements ()
 Delete the 1d interface elements. More...
 
void deform_free_surface (const double &epsilon, const unsigned &n_periods)
 Deform the mesh/free surface to a prescribed function. More...
 
void fix_pressure (const unsigned &e, const unsigned &pdof, const double &pvalue)
 Fix pressure in element e at pressure dof pdof and set to pvalue. More...
 
void actions_before_newton_convergence_check ()
 Spine heights/lengths are unknowns in the problem so their values get corrected during each Newton step. However, changing their value does not automatically change the nodal positions, so we need to update all of them here. More...
 
void actions_before_newton_solve ()
 No actions required before solve step. More...
 
void actions_after_newton_solve ()
 Update after solve can remain empty, because the update is performed automatically after every Newton step. More...
 
void deform_free_surface (const double &epsilon, const unsigned &n_periods)
 Deform the mesh/free surface to a prescribed function. More...
 
void fix_pressure (const unsigned &e, const unsigned &pdof, const double &pvalue)
 Fix pressure in element e at pressure dof pdof and set to pvalue. More...
 

Private Attributes

ElasticRefineableTwoLayerMesh
< ELEMENT > * 
Bulk_mesh_pt
 Pointer to the (specific) "bulk" mesh. More...
 
Mesh * Surface_mesh_pt
 Pointer to the "surface" mesh. More...
 
ConstitutiveLaw * Constitutive_law_pt
 
double Lx
 Width of domain. More...
 
ofstream Trace_file
 Trace file. More...
 
TwoLayerSpineMesh< ELEMENT > * Bulk_mesh_pt
 Pointer to the specific bulk mesh. More...
 

Detailed Description

template<class ELEMENT, class TIMESTEPPER>
class InterfaceProblem< ELEMENT, TIMESTEPPER >

Two fluid interface problem in a rectangular domain which is periodic in the x direction

Definition at line 177 of file elastic_two_layer_interface.cc.

Constructor & Destructor Documentation

template<class ELEMENT , class TIMESTEPPER >
InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem ( )
template<class ELEMENT, class TIMESTEPPER>
InterfaceProblem< ELEMENT, TIMESTEPPER >::~InterfaceProblem ( )
inline

Destructor (empty)

Definition at line 186 of file elastic_two_layer_interface.cc.

template<class ELEMENT , class TIMESTEPPER >
InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem ( const unsigned &  n_x,
const unsigned &  n_y1,
const unsigned &  n_y2,
const double &  l_x,
const double &  h1,
const double &  h2 
)
template<class ELEMENT, class TIMESTEPPER>
InterfaceProblem< ELEMENT, TIMESTEPPER >::~InterfaceProblem ( )
inline

Destructor (empty)

Definition at line 109 of file spine_two_layer_interface.cc.

Member Function Documentation

template<class ELEMENT , class TIMESTEPPER >
void InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_after_adapt ( )
private

Rebuild the mesh of interface elements after adapting the bulk mesh.

Definition at line 522 of file elastic_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_after_newton_solve ( )
inlineprivate

Update after solve can remain empty, because the update is performed automatically after every Newton step.

Definition at line 139 of file spine_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_after_newton_solve ( )
inlineprivate

No actions required after solve step.

Definition at line 206 of file elastic_two_layer_interface.cc.

template<class ELEMENT , class TIMESTEPPER >
void InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_before_adapt ( )
private

Strip off the interface elements before adapting the bulk mesh.

Definition at line 506 of file elastic_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_before_implicit_timestep ( )
inlineprivate

Actions before the timestep: For maximum stability, reset the current nodal positions to be the "stress-free" ones.

Definition at line 210 of file elastic_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_before_newton_convergence_check ( )
inlineprivate

Spine heights/lengths are unknowns in the problem so their values get corrected during each Newton step. However, changing their value does not automatically change the nodal positions, so we need to update all of them here.

Definition at line 129 of file spine_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_before_newton_solve ( )
inlineprivate

No actions required before solve step.

Definition at line 135 of file spine_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_before_newton_solve ( )
inlineprivate

No actions required before solve step.

Definition at line 203 of file elastic_two_layer_interface.cc.

template<class ELEMENT , class TIMESTEPPER >
void InterfaceProblem< ELEMENT, TIMESTEPPER >::create_interface_elements ( )
private

Create the 1d interface elements.

Create interface elements between the two fluids in the mesh pointed to by Bulk_mesh_pt and add the elements to the Mesh object pointed to by Surface_mesh_pt.

Definition at line 562 of file elastic_two_layer_interface.cc.

References Global_Physical_Variables::Ca, Global_Physical_Variables::St, and Global_Physical_Variables::Viscosity_Ratio.

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::deform_free_surface ( const double &  epsilon,
const unsigned &  n_periods 
)
private

Deform the mesh/free surface to a prescribed function.

template<class ELEMENT , class TIMESTEPPER >
void InterfaceProblem< ELEMENT, TIMESTEPPER >::deform_free_surface ( const double &  epsilon,
const unsigned &  n_periods 
)
private

Deform the mesh/free surface to a prescribed function.

Definition at line 646 of file elastic_two_layer_interface.cc.

template<class ELEMENT , class TIMESTEPPER >
void InterfaceProblem< ELEMENT, TIMESTEPPER >::delete_interface_elements ( )
private

Delete the 1d interface elements.

Delete the interface elements and wipe the surface mesh.

Definition at line 623 of file elastic_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::doc_solution ( DocInfo &  doc_info)

Doc the solution.

template<class ELEMENT , class TIMESTEPPER >
void InterfaceProblem< ELEMENT, TIMESTEPPER >::doc_solution ( DocInfo &  doc_info)

Document the solution.

Doc the solution.

Definition at line 674 of file elastic_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::fix_pressure ( const unsigned &  e,
const unsigned &  pdof,
const double &  pvalue 
)
inlineprivate

Fix pressure in element e at pressure dof pdof and set to pvalue.

Definition at line 145 of file spine_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::fix_pressure ( const unsigned &  e,
const unsigned &  pdof,
const double &  pvalue 
)
inlineprivate

Fix pressure in element e at pressure dof pdof and set to pvalue.

Definition at line 231 of file elastic_two_layer_interface.cc.

Referenced by InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem().

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::set_boundary_conditions ( )

Set boundary conditions.

template<class ELEMENT , class TIMESTEPPER >
void InterfaceProblem< ELEMENT, TIMESTEPPER >::set_boundary_conditions ( )

Set boundary conditions.

Set boundary conditions: Set both velocity components to zero on the top and bottom (solid) walls and the horizontal component only to zero on the side (periodic) boundaries.

Definition at line 471 of file elastic_two_layer_interface.cc.

Referenced by InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem().

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::set_initial_condition ( )

Set initial conditions.

template<class ELEMENT , class TIMESTEPPER >
void InterfaceProblem< ELEMENT, TIMESTEPPER >::set_initial_condition ( )

Set initial conditions.

Set initial conditions: Set all nodal velocities to zero and initialise the previous velocities and nodal positions to correspond to an impulsive start.

Definition at line 441 of file elastic_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
void InterfaceProblem< ELEMENT, TIMESTEPPER >::unsteady_run ( const double &  t_max,
const double &  dt 
)

Do unsteady run up to maximum time t_max with given timestep dt.

template<class ELEMENT , class TIMESTEPPER >
void InterfaceProblem< ELEMENT, TIMESTEPPER >::unsteady_run ( const double &  t_max,
const double &  dt 
)

Do unsteady run up to maximum time t_max with given timestep dt.

Perform run up to specified time t_max with given timestep dt.

Definition at line 727 of file elastic_two_layer_interface.cc.

Referenced by main().

Member Data Documentation

template<class ELEMENT, class TIMESTEPPER>
TwoLayerSpineMesh<ELEMENT>* InterfaceProblem< ELEMENT, TIMESTEPPER >::Bulk_mesh_pt
private

Pointer to the specific bulk mesh.

Definition at line 161 of file spine_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
ElasticRefineableTwoLayerMesh<ELEMENT>* InterfaceProblem< ELEMENT, TIMESTEPPER >::Bulk_mesh_pt
private

Pointer to the (specific) "bulk" mesh.

Definition at line 241 of file elastic_two_layer_interface.cc.

Referenced by InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem().

template<class ELEMENT, class TIMESTEPPER>
ConstitutiveLaw* InterfaceProblem< ELEMENT, TIMESTEPPER >::Constitutive_law_pt
private

Definition at line 247 of file elastic_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
double InterfaceProblem< ELEMENT, TIMESTEPPER >::Lx
private

Width of domain.

Definition at line 250 of file elastic_two_layer_interface.cc.

template<class ELEMENT, class TIMESTEPPER>
Mesh * InterfaceProblem< ELEMENT, TIMESTEPPER >::Surface_mesh_pt
private

Pointer to the "surface" mesh.

Pointer to the surface mesh.

Definition at line 244 of file elastic_two_layer_interface.cc.

Referenced by InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem().

template<class ELEMENT, class TIMESTEPPER>
ofstream InterfaceProblem< ELEMENT, TIMESTEPPER >::Trace_file
private

Trace file.

Definition at line 253 of file elastic_two_layer_interface.cc.


The documentation for this class was generated from the following files: