抄録
In this paper, we develop new discrete method which applied the approach of the hybrid model to a seepage flow problem with phreatic surface. First, the formulation which introduced the potential continuity by penalty function in general weighting residual procedure is shown. Compatibility of the potential on the element boundary edge is approximately introduced using the penalty function. In this method, an independent potential field is assumed for each element, and the head and the hydraulic gradient of each element are used as the degrees of freedom. Therefore, the shape of the element is not limited and the node does not need to have the degree of freedom and is used only for recognition of the domain shape. Then, we developed the technique in which only the integral region changed without initial mesh division in order to decide the phreatic surface. The equivalent solution with the analytical solution is obtained when it examined the accuracy of numerical solution using this technique.
