Viscoelastic artificial boundaries are widely adopted in numerical simulations of wave propagation problems.When explicit time-domain integration algorithms are used,the stability condition of the boundary domain is s...Viscoelastic artificial boundaries are widely adopted in numerical simulations of wave propagation problems.When explicit time-domain integration algorithms are used,the stability condition of the boundary domain is stricter than that of the internal region due to the influence of the damping and stiffness of an viscoelastic artificial boundary.The lack of a clear and practical stability criterion for this problem,however,affects the reasonable selection of an integral time step when using viscoelastic artificial boundaries.In this study,we investigate the stability conditions of explicit integration algorithms when using three-dimensional(3D)viscoelastic artificial boundaries through an analysis method based on a local subsystem.Several boundary subsystems that can represent localized characteristics of a complete numerical model are established,and their analytical stability conditions are derived from and further compared to one another.The stability of the complete model is controlled by the corner regions,and thus,the global stability criterion for the numerical model with viscoelastic artificial boundaries is obtained.Next,by analyzing the impact of different factors on stability conditions,we recommend a stability coefficient for practically estimating the maximum stable integral time step in the dynamic analysis when using 3D viscoelastic artificial boundaries.展开更多
Groundwater flows play a key role in the recharge of aquifers, the transport of solutes through subsurface systems or the control of surface runoff. Predicting these processes requires the use of groundwater models wi...Groundwater flows play a key role in the recharge of aquifers, the transport of solutes through subsurface systems or the control of surface runoff. Predicting these processes requires the use of groundwater models with their applicability directly linked to their accuracy and computational efficiency. In this paper, we present a new method to model water dynamics in variably- saturated porous media. Our model is based on a fully-explicit discontinuous-Galerkin formulation of the 3D Richards equation, which shows a perfect scaling on parallel architectures. We make use of an adapted jump penalty term for the discontinuous-Galerkin scheme and of a slope limiter algorithm to produce oscillation-free exactly conservative solutions. We show that such an approach is particularly well suited to infiltration fronts. The model results are in good agreement with the reference model Hydrus-lD and seem promising for large scale applications involving a coarse representation of saturated soil.展开更多
基金National Natural Science Foundation of China under Grant Nos.52108458 and U1839201China National Postdoctoral Program of Innovative Talents under Grant No.BX20200192+1 种基金Shuimu Tsinghua Scholar Program under Grant No.2020SM005National Key Research and Development Program of China under Grant No.2018YFC1504305。
文摘Viscoelastic artificial boundaries are widely adopted in numerical simulations of wave propagation problems.When explicit time-domain integration algorithms are used,the stability condition of the boundary domain is stricter than that of the internal region due to the influence of the damping and stiffness of an viscoelastic artificial boundary.The lack of a clear and practical stability criterion for this problem,however,affects the reasonable selection of an integral time step when using viscoelastic artificial boundaries.In this study,we investigate the stability conditions of explicit integration algorithms when using three-dimensional(3D)viscoelastic artificial boundaries through an analysis method based on a local subsystem.Several boundary subsystems that can represent localized characteristics of a complete numerical model are established,and their analytical stability conditions are derived from and further compared to one another.The stability of the complete model is controlled by the corner regions,and thus,the global stability criterion for the numerical model with viscoelastic artificial boundaries is obtained.Next,by analyzing the impact of different factors on stability conditions,we recommend a stability coefficient for practically estimating the maximum stable integral time step in the dynamic analysis when using 3D viscoelastic artificial boundaries.
基金funded by the Fond de la Recherche Scientifique de Belgique (FRSFNRS)
文摘Groundwater flows play a key role in the recharge of aquifers, the transport of solutes through subsurface systems or the control of surface runoff. Predicting these processes requires the use of groundwater models with their applicability directly linked to their accuracy and computational efficiency. In this paper, we present a new method to model water dynamics in variably- saturated porous media. Our model is based on a fully-explicit discontinuous-Galerkin formulation of the 3D Richards equation, which shows a perfect scaling on parallel architectures. We make use of an adapted jump penalty term for the discontinuous-Galerkin scheme and of a slope limiter algorithm to produce oscillation-free exactly conservative solutions. We show that such an approach is particularly well suited to infiltration fronts. The model results are in good agreement with the reference model Hydrus-lD and seem promising for large scale applications involving a coarse representation of saturated soil.