共查询到20条相似文献，搜索用时 734 毫秒
1.
为了研究污染物对流扩散方程中不同时空格式的适用性,针对对流扩散方程的一维﹑二维和三维3种情况,分别建立了预报-校正的有限差分数值模型。在时间步进格式上分别采用了Crank-Nicolson格式或混合4阶Adams-Bashforth-Moulton格式,对对流项分别采用2阶精度或4阶精度,对扩散项采用了2阶精度。利用建立的数值模型求解了经典的污染物浓度场对流扩散,通过数值解与解析解的比较讨论了不同时空格式对数值模型计算结果的影响。结果表明:对空间一次导数采用4阶精度可以避免采用2阶精度带来的误差。采用混合4阶Adams-Bashforth-Moulton格式或Crank-Nicolson格式数值计算结果均与解析解吻合程度较好,但对于数组为[40,40,40]的三维对流扩散问题,前者比后者省时20.7%。 相似文献
2.
W. E. Rogers J. M. Kaihatu H. A. H. Petit N. Booij L. H. Holthuijsen 《Ocean Engineering》2002,29(11):1357-1390
The numerical schemes for the geographic propagation of random, short-crested, wind-generated waves in third-generation wave models are either unconditionally stable or only conditionally stable. Having an unconditionally stable scheme gives greater freedom in choosing the time step (for given space steps). The third-generation wave model SWAN (“Simulated WAves Nearshore”, Booij et al., 1999) has been implemented with this type of scheme. This model uses a first order, upwind, implicit numerical scheme for geographic propagation. The scheme can be employed for both stationary (typically small scale) and nonstationary (i.e. time-stepping) computations. Though robust, this first order scheme is very diffusive. This degrades the accuracy of the model in a number of situations, including most model applications at larger scales. The authors reduce the diffusiveness of the model by replacing the existing numerical scheme with two alternative higher order schemes, a scheme that is intended for stationary, small-scale computations, and a scheme that is most appropriate for nonstationary computations. Examples representative of both large-scale and small-scale applications are presented. The alternative schemes are shown to be much less diffusive than the original scheme while retaining the implicit character of the particular SWAN set-up. The additional computational burden of the stationary alternative scheme is negligible, and the expense of the nonstationary alternative scheme is comparable to those used by other third generation wave models. To further accommodate large-scale applications of SWAN, the model is reformulated in terms of spherical coordinates rather than the original Cartesian coordinates. Thus the modified model can calculate wave energy propagation accurately and efficiently at any scale varying from laboratory dimensions (spatial scale O(10 m) with resolution O(0.1 m)), to near-shore coastal dimension (spatial scale O(10 km) with resolution O(100 m)) to oceanic dimensions (spatial scale O(10 000 km) with resolution O(100 km). 相似文献
3.
The typical equation for bed level change in sediment transport in river, estuary and near shore systems is based on conservation of sediment mass. It is generally a nonlinear conservation equation for bed level. The physics here are similar to shallow water wave equations and gas dynamics equation which will develop shock waves in many circumstances. Many state-of-art morphological models use classical lower order Lax–Wendroff or modified Lax–Wendroff schemes for morphology which are not very stable for long time sediment transport processes simulation. Filtering or artificial diffusion are often added to achieve stability. In this paper, several shock capturing schemes are discussed for simulating bed level change with different accuracy and stability behaviors. The conclusion is in favor of a fifth order Euler-WENO scheme which is introduced to sediment transport simulations here over other schemes. The Euler-WENO scheme is shown to have significant advantages over schemes with artificial viscosity and filtering processes, hence is highly recommended especially for phase-resolving sediment transport models. 相似文献
4.
A novel computational approach is introduced for the efficient computation of equilibrium solutions of seasonally forced ocean biogeochemical models. The essential idea is to formulate the problem as a large system of nonlinear algebraic equations to be solved with a class of methods known as matrix-free Newton–Krylov (MFNK). MFNK is a combination of Newton-type methods for superlinearly convergent solution of nonlinear equations, and Krylov subspace methods for solving the Newton correction equations. The basic link between the two methods is the Jacobian-vector product, which may be probed approximately without forming and storing the elements of the true Jacobian. To render this approach practical for global models with O(106) degrees of freedom, a flexible preconditioning strategy is developed. The result is an essentially “black-box” numerical scheme than can be applied to most existing biogeochemical models. The method is illustrated by applying it to find the equilibrium solutions of two realistic biogeochemical problems. Compared with the conventional approach of direct time integration, the preconditioned-MFNK scheme is shown to be roughly two orders of magnitude more efficient. Several potential refinements of the basic algorithm that may yield further performance gains are discussed. The numerical scheme described here addresses a fundamental challenge to using ocean biogeochemical models more effectively. 相似文献
5.
《Coastal Engineering》2005,52(2):139-149
An improved parabolic water wave transformation model is developed based on generalized [1/1] Padé approximation. For forward scattered waves, the parabolic equation is solved using a marching scheme. The values of wave angles are calculated after the solution of each line; so that better [1/1] generalized Padé approximation is performed. The nonlinear effects are included using a modified dispersion equation. The model is easy to use and performs very well for complex bathymetry. The model is tested for cases of wave angles up to 70°. The numerical results show that for large wave angles, the new parabolic model is better than all the existing parabolic models based on rational approximation. 相似文献
6.
Parabolic Approximation of the Weakly Nonlinear Mild Slope Equation with Bottom Friction 总被引：2，自引：0，他引：2
This paper presents a refined parabolic approximation model of the mild slope equation to simu-late the combination of water wave refraction and diffraction in the large coastal region.The bottom frictionand weakly nonlinear term are included in the model.The difference equation is established with the Crank-Nicolson scheme.The numerical test shows that some numerical prediction results will be inaccurate in com-plicated topography without considering weak nonlinearity;the bottom friction will make wave height damp-ing and it can not be neglected for calculation of wave field in large areas. 相似文献
7.
This paper presents the numerical solution of a new nonlinear mild-slope equation governing waves with different frequency components propagating in a region of varying water depth. There are two new nonlinear equations. The linear part of the equations is the mild-slope equation, and one of the models has the same non-linearity as the Boussinesq equations. The new equations are directly applicable to the problems of nonlinear wave-wave interactions over variable depth. The equations are first simplified with the parabolic approximation, and then solved numerically with a finite difference method. The Crank-Nicolson method is used to discretize the models. The numerical models are applied to a set of published experimental cases, which are nonlinear combined refraction-diffraction with generation of higher harmonic waves. Comparison of the results shows that the present models generally predict the measurements better than other nonlinear numerical models which have been applied to the data set. 相似文献
8.
ADCIRC, a finite element circulation model for shelves, coasts and estuaries, will be used for variational data assimilation. The nonlinear Euler–Lagrange (EL) problem will be solved using the iterated indirect representer algorithm. This algorithm makes such large, nonlinear but functionally smooth optimization problems feasible by iterating on linear approximations of the nonlinear problem (Picard iterations) and by making preconditioned searches in the “data subspace” at each iterate. Before solving the nonlinear EL using such Picard iterations, it essential that the iteration scheme be carefully examined within the framework of the nonassimilative or forward problem.The purpose of this paper is (1) to detail a Picard iteration procedure for ADCIRC, including the problematic bottom friction term; (2) to examine the ability of the iteration scheme to recover the nonlinear forward solution from deficient background fields; and (3) to present a study of different interpolation methods for reducing the memory/disk requirements of the iteration scheme. The iteration scheme is shown to be quite robust in its ability to recover the nonlinear solution from a variety of deficient background fields. A new cubic Hermitian interpolation method is shown to be a more effective alternative to standard linear interpolation for reducing memory/disk requirements, especially for high frequency overtides. 相似文献
9.
Dynamic responses of structures due to earthquake excitation are the important problems in engineering, thus, the information concerned is plenty. However, most of the literature is relating to the discrete methods, particularly to the finite element method (FEM), and the one relating to the method combining both the “continuous” and “discrete” models is rare. The objective of this paper is to provide some information in this respect. First, the analytical solution for the natural frequencies and normal mode shapes of a “continuous” tower, without contacting water (or “dry” tower), carrying an eccentric tip mass possessing rotary inertia is determined. Next, the partial differential equation of motion for the forced vibration of the tower, contacting water (or “wet” tower), subjected to support excitation is transformed into a matrix equation by using the last natural frequencies and normal modes shape of the freely vibrating dry tower. Finally, the numerical integration method is used to solve the matrix equation to yield the seismic response of the wet tower. In theory, the mode superposition method is correct only if the total number of modes considered approaches infinity, however, numerical results of this paper reveal that superposition of only the lowest six modes will yield excellent results to be very close to the corresponding ones obtained from the conventional FEM. For this reason, the CPU time required by the presented approach is less than 5% of that required by the conventional FEM. 相似文献
10.
研究并行算法解决应用并行计算机完成规模尽可能大的偏微分方程的数值求解问题。利用Hopf-Cole变换,将一维非线性Burgers方程转化为线性扩散方程,基于第二类Saul’yev型非对称格式和Crank-Nicolson格式对扩散方程进行差分离散,建立解Burgers方程的交替分段并行差分格式,并讨论该方法的稳定性,给出了数值算例。此算法把剖分节点分成若干组,在每组上构造能够独立求解的差分方程,因此具有并行本性,适合在高性能多处理器的并行计算机上使用。数值试验的结果表明此方法是有效的,且有较高的精度。 相似文献
11.
12.
A numerical solution is developed to investigate the generation and propagation of small-amplitude water waves in a semi-infinite rectangular wave basin. The three-dimensional wave field is produced by the prescribed “snake-like” motion of an array of segmented wave generators located along the wall at one end of the tank. The solution technique is based on the boundary element approach and uses an appropriate three-dimensional Green function which explicitly satisfies the tank-wall boundary conditions. The Green function and its derivatives which appear in the integral equation formulation can be shown to be slowly convergent when the source and field points are in close proximity. Therefore, when computing the velocity potentials on the wave generators, the source points are chosen outside the fluid domain, thereby ensuring the rapid convergence of these functions and rendering the integral equations non-singular. Numerical results are shown which illustrate the influence of the various wavemaker and basin parameters on the generated wave field. Finally, the complete wave field produced by the diffraction of oblique waves by a vertical circular cylinder in a basin is presented. 相似文献
13.
We point out one problem of the grid advection schemes when used in wave models in coastal areas. The deficiency of the schemes is investigated by means of the ‘third' generation WAM wave model, in which the wave energy is advected by a first order upwind scheme. Two similar, alternative modifications of this scheme are analyzed, the second of which is shown to solve most of the problems encountered with advection along the co-ordinate axes. 相似文献
14.
Realistic representation of sea ice in ocean models involves the use of a non-linear free-surface, a real freshwater flux and observance of requisite conservation laws. We show here that these properties can be achieved in practice through use of a rescaled vertical coordinate “z*” in z-coordinate models that allows one to follow undulations in the free-surface under sea ice loading. In particular, the adoption of “z*” avoids the difficult issue of vanishing levels under thick ice.Details of the implementation within MITgcm are provided. A high resolution global ocean sea ice simulation illustrates the robustness of the z* formulation and reveals a source of oceanic variability associated with sea ice dynamics and ice-loading effects. The use of the z* coordinate allows one to achieve perfect conservation of fresh water, heat and salt, as shown in extended integration of coupled ocean sea ice atmospheric model. 相似文献
15.
A limitation of many global climate models with explicit finite-difference numerics is the timestep restriction caused by the decrease in cell size associated with the convergence of meridians near the poles. To keep the longitudinal width of model cells as uniform as possible, we apply a “reduced” grid to a three-dimensional primitive equation ocean-climate model. With this grid the number of cells in the longitudinal direction is reduced at high latitudes. The grid consists of subgrids which interact at interfaces along their northern and southern boundaries, where the resolution changes by a factor of three. We extend the finite-difference techniques to these interfaces, focusing on the conservation required to perform long time integrations, while preserving the staggered spatial arrangement of variables and the numerics used on subgrids. The common alternative used to reduce the timestep restriction caused by the spherical grid is the filtering of high-frequency modes from the high-latitude solution. The reduced grid allows an increased timestep while eliminating the need for filtering and reduces execution time per model step by roughly 20%. We implement the reduced grid model for parallel computer architectures with two-dimensional domain decomposition and message passing, with speedup results similar to those of the original model. We present results of model runs showing small effects on the solution and sizable improvements to the execution time. 相似文献
16.
Multidimensional Positive Definite Advection Transport Algorithm (MPDATA) and Monotonic Upstream Schemes for Conservation Laws (MUSCL) are compared and evaluated in a biogeochemical-like tracer model. Based on the MUSCL scheme, an adapted scheme for σ-coordinate models (MUSCL_VD) is presented. These schemes are first implemented in idealized experiments configured for a closed channel with and without a seamount, and conservation is confirmed. MUSCL_VD and MUSCL generate the same results in the flat bottom case, while MUSCL_VD is the only scheme that avoids unphysical values in the seamount experiments. When implemented in a three-dimensional model of the Gulf of Maine (GoM) with realistic topography and real-time forcing, the modeled distributions of an inert tracer concentration are sensitive to the advection schemes. The experiment with MUSCL shows persistent overshoots. MPDATA and MUSCL_VD produce similar spatial distributions and seasonal cycles with MUSCL_VD being slightly more diffusive in the horizontal and less so in the vertical. These two schemes are about equally efficient in computation. 相似文献
17.
《Coastal Engineering》2005,52(5):391-407
A numerical solver is presented of the modified time-independent mild-slope equation, which incorporates energy dissipation. Using a second-order parabolic approximation, the following external boundary conditions are modelled: open and fully transmitting to both incoming and outgoing waves; partially reflecting, and; fully absorbing. Discretisation of the governing equation and boundary conditions is by means of a second-order accurate central difference scheme. The resulting sparse-banded matrix is solved using an inexpensive banded solver with Gaussian elimination. The numerical predictions are in excellent agreement with the analytical solution for the interaction of non-breaking waves with an array of vertical surface-piercing circular cylinders on a horizontal bed. Results are compared with those for the same array on various seabed topographies. The model is robust and can be used for wave propagation in complex geometries. It has fewer restrictions associated with wave obliqueness at boundaries than traditional models based on the mild-slope equation. 相似文献
18.
David Prandle 《Progress in Oceanography》1991,27(3-4)
Since 1984 the OSCR HF Radar system has been used in over 50 deployments to measure near-shore surface currents for both scientific and engineering applications. The enhanced scope, resolution and accuracy of these measurements have yielded new insights into the tidal, wind and density driven dynamics of the near-shore zone.Tidal current ellipses obtained from these radar measurements have been shown to be in good aggrement with values calculated by numerical models both for the predominant constituents and also for higher harmonics. Coherent patterns of wind-forced currents ahve been determined with strong evidence of a “slab-like” surface response. In one deployment, with offshore winds blowing over relatively deep water, this “slab” rotated clockwise at near-inertial frequency. Strong (up to 20cm s−1), persistent surface residual currents are commonly observed, these are almost certainly generated by (small) horizontal density gradients. These observed surface residuals provide ideal data for rigorous testing of 3-D numerical models.With a threatened rise in sea level, HF Radar is well-suited for observing the expected changes in the dynamics of near-shore regions. Continuing development of these radar systems offers exciting prospects of remote sensing of both surface waves and currents. Future applications may extend beyond the near-shore region to measurements along the shelf-edge, in oceanic gyres and for “beach-processes”. 相似文献
19.
解四阶拟线性波动方程的一类二阶差分格式 总被引：2，自引：1，他引：2
建立了解一类四阶拟线性耗散、色散波动方程初边值问题的Crank-Nicolson差分方法，并结合外推的技巧，给出了1个线性化方法；证明了差分解的存在唯一性；用能量估计的方法证明了此格式的二阶收敛性和无条件稳定性；给出了一些数值结果。 相似文献
20.
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. 相似文献