nd = 2
T = opts.T
nDTout = opts.nDTout
L = (75.0, 30.0)
g = 9.81
# PARAMETERS #
mannings = 0.02
domainRect = RectangularDomain(L=L)
if opts.structured:
domain = domainRect
else:
domainRect.writePoly("hump")
domain = PlanarStraightLineGraphDomain("hump")
domain.boundaryTags = domainRect.boundaryTags
#This is relevant just when use_second_order_NonFlatB_with_EV_stabilization=True
cE = 1
LUMPED_MASS_MATRIX = 0
bt = domain.boundaryTags
bt['front'] = bt['bottom']
bt['back'] = bt['top']
domain.writePoly("tank2d")
######################
##### BATHYMETRY #####
######################
def bathymetry_function(X):
nd=2
T=opts.T
nDTout=opts.nDTout
L=(75.0,30.0)
g = 9.81
# PARAMETERS #
mannings=0.02
domainRect = RectangularDomain(L=L)
if opts.structured:
domain=domainRect
else:
domainRect.writePoly("hump")
domain = PlanarStraightLineGraphDomain("hump")
domain.boundaryTags = domainRect.boundaryTags
#This is relevant just when use_second_order_NonFlatB_with_EV_stabilization=True
cE=1
LUMPED_MASS_MATRIX=0
bt = domain.boundaryTags
bt['front'] = bt['bottom']
bt['back'] = bt['top']
domain.writePoly("tank2d")
######################
##### BATHYMETRY #####
######################
def bathymetry_function(X):
import numpy as np
# DOMAIN AND MESH #
###################
L = (75.0, 30.0)
refinement = opts.refinement
rectangle = RectangularDomain(L=L)
# CREATE REFINEMENT #
nnx0 = 6
nnx = (nnx0 - 1) * (2**refinement) + 1
nny = old_div((nnx - 1), 2) + 1
he = old_div(L[0], float(nnx - 1))
if opts.structured:
domain = rectangle
else:
rectangle.writePoly("dam3bumps")
domain = PlanarStraightLineGraphDomain(fileprefix="dam3bumps")
domain.MeshOptions.triangleOptions = "pAq30Dena%f" % (0.5 * opts.he**2, )
nnx = None
nny = None
domain.MeshOptions.nnx = nnx
domain.MeshOptions.nny = nny
######################
##### BATHYMETRY #####
######################
def bathymetry_function(X):
x = X[0]
y = X[1]
def create_domain2D():
return PlanarStraightLineGraphDomain()
runCFL = 0.25
useSuperlu = True
triangleFlag = 1
reflecting_BCs = opts.reflecting_BCs
##################
##### DOMAIN #####
##################
nd = 2
domainRect = RectangularDomain(L=L)
if opts.structured:
domain = domainRect
else:
pathfile = os.path.dirname(os.path.abspath(__file__)) + "/" + "hump"
#domainRect.writePoly(pathfile)
domain = PlanarStraightLineGraphDomain(pathfile)
domain.boundaryTags = domainRect.boundaryTags
domain.polyfile = pathfile
genMesh = False
bt = domain.boundaryTags
bt['front'] = bt['bottom']
bt['back'] = bt['top']
if opts.structured:
domain.writePoly("tank2d")
################
##### MESH #####
################
nnx0 = 6
# DOMAIN AND MESH #
###################
L = (48.8, 26.5) # this is length in x direction and y direction
refinement = opts.refinement
rectangle = RectangularDomain(L=L, x=[0, -13.25, 0])
# CREATE REFINEMENT #
nnx0 = 6
nnx = (nnx0 - 1) * (2**refinement) + 1
nny = old_div((nnx - 1), 2) + 1
he = old_div(L[0], float(nnx - 1))
if opts.structured:
domain = rectangle
else:
rectangle.writePoly("reef")
domain = PlanarStraightLineGraphDomain(fileprefix="reef")
domain.MeshOptions.triangleOptions = "pAq30Dena%f" % (0.5 * opts.he**2, )
nnx = None
nny = None
###############################
# CONSTANTS NEEDED FOR SETUP #
###############################
g = 9.81
# stuff for solitary wave
h0 = 0.78
alpha = 0.4 # 0.5 * h0
xs = 5.0
r = np.sqrt(old_div(3. * alpha, (4. * h0**2 * (h0 + alpha))))
c = np.sqrt(g * (h0 + alpha))
###################
L = (25.0, 2.0)
refinement = opts.refinement
rectangle = RectangularDomain(L=L, x=[0, 0, 0])
# CREATE REFINEMENT #
nnx0 = 6
nnx = (nnx0 - 1) * (2**refinement) + 1
nny = old_div((nnx - 1), 10) + 1
he = old_div(L[0], float(nnx - 1))
triangleOptions = "pAq30Dena%f" % (0.5 * he**2, )
if opts.structured:
domain = rectangle
else:
rectangle.writePoly("solitary")
domain = PlanarStraightLineGraphDomain(fileprefix="solitary")
domain.MeshOptions.triangleOptions = "pAq30Dena%f" % (0.5 * opts.he**2, )
nnx = None
nny = None
domain.MeshOptions.nnx = nnx
domain.MeshOptions.nny = nny
##################################
# SOLITARY WAVE FUCTION AND BATH #
##################################
g = 9.81
# constants for solitary wave
h0 = 0.5
alpha = 0.2 * h0
xs = 5.0
r = np.sqrt(old_div(3. * alpha, (4. * h0**2 * (h0 + alpha))))
###################
L = (25.0, 1.0)
refinement = opts.refinement
rectangle = RectangularDomain(L=L, x=[0, 0, 0])
# CREATE REFINEMENT #
nnx0 = 6
nnx = (nnx0 - 1) * (2**refinement) + 1
nny = old_div((nnx - 1), 10) + 1
he = old_div(L[0], float(nnx - 1))
triangleOptions = "pAq30Dena%f" % (0.5 * he**2, )
if opts.structured:
domain = rectangle
else:
rectangle.writePoly("bump")
domain = PlanarStraightLineGraphDomain(fileprefix="bump")
domain.MeshOptions.triangleOptions = "pAq30Dena%f" % (0.5 * opts.he**2, )
nnx = None
nny = None
##########################################
# DEFINE INITIAL CONSTANTS AND FUNCTIONS #
##########################################
g = 9.81
# constants for transcritical bump
hL = 0.28205279813802181
q_in = 0.18
def bathymetry_function(X):
###################
# DOMAIN AND MESH #
###################
L = (20.0, 20.0) # this is length in x direction and y direction
refinement = opts.refinement
rectangle = RectangularDomain(L=L)
# CREATE REFINEMENT #
nnx = 50
nny = 50
he = old_div(L[0], float(nnx - 1))
if opts.structured:
domain = rectangle
else:
rectangle.writePoly("circular_damBreak")
domain = PlanarStraightLineGraphDomain(fileprefix="circular_damBreak")
domain.MeshOptions.triangleOptions = "pAq30Dena%f" % (0.5 * opts.he**2, )
nnx = None
nny = None
######################
##### BATHYMETRY #####
######################
def bathymetry_function(X):
x = X[0]
return x * 0.
##############################
# DOMAIN AND MESH #
###################
L = (15.0, 2.0) # this is domain length in x direction and y direction
refinement = opts.refinement
rectangle = RectangularDomain(L=L)
# CREATE REFINEMENT #
nnx0 = 6
nnx = (nnx0 - 1) * (2**refinement) + 1
nny = old_div((nnx - 1), 10) + 1
he = old_div(L[0], float(nnx - 1))
if opts.structured:
domain = rectangle
else:
rectangle.writePoly("seawall")
domain = PlanarStraightLineGraphDomain(fileprefix="seawall")
domain.MeshOptions.triangleOptions = "pAq30Dena%f" % (0.5 * opts.he**2,)
nnx = None
nny = None
###############################
# CONSTANTS NEEDED FOR SETUP #
###############################
g = 9.81
h0 = 0.2 # water depth
a = 0.35 # relative amplitude
k_wavenumber = np.sqrt(3.0 * a / (4.0 * h0**3)) # wavenumber
z = np.sqrt(3.0 * a * h0) / (2.0 * h0 * np.sqrt(h0 * (1.0 + a))) # width
c = np.sqrt(g * (1.0 + a) * h0) # wave celerity
x0 = 5.9 # initial location of solitary wave
# DOMAIN AND MESH #
###################
L = (25.0, 30.0) # this is length in x direction and y direction
refinement = opts.refinement
rectangle = RectangularDomain(L=L)
# CREATE REFINEMENT #
nnx0 = 6
nnx = (nnx0 - 1) * (2**refinement) + 1
nny = old_div((nnx - 1), 2) + 1
he = old_div(L[0], float(nnx - 1))
if opts.structured:
domain = rectangle
else:
rectangle.writePoly("island")
domain = PlanarStraightLineGraphDomain(fileprefix="island")
domain.MeshOptions.triangleOptions = "pAq30Dena%f" % (0.5 * opts.he**2, )
nnx = None
nny = None
###############################
# CONSTANTS NEEDED FOR SETUP #
###############################
g = 9.81
# stuff for solitary wave
h0 = 0.32
alpha = 0.181 * h0
xs = 7.56023
if opts.which_case == 1:
alpha = 0.091
g = 9.81
# PARAMETERS #
mannings = 0.02
# DOMAIN #
vertices = [
[0.0, 0.0], #0
[75.0, 0.0], #1
[90.0, 15.0], #2
[75.0, 30.0], #3
[0.0, 30.0]
] #4
segments = [[0, 1], [1, 2], [2, 3], [3, 4], [4, 0]]
domain = PlanarStraightLineGraphDomain(vertices=vertices, segments=segments)
domain.writePoly("tank2d")
#This is relevant just when use_second_order_NonFlatB_with_EV_stabilization=True
cE = 1
LUMPED_MASS_MATRIX = 0
######################
##### BATHYMETRY #####
######################
def bathymetry_function(X):
import numpy as np
x = X[0]
y = X[1]
bump1 = 1 - 1. / 8 * np.sqrt((x - 30)**2 + (y - 6)**2)
###################
L = (30.0, 25.0) # this is length in x direction and y direction
refinement = opts.refinement
rectangle = RectangularDomain(L=L)
# CREATE REFINEMENT #
nnx0 = 6
nnx = (nnx0 - 1) * (2**refinement) + 1
nny = old_div((nnx - 1), 2) + 1
he = old_div(L[0], float(nnx - 1))
triangleOptions = "pAq30Dena%f" % (0.5 * he**2, )
if opts.structured:
domain = rectangle
else:
rectangle.writePoly("well_balancing")
domain = PlanarStraightLineGraphDomain(fileprefix="well_balancing")
domain.MeshOptions.triangleOptions = "pAq30Dena%f" % (0.5 * opts.he**2, )
nnx = None
nny = None
###############################
# CONSTANTS NEEDED FOR SETUP #
###############################
g = 9.81
h0 = 1.0
# stuff for cone bathymetry
htop = 0.625
rcone = 3.6
scone = 4.0
hcone = 0.9
# PARAMETERS #
mannings=0.02
# DOMAIN #
vertices=[[0.0, 0.0], #0
[75.0, 0.0], #1
[90.0, 15.0], #2
[75.0, 30.0], #3
[0.0, 30.0]] #4
segments=[[0, 1],
[1, 2],
[2, 3],
[3, 4],
[4, 0]]
domain = PlanarStraightLineGraphDomain(vertices=vertices,
segments=segments)
domain.writePoly("tank2d")
#This is relevant just when use_second_order_NonFlatB_with_EV_stabilization=True
cE=1
LUMPED_MASS_MATRIX=0
######################
##### BATHYMETRY #####
######################
def bathymetry_function(X):
import numpy as np
x = X[0]
y = X[1]
bump1 = 1-1./8*np.sqrt((x-30)**2+(y-6)**2)
bump2 = 1-1./8*np.sqrt((x-30)**2+(y-24)**2)
reflecting_BCs = opts.reflecting_BCs
##################
##### DOMAIN #####
##################
nd = 2
Lx = 90
vertices = [
[0.0, 0.0], #0
[75.0, 0.0], #1
[90.0, 15.0], #2
[75.0, 30.0], #3
[0.0, 30.0]
] #4
segments = [[0, 1], [1, 2], [2, 3], [3, 4], [4, 0]]
domain = PlanarStraightLineGraphDomain(vertices=vertices, segments=segments)
#domain.writePoly("tank2d")
genMesh = False
pathfile = os.path.dirname(os.path.abspath(__file__)) + "/" + "tank2d"
domain.polyfile = pathfile
#domain.writePoly(domain.polyfile)
################
##### MESH #####
################
nnx0 = 6
nnx = (nnx0 - 1) * (2**refinement) + 1
nny = old_div((nnx - 1), 2) + 1
nnz = 1
he = old_div(Lx, float(nnx - 1))