Inheritance diagram for InterfaceProblem< ELEMENT, TIMESTEPPER >:

Public Member Functions
	InterfaceProblem (const unsigned &n_x, const unsigned &n_y, const double &l_x, const double &h)
	Constructor for single fluid free surface 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)
	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_y, const double &l_x, const double &h)

	~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...

Public Attributes
SingleLayerSpineMesh< ELEMENT > *	Bulk_mesh_pt
	The bulk fluid mesh, complete with spines and update information. 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	deform_free_surface (const double &epsilon, const unsigned &n_periods)
	Deform the mesh/free surface to a prescribed function. 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...

Private Attributes
ElasticRectangularQuadMesh< 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...

Detailed Description

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

Single fluid free surface problem in a rectangular domain which is periodic in the x direction

Definition at line 96 of file elastic_single_layer.cc.

Constructor & Destructor Documentation

◆ InterfaceProblem() [1/2]

template<class ELEMENT , class TIMESTEPPER >

InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem	(	const unsigned &	n_x,
		const unsigned &	n_y,
		const double &	l_x,
		const double &	h
	)

Constructor for single fluid free surface problem.

Constructor: Pass the number of elements and the lengths of the domain in the x and y directions (h is the height of the fluid layer i.e. the length of the domain in the y direction)

Definition at line 164 of file elastic_single_layer.cc.

References InterfaceProblem< ELEMENT, TIMESTEPPER >::Bulk_mesh_pt, Global_Physical_Variables::Ca, InterfaceProblem< ELEMENT, TIMESTEPPER >::Constitutive_law_pt, Global_Physical_Variables::G(), Global_Physical_Variables::Nu, Global_Physical_Variables::Re, Global_Physical_Variables::ReInvFr, Global_Physical_Variables::ReSt, InterfaceProblem< ELEMENT, TIMESTEPPER >::set_boundary_conditions(), Global_Physical_Variables::St, and InterfaceProblem< ELEMENT, TIMESTEPPER >::Surface_mesh_pt.

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

◆ ~InterfaceProblem() [1/2]

template<class ELEMENT, class TIMESTEPPER>

InterfaceProblem< ELEMENT, TIMESTEPPER >::~InterfaceProblem ( )

inline

Destructor (empty)

Definition at line 108 of file elastic_single_layer.cc.

◆ InterfaceProblem() [2/2]

template<class ELEMENT, class TIMESTEPPER>

InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem	(	const unsigned &	n_x,
		const unsigned &	n_y,
		const double &	l_x,
		const double &	h
	)

Constructor: Pass the number of elements and the lengths of the domain in the x and y directions (h is the height of the fluid layer i.e. the length of the domain in the y direction)

◆ ~InterfaceProblem() [2/2]

template<class ELEMENT, class TIMESTEPPER>

InterfaceProblem< ELEMENT, TIMESTEPPER >::~InterfaceProblem ( )

inline

Destructor (empty)

Definition at line 99 of file spine_single_layer.cc.

Member Function Documentation

◆ actions_after_newton_solve() [1/2]

template<class ELEMENT, class TIMESTEPPER>

void InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_after_newton_solve ( )

inlineprivate

No actions required after solve step.

Definition at line 128 of file elastic_single_layer.cc.

◆ actions_after_newton_solve() [2/2]

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 135 of file spine_single_layer.cc.

References Global_Physical_Variables::Ca, InterfaceProblem< ELEMENT, TIMESTEPPER >::deform_free_surface(), InterfaceProblem< ELEMENT, TIMESTEPPER >::doc_solution(), Global_Physical_Variables::G(), InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem(), Global_Physical_Variables::Re, Global_Physical_Variables::ReInvFr, Global_Physical_Variables::ReSt, InterfaceProblem< ELEMENT, TIMESTEPPER >::set_boundary_conditions(), InterfaceProblem< ELEMENT, TIMESTEPPER >::set_initial_condition(), Global_Physical_Variables::St, and InterfaceProblem< ELEMENT, TIMESTEPPER >::unsteady_run().

◆ actions_before_implicit_timestep()

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 132 of file elastic_single_layer.cc.

◆ actions_before_newton_convergence_check()

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 125 of file spine_single_layer.cc.

◆ actions_before_newton_solve() [1/2]

template<class ELEMENT, class TIMESTEPPER>

void InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_before_newton_solve ( )

inlineprivate

No actions required before solve step.

Definition at line 125 of file elastic_single_layer.cc.

◆ actions_before_newton_solve() [2/2]

template<class ELEMENT, class TIMESTEPPER>

void InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_before_newton_solve ( )

inlineprivate

No actions required before solve step.

Definition at line 131 of file spine_single_layer.cc.

◆ deform_free_surface() [1/2]

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 408 of file elastic_single_layer.cc.

References InterfaceProblem< ELEMENT, TIMESTEPPER >::Bulk_mesh_pt, and InterfaceProblem< ELEMENT, TIMESTEPPER >::Lx.

Referenced by InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_after_newton_solve(), InterfaceProblem< ELEMENT, TIMESTEPPER >::set_boundary_conditions(), and InterfaceProblem< ELEMENT, TIMESTEPPER >::unsteady_run().

◆ deform_free_surface() [2/2]

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.

◆ doc_solution() [1/2]

template<class ELEMENT, class TIMESTEPPER>

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

Doc the solution.

◆ doc_solution() [2/2]

template<class ELEMENT , class TIMESTEPPER >

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

Document the solution.

Doc the solution.

Definition at line 436 of file elastic_single_layer.cc.

References InterfaceProblem< ELEMENT, TIMESTEPPER >::Bulk_mesh_pt, InterfaceProblem< ELEMENT, TIMESTEPPER >::Surface_mesh_pt, InterfaceProblem< ELEMENT, TIMESTEPPER >::Trace_file, and InterfaceProblem< ELEMENT, TIMESTEPPER >::unsteady_run().

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

◆ set_boundary_conditions() [1/2]

template<class ELEMENT, class TIMESTEPPER>

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

Set boundary conditions.

◆ set_boundary_conditions() [2/2]

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 bottom (solid) wall and the horizontal component only to zero on the side (periodic) boundaries.

Definition at line 370 of file elastic_single_layer.cc.

References InterfaceProblem< ELEMENT, TIMESTEPPER >::Bulk_mesh_pt, and InterfaceProblem< ELEMENT, TIMESTEPPER >::deform_free_surface().

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

◆ set_initial_condition() [1/2]

template<class ELEMENT, class TIMESTEPPER>

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

Set initial conditions.

◆ set_initial_condition() [2/2]

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 340 of file elastic_single_layer.cc.

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

◆ unsteady_run() [1/2]

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.

◆ unsteady_run() [2/2]

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 489 of file elastic_single_layer.cc.

References InterfaceProblem< ELEMENT, TIMESTEPPER >::deform_free_surface(), InterfaceProblem< ELEMENT, TIMESTEPPER >::doc_solution(), InterfaceProblem< ELEMENT, TIMESTEPPER >::set_initial_condition(), and InterfaceProblem< ELEMENT, TIMESTEPPER >::Trace_file.

Referenced by InterfaceProblem< ELEMENT, TIMESTEPPER >::actions_after_newton_solve(), InterfaceProblem< ELEMENT, TIMESTEPPER >::doc_solution(), and main().

Member Data Documentation

◆ Bulk_mesh_pt [1/2]

template<class ELEMENT, class TIMESTEPPER>

SingleLayerSpineMesh<ELEMENT>* InterfaceProblem< ELEMENT, TIMESTEPPER >::Bulk_mesh_pt

The bulk fluid mesh, complete with spines and update information.

Definition at line 114 of file spine_single_layer.cc.

◆ Bulk_mesh_pt [2/2]

template<class ELEMENT, class TIMESTEPPER>

ElasticRectangularQuadMesh<ELEMENT>* InterfaceProblem< ELEMENT, TIMESTEPPER >::Bulk_mesh_pt

private

Pointer to the (specific) "bulk" mesh.

Definition at line 141 of file elastic_single_layer.cc.

Referenced by InterfaceProblem< ELEMENT, TIMESTEPPER >::deform_free_surface(), InterfaceProblem< ELEMENT, TIMESTEPPER >::doc_solution(), InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem(), and InterfaceProblem< ELEMENT, TIMESTEPPER >::set_boundary_conditions().

◆ Constitutive_law_pt

template<class ELEMENT, class TIMESTEPPER>

ConstitutiveLaw* InterfaceProblem< ELEMENT, TIMESTEPPER >::Constitutive_law_pt

private

Definition at line 147 of file elastic_single_layer.cc.

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

◆ Lx

template<class ELEMENT, class TIMESTEPPER>

double InterfaceProblem< ELEMENT, TIMESTEPPER >::Lx

private

Width of domain.

Definition at line 150 of file elastic_single_layer.cc.

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

◆ Surface_mesh_pt

template<class ELEMENT, class TIMESTEPPER>

Mesh * InterfaceProblem< ELEMENT, TIMESTEPPER >::Surface_mesh_pt

private

Pointer to the "surface" mesh.

The mesh that contains the free surface elements.

Definition at line 144 of file elastic_single_layer.cc.

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

◆ Trace_file

template<class ELEMENT, class TIMESTEPPER>

ofstream InterfaceProblem< ELEMENT, TIMESTEPPER >::Trace_file

private

Trace file.

Definition at line 153 of file elastic_single_layer.cc.

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

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

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Detailed Description

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

Constructor & Destructor Documentation

◆ InterfaceProblem() [1/2]

◆ ~InterfaceProblem() [1/2]

◆ InterfaceProblem() [2/2]

◆ ~InterfaceProblem() [2/2]

Member Function Documentation

◆ actions_after_newton_solve() [1/2]

◆ actions_after_newton_solve() [2/2]

◆ actions_before_implicit_timestep()

◆ actions_before_newton_convergence_check()

◆ actions_before_newton_solve() [1/2]

◆ actions_before_newton_solve() [2/2]

◆ deform_free_surface() [1/2]

◆ deform_free_surface() [2/2]

◆ doc_solution() [1/2]

◆ doc_solution() [2/2]

◆ set_boundary_conditions() [1/2]

◆ set_boundary_conditions() [2/2]

◆ set_initial_condition() [1/2]

◆ set_initial_condition() [2/2]

◆ unsteady_run() [1/2]

◆ unsteady_run() [2/2]

Member Data Documentation

◆ Bulk_mesh_pt [1/2]

◆ Bulk_mesh_pt [2/2]

◆ Constitutive_law_pt

◆ Lx

◆ Surface_mesh_pt

◆ Trace_file

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