The thermodynamic structure atop a penetrating convective thunderstorm

The thermodynamic structure on top of a numerically simulated severe storm is examined to explain the satellite observed plume formation above thunderstorm anvils. The same mechanism also explains the formation of jumping cirrus observed by Fujita on board of a research aircraft. A three-dimensional, non-hydrostatic cloud model is used to perform numerical simulation of a supercell that occurred in Montana in 1981. Analysis of the model results shows that both the plume and the jumping cirrus phenomena are produced by the high instability and breaking of the gravity waves excited by the strong convection inside the storm. These mechanisms dramatically enhance the turbulent diffusion process and cause some moisture to detach from the storm cloud and jump into the stratosphere. The thermodynamic structure in terms of the potential temperature isotherms above the simulated thunderstorm is examined to reveal the instability and wave breaking structure. The plumes and jumping cirrus phenomena represent an irreversible transport mechanism of materials from the troposphere to the stratosphere that may have global climatic implications.     

Triple diagram method for the prediction of wave height and period

Many formulations have been developed so far to predict the wave height and period from fetch length and wind blowing duration for a constant wind speed. This study aimed to predict wave parameters from fetch length and meteorological factors by using triple diagram methodology based on Kriging principles. Proposed model results were compared with Joint North Sea Wave Project (JONSWAP) model which is used so commonly in the ocean and coastal engineering studies. For the implementation of the methodology hourly wave and wind data were obtained from a buoy located in Lake Ontario. Numerical and graphical comparisons demonstrated that the proposed method outperforms the classical formulation.     

Short-term shoreline changes due to cross-shore structures: a one-line numerical model

The proposed numerical model simulates the short-term temporal changes in shoreline position due to a structure interrupting the longshore sediment flux. The impacts of both the groin-type construction and underwater trench of arbitrary orientation relative to the shore are discussed. In order to estimate the sediment mass trapped by the structure, a submodel of the longshore sediment transport induced by a random wave field is developed. The contribution of the surface roller in momentum balance as well as in sediment suspension is included. The shoreline changes are computed from the equation deduced from the mass conservation. The perturbations in the longshore sediment discharge caused by a structure are assumed to concentrate within some boundary area of which the spatial scale is proportional to the structure's length until the latter is exceeded by the width of the sediment flux. It is shown in particular that the total effect of a long trench (channel) and a pier in its nearshore part results in general shoreline recession except for the vicinity of a pier. The model is tested against the laboratory data of Baidei et al. (1994) and applied to the Baidara Bay coast (Kara Sea) where a pipeline would be designed.     

Analytical Study on Wave Diffraction from a Vertical Circular Cylinder in Front of Orthogonal Vertical Walls

In this paper, the principle of mirror image is used to transform the problem of wave diffraction from a circular cylinder in front of orthogonal vertical walls into the problem of diffraction of four symmetric incident waves from four symmetrically arranged circular cylinders, and then the eigenfunction expansion of velocity potential and Grafs addition theorem are used to give the analytical solution to the wave diffraction problem. The relation of the total wave force on cylinder to the distance between the cylinder and orthogonal vertical walls and the incidence angle of wave is also studied by numerical computation.     

Transition from wavelets to ripples in a laboratory flume with a diverging channel

An experimental investigation on the initiation and development of bed forms on a bed of fine silica sand was conducted under alluvial flow conditions in a laboratory flume with a diverging channel. The main aims of the study were to assess： i） the steepness of bed forms in the transition stage of development; and ii） the threshold height of wavelets （ηt） that triggered the start of ripple development. Detailed bed profile measurements were carried out using an acoustic Doppler probe, traversed longitudinally over the sediment bed at various experimentation times. The bed form dimensions were extracted from such bed profile records and analysed for the wavelet, transition and equilibrium stages. It was found that the steepness of ripples in the transition and equilibrium stages were similar, confirming predictions of previous mathematical model simulations. A lognormal distribution fitted the wavelet length data. The wavelet threshold height was estimated as ηt ≈ 7 mm, or ηt≈ 80 in wall units. Such a height magnitude suggested that ripple development could be triggered by the wavelets reaching the outer flow zone of a turbulent boundary layer. The ηt value obtained corresponded generally to the intersection point between two predictive equations for bed form dimensions. A formulation was developed to predict ηt as a function of the sediment grain size, which was confirmed for the fine sand used in this study.     

Modelling linear wave transformation induced by dissipative structures—Random waves

The relevant theory is presented and numerical results are compared with the analytical solution for the interaction of non-breaking waves with an array of vertical porous circular cylinders on a horizontal bed. The extension to the cases of unidirectional and multidirectional waves is obtained by means of a transfer function. The influence of the mechanical properties of porous structures and wave irregularity on wave transformation is analysed. Results for unidirectional and multidirectional wave spectra are compared to those obtained for regular waves. The model presented reproduces well the analytical results and provides a tool for analysing several engineering problems.     

Wave interaction with a perforated wall breakwater with a submerged horizontal porous plate

This study examines the hydrodynamic performance of a new perforated-wall breakwater. The breakwater consists of a perforated front wall, a solid back wall and a submerged horizontal porous plate installed between them. The horizontal porous plate enhances the stability and wave-absorbing capacity of the structure. An analytical solution based on linear potential theory is developed for the interaction of water waves with the new proposed breakwater. According to the division of the structure, the whole fluid domain is divided into three sub-domains, and the velocity potential in each domain is obtained using the matched eigenfunction method. Then the reflection coefficient and the wave forces and moments on the perforated front wall and the submerged horizontal porous plate are calculated. The numerical results obtained for limiting cases are exactly the same as previous predictions for a perforated-wall breakwater with a submerged horizontal solid plate [Yip, T.L., Chwang, A.T., 2000. Perforated wall breakwater with internal horiontal plate. Journal of Engineering Mechanics ASCE 126 (5), 533–538] and a vertical wall with a submerged horizontal porous plate [Wu, J.H., Wan, Z.P., Fang, Y., 1998. Wave reflection by a vertical wall with a horizontal submerged porous plate. Ocean Engineering 25 (9), 767–779]. Numerical results show that with suitable geometric porosity of the front wall and horizontal plate, the reflection coefficient will be always rather small if the relative wave absorbing chamber width (distance between the front and back walls versus incident wavelength) exceeds a certain small value. In addition, the wave force and moment on the horizontal plate decrease significantly with the increase of the plate porosity.     

An empirical formula for friction coefficient of a perforated wall with vertical slits

The friction coefficient in the permeability parameter of a perforated wall has been estimated on the basis of a best fit between measured and predicted values of such hydrodynamic coefficients as reflection and transmission coefficients. In the present study, an empirical formula for the friction coefficient is proposed in terms of known variables, i.e., the porosity and thickness of the perforated wall and the water depth. This enables direct estimation of the friction coefficient without invoking a best fit procedure. To obtain the empirical formula, hydraulic experiments are carried out, the results of which are used along with other researchers' results. The proposed formula is used to predict the reflection and transmission coefficients of various types of structures including a perforated wall. The concurrence between the experimental data and calculated results is good, verifying the appropriateness of the proposed formula. It is also shown that the proposed formula can be used for irregular waves as well.     

On crustal corrections in surface wave tomography

Mantle models from surface waves rely on good crustal corrections. We investigated how far ray theoretical and finite frequency approximations can predict crustal corrections for fundamental mode surface waves. Using a spectral element method, we calculated synthetic seismograms in transversely isotropic PREM and in the 3-D crustal model Crust2.0 on top of PREM, and measured the corresponding time-shifts as a function of period. We then applied phase corrections to the PREM seismograms using ray theory and finite frequency theory with exact local phase velocity perturbations from Crust2.0 and looked at the residual time-shifts. After crustal corrections, residuals fall within the uncertainty of measured phase velocities for periods longer than 60 and 80 s for Rayleigh and Love waves, respectively. Rayleigh and Love waves are affected in a highly non-linear way by the crustal type. Oceanic crust affects Love waves stronger, while Rayleigh waves change most in continental crust. As a consequence, we find that the imperfect crustal corrections could have a large impact on our inferences of radial anisotropy. If we want to map anisotropy correctly, we should invert simultaneously for mantle and crust. The latter can only be achieved by using perturbation theory from a good 3-D starting model, or implementing full non-linearity from a 1-D starting model.     

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

目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推导出准压缩波qP、准剪切波qSV和纯剪切波SH的复速度、相速度、能量速度以及品质因子的解析表达式.对模型的正确性进行了数值验证,并分析了qP,qSV和SH波在介质中的传播特性.数值试验结果表明：本模型能够实现理想的恒定Q行为,表现了品质因子和速度的各向异性特征,显示出黏弹性增强将导致能量速度和相速度的频散曲线变化剧烈;速度和衰减各向异性参数与传播角度之间的耦合效应对qP,qSV和SH波的速度和能量影响明显;qP,qSV和SH波的频散曲线和波前面随着衰减各向异性强度的改变发生显著变化,其中耦合在一起的qP和qSV波变化趋势相同,而SH波与它们呈现相反的变化规律.本研究为从常Q模型角度分析地震波在衰减各向异性黏弹性介质中的传播特征奠定了理论基础.