共查询到20条相似文献，搜索用时 218 毫秒
1.
边界拟合坐标下流场计算模式 总被引：1，自引：0，他引：1
流场数值模拟中，采用曲线坐标法准确模拟计算水域自然边界情形是当今研究的热点之一。本文借助于自然界中某些物理现象中其函数符合椭圆型且必然正交这一特性，建立了曲线坐标系，弥补在笛卡尔坐标系下二维流场模拟的边界贴合不足。应用本模式对我国某一河口水域的潮位和流场进行计算，结果表明；该模式具有物理意义清晰、稳定性良好、精度高等特点，模式对实际工程中流场模拟具一定的应用价值。 相似文献
2.
V. B. Zalesnyi A. V. Gusev V. I. Agoshkov 《Izvestiya Atmospheric and Oceanic Physics》2016,52(3):277-293
We present a numerical model of Black Sea circulation based on primitive equations with improved spatial resolution in the coastal zone. The model equations are formulated in a two-pole orthogonal coordinate system with arbitrary locations of the poles and a vertical σ coordinate. Increased horizontal resolution is gained by displacing the pole into the vicinity of the separated subdomain. The problem is solved over a grid with a variable step. The northern coordinate pole is displaced to the vicinity of Gelendzhik; the grid step varies from 150 m in the coastal zone to 4.6 km in the main basin. We simulated the fields of currents, sea level, temperature, and salinity under the given atmospheric forcing in 2007. The model is capable of reproducing the large-scale Black Sea circulation and submesoscale variations in the coastal currents. 相似文献
3.
《Coastal Engineering》2005,52(6):513-533
Using the perturbation method, a time dependent parabolic equation is developed based on the elliptic mild slope equation with dissipation term. With the time dependent parabolic equation employed as the governing equation, a numerical model for wave propagation including dissipation term in water of slowly varying topography is presented in curvilinear coordinates. In the model, the self-adaptive grid generation method is employed to generate a boundary-fitted and varying spacing mesh. The numerical tests show that the effects of dissipation term should be taken into account if the distance of wave propagation is large, and that the outgoing boundary conditions can be treated more effectively by introduction of the dissipation term into the numerical model. The numerical model is able to give good results of simulating wave propagation for waters of complicatedly boundaries and effectively predict physical processes of wave propagation. Moreover, the errors of the analytical solution deduced by Kirby et al. (1994) [Kirby, J.T., Dalrymple, R.A., Kabu, H., 1994. Parabolic approximation for water waves in conformal coordinate systems. Coastal Engineering 23, 185–213.] from the small-angle parabolic approximation of the mild-slope equation for the case of waves between diverging breakwaters in a polar coordinate system are corrected. 相似文献
4.
数值模式与统计模型相耦合的近岸海浪预报方法 总被引：2，自引：2，他引：0
针对数值模式和统计模型预报近岸海浪存在的局限性,构建了数值模式和统计模型相耦合的近岸海浪预报框架,在模式计算格点和近岸预报目标点之间定义一个海浪能量密度谱传递系数,通过经验正交函数分解和卡尔曼滤波方法建立传递系数的统计预报模型并与数值模式进行耦合。经过对近岸波浪观测站1a的预报试验表明:该方法能够提高近岸海浪有效波高预报精度,有效波高的均方根误差降低了约0.16m,平均相对误差降低约9%。进一步试验和分析发现,该方法的预报有效时间小于24h,将海浪能量密度谱经过分解后得到的基本模态反映了近岸波侯的主要特征,海浪能量密度谱传递系数的变化体现了波侯的季节变化特点。 相似文献
5.
The slightly compressible flow formulation is applied to the free-surface, three-dimensional turbulent flow around a Wigley hull. Two turbulence models (large eddy simulation and Baldwin–Lomax) are used and compared. The simulation conditions are the ones for which experimental and numerical results exist. The computational grid is built using an algebraic grid generator with the model fixed in space. The codes use the interface-capturing technique for computing the free-surface displacements and the Beam and Warming scheme for marching in time the numerical model. The results compare well with the experimental data available. 相似文献
6.
基于2000年8月在南海调查航次得到的水文资料,首次采用广义随底坐标形式的改进POM模式对南海夏季环流进行了数值研究.用正交曲线性水平网格覆盖观测区域,在垂向上对近表海面层次采用近似z坐标,而近底层则为随底坐标.在计算海区实际地形及假设的水平均匀而垂直层化的密度分布下,实施的两个数值计算试验表明,本模式采用的垂直坐标方案比传统的σ坐标方案优越,随底坐标模式因压力梯度项在起伏地形下产生的系统计算误差将变得十分的微小.在南海2000年夏季环流的实际计算中,首先对观测资料进行了60d的诊断计算,然后在诊断已得到的动力场结果基础上,又进行了10d左右的预报运行得到半诊断结果.从计算结果来看,它依赖于参数Cvis与Cdif的选择,特别是参数Cvis,文中取值为Cvis=Cdif=008.比较诊断与半诊断两个计算过程的结果,它们在定性上较为一致,在定量上有些差别.这是因为半诊断计算的方法对密度场作适当的动力调整,使其与地形、风场等更加匹配.在大尺度环流结构不受影响的情况下,尽可能地消除了小尺度噪声,可使计算得到的流场更为清晰.2000年8月南海计算区域环流的最大特点是多涡结构,其中有些反气旋暖涡和气旋式冷涡相间分布.在越南东南海域自表层至1000m水层稳定存在着一个显著的反气旋暖涡,其中心位置在11°51'N,112°07'E(诊断计算),水平尺度约为300km.此暖涡以东存在一个气旋式冷涡,这两个冷、暖涡是研究海区夏季环流的重要环流特征之一.在计算区域东北部夏季环流以反气旋环流系统为主;在计算区域东南部夏季环流以气旋系统为主;南海夏季环流分布,明显出现西部强化特征. 相似文献
7.
Gian Mario BeltramiRiccardo Briganti 《Ocean Engineering》2003,30(7):833-853
This paper presents a simple and straightforward method for carrying out the direct numerical solution of the eigenvalue problem associated to the homogeneous linear shallow-water equations expressed using orthogonal curvilinear coordinates, when ‘adiabatic’ boundary conditions apply. These equations, together with the boundary conditions, define a self-adjoint problem in the continuum. The method presented here, which is thought for calculating the 2-D theoretical gravity modes of both natural and artificial basins, relies on a change of basis of the dependent variable vector. This preliminary transformation makes it, in fact, possible to formulate two different numerical approaches which guarantee the self-adjoint property of the discrete form of the system consisting of the governing equations and the boundary conditions. The method is tested using a square and a fully circular domain, both of which allow comparisons with well-known analytical and numerical solutions. Discretizing the physical domain of a fully circular basin by a cylindrical coordinate grid makes it possible to show the actual efficiency of the method in calculating the theoretical gravity modes of basins discretized by a boundary-following coordinate grid which allows laterally variable resolution. 相似文献
8.
Hemming A. Schäffer 《Coastal Engineering》2009,56(5-6):517-533
The transformation of irrotational surface gravity waves in an inviscid fluid can be studied by time stepping the kinematic and dynamic surface boundary conditions. This requires a closure providing the normal surface particle velocity in terms of the surface velocity potential or its tangential derivative. A convolution integral giving this closure as an explicit expression is derived for linear 1D waves over a mildly sloping bottom. The model has exact linear dispersion and shoaling properties. A discrete numerical model is developed for a spatially staggered uniform grid. The model involves a spatial derivative which is discretized by an arbitrary-order finite-difference scheme. Error control is attained by solving the discrete dispersion relation a priori and model results make a perfect match to this prediction. A procedure is developed by which the computational effort is minimized for a specific physical problem while adapting the numerical parameters under the constraint of a predefined tolerance of damping and dispersion error. Two computational examples show that accurate irregular-wave transformation on the kilometre scale can be computed in seconds. Thus, the method makes up a highly efficient basis for a forthcoming extension that includes nonlinearity at arbitrary order. The relation to Boussinesq equations, mild-slope wave equations, boundary integral equations and spectral methods is briefly discussed. 相似文献
9.
为了能够精确地拟合天然水域不规则的水陆边界,并能保持有限差分数值计算在矩形网格中进行,本文利用Thompson的数值网格生成技术,在长江口南支七丫口至横沙河段设计了一个椭圆型边界拟合坐标系,并在这个坐标系下建立了平面二维潮流数学模型。通过数值求解此数学模型,实现了计算区域内现状潮流的一个数值模拟,模拟结果表明,模式的设计是成功的。 相似文献
10.
《Ocean Modelling》1999,1(2-4):71-80
Ocean general circulation models (OGCMs) which represent the governing equations on a finite difference grid require shorter time steps with increasing resolution. Thus, until now, in the absence of filtering, the time step length has been determined by the smallest grid spacing within the model domain. Here we present a method for reducing the time step length (and increasing the number of time steps taken) at selected points in the grid, so as to minimise the computational cost of integrating the OGCM, whilst achieving numerical stability throughout the model domain without filtering. This variable time stepping method can be used to overcome numerical constraints associated with the convergence of longitude–latitude grids at the poles, and also to allow efficient integration of model domains with variable resolution. Examples of the computational saving are given. 相似文献
11.
Parametric vertical coordinate formulation for multiscale, Boussinesq, and non-Boussinesq ocean modeling 总被引：2，自引：1，他引：2
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. 相似文献
12.
A 2‐time‐level finite difference atmospheric general circulation model based on the semi‐Lagrangian advection of pseudo potential vorticity (which becomes potential vorticity in that part of the domain where the hybrid vertical coordinate becomes isentropic) has been formulated. At low levels, the hybrid vertical coordinate is terrain following. The problem of isentropic potential vorticity possibly becoming ill‐defined in the regions of planetary boundary layer is thus circumvented. The divergence equation is a companion to the (pseudo) potential vorticity equation and the model is thus called a PV‐D model. Many features of a previously developed shallow water PV‐D model are carried over: a modification of the PV equation needed to give computational stability of long Rossby waves; a semi‐Lagrangian semi‐implicit treatment of both the linear and the nonlinear terms; the use of an unstaggered grid in the horizontal; the use of a nonlinear multigrid technique to solve the nonlinear implicit equations. A linear numerical stability analysis of the model's gravity–inertia waves indicates that the potential temperature needs to be separated into horizontal mean and perturbation parts. This allows an implicit treatment of the vertical advection associated with the mean in the thermodynamic equation. Numerical experiments with developing baroclinic waves have been carried out and give realistic results. 相似文献
13.
14.
A BFG model for calculation of tidal current and diffusion of pollutants in nearshore areas 总被引：4，自引：0，他引：4
ＡＢＦＧｍｏｄｅｌｆｏｒｃａｌｃｕｌａｔｉｏｎｏｆｔｉｄａｌｃｕｒｒｅｎｔａｎｄｄｉｆｆｕｓｉｏｎｏｆｐｏｌｌｕｔａｎｔｓｉｎｎｅａｒｓｈｏｒｅａｒｅａｓ￥ＳｈｉＦｅｎｇｙａｎ；ａｎｄＺｈｅｎｇ；Ｌｉａｎｙｕａｎ（ＳｔａｔｅＫｅｙＬａｂｏｒａｔｏｒｙ... 相似文献
15.
Hastings F.D. Schneider J.B. Broschat S.L. Thorsos E.I. 《Oceanic Engineering, IEEE Journal of》2001,26(1):94-101
A finite-difference time-domain (FDTD) method for scattering by one-dimensional, rough fluid-fluid interfaces is presented, modifications to the traditional FDTD algorithm are implemented which yield greater accuracy at lower computational cost. These modifications include use of a conformal technique, in which the grid conforms locally to the interface, and a correction for the numerical dispersion inherent to the FDTD algorithm, Numerical results are presented for fluid-fluid cases modeling water-sediment interfaces. Two different roughness spectra, the single-scale Gaussian roughness spectrum and a multiscale modified power-law spectrum, are used. The Gaussian results are calculated as a function of the dimensionless parameters kh and kl, where k is the wavenumber in water, h is the rms surface height, and l is the surface correlation length. For the modified power-law spectrum, statistical parameters consistent with an insonification frequency of 7.5 kHz are used. Results are compared with those obtained using an integral equation technique both for scattering from single-surface realizations and for Monte Carlo averages of scattering from an ensemble of surface realizations. Scattering strengths are calculated as a function of scattering angle for an incident angle of 70° (20° grazing). The results agree well over all scattering angles for the cases examined 相似文献
16.
Dmitry S. Dukhovskoy Steven L. Morey Paul J. Martin James J. O&#x;Brien Cortis Cooper 《Ocean Modelling》2009,28(4):250-265
Recent observations over the Sigsbee Escarpment in the Gulf of Mexico have revealed extremely energetic deep currents (near 1 m s−1), which are trapped along the escarpment. Both scientific interest and engineering needs demand dynamical understanding of these extreme events, and can benefit from a numerical model designed to complement observational and theoretical investigations in this region of complicated topography. The primary objective of this study is to develop a modeling methodology capable of simulating these physical processes and apply the model to the Sigsbee Escarpment region. The very steep slope of the Sigsbee Escarpment (0.05–0.1) limits the application of ocean models with traditional terrain-following (sigma) vertical coordinates, which may represent the very complicated topography in the region adequately, can result in large truncation errors during calculation of the horizontal pressure gradient. A new vertical coordinate system, termed a vanishing quasi-sigma coordinate, is implemented in the Navy Coastal Ocean Model for application to the Sigsbee Escarpment region. Vertical coordinate surfaces for this grid have noticeably gentler slopes than a traditional sigma grid, while still following the terrain near the ocean bottom. The new vertical grid is tested with a suite of numerical experiments and compared to a classical sigma-layer model. The numerical error is substantially reduced in the model with the new vertical grid. A one-year, realistic, numerical simulation is performed to simulate strong, deep currents over the Escarpment using a very-high-resolution nested modeling approach. The model results are analyzed to demonstrate that the deep-ocean currents in the simulation replicate the prominent dynamical features of the observed intense currents in the region. 相似文献
17.
基于推板造波理论和摇板造波理论,在Open FOAM平台上采用重叠网格技术建立黏性数值波浪水槽,并使用一种结合SIMPLE算法和PISO算法的PIMPLE算法对数值模型进行求解。利用开发的数值模型通过数值收敛性测试和网格独立性测试分别重点研究了时间步长、库朗数和网格尺寸对数值精度和计算效率的影响。并对比研究了此数值模型分别嵌入层流模型和湍流模型的计算精度和计算效率。实现的规则波和二阶有限振幅波与理论结果和试验结果吻合,验证了此黏性数值波浪水槽的造波和主动消波功能。基于二维数值波浪水槽,进一步研究了三维数值造波,数值计算结果与理论结果吻合良好。研究结果不仅验证了重叠网格在二维和三维两相流体域中求解运动物体与流场交互的可靠性和正确性,而且为使用此黏性数值波浪水槽解决更复杂的海洋工程问题提供了依据。 相似文献
18.
A new adaptive Cartesian-grid for the CIP (constrained interpolation profile) method is proposed and applied to two-dimensional numerical simulations of violent free-surface flows. The CCUP (CIP combined and unified procedure) method is employed and combined with this adaptive Cartesian-grid for robust and efficient computation. This adaptive grid is capable of tracking regions where flows vary violently, and a much finer grid is then concentrated automatically on these regions to adapt to the violent changing of the flow. Unlike the abacus-like Soroban grid which is an adaptive meshless grid with complicated algorithms and inefficiency of evaluation of frequently computed spatial derivatives, the present approach not only simplifies computational algorithm but also enhances efficiency of frequently-computed spatial derivatives. It is also different from most of the remeshing schemes that no additional CPU-time for the value-mapping from the old grid to the new grid is taken in this adaptive grid system provided that the advection velocity is interpolated, since the value-mapping process is accomplished simultaneously within the advection process. To validate the accuracy and efficiency of this newly-proposed CFD model, several two-dimensional benchmark problems are performed, and the results are compared with experimental measurements and other published numerical results. Numerical simulations show that the proposed numerical model is robust, accurate, and efficient for strongly nonlinear free-surface flows. 相似文献
19.
The paper develops and analyzes two fully nonlinear boundary conditions that incorporate the motion of the shoreline in nonlinear time domain nearshore models. A moving shoreline essentially means the computational domain is changing with the solution of the flow. The problem is solved in two steps. The first is to establish an equation that determines the motion of the shoreline based on the local momentum balance. The second is to develop and implement into a shoreline model the capability of accommodating a changing computational domain. The two models represent two different ways of addressing this step: one is to track the position of the shoreline in a fixed grid by establishing a special shoreline point which generally is not a fixed grid point. The second is by a coordinate transformation that maps the changing domain onto a fixed domain and solves the basic equations in the mapped domain. The two shoreline conditions are tested against three known solution for nonlinear shoreline motion. Two are the 1-D solutions to the nonlinear shallow water (NSW) equations by Carrier and Greenspan [J. Fluid Mech. 4 (1958) 97], one representing the response to a transient change in the offshore water level, the other the motion due to a periodic standing wave, both on slopes steep enough to allow full reflection. The third is the 2-D horizontal (2DH) computational solution by Zelt [Coast. Eng. 15 (1991) 205] for the run-up of a solitary wave on a cusped beach. In all cases, both models are shown to behave well and give high accuracy results for suitably chosen grid and time spacings. 相似文献
20.
An alternative ‘compressed’ coordinate system for representing finite difference grids in ocean models is presented and compared with the traditional Cartesian method. The alternative method represents any arbitrary three-dimensional domain as a one-dimensional vector. Advantages of using this method include the exclusion of all land ‘dry’ points from the computational grid, leading to savings in the memory used and increased execution speed where domains have a low ratio of wet to total computational cells. The Cartesian and compressed models are applied to a variety of domains to assess the performance of each in terms of memory use and computational speed. 相似文献