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.     

略论目标规划及交互方法在小康示范村规划中的实践--以西藏贡嘎县沃拉村为例

建设全面小康社会是我国本世纪头20年的核心内容,目前已经有越来越多的人加入到这一研究领域中,该领域的规划研究工作正在逐渐加强。本文以西藏贡嘎县沃拉村为例,介绍了目标规划及交互方法小康示范村规划中的实践,讨论了规划的过程、建模方法及模型的精度等问题,认为在研究人员的指导下,由村民学习并建模的方法是今后小康村规划的发展方向。     

黄土动力学数据处理软件的研制

从数据分析和处理流程和功能模块设计方面讨论了应用Visual Basic进行土动力学数据处理系统开发时的问题及可行的解决方案，介绍了编制的土动力学数据处理系统的主要功能。     

风机基础TLP平台在波浪中运动性能的CFD数值分析

随着海上风能的开发向深水发展,支撑风机的载体平台越来越受到关注。在经济性与安全性、稳定性的多重要求下,张力腿平台（TLP）在海洋风能资源的开发中体现出了重要地位。采用基于开源平台OpenFOAM开发的计算流体动力学（CFD）水动力学求解器naoe-FOAM-SJTU对一座处于中等水深下的风机基础水下TLP（STLP）的运动响应进行了数值模拟与研究。文中使用弹簧锚链模型模拟STLP的垂向系泊锁链系统,模拟该平台在不同波浪环境下的运动响应情况。首先将STLP单自由度自由衰减CFD模拟结果与已有全耦合时域分析结果进行对比,验证了naoe-FOAM-SJTU求解器及使用弹簧模型模拟STLP系泊系统的准确性与可靠性。随后在考虑非线性波浪载荷的情况下研究极端海况下与一般作业海况下STLP的运动响应情况,计算工况中的风机基础所受弯矩及锚链受力情况,并详细展示流场、速度场信息,分析高阶波浪成分、不同海况等条件对于STLP运动性能的影响。研究结果表明,TLP在中等水深中具有良好的运动性能,naoe-FOAM-SJTU求解器可以有效模拟水中生产平台在波浪环境下的水动力问题,并可以对整个流场进行可视化展示与分析。     

Soft computing approach for real-time estimation of missing wave heights

This paper presents soft computing approach for estimation of missing wave heights at a particular location on a real-time basis using wave heights at other locations. Six such buoy networks are developed in Eastern Gulf of Mexico using soft computing techniques of Artificial Neural Networks (ANN) and Genetic Programming (GP). Wave heights at five stations are used to estimate wave height at the sixth station. Though ANN is now an established tool in time series analysis, use of GP in the field of time series forecasting/analysis particularly in the area of Ocean Engineering is relatively new and needs to be explored further. Both ANN and GP approach perform well in terms of accuracy of estimation as evident from values of various statistical parameters employed. The GP models work better in case of extreme events. Results of both approaches are also compared with the performance of large-scale continuous wave modeling/forecasting system WAVEWATCH III. The models are also applied on real time basis for 3 months in the year 2007. A software is developed using evolved GP codes (C++) as back end with Visual Basic as the Front End tool for real-time application of wave estimation model.     

Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

线性规划在徐州地下水资源开发中的应用

本文从地下水系统的观点出发,对徐州地区中寒武统—中奥陶统碳酸盐岩地下水系统进行了分析。并应用线性规划的原理和方法对徐州市西郊四个水源地地下水进行了水文地质计算。在规定水位允许降深的约束条件下,求得至2000年,该区地下水最大开采量为14万米~3/天及各水源地相应取水量的最佳分配。文中采用特定条件图解法计算多边界、多主井抽水干扰影响的方法,具有通用意义。     

投影寻踪模型在地下水水质评价中的应用

文中分析了常用的水体质量评价模型的缺陷,提出了基于序列二次规划算法SQP(Sequentialquadraticprogramming)的投影寻踪评价模型,并利用多项式插值的方法进行水质分类,取得了符合客观实际的评价结果。某水源地地下水质量评价结果表明:用投影寻踪方法可以揭示高维数据的结构特征;所建立的评价模型不仅简单、直观和易于理解,而且结果合理。为水源地的保护和合理利用提供了决策依据。     

非线性最小二乘参数平差迭代算法

在非线性最小二乘问题现有的3类主要算法－－高斯－牛顿法、阻尼最小二乘法和最小二乘的拟牛顿法的基础上，引入了综合性能更优的非线性规划的SQPM（序列二次规划法）算法，并且为进一步提高SQPM算法迭代的收敛性，对其步长策略进行了改进。改进的SQPM算法成为无需精确计算参数概略值的非线性最小二乘参数平差的实用和有效算法。     

非结构网格上的三维浅水流动数值模型

针对当前复杂环境水流模拟的需求,建立了新型的基于特征型高分辨率数值算法的三维非结构网格浅水动力模型。模型采用有限体积法离散sigma坐标下的三维浅水方程,运用Roe黎曼近似解评估水平界面通量。模型网格拟合边界能力强,可根据需要局部加密;格式数值性能优良,具有守恒性、单调迎风性、高数值分辨率等特性。同时,应用干湿判别法处理动边界,以适应浅滩地形漫/露过程模拟的需要。封闭水池内部风生环流、干河床上溃坝过程和长江口实际潮流场的模拟从不同侧面展示了模型的特点,结果表明它能够准确地预测水流的三维流动结构,而且计算简单高效,具有良好的数值稳定性。