水平非均匀性蒸发波导诊断及其对雷达探测的影响

蒸发波导是发生在海气边界层的一种异常折射现象,因为其分布广、发生概率大,所以被认为是对海上电子装备影响最为显著的波导类型。然而由于其形成机制复杂,且在近岸地区存在水平不均匀性,使得目前非均匀蒸发波导的诊断及其应用还未能落实到实际工作中。针对这一现状,首先利用G L Geernaert的方法修正了Monin-Obukhov相似理论,将其扩展到海洋大气表面边界层不均匀条件下;其次在Babin模式的基础上引入张强普适函数的非线性修正因子与阵性风速,从而将蒸发波导诊断模式的适用范围拓展到近海沿岸地区和甚低风速条件下。并在此基础上研究了蒸发波导水平非均匀性对雷达探测的影响,得到了水平非均匀蒸发波导能够改变均匀波导环境下雷达的探测距离及其盲区的分布。     

底部呈二维缓变情况下超大型浮体水弹性响应的两种计算方法比较

超大型浮体在海洋资源开发和海洋空间利用方面有重要应用前景.非均匀海洋环境中的水弹性响应是其应用中的一个重要问题.在近海中最典型的非均匀海洋环境当属由于底部变化引起的非均匀现象.本文分别采用多重尺度法（零阶近似）和常规的有限水深势流格林函数边界积分法,对底部呈二维缓变情况下超大型浮体的水弹性响应问题进行了研究和对比,并与实验工况进行了对照.两种方法与试验结果吻合较好,证明非均匀海洋环境确实对超大型浮体的水弹性响应具有一定的影响.     

Effect of the Planetary Boundary Layer Horizontal Inhomogeneity in the Two-Dimensional Fluid Motion Model

The generalized two-dimensional vortex equation is derived for an incompressible viscous fluid in a rotating system for a vertically averaged flow taking into account the variability of the boundary layer characteristics. The resulting equation contains parameters and their spatial derivatives determined by the second moments of functions describing the vertical profiles of the flow components. Numerical experiments demonstrate the influence of the boundary-layer horizontal inhomogeneity on the evolution of the vorticity field of a pair of atmospheric vortices.     

波速线性变化介质中地震波传播特征

根据变速介质中地震波传播所满足的变系数波动方程,推导出波速线性变化介质中一维波动方程的解析解,讨论了波速不均匀性对地震波传播的影响：地震波振幅衰减大小与介质速度梯度有关;频率转折点两侧地震波传播性质;频率与地震波波速关系。     

非均一相互作用能对超薄膜生长影响的Monte Carlo模拟研究

利用Monte Carlo (MC)模拟技术研究了非均一的吸附原子与基底相互作用能在一定的生长条件下对超薄膜生长过程的影响.非均一相互作用能是由基底表面原子在垂直和水平方向上实际位置与理想晶格原子位置的偏差所造成.本文用高斯分布来表示这种非均一相互作用能.模拟结果表明:非均一相互作用能对超薄膜的生长过程及薄膜的形貌有显著的影响.这种影响同时受到生长条件的限制,在中等温度时相互作用能的非均一性对岛的个数、平均大小的影响最显著;温度的增加在一定程度上可抵御相互作用能的非均一性对薄膜生长的影响.     

Seabed shear stress and bedload transport due to asymmetric and skewed waves

A simple conceptual formulation to compute seabed shear stress due to asymmetric and skewed waves is presented. This formulation generalizes the sinusoidal wave case and uses a variable friction factor to describe the physics of the boundary layer and to parameterize the effects of wave shape. Predictions of bed shear stresses agree with numerical computations using a standard boundary layer model with a k–ε turbulence closure. The bed shear stress formulation is combined with a Meyer-Peter and Müller-type formula to predict sheet flow bedload transport under asymmetric and skewed waves for a horizontal or sloping bed. The predictions agree with oscillatory water tunnel measurements from the literature.     

Bottom friction and its effects on periodic long wave propagation

A new set of Boussinesq-type equations describing the free surface evolution and the corresponding depth-integrated horizontal velocity is derived with the bottom boundary layer effects included. Inside the boundary layer the eddy viscosity gradient model is employed to characterize Reynolds stresses and the eddy viscosity is further approximated as a linear function of the distance measured from the seafloor. Boundary-layer velocities are coupled with the irrotational velocity in the core region through boundary conditions. The leading order boundary layer effects on wave propagation appear in the depth-integrated continuity equation to account for the velocity deficit inside the boundary layer. This formulation is different from the conventional approach in which a bottom stress term is inserted in the momentum equation. An iterative scheme is developed to solve the new model equations for the free surface elevation, depth-integrated velocity, the bottom stress, the boundary layer thickness and the magnitude of the turbulent eddy viscosity. A numerical example for the evolution of periodic waves propagating in one-dimensional channel is discussed to illustrate the numerical procedure and physics involved. The differences between the conventional approach and the present formulation are discussed in terms of the bottom frictional stress and the free surface profiles.     

Preliminary analysis of the applicability of adiabatic modes to inverting synthetic acoustic data in shallow water over a sloping sea floor

The May 2001 Geoacoustic Inversion Techniques Workshop provided synthetic transmission loss (TL) data for four cases with range-dependent shallow-water all-liquid environments. In two of these cases ("0" and "1"), the sea floor has constant slope and the geoacoustic model (GAM) is range independent. Cost functions have been computed using a new adiabatic-mode TL algorithm (which uses an exact velocity boundary condition at the sloping sea floor), as one parameter in the GAM is varied. Two frequencies (80 and 220 Hz) were selected. In case 0, the sea-floor slope is 0.0183 and the GAM comprises an inhomogeneous layer over a basement. The sea-floor sound-speed was selected as the variable parameter. The resulting cost minima at 80 and 220 Hz are displaced from the actual sound speed by 2.3 and 3.4 m/s, respectively. In case 1, the sea-floor slope is 0.012 and the GAM comprises one homogenous layer, five inhomogeneous layers, and a basement. The selected parameter was the sound-speed in the homogeneous layer. The corresponding cost minima are displaced by -1.2 and +1.1 m/s. The relative values of these four errors indicate that mode coupling increases with sea-floor slope and that there may be a dependence on frequency at the greater slope.     

浅化波浪层流边界层流速分布特性的数值分析

建立了同时考虑波致雷诺应力和时均水平压强梯度影响的二阶波浪边界层数学模型,模型计算得到的浅化波浪层流边界层内瞬时流速剖面、振荡速度幅值和时均流速剖面均与水槽实验数据吻合较好,在此基础上探讨了浅化波浪边界层流速分布特性及其影响机制。随着波浪的浅化变形,边界层内时均流速剖面"底部向岸、上部离岸"的变化特征越来越明显。这是二阶对流项引起的波致雷诺应力和离岸回流引起的时均水平压强梯度共同作用的结果,在床面附近由波致雷诺应力占主导作用并趋于引起向岸流动,在上部区域由时均水平压强梯度占主导作用并趋于引起离岸流动。     

海上风的边界层模式

由于航海、海上开发作业等对海洋上风和海浪的预报提出越来越高的要求,而海浪、风暴潮等海洋水文要素的数值计算和预报,迫切需要解决海洋上风场的精确计算。但是,复杂的海面结构,大气稳定度的影响以及风、浪之间动量的交换等,使海上风的理论计算遇到很多困难,至今大部分工作是依靠统计方法。利用天气预报的形势场计算地转风或梯度风,以及它与海面摩擦、大气稳定度的经验订正关系。     

Longitudinal dispersion due to the boundary layer in an oscillatory current: Theoretical analysis in the case of an instantaneous line source

The longitudinal dispersion effect due to the boundary layer formed by a tidal oscillatory current is examined theoretically. This analysis reveals the process whereby the dispersion coefficient becomes steady after the release of the diffusing substance. Though the dispersion due to an oscillatory current has so far been investigated mostly in the case of a linear velocity profile, the following result was found by taking account of the boundary layer in the oscillatory current. The depth-averaged dispersion coefficient for the case of a current having a boundary layer can be a few times larger than in the case of a linear velocity profile when the characteristic mixing time is long; the phase lead in the boundary layer induces nearly 20 percent of the longitudinal dispersion effect.     

Permeability in Rotliegend gas sandstones to gas and brine as predicted from NMR,mercury injection and image analysis

Permeability characterisation of low permeability, clay-rich gas sandstones is part of production forecasting and reservoir management. The physically based Kozeny (1927) equation linking permeability with porosity and pore size is derived for a porous medium with a homogeneous pore size, whereas the pore sizes in tight sandstones can range from nm to μm. Nuclear magnetic resonance (NMR) transverse relaxation was used to estimate a pore size distribution for 63 samples of Rotliegend sandstone. The surface relaxation parameter required to relate NMR to pore size is estimated by combination of NMR and mercury injection data. To estimate which pores control permeability to gas, gas permeability was calculated for each pore size increment by using the Kozeny equation. Permeability to brine is modelled by assuming a bound water layer on the mineral pore interface. The measured brine permeabilities are lower than predicted based on bound water alone for these illite rich samples. Based on the fibrous textures of illite as visible in electron microscopy we speculate that these may contribute to a lower brine permeability.     

Nonlinear response of membranes to ocean waves using boundary and finite elements

The behavior of a highly deformable membrane to ocean waves was studied by coupling a nonlinear boundary element model of the fluid domain to a nonlinear finite element model of the membrane. The hydrodynamic loadings induced by water waves are computed assuming large body hydrodynamics and ideal fluid flow and then solving the transient diffraction/radiation problem. Either linear waves or finite amplitude waves can be assumed in the model and thus the nonlinear kinematic and dynamic free surface boundary conditions are solved iteratively. The nonlinear nature of the boundary condition requires a time domain solution. To implicitly include time in the governing field equation, Volterra's method was used. The approach is the same as the typical boundary element method for a fluid domain where the governing field equation is the starting point. The difference is that in Volterra's method the time derivative of the governing field equation becomes the starting point.The boundary element model was then coupled through an iterative process to a finite element model of membrane structures. The coupled model predicts the nonlinear interaction of nonlinear water waves with highly deformable bodies. To verify the coupled model a large scale test was conducted in the OH Hinsdale wave Research Laboratory at Oregon State University on a 3-ft-diameter fabric cylinder submerged in the wave tank. The model data verified the numerical prediction of the structure displacements and of the changes in the wave field.The boundary element model is an ideal modeling technique for modeling the fluid domain when the governing field equations is the Laplace equation. In this case the nonlinear boundary element model was coupled with a finite element model of membrane structures, but the model could have been coupled with other finite element models of more rigid structures, such as a pontoon floating breakwater.     

潮滩水流加减速阶段边界层流速剖面的比较

江苏大丰地区潮滩由于水深较浅,潮流、波浪等动力较强,整个水层可视为边界层,其主体部分是对数层,即水流流速在垂向上呈对数分布。在潮流的加减速阶段,流速剖面将可能偏离对数分布,从而使对数剖面法计算出来的边界层参数造成误差。使用MIDAS-400用户化数据采集系统,在大丰潮滩获得了多层流速、浊度等同步高频观测数据,基于修正后的von Karman-Prandtl模型对u-lnz进行回归分析、数据内部一致性分析来定义流速对数剖面并与未修正前经典理论得到的边界层参数进行比较。分析结果表明,修正后的流速剖面更符合实际情况,边界层参数除了受水流加速度的影响外,还和沙纹等因素有关。另外,边界层参数的变化量与特征加速度负相关。     

On acoustic N-wave reflections from atmospheric layered inhomogeneities

The nonlinear propagation of acoustic pulses from a point source of an explosive character (surface explosion or volcano) throughout the atmosphere with stratified wind-velocity and temperature inhomogeneities is studied. The nonlinear distortions of acoustic pulse and its transformation into an N-wave during its propagation to the upper atmosphere are analyzed in the context of a modified Burgers’ equation which takes into account a geometric ray-tube divergence simultaneously with an increase in both nonlinear and dissipative effects with height due to a decrease in atmospheric density. The problem of reflection of a spherical N-wave from an atmospheric inhomogeneous layer with model vertical wind-velocity and temperature fluctuations having a vertical spectrum that is close to that observed within the middle atmosphere is considered. The relation between the parameters (form, length, frequency spectrum, and intensity) of signals reflected from an atmospheric inhomogeneous layer and the parameters of the atmospheric fine layered structure at reflection heights is analyzed. The theoretically predicted forms of signals reflected from stratified inhomogeneities within the stratosphere and the lower thermosphere are compared to the observed typical forms of both stratospheric and thermospheric arrivals from surface explosions and volcanoes in the zones of an acoustic shadow.     

Experimental and numerical investigation of viscous effects on solitary wave propagation in a wave tank

Two-dimensional depth-averaged Boussinesq-type equations were presented with the consideration of slowly varying bathymetry and effects of bottom viscous boundary layer. These Boussinesq-type equations were written in terms of the horizontal velocity components evaluated at an arbitrary elevation in the water depth and the free surface displacement. The leading order effects of the bottom boundary layer were represented by a convolution integral in the depth-integrated continuity equation. To test the validity of the theory, a set of laboratory experiments was performed to measure the viscous damping and shoaling of a solitary wave propagating in a wave tank. The time histories of the free surface profiles were measured at several locations along the centerline of the flume. To compare these laboratory data with theoretical results, the two-dimensional Boussinesq-type equations were integrated across the wave tank, resulting in a set of one-dimensional equations, while the side-wall boundary layers were properly considered. The agreement between the experimental data and numerical results was very good. The bottom shear stress formula was also given and its impact on the sediment transport rate was discussed.     

Modeling a neutrally stratified turbulent air flow over a horizontal rough surface

The neutrally stratified boundary layer over a smooth rough surface is consider. The turbulent flow is simulated using a finite-difference eddy-resolving model of the atmospheric boundary layer (ABL). The model includes different turbulence closure schemes and numerical approximations for advection components of the momentum balance equation. We investigate the quality of reproduction of spectral characteristics of the turbulent flow and the model’s capabilities to reproduce the observed profile of mean wind velocity near the rough surface. It is shown that the best result is obtained by coupling a numerical scheme of higher order of accuracy with a mixed closure scheme based on an adaptive estimation of the mixing length for subgrid-scale fluctuations. Here, we are able to reproduce the asymptotics of the fluctuation spectrum of the longitudinal component of wind velocity near the surface and within the boundary layer as well as the logarithmic profile of mean velocity near the surface.     

A theory of semigeostrophic gravity waves and its application to the intrusion of a density current along a coast

Semigeostrophic gravity waves associated with a coastal boundary current, which has finite and uniform potential vorticity and is bounded away from the coastline by a density front on the ocean surface, are investigated. It is shown that the semigeostrophic coastal current has two waves which are named here the Semigeostrophic Coastal Wave (SCW) and the Semigeostrophic Frontal Wave (SFW). The SCW becomes an elementary Kelvin wave at some limit while the SFW is caused by the existence of the surface density front. The SCW appears mainly as variations in the upper layer depth at the coast and as alongshore velocity at the density front. On the other hand, the SFW appears mainly as variations in the width of the current. When the weak nonlinearity and ageostrophic effect are included, these semigeostrophic gravity waves satisfy the Kortweg- de Vries equation, which suggests that the local changes in the width and/or velocity of the semigeostrophic coastal current propagate as wave-like disturbances.     

Full-range equation for wave boundary layer thickness

Full-range equation covering all the flow regimes in a wave boundary layer is proposed for the boundary layer thickness. The results are compared with the available experimental data and good agreement has been found. In case of wave boundary layers, there are three definitions of boundary layer thickness in use. Therefore, the full-range equation is derived for three of the definitions. The findings of this study may be useful in calculating suspended sediment transport in coastal environments and studying wave–current combined motion.