Removing irregular frequencies by a partial discontinuous higher order boundary element method

The phenomenon of irregular frequencies is a puzzle in the course of calculating the interaction of waves and structures by the boundary element method. To remove the irregular frequencies, the modified integral domain method is adopted, and continuous higher order elements and partial discontinuous higher order elements are used for discretization. By these means, the effects of the irregular frequencies are effectively removed. Effective strategies have been adopted to deal with singular integrals and nearly singular integrals at different situations. The numerical results of the horizontal wave force on a uniform cylinder in the first order and second order diffraction problems show that the present method has a good validity. At the same time, the influence of collocation parameter on accuracy of numerical results is examined in detail.     

Application of a 2D harmonic polynomial cell (HPC) method to singular flows and lifting problems

Further developments and applications of the 2D harmonic polynomial cell (HPC) method proposed by Shao and Faltinsen [22] are presented. First, a local potential flow solution coupled with the HPC method and based on the domain decomposition strategy is proposed to cope with singular potential flow characteristics at sharp corners fully submerged in a fluid. The results are verified by comparing them with the analytical added mass of a double-wedge in infinite fluid. The effect of the singular potential flow is not dominant for added mass and damping, but the error is non-negligible when calculating mean wave loads using direct pressure integration. Then, the double-layer nodes technique is used to simulate a thin free shear layer shed from lifting bodies, across which the velocity potential is discontinuous. The results are verified by comparing them with analytical results for steady and unsteady lifting problems of a flat plate in infinite fluid. The latter includes the Wagner problem and the Theodorsen functions. Satisfactory agreement with other numerical results is documented for steady linear flow past a foil and beneath the free surface.     

Hydrostatic stability calculations by pressure integration

It is the purpose of this paper to present a new method for the calculation of hydrostatic properties of intact and damaged ship hulls and other floating structures. The geometry of the floating structure is modelled as a set of compartments, bounded by flat panels. Hydrostatic properties are derived by pressure integration on these panel elements. By application of Greens integration theorem the area integration is transformed into line integrals around the contour of each element. The line integrals can be expressed analytically such that the final result is that all the hydrostatic properties can be determined as a summation of easily evaluated expressions. This calculation procedure has the advantage of being well suited for computer calculation, and is easily applied to truly arbitrarily shaped floating structures. Another advantage is that it is based on a geometrical model which is equally suited for finite element strength analysis and hydrodynamic calculations based on a panel method. The application of the procedure is demonstrated by an example.     

下一代大气模式中的数值方法综述

数值预报是逐日天气预报、气候预测和气象防灾减灾的核心科技支撑。为进一步提高预报预测的准确度和延长预见期,甚高分辨率、多圈层耦合、多尺度嵌套、多尺度集合、数值地球系统模拟技术等是下一代数值预报的重要发展方向。异构众核高性能计算机和E级计算的高速发展为这一发展提供了契机,但也对现有业务数值预报中采用的数值方法提出了挑战。此文仅对国内外下一代大气模式涉及到的数值方法进行综述,着重于数值算法、准均匀球面网格和时间积分方案等3个方面,期望为相关研究者提供参考。     

Second-order Taylor expansion boundary element method for the second-order wave radiation problem

A novel Boundary Element Method (BEM) named the second-order Taylor Expansion Boundary Element Method (the 2nd order TEBEM) is developed for the solution of the second-order wave radiation velocity potential and sum-frequency wave loads for floating bodies. The radiation condition is enforced by a hybrid method of the multi-transmitting formula and damping zone. For the interior domain problem of a cube and a sphere, numerical results demonstrate that the 2nd order TEBEM can accurately solve the first and second-order gradients of velocity potential on the no-smoothed and smoothed boundary compared to the low-order BEM. The double frequency forces acting on the truncated cylinder are calculated under finite water depth. The agreement between the 2nd order TEBEM and others' numerical results is good. Moreover, all of the singular integrals in the 2nd order TEBEM can be solved analytically, so its implementation is much easier compared to the high-order BEM.     

Variational method for computing surface acoustic pressure on vibrating bodies, applied to transversely oscillating disks

A variational principle derived from the Kirchhoff-Helmholtz integral relation can be applied to acoustic radiation and diffraction problems. An illustrative example discussed here is that of sound radiation from a flat rigid circular disk in transverse oscillation. The variational formulation has the surface pressure as the unknown variable, with the normal velocity of the surface taken as given. The Rayleigh-Ritz method used in determining approximate solutions in terms of truncated expansions of basis functions encounters some numerical problems in the evaluation of integrals with singular integrands. The integrands are nevertheless integrable and techniques are described for handling the singularities. Another potential source of difficulty is that the tangential derivative of the surface pressure for the exact solution must be infinite at the edge of the disk. One makes use of prior knowledge of such a fact by using basis functions with the correct dependence on radial distance near the disk edge. Because basis functions in the Rayleigh-Ritz procedure have been selected with the aid of prior insight into the nature of the true solution, accurate results are obtained with a relatively small number of basis functions. The numerical solutions agree well with results calculated by Leitner in 1949.     

A Higher Order BEM for Wave-Current Action on Structures ——Direct Computation of Free-Term Coefficient and CPV Integrals

1 .Introduction Research on wave action onforward moving bodies or wave-current action on structures is impor-tant in offshore and coastal engineering.For bodies withsimple geometry,analytic solutions have beendeveloped . Matsuiet al .(1991) and Emmerhoff and Sclavounos (1992) derived analytic solutionsforuniformcylinders infinite andinfinite water depth.Bao and Kinoshita (1993) expendedthe theorytotruncated cylinders .In engineering practice ,structures are complexand numerical methods have…     

Simulating the influence of an atmospheric fine inhomogeneous structure on long-range propagation of pulsed acoustic signals

The results of simulating the influence of an atmospheric fine structure on the characteristics of acoustic signals propagating throughout the atmosphere for long distances from their sources are presented. A numerical model of an atmospheric fine inhomogeneous structure within the height range z = 20…120 km is proposed to perform calculations. This model and its numerical parameters are based on the current notions of the formation of an atmospheric fine structure due to internal gravity waves. The numerical calculations were performed using the parabolic-equation method. A spatial structure of the acoustic field and the structure of an acoustic signal at long distances from a pulsed source were calculated. It is shown that the presence of an atmospheric fine structure results in a scattering of acoustic signals and their recording in the geometric shadow region. The results of calculations of signal forms are in a satisfactory agreement with data on signals recorded in the geometric shadow region which is formed at a distance of about 300 km from an experimental explosion.     

An investigation into parametric roll resonance in regular waves using a partly non-linear numerical model

A partly non-linear time-domain numerical model is used for the prediction of parametric roll resonance in regular waves. The ship is assumed to be a system with four degrees of freedom, namely, sway, heave, roll and pitch. The non-linear incident wave and hydrostatic restoring forces/moments are evaluated considering the instantaneous wetted surface whereas the hydrodynamic forces and moments, including diffraction, are expressed in terms of convolution integrals based on the mean wetted surface. The model also accounts for non-potential roll damping expressed in an equivalent linearised form. Finally, the coupled equations of motion are solved in the time-domain referenced to a body fixed axis system.This method is applied to a range of hull forms, a post-Panamax C11 class containership, a transom stern Trawler and the ITTC-A1 containership, all travelling in regular waves. Obtained results are validated by comparison with numerical/experimental data available in the literature. A thorough investigation into the influence of the inclusion of sway motion is conducted. In addition, for the ITTC-A1 containership, an investigation is carried out into the influence of tuning the numerical model by modifying the numerical roll added inertia to match that obtained from roll decay curves.     

The asymmetric elastic wavefields in a model comprising a liquid layer overlying an anisotropic solid seabed due to an arbitrary source within the solid

The asymmetric three-dimensional radiation pattern and resultant elastodynamic response of stress waves in a model comprising a compressible water column overlying a transversely isotropic seabed in which a time-harmonic source acts is theoretically investigated. The use of potential functions, the Hankel transform, and a Fourier series expansion are adopted to deal with the equations of motion for both media. Closed-form integral expressions are developed for the potentials and the stress/displacement components. The expressions and introduced procedure are sufficiently flexible to incorporate various types of source loads. To evaluate the field quantities, the residue method and a robust integration scheme are utilized to handle the poles and branch points within the integrand. Any possible number of dispersive propagation modes are taken into account in the integral evaluation. The deduced velocity dispersion curves depict the characteristics of the various modes. They also indicate the existing singular points (poles) for a specific dimensionless frequency and the surface wave type associated with each pole. Numerical results are presented for the hydrodynamic pressure and displacement in the liquid layer and stress and displacement components in the solid seabed due to distributed and concentrated source excitations. The formulation and the numerical scheme are valid for calculating the wavefield anywhere within the model including both far- and near-field effects. The sensitivity of the results to different parameters is also analyzed. Both analytical and numerical comparisons with existing solutions for simpler cases are made to confirm the validity of the results. The results are especially useful in seismic hazard assessment of submarine earthquakes, landslides, and tsunamis. They can also be extended to deal with the fluid-solid-structure interaction problems.     

三维斜压台风模式Ⅰ.数值方法

一种多重移动套网格斜压台风模式已被应用于海洋环境数值预报.模式采用σ坐标系的原始方程组作为控制方程.现在用于国家海洋环境预报中心的模式垂直方向为非等距4层,水平方向为Arakawa B型格式,所采用的差分格式满足动量和能量守恒原理.模式控制方程组分离成平流过程和适应过程二组方程,并根据大气运动不同过程的特性,分别采用不同的时间步长和不同的积分方法.预报和后报结果显示该数值方法不仅可以缩短机时,而且可以得到稳定的预报结果.     

三维斜压台风模式I.数值方法

一种多重移动套网格斜压台风模式已被应用于海洋环境数值预报。模式采用坐标系的原始方程组作为控制方程。现在用于国家海洋环境预报中心的模式垂直方向为非等距４层,水平方向为ＡｒａｋａｗａＢ型格式,所采用的差分格式满足动量和能量守恒原理。模式控制方程组分离成平流过程和适应过程二组方程,并根据大气运动不同过程的特性,分别采用不同的时间步长和不同的积分方法。预报和后报结果显示该数值方法不仅可以缩短机时,而且可以得到稳定的预报结果。     

Geoacoustic Inversion via Genetic Algorithm and Its Application to Manganese Sediment Identification

An acoustic inversion method using a wide-band signal and two near field receivers is proposed and applied to multiple layered seabed models including a manganese sediment. The inversion problem can be formulated into a probabilistic model comprised of signals, a forward model, and additive noise. The forward model simulates wide-band signals, such as chirp signals, and is chosen to be the source-waveletconvolution plane wave modeling method. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the possible numerical ranges for a priori uniform distribution is based. The genetic algorithm is applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L 2 norm of the difference between measured and modeled signals. Not only the marginal pdf but also its statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm.     

Geoacoustic Inversion via Genetic Algorithm and Its Application to Manganese Sediment Identification

Abstract

An acoustic inversion method using a wide-band signal and two near field receivers is proposed and applied to multiple layered seabed models including a manganese sediment. The inversion problem can be formulated into a probabilistic model comprised of signals, a forward model, and additive noise. The forward model simulates wide-band signals, such as chirp signals, and is chosen to be the source-wavelet-convolution plane wave modeling method. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the possible numerical ranges for a priori uniform distribution is based. The genetic algorithm is applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L₂ norm of the difference between measured and modeled signals. Not only the marginal pdf but also its statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm.     

Spar平台垂荡板水动力系数计算与分析

Spar平台带有厚度很薄的垂荡板,使用边界元方法计算时存在着准奇异积分的数值计算问题。采用自适应高斯积分法对其准奇异积分问题进行处理,较少单元即可达到很高精度。在势流理论下对Spar平台垂荡板的水动力系数进行了计算,并与实验结果进行了对比,发现真实流体中垂荡板的附加质量约为势流理论下附加质量计算结果的1.1倍,而辐射阻尼在总阻尼中所占的百分比很小,此时粘性阻尼成为垂荡板的主要阻尼来源,在工程设计中不可忽略。     

Using a mesoscale ensemble to predict forecast error and perform targeted observation

Using NCEP short range ensemble forecast（SREF） system,demonstrated two fundamental on-going evolutions in numerical weather prediction（NWP） are through ensemble methodology.One evolution is the shift from traditional single-value deterministic forecast to flow-dependent（not statistical） probabilistic forecast to address forecast uncertainty.Another is from a one-way observation-prediction system shifting to an interactive two-way observation-prediction system to increase predictability of a weather system.In the first part,how ensemble spread from NCEP SREF predicting ensemble-mean forecast error was evaluated over a period of about a month.The result shows that the current capability of predicting forecast error by the 21-member NCEP SREF has reached to a similar or even higher level than that of current state-of-the-art NWP models in predicting precipitation,e.g.,the spatial correlation between ensemble spread and absolute forecast error has reached 0.5 or higher at 87 h（3.5 d） lead time on average for some meteorological variables.This demonstrates that the current operational ensemble system has already had preliminary capability of predicting the forecast error with usable skill,which is a remarkable achievement as of today.Given the good spread-skill relation,the probability derived from the ensemble was also statistically reliable,which is the most important feature a useful probabilistic forecast should have.The second part of this research tested an ensemble-based interactive targeting（E-BIT） method.Unlike other mathematically-calculated objective approaches,this method is subjective or human interactive based on information from an ensemble of forecasts.A numerical simulation study was performed to eight real atmospheric cases with a 10-member,bred vector-based mesoscale ensemble using the NCEP regional spectral model（RSM,a sub-component of NCEP SREF） to prove the concept of this E-BIT method.The method seems to work most effective for basic atmospheric state variables,moderately effective for convective instabilities and least effective for precipitations.Precipitation is a complex result of many factors and,therefore,a more challenging field to be improved by targeted observation.     

Stabilization of a linearly unstable wave participating in a three-wave resonant interaction

An analytical study of the influence of three-wave resonant interactions on the evolution of unstable wave disturbances is presented in the Kelvin-Helmholtz model. These results may be of interest in analyzing the dynamics of disturbances at the ocean-atmosphere interface and in two-layer flows which arise in the ocean and are characterized by large gradients of flow velocity at the boundary of layers. In the case under consideration, the instability arises when eigenfrequencies coincide in the framework of a single mode and the instability is algebraic. The amplitudes of the two other interacting stable waves are assumed to be small compared to the amplitude of the third, unstable, mode. The system of amplitude equations for this case is investigated using the WKB method. As a result, we obtain the formulas coupling the solutions for the time before and after a transition through a singular point, where the amplitude of the linearly unstable wave has a local minimum. These formulas give the rule of transformation of the parameter that characterizes a phase shift between fast and slow modes and determines the behavior of the system. It is shown that, in a transition through a singular point, this parameter changes randomly. As long as the parameter is positive, the amplitude of the linearly unstable wave remains limited and oscillates stochastically. In a transition of the parameter through zero, we exit the stabilization region and have an infinite growth of amplitude. The transition into the instability region is random. However, if the time interval where the amplitude remains limited is large enough, the scenario of the behavior of the system we have obtained can be treated as the partial stabilization of instability. The results make it possible for us to investigate the stochasticity caused by the nonlinear interaction of gravity-capillary waves in a two-layer model of a shear flow. These results are also of interest in analyzing secondary flows in laboratory facilities modeling the ocean and atmospheric processes.     

Separation of regular waves by a transfer function method

This paper presents a transfer function method (TFM) which can separate a regular wave field into incident and reflected waves based on the linear wave theory. The TFM uses specific transfer functions and corresponding convolution integrals to separate time series data measured in a combined partial standing wave system into incident and reflected waves. After this separation, estimation of the reflection coefficient becomes very easy. All manipulations have been performed in time domain. Furthermore, this method does not involve the calculation of wave heights and/or phase differences. The present method is demonstrated through numerical sample and physical model experiments carried out in a wave flume. Compared with other methods, the TFM gives much better estimates of the incident wave heights for physical model experiments in this study.     

Numerical investigation of the influence of the Caucasus relief on the distribution of hydrometeorological fields

The influence of the Caucasus relief on temporal and spatial variations in meteorological fields is investigated through the numerical integration of nonadiabatic, nonstationary, and nonlinear equations of atmospheric hydrothermodynamics. A dimensionless vertical coordinate is used. In numerical experiments, this coordinate makes it possible to approximate a real relief at a high accuracy and to take into account specific features of a complex orography of the region. The cases of flowing around the relief by large-scale westerly, easterly, northerly, and southerly background winds are modeled. The relief in the lower troposphere is shown to considerably change the spatial structures of large-scale fields of meteorological elements and to form the mesoscale structures of wind, temperature, pressure, and cloudiness that are characteristic of the Caucasus.     

Objective analysis of the structure of three-dimensional atmospheric fronts

A computational method for diagnosing three-dimensional atmospheric fronts from temperature, wind, and geopotential fields on a three-dimensional regular grid is proposed. The criterion, which serves for the diagnosis of atmospheric fronts, is discussed. The weights of the input information about the mentioned fields are optimized based on the maximal difference between the correlation functions for (a) pairs of particles separated by the front and (b) pairs from one synoptic mass. These weights were different for different baric levels. The correlation functions and the optimization of weights were estimated on the basis of the archive of fields of the NCEP objective analysis on the half-degree latitude-longitude grid and data from aerological observations. The results of numerical experiments on the construction of atmospheric fronts are presented. Applying the described method to fields predicted for a term of up to 36 h showed that errors in the prognostic models introduce a relatively weak distortion into the geometry of atmospheric fronts.