Resting lake scenario for the shallow water equations. More...

Inheritance diagram for scenarios.swe_resting_lake.SWERestingLake:

Collaboration diagram for scenarios.swe_resting_lake.SWERestingLake:

Public Member Functions
	__init__ (self, initial_conditions="sinusoidal")
	initial_conditions (self)
	analytical_solution (self)
Public Member Functions inherited from scenarios.scenario.Scenario
	boundary_conditions (self)
	refinement_criterion (self)
	set_global_simulation_parameters (self, project)

Protected Member Functions
	_initial_conditions (self, Q)

Protected Attributes
	_scenario = initial_conditions

Additional Inherited Members
Static Protected Attributes inherited from scenarios.scenario.Scenario
	_equation = None
int	_dimensions = 2
float	_end_time = 1.0
float	_plot_dt = 0.1
float	_offset = 0.0
float	_domain_size = 1.0
bool	_periodic_bc = False

Detailed Description

Resting lake scenario for the shallow water equations.

The real water height H as defined by the sum of the water column h and the bathymetry b is constant over the entire domain, meaning that there should be no changes on the entire domain, but because we use the sum of the derivatives of h and b (h' + b') instead of the derivative of the sum (h + b)' some rounding errors will creep in, which causes unphysical waves to appear. As such this scenario is nice for testing how large these unphysical waves are for a given algorithm, and how stable the algorithm is, i.e. does it dampen out these waves or do they oscillate out of control.

Definition at line 11 of file swe_resting_lake.py.