Bt = Transmissive_boundary(domain)
#Constant inflow
Bd = Dirichlet_boundary([inflow_stage, 0.0, 0.0])
#Time dependent inflow
from math import sin, pi
Bw = Time_boundary(domain=domain,
f=lambda x: [(1 + sin(x*pi/4))*\
(inflow_stage*(sin(2.5*x*pi)+0.7)),0,0])
#Set boundary conditions
domain.set_boundary({'left': Bd, 'right': Br, 'bottom': Br, 'top': Br})
######################
#Initial condition
#
log.critical('Initial condition')
domain.set_quantity('stage', Z)
#Evolve
for t in domain.evolve(yieldstep = 0.1, finaltime = 30):
domain.write_time(track_speeds=True)
domain.write_boundary_statistics(['stage'],'left')
log.critical('Done')
# Boundary conditions
#
log.critical('Boundaries')
Br = Reflective_boundary(domain)
Bt = Transmissive_boundary(domain)
#Constant inflow
Bd = Dirichlet_boundary([inflow_stage, 0.0, 0.0])
#Time dependent inflow
from math import sin, pi
Bw = Time_boundary(domain=domain,
f=lambda x: [(1 + sin(x*pi/4))*\
(inflow_stage*(sin(2.5*x*pi)+0.7)),0,0])
#Set boundary conditions
domain.set_boundary({'left': Bd, 'right': Br, 'bottom': Br, 'top': Br})
######################
#Initial condition
#
log.critical('Initial condition')
domain.set_quantity('stage', Z)
#Evolve
for t in domain.evolve(yieldstep=0.1, finaltime=30):
domain.write_time(track_speeds=True)
domain.write_boundary_statistics(['stage'], 'left')
log.critical('Done')
END\n\
2 ?? ??\n\
10.0 80.0\n\
10.0 90.0\n\
20.0 90.0\n\
10.0 80.0\n\
END\n\
END\n")
file.close()
# import the ungenerate file
m.import_ungenerate_file(fileName)
os.remove(fileName)
m.generate_mesh(maximum_triangle_area=max_area,verbose=False)
mesh_filename = "mesh.tsh"
m.export_mesh_file(mesh_filename)
# Run a simulation on the mesh
domain = Domain(mesh_filename, use_cache = False)
Br = Reflective_boundary(domain)
Bd = Dirichlet_boundary([3,0,0])
domain.set_boundary( {'wall': Br, 'wave': Bd} )
yieldstep = 0.1
finaltime = 20
for t in domain.evolve(yieldstep, finaltime):
domain.write_time()
