Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

IMPLICIT NONLINEAR NORMAL MODE INITIALIZATION:A MULTIGRID APPROACH

The implicit nonlinear normal mode initialization (INMI) is applied to a tropical limited area shallow water modelin spherical coordinates.The boundary condition for the INMI scheme is based on the boundary formulation of themodel.The INMI scheme is found to be very efficient in suppressing spurious gravity wave oscillation and providing awell balanced initial data set for the model.The INMI scheme involves solving a number of elliptic type equations withvarying complexity.and hence an efficient numerical technique is required for solving such equations.In order to makeINMI computationally more attractive,we are employing the multigrid method for solving all the elliptic type equationsin the INMI scheme.The numerical procedures for the development of such multigrid solvers are briefly described.     

Efficient block preconditioners for the coupled equations of pressure and deformation in highly discontinuous media

Large‐scale simulations of flow in deformable porous media require efficient iterative methods for solving the involved systems of linear algebraic equations. Construction of efficient iterative methods is particularly challenging in problems with large jumps in material properties, which is often the case in geological applications, such as basin evolution at regional scales. The success of iterative methods for this type of problems depends strongly on finding effective preconditioners. This paper investigates how the block‐structured matrix system arising from single‐phase flow in elastic porous media should be preconditioned, in particular for highly discontinuous permeability and significant jumps in elastic properties. The most promising preconditioner combines algebraic multigrid with a Schur complement‐based exact block decomposition. The paper compares numerous block preconditioners with the aim of providing guidelines on how to formulate efficient preconditioners. Copyright © 2010 John Wiley & Sons, Ltd.     

多重网格法在MM5自适应模式中的应用

将多重网格法应用于MM5中来加速自适应网格的生成,并用这种方法来模拟2002年6月17-18日发生在江南和华南的暴雨过程。试验表明：在收敛精度要求高的情况下自适应模式模拟结果的精度高于均匀网格模式,且多重网格法有助于提高运算效率。     

直流电阻率法2.5维正演的外推瀑布式多重网格法

引入外推瀑布式多重网格法(EXCMG)求解2.5维直流电阻率有限元计算形成的大型稀疏线性方程组,结合基于地址矩阵的压缩存贮方式以及最优化离散波数,使得2.5维电阻率正演程序的计算速度大大提高而内存需求大大减小.研究结果表明:EXCMG法的收敛速度与网格尺寸无关,计算速度明显优于不完全Cholesky共轭梯度(ICCG)...     

基于粗细网格的有限元并行分析方法

并行计算己成为求解大规模岩土工程问题的一种强大趋势。探讨了粗细网格与预处理共轭梯度法结合的并行有限元算法。从多重网格刚度矩阵推得有效的预处理子。该算法在工作站机群上实现。用地基处理时土体强夯的数值模拟分析进行了数值测试,对其并行性能进行了详细分析。计算结果表明：该算法具有良好的并行加速比和效率,是一种有效的并行算法。     

A PRELIMINARY STUDY ON THE 3DVAR ASSIMILATION OF THE AMSU-A DATA IN SPACE-TIME MULTISCALE ANALYSIS SYSTEM

Assimilating satellite radiances into Numerical Weather Prediction (NWP) models has become an important approach to increase the accuracy of numerical weather forecasting. In this study, the assimilation technique scheme was employed in NOAA’s STMAS (Space-Time Multiscale Analysis System) to assimilate AMSU-A radiances data. Channel selection sensitivity experiments were conducted on assimilated satellite data in the first place. Then, real case analysis of AMSU-A data assimilation was performed. The analysis results showed that, following assimilating of AMSU-A channels 5–11 in STMAS, the objective function quickly converged, and the channel vertical response was consistent with the AMSU-A weighting function distribution, which suggests that the channels can be used in the assimilation of satellite data in STMAS. With the case of the Typhoon Morakot in Taiwan Island in August 2009 as an example, experiments on assimilated and unassimilated AMSU-A radiances data were designed to analyze the impact of the assimilation of satellite data on STMAS. The results demonstrated that assimilation of AMSU-A data provided more accurate prediction of the precipitation region and intensity, and especially, it improved the 0–6h precipitation forecast significantly.     

Strategies for Reducing Computing Time of Nuclear Waste Management Simulations Using the PORFLOW ™ Software

PORFLOW is a CFD software which is dedicated to groundwater flow and nuclear waste management simulations. PORFLOW is developed by the ACRi company located in Bel Air (California, USA) (http://www.acricfd.com) with a subsidiary in Sophia Antipolis (France) (http://www.acri.fr). In this paper, we report on our recent efforts to improve the performances of PORFLOW by focusing our attention on the solution of the large sparse linear systems resulting from the space and time discretizations of flow and transport equations. This work is illustrated by applying PORFLOW to the calculation of the COUPLEX1 and part of the COUPLEX2 problems.     

Couplex Benchmark Computations Obtained with the Software Toolbox UG

This paper describes the numerical results for the COUPLEX benchmark obtained with the simulation software UG using vertex centered finite volume and higher order discontinuous Galerkin schemes. Multigrid solvers on unstructured grids, local mesh refinement and parallel computation are employed to yield very accurate solutions. Since the full range of results required in the benchmarks is too large to be displayed in this paper we focus on the comparison of discretization schemes, assessment of numerical errors and the presentation of parallel computations.     

多重网格方法在资料同化中的应用

多重网格方法通过在疏密不同的网格层上进行迭代，以平滑不同频率的误差分量。具有收敛速度快、精度高等优点，是现有计算工作量最少的算法之一。文章介绍了多重网格方法的原理和具体实现方案，将其用于卫星云图资料的变分同化中，并与超松弛迭代法作了比较。结果表明采用多重网格法，有助于提高同化效率。