On the existence of atmospheric attractors

In this paper, we consider the long-time dynamics for the primitive equations of large-scale dry at- mosphere. First, by energy methods, we obtain the existence and uniqueness of global strong solu- tions of the problem. Second, by studying the long-time behavior of strong solutions, we construct a global attractor which captures all the trajectories.     

利用PETSc和FEPG编制海啸数值模拟程序的研究

海啸的数值模拟是海啸研究的一个重要领域, 它对于帮助理解海啸的基本物理特性和预防减灾具有重要意义。 海啸数值模拟程序的编制是一项繁杂的工作, 该文介绍了利用两种功能强大的通用软件(PETSc和FEPG)来进行海啸数值模拟程序编写的研究。 PETSc和FEPG采用有限差分、 有限元、 有限体积等多种离散方法, 可以对数值问题给出稳定的求解。 该文采用了直角坐标系下的非线性浅水波方程作为海啸波传播的控制方程进行离散求解, 并将其结果与TUNAMI N1模型进行了比较, 表明这两种方法方便而且有效。     

A survey of numerical methods for stochastic differential equations

The development of numerical methods for stochastic differential equations has intensified over the past decade. The earliest methods were usually heuristic adaptations of deterministic methods, but were found to have limited accuracy regardless of the order of the original scheme. A stochastic counterpart of the Taylor formula now provides a framework for the systematic investigation of numerical methods for stochastic differential equations. It suggests numerical schemes, which involve multiple stochastic integrals, of higher order of convergence. We shall survey the literature on these and on the earlier schemes in this paper. Our discussion will focus on diffusion processes, but we shall also indicate the extensions needed to handle processes with jump components. In particular, we shall classify the schemes according to strong or weak convergence criteria, depending on whether the approximation of the sample paths or of the probability distribution is of main interest.     

The SOMS model and its application to Lake Neuchâtel

A three-dimensional numerical circulation model (SOMS) based on primitive equations is described. The algorithm, by which Coriolis and vertical diffusion terms are treated implicitly while mass is still conserved exactly (algebraically), is discussed in detail. The model is applied to Lake Neuchâtel (Switzerland), to determine the general circulation under influence of the most prevailing wind.     

Development of the instantaneous unit hydrograph using stochastic differential equations

Recognizing that simple watershed conceptual models such as the Nash cascade ofn equal linear reservoirs continue to be reasonable means to approximate the Instantaneous Unit Hydrograph (IUH), it is natural to accept that random errors generated by climatological variability of data used in fitting an imprecise conceptual model will produce an IUH which is random itself. It is desirable to define the random properties of the IUH in a watershed in order to have a more realistic hydrologic application of this important function. Since in this case the IUH results from a series of differential equations where one or more of the uncertain parameters is treated in stochastic terms, then the statistical properties of the IUH are best described by the solution of the corresponding Stochastic Differential Equations (SDE's). This article attempts to present a methodology to derive the IUH in a small watershed by combining a classical conceptual model with the theory of SDE's. The procedure is illustrated with the application to the Middle Thames River, Ontario, Canada, and the model is verified by the comparison of the simulated statistical measures of the IUH with the corresponding observed ones with good agreement.     

Dynamic faulting under rate-dependent friction

We discuss the effects of rate-dependent friction on the propagation of seismic rupture on active faults. Several physicists using Burridge and Knopoff's box and spring model of faulting have proposed that fault complexity may arise from the spontaneous development of a self-similar stress distribution on the fault plane. If this model proves to be correct, it has important consequences for the origin of the complexity of seismic sources. In order to test these ideas on a more realistic earthquake model, we developed a new boundary integral equation method for studying rupture propagation along an antiplane fault in the presence of nonlinear rate-dependent friction. We study rupture dynamics of models with single and twin asperities. In our models, asperities are places on the fault with a higher value of prestress. Othewise all fault parameters are homogeneous. We show that for models with such asperities, a slip velocity weakening friction leads to the propagation of supersonic healing phases and to the spontaneous arrest of fracture if the prestress outside the asperities is low enough. For models with asperities, we can also observe narrow slip velocity pulses, qualitatively similar to the so-called Heaton pulses observed in some earthquake accelerograms. We also observe a complex distribution of stress after the rupture that depends on details of the initial distribution of asperities and on the details of the friction law.     

基于精确Zoeppritz方程的储层含油气性预测方法

流体指示因子和泊松比作为重要的指标参数在储层含油气性预测中发挥着至关重要的作用,大量学者开展了这两个参数的直接反演研究.然而,现有反演方法主要是以精确Zoeppritz方程的近似公式为正演方程,近似公式诸多的假设条件及较低的计算精度极大地限制了这类方法在复杂储层的应用效果.因此,为了提升储层含油气性预测精度,文中提出了一种新的基于精确Zoeppritz方程的流体因子和泊松比反演方法.首先,借助敏感流体因子定量分析法对现有流体指示因子进行优选,并将传统形式的精确Zoeppritz方程改写为包含该优选流体指示因子和泊松比的新形式.然后,基于新方程构建贝叶斯理论框架下的非线性反演目标函数.同时,为了进一步提升流体因子和泊松比对储层的刻画精度,在假设背景先验模型服从高斯分布的同时引入服从微分拉普拉斯分布的块约束项.最后,借助泰勒级数展开对上述非线性目标函数进行求解.合成数据和油田数据验证结果表明新方法能够稳定合理地估计流体指示因子和泊松比,且精度远高于基于近似公式的传统方法.此外,实际数据测试表明新方法能够有效提升储层含油气性预测精度,降低不确定性.

     

The tripole: A new coherent vortex structure of incompressible two-dimensional flows

Abstract

Using a contour dynamical algorithm, we have found rotating tripolar V-state solutions for the inviscid Euler equations in two-dimensions. We have studied their geometry as a function of their physical parameters. Their stability was investigated with the aid of contour surgery, and most of the states were found to be stable. Under finite-amplitude perturbations, tripoles are shown to either fission into two asymmetric dipoles or to evolve into a shielded axisymmetric vortex, demonstrating the existence of two new ‘‘reversible transitions'’ between topologically distinct coherent vortex structures. These dynamical results are confirmed by pseudo-spectral simulations, with which we also show how continuous tripolar long-lived coherent vortex structures can be generated in a variety of ways.     

港湾中海啸流的模拟及其时空演变特征对输入条件不确定性的响应分析

一直以来,海啸波特征作为表征海啸潜在破坏性的参数指标得到了广泛应用,特别是针对近场极端海啸事件造成的灾害来说,这种表征具有较好的适用性.然而总结分析历史海啸事件造成的损失发现：在远场近岸及港湾系统中,海啸诱导的强流却是造成损失的主要原因.陆架或港湾振荡导致海啸波幅快速升降诱发强流,可能促使港工设施受到威胁及损害,进而对海啸预警服务及海事应急管理提出了新的挑战.因此,全面理解与评估海啸在港湾中诱发的灾害特征,探索港湾中海啸流的数值模拟方法,发展针对港湾尺度的海啸预警服务指导产品尤为迫切.受限于海啸流验证数据的缺乏及准确模拟海啸流技术方法的诸多不确定性,大部分海啸数值模拟研究工作主要是针对水位特征的研究及验证,可能导致对港湾中海啸灾害危险性认识的曲解与低估.本研究基于非线性浅水方程,针对夏威夷群岛三个典型港湾建立了精细化海啸数值模型（空间分辨率达到10 m）,并联合有限断层破裂模型计算分析了日本东北地震海啸在三个港湾及其邻近区域的海啸特征,波、流计算结果与实测结果吻合较好,精细化的海啸港湾模型模拟结果可信.模拟发现港湾中较小的波幅,同样可以产生强流.综合分析日本东北地震海啸波、流特征对输入条件不确定性的响应结果发现：港湾中海啸波-流能量的空间分布特征差异较大,这与港湾系统中海啸波的驻波特性相关;相比海啸波幅空间特征,海啸流特征具有更强的空间敏感性;海啸流时空分布特征对输入条件的不确定性响应比海啸波幅对这些不确定性的响应更强,海啸流的模拟与预报更有挑战性;不确定性对海啸流计算精度的影响会进一步传导放大港湾海啸流危险性的评估及对港工设施产生的应力作用的误差,合理的输入条件对海啸流的精确模拟至关重要.最后,希望通过本文的研究可以从海啸波-流特征角度更加全面认识近岸海啸灾害特征,拓展海啸预警服务的广度与深度,从而为灾害应急管理部门提供更加科学合理的辅助决策产品.

     

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.