Analytic solutions

class freshkiss3d.AnalyticSol
Parameters
  • triangular_mesh (freshkiss3d.mesh.TriangularMesh) –

  • scheduler (freshkiss3d.Scheduler) – Analytic solution calculus scheduler (default = fk.AlwaysScheduler)

  • compute_error (bool) – Error computation selector (default = True)

  • error_type (string ('L2', 'L1')) – Error type selector (default = ‘L2’)

  • error_output (string ('none', 'terminal', 'txt')) – Error output type (default = ‘none’)

topography

Ground elevation

Type

array, float [NC]

H

Analatic water height

Type

array, float [NC]

elevation

z coordonate of the free surface = H + topography

Type

array, float [NC]

U

x analitic velocity

Type

array, float [NC,1]

V

y analitic velocity

Type

array, float [NC,1]

error_H
Type

float

error_U
Type

float

error_V
Type

float

class freshkiss3d.Bump

Bumps 1d stationary analytical solutions:

Parameters
  • args (AnalyticSol) –

  • case (string) – Flow type

  • q_in (float) – inlet fluvial flowrate

  • h_out (float) – outlet fluvial height

  • free_surface_0 (float) – Initial free surface

  • L (float) – Domain length (default: 25.)

  • x_b (float) – center of the bump (default: 10.)

  • topo (function) – Topography

z_m

Max elevation of topography

Type

float

h_c

Critical height

Type

float

compute_solution()

Computes the analytic solution at a given time.

Analytical solution can be found in: https://hal.archives-ouvertes.fr/hal-00628246/document section 3.1 (Bumps)

Parameters

time – Simulation time (usually Simutime.Time)

Examples using freshkiss3d.Bump

Bump

Bump

Bump
Bump mesh convergence

Bump mesh convergence

Bump mesh convergence
class freshkiss3d.DamBreak

Dam breaks 1d analytical solutions:

Parameters
  • args (AnalyticSol) –

  • h_l (float) – Left side height (Riemann problem)

  • h_r (float) – Right side height (Riemann problem)

  • x_d (float) – Dam position (default: 5.)

  • L (float) – Domain length (default: 10.)

compute_solution()

Computes the analytic solution at a given time.

Analytical solution can be found in: https://hal.archives-ouvertes.fr/hal-00628246/document section 4.1 (Dam breaks)

Parameters

time (float) – Simulation time (usually Simutime.Time)

Examples using freshkiss3d.DamBreak

Dam break

Dam break

Dam break
class freshkiss3d.Thacker2D

Thacker 2d analytical solution:

Parameters
compute_solution()

Computes the analytic solution at a given time.

Analytical solution can be found in: https://hal.archives-ouvertes.fr/hal-00628246/document section 4.2.2 (Planar surface in a paraboloid)

Parameters

time (float) – Simulation time (usually Simutime.Time)

Examples using freshkiss3d.Thacker2D

Thacker 2d

Thacker 2d

Thacker 2d
Thacker2d mesh convergence

Thacker2d mesh convergence

Thacker2d mesh convergence
Analytical solutions

Analytical solutions

Analytical solutions