Changeset 8005


Ignore:
Timestamp:
Sep 4, 2010, 9:47:43 PM (9 years ago)
Author:
habili
Message:

using losses

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_core/source/anuga/structures/testing_wide_bridge.py

    r7995 r8005  
    166166domain.forcing_terms.append(culvert_energy)
    167167"""
    168 from anuga.structures.boyd_box_operator import Boyd_box_operator
     168#from anuga.structures.boyd_box_operator import Boyd_box_operator
    169169#culvert0 = Culvert_operator(domain,
    170170#                            end_point0=[40.0, 75.0],
     
    180180#------------------------------------------------------------------------------
    181181culverts = []
    182 number_of_culverts = 1
     182number_of_culverts = 1 
    183183for i in range(number_of_culverts):
    184184    culvert_width = 50.0/number_of_culverts
     
    186186    ep0 = [40.0, y]
    187187    ep1 = [50.0, y]
    188     culverts.append(Boyd_box_operator(domain,
     188    culverts.append(anuga.Boyd_box_operator(domain,
    189189                            end_point0=ep0,
    190190                            end_point1=ep1,
    191                             width=culvert_width,
    192                             height=10.0,
     191                            losses=1.5,
     192                            width=3.658,
     193                            height=3.658,
    193194                            apron=5.0,
    194195                            use_momentum_jet=True,
    195                             use_velocity_head=False,
     196                            use_velocity_head=True,
    196197                            manning=0.013,
     198                            description='bridge culvert',
    197199                            verbose=False))
    198200
     
    244246#Btds = Time_boundary(domain, lambda t: [0.0+ 0.0075*(1+num.sin(2*pi*(t-4)/20)), 0.0, 0.0])
    245247
    246 Btus = anuga.Dirichlet_boundary([18.0, 0, 0])           # Outflow water at 5.0
    247 Btds = anuga.Dirichlet_boundary([0.0, 0, 0])           # Outflow water at 1.0
     248Btus = anuga.Dirichlet_boundary([20.0, 0, 0])           # Outflow water at 5.0
     249Btds = anuga.Dirichlet_boundary([19.0, 0, 0])           # Outflow water at 1.0
    248250domain.set_boundary({'left': Btus, 'right': Btds, 'top': Br, 'bottom': Br})
    249251
     
    253255#------------------------------------------------------------------------------
    254256
    255 for t in domain.evolve(yieldstep = 1, finaltime = 100):
     257for t in domain.evolve(yieldstep = 1, finaltime = 10):
    256258    print domain.timestepping_statistics()
    257259    print domain.volumetric_balance_statistics()
    258     for culvert in culverts:
     260    for i, culvert in enumerate(culverts):
     261        print 'culvert: ', i
    259262        print culvert.structure_statistics()
    260263   
Note: See TracChangeset for help on using the changeset viewer.