In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,t...In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,the original Boussinesq system is transformed into an equivalent one.Then we discretize it using the second-order backward di erentiation formula(BDF2)and Crank-Nicolson(CN)to obtain two second-order time-advanced schemes.In both numerical schemes,a pressure-correction method is employed to decouple the velocity and pressure.These two schemes possess the desired property that they can be fully decoupled with satisfying unconditional stability.We rigorously prove both the unconditional stability and unique solvability of the discrete schemes.Furthermore,we provide detailed implementations of the decoupling procedures.Finally,various 2D numerical simulations are performed to verify the accuracy and energy stability of the proposed schemes.展开更多
In order to study the water flow in the drainage layer of highway under steady-state condition, one-dimensional （1D） Boussinesq equation-based model with Dupuit-Forchheimer assumption was established and the semi-an...In order to study the water flow in the drainage layer of highway under steady-state condition, one-dimensional （1D） Boussinesq equation-based model with Dupuit-Forchheimer assumption was established and the semi-analytical solutions to predict the water-table height were presented. In order to validate the model, a two-dimensional （2D） saturated flow model based on the Laplace equation was applied for the purpose of the model comparison. The water-table elevations predicted by ID Boussinesq equation-based model and 2D Laplace equation-based model match each other well, which indicates that the horizontal flow in drainage layer is dominated. Also, it validates the 1D Boussinesq equation-based model which can be applied to predict the water-table elevation in drainage layer. Further, the analysis was conducted to examine the effect of infiltration rate, hydraulic conductivity and slope of drainage layer on the water-table elevation. The results show that water-table falls down as the ratio of Is to K decreases and the slope increases. If the aquifer becomes confined by the top of drainage layer due to the larger ratio of Is to K or smaller slope, the solution presented in this work can also be applied to approximate the water-table elevation in unconfined sub-section as well as hydraulic head in the confined sub-section.展开更多
基金Supported by Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)+2 种基金Basic Research Plan of Shanxi Province(202203021211129)Shanxi Province Natural Science Research(202203021212249)Special/Youth Foundation of Taiyuan University of Technology(2022QN101)。
文摘In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,the original Boussinesq system is transformed into an equivalent one.Then we discretize it using the second-order backward di erentiation formula(BDF2)and Crank-Nicolson(CN)to obtain two second-order time-advanced schemes.In both numerical schemes,a pressure-correction method is employed to decouple the velocity and pressure.These two schemes possess the desired property that they can be fully decoupled with satisfying unconditional stability.We rigorously prove both the unconditional stability and unique solvability of the discrete schemes.Furthermore,we provide detailed implementations of the decoupling procedures.Finally,various 2D numerical simulations are performed to verify the accuracy and energy stability of the proposed schemes.
基金Project(511114) supported by the Natural Science Foundation of Hainan Province, ChinaProject(2009YBFZ05) supported by Postgraduate Award of Central South University, China+1 种基金Project(200731) supported by Traffic Technology Fund of Hunan Province, ChinaProject(2008BAG10B01) supported by the National Key Technology R&D Program, China
文摘In order to study the water flow in the drainage layer of highway under steady-state condition, one-dimensional （1D） Boussinesq equation-based model with Dupuit-Forchheimer assumption was established and the semi-analytical solutions to predict the water-table height were presented. In order to validate the model, a two-dimensional （2D） saturated flow model based on the Laplace equation was applied for the purpose of the model comparison. The water-table elevations predicted by ID Boussinesq equation-based model and 2D Laplace equation-based model match each other well, which indicates that the horizontal flow in drainage layer is dominated. Also, it validates the 1D Boussinesq equation-based model which can be applied to predict the water-table elevation in drainage layer. Further, the analysis was conducted to examine the effect of infiltration rate, hydraulic conductivity and slope of drainage layer on the water-table elevation. The results show that water-table falls down as the ratio of Is to K decreases and the slope increases. If the aquifer becomes confined by the top of drainage layer due to the larger ratio of Is to K or smaller slope, the solution presented in this work can also be applied to approximate the water-table elevation in unconfined sub-section as well as hydraulic head in the confined sub-section.
