水下滑翔器主体外形优化设计

对水下滑翔器主体外形的优化设计问题进行研究.首先,在考虑内部机械结构限制的条件下,建立了主体外形优化数学模型,并采用CFX计算了若干主体形状的绕流阻力.在此基础上建立了映射主体外形的尺寸参数与绕流阻力的BP神经网络.然后将建立的神经网络作为优化设计问题的目标函数,采用坐标轮换法对滑翔器主体外形进行了优化设计,确定了尺寸参数的最优解.     

克里金内插法实现坐标转换的应用研究

以陕西地区2000国家大地坐标系(CGCS2000)和旧1954年北京大地坐标系转换为例,探讨了利用普通克里金内插法进行坐标转换的实现方法,比较了采用不同变异函数和空间数据窗口的转换精度。试算结果表明:大地控制网的变形与方向无关;采用适当的变异函数模型比采用简单的线性变异函数模型精度提高约40%且空间数据窗口可大幅度缩小;基于普通克里金内插方法实现陕西地区坐标转换精度可达±0.092m。     

Research on Design Method of the Full Form Ship with Minimum Thrust Deduction Factor

In the preliminary design stage of the full form ships, in order to obtain a hull form with low resistance and maximum propulsion efficiency, an optimization design program for a full form ship with the minimum thrust deduction factor has been developed, which combined the potential flow theory and boundary layer theory with the optimization technique. In the optimization process, the Sequential Unconstrained Minimization Technique(SUMT) interior point method of Nonlinear Programming(NLP) was proposed with the minimum thrust deduction factor as the objective function. An appropriate displacement is a basic constraint condition, and the boundary layer separation is an additional one. The parameters of the hull form modification function are used as design variables. At last, the numerical optimization example for lines of after-body of 50000 DWT product oil tanker was provided, which indicated that the propulsion efficiency was improved distinctly by this optimal design method.     

Objective array design for three-dimensional temperature and salinity observation: Application to the South China Sea

In this study, a moored array optimization tool (MAOT) was developed and applied to the South China Sea (SCS) with a focus on three-dimensional temperature and salinity observations. Application of the MAOT involves two steps: (1) deriving a set of optimal arrays that are independent of each other for different variables at different depths based on an empirical orthogonal function method, and (2) consolidating these arrays using a K-center clustering algorithm. Compared with the assumed initial array consisting of 17 mooring sites located on a 3°×3° horizontal grid, the consolidated array improved the observing ability for three-dimensional temperature and salinity in the SCS with optimization efficiencies of 19.03% and 21.38%, respectively. Experiments with an increased number of moored sites showed that the most cost-effective option is a total of 20 moorings, improving the observing ability with optimization efficiencies up to 26.54% for temperature and 27.25% for salinity. The design of an objective array relies on the ocean phenomenon of interest and its spatial and temporal scales. In this study, we focus on basin-scale variations in temperature and salinity in the SCS, and thus our consolidated array may not well resolve mesoscale processes. The MAOT can be extended to include other variables and multi-scale variability and can be applied to other regions.     

Intercomparison of three South China Sea circulation models

1IntroductionTheSouthChinaSeaisthelargesttropicalmarginaldeepsealocatingbetweenthewesternPacificOceanandtheeasternIndianOcean.AsapartofAsia-Australiamaritimecontinent,monsoonisaprimaryfactorforcingtheSouthChinaSeaCurrent(SCSC)variation.Drivenbynortheasterlymonsooninwinterandsouth-westerlymonsooninsummer,respectively,theSCSCbehavesacyclonicgyreandananticy-clonicgyre,correspondingly(Wyrtki,1961;Xuetal.,1982).Owingtotheshortageandexpen-sivenessofdirectobservationsintheSCS,fur-therunder…     

Parametric Geometric Model and Hydrodynamic Shape Optimization of A Flying-Wing Structure Underwater Glider

Combining high precision numerical analysis methods with optimization algorithms to make a systematic exploration of a design space has become an important topic in the modern design methods. During the design process of an underwater glider''s flying-wing structure, a surrogate model is introduced to decrease the computation time for a high precision analysis. By these means, the contradiction between precision and efficiency is solved effectively. Based on the parametric geometry modeling, mesh generation and computational fluid dynamics analysis, a surrogate model is constructed by adopting the design of experiment (DOE) theory to solve the multi-objects design optimization problem of the underwater glider. The procedure of a surrogate model construction is presented, and the Gaussian kernel function is specifically discussed. The Particle Swarm Optimization (PSO) algorithm is applied to hydrodynamic design optimization. The hydrodynamic performance of the optimized flying-wing structure underwater glider increases by 9.1%.     

基于拓扑优化与多目标优化的集装箱相似畸变模型设计

使用相似原理构建缩尺模型是探究大型结构性能的最基本和最优途径。基于PATRAN软件的PCL语言与ISIGHT软件,提出一种联合静动力学、拓扑优化与多目标优化的集装箱相似畸变模型构建方法。利用拓扑优化的SIMP法(Solid Isotropic Material with Penalization)得到缩尺集装箱墙壁结构的最优材料分布,既保证了结构性能又降低了缩尺模型的重量;基于全局多目标算法建立优化模型,采用有限元方法对集装箱性能进行分析,以结构刚度和一阶模态频率为目标函数进行寻优,得到相似模型的最优设计参数。最终利用有限元分析验证该模型的合理性。结果表明:该方法具有一定的通用性,能为其它类似结构的相似模型构建提供参考。     

A hybrid coordinate ocean model for shelf sea simulation

The general circulation in the North Sea and Skagerrak is simulated using the hybrid coordinate ocean model (HYCOM). Although HYCOM was originally developed for simulations of the open ocean, it has a design which should make it applicable also for coastal and shallow shelf seas. Thus, the objective of this study has been to examine the skills of the present version of HYCOM in a coastal shelf application, and to identify the areas where HYCOM needs to be further developed. To demonstrate the capability of the vertical coordinate in HYCOM, three experiments with different configurations of the vertical coordinate were carried out. In general, the results from these experiments compares quite well with in situ and satellite data, and the water masses and the general circulation in the North Sea and Skagerrak is reproduced in the simulations. Differences between the three experiments are small compared to other errors, which are related to a combined effect of model setup and properties of the vertical mixing scheme. Hence, it is difficult to quantify which vertical coordinate configuration works best for the coastal region. It is concluded that HYCOM can be used for simulations of coastal and shelf seas, and further suggestions for improving the model results are given. Since HYCOM also works well in open ocean and basin scale simulations, it may allow for a realistic modelling of the transition region between the open ocean and coastal shelf seas.     

多波束回波强度位置归算及图像化表示

多波束在测深的同时也记录回波强度数据。在分析多波束回波强度记录方式的基础上,主要针对波束内回波强度的定位和数据优化等问题进行了探讨,确定了回波强度位置归算方法。最后,以Simrad EM1002多波束实测的回波强度数据为例,进一步阐述了在地理坐标框架下回波强度的图像化表示方法。     

基于改进量子粒子群算法的AUV路径规划研究

针对海洋环境下自主水下机器人（AUV）的路径规划问题,提出了一种基于框架四叉树的改进量子粒子群算法（QPSO）,首先使用框架四叉树的方法对障碍物建模,该方法提高了建模的精度且对后续算法的效率也有极大的改进,之后设计改进的量子粒子群算法,并且结合水下环境的特殊性设计适应度函数,综合考虑航线路径长度、偏转角度以及海流影响,使得算法可以在水下环境中寻得能耗最短的解路径。最后通过仿真试验验证,相比于传统的栅格法和粒子群算法,改进量子粒子群算法的运算时间更短,收敛速度更快,其独特的适应度函数可以使AUV能更好适应水下多变的环境,且能利用海流设计能耗更小的路径,具有很大的实用价值。     

空间直角坐标系统转换的抗差算法研究

为了解决空间直角坐标系间相互转换时,由于公共点的坐标精度差直接导致坐标转换精度差的问题.探讨利用抗差估计理论解决公共点坐标粗差的影响,结合具体的坐标转换算例,分别采用Tukey、IGG1及IGG3三种权函数进行迭代计算,完成坐标转换.结果表明,抗差估计原理用于空间直角坐标系间的相互转换,可以降低公共点粗差的影响,提高坐...     

Parametric vertical coordinate formulation for multiscale, Boussinesq, and non-Boussinesq ocean modeling

Two physical parameters are introduced into the basic ocean equations to generalize numerical ocean models for various vertical coordinate systems and their hybrid features. The two parameters are formulated by combining three techniques: the arbitrary vertical coordinate system of Kasahara [Kasahara, A., 1974. Various vertical coordinate systems used for numerical weather prediction. Mon. Weather Rev. 102, 509–522], the Jacobian pressure gradient formulation of Song [Song, Y.T., 1998. A general pressure gradient formation for ocean models. Part I: Scheme design and diagnostic analysis. Mon. Weather Rev. 126 (12), 3213–3230], and a newly introduced parametric function that permits both Boussinesq (volume-conserving) and non-Boussinesq (mass-conserving) conditions. Based on this new formulation, a generalized modeling approach is proposed. Several representative oceanographic problems with different scales and characteristics––coastal canyon, seamount topography, non-Boussinesq Pacific Ocean with nested eastern Tropics, and a global ocean model––have been used to demonstrate the model’s capabilities for multiscale applications. The inclusion of non-Boussinesq physics in the topography-following ocean model does not incur computational expense, but more faithfully represents satellite-observed ocean-bottom-pressure data. Such a generalized modeling approach is expected to benefit oceanographers in solving multiscale ocean-related problems by using various coordinate systems on the same numerical platform.     

Hull form reliability-based robust design optimization combining polynomial chaos expansion and maximum entropy method

In the ship design optimization process, the neglect of the unavoidable uncertainty of parameters in the actual navigation and experimental observations, may lead to a bad result with some hidden dangers in practical applications. Considering the influence of the uncertainty, a new and effective hull form reliability-based robust design optimization (RBRDO) framework, including the hull form modification module and RBRDO module, has been developed and tested in the present work. Radial basis function method is utilized as the parametric hull surface modification technique to generate a series of smooth hull forms while combining polynomial chaos expansions (PCE) method with maximum entropy method (MEM) to conduct the uncertainty analysis for the prediction of the mean and the standard deviation of the objective and the failure probability of constraints. To verify the validity of the method, hull form design optimization of the bow of KCS model is implemented under the influence of the uncertainty. Numerical results indicate that the proposed RBRDO framework is effective compared with traditional Monte Carlo method. Meanwhile, compared with traditional DO case, RBRDO case has higher adaptability to the environmental uncertainty with the lower failure probability, which ensure the robustness and reliability of the optimal hull form.     

Stationary Rossby waves and shocks on the Sverdrup coordinate

Gyre-scale frontal structures, e.g., the Subtropical Front, have been well documented in the oceans. Although the generation mechanism of such a front remains unclear, it may be ascribed to the steepening of nonlinear planetary waves as discussed by Dewar (1992). To discuss the stationary wave characteristics and their shock formation in a two-layer wind-driven gyre, the present paper introduces a new coordinate, referred to as theSverdrup coordinate here, in which the Sverdrup function is used instead of the longitudinal coordinate, and especially, investigates the possibility of the frontgenesis caused by the Rossby waves emanating from the western boundary region. On theSverdrup coordinate, since the advection by the Sverdrup flow is in the direction normal to that of the Rossby wave propagation, the system becomes much simpler than that on the (x,y)-coordinate, and the solution to this coordinate does not explicitly depend on the distribution of the Ekman pumping, i.e., we can treat cases, in which the Ekman pumping is a function of bothx andy, in a similar way to the case with zonally uniform Ekman pumping.     

Formulation, implementation and examination of vertical coordinate choices in the Global Navy Coastal Ocean Model (NCOM)

A 1/8° global version of the Navy Coastal Ocean Model (NCOM) is described with details of its formulation, implementation, and configuration of the vertical coordinate. NCOM is a baroclinic, hydrostatic, Boussinesq, free-surface ocean model that allows its vertical coordinate to consist of σ coordinates for the upper layers and z-levels below a user-specified depth. This flexibility allows implementation of a hybrid σ–z coordinate system that is expected to mitigate some of the weaknesses that can be associated with either pure coordinate option. For the global NCOM application, the σ–z coordinate is used to allow terrain-following σ coordinates in the upper ocean, providing better resolution and topographic fidelity in shelf regions where flow is most sensitive to its representation. Including z coordinates for deeper regions efficiently maintains high near-surface vertical resolution in the open ocean. Investigation into the impact of the selected coordinate system focuses on differences between atmospherically-forced free-running (no assimilation) global solutions using σ–z and pure z coordinates. Comparisons with independent temperature observations indicate that global NCOM using the σ–z coordinate has improved skill relative to its z coordinate implementation. Among other metrics, we show that in comparison with time series of surface temperature from National Oceanic Data Center (NODC) buoys, mostly located in coastal regions, root mean squared differences (RMSD) improved for 63% and correlation improved for 71% of the stations when σ–z coordinates were used instead of pure z. For the exclusively open-ocean Tropical Atmosphere-Ocean (TAO) buoys, differences between the simulations were small, with the σ–z showing smaller RMSD for 45% of the stations and higher correlation for 65% of the stations. Additional comparisons using temperature profile observations further confirm a tendency for improved performance using the hybrid σ–z coordinates.     

Free fall water entry of a wedge tank into calm water in three degrees of freedom

The hydrodynamic problem of a two dimensional wedge tank filled with liquid entering a calm water surface is analysed based on the incompressible velocity potential theory. The motion effect of inner liquid on the entry process is investigated through comparison with the result containing equivalent solid mass or the liquid being frozen. The problem is solved through the boundary element method in the time domain. Two separated computational regions are constructed. One is the inner domain for the internal liquid, and the other is the outer open domain for the open water. The former is solved in the physical coordinate system, and the latter is solved in a stretched coordinate system. The solutions of two separated domains are connected through the motion of the body. The auxiliary function method is extended to decouple the nonlinear mutual dependence between fluid loads from two separated domains and the body motion. Detailed results for wedge motion, external impact pressure and free surface, and for internal pressure, free surface deformation and liquid motion are provided. Through comparison with the results of a wedge tank with frozen ice, in-depth discussion on the effect of the inner liquid is provided.     

Research on Bulbous Bow Optimization Based on the Improved PSO Algorithm

In order to reduce the total resistance of a hull, an optimization framework for the bulbous bow optimization was presented. The total resistance in calm water was selected as the objective function, and the overset mesh technique was used for mesh generation. RANS method was used to calculate the total resistance of the hull. In order to improve the efficiency and smoothness of the geometric reconstruction, the arbitrary shape deformation (ASD) technique was introduced to change the shape of the bulbous bow. To improve the global search ability of the particle swarm optimization (PSO) algorithm, an improved particle swarm optimization (IPSO) algorithm was proposed to set up the optimization model. After a series of optimization analyses, the optimal hull form was found. It can be concluded that the simulation based design framework built in this paper is a promising method for bulbous bow optimization.     

A New Maximum Entropy Probability Function for the Surface Elevation of Nonlinear Sea Waves

Based on the maximum entropy principle a new probability density function （PDF） f（x） for the surface elevation of nonlinear sea waves, X, is derived through performing a coordinate transform of X and solving a variation problem subject to three constraint conditions of f（ x ）. Compared with the maximum entropy PDFs presented previously, the new PDF has the following merits： （1） it has four parameters to be determined and hence can give more refined fit to observed data and has wider suitability for nonlinear waves in different conditions; （2） these parameters are expressed in terms of distribution moments of X in a relatively simple form and hence are easy to be determined from observed data; （3） the PDF is free of the restriction of weak nonlinearity and possible to be used for sea waves in complicated conditions, such as those in shallow waters with complicated topography; and （4） the PDF is simple in form and hence convenient for theoretical and practical uses. l.aboratory wind-wave experiments have been conducted to test the competence of the new PDF for the surface elevation of nonlinear waves. The experimental results manifest that the new PDF gives somewhat better fit to the laboratory wind-wave data than the well-known Gram-Charlier PDF and beta PDF.     

从地方坐标系到2000国家大地坐标系的转换方法

探讨了我国原有地方坐标系与CGCS2000坐标系的定义差别,在我国相关学者研究成果的基础上,设计一种基于布尔莎七参数转换模型的坐标系统转换方法,可以实现除高等级大地控制点之外的各种地理信息数据从原有地方坐标到CGCS2000坐标的快速转换。     

Reliability-based design optimization of an FPSO riser support using moving least squares response surface meta-models

The paper deals with the reliability-based design optimization (RBDO) of a riser support installed on a floating production storage and offloading (FPSO) unit under operation, extreme, damaged, and one line failure cases and installation loading conditions. The optimization problem is formulated such that probabilistic thickness variables described with random characteristics are determined by minimizing the weight of the riser support structure subjected to stress constraints for the given target reliability. The initial design model is generated based on actual FPSO riser support specifications. The finite element analysis is conducted using NASTRAN, and the probabilistic optimal solutions are obtained via the moving least squares method in the context of RBDO using a response surface meta-model. For the meta-modeling of the inequality constraint functions of stresses, a constraint-feasible moving least squares method (CF-MLSM) is adopted in the present study. The CF-MLSM has been shown to ensure constraint feasibility regardless of the nonlinearity of the constraint function, the feasible bounds, and the random characteristics during the meta-model-based RBDO process. The solution results from the proposed RBDO strategy present improved design performances under various riser operating conditions.