Incorporating uncertainty into the estimation of biological reference points for a spiny lobster (Panulirus penicillatus) fishery

The status of a fishery is often defined as the probability of fishing mortality rate exceeding a perilous level for long‐term sustainability. Lobster stock assessments are often subject to large uncertainty in input data and high levels of natural variability in lobster life history processes, which calls for incorporating uncertainty associated with both indicator and management reference points in an evaluation of biological risk of overfishing. Using a Monte Carlo simulation approach, we evaluated the impacts of uncertainty in modelling on the determination of the status of the Taitung spiny lobster (Panulirus penicillatus) fishery (Taiwan), which has not been quantitatively determined despite its commercial importance. The commonly used biological reference points derived from the per recruit model (F _0.1 the fishing mortality rate where the slope of the curve of yield‐per‐recruit model is 10% of the maximum slope and F _4Q%, the fishing mortality rate that reduces the expected egg production for a cohort of female lobsters to 40% of that produced in the absence of a fishery of the egg‐per‐recruit model) were influenced by uncertainties associated with lobster life history and fishery parameters. A large uncertainty in the current fishing mortality rate (F _cnr) and estimates of biological reference points (F _BRPs) increased the uncertainty in determining the risk of overexploitation throughout the confidence levels of the stochastic decision‐making framework. This simulation study suggests that the target reference point of F _40% is less sensitive to the input parameters’ uncertainty than F _0.1 We suggest a further evaluation of other F‐based references points and development of biomass‐based reference points before final selection and implementation for the management of the Taitung lobster fishery.     

颗粒材料数值样本的坐标排序生成技术

颗粒材料离散颗粒模型的数值模拟结果与颗粒材料的数值样本密切相关,随着离散单元在颗粒材料数值模拟领域的广泛应用,颗粒材料的数值样本生成技术日益受到重视。基于RSA模型研究如何使随机生成的颗粒材料更密实,对均匀颗粒而言亦即如何在指定区域内生成更多的颗粒,讨论了4类修正方案,并建议了一种基于坐标排序的样本生成技术。研究表明,在传统的颗粒体随机生成技术基础上,通过对随机生成的x坐标序列或y坐标序列进行排序,可使生成的颗粒材料数值样本更密实。     

Application of an improved semi-Lagrangian procedure to fully-nonlinear simulation of sloshing in non-wall-sided tanks

The semi-Lagrangian procedure is widely used for updating the fully-nonlinear free surface in the time domain. However, this procedure is only available to cases when the body surface is vertical near the waterline. Present study introduces an improved semi-Lagrangian procedure which removes this ‘vertical-wall’ limitation. Coupling with the boundary element method, the improved semi-Lagrangian procedure is applied to the simulation of fully-nonlinear sloshing waves in non-wall-sided tanks. From the result comparison with the open source CFD software OpenFOAM, it is confirmed that this numerical scheme could guarantee a sufficient accuracy. Further series studies on 2D and 3D fully-nonlinear sloshing waves in wedged tanks are performed. Featured phenomena are observed which are distinct from those in wall-sided tanks.     

工程海域波浪特征分析方法比较

通过选取两个典型海外工程项目案例,分别采用波浪特征参数法和波浪谱法对工程海域波浪条件进行分析,并比较这两种方法在作业波浪条件、波浪特征认知程度上的差异,以及对海港水域平面布置产生的影响。对于波浪能量在频域及方向上存在多峰谱情况的海域,基于全谱的波浪特征参数系列数据有可能出现以下两种情况:1)主峰波向掩盖次峰波向;2)主峰波向错配次峰周期。对于上述两种情况,若仅采用波浪特征参数法对工程海域作业波浪条件进行分析,并指导海港水域平面布置,则极可能忽略或误判个别敏感波向、关键频域段的部分波浪能量,从而误导设计,为后续项目运营埋下隐患。对于以上情形,宜采用波浪谱法深入分析工程海域波浪能量在整个频域及各个方向上的详细分布及构成,并对关键波向、关键频域段的波浪能量进行量化,评估其对海港水域平面布置的影响,从而最大程度地避免由于波浪输入条件误判带来的风险。     

An Iteration Method for Computation of Flexible Fender Piles

—In the designing of the flexible fender pile,M-method is usually adopted.But M-method isusable only in the case that the value of the force(moment M and shear force Q)exerted on the pile abovethe soil is known.In practical cases,only the percussive energy of the ship can be found,while the exactvalue of the percussive force cannot be computed directly.In this paper,a method is introduced which canbe used to calculate the inner force of flexible fender piles by a computer in these cases.First,according tothe known docking dynamic energy of a ship,the absorbed energy can be found from the deformation ofthe pile and the rubber fender respectively,at the same time the value of the force acting on the fender pilesand the reaction of fenders can be calculated.Then the inner force and displacement of the fender piles canalso be calculated further by M-method.The whole computing process is realized by a computer programand the results obtained are reasonable.In this paper,the computation,the iteration method a     

A finite differences formulation for the linear and nonlinear dynamics of 2D catenary risers

A finite differences (FD) solution method is proposed for the numerical treatment of the dynamic equilibrium problem of 2D catenary risers. The method is based on the so-called Box approximation, which in the scope of the present contribution is applied to the complete nonlinear model as well as to the reduced linearized formulation. The application of the Box method transforms the original governing systems into convenient sets of algebraic equations, which in turn are solved efficiently by the relaxation method. Extensive numerical calculations are presented that describe the dynamic behaviour of the structure and evaluate the amplification in loading due to the dynamic components. The effect of the geometric nonlinearities is assessed through comparative calculations that concern both mathematical formulations examined in the present, i.e. the complete nonlinear, and the reduced linearized model. Special attention is paid to the heave excitations as they amplify significantly the magnitudes of the loading components.     

Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.     

A laboratory study of directional spectra with maximum likelihood method-Ⅰ Developing wind wave

measurements by a circular array consisting of 18 wave gauges in a large wave tank, directional spectra of wind-generated waves in deep water are systematically determined by using maximum likelihood method The investigations reveal that the angular spreading of the wave energy is consistent with cos^2s(θ/2) proposed by Longuet-Higgins et al.(1963, Ocean Wave Spectra, 111-136), if the bimodal distributions of wave energy are not taken into account. Bimodality occurring on higher frequency than peak frequency is too rare to affect our whole resalts. Surprisingly, a much broader directional spreading than that of the field, which is interpreted by the strongly nonlinear energy transfer because of the very young waves in laboratory, is found. The parameter s depends on frequency in the same way as observed by Mitsuyasu et al.(1975, Journal of Physical Oceanography, 5, 750-760) and Hasselmann et al.(1980, Journal of Physical Oceanography, 10, 1264-1280) in the field, and the relationship between the four nondimensional parameters s_m, f_o, b₁ and b₂, determining the directional width, and ū₁₀ (corresponding to the inverse of wave age) are given respectively. The observed distributions are found to agree well with the suggestion of Donelan et al.(1985, Philosophical Transaction of Royal Society of London, A315, 509-562) when applied to field waves.     

支持向量机在海洋随机资料分析中的应用

在随机海洋资料的分析中,需要从大量的实际数据中提取出有用的信息来建立或检验各种海洋理论模型,因此我们要正确运用适当的手段来分析处理各种实际数据。文中简要介绍了数据挖掘中一项新的技术支持向量机,着重介绍了最优化问题的求解算法,以及如何应用于海洋资料的分析处理。