A simulation tool for high data-rate acoustic communication in ashallow-water, time-varying channel

The paper discusses the development of a simulation tool to model high data-rate acoustic communication in shallow water. The simulation tool is able to generate synthetic time series of signals received at a transducer array after transmission across a shallow-water communication channel. The simulation tool is suitable for testing advanced signal processing techniques for message recovery. A channel model has been developed based on the physical aspects of the acoustic channel. Special emphasis has been given to fluctuations of the signal transmission caused by time-varying multipath effects. At shorter ranges, the temporal variations are dominated by acoustic scattering from the moving sea surface. Therefore, the channel model produces a coherence function which may be interpreted as a time-varying reflection coefficient for the surface scattered acoustical path. A static, range-independent ray model identifies the significant multipaths, and the surface path is modulated with the time-varying reflection coefficient. The advantages and limitations of the channel model are discussed and assumptions necessary to overcome the limitations are emphasised. Based on the assumptions, an algorithm has been developed and implemented to model how a binary message will be modulated when transmitted by a transducer, is distorted in the channel and finally is received by a transducer array     

Limiting heave response of a wave energy device by draft adjustment with upper surface immersion

Many studies have been published concerning the influence of the immersed shape (in still water) of a floating body on its response and power capture from ocean waves. With a few notable exceptions, much of this analysis has assumed small amplitude motion and linear models have been employed to predict response. The form of the upper surface of such a body has received little attention. Here, we show how the shape of the upper (top) surface of a floating body can be designed to ensure that the response amplitude of the body is within a specified value. This is of considerable importance to the survivability of wave energy devices. The approach used is to achieve a large increase of both natural period and hydrodynamic damping for only a small change of float mass. These two factors impose a hydrodynamic limit on the displacement which may be exploited to avoid the ‘end-stop’ problem often encountered in wave device design. To demonstrate the change of response, experimental measurements are presented of the response of an axisymmetric float with rounded base and conical upper surface with rounded perimeter due to a range of regular, irregular and focused wave conditions. Power extraction is not considered since the mechanically undamped response represents the worst case. In contrast to a simple, straight-sided axisymmetric float, a smaller change of mass is required to satisfy a particular response amplitude limit. Although a significant reduction is not expected, hydrodynamic damping may reduce with increasing physical scale, and this remains to be quantified.     

Functional responses for zooplankton feeding on multiple resources: a review of assumptions and biological dynamics

Modelers often need to quantify the rates at which zooplankton consume a variety of species, size classes and trophic types. Implicit in the equations used to describe the multiple resource functional response (i.e. how nutritional intake varies with resource densities) are assumptions that are not often stated, let alone tested. This is problematic because models are sensitive to the details of these formulations. Here, we enable modelers to make more informed decisions by providing them with a new framework for considering zooplankton feeding on multiple resources. We define a new classification of multiple resource responses that is based on preference, selection and switching, and we develop a set of mathematical diagnostics that elucidate model assumptions. We use these tools to evaluate the assumptions and biological dynamics inherent in published multiple resource responses. These models are shown to simulate different resource preferences, implied single resource responses, changes in intake with changing resource densities, nutritional benefits of generalism, and nutritional costs of selection. Certain formulations are further shown to exhibit anomalous dynamics such as negative switching and sub-optimal feeding. Such varied responses can have vastly different ecological consequences for both zooplankton and their resources; inappropriate choices may incorrectly quantify biologically-mediated fluxes and predict spurious dynamics. We discuss how our classes and diagnostics can help constrain parameters, interpret behaviors, and identify limitations to a formulation's applicability for both regional (e.g. High-Nitrate-Low-Chlorophyll regions comprising large areas of the Pacific) and large-scale applications (e.g. global biogeochemical or climate change models). Strategies for assessing uncertainty and for using the mathematics to guide future experimental investigations are also discussed.     

On the long-term response of marine structures

This paper addresses some important issues related to the estimation of long-term extreme responses of marine structures. Several convolution models to establish the long-term distribution of a marine structure response parameter are available in the literature. These methods are typically based either on all short-term peaks, all extreme short-term peaks or all short-term upcrossing rates. The main assumptions and simplifications of the five models most usually found in the literature are discussed in this paper. A linear single-degree-of-freedom (SDOF) system along with a bi-lognormal probability model for significant wave heights and zero-crossing wave periods have been used for numerical tests. An improved approach to efficiently evaluate the long-term convolution integrals is also proposed in this paper. It is shown that a combination of the Inverse First Order Reliability Method (IFORM) and an Importance Sampling Monte Carlo Simulation (ISMCS) approach can be used to obtain a very good result for the exact solution of long-term integrals.     

The threshold feeding response of microzooplankton within Pacific high-nitrate low-chlorophyll ecosystem models under steady and variable iron input

The equatorial Pacific is an HNLC (High-Nitrate Low-Chlorophyll) region. Modeling and in-situ process studies have confirmed the importance of microzooplankton grazing in this ecosystem. Unfortunately, both the parameters and functions representing microzooplankton grazing within current ecosystem models are poorly constrained. We used a simple 4-component food web model to test the assumption that a lower grazing threshold, which is common in many models, is necessary to achieve the HNLC condition. Without the grazing threshold, the model did not reproduce the HNLC condition. However, by raising the half-saturation constant within the microzooplankton functional response with no threshold, it was possible to reproduce the critical dynamics of the HNLC condition under both steady and moderate seasonal variability in nutrient input. It was also possible to reproduce the HNLC system using a sigmoidal functional response for the microzooplankton, with results somewhere between the other two forms of the model, although this version had the highest sensitivity to changes in its parameters. The three models predicted similar phytoplankton biomass and primary productivity under steady nutrient input, but diverge in these metrics as the amplitude of nutrient input variability increases. These three functional responses also imply certain important differences in the microzooplankton community. Whereas the threshold model had the least sensitivity to parameter choice, the high half-saturation constant, no-threshold model may actually be a better approximation when modeling a community of grazers. Ecosystem models that predict carbon production and export in HNLC regions can be very sensitive to assumptions concerning microzooplankton grazing; future studies need to concentrate on the functional responses of microzooplankton before these models can be used for predicting fluxes in times or regions where forcing is beyond that used to constrain the original model.     

Incorporating irregular nonlinear waves in coupled simulation and reliability studies of offshore wind turbines

Design of an offshore wind turbine requires estimation of loads on its rotor, tower and supporting structure. These loads are obtained by time-domain simulations of the coupled aero-servo-hydro-elastic model of the wind turbine. Accuracy of predicted loads depends on assumptions made in the simulation models employed, both for the turbine and for the input wind and wave conditions. Currently, waves are simulated using a linear irregular wave theory that is not appropriate for nonlinear waves, which are even more pronounced in shallow water depths where wind farms are typically sited. The present study investigates the use of irregular nonlinear (second-order) waves for estimating loads on the support structure (monopile) of an offshore wind turbine. We present the theory for the irregular nonlinear model and incorporate it in the commonly used wind turbine simulation software, FAST, which had been developed by National Renewable Energy Laboratory (NREL), but which had the modeling capability only for irregular linear waves. We use an efficient algorithm for computation of nonlinear wave elevation and kinematics, so that a large number of time-domain simulations, which are required for prediction of long-term loads using statistical extrapolation, can easily be performed. To illustrate the influence of the alternative wave models, we compute loads at the base of the monopile of the NREL 5MW baseline wind turbine model using linear and nonlinear irregular wave models. We show that for a given environmental condition (i.e., the mean wind speed and the significant wave height), extreme loads are larger when computed using the nonlinear wave model. We finally compute long-term loads, which are required for a design load case according to the International Electrotechnical Commission guidelines, using the inverse first-order reliability method. We discuss a convergence criteria that may be used to predict accurate 20-year loads and discuss wind versus wave dominance in the load prediction. We show that 20-year long-term loads can be significantly higher when the nonlinear wave model is used.     

Modelling of multi-bodies in close proximity under water waves—Fluid forces on floating bodies

This work presents two-dimensional numerical results of the dependence of wave forces of multiple floating bodies in close proximity on the incident wave frequency, gap width, body draft, body breadth and body number based on both viscous fluid and potential flow models. The numerical models were validated by the available experimental data of fluid oscillation in narrow gaps. Numerical investigations show that the large amplitude responses of horizontal and vertical wave forces appear around the fluid resonant frequencies. The convectional potential flow model is observed to un-physically overestimate the magnitudes of wave forces as the fluid resonance takes place. By introducing artificial damping term with appropriate damping coefficients μ∈[0.4, 0.5], the potential flow model may work as well as the viscous fluid model, which agree with the damping coefficients used in our previous work for the predication of wave height under gap resonance. In addition, the numerical results of viscous fluid model suggest that the horizontal wave force is highly dependent on the water level difference between the opposite sides of an individual body and the overall horizontal wave force on the floating system is generally smaller than the summation of wave force on each body.     

一个区域海洋——海冰——大气耦合模式中的北极气候偏差：年结果检验

The Coupling of three model components, WRF/PCE （polar climate extension version of weather research and forecasting model （WRF））, ROMS （regional ocean modeling system）, and CICE （community ice code）, has been implemented, and the regional atmosphere-ocean-sea ice coupled model named WRF/PCE- ROMS-CICE has been validated against ERA-interim reanalysis data sets for 1989. To better understand the reasons that generate model biases, the WRF/PCE-ROMS-CICE results were compared with those of its components, the WRF/PCE and the ROMS-CICE. There are cold biases in surface air temperature （SAT） over the Arctic Ocean, which contribute to the sea ice concentration （SIC） and sea surface temperature （SST） biases in the results of the WRF/PCE-ROMS-CICE. The cold SAT biases also appear in results of the atmo- spheric component with a mild temperature in winter and similar temperature in summer. Compared to results from the WRF/PCE, due to influences of different distributions of the SIC and the SST and inclusion of interactions of air-sea-sea ice in the WRF/PCE-ROMS-CICE, the simulated SAT has new features. These influences also lead to apparent differences at higher levels of the atmosphere, which can be thought as responses to biases in the SST and sea ice extent. There are similar atmospheric responses in feature of distribution to sea ice biases at 700 and 500 hPa, and the strength of responses weakens when the pressure decreases in January. The atmospheric responses in July reach up to 200 hPa. There are surplus sea ice ex- tents in the Greenland Sea, the Barents Sea, the Davis Strait and the Chukchi Sea in winter and in the Beau- fort Sea, the Chukchi Sea, the East Siberian Sea and the Laptev Sea in summer in the ROMS-CICE. These differences in the SIC distribution can all be explained by those in the SST distributions. These features in the simulated SST and SIC from ROMS-CICE also appear in the WRF/PCE-ROMS-CICE. It is shown that the performance of the WRF/PCE-ROMS-CICE is determined to a l     

Ekman模型时间变化风应力系数的伴随参数估计

本文基于时间分布参数设置,利用伴随同化方法,反演了Ekman模型中随时间变化的风应力拖曳系数,并在孪生实验和实际实验中对该方法进行了验证。在孪生实验中,研究了参数反演结果对不同影响因素的响应,包括:风速分布、风应力系数分布、风应力系数初始猜测值、风应力系数独立变量个数、观测数据误差和观测的深度。孪生实验结果验证了伴随同化方法反演Ekman模型中时变风应力系数的有效性,具体包括如下五个方面结论:1)不同风速分布下均能成功反演出不同风应力拖曳系数分布; 2)反演结果对初始猜测值较为敏感,风应力系数初始猜测值越接近给定值,反演结果越好;3)风应力系数独立点个数的选取会显著影响反演结果,合理的选择有利于提高反演效率及减小观测数据误差;4)观测误差能够影响反演结果,观测数据误差在20%以下时能取得合理的反演结果; 5)反演结果对观测数据的表层和次表层流速更为敏感,这是由Ekman流的物理性质决定的。实际实验,利用百慕大锚系试验平台的风速和流速数据,去除周期性潮流和地转流成分后得到Ekman流成分,并作为观测输入到该同化模型,反演出了适用于该区域和该时段的随时间变化的风应力系数。通过比较模拟流速和观测流速,证明利用伴随同化方法能从实测数据中反演出合理的时变风应力系数,对于海洋模型风应力系数的确定是一项有益的尝试。     

CAR模型的阶与时滞的强一致估计

研究随机线性反馈控制系统的结构辨识问题。在已知时滞的下界和模型阶的上界的假定下，通过使修改的Ｂａｙｅｓｉａｎ信息准则最小化，推导出由多输入多输出ＣＡＮ模型描述的系统的未知阶与时滞的估计算法，证明了算法是强一致收敛的，且能在有限步内达到其模型结构参数的真值。讨论了当模型的参数矩阵不满秩时减弱条件Ｈ’ｓ的强一致估计算法。     

Identification-based simplified model of large container ships using support vector machines and artificial bee colony algorithm

The 6 degrees of freedom (DOF) model with a high degree of complexity for capturing ship dynamics is generally able to track the nonlinear and coupling dynamics of ships. However, the 6 DOF model makes challenges in estimating model coefficients and designing the model-based control. Therefore, simplified ship dynamic models within allowed accuracy are essential. This paper simplified the 6 DOF nonlinear dynamic model of ships into two decoupled models including the speed model and the steering model through reasonable assumptions. Those models were tested through maneuvering simulations of a container ship with a 4 DOF dynamic model. Support vector machines (SVM) optimized by the artificial bee colony algorithm (ABC) was used to identify parameters of speed and steering models by analyzing the rudder angle, propeller shaft speed, surge and sway velocities, and yaw rate from simulated data extracted from a series of maneuvers made by the container ship. Comparisons with the first order linear and nonlinear Nomoto models show that the simplified nonlinear steering model can capture more complicated dynamics and performs better. Additionally, comparisons among three different parameter identification methods demonstrate similar identification results but the different performance involving the applicability and effectiveness. SVM optimized by ABC is relatively convenient and effective for parameter identification of ship simplified dynamic models.     

Boundary element analysis of liquid sloshing characteristics in axisymmetric tanks with various porous baffles

This paper aims to investigate the effects of the porous baffles on the suppression of sloshing for the tanks with axisymmetric geometries under lateral excitation. Based on the assumptions of inviscid, irrotational, incompressible liquid and small amplitude sloshing, an axisymmetric boundary element method (BEM) for 3D Laplace equation is derived by using the Green's theorem together with the weighted residual method. And a zoning method is employed to model fluid domain in the tanks with complex porous baffles. Meanwhile, the porous baffles are treated motioning together with the tanks, and the velocity across the porous baffle is assumed to be linearly proportional to the pressure gradient between each side of the porous baffle. And the mechanism of suppressing the sloshing response is mainly the energy dissipation of the fluid passing through the porous baffle. Moreover, the linear free surface boundary conditions are also used to solve the governing equations. Compared with other numerical methods, the most prominent advantage of the BEM in solving axisymmetric potential problem is that only the boundaries of half the cross-section instead of the entire problem domain should be discretized, which can cut down large amount of memory and time costs. The present method is verified by comparing the numerical results with the existing literatures, and excellent agreements are obtained. Meanwhile, the proposed models are applied to investigate the effects of the porous baffles on sloshing response in circular cylindrical, annular cylindrical and conical tanks. The effects of the porous baffle length, porous-effect parameter, installation angle and baffle height on the sloshing force, natural frequency and surface elevation are studied. Additionally, some typical sloshing pressure distributions, velocity potential contours and velocity fields are plotted. The results show that swirls at the tips of the baffles can be observed in many cases, and the top-mounted porous baffle makes more significant suppression effects on sloshing response than that of bottom-mounted porous baffle, while increasing the number of ring porous baffles can achieve better restraint effects on sloshing response. And increasing the baffle length of the horizontal wall-mounted ring porous baffle can significantly decrease the sloshing frequencies, as well as the first non-dimensional natural frequency decreases with decrease in porous-effect parameter of the coaxial porous baffle. In addition, remarkable effects on sloshing can be obtained when reasonable designed by selecting the optimal porous-effect parameter, installation angle and baffle height. And this paper can be a useful guide for the seismic design and analysis of many actual liquid storage tanks (such as the Advanced Passive PWR, large water cooling tower, etc.).     

Ecosystem modelling in the southern Benguela: comparisons of Atlantis,Ecopath with Ecosim,and OSMOSE under fishing scenarios

Ecosystem-based management of marine fisheries requires the use of simulation modelling to investigate the system-level impact of candidate fisheries management strategies. However, testing of fundamental assumptions such as system structure or process formulations is rarely done. In this study, we compare the output of three different ecosystem models (Atlantis, Ecopath with Ecosim, and OSMOSE) applied to the same ecosystem (the southern Benguela), to explore which ecosystem effects of fishing are most sensitive to model uncertainty. We subjected the models to two contrasting fishing pressure scenarios, applying high fishing pressure to either small pelagic fish or to adult hake. We compared the resulting model behaviour at a system level, and also at the level of model groups. We analysed the outputs in terms of various commonly used ecosystem indicators, and found some similarities in the overall behaviour of the models, despite major differences in model formulation and assumptions. Direction of change in system-level indicators was consistent for all models under the hake pressure scenario, although discrepancies emerged under the small-pelagic-fish scenario. Studying biomass response of individual model groups was key to understanding more integrated system-level metrics. All three models are based on existing knowledge of the system, and the convergence of model results increases confidence in the robustness of the model outputs. Points of divergence in the model results suggest important areas of future study. The use of feeding guilds to provide indicators for fish species at an aggregated level was explored, and proved to be an interesting alternative to aggregation by trophic level.     

误差方差二次型估计模型的转化

研究一般的回归模型中误差方差的二次型估计的容许性，研究方法是模型的整体转化和局部转化，结果有：（１）二次约束下的线性模型等价于相应的无约束的线性模型。（２）线性（齐次或非齐次）等式约束下的线性模型等价于某个无约束的线性模型。（３）单个非齐次不等式约束下的线性模型等价于某个无约束的线性模型。（４）通过例子证明了多个线性不等式约束的线性模型不能等价于某个无约束的线性模型。（５）某类非齐次二次型估计的容许性等价于相应的齐次二次型估计的容许性     

非线性随机系统的参数辨识及动态预报

鉴于含有加性噪声的指数模型描述了一类重要的非线性随机系统。本文给出这类系统的参数递推辨识算法，克服了迭代算法不能在线运行、需反复矩阵求逆的不足。当系统时变时，还采用了虚拟噪声技术来补偿因参数时变引起的建模误差，从而改善动态预报器的性能。应用这种方法，对油田采油井、注水井的套管损坏情况进行了多步动态预报，其精度令人满意     

Wave measurements by pressure transducers using artificial neural networks

Underwater ultrasonic acoustic transducers are frequently used in ocean wave measurements as they measure surface level using acoustic waves. However, their effectiveness can be severely affected in rough sea conditions, when bubbles generated by breaking waves interfere with their acoustic signals. When the seas are rough, one therefore often has to rely on a pressure transducer, which is generally used as a back-up for the acoustic wave gauge. A pressure transfer function is then used to obtain the surface wave information. Alternatively, the present study employed an artificial neural network to convert the pressure signal into significant wave height, significant wave period, maximum wave height, and spectral peakedness parameter using data obtained from various water depths. The results showed that, for water depths greater than 20 m, the wave parameters obtained from the artificial neural network were significantly closer to those obtained by the acoustic measurements than those obtained by using a linear pressure transfer function. Moreover, for a given water depth, the wave heights estimated by the network model from pressure data were not as good as those estimated by linear wave theory for large wave heights (above a 4 m significant wave height in this study). This can be improved if the training data set has more records with large wave heights.     

Limitations of the streamtube model

A common feature of one-dimensional streamtube models is that the pressure gradient includes only the effect of the density difference between the plume and the ambient water on a sloping bottom. Variations of the streamtube height and density along the flow path are thus neglected. In the present paper a streamtube model is used to highlight the assumptions underlying the simplification and to obtain the limitations this prerequisite presents. The results indicate that the parameter range for which the simplification is applicable is largest when the initial outflow is parallel to the isobaths. More specific limits are given for three major outflows — the Mediterranean, the Faroe Bank Channel and the Denmark Strait outflow.     

Combined use of dimensional analysis and modern experimental design methodologies in hydrodynamics experiments

In this paper, a combined use of dimensional analysis (DA) and modern statistical design of experiment (DOE) methodologies is proposed for a hydrodynamics experiment where there are a large number of variables. While DA is well-known, DOE is still unfamiliar to most ocean engineers although it has been shown to be useful in many engineering and non-engineering applications. To introduce and illustrate the method, a study concerning the thrust of a propeller is considered. Fourteen variables are involved in the problem and after dimensional analysis this reduces to 11 dimensionless parameters. Then, a two-level fractional factorial design was used to screen out parameters that do not significantly contribute to explaining the dependent dimensionless parameter. With the remaining five statistically significant dimensionless parameters, various response surface methodologies (RSM) were used to obtain a functional relationship between the dependent dimensionless thrust coefficient, and the five dimensionless parameters. The final model was found to be of reasonable accuracy when tested against results not used to develop the model. The methodologies presented in the paper can be similarly applied to systems with a large number of control variables to systematically derive approximate mathematical models to predict the responses of the system economically and accurately.     

我国主要海洋产业发展规模预测分析

通过收集整理1952—2002年50a的海洋经济统计数据,结合计量经济学手段进行时间序列分析。依样本空间和分析对象的不同,运用趋势外推法分别建立了线性模型、抛物线模型、生长曲线模型3种预测外延模型结构。并对预测外延模型结构进行了数据拟合与相关参数分析,总结出了适合预测海洋经济发展的线性模型和抛物线模型。根据海洋渔业、海洋油气业、海滨砂矿业、海洋盐业、海洋船舶工业、海洋交通运输业、滨海国际旅游业等7个主要海洋产业各自发展特点,用预测模型分别进行拟合,建立适合各海洋产业的预测模型。最终对2010年全国6个主要海洋产业、分产业总产值进行了预测。     

线性时变大系统的分散最优控制与镇定

本文提出一种关于线性时变大系统最优控制的二次型性能指标分解方法。据此研究线性时变大系统局部分散最优控制的逆问题．得到大系统关于某二次型性能指标存在最优状态反馈控制律的条件，同时它也是大系统局部分散镇定的一个新的充分条件．