Efficient flow and transport simulations in reconstructed 3D pore geometries

Upscaling pore-scale processes into macroscopic quantities such as hydrodynamic dispersion is still not a straightforward matter for porous media with complex pore space geometries. Recently it has become possible to obtain very realistic 3D geometries for the pore system of real rocks using either numerical reconstruction or micro-CT measurements. In this work, we present a finite element–finite volume simulation method for modeling single-phase fluid flow and solute transport in experimentally obtained 3D pore geometries. Algebraic multigrid techniques and parallelization allow us to solve the Stokes and advection–diffusion equations on large meshes with several millions of elements. We apply this method in a proof-of-concept study of a digitized Fontainebleau sandstone sample. We use the calculated velocity to simulate pore-scale solute transport and diffusion. From this, we are able to calculate the a priori emergent macroscopic hydrodynamic dispersion coefficient of the porous medium for a given molecular diffusion D_m of the solute species. By performing this calculation at a range of flow rates, we can correctly predict all of the observed flow regimes from diffusion dominated to convection dominated.     

An Algebraic Method to Compute the Critical Points of the Distance Function Between Two Keplerian Orbits

We describe an efficient algorithm to compute all the critical points of the distance function between two Keplerian orbits (either bounded or unbounded) with a common focus. The critical values of this function are important for different purposes, for example to evaluate the risk of collisions of asteroids or comets with the Solar system planets. Our algorithm is based on the algebraic elimination theory: through the computation of the resultant of two bivariate polynomials, we find a 16th degree univariate polynomial whose real roots give us one component of the critical points. We discuss also some degenerate cases and show several examples, involving the orbits of the known asteroids and comets.      

天气雷达测定区域降水量方法的改进与比较

讨论变分校准法用于雷达－雨量计系统联合探测降水。由雷达反射率因子Ｚ和地面降水强度Ｉ实时地获得最优Ｚ－Ｉ关系，在求解欧拉方程时采用多重网格法，不仅可提高计算结果的精度，还可大大提高计算速度。     

Review of the algebraic linear methods and parallel implementation in numerical simulation of groundwater flow

The desire to increase spatial and temporal resolution in modeling groundwater system has led to the requirement for intensive computational ability and large memory space. In the course of satisfying such requirement, parallel computing has played a core role over the past several decades. This paper reviews the parallel algebraic linear solution methods and the parallel implementation technologies for groundwater simulation. This work is carried out to provide guidance to enable modelers of groundwater systems to make sensible choices when developing solution methods based upon the current state of knowledge in parallel computing.     

System of Multigrid Nonlinear Least-squares Four-dimensional Variational Data Assimilation for Numerical Weather Prediction (SNAP): System Formulation and Preliminary Evaluation

A new forecasting system—the System of Multigrid Nonlinear Least-squares Four-dimensional Variational(NLS-4DVar) Data Assimilation for Numerical Weather Prediction(SNAP)—was established by building upon the multigrid NLS-4DVar data assimilation scheme, the operational Gridpoint Statistical Interpolation(GSI)-based data-processing and observation operators, and the widely used Weather Research and Forecasting numerical model. Drawing upon lessons learned from the superiority of the operational GS...