共查询到20条相似文献，搜索用时 203 毫秒
1.
研究并行算法解决应用并行计算机完成规模尽可能大的偏微分方程的数值求解问题。利用Hopf-Cole变换,将一维非线性Burgers方程转化为线性扩散方程,基于第二类Saul’yev型非对称格式和Crank-Nicolson格式对扩散方程进行差分离散,建立解Burgers方程的交替分段并行差分格式,并讨论该方法的稳定性,给出了数值算例。此算法把剖分节点分成若干组,在每组上构造能够独立求解的差分方程,因此具有并行本性,适合在高性能多处理器的并行计算机上使用。数值试验的结果表明此方法是有效的,且有较高的精度。 相似文献
2.
An improved coupling of numerical and physical models for simulating 2D wave propagation is developed in this paper. In the proposed model, an unstructured finite element model （FEM） based Boussinesq equations is applied for the numerical wave simulation, and a 2D piston-type wavemaker is used for the physical wave generation. An innovative scheme combining fourth-order Lagrange interpolation and Runge-Kutta scheme is described for solving the coupling equation. A Transfer function modulation method is presented to minimize the errors induced from the hydrodynamic invalidity of the coupling model and/or the mechanical capability of the wavemaker in area where nonlinearities or dispersion predominate. The overall performance and applicability of the coupling model has been experimentally validated by accounting for both regular and irregular waves and varying bathymetry. Experimental results show that the proposed numerical scheme and transfer function modulation method are efficient for the data transfer from the numerical model to the physical model up to a deterministic level. 相似文献
3.
An improved coupling of numerical and physical models for simulating 2D wave propagation is developed in this paper. In the proposed model, an unstructured finite element model (FEM) based Boussinesq equations is applied for the numerical wave simulation, and a 2D piston-type wavemaker is used for the physical wave generation. An innovative scheme combining fourth-order Lagrange interpolation and Runge-Kutta scheme is described for solving the coupling equation. A Transfer function modulation method is presented to minimize the errors induced from the hydrodynamic invalidity of the coupling model and/or the mechanical capability of the wavemaker in area where nonlinearities or dispersion predominate. The overall performance and applicability of the coupling model has been experimentally validated by accounting for both regular and irregular waves and varying bathymetry. Experimental results show that the proposed numerical scheme and transfer function modulation method are efficient for the data transfer from the numerical model to the physical model up to a deterministic level. 相似文献
4.
Efficient Numerical Solution of the Modified Mild-Slope Equation 总被引：12，自引：1，他引：11
An efficient numerical model for wave refraction,diffraction and reflection is presented in thispaper.In the model,the modified time-dependent mild-slope equation is transformed into an evolutionequation and an improved ADI method involving a relaxation factor is adopted to solve it.The methodhas the advantage of improving the numerical stability and convergence rate by properly determining therelaxation factor.The range of the relaxation factor making the differential scheme unconditionally stableis determined by stability analysis.Several verifications are performed to examine the accuracy of the pres-ent model.The numerical results coincide with the analytic solutions or experimental data very well,andthe computer time is reduced. 相似文献
5.
Presented here is a compact explicit difference scheme of high accuracy for solving the extended Boussinesq equations.For time discretization,a three-stage explicit Runge-Kutta method with TVD property is used at predicting stage,a cubic spline function is adopted at correcting stage,which made the time discretization accuracy up to fourth order;For spatial discretization,a three-point explicit compact difference scheme with arbitrary order accuracy is employed.The extended Boussinesq equations derived by Beji and Nadaoka are solved by the proposed scheme.The numerical results agree well with the experimental data.At the same time,the comparisons of the two numerical results between the present scheme and low accuracy difference method are made,which further show the necessity of using high accuracy scheme to solve the extended Boussinesq equations.As a valid sample,the wave propagation on the rectangular step is formulated by the present scheme,the modelled results are in better agreement with the experimental data than those of Kittitanasuan. 相似文献
6.
David F. Paskausky 《Journal of Oceanography》1978,34(5):233-236
The diffusive temperature equation is solved analytically and numerically for a simple oceanic case. Errors introduced by various numerical schemes are determined by comparison with the analytical solution. Less error is introduced by the standard differencing scheme than by the DuFort-Frankel scheme. A variable grid scheme produces the least error and the greatest savings in computer core and time. Usefulness of various numerical schemes is situation dependent. Any numerical forecasting technique using initialization data and frequent data updating will need to be verified with adequate field data before it can be considered reliable and used to predict short term (1 to 7 days) events in the ocean. 相似文献
7.
一个两时间层分裂显格式海洋环流模式(MASNUM)及其检验 总被引：1，自引：0，他引：1
HAN Lei 《海洋学报(英文版)》2014,33(11):11-35
A two-time-level, three-dimensional numerical ocean circulation model（named MASNUM） was established with a two-level, single-step Eulerian forward-backward time-differencing scheme. A mathematical model of large-scale oceanic motions was based on the terrain-following coordinated, Boussinesq, Reynolds-averaged primitive equations of ocean dynamics. A simple but very practical Eulerian forward-backward method was adopted to replace the most preferred leapfrog scheme as the time-differencing method for both barotropic and baroclinic modes. The forward-backward method is of second-order of accuracy, computationally efficient by requiring only one function evaluation per time step, and free of the computational mode inherent in the three-level schemes. This method is superior to the leapfrog scheme in that the maximum time step of stability is twice as large as that of the leapfrog scheme in staggered meshes thus the computational efficiency could be doubled. A spatial smoothing method was introduced to control the nonlinear instability in the numerical integration. An ideal numerical experiment simulating the propagation of the equatorial Rossby soliton was performed to test the amplitude and phase error of this new model. The performance of this circulation model was further verified with a regional（northwest Pacific） and a quasi-global（global ocean simulation with the Arctic Ocean excluded） simulation experiments. These two numerical experiments show fairly good agreement with the observations. The maximum time step of stability in these two experiments were also investigated and compared between this model and that model which adopts the leapfrog scheme. 相似文献
8.
This work presents a new approach for simulating the random waves in viscous fluids and the associated bottom shear stresses. By generating the incident random waves in a numerical wave flume and solving the unsteady two-dimensional Navier-Stokes equations and the fully nonlinear free surface boundaiy conditions for the fluid flows in the flume, the viscous flows and laminar bottom shear stresses induced by random waves axe determined. The deterministic spectral amplitude method implemented by use of the fast Fourier transform algorithm was adopted to generate the incident random waves. The accuracy of the numerical scheme is confirmed by comparing the predicted wave spectrum with the target spectrum and by comparing the nanlerical transfer function between the shear stress and the surface elevation with the theoretical transfer function. The maximum bottom shear stress caused by random waves, computed by this wave model, is compared with that obtained by Myrhaug＇ s model （1995）. The transfer function method is also employed to determine the maximum shear stress, and is proved accurate. 相似文献
9.
In this paper, we present a numerical procedure for solving a 2‐dimensional, compressible, and nonhydrostatic system of equations. A forward‐backward integration scheme is applied to treat high‐frequency and internal gravity waves explicitly. The numerical procedure is shown to be neutral in time as long as a Courant–Friedrichs–Lewy criterion is met. Compared to the leap‐frog‐scheme most models use, this method involves only two time steps, which requires less memory and is also free from unstable computational modes. Hence, a time‐filter is not needed. Advection and diffusion terms are calculated with a time step longer than sound‐wave related terms, so that extensive computer time can be saved. In addition, a new numerical procedure for the free‐slip bottom boundary condition is developed to avoid using inaccurate one‐sided finite difference of pressure in the surface horizontal momentum equation when the terrain effect is considered. We have demonstrated the accuracy and stability of this new model in both linear and nonlinear situations. In linear mountain wave simulations, the model results match the corresponding analytical solution very closely for all three cases presented in this paper. The analytical streamlines for uniform flow over a narrow mountain range were obtained through numerical integration of Queney's mathematical solution. It was found Queney's original diagram is not very accurate. The diagram had to be redrawn before it was used to verify our model results. For nonlinear tests, we simulated the famous 1972 Boulder windstorm and a bubble convection in an isentropic enviroment. Although there are no analytical solutions for the two nonlinear tests, the model results are shown to be very robust in terms of spatial resolution, lateral boundary conditions, and the use of the time-split scheme. 相似文献
10.
This paper aims to validate a numerical seakeeping code based on a 3D Rankine panel method by comparing its results with experimental data. Particularly, the motion response and hull-girder loads on a real modern ship, a 6500 TEU containership, are considered in this validation study. The method of solution is a 3D Rankine panel method which adopts B-spline basis function in the time domain. The numerical code is based on the weakly nonlinear scheme which considers nonlinear Froude-Krylov and restoring forces. The main focus of this study is given to investigate the nonlinear characteristics of wave-induced loads, and to validate this present scheme for industrial use in the range of low Froude number. The comparisons show that the nonlinear motions and hull-girder loads, computed by the present numerical code, have good overall agreements with experimental results. It is found that, for the better accuracy of computational results, particularly at extreme waves in oblique seas, the careful treatment of soft-spring (or compatible) system is recommended to the control of non-restoring motions such as surge, sway, and yaw. 相似文献
11.
The problem of static three-dimensional, nonlinear, large deformation response of a marine riser is formulated within small strain theory and solved numerically. This type of analysis is necessary, for the new generation of drilling and production risers. The mathematical model takes properly into account the effects of internal and external pressure and the complete nonlinear boundary conditions, without linearizing the follower forces. The extensibility or inextensibility condition is used as the constitutive relation in the tangential direction. Torsion and bending are coupled. The external load and the boundary conditions are deformation dependent. A solution method is developed based on an incremental finite element algorithm, which involves a prediction-correction scheme. In the correction phase deformation dependent quantities are updated. The extensibility or inextensibility condition is used to reduce the degrees of freedom of the system. The numerical results of the developed computer code compare very well with available semi-analytical and numerical solutions. Three numerical applications are used to demonstrate the importance of large deformation, nonlinear and three-dimensional analyses. 相似文献
12.
海岸河口三维潮流数学模型 总被引：1，自引：0，他引：1
本文采用有限元法建立了适用于海岸河口浅水地区的三维潮流数学模型，垂向采用绝对分层坐标系统，将整个水柱分成若干层，在每层内通过垂向积分平均，将三维问题简化为多个平面二维问题。在求解有限元方程中，引入集中质量矩阵技术，在时间上采用两步ＬａｘＷｅｎｄｒｏｆｆ格式，使有限元方程直接以显式解出，不需联立求解，节省了大量的计算时间和计算内存，通过模拟表面切应力作用在矩形水地上而引起的水流，计算结果与分析解比较一致，并将本模型应用到香港维多利亚水道中，计算结果与实测值亦符合较好，证明本模型是一个实用而有效的三维潮流数学模型。 相似文献
13.
14.
Computation of wave-making resistance of a catamaran in deep water using a potential-based panel method 总被引：1，自引：0，他引：1
This paper presents a potential-based boundary element method for solving a nonlinear free-surface flow problem for a Wigley catamaran moving with a uniform speed in deep water. Since the interior flow of each monohull of the catamaran is different from the exterior flow, both monohulls must be considered as lifting bodies. The pressure Kutta condition is imposed at the trailing-edge of the lifting body by determining the dipole distribution, which generates required circulation on the lifting part. The effects of wave interference and hull separation on the hydrodynamic characteristics of the catamaran hull are analyzed and the validity of the computer scheme is examined by comparing the wave resistance with the numerical results of others. The present method could be a useful design tool for screening the suitable combinations of hull parameters and hull spacing at the preliminary design stage of catamaran hull. 相似文献
15.
将跳点法的基本思想应用于求解二维水位抛物型方程，从而构造了一种用于风暴潮数值预报的高效率差分格式。文中给出了这一格式的稳定性证明，并通过对实际风暴潮过程的模拟，证实了这一格式的稳定性和时效性 相似文献
16.
A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed. 相似文献
17.
A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropr... 相似文献
18.
The present study develops a numerical model of the two-dimensional fully nonlinear shallow water equations （NSWE） for the wave run-up on a beach. The finite volume method （FVM） is used to solve the equations, and a second-order explicit scheme is developed to improve the computation efficiency. The numerical fluxes are obtained by the two dimensional Roe＇ s flux function to overcome the errors caused by the use of one dimensional fluxes in dimension splitting methods. The high-resolution Godunov-type TVD upwind scheme is employed and a second-order accuracy is achieved based on monotonic upstream schemes for conservation laws （MUSCL） variable extrapolation; a nonlinear limiter is applied to prevent unwanted spurious oscillation. A simple but efficient technique is adopted to deal with the moving shoreline boundary. The verification of the solution technique is carried out by comparing the model output with documented results and it shows that the solution technique is robust. 相似文献
19.
数值模式与统计模型相耦合的近岸海浪预报方法 总被引：2，自引：2，他引：0
针对数值模式和统计模型预报近岸海浪存在的局限性,构建了数值模式和统计模型相耦合的近岸海浪预报框架,在模式计算格点和近岸预报目标点之间定义一个海浪能量密度谱传递系数,通过经验正交函数分解和卡尔曼滤波方法建立传递系数的统计预报模型并与数值模式进行耦合。经过对近岸波浪观测站1a的预报试验表明:该方法能够提高近岸海浪有效波高预报精度,有效波高的均方根误差降低了约0.16m,平均相对误差降低约9%。进一步试验和分析发现,该方法的预报有效时间小于24h,将海浪能量密度谱经过分解后得到的基本模态反映了近岸波侯的主要特征,海浪能量密度谱传递系数的变化体现了波侯的季节变化特点。 相似文献
20.
为建立高效的Boussinesq类水波数值模型,提出了一种新型的、基于有限差分和有限体积方法的混合数值格式。针对守恒形式的一维控制方程,在等间距矩形控制体内对其进行积分并离散,采用有限体积方法计算界面数值通量,剩余源项采用有限差分方法计算。其中,采用MUSTA格式并结合高精度状态插值方法计算控制体界面数值通量。时间积分则采用具有TVD性质的三阶龙格-库塔多步积分法进行。除验证模型外,重点对MUSTA格式和广泛使用的HLL格式进行了比较。结果表明,MUSTA格式可用于Boussinesq类水波方程数值求解,综合考虑数值精度、计算效率、程序编制和实际应用这几个方面,其较HLL格式更具有优势。 相似文献