/////////////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////////////// More...
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: Pass the number of elements and the width of the domain in the x direction. Also pass the number of elements in the y direction of the bottom (fluid 1) and top (fluid 2) layers, along with the heights of both layers. 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 173 of file elastic_two_layer_interface.cc.
Constructor & Destructor Documentation
◆ InterfaceProblem() [1/2]
| InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem |
Constructor.
Constructor for two fluid interface problem.
Definition at line 259 of file elastic_two_layer_interface.cc.
References Problem_Parameter::Constitutive_law_pt, Global_Physical_Variables::Density_Ratio, Global_Physical_Variables::G(), Global_Physical_Variables::Nu, Global_Physical_Variables::Re, Global_Physical_Variables::ReInvFr, Global_Physical_Variables::ReSt, and Global_Physical_Variables::Viscosity_Ratio.
◆ ~InterfaceProblem() [1/2]
|
inline |
Destructor (empty)
Definition at line 182 of file elastic_two_layer_interface.cc.
◆ InterfaceProblem() [2/2]
| 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 | ||
| ) |
Constructor: Pass the number of elements and the width of the domain in the x direction. Also pass the number of elements in the y direction of the bottom (fluid 1) and top (fluid 2) layers, along with the heights of both layers.
Constructor for two fluid interface problem.
Definition at line 170 of file spine_two_layer_interface.cc.
References InterfaceProblem< ELEMENT, TIMESTEPPER >::Bulk_mesh_pt, Global_Physical_Variables::Ca, Global_Physical_Variables::Density_Ratio, InterfaceProblem< ELEMENT, TIMESTEPPER >::fix_pressure(), Global_Physical_Variables::G(), Global_Physical_Variables::Re, Global_Physical_Variables::ReInvFr, Global_Physical_Variables::ReSt, InterfaceProblem< ELEMENT, TIMESTEPPER >::set_boundary_conditions(), Global_Physical_Variables::St, InterfaceProblem< ELEMENT, TIMESTEPPER >::Surface_mesh_pt, and Global_Physical_Variables::Viscosity_Ratio.
◆ ~InterfaceProblem() [2/2]
|
inline |
Destructor (empty)
Definition at line 105 of file spine_two_layer_interface.cc.
Member Function Documentation
◆ actions_after_adapt()
|
private |
Rebuild the mesh of interface elements after adapting the bulk mesh.
Definition at line 518 of file elastic_two_layer_interface.cc.
◆ actions_after_newton_solve() [1/2]
|
inlineprivate |
No actions required after solve step.
Definition at line 202 of file elastic_two_layer_interface.cc.
◆ actions_after_newton_solve() [2/2]
|
inlineprivate |
Update after solve can remain empty, because the update is performed automatically after every Newton step.
Definition at line 135 of file spine_two_layer_interface.cc.
◆ actions_before_adapt()
|
private |
Strip off the interface elements before adapting the bulk mesh.
Definition at line 502 of file elastic_two_layer_interface.cc.
◆ 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 206 of file elastic_two_layer_interface.cc.
◆ 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 125 of file spine_two_layer_interface.cc.
◆ actions_before_newton_solve() [1/2]
|
inlineprivate |
No actions required before solve step.
Definition at line 199 of file elastic_two_layer_interface.cc.
◆ actions_before_newton_solve() [2/2]
|
inlineprivate |
No actions required before solve step.
Definition at line 131 of file spine_two_layer_interface.cc.
◆ 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 558 of file elastic_two_layer_interface.cc.
References Global_Physical_Variables::Ca, Global_Physical_Variables::St, and Global_Physical_Variables::Viscosity_Ratio.
◆ deform_free_surface() [1/2]
|
private |
Deform the mesh/free surface to a prescribed function.
Definition at line 641 of file elastic_two_layer_interface.cc.
◆ deform_free_surface() [2/2]
|
private |
Deform the mesh/free surface to a prescribed function.
◆ delete_interface_elements()
|
private |
Delete the 1d interface elements.
Delete the interface elements and wipe the surface mesh.
Definition at line 619 of file elastic_two_layer_interface.cc.
◆ doc_solution() [1/2]
| void InterfaceProblem< ELEMENT, TIMESTEPPER >::doc_solution | ( | DocInfo & | doc_info | ) |
Document the solution.
Doc the solution.
Definition at line 670 of file elastic_two_layer_interface.cc.
References Problem_Parameter::Trace_file.
◆ doc_solution() [2/2]
| void InterfaceProblem< ELEMENT, TIMESTEPPER >::doc_solution | ( | DocInfo & | doc_info | ) |
Doc the solution.
◆ fix_pressure() [1/2]
|
inlineprivate |
Fix pressure in element e at pressure dof pdof and set to pvalue.
Definition at line 227 of file elastic_two_layer_interface.cc.
Referenced by InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem().
◆ fix_pressure() [2/2]
|
inlineprivate |
Fix pressure in element e at pressure dof pdof and set to pvalue.
Definition at line 141 of file spine_two_layer_interface.cc.
◆ set_boundary_conditions() [1/2]
| 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 467 of file elastic_two_layer_interface.cc.
Referenced by InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem().
◆ set_boundary_conditions() [2/2]
| void InterfaceProblem< ELEMENT, TIMESTEPPER >::set_boundary_conditions | ( | ) |
Set boundary conditions.
◆ set_initial_condition() [1/2]
| 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 437 of file elastic_two_layer_interface.cc.
◆ set_initial_condition() [2/2]
| void InterfaceProblem< ELEMENT, TIMESTEPPER >::set_initial_condition | ( | ) |
Set initial conditions.
◆ unsteady_run() [1/2]
| 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 722 of file elastic_two_layer_interface.cc.
References Problem_Parameter::Trace_file.
Referenced by main().
◆ unsteady_run() [2/2]
| 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.
Member Data Documentation
◆ Bulk_mesh_pt [1/2]
|
private |
Pointer to the (specific) "bulk" mesh.
Definition at line 237 of file elastic_two_layer_interface.cc.
Referenced by InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem().
◆ Bulk_mesh_pt [2/2]
|
private |
Pointer to the specific bulk mesh.
Definition at line 157 of file spine_two_layer_interface.cc.
◆ Constitutive_law_pt
|
private |
Definition at line 243 of file elastic_two_layer_interface.cc.
◆ Lx
|
private |
Width of domain.
Definition at line 246 of file elastic_two_layer_interface.cc.
◆ Surface_mesh_pt
|
private |
Pointer to the "surface" mesh.
Pointer to the surface mesh.
Definition at line 240 of file elastic_two_layer_interface.cc.
Referenced by InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem().
◆ Trace_file
|
private |
Trace file.
Definition at line 249 of file elastic_two_layer_interface.cc.
The documentation for this class was generated from the following files: