A temporally second-order accurate Godunov-type scheme for solving the extended Boussinesq equations     

J.B. Shiach  C.G. Mingham 《Coastal Engineering》2009

A numerical scheme for solving the class of extended Boussinesq equations is presented. Unlike previous schemes, where the governing equations are integrated through time using a fourth-order method, a second-order Godunov-type scheme is used thus saving storage and computational resources. The spatial derivatives are discretised using a combination of finite-volume and finite-difference methods. A fourth-order MUSCL reconstruction technique is used to compute the values at the cell interfaces for use in the local Riemann problems, whilst the bed source and dispersion terms are discretised using centred finite-differences of up to fourth-order accuracy. Numerical results show that the class of extended Boussinesq equations can be accurately solved without the need for a fourth-order time discretisation, thus improving the computational speed of Boussinesq-type numerical models. The numerical scheme has been applied to model a number of standard test cases for the extended Boussinesq equations and comparisons made to physical wave flume experiments.  相似文献   

船载物探电子设备的接地问题及解决方案     

韦成龙  罗文造  王立明  杨蜀冀 《海洋技术学报》2009,28(1)

船栽物探电子设备,是一个电气结合复杂的数据采集记录系统.各系统电子线路因船舶配电方式、工作环境的特殊性,其接地很重要,但容易被忽视.而接地处理不当,在使用时会受到噪声干扰、振荡;或造成仪器的输出榆入电路损坏,不能正常工作.更严重地甚至机壳漏电,危及人身安全.船栽物探电子设备接地问题.可分为"接地保护"和"接地故障"."接地保护"是将电气设备的金属外壳、支架和电缆的护套与大地等单位的金属船体(船壳)作永久性良好的电气连接,它分为保护接地和功能接地."接地故障"是指因设备不正确接地或由于绝缘层的破损、老化及环境恶劣引起绝缘电阻低下或者接地干扰等问题.文中从船载物探电子设备的使用特性出发.讨论在实际接地处理中出现问题的原因,分析了各种接地干扰的原理.并提出了相应的解决方案,对于船栽物探电子设备的实际安装应用有一定的指导意义.  相似文献   

双偏振雷达差分传播相移的五种滤波方法对比分析     

魏庆  ;胡志群  ;刘黎平 《成都信息工程学院学报》2014,(6):596-602

为了对双线偏振雷达的差分传播相移ΦDP进行滤波降噪处理,提高双线偏振雷达识别粒子、估测降水的能力。通过对C波段双线偏振雷达差分传播相移分别采用滑动平均、中值滤波、FIR滤波、卡尔曼滤波和小波分析5种不同方法地介绍,利用中国气象科学研究院灾害天气国家重点实验室移动式C波段双线偏振雷达(POLC)2013年外场试验观测资料,对比分析了经5种方法处理后的差分传播相移ΦDP径向数据,及其分别在降水估测中应用的效果。结果表明:采用FIR滤波方法、卡尔曼滤波方法和小波分析方法对ΦDP径向数据进行处理的效果比传统的滑动平均和中值滤波方法好;其中小波分析处理的ΦDP径向数据波动最小,再经过变距离法拟合的KDP数据的负值能够得到最有效地抑制,最后应用到R(KDP)降水估测中的精度最高。  相似文献   

Precise, automatic and fast method for vanishing point detection   总被引：1，自引：0，他引：1  

Mahzad  Kalantari  Franck  Jung  Jeanpierre  Guedon 《The Photogrammetric Record》2009,24(127):246-263

A new automated approach for vanishing point detection in images of urban scenes is described. This method is based on the theorem of Thales. The main contribution of this paper is the automatic and simultaneous detection of all vanishing points of the image, achieved by converting this problem into the detection of circles in a complex cloud of points, in which each point corresponds to a segment and is associated with an uncertainty. This extraction of circles uses a RANSAC method, modified to improve its speed by using accumulation techniques (Hough transform or otherwise). This robust estimation is then refined by least squares error propagation using the individual variances of each segment. The algorithm is robust, its accuracy is optimised and it is entirely automatic. Its successful operation has been tested on a large number of images of varied urban scenes.  相似文献   

局部角度域波传播步进算法研究   总被引：4，自引：2，他引：2       下载免费PDF全文

高静怀  周艳辉  毛剑  陈文超  吴如山  李幼铭 《地球物理学报》2007,50(1):248-259

本文从非均匀介质中波动方程出发,提出了基于一般标架的相空间（局部角度域）波传播的步进算法. 该方法在构造单程波的步进算法时,在选择标架或正交基等方面有更大的自由度. 我们以不随频率及深度变化且具有变尺度特性的Gabor_Daubechies紧标架为例,给出了单程波传输算子的具体形式及相应的波场步进算法;详细讨论了基于Gabor_Daubechies标架的传输算子的高频渐近展开问题,得出了在高频、小传输步长条件下传输算子的近似解析表达式,并给出使用条件. 通过模型算例,比较了精确传输算子与高频近似传输算子用于非均匀介质中波传播的结果,说明在一定条件下由两者得出的波场几乎是相同的.  相似文献   

基于分数阶时间导数的黏弹性衰减VTI介质中平面波传播     

程时俊  毛伟建  欧阳威 《地球物理学报》2021,64(3):965-981

目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推...  相似文献   

Modeling 3D thermal fracture propagation by transient cooling using virtual multidimensional internal bonds          下载免费PDF全文

K. Huang  A. Ghassemi 《国际地质力学数值与分析法杂志》2016,40(17):2293-2311

Thermal fracturing can play an important role in development of unconventional petroleum and geothermal resources. Thermal fractures can result from the nonlinear deformation of the rock in response to thermal stress related to cold water injection as well as heating. Before the rock reaches the final failure stage, material softening and bulk modulus degradation can cause changes in the thermo‐mechanical properties of the solid. In order to capture this aspect of the rock fracture, a virtual multidimensional internal bond‐based thermo‐mechanical model is derived to track elastic, softening, and the failure stages of the rock in response to the temporal changes of its temperature field. The variations in thermo‐mechanical properties of the rock are derived from a nonlinear constitutive model. To represent the thermo‐mechanical behavior of pre‐existing fractures, the element partition method is employed. Using the model, numerical simulation of 3D thermal fracture propagation in brittle rock is carried out. Results of numerical simulations provide evidence of model verification and illustrate nonlinear thermal response and fracture development in rock under uniform cooling. In addition, fracture coalescence in a cluster of fractures under thermal stress is illustrated, and the process of thermal fracturing from a wellbore is captured. Results underscore the importance of thermal stress in reservoir stimulation and show the effectiveness of the model to predict 3D thermal fracturing. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Variability of temperature,salinity and tidally‐averaged density in the middle estuary of the St Lawrence     

Langley R. Muir 《大气与海洋》2013,51(4):320-336

Abstract

Surveys in the Middle Estuary of the St Lawrence have yielded a data base consisting of more than 15,000 T‐S pairs distributed over 62 13‐h profiling stations. Although the T‐S curves at each station are remarkably linear, the variability of the slopes and intercepts of the lines is considerable. The means and standard deviations of the temperature and salinity at each individual station are not explicable in terms of linear combinations of the parameters for location in the Estuary, the upstream water properties, the phase of the spring‐neap cycle and the tidal energies.

It is shown that the tidally‐averaged density structure is separable into horizontal and vertical components and that its vertical variation over the whole Estuary may be explained by any one of three different functional forms. However, its horizontal variation is not explicable in terms of linear combinations of the parameters mentioned in the paragraph above.

Plots of the horizontal variations in temperature, salinity or density may only be meaningful if the data are collected synoptically, and even then cannot be considered to be accurate over time‐scales longer than one tidal cycle.  相似文献   

A ‘string function’ for describing the propagation of large-scale potential energy anomalies in a rotating fluid     

R. H. Tyler  R. Käse 《地球物理与天体物理流体动力学》2013,107(1-2):31-64

  相似文献   

Shear wave propagation along infinite slopes: A theoretically based numerical study     

C. di Prisco  M. Pastor  F. Pisanò 《国际地质力学数值与分析法杂志》2012,36(5):619-642

In this paper, the seismic response of ‘infinitely’ long slopes is numerically analysed via the formulation of a 1D analytical/numerical model, in which the soil mechanical behaviour is assumed to be elasto‐perfectly viscoplastic and simple shear (SS) kinematical constraints are imposed. In order to simplify the problem, a theoretically based procedure to set up a fully 1D shear constitutive model is defined, within which the mechanical response of a multiaxial relationship is condensed. The use of a 1D shear constitutive model is aimed at reducing the number of unknowns and, therefore, the computational costs. In particular, the case of the Mohr–Coulomb yield criterion is considered, while an enhanced Taylor–Galerkin finite element algorithm is employed to simulate the seismic wave propagation within the soil stratum. The proposed ‘condensation’/calibration procedure captures both the ‘pseudo’‐hardening pre‐failure behaviour and the influence of dilation on the occurrence of strain‐localization, which characterize, under SS conditions, the static response of virgin perfectly plastic soils. The effectiveness of the conceived method is shown with reference to freshly deposited deposits, while, in the case of highly overconsolidated strata, some difficulties arise because of the brittle behaviour induced both by unloading and non‐associativeness. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献