预修正快速傅里叶变换高阶边界元方法在波浪对物体作用问题中的应用

边界元方法被广泛应用于波浪对海上婕筑物作用领域,但由于传统边界元方法的存储量和计算量均为未知量的平方量级,很难满足大范围多未知量问题的计算需要.本文基于高阶边界元方法,应用预修正快速傅里叶变换方法,使计算量与存储量分别降低至O(NlogN)和O(N)量级,并得到一个连续的压强分布以适应结构设计的要求,同时可以通过使用满...     

预修正快速傅里叶变换方法在波浪绕射和辐射问题中应用

在波浪对海上建筑物绕射和辐射问题的计算中,面元法被广泛使用,但由于传统面元法的存储量和计算量均为未知量的平方量级,很难满足大范围多未知量问题的计算需要。采用预修正快速傅里叶变换方法(pFFT方法),使计算量与存储量都降低到未知量的线性量级。以淹没圆球与漂浮圆柱两个典型算例为基础,通过不同未知量时pFFT方法与传统面元法的计算量与存储量的对比,以及pFFT方法自身各步骤计算时间的对比,研究了不同网格方案的选取对pFFT方法计算量和存储量的影响,推荐根据未知量个数采用计算时间最小化原则选取pFFT网格参数。     

B样条面元法在波浪载荷计算中若干问题研究

对基于B样条面元法波浪载荷计算、载荷加载及速度势的二阶偏导数计算中存在若干问题的解决办法进行了探讨.对于加载时压力点的参数值计算,给出一种基于最小二乘法的参数值反算方法;对于二阶偏导数计算,给出一种基于B样条面元法求解速度势二阶导数的隐式方法,并通过对一无限水域中椭球潜体定常速度势二阶导数的计算对该方法进行了验证.此外,对B样条面元法的计算效率进行了讨论,并对对称性在B样条面元法中的应用以及其对计算效率的提高进行了说明.所介绍的方法简单实用,对于发展和推广B样条面元法在波浪载荷计算中的应用具有一定的价值.     

基于PCL语言的浮体湿表面B样条面元生成方法

提出一种基于PCL语言生成浮体湿表面B样条面元的方法.介绍利用PCL语言生成浮体湿表面型值数据点,对型值数据点进行规范参数化,反算B样条面元控制点,以及生成由控制点计算得到的用于检查B面元质量的加密曲面网格这四个主要模块的程序流程,并对编写该程序时需要用到的PCL函数及注意事项进行说明.该方法充分利用MSC.Patran强大的图形处理及三维显示功能,可以准确生成任意形状浮体湿表面的B样条面元,为获得准确的高阶水动力结果奠定了基础.     

信号快速带通滤波方法在海洋信号分析中的应用

借助于离散"解析"信号提出了一种对信号进行快速滤波的方法.与传统方法比较,手续简单,快速有效,尤其适用复杂多变的海洋长信号的滤波.以2个典型信号的带通滤波为例,验证了快速滤波方法的有效性,2个复杂海洋信号的滤波结果表明,此滤波方法在海洋信号分析中有应用意义.     

海水入侵地区咸水与淡水分界面计算──B样条函数方法的应用

提出用Ｂ样条函数方法求解咸、淡水分界面的原理、方法和步骤，指出该方法原理简明、易于实现、计算可靠。作为实例，将此方法应用于莱州市滨海地区。     

基于B—样条的边界元方法及其在波浪力计算中的应用

就波浪与结构物相互作用问题,应用Ｂ－样条函数对物体表面上的速度势进行了数值离散,物体表面的控制参数通过最小二乘原理确定,速度势控制参数通过Ｇａｌｅｒｋｉｎ理论确定。本理论给 其导数在整个光滑物面上都是连续的。采用此种方法所计算的数值解十分准确,并且收敛得很快。     

对流扩散方程欧拉──拉格朗日方法的B─样条插值格式

利用Ｂ－样条函数．文中发展了一个对流扩散方程欧拉－拉格朗日方法的插值格式．改善了其数值频散性。算例表明，此格式优于一次和二次插值格式。文中给出了Ｂ－样条插值格式的二维及三维推广。     

基于薄板样条函数的HY-1B卫星 L1B 数据几何校正方法研究

将HY-1BL1B产品中已有的经纬度数据抽样作为控制点,根据控制点的图像坐标和经纬度坐标间的对应关系,基于薄板样条函数的坐标变换关系对HY-1BL1B数据进行了几何校正,并与陆地矢量边界进行了比较,结果表明该方法能较好地完成HY-1B卫星L1B数据的几何校正。     

基于B样条函数的不定积分变换法在磁力资料处理中的应用

磁异常的不定积分变换法是1种提高磁异常处理及解释能力的方法.提出利用3次B样条插值函数对等距节点离散磁异常进行不定积分变换,实现了对实测数据处理和不定积分参量图的应用,验证了该方法的正确性和可行性.选取胶州湾湾口区2条典型测线,对磁异常数据进行变换得到不定积分参量图,将实测磁异常与模型磁异常的剖面图、参量图、不定积分参量图综合对比得出,不定积分参量图在判断地质体参数值方面有着明显的优势,利用3种图形综合解释可进行定量解释并减少解释多解性.     

A B-spline based higher order panel method for analysis of steady flow around marine propellers

A higher order panel method based on B-spline representation for both the geometry and the solution is developed for the analysis of steady flow around marine propellers. The self-influence functions due to the normal dipole and the source are desingularized through the quadratic transformation, and then shown to be evaluated using conventional numerical quadrature. By selecting a proper order for numerical quadrature, the accuracy of the present method can be increased to the machine limit. The far- and near-field influences are shown to be evaluated based on the same far-field approximation, but the near-field solution requires subdividing the panels into smaller subpanels continuously, which can be effectively implemented due to the B-spline representation of the geometry. A null pressure jump Kutta condition at the trailing edge is found to be effective in stabilizing the solution process and in predicting the correct solution. Numerical experiments indicate that the present method is robust and predicts the pressure distribution on the blade surface, including very close to the tip and trailing edge regions, with far fewer panels than existing low-order panel methods.     

NURBS高阶面元法在多体水动力分析中应用

当位于海上的多个浮体距离相近时,它们间将产生显著的水动力相互作用。使用NURBS(非均匀有理B样条)高阶面元法进行多体水动力分析,速度势用任意阶的B样条函数近似,物面用NURBS曲面表示,由于水动力与几何表示相分离,就可以方便地实现计算机辅助设计与水动力分析的集成。对截断圆柱与椭球间的水动力相互作用和半潜式平台与箱型驳船间的水动力相互作用的计算显示此种方法可以对多体间的水动力相互作用给出准确的分析,计算时间和占用内存大幅减少。     

Downward continuation of aeromagnetic data to the marine surface

The existing methods of downward continuation of potential field cannot be used to continue the aeromagnetic data to the marine surface because of the limited continuation distance. An iteration method for the downward continuation of potential field with a much larger continuation distance has been developed, which can continue the aeromagnetic data to the marine surface and get the marine - magnetic chart with the same scale as the aeromagnetic data. This downward continuation method will greatly raise the ef- ficiency of marine - magnetic investigation. The principle of the iteration method is presented. The method is demonstrated on synthetic models and real aeromagnetic data. Also, the error brought by continuation is discussed. The efficiency of the iteration method for the downward continuation of potential field is compared with the fast fourier transform （FFT） method, the former is much better than the latter.     

A 3D potential-based and desingularized high order panel method

In this paper, a novel high order panel method based on doublet distribution and Gaussian quadrature was adopted to deal with the potential flow problem. In the geometry representation we employed both the exact surface and NURBS surface form to construct the surface panel. These data were calculated directly from the mathematical shape definition. Furthermore, no fixed order of doublet density distribution was assumed on each panel. Not only the number of panels could be chosen, but also the Gaussian order of each panel. The numerical results for sphere, ellipsoid and Wigley hull demonstrated here indicated that the present method was adapted to the potential flow problem. Moreover, the NURBS surface geometry representation was capable of further potential flow optimal calculation.     

数字滤波技术在潮汐观测中的应用研究

潮汐是重要的海洋物理要素。瞬时变化的海面高度信息除了包含潮汐信息外,还包括波浪等高频干扰信息,因此需要对干扰信息进行滤除。如何从获取到的海面瞬时高度变化中提取出潮汐信息就成为一个关键问题。海洋学中通常是通过低通滤波来获取低频的潮汐信息。通过分析常用的3种低通数字滤波的方法,即滑动平均法、快速傅里叶变换法和小波变换法,对3种低通滤波方法应用于潮汐数据处理的效果进行了比较分析,得出滑动平均法在潮汐信息提取中具有较高应用价值的结论。     

A B-spline-based method for radiation and diffraction problems

An open uniform B-spline-based panel method is developed for solution of potential flow problems. In this method, both geometry as well as the field variables are represented by the same open uniform B-spline basis function. The method is initially applied for the radiation problem in unbounded fluid. Computed results for a spheroid of different aspect ratio are found to be in excellent agreement with analytical results. The method is then applied for diffraction problem formulated based on the transient (time-domain) Green's function. Computed results for a hemisphere and Wigley hull are compared with published results and the comparison shows good agreement.     

基于频域Rankine面元法的船舶耐波性计算和分析

在频域中分别使用Rankine面元法和去奇点Rankine面元法（desingularized-Rankine panel method，简称DRPM）快速求解有航速船舶耐波性问题。使用两种不同的线性化方式：叠模流线性化法（double-body flow，简称DB）和均匀流线性化方法（Neumann-Kelvin，简称NK）计算了Wigley I、Wigley III和S175在有航速时船舶的水动力系数、波浪激励力和运动响应。两种线性化方法的计算结果与试验结果进行了比较，研究表明：使用Rankine面元法和去奇点Rankine面元法计算的结果相差不大，且与NK法相比使用DB法计算的结果和试验值更吻合，尤其是非对角线上交叉耦合水动力系数A35,B35,A53,B53的计算结果。运动响应对于边界条件线性化方式比较敏感，尤其是垂荡运动，在船共振频率附近，运动受到的影响最明显。     

Design of ship's bow form by potential-based panel method

A surface panel method treating a boundary-value problem of the Dirichlet type with the free surface is presented to design a three-dimensional body corresponding to a prescribed pressure distribution. The free surface boundary condition is linearized with respect to the oncoming flow, and computed by four-point finite difference scheme. Sample designs for submerged spheroids and Wigley hull are carried out to demonstrate the stable convergence, the effectiveness and the robustness of the method. The design of a 5500TEU container carrier is performed with respect to reduction of the wave resistance. To reduce the wave resistance, calculated pressure on the hull surface is modified to have the lower fluctuation, and is applied as a Dirichlet type dynamic boundary condition on the hull surface. The designed hull form is verified to have the lower wave resistance than the initial one not only by computation but also by experiment.     

Flow around circular and noncircular cylinders by a low-order panel method

Time-dependent cross-flow was studied around cylinders with circular and noncircular cross-sections. The numerical approach for the analysis was a low-order panel method based on constant source and dipole values along each panel. The method was previously used successfully for several applications, such as calculation of the added mass and damping coefficients. In simulating the viscous time-dependent flow around the cylinder, the time-dependent wake feature of the code was used. For the circular and D-cylinders, the results agreed well with the experiments. Suggestions for improving the results for T-cylinders with angle of attack are included.