An efficient discontinuous Galerkin method for advective transport in porous media

We consider a discontinuous Galerkin scheme for computing transport in heterogeneous media. An efficient solution of the resulting linear system of equations is possible by taking advantage of a priori knowledge of the direction of flow. By arranging the elements in a suitable sequence, one does not need to assemble the full system and may compute the solution in an element-by-element fashion. We demonstrate this procedure on boundary-value problems for tracer transport and time-of-flight.     

企业研究的时间地理学框架 --兼论泰勒模式的时间地理学解释

企业空间扩张是公司地理学研究的重要内容。文章引入关注时空坐标系下对个人或群体的连续活动以及所面临的制约条件进行独到分析的时间地理学概念和理论体系，通过“时空棱柱”、“制约”等核心概念扩张性的应用，对企业空间扩张的代表性理论－泰勒模式做出新的诠释，并以此为基础尝试提出企业研究的时间地理学框架，展望了时间地理学在有关企业家研究、创新扩散研究、企业时间管理、企业行为的制约体系等方面的应用前景。     

Griffith Taylor's ‘palimpsest theory’ of the dry valleys of Victoria land,antarctica

Griffith Taylor was the first geomorphologist to work in the Dry Valleys of southern Victoria Land, Antarctica. Following his field work in February 1911, he proposed a multistage model in which earlier cirque erosion was later swamped by expanding outlet glaciers. Subsequently these glaciers retreated leaving the present form of the valleys. The topography retained the imprint of each episode, hence his name ‘palimpsest theory’. I summarise later research and compare Taylor's theory with current views.     

用无网格伽辽金法求解稳定地下水流问题

详细介绍了无网格伽辽金法（EFGM）基本原理,并将其应用于非均质多孔介质中的稳定地下水流问题,用具体算例将无网格伽辽金法计算结果与传统有限元法（LFEM）计算结果作比较,计算表明无网格伽辽金法具有较高的精度。     

Higher-order approximation techniques for estimating stochastic parameter of a sediment transport model

Higher-order approximation techniques for estimating stochastic parameter of the non-homogeneous Poisson (NHP) model are presented. The NHP model is characterized by a two-parameter cumulative probability distribution function (CDF) of sediment displacement. Those two parameters are the temporal and spatial intensity functions, physically representing the inverse of the average rest period and step length of sediment particles, respectively. Difficulty of estimating the parameters has, however, restricted the applications of the NHP model. The approximation techniques are proposed to address such problem. The basic idea of the method is to approximate a model involving stochastic parameters by Taylor series expansion. The expansion preserves certain higher-order terms of interest. Using the experimental (laboratory or field) data, one can determine the model parameters through a system of equations that are simplified by the approximation technique. The parameters so determined are used to predict the cumulative distribution of sediment displacement. The second-order approximation leads to a significant reduction of the CDF error (of the order of 47%) compared to the first-order approximation. Error analysis is performed to evaluate the accuracy of the first- and second-order approximations with respect to the experimental data. The higher-order approximations provide better estimations of the sediment transport and deposition that are critical factors for such environment as spawning gravel-bed.     

On the evaluation of time-domain Green function

An analytical method has been developed to evaluate the wave part of the time-domain Green function and its derivatives. Based on Taylor series expansion, the Green function is obtained by solving a fourth-order ordinary differential equation. The method accelerates the convergence of the summation of an infinite series in the numerical computation. The accuracy of this method was demonstrated by its comparison with other method and its application to solve the radiation problem of a floating hemisphere using a panel-free method. The computed hydrodynamic coefficients agree well with the analytical solutions.     

Application of a coupled discontinuous–continuous Galerkin finite element shallow water model to coastal ocean dynamics

A coupled discontinuous–continuous Galerkin (DG–CG) shallow water model is compared to a continuous Galerkin generalized wave-continuity equation (GWCE) based model for the coastal ocean, whereby local mass imbalance typical of GWCE-based solutions is eliminated using the coupled DG–CG approach. Two mass imbalance indicators for the GWCE-based model are presented and analyzed. The indicators motivate discussion on the suitability of using a GWCE-based model versus the locally conservative coupled DG–CG model. Both realistic and idealized test problems for tide, wind, and wave-driven circulation form the basis of the study. For the problems studied, coupled DG–CG solutions retain the robustness of well-documented solutions from GWCE-based models and also capture the dynamics driven by small-scale, highly advective processes which are problematic for GWCE-based models. Issues associated with the coupled DG–CG model are explored, including increased cost due to increased degrees of freedom, the necessary application of slope limiters, as well as the actual coupling process.     

Paleomagnetism of ordovician lamprophyres from Taylor Valley,Victoria Land,Antarctica

Paleomagnetic study of specimens from four lamprophyric dykes on the Kukri Hills, Taylor Valley (77.64°S, 163.35°E) has yielded a primary mean direction of magnetization ofD=222.6 andI=+0.6 with ₉₅=10.9° after AF cleaning. The magnetization of five other dykes and of the amphibolitic basement was either unstable or not fully reliable.The corresponding pole position lies at 9.3°S and 26.7°E and confirms the previous results from Lower Ordovician rocks from distant areas of East Antarctica.A Lower Ordovician mean pole position recalculated from valid data lies at 17°S, 21°E.     

Some results comparing Monte Carlo simulation and first order Taylor series approximation for steady groundwater flow

The expected head and standard deviation of the head from the first order Taylor series approximation is compared to Monte Carlo simulation, for steady flow in a confined aquifer with transmissivity as a random variable. Emphasis is on the effect of changes in the covariance structure of the transmissivity, and pumping rates, on the errors in the first order Taylor series approximation. The accuracy of the first order Taylor series approximation is found to be particularly sensitive to pumping rates. With significant pumping the approximation is found to under estimate both the expected drawdown and head variance, and the error increases as the pumping rate increases. This can lead to large errors in probability constraints based on moments from the first order Taylor series approximation.     

Stability of buoyancy opposed mixed convection in a vertical channel and its dependence on permeability

Non-Darcy mixed convective flow of water due to external pressure gradient and buoyancy opposed forces are considered in a vertical channel filled with porous medium, which can be either isotropic or anisotropic. The linear theory of stability analysis has been used to numerically investigate the dependence of the transition behavior of the fully developed basic flow on the permeability of the medium. Numerical experiments indicate that mainly two main instability modes appear: Rayleigh–Taylor (R–T) and buoyant instability. For Darcy numbers (Da) ?10⁻⁹, R–T instability dominates within the entire Reynolds number (Re) range considered here. It was also found that for the same Re, the fully developed base flow is highly unstable (stable) for porous media with high (low) permeability. Further, it was seen that the disturbance isotherm cells migrate from the channel walls toward the centerline when permeability is reduced. Reducing the permeability by one order of magnitude (corresponding to a decrease of Darcy number from 10⁻⁶ to 10⁻⁷) increases base flow stability approximately 20-fold. For higher Reynolds numbers, buoyant, mixed and shear instability of the basic flow were found when Da was increased from 10⁻⁷ to 10⁻³. However, for cases in which permeability and porosity behaved as suggested by Carman–Kozeny relation (CKR), buoyant stability was the only mode of instability. Critical values of the Rayleigh (Ra) and Darcy (Da) numbers in the R–T mode of instability were related to each other by the hyperbolic function RaDa = −2.465.