“96·8”特大暴雨和中尺度系统发展结构的非静力数值模拟

1996年8月3～5日(“96     

关于控制迭代过程的结束条件

迭代法属于逐次逼近法,应用很广。它是求解函数方程、方程组和微分方程等的一种基本而重要的数值方法。在具体求解时,需要给出控制迭代过程结束的条件。现有的某些计算方法(或数值分析)的教材和书,在论述求解函数方程的计算方案、流程图(或框图)及算法语言程序时,常常将此条件取为相邻两次迭代值xk、xk-1的差小于给定的准确度本文指出这个条件是不完全的,讨论了该条件的正确表达式,由误差估计式给出了控制迭代过程结束条件的五种形式。     

台风"云娜"的热动力结构模拟试验研究

对2004年8月登陆我国的台风"云娜"的热动力结构及其演变进行了数值模拟。试验从大尺度环境场、台风中尺度动力和热力场三维结构、台风云雨特征等多角度,对模拟结果与实际观测进行了对比。结果表明,模式不仅较准确地再现了台风登陆前24h移动和演变的全过程,在移向、移速和登陆时间、地点上与实况十分接近,而且对台风系统中尺度结构、云雨特征的模拟与分析场及卫星观测也有较好的一致性。采用的模拟工具是新开发的双重嵌套AREM模式,该模式是静力模式,但对台风仍有较强的模拟能力。此外,本试验无任何人为台风模型的嵌入是其中的一个特点,从而增加了试验的可信度。     

http://www.sciencedirect.com/science/article/pii/S1674987111000508

This paper presents the results of a set of numerical models focussing on structural controls on hydrothermal mineralization.We first give an overview of natural phenomena of structurally-controlled ore formation and the background theory and mechanisms for such controls. We then provide the results of a group of simple 2D numerical models validated through comparison with Cu-vein structure observed near the Shilu Copper deposit(Yangchun,Guangdong Province, China) and finally a case study of 3D numerical modelling applied to the Hodgkinson Province in North Queensland(Australia).Two modelling approaches,discrete deformation modelling and continuum coupled deformation and fluid flow modelling,are involved.The 2D model-derived patterns are remarkably consistent with the Cu-vein structure from the Shilu Copper deposit,and show that both modelling approaches can realistically simulate the mechanical behaviours of shear and dilatant fractures.The continuum coupled deformation and fluid flow model indicates that pattern of the Cuveins near the Shilu deposit is the result of shear strain localization,development of dilation and fluid focussing into the dilatant fracture segments.The 3D case-study models(with deformation and fluid flow coupling) on the Hodgkinson Province generated a number of potential gold mineralization targets.     

库车坳陷克拉2气田异常流体压力演化史

根据克拉2气田现今实测地层压力分析了异常高压发育的特征,并通过流体包裹体测试求取包裹体形成时期的古压力。计算结果表明:下白垩统和古近系库姆格列木群砂岩储层在距今8～6Ma(康村组沉积期末)时,古压力系数大致为1.09～1.18,发育过剩压力为4.6～9.6MPa的异常高压。在此基础上,以古、今压力数据为约束条件,利用数值模拟技术定量恢复出下白垩统与侏罗系异常流体压力的演化历史。研究发现:下白垩统储层流体压力经历了沉积型异常流体压力形成发展(白垩纪-库车组沉积期末)和沉积型异常流体压力萎缩—构造挤压型异常流体压力孕育(库车组沉积期末至今)两个阶段。康村组沉积早期沉积型超压开始出现,库车组沉积期末达到顶峰,过剩压力达40～50MPa;此后,在沉积型超压萎缩—构造挤压型超压孕育过程的耦合下,过剩压力曾降至25～30MPa,更新世起超压略有升高。侏罗系流体压力演化较为简单,吉迪克组沉积期开始出现超压,至库车组沉积期末达到高峰,此后异常高压逐渐萎缩。在克拉2气田的主要成藏时期(距今3～1Ma),烃源岩较储层具有更高的异常压力,源—储压力差为10～30MPa。     

通量和平流形式预报方程的数值预报试验

本文利用1982年5月8日和1979年6月28日两个实例,进行了四种方案的48小时数值预报试验,比较了它们的预报质量及彼此间的差异。四种方案是:经典通量、通量修正、经典平流和平流修正。试验表明,在有地形的数值预报模式中,采用平流形式的预报方程组比采用通量形式的为好,既可有较高的计算精度,又可节省计算时间。     

Effect of bed surface roughness on longitudinal dispersion in artificial open channels

In several empirical and modelling studies on river hydraulics, dispersion was negatively correlated to surface roughness. In this study, it was aimed to investigate the influence of surface roughness on longitudinal dispersion under controlled conditions. In artificial flow channels with a length of 104 m, tracer experiments with variations in channel bed material were performed. By use of measured tracer breakthrough curves, average flow velocity, mean longitudinal dispersion, and mean longitudinal dispersivity were calculated. Longitudinal dispersion coefficients ranged from 0·018 m² s^?1 in channels with smooth bed surface up to 0·209 m² s^?1 in channels with coarse gravel as bed material. Longitudinal dispersion was linearly related to mean flow velocity. Accordingly, longitudinal dispersivities ranged between 0·152 ± 0·017 m in channels with smooth bed surface and 0·584 ± 0·015 m in identical channels with a coarse gravel substrate. Grain size and surface roughness of the channel bed were found to correlate positively to longitudinal dispersion. This finding contradicts several existing relations between surface roughness and dispersion. Future studies should include further variation in surface roughness to derive a better‐founded empirical equation forecasting longitudinal dispersion from surface roughness. Copyright © 2011 John Wiley & Sons, Ltd.     

Universal hydrofracturing algorithm for shear-thinning fluids: Particle velocity based simulation

A universal particle velocity based algorithm for simulating hydraulic fractures, previously proposed for Newtonian fluids, is extended to the class of shear-thinning fluids. The scheme is not limited to any particular elasticity operator or crack propagation regime. The computations are based on two dependent variables: the crack opening and the reduced particle velocity. The application of the latter facilitates utilization of the local condition of Stefan type (speed equation) to trace the fracture front. The condition is given in a general explicit form which relates the crack propagation speed (and the crack length) to the solution tip asymptotics. The utilization of a modular structure, and the adaptive character of its basic blocks, result in a flexible numerical scheme. The computational accuracy of the proposed algorithm is validated against a number of analytical benchmark solutions.     

Modelling the effect of wave overtopping on nearshore hydrodynamics and morphodynamics around shore-parallel breakwaters

Wave overtopping nearshore coastal structures, such as shore-parallel breakwaters, can significantly alter the current circulation and sediment transport patterns around the structures, which in turn affects the formation of tombolos and salients in the nearshore area. This paper describes the implementation of a wave overtopping module into an existing depth-averaged coastal morphological mode: COAST2D and model applications to investigate the effect of wave overtopping on the hydrodynamics and morphodynamics around a group of shore-parallel breakwaters. The hydrodynamic aspects of the model were validated against a series of laboratory conditions. The model was then applied to a study site at Sea Palling, Norfolk, UK, where 9 shore-parallel segmented breakwaters including 4 surface-piercing and 5 low-crested breakwaters are present, for the storm conditions in Nov 2006. The model results were compared with laboratory data and field measurements, showing a good agreement on both hydrodynamics and morphological changes. Further analysis of wave overtopping effect on the nearshore hydrodynamics and morphodynamics reveals that wave overtopping has significant impacts on the nearshore circulation, sediment transport and the resulting morphological changes within such a complex breakwater scheme under the storm and macro-tide conditions. The results indicate the importance of including the wave overtopping in modelling nearshore morphodynamics with the presence of coastal structures.