共查询到20条相似文献，搜索用时 106 毫秒
1.
A high-order Rankine panel method based on Non-Uniform Rational B-Spline (NURBS) is developed for solving the three-dimensional radiation and diffraction problems with forward speed. A NURBS surface is used to precisely represent the body geometry. Velocity potential on the body surface is described by B-spline after the source density distribution on the boundary surface is determined. A collocation approach is applied to numerical computation and the integral equations are evaluated by applying Gauss–Legendre quadrature. The mj-terms are evaluated by a desingularized method which utilizes NURBS technique. In order to verify the method proposed, it is firstly applied to the unbounded flow problem of a sphere and spheroids. The numerical results are found to be in good agreement with analytical solutions. Then the method is used to solve the radiation and diffraction problems of a sphere and the diffraction problem of a spheroid moving with a forward speed beneath the free surface in frequency domain. The numerical results are satisfactory in comparison with the published analytical results and experimental results. 相似文献
2.
A panel-free method (PFM), based on the desingularized Green’s formulae proposed by Landweber and Macagno, has been developed to solve the radiation problem of a floating body in the time domain. The velocity potential due to a non-impulsive velocity is obtained by solving the boundary integral equation in terms of source strength distribution. The singularity in the Rankine source term of the time-dependent Green function is removed. The geometry of a body surface is mathematically represented by NURBS surfaces. The integral equation can be globally discretized over the body surface by Gaussian quadratures. No assumption is needed for certain degree of approximation of distributed source strength on the body surface. The accuracy of PFM was demonstrated by its application to a classical problem of uniform flow past a sphere. The response function of a hemisphere at zero speed was then computed by PFM. The computed response function, added-mass and damping coefficients are compared with other published results. 相似文献
3.
提出求解地下流稳定问题的一种新方法——样条函数法,该法较其它算法有计算精度高,程序易于编制,边界处理灵活等优点。 相似文献
4.
Gennady A. Kivman 《Progress in Oceanography》1997,40(1-4)
Estimation of the open-boundary inputs solving a weak constraint variational formulation for an Arctic tide model is considered as an ill-posed problem in the sense that the solution is very sensitive to the data noise and to grid size. Mathematically, spatial discretization of a cost function to be minimized and penalization of normal flow through the open boundary act as regularization of the problem. An heuristic choosing rule for the regularization parameter is applied to assess a suitable spatial resolution and the weight referred to the open boundary penalty. It is shown that these provide a better fit of the solution to a control data set compared with a finer grid, the value of the energy flux through the open boundary being in agreement with other model estimates. The M2 solution obtained is much closer to the control data than other modern solutions while the accuracy of the simulated K1 constituent is within the same error level. The tidal maps for these waves exhibit certain distinctions in comparison with other charts. 相似文献
5.
A computational method has been developed to predict the hydrodynamic performance of the propeller–rudder systems (PRS) and azimuthing podded drive (AZIPOD) systems. The method employs a vortex-based lifting theory for the propeller and the potential surface panel method for the steering system. Three propeller models along with three steering systems (rudder and strut, flap and pod (SFP)) are implemented in the present calculations for the cases of uniform and non-uniform conditions. Computed velocity components show good agreement with the experimental measurements behind a propeller with or without the rudder. Calculated thrust, torque and lift also agree well with the experimental results. Computations are also performed for an AZIPOD system in order to obtain the pressure distributions on the SFP, and the hydrodynamic performance (thrust, torque and lift coefficients). The present method is useful for examining the performance of the PRS and AZIPOD systems in the hope of estimating the propulsion and the maneuverability characteristics of the marine vehicles more accurately. 相似文献
6.
The title problem is solved in the case where uniform applied loading is present at the plate outer boundary. Two independent solutions are obtained: the optimized Rayleigh–Ritz method and the finite element algorithmic procedure. Good engineering agreement is shown to exist. The proposed analytical procedure is convenient from a designer's viewpoint. Admittedly exact analytical solutions can be found in terms of Bessel functions, the procedure being rather complex from an engineering viewpoint. 相似文献
7.
The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from lO0 to lO00. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem. 相似文献
8.
一种基于Roe格式的有限体积法在二维溃坝问题中的应用 总被引：3，自引：0，他引：3
溃坝属于典型的非恒定含间断的浅水问题。应用有限体积法离散二维浅水控制方程的守恒型方程组,将基于近似黎曼解的Roe格式用于数值计算溃坝问题,并利用MUSCL方法构造二阶空间积分格式和预测-校正二步法构造二阶时间格式,从而使数值解的整体达到二阶,提高了精度。文中算法在一维溃坝的Stoker问题的数值结果与解析解进行对比,结果证明了此方法的可行性。应用此方法在二维溃坝问题上的结果,说明了此算法可有效模拟溃坝水流的演进过程。 相似文献
9.
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. 相似文献
10.
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 L2 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. 相似文献
11.
Domain decomposition and matching for time-domain analysis of motions of ships advancing in head sea
A domain decomposition and matching method in the time-domain is outlined for simulating the motions of ships advancing in waves. The flow field is decomposed into inner and outer domains by an imaginary control surface, and the Rankine source method is applied to the inner domain while the transient Green function method is used in the outer domain. Two initial boundary value problems are matched on the control surface. The corresponding numerical codes are developed, and the added masses, wave exciting forces and ship motions advancing in head sea for Series 60 ship and S175 containership, are presented and verified. A good agreement has been obtained when the numerical results are compared with the experimental data and other references. It shows that the present method is more efficient because of the panel discretization only in the inner domain during the numerical calculation, and good numerical stability is proved to avoid divergence problem regarding ships with flare. 相似文献
12.
Iskender Sahin Assistant Professor Allen H. Magnuson Asssciate Professor 《Ocean Engineering》1984,11(5):451-461
A modified source-and-dipole type singularity panel method is proposed to calculate the flow properties for an oscillating arbitrary body in the presence of a free surface. The technique is based on Green's identity whereby the boundary value problem is expressed as a boundary integral equation which is solved numerically. The free-space Green function is used in the integral equation. To demonstrate the feasibility of the method, the problem of a pulsating submerged line source under a free surface is treated and results are compared with the exact solution.An excellent agreement with the theory is obtained for panel density of about ten panels per wavelength and paneled water surface length of two wavelengths with very low computing times, indicating the feasibility of the method for unsteady water wave problems. 相似文献
13.
The paper provides a detailed analysis for the second-order diffraction of monochromatic waves. For the second-order potential on the free surface, the paper proposed a forward prediction method for computing the integration on the free surface. By this method we only need to run the infinity integration on the free surface directly for a few points; a one-step quadrature can then be applied successively outward from the body for potentials at other points. For wave diffraction from a body of revolution with a vertical axis, the paper derives a new integral equation, which can cancel the leading singularity in the derivative of ring Green's functions automatically. To obtain accurate results, different approaches are also used to deal with singularities in the ring Green's functions in the integration on both the body surface and free surface. The method has been implemented for bodies of revolution with vertical axes, but the theory is also available for arbitrary bodies.A numerical examination is made to validate the numerical code by comparing the second-order forces and moments on uniform and truncated cylinders and second-order diffraction potentials on the free surface with some published results. The comparison shows that the present results are in good agreement with those published. The method is also used to compute the second-order wave elevation around uniform and truncated cylinders. 相似文献
14.
Numerical solutions for the hydroelastic problems of bodies are studied directly in the time domain using Neumann–Kelvin formulation. In the hydrodynamic part of problem, the exact initial boundary value problem is linearized using the free stream as a basis flow, replaced by the boundary integral equation applying Green theorem over the transient free surface Green function. The resultant boundary integral equation is discretized using quadrilateral elements over which the value of the potential is assumed to be constant and solved using the trapezoidal rule to integrate the memory or convolution part in time. In the structure part of the problem, the finite element method is used to solve the hydroelastic problem. The Mindlin plate as a bending element, which includes transverse shear effect and rotary inertia effect are used. The present numerical results show acceptable agreement with experimental, analytical, and other published numerical results. 相似文献
15.
1 IntroductionThe shallow water equations (SWE) are frequent-ly used as a mathematical model for water flows incoastal areas, lakes, estuaries, etc. Thus, they are animportant tool to simulate a variety of problems relat-ed to coastal engineering, environment, ecology, etc.(Bermúdez et al., 1998). On the basis of solving theone-dimensional (1D) SWE, Hu et al. (2000) have de-veloped a model capable of simulating storm wavespropagating in the coastal surf zone and overtopping asea wall. Ano… 相似文献
16.
On the basis of the reproducing kernel particle method (RKPM), a new meshless method, which is called the complex variable reproducing kernel particle method (CVRKPM), for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM. 相似文献
17.
A method is introduced to calculate and to account for the uncertainties in the predictions of oil spill trajectories using a classic oil spill model. The method considers the output of the oil spill model as a function of random variables, which are the input parameters, and calculates the standard deviation of the output results which gives a measure of the uncertainty of the model given the uncertainties of the input parameters.Instead of a single trajectory that is calculated by the oil spill model using the mean values of the parameters, a band of trajectories can be defined when various simulations are done taking into account the uncertainties of the input parameters. This band of trajectories defines envelopes of the trajectories that are likely to be followed by the spill given the uncertainties of the input.The method is applied to an oil spill that occurred in open sea near Madeira Islands, in the Atlantic Ocean, in December 1989. The simulations allow the understanding of how a change in the wind direction drove the spill towards the Islands.The envelope of likely trajectories that is obtained with the uncertainty modelling shows a band of trajectories that is in better agreement with the observations than the single trajectory simulated by the oil spill model, based on mean parameters. 相似文献
18.
19.
The flow relaxation scheme (FRS) is tested as an open boundary condition (OBC). This scheme was originally designed to relax external solutions (ES) from a large area model towards solutions in a limited area model with a fine mesh. When one uses the FRS as a pure OBC; i.e. the ES are unknown and set equal to zero, the FRS method degenerates to a sponge type OBC. Since a sponge type OBC does not work well in all cases, for instance in the case of an alongshore uniform wind; Chapman (1985) and Røed and Cooper (1987), we introduced local wind induced solutions as a part of the ES.When we also included the tide in the ES, the tests showed that the FRS method handled the tidal input problem very well. Even in the case where noise was added to the tide in the ES, the results were good.The major drawbacks for the FRS method are the required extention of the computer-space and -time. 相似文献
20.
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… 相似文献