def turbine_drag_function():
   from firedrake_fluids.shallow_water import ShallowWater
   
   sw = ShallowWater(path=os.path.join(cwd, "bump.swml"))
   array = sw.get_turbine_array()
   bump = assemble(array.turbine_drag*dx)

   sw = ShallowWater(path=os.path.join(cwd, "tophat.swml"))
   array = sw.get_turbine_array()
   tophat = assemble(array.turbine_drag*dx)
   
   return bump, tophat
Example #2
0
def turbine_drag_function():
    from firedrake_fluids.shallow_water import ShallowWater

    sw = ShallowWater(path=os.path.join(cwd, "bump.swml"))
    array = sw.get_turbine_array()
    bump = assemble(array.turbine_drag() * dx)

    sw = ShallowWater(path=os.path.join(cwd, "tophat.swml"))
    array = sw.get_turbine_array()
    tophat = assemble(array.turbine_drag() * dx)

    return bump, tophat