Development of a cross-section based streamflow routing package for MODFLOW

It is challenging to simulate stream-aquifer interactions for the wide channel streams with the existing stream routing packages of MODFLOW. To overcome this limitation, a Cross-Section streamflow Routing (CSR) package is developed to simulate the streamflow and the interaction between streams and aquifers for the stream with a width larger than the MODFLOW grid size. In the CSR package, streams are divided into stream segments which are formed by two consecutive cross-sections. A cross-section is described by a number of streambed points that determine the geometry and hydraulic properties of the streambed. The stream water depth and streamflow at the cross-sections are related by the Single Channel method, the Divided Channel method, a data table or a power function. A rapid algorithm is used to compute the submerged area of the MODFLOW grid. The streambed conductance of a grid cell is computed based on its submerged area, streambed hydraulic conductivity and thickness. Stream-aquifer seepage is subsequently estimated as the product of the streambed conductance and difference between the stream stage and groundwater hydraulic head. Stream-aquifer seepage is treated as lateral flow in the streamflow routing computation with the Muskingum-Cunge method or mass conservation method. A hypothetical problem is established to test the capabilities of the CSR package with steady- and transient-state models. The results compare favorably with the SFR2 package and the HEC-RAS model. However, significant difference in flood wave attenuation is observed between the CSR package and the SFR2 package. It proves that the CSR package is capable of simulating the variation of stream-aquifer interactions in both space and time efficiently. The CSR package represents a certain improvement over previous MODFLOW streamflow packages by providing the efficient cross-section based computation and the unique capability of simulating streambed heterogeneity in longitudinal and transverse directions.     

Building a geodatabase for mapping hydrogeological features and 3D modeling of groundwater systems: Application to the Saguenay-Lac-St.-Jean region, Canada

Understanding and managing groundwater resources require the integration of a large amount of high-quality data from a variety of sources. Due to the limitations in accessing information related to groundwater and subsurface conditions, the gathering of available existing information is of crucial importance when conducting a successful hydrogeological study. Here, we present an approach for the development of an exhaustive and comprehensive groundwater database through (1) the gathering of relevant sources of information relating to groundwater, and (2) the application of a quality control process in order to screen the data for accuracy and quality. This hydrogeological database is then implemented within a GIS (geographic information system) framework coupled to a Relational Database Management System (RDBMS) as a personal geodatabase (ESRI format) GIS technology. Once established, the spatial database allows a user to request the relevant data required for a specific hydrogeological study. In addition, stratigraphic data stored within the spatial database may be utilized for constructing 3D subsurface hydrostructural models. In order to achieve this objective, the software Arc Hydro Groundwater combined with the ArcGIS spatial database is shown to be appropriate for the 3D structural representation of aquifers (groundwater reservoirs). The innovative contribution of this approach in building 3D hydrostructural subsurface models from a spatial database resides in simplifying the required step-by-step processes by considering a unified compatible combination of “RDBMS-ArcGIS-Arc Hydro Groundwater” technologies. The proposed methodology is illustrated using data from an ongoing project aimed at developing an inventory of the groundwater resources of the Saguenay-Lac-Saint-Jean region, Quebec (Canada).     

Unconstrained periodic boundary conditions for solid state elasticity

We introduce a method to implement dynamics on an elastic lattice without imposing constraints via boundary or loading conditions. Using this method we are able to examine fracture processes in two-dimensional systems previously inaccessible for reliable computer simulations. We show the validity of the method by benchmarking and report a few preliminary results.     

Spectral and pseudo-spectral methods for parabolic problems with non periodic boundary conditions

The advection-diffusion equation is approximated by Chebyshev and Legendre spectral and pseudo-spectral methods. Stability results in the energy norm and error estimates in terms of the discretization parameter and of the regularity of the solution in weighted Sobolev norms are presented.     

叉指式加速度计的一种仿真法分析

微电子机械系统的快速发展对其建模仿真也提出了较高要求。首先介绍了一种微电子机械系统原理级描述仿真方法,它基于成熟的微机械元件库,提供微器件的原理级描述,使得仿真更简单、更高效,但目前其适用范围有限。然后介绍了用这种方法分析叉指式微电容加速度计的结果。     

A new fourth-order difference scheme for solving an N-carrier system with Neumann boundary conditions

In this paper, a numerical method is developed to solve an N-carrier system with Neumann boundary conditions. First, we apply the compact finite difference scheme of fourth order for discretizing spatial derivatives at the interior points. Then, we develop a new combined compact finite difference scheme for the boundary, which also has fourth-order accuracy. Lastly, by using a Padé approximation method for the resulting linear system of ordinary differential equations, a new compact finite difference scheme is obtained. The present scheme has second-order accuracy in time direction and fourth-order accuracy in space direction. It is shown that the scheme is unconditionally stable. The present scheme is tested by two numerical examples, which show that the convergence rate with respect to the spatial variable from the new scheme is higher and the solution is much more accurate when compared with those obtained by using other previous methods.     

A model of seepage field in the tailings dam considering the chemical clogging process

The radial collector well, an important water drainage construction, has been widely applied to the tailings dam. Chemical clogging frequently occurs around the vertical shaft in radial collector well due to enough dissolved oxygen and some heavy metals in groundwater flow of tailings dam. Considering the contribution of water discharge from both vertical shaft and horizontal screen laterals and chemical clogging occurring around vertical shaft well, a new model was developed on the basis of Multi-Node Well (MNW2) package of MODFLOW. Moreover, two cases were calculated by the newly developed model. The results indicate that the model considering chemical clogging occurring around the vertical shaft well is reasonable. Owing to the decrease in hydraulic conductivity caused by chemical clogging, the groundwater level in dam body increases constantly and water discharge of radial collector well declines by 10-15%. For ordinary vertical well, it decreases by 30%. Therefore, chemical clogging occurring around radial collector well can arouse increases of groundwater level, and influence dambody safety.     

A comprehensive modeling environment for the simulation of groundwater flow and transport

A comprehensive graphical modeling environment has been developed to address the needs of the computer simulation of groundwater flow and transport. The Department of Defense Groundwater Modeling Systems (GMS), developed at the Engineering Computer Graphics Laboratory at Brigham Young University, is part of a multi-year project funded through the Department of Defense, Department of Energy and Environmental Protection Agency. GMS is a graphically based software tool providing facility through all aspects of the groundwater flow and transport modeling process. Facilities include geometric modeling of hydrostratigraphy, two- and three-dimensional mesh generation, graphically based model input for specific flow and transport codes, interpolation and geostatistics as well as complete three-dimensional scientific visualization.     

An R package for simulating metapopulation dynamics and range expansion under environmental change

The metapopulation paradigm is central in ecology and conservation biology to understand the dynamics of spatially-structured populations in fragmented landscapes. Metapopulations are often studied using simulation modelling, and there is an increasing demand of user-friendly software tools to simulate metapopulation responses to environmental change. Here we describe the MetaLandSim R package, which integrates ideas from metapopulation and graph theories to simulate the dynamics of real and virtual metapopulations. The package offers tools to (i) estimate metapopulation parameters from empirical data, (ii) to predict variation in patch occupancy over time in static and dynamic landscapes, either real or virtual, and (iii) to quantify the patterns and speed of metapopulation expansion into empty landscapes. MetaLandSim thus provides detailed information on metapopulation processes, which can be easily combined with land use and climate change scenarios to predict metapopulation dynamics and range expansion for a variety of taxa and ecological systems.     

A model integration framework for linking SWAT and MODFLOW

Assessment of long-term anthropogenic impacts on agro-ecosystems requires comprehensive modelling capabilities to simulate water interactions between the surface and groundwater domains. To address this need, a modelling framework, called “SWATmf”, was developed to link and integrate the Soil Water Assessment Tool (SWAT), a widely used surface watershed model with the MODFLOW, a groundwater model. The SWATmf is designed to serve as a project manager, builder, and model performance evaluator, and to facilitate dynamic interactions between surface and groundwater domains at the watershed scale, thus providing a platform for simulating surface and groundwater interactions. Using datasets from the Fort Cobb Reservoir experimental watershed (located in Oklahoma, USA), the SWATmf to facilitate linkage and dynamic simulation of SWAT and MODFLOW models. Simulated streamflow and groundwater levels generally agreed with observations trends showing that the SWATmf can be used for simulating surface and groundwater interactions.     

Creases and boundary conditions for subdivision curves

Our goal is to find subdivision rules at creases in arbitrary degree subdivision for piece-wise polynomial curves, but without introducing new control points e.g. by knot insertion. Crease rules are well understood for low degree (cubic and lower) curves. We compare three main approaches: knot insertion, ghost points, and modifying subdivision rules. While knot insertion and ghost points work for arbitrary degrees for B-splines, these methods introduce unnecessary (ghost) control points.The situation is not so simple in modifying subdivision rules. Based on subdivision and subspace selection matrices, a novel approach to finding boundary and sharp subdivision rules that generalises to any degree is presented. Our approach leads to new higher-degree polynomial subdivision schemes with crease control without introducing new control points.     

A new algorithm for scheduling periodic, real-time tasks

We consider the problem of preemptively scheduling a set of periodic, real-time tasks on a multiprocessor computer system. We give a new scheduling algorithm, the so-called Slack-Time Algorithm, and show that it is more effective than the known Deadline Algorithm. We also give an (exponential-time) algorithm to decide if a task system is schedulable by the Slack-Time or the Deadline Algorithm. The same algorithm can also be used to decide if a task system is schedulable by any given fixed-priority scheduling algorithm. This resolves an open question posed by Leung and Whitehead. Finally, it is shown that the problem of deciding if a task system is schedulable by the Slack-Time, the Deadline, or any given fixed-priority scheduling algorithm is co-NP-hard for each fixedm.     

Implementation of periodic boundary conditions for loading of mechanical metamaterials and other complex geometric microstructures using finite element analysis

The implementation of periodic boundary conditions (PBCs) is one of the most important and difficult steps in the computational analysis of structures and materials. This is especially true in cases such as mechanical metamaterials which typically possess intricate geometries and designs which makes finding and implementing the correct PBCs a difficult challenge. In this work, we analyze one of the most common PBCs implementation technique, as well as implement and validate an alternative generic method which is suitable to simulate any possible 2D microstructural geometry with a quadrilateral unit cell regardless of symmetry and mode of deformation. A detailed schematic of how both these methods can be employed to study 3D systems is also presented.

     

Well-posedness in the generalized sense for boundary layer suppressing boundary conditions

The time-dependent Navier-Stokes equation for incompressible fluid flow together with new boundary layer suppressing boundary conditions for open boundaries is investigated. In these new boundary conditions one typically prescribes a high-order derivative of some of the dependent variables. We prove that these boundary conditions give rise to a problem that is well posed in the generalized sense. This means that there exists a unique smooth solution of the linearized problem and that this solution can be estimated by data.     

A parallel two-level method for simulating blood flows in branching arteries with the resistive boundary condition

Computer modeling of blood flows in the arteries is an important and very challenging problem. In order to understand, computationally, the sophisticated hemodynamics in the arteries, it is essential to couple the fluid flow and the elastic wall structure effectively and specify physiologically realistic boundary conditions. The computation is expensive and the parallel scalability of the solution algorithm is a key issue of the simulation. In this paper, we introduce and study a parallel two-level Newton–Krylov–Schwarz method for simulating blood flows in compliant branching arteries by using a fully coupled system of linear elasticity equation and incompressible Navier–Stokes equations with the resistive boundary condition. We first focus on the accuracy of the resistive boundary condition by comparing it with the standard pressure type boundary condition. We then show the parallel scalability results of the two-level approach obtained on a supercomputer with a large number of processors and on problems with millions of unknowns.     

A new adaptive control for periodic tracking/disturbance rejection

This paper proposes a new adaptive feedforward cancellation (AFC) control that achieves periodic tracking and/or periodic disturbance rejection. The new control design is a direct scheme in the sense that it adaptively updates the desired control without estimating the unknown disturbance. The proposed new control has several advantages. First, its adaptation gain can be arbitrarily chosen without upsetting the system stability. Second, it can be applied to not only minimum‐phase systems, but also non‐minimum phase systems. Finally, it is shown that the proposed AFC control is independent of where the disturbance enters the system. Copyright © 2011 John Wiley and Sons Asia Pte Ltd and Chinese Automatic Control Society     

A new fair scheduling algorithm for periodic tasks on multiprocessors

We present a new scheduling algorithm, called PL that is work-conserving and in terms of schedulability, optimal on multiprocessors for a synchronous periodic task set. The PL algorithm is a laxity based algorithm and ensures execution of a task with approximate proportional fairness at each task's period. Existing optimal algorithms on multiprocessors may cause excessive scheduling decisions and preemptions or may not be applied in a discrete environment. The proposed algorithm can be applied in a discrete environment and reduce the number of scheduling decisions and preemptions compared with a Pfair algorithm.     

A new model for Brownian force and the application to simulating nanofluid flow

A new expression is proposed to simulate Brownian force based on the experimental measurement results of Brownian motion, which follows white Gaussian noise process. As the time t → 0 and the particle density is equal to the fluid density, the new expression approaches the classical formula of the model used by many researchers. The modified model is validated by theoretical and experimental data. On the other hand, as it origins from the unbalanced force exerted by surrounding fluid molecules, the drag analogy force model is constructed describing the Brownian force, which depends on size-related statistical velocity. Thus, a different expression for the Langevin equation is presented. The present model is applied in simulating flow and heat transfer in a channel utilizing alumina–water nanofluid. Navier–Stokes equations with modified source terms for the continuous flow have been discretized using finite element method. The velocities and temperatures of nanoparticles are determined in the Lagrangian reference frame. The simulation results show that the distribution of nanoparticles inside the channel is obviously unsteady and nonuniform. The fluid velocity and temperature profiles show significant fluctuation feature at low Reynolds numbers (Re). The impact of Brownian motion on the fluid flow is analyzed quantitatively. We have found that for Re < 0.06, the affected intensity increases rapidly.