Multiscale finite element method with JUMP vector for simulating Darcy velocity at the interface of media
A JUMP vector multiscale finite element method(MSFEM-J)for simulating Darcy flow velocity at the interface of media is proposed,which can solve two shortcomings of current methods such as finite element meth-od when simulating groundwater problems with interfaces.First,the current method needs to accurately dissect and characterize the heterogeneity of media to ensure the accuracy of solution,and consumes huge computational resources;Secondly,the current numerical method cannot ensure that the Darcy velocity at the interface of the medium obeys the refraction law,that is,the normal velocity is continuous and the tangential velocity is propor-tional to the permeability coefficient.MSFEM-J first significantly improves the calculation efficiency of the groundwater head and Darcy velocity processes through the multi-scale finite element(MSFEM)basis function;Then,the JUMP vector is constructed based on the refraction law,and the right end term of the weak side node of the interface is iteratively corrected until convergence,so that the Darcy velocity at the interface meets the re-fraction law.The numerical experiments show that the Darcy velocity simulated by MSFEM-J in this paper con-forms to the refraction law at the interface of the medium,and has continuity at other locations in the study area.The accuracy of this method is higher than that of the current method of fine dissection,and its calculation con-sumption is only about 2％of that of the same method.
groundwater numerical simulationmultiscale finite element methodDarcy velocitymedium inter-faceJUMP vector
万方数据
维普
