共查询到20条相似文献，搜索用时 0 毫秒
1.
刘儒勋 《应用数学和力学(英文版)》1992,13(3):233-240
In this paper a new approach for designing upwind type schemes-the characterizing-integral method and its applied skills are introduced, The method is simple, convenient and eff ective. And the method isn't only limited to conservation laws unlike other methods and may be easily extended to multi-dimension problems. Furthermore, the numerical dissipation of the method can be flexibly regulated, so that it is especially suitable for solving various discontinuity problems. The paper shows us how to use this approach to simulate deformation and breaking of a nonlinear shallow water wave on a gentle slope, and to compute two-dimensional dam failure problem. 相似文献
2.
Liu Ru-xun 《应用数学和力学(英文版)》1992,13(3):233-240
In this paper a new approach for designing upwind type schemes-the characterizing-integral method and its applied skills are introduced. The method is simple, convenient and eff ective. And the method isn ’t only limited to conservation laws unlike other methods and maybe easily extended to multi-dimension problems. Furthermore, the numerical dissipation of the method can be flexibly regulated, so that it is especially suitable for solving various discontinuity problems.The paper shows us now to use this approach to simulate deformation and breaking of a nonlinear shallow water wave on a gentle slope, and to compute two-dimensional dam failure problem. 相似文献
3.
In this paper, a new numerical model, based on a set of non‐linear shallow water equations is developed for the simulation of the formation and evolution of tidal bore in the Hangzhou Gulf and Qiantangjiang river of China. The numerical method and boundary conditions are described in detail. The method is validated against analytical solutions and experimental data. Simulation of the actual tidal bore in Hangzhou Gulf and its propagation in the Qiantangjiang river are performed. Numerical results show that this proposed method is effective for the prediction of tidal bore and current flow at the entrance of Qiantangjiang river. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
4.
In this paper, a new vector‐filtering criterion for dynamic modes selection is proposed that is able to extract dynamically relevant flow features from dynamic mode decomposition of time‐resolved experimental or numerical data. We employ a novel modes selection criterion in parallel with the classic selection based on modes amplitudes, in order to analyze which of these procedures better highlight the coherent structures of the flow dynamics. Numerical tests are performed on two distinct problems. The efficiency of the proposed criterion is proved in retaining the most influential modes and reducing the size of the dynamic mode decomposition model. By applying the proposed filtering mode technique, the flow reconstruction error is shown to be significantly reduced. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
5.
A review of wetting and drying (WD) algorithms used by contemporary numerical models based on the shallow water equations is presented. The numerical models reviewed employ WD algorithms that fall into four general frameworks: (1) Specifying a thin film of fluid over the entire domain; (2) checking if an element or node is wet, dry or potentially one of the two, and subsequently adding or removing elements from the computational domain; (3) linearly extrapolating the fluid depth onto a dry element and its nodes from nearby wet elements and computing the velocities; and (4) allowing the water surface to extend below the topographic ground surface. This review presents the benefits and drawbacks in terms of accuracy, robustness, computational efficiency, and conservation properties. The WD algorithms also tend to be highly tailored to the numerical model they serve and therefore difficult to generalize. Furthermore, the lack of temporally and spatially defined validation data has hampered comparisons of the models in terms of their ability to simulate WD over real domains. A short discussion of this topic is included in the conclusion. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
6.
The approximation of reduced linear evolution operator (propagator) via dynamic mode decomposition (DMD) is addressed for both linear and nonlinear events. The 2D unsteady supersonic underexpanded jet, impinging the flat plate in nonlinear oscillating mode, is used as the first test problem for both modes. Large memory savings for the propagator approximation are demonstrated. Corresponding prospects for the estimation of receptivity and singular vectors are discussed. The shallow water equations are used as the second large‐scale test problem. Excellent results are obtained for the proposed optimized DMD method of the shallow water equations when compared with recent POD‐based/discrete empirical interpolation‐based model reduction results in the literature. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
7.
Depth‐averaged non‐hydrostatic extension for shallow water equations with quadratic vertical pressure profile: equivalence to Boussinesq‐type equations 
下载免费PDF全文
下载免费PDF全文 We reformulate the depth‐averaged non‐hydrostatic extension for shallow water equations to show equivalence with well‐known Boussinesq‐type equations. For this purpose, we introduce two scalars representing the vertical profile of the non‐hydrostatic pressure. A specific quadratic vertical profile yields equivalence to the Serre equations, for which only one scalar in the traditional equation system needs to be modified. Equivalence can also be demonstrated with other Boussinesq‐type equations from the literature when considering variable depth, but then the non‐hydrostatic extension involves mixed space–time derivatives. In case of constant bathymetries, the non‐hydrostatic extension is another way to circumvent mixed space–time derivatives arising in Boussinesq‐type equations. On the other hand, we show that there is no equivalence when using the traditionally assumed linear vertical pressure profile. Linear dispersion and asymptotic analysis as well as numerical test cases show the advantages of the quadratic compared with the linear vertical non‐hydrostatic pressure profile in the depth‐averaged non‐hydrostatic extension for shallow water equations. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
8.
9.
Time domain simulation of the interaction between offshore structures and irregular waves in shallow water becomes a focus due to significant increase of liquefied natural gas (LNG) terminals. To obtain the time series of irregular waves in shallow water, a numerical wave tank is developed by using the meshless method for simulation of 2D nonlinear irregular waves propagating from deep water to shallow water. Using the fundamental solution of Laplace equation as the radial basis function (RBF) and locating the source points outside the computational domain, the problem of water wave propagation is solved by collocation of boundary points. In order to improve the computation stability, both the incident wave elevation and velocity potential are applied to the wave generation. A sponge damping layer combined with the Sommerfeld radiation condition is used on the radiation boundary. The present model is applied to simulate the propagation of regular and irregular waves. The numerical results are validated by analytical solutions and experimental data and good agreements are observed. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
10.
In this paper, we introduce a finite‐volume kinetic BGK scheme and its applications to the study of roll and solitary waves. The current scheme is based on the numerical solution of the gas‐kinetic Bhatnagar–Gross–Krook model in the flux evaluation across each cell interface. An intrinsic connection between the BGK model and time‐dependent, non‐linear, non‐homogeneous shallow‐water equations enables us to solve shallow‐water equations automatically with our kinetic scheme. The analytical solution, experimental measurements, and numerical calculations for problems associated with roll‐waves down an inclined open channel and solitary waves incident on a sloped beach are also presented. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
11.
12.
Long waves generated by a moving atmospheric pressure distribution, associated with a storm, in coastal region are investigated numerically. For simplicity the moving atmospheric pressure is assumed to be moving only in the alongshore direction and the beach slope is assumed to be a constant in the on-offshore direction. By solving the linear shallow water equations we obtain numerical solutions for a wide range of physical parameters, including storm size (2a), storm speed (U), and beach slope (α). Based on the numerical results, it is determined that edge wave packets are generated if the storm speed is equal to or greater than the critical velocity, Ucr, which is defined as the phase speed of the fundamental edge wave mode whose wavelength is scaled by the width of the storm size. The length and the location of the positively moving edge wave packet is roughly Ut/2 ≤ y ≤ Ut, where y is in the alongshore direction and t is the time. Once the edge wave packet is generated, the wavelength is the same as that of the fundamental edge wave mode corresponding to the storm speed and is independent of the storm size, which can, however, affect the wave amplitude. When the storm speed is less than the critical velocity, the primary surface signature is a depression directly correlated to the atmospheric pressure distribution. 相似文献
13.
A multidimensional discretisation of the shallow water equations governing unsteady free-surface flow is proposed. The method, based on a residual distribution discretisation, relies on a characteristic eigenvector decomposition of each cell residual, and the use of appropriate distribution schemes. For uncoupled equations, multidimensional convection schemes on compact stencils are used, while for coupled equations, either system distribution schemes such as the Lax–Wendroff scheme or scalar schemes may be used. For steady subcritical flows, the equations can be partially diagonalised into a purely convective equation of hyperbolic nature, and a set of coupled equations of elliptic nature. The multidimensional discretisation, which is second-order-accurate at steady state, is shown to be superior to the standard Lax–Wendroff discretisation. For steady supercritical flows, the equations can be fully diagonalised into a set of convective equations corresponding to the steady state characteristics. Discontinuities such as hydraulic jumps, are captured in a sharp and non-oscillatory way. For unsteady flows, the characteristic equations remain coupled. An appropriate treatment of the coupling terms allows the discretisation of these equations at the scalar level. Although presently only first-order-accurate in space and time, the classical dam-break problem demonstrates the validity of the approach. © 1998 John Wiley & Sons, Ltd. 相似文献
14.
K. N. Gavrilova 《Journal of Applied Mechanics and Technical Physics》2004,45(1):37-44
A hyperbolic model of a shallow water flow is considered with allowance for nonlinear and dispersion effects. The structure of traveling waves above a flat bottom is studied. Stability of small disturbances of a homogeneous flow and development of instability of a nonstationary flow above an inclined bottom are analyzed. 相似文献
15.
This paper presents a parametric finite‐difference scheme concerning the numerical solution of the one‐dimensional Boussinesq‐type set of equations, as they were introduced by Peregrine (J. Fluid Mech. 1967; 27 (4)) in the case of waves relatively long with small amplitudes in water of varying depth. The proposed method, which can be considered as a generalization of the Crank‐Nickolson method, aims to investigate alternative approaches in order to improve the accuracy of analogous methods known from bibliography. The resulting linear finite‐difference scheme, which is analysed for stability using the Fourier method, has been applied successfully to a problem used by Beji and Battjes (Coastal Eng. 1994; 23 : 1–16), giving numerical results which are in good agreement with the corresponding results given by MIKE 21 BW (User Guide. In: MIKE 21, Wave Modelling, User Guide. 2002; 271–392) developed by DHI Software. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
16.
This paper introduces a sparse matrix discrete interpolation method to effectively compute matrix approximations in the reduced order modeling framework. The sparse algorithm developed herein relies on the discrete empirical interpolation method and uses only samples of the nonzero entries of the matrix series. The proposed approach can approximate very large matrices, unlike the current matrix discrete empirical interpolation method, which is limited by its large computational memory requirements. The empirical interpolation indices obtained by the sparse algorithm slightly differ from the ones computed by the matrix discrete empirical interpolation method as a consequence of the singular vectors round‐off errors introduced by the economy or full singular value decomposition (SVD) algorithms when applied to the full matrix snapshots. When appropriately padded with zeros, the economy SVD factorization of the nonzero elements of the snapshots matrix is a valid economy SVD for the full snapshots matrix. Numerical experiments are performed with the 1D Burgers and 2D shallow water equations test problems where the quadratic reduced nonlinearities are computed via tensorial calculus. The sparse matrix approximation strategy is compared against five existing methods for computing reduced Jacobians: (i) matrix discrete empirical interpolation method, (ii) discrete empirical interpolation method, (iii) tensorial calculus, (iv) full Jacobian projection onto the reduced basis subspace, and (v) directional derivatives of the model along the reduced basis functions. The sparse matrix method outperforms all other algorithms. The use of traditional matrix discrete empirical interpolation method is not possible for very large dimensions because of its excessive memory requirements. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
17.
The present paper makes use of a wave equation formulation of the primitive shallow water equations to simulate one-dimensional free surface flow. A numerical formulation of the boundary element method is then developed to solve the wave continuity equation using a time-dependent fundamental solution, while an explicit finite difference scheme is used to derive velocities from the primitive momentum equation. One-dimensional free surface flows in open channels are treated and the results compared with analytical and numerical solutions. © 1997 John Wiley & Sons, Ltd. 相似文献
18.
Wave equation models currently discretize the generalized wave continuity equation with a three‐time‐level scheme centered at k and the momentum equation with a two‐time‐level scheme centered at k+1/2; non‐linear terms are evaluated explicitly. However in highly non‐linear applications, the algorithm becomes unstable at even moderate Courant numbers. This paper examines an implicit treatment of the non‐linear terms using an iterative time‐marching algorithm. Depending on the domain, results from one‐dimensional experiments show up to a tenfold increase in stability and temporal accuracy. The sensitivity of stability to variations in the G‐parameter (a numerical weighting parameter in the generalized wave continuity equation) was examined; results show that the greatest increase in stability occurs with G/τ=2–50. In the one‐dimensional experiments, three different types of node spacing techniques—constant, variable, and LTEA (Localized Truncation Error Analysis)—were examined; stability is positively correlated to the uniformity of the node spacing. Lastly, a scaling analysis demonstrates that the magnitudes of the non‐linear terms are positively correlated to those that most influence stability, particularly the term containing the G‐parameter. It is evident that the new algorithm improves stability and temporal accuracy in a cost‐effective manner. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
19.
20.
The objective of this paper is to present a methodology of using a two‐step split‐operator approach for solving the shallow water flow equations in terms of an orthogonal curvilinear co‐ordinate system. This approach is in fact one kind of the so‐called fractional step method that has been popularly used for computations of dynamic flow. By following that the momentum equations are decomposed into two portions, the computation procedure involves two steps. The first step (dispersion step) is to compute the provisional velocity in the momentum equation without the pressure gradient. The second step (propagation step) is to correct the provisional velocity by considering a divergence‐free velocity field, including the effect of the pressure gradient. This newly proposed method, other than the conventional split‐operator methods, such as the projection method, considers the effects of pressure gradient and bed friction in the second step. The advantage of this treatment is that it increases flexibility, efficiency and applicability of numerical simulation for various hydraulic problems. Four cases, including back‐water flow, reverse flow, circular basin flow and unsteady flow, have been demonstrated to show the accuracy and practical application of the method. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献