初始倾斜跃层所致内波的数值研究

The pycnocline in a closed domain is tilted by external wind forcing and tends to restore to a level posi- tion when the wind falls. An internal seiche oscillation exhibits if the forcing is weak, otherwise internal surge and internal solitary waves emerge, which serve as a link to cascade energy to small-scale processes. A two-dimensional non-hydrostatic code with a turbulence closure model is constructed to extend previous laboratory studies. The model could reproduce all the key phenomena observed in the corresponding labo- ratory experiments. The model results further serve as a comprehensive and reliable data set for an in-depth understanding of the related dynamical process. The comparative analyses indicate that nonlinear term favors the generation of internal surge and subsequent internal solitary waves, and the linear model predicts the general trend reasonably well. The vertical boundary can approximately reflect all the incoming waves, while the slope boundary serves as an area for small-scale internal wave breaking and energy dissipation. The temporal evolutions of domain integrated kinetic and potential energy are also analyzed, and the results indicate that about 20% of the initial available potential energy is lost during the first internal wave breaking process. Some numerical tactics such as grid topology and model initialization are also briefly discussed.     

海洋沉积样品~(14)C测年差异来源及其结果解释

~(14)C测年已为海洋研究工作广泛应用,但由于时有不合理取样或在运用数据时没有充分考虑各种海洋环境和动力因素,因而产生无法使用数据作合理解释的现象。本文试从海洋样品本身和海洋各种营力作用分析各种数据差异的来源,认为海洋环境特别是海岸带环境的复杂性是造成海洋样品~(14)C测年数据差异的主要原因,它可能影响样品的代表性或使数据解释时得出错误结论.并指出在运用年代结果时一定要充分考虑地质地理和各种营力作用因素。     

On the divergence-form finite-difference approximation to the momentum advection in curvilinear orthogonal coordinates

The traditional finite-difference schemes for the dynamical equations in curvilinear orthogonal coordinates have a basic flaw: they conserve only energy, but not momentum. A finite-difference approximation to these equations is suggested that conserves both energy and momentum.     

Dynamic Properties and Energy Conversion Efficiency of A Floating Multi-Body Wave Energy Converter

The present study proposed a floating multi-body wave energy converter composed of a floating central platform, multiple oscillating bodies and multiple actuating arms. The relative motions between the oscillating bodies and the floating central platform capture multi-point wave energy simultaneously. The converter was simplified as a forced vibration system with three degrees of freedom, namely two heave motions and one rotational motion. The expressions of the amplitude-frequency response and the wave energy capture width were deduced from the motion equations of the converter. Based on the built mathematical model, the effects of the PTO damping coefficient, the PTO elastic coefficient, the connection length between the oscillating body and central platform, and the total number of oscillating bodies on the performance of the wave energy converter were investigated. Numerical results indicate that the dynamical properties and the energy conversion efficiency are related not only to the incident wave circle frequency but also to the converter's physical parameters and interior PTO coefficients. By adjusting the connection length, higher wave energy absorption efficiencies can be obtained. More oscillating bodies installed result in more stable floating central platform and higher wave energy conversion efficiency.     

Finite element modeling of heating phenomena of cracks excited by high-intensity ultrasonic pulses

A three-dimensional thermo-mechanical coupled finite element model is built up to simulate the phenomena of dynamical contact and frictional heating of crack faces when the plate containing the crack is excited by high-intensity ultrasonic pulses. In the finite element model, the high-power ultrasonic transducer is modeled by using a piezoelectric thermal-analogy method, and the dynamical interaction between both crack faces is modeled using a contact-impact theory. In the simulations, the frictional heating taking place at the crack faces is quantitatively calculated by using finite element thermal-structural coupling analysis, especially, the influences of acoustic chaos to plate vibration and crack heating are calculated and analysed in detail. Meanwhile, the related ultrasonic infrared images are also obtained experimentally, and the theoretical simulation results are in agreement with that of the experiments. The results show that, by using the theoretical method, a good simulation of dynamic interaction and friction heating process of the crack faces under non-chaotic or chaotic sound excitation can be obtained.     

有限区域海面风场数值模式的建立

目前,国内外大尺度天气预报数值模式已经取得很大成就,如在预报地面气旋系统的移动和发展等方面。然而,这些模式的铅直分辨率一般不高,通常不包括边界层的动力过程,因而,在预报海面风场上,数值天气预报模式还没有令人满意的效果。目前所用的边界层预报模式尽管具有较高的水平分辨率和垂直分辨率,但要求大量的计算和资料,很难用于日常业务预报;另一种边界层诊断模式,用到的风场数据是气象工作者从天气形势(大尺度)分析中,通过差值等方法获得的,缺乏针对性,而且精度较低(WMO,1988;WMO,1990)。 作者认为,适于海洋要素预报的风场模式应是定位于有限区域,或称中尺度模式,并根据台站的具体资料来源和工作状况力争建立一套适于海洋要素预报的中尺度风场数值模式。鉴于以上,作者建立了一个适于海洋要素预报的有限区域海面风场数值模式。模式采用符合动力及热力学条件的简化方程组对大尺度风场进行加密,通过数值模拟得到适于海洋要素预报的有限区域海面风场。并将该模式运用于渤海这一有限区域,获得了良好的效果,说明该模式对海洋要素预报具有很好的应用前景。     

Stochastic systems with cross-correlated Gaussian white noises

This paper theoretically investigates three stochastic systems with cross-correlation Gaussian white noises. Both steady state properties of the stochastic nonlinear systems and the nonequilibrium transitions induced by the cross-correlated noises are studied. The stationary solutions of the Fokker-Planck equation for three specific examples are analysed. It is shown explicitly that the cross-correlation of white noises can induce nonequilibrium transitions.     

特殊动力作用在海洋生物中的研究进展

由摄食而引起的特殊动力作用(SDA)使动物在处理和转化食物时需要额外消耗能量。阐述了SDA的含义及其研究参数,讨论了目前国内外对SDA研究的主要途径及不同的研究重点,探讨了影响SDA耗能的主要因素,指出了SDA与动物生长之间的关系。从两方面的研究角度考虑,展望了研究SDA的前景和意义。     

粤东甲子海域潮波异常和南海北部潮波的传播

从南海潮波数值模拟潮能通量的结果中,勾勒出南海北部潮波的传播路径,显示出众多的路径分支。分析表明南海北部陆架包括北部湾海区,仅获得太平洋传入潮能的一小部分,因而南海北部沿岸的潮汐潮流都表现得比较弱。由于复杂的地形,导致了潮波传播路径指向不同方向。文章通过实测资料以及潮波的传播路径,特别讨论了粤东甲子站附近的潮波异常现象,并指出流向台湾海峡南口的分支和流向珠江口、广州湾的分支之间的潮能辐散,可能是造成甲子站附近潮性系数特别大以及潮汐潮流性质迥异现象的重要原因之一。     

黑潮延伸体海域典型涡旋的次中尺度特征分析

本文基于卫星遥感资料和高分辨率ROMS(Regional Ocean Modeling System)数值模拟结果, 对黑潮延伸体海域典型中尺度涡旋的次中尺度特征进行了探讨。卫星观测和模拟结果显示, 黑潮延伸体涡旋海域伴随着活跃的次中尺度现象。涡旋演变与多尺度能量分析结果表明, 涡旋海域次中尺度动能的强弱与涡旋海域地转流动能有着密切联系, 锋生可能是涡旋边缘次中尺度动能增强的重要机制。次中尺度现象在中尺度涡旋海域具有沿地转流方向的复杂涡丝状结构特征, 意味着涡旋边缘较强的水平浮力梯度和地转流侧向剪切为次中尺度过程形成与发展提供了有利条件。此外, 垂向结构分析表明, 次中尺度过程能引起较大的垂向速度, 最大可达100m·day^-1, 该垂向速度可以影响至混合层下200m深度处, 对海洋内部的垂向物质能量交换、海—气相互作用等有着重要的影响。     

Theoretical models of the height of the atmospheric boundary layer and turbulent entrainment at its upper boundary

The planetary boundary layer (PBL), which directly interacts with the underlying surface, differs significantly in its nature from the low-turbulent stably stratified free atmosphere. Fluctuations of the Earth’s surface heat balance immediately affect the PBL and assimilate there owing to the effective mechanism of turbulent heat transfer. In this case the upper boundary of the PBL plays the role of a cover, preventing the direct penetration of thermal effects and contaminants into an overlying atmospheric layer. In view of this, air pollution is especially dangerous when associated with shallow PBL. In addition, local peculiarities of climate change are mainly determined by the PBL height due to the high sensitivity of thin stably stratified PBLs to the thermal effects. Deep convective PBLs are not very sensitive to weak thermal effects, but they significantly affect the formation of convective cloudiness and the climate system as a whole by means of the turbulent entrainment of the thermal energy, humidity, aerosols, and other admixtures through the upper boundary. The PBL height and turbulent entrainment must be calculated when simulating and forecasting air pollution, abnormal frosts and heat, and other hazardous phenomena. In this paper we discuss the state-of-the-art knowledge in the area of PBL height simulation and suggest a new model of turbulent entrainment for convective PBLs.     

南海西北部上层海洋对台风“杜苏芮”的响应分析

基于锚碇观测资料,本文分析了南海西北部陆坡区上层海洋对台风“杜苏芮”的动力学和热力学响应特征。在动力学响应方面,台风“杜苏芮”期间上层流速显著增强,混合层纬向流速可达1.20 m/s;“杜苏芮”经过后上层海水运动以近惯性振荡为主(流向顺时针旋转周期在36～40 h之间)。近惯性能量在垂向分布上存在两个高值中心,分别位于混合层和温跃层深度上。近惯性能量耗散过程的e折时间尺度约为3.7 d,我们认为能量的向下传播在局地近惯性能量衰减过程中起主要作用。对能量谱的分析表明,“杜苏芮”作用期间近惯性频率能量相对于其作用前增大了约29.4倍,而全日和半日频率(K₁和M₂)能量有所减弱。此外,能量谱显示近惯性频率存在明显的“蓝移”现象,即对于纬向和经向流速分量在400 m以浅平均的近惯性振荡频率分别为1.167 f₀和1.170 f₀(f₀为局地惯性频率)。蓝移与近惯性内波的向下传播及正的相对涡度的输入有关。在热力学响应方面,上层海洋在台风的搅拌作用下,40～250 m深度均出现较小增温,最大增温幅度接近1°C;此外70 m以浅盐度的降低可能与台风过境时的降水相关,而Ekman抽吸引起的上升流则可能对70～100 m深度盐度的升高具有重要作用。     

Bending localization of nitrous oxide under anharmonicity and Fermi coupling: the dynamical potential approach

This paper studies the vibrational nonlinear dynamics of nitrous oxide with Fermi coupling between the symmetric stretching and bending coordinates by classical dynamical potential approach. This is a global approach in the sense that the overall dynamics is evidenced by the classical nonlinear variables such as the fixed points and the focus are on a set of levels instead of individual ones. The dynamics of nitrous oxide is demonstrated to be not so much dependent on the excitation energy. Moreover, the localized bending mode is shown to be ubiquitous in all the energy range studied.     

A combined wind and wave energy-converter concept in survival mode: Numerical and experimental study in regular waves with a focus on water entry and exit

A combined wind and wave energy converter concept, named STC concept was proposed. Model tests were performed in terms of operational and survival modes. Water entry and exit phenomena as well as green water on deck were observed during the survivability model tests. In this paper, a nonlinear numerical model based on a blended station-keeping potential-flow solver with a local impact solution for bottom slamming events and an approximated model for the water shipped on the deck is proposed to simulate these nonlinear phenomena. Physical investigation of the water entry and exit process was firstly carried out and uncertainty analysis of the model test results were performed. Numerical comparisons between the nonlinear solver and model test results are then performed in terms of mean, wave frequency and double wave frequency motion response components. The slamming and green water involved in the water entry process are specially investigated, in terms of the physical evolution and the effects on the dynamic motion responses. The validation work on the occurrence of slamming and water on deck as well as the slamming pressure are performed.     

Sub-inertial oscillations in the Black Sea generated by the semidiurnal tidal potential

The Black Sea shelf is a region of intense manifestation of various dynamical processes. Under the influence of different natural forces, eddy-wave phenomena develop here, which influence the general circulation of sea waters, biological productivity, and the condition of the engineering structures. Modern numerical models allow us to simulate and analyze the processes of the joint dynamics of marine circulation and large-scale waves. In this work, we study the spatiotemporal spectral characteristics of the sea level and velocity fluctuations formed due to atmospheric forcing and tidal potential. The hydrophysical fields are calculated using the Institute of Numerical Mathematics, Russian Academy of Sciences (INM RAS), σ model based on primitive equations. We use the CORE data as atmospheric forcing at the sea surface; the tidal potential is described by the semidiurnal lunar constituent M₂. Analyzing the simulation results makes it possible to emphasize that accounting for the semidiurnal tidal potential not only improves the accuracy of the sea-level calculation at coastal stations, but also generates subinertial baroclinic oscillations previously found in the Black Sea from the data of in situ observations.     

Calculation of energy of the barotropic-baroclinic interaction in the strait of Bab el Mandeb

Equations for the kinetic, potential, and total energy in the curvilinear boundary-fitted coordinates are considered in the strait area and adjacent regions of the Red Sea and Gulf of Aden. The budgets of different forms of energy are presented, making it possible to assess the correctness of the formulation of the three-dimensional boundary-value problem, the accuracy of the calculation of the barotropic-baroclinic interaction in an energy norm, and the relationship between the components of energy balance. The energy fluxes in the boundary sections of the strait are estimated for the interaction of thermohaline currents with main tidal harmonics and with the total tide, as well as for the amplitude of fluctuations in boundary energy fluxes in a synodic month. The results illustrating the behavior of some dynamical characteristics and components of energy balance at vertical and horizontal sections of the strait are given.     

A note on a critical wind speed for air–sea boundary processes

Wind and wind-generated waves were measured in a wind-wave tank. A clear transition was found in the relation between the wind speed U ₁₀ and the wind friction velocity u _* near u _* = 0.2 m/s, where U ₁₀ is the wind speed at 10 m height extrapolated from the measured wind profile in a logarithmic layer, and u _* = 0.2 m/s corresponds roughly to U ₁₀ = 8 m/s in the present measurement. Quite a similar transition was found in the relation between the spectral density of high frequency wind waves and u _*. These results suggest the existence of the critical wind speed for air–sea boundary processes, which was proposed by Munk (J Marine Res 6:203–218, 1947) more than half a century ago. His original idea of the critical wind speed was based on the discontinuities in such phenomena as white caps, wind stress, and evaporation, which commonly appear at a wind speed near 7 m/s. On the basis of the results of our present study and those of earlier studies, we discuss the phenomena which are relevant to the critical wind speed for the air–sea boundary processes. The conclusion is that the critical wind speed exists and it is attributed to the start of wave breaking rather than the Kelvin–Helmholtz instability, but the air–sea boundary processes are not discontinuous at a particular wind speed; because of the stochastic nature of breaking waves, the changes occur over a range of wind speeds. Detailed discussions are presented on the dynamical processes associated with the critical wind speed such as wind-induced change of sea surface roughness and high frequency wave spectrum. Future studies are required, however, to clarify the dynamical processes quantitatively. In particular, there is a need to further examine the gradual change of breaking patterns of wind waves with the increase of wind speed, and the associated change of the structure of the wind over wind waves, such as separation of the airflow at the crest of wind waves, the turbulent stress, and wave-induced stress. Studies on the dynamical structure of the high frequency wave spectrum are also needed.     

黑潮延伸体海域次中尺度过程的季节变化研究

近年来的现场观测和理论研究发现, 次中尺度现象广泛存在于上层海洋, 其产生与锋生作用及混合层斜压不稳定存在密切联系。本文利用高分辨率的数值模拟结果并结合动力学及能量诊断分析, 对黑潮延伸体海域次中尺度过程的季节变化进行了探讨。探讨结果表明, 黑潮延伸体海域次中尺度过程具有冬季最强, 春季和秋季次之, 夏季最弱的显著季节变化特征。基于冬、夏季次中尺度能量源的诊断可以看到, 这些季节变化特征主要与上层海洋的斜压不稳定和锋生作用有关。冬季, 黑潮延伸体海域的中尺度能量较弱, 但次中尺度过程在季节尺度上表现最为活跃, 这主要与混合层斜压不稳定的作用有关; 夏季, 黑潮延伸体海域的混合层较浅, 次中尺度过程较弱, 但中尺度涡旋活跃, 中尺度流场变形引起的锋生作用对夏季次中尺度现象的产生具有重要影响。在次中尺度能量的季节变化方面, 冬季次中尺度过程从中尺度过程汲取能量的速率远高于夏季, 这是冬季次中尺度过程比夏季更为活跃的主要原因。本文研究结果有助于加深对黑潮延伸体海域次中尺度过程季节性变化及其动力机制的理解。     

An accurate four-quadrant nonlinear dynamical model for marinethrusters: theory and experimental validation

This paper reports two specific improvements in the finite-dimensional nonlinear dynamical modeling of marine thrusters. Previously reported four-quadrant models have employed thin airfoil theory considering only axial fluid flow and using sinusoidal lift/drag curves. First, we present a thruster model incorporating the effects of rotational fluid velocity and inertia on thruster response. Second, we report a novel method for experimentally determining nonsinusoidal lift/drag curves. The model parameters are identified using experimental thruster data (force, torque, and fluid velocity). The models are evaluated by comparing experimental performance data with numerical model simulations. The data indicates that thruster models incorporating both reported enhancements provide superior accuracy in both transient and steady-state responses