We present a compact, high-order Richards’ equation solver using a local discontinuous Galerkin finite element method in space and a dual-time stepping method in time. Dual-time stepping methods convert a transient problem to a steady state problem, enabling direct evaluation of residual terms and resolve implicit equations in a step-wise manner keeping the method compact and amenable to parallel computing. Verification of our solver against an analytical solution shows high-order error convergence and demonstrates the solvers ability to maintain high accuracy using low spatial resolution; the method is robust and accurately resolves numerical solutions with time steps that are much larger than what is normally required for lower-order implicit schemes. Resilience of our solver (in terms of nonlinear convergence) is demonstrated in ponded infiltration into homogeneous and layered soils, for which HYDRUS-1D solutions are used as qualitative references to gauge performance of two slope limiting schemes.
Landslides - The multiple-barrier mitigation strategy in the debris-flow source area is an effective approach to inhibit debris-flow entrainment and scale amplification along the flow path.... 相似文献
The northern Wuyi area, which is located in the northern Wuyi metallogenic belt, has superior mineralization conditions. The Pingxiang–Guangfeng–Jiangshan–Shaoxing fault (PSF) extends across the whole region regardless of whether or how the PSF relates to the near-surface mineralization. We carried out an MT survey in the region and obtained a reliable 2D model of the crustal electrical structure to a depth of 30 km. In the resistivity model, we inferred that a continuous high conductivity belt that ranges from the shallow to deep crust is a part of the PSF. Then, we estimated the fluid content and pressure gradient to identify the deep sources of fluid as well as its pattern of motion pattern. Finally, we proposed a model for the deep metallogenic migration processes that combines geological data, fluid content data, pressure gradient data, and the subsurface resistivity model. The model analysis showed that the Jiangnan orogenic belt and the Cathaysia block formed the PSF during the process of com. The deep fluid migrated upward through the PSF to the shallow crust. Therefore, we believe that the PSF is an ore-forming fluid migration channel and that it laid the material basis for large-scale mineralization in the shallow crust. 相似文献