Search

Abstract: A groundwater flow model is typically used to provide the flow field for conducting groundwater solute transport simulations. The advection term of the mass conserved formulation for groundwater transport assumes that the flow field is perfectly balanced and that all water flowing into a numerical grid cell is exactly balanced by outflows after accounting for sources/sinks or internal storage. However, in many complicated regional or site‐scale models, there may be localized flow balance errors that may be difficult to eliminate through tighter flow convergence tolerances due to simulation time constraints or numerical limits on convergence tolerances. Thus, if water is erroneously gained or lost within a grid cell during the flow computation, the solutes within it will also be numerically affected in the associated transport simulation. Transport solutions neglect this error in groundwater flow as the transport equations that are solved assume no error in flow. This flow imbalance error can however have consequences on the transport solution ranging from unnoticeable errors in the resulting concentrations to spurious oscillations that can grow in time and hinder further solution. An approach has been suggested here, to explicitly handle these flow imbalances during mass conserved advective transport computations and report them in the corresponding transport mass balance output, as corrections that are needed to handle errors originating in the flow solution. Example problems are provided to explain the concepts and demonstrate the impacts. [ABSTRACT FROM AUTHOR]

Abstract: As an alternative to grid refinement, the concept of a ghost node, which was developed for nested grid applications, has been extended towards improving sub-grid scale accuracy of flow to conduits, wells, rivers or other boundary features that interact with a finite-difference groundwater flow model. The formulation is presented for correcting the regular finite-difference groundwater flow equations for confined and unconfined cases, with or without Newton Raphson linearization of the nonlinearities, to include the Ghost Node Correction (GNC) for location displacement. The correction may be applied on the right-hand side vector for a symmetric finite-difference Picard implementation, or on the left-hand side matrix for an implicit but asymmetric implementation. The finite-difference matrix connectivity structure may be maintained for an implicit implementation by only selecting contributing nodes that are a part of the finite-difference connectivity. Proof of concept example problems are provided to demonstrate the improved accuracy that may be achieved through sub-grid scale corrections using the GNC schemes. [Copyright &y& Elsevier]

A diverse group of problems requires quantification of the entire hydrologic cycle by the integrated simulation of water flow in the surface and subsurface regimes. In a transient integrated simulation of the water cycle, the time step size is a key factor in controlling the solution accuracy and the simulation efficiency for a given spatial discretization. In general, if the time step size is sufficiently small, the resulting solution becomes more accurate but with higher computational cost. Thus, to maintain an acceptable level of solution accuracy in the entire simulation domain, the time step size is restricted by the relatively rapid responses in the surface flow regime. As the relatively rapid responses are typically limited to a small portion of the surface domain compared with the groundwater system, a large portion of the domain tends to be temporally overdiscretized. The implicit subtime stepping approach described here can apply smaller subtime steps only to the subdomain where the accuracy requirements are needed. In this work, generalized formulations for implicit subtime stepping in the numerical solution of the nonlinear coupled surface-subsurface equations were derived and implemented into the integrated model HydroGeoSphere. Application to several problems showed that implicit subtime stepping can significantly improve the simulation efficiency with minimal loss in accuracy. The methodology was successfully applied to enhance the computational efficiency of an integrated flow simulation in the San Joaquin Valley, California, where the characteristic response time near surface drainage streams is orders of magnitude shorter than in the groundwater regime. [ABSTRACT FROM AUTHOR]