A numerical model of hydraulic resistance by hanging aquaculture facilities is developed and applied to a model basin and a field. A drag stress term formulated by the quadratic law of drag force is introduced Tn the equations of motion for a two-dimensional depth-averaged flow. In the model basin, numerical experiments ave tarried out for the various shape of obstructions, string density and layout of facilities etc.. The flow pattern around the facilities is affected sensitively by the density of string and the layout of facilities. On the other hand, the velocity decay due to the hanging oyster aquaculture facilities is observed in Kamak bay, where the maximum velcocity decay rate is $25\%$ in spring tide. The model is also applied to the field, Kamak bay. The velocity decay rate in the model is comparable with the field measurement data. The velocity decreases around the down-stream area of the facilities, .hut it increases in the other region. The water elevation decreases during the flood and it increases during the ebb.