首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 328 毫秒
1.
1.IN~DUCTIONTurbiditycurrentisoneclassofflowsnameddensitycurrentorgravitycurrent(therHunterRouse(Yih(1980)),whichdePictstheintmsionofheaVyfluidintoalighterone.Usually,thedensitydifferencebetWeentWonuidisrelativelysmallandmixingacrosstheimerfaceoccurs.ThedrivingforceofdensitycurrentsisnotdensitydifferenceitselfbutthedifferenceinspeCmcweights.Turbiditycurrentisnamedwhenthedensitydifferenceisespeciallycausedbysuspendedfinesedimentparticles.Sincesediment-ladenflowcaninteraCtwiththelowerbou…  相似文献   

2.
本文将DRP/opt MacCormack有限差分格式用于模拟二维各向异性介质中的地震波传播.DRP/opt MacCormack是一种同位网格下的差分格式,避免了传统的交错网格在计算各向异性问题时由于变量插值而导致的误差.而且相对于低阶同位网格差分格式,它具有低色散、低耗散的优点.此格式将中心差分算子分成前向和后向两个空间单边差分,然后在4-6步Runge-Kutta时间积分中使用单边差分组合.在具有垂直对称轴的横向各向同性(VTI)模型下,通过对比DRP/opt MacCormack有限差分和谱元方法的模拟结果,验证了前者具有很高的精度和稳定性.由于实际地质条件下TI介质的对称轴通常是倾斜的(TTI),本文在二维三分量框架下模拟TTI介质中的地震波场.结果显示横波分裂和切平面/反平面运动耦合的特征.数值实验表明DRP/opt MacCormack是一种有效的研究各向异性介质中地震波传播规律的差分格式.  相似文献   

3.
In this paper, we deduced the corresponding first-order velocity–stress equation for curvilinear coordinates from the first-order velocity–stress equation based on the modified Biot/squirt model for a two-dimensional two-phase medium. The equations are then numerically solved by an optimized high-order non-staggered finite difference scheme, that is, the dispersion relation preserving/optimization MacCormack scheme. To implement undulating free-surface topography, we derive an analytical relationship between the derivatives of the particle velocity components and use the compact finite-difference scheme plus a traction-image method. In the undulating free surface and the undulating subsurface interface of two-phase medium, the complex reflected wave and transmitted wave can be clearly recognized in the numerical simulation results. The simulation results show that the curvilinear-grid finite-difference method, which uses a body-conforming grid to describe the undulating surface, can accurately reduce the numerical scattering effect of seismic wave propagation caused by the use of ladder-shaped grid to fit the surfaces when undulating topography is present in a two-phase isotropic medium.  相似文献   

4.
求解弹性波方程的辛RKN格式   总被引:2,自引:2,他引:0       下载免费PDF全文
将弹性波方程变换至Hamilton体系,构造适用于弹性波模拟的高效显式二阶辛Runge-Kutta-Nystrm(RKN)格式,运用根数理论得到此格式的阶条件方程组.通过给定系数的限定条件,得到方程的对称解.为了使时间离散误差达到极小,提出数值频率与真实频率比较,通过Taylor展开,得到关于辛系数的限定方程,求解方程组得到最小频散辛RKN格式.对比分析时间演进方程的稳定性,得到使库朗数达到极大值的限定方程,求解方程组得到最稳定辛RKN格式.发现此两种格式为同一格式.新得到的辛RKN格式不依赖于空间离散方法,为了对比的需要,选取有限差分法进行空间离散.在频散、稳定性分析中,与常见辛格式对比,从理论上分析了本文提出的格式在数值频散压制、稳定性提升等方面的优势,数值实验进一步证实了理论分析的正确性.  相似文献   

5.
We derive a governing second-order acoustic wave equation in the time domain with a perfectly matched layer absorbing boundary condition for general inhomogeneous media. Besides, a new scheme to solve the perfectly matched layer equation for absorbing reflections from the model boundaries based on the rapid expansion method is proposed. The suggested scheme can be easily applied to a wide class of wave equations and numerical methods for seismic modelling. The absorbing boundary condition method is formulated based on the split perfectly matched layer method and we employ the rapid expansion method to solve the derived new perfectly matched layer equation. The use of the rapid expansion method allows us to extrapolate wavefields with a time step larger than the ones commonly used by traditional finite-difference schemes in a stable way and free of dispersion noise. Furthermore, in order to demonstrate the efficiency and applicability of the proposed perfectly matched layer scheme, numerical modelling examples are also presented. The numerical results obtained with the put forward perfectly matched layer scheme are compared with results from traditional attenuation absorbing boundary conditions and enlarged models as well. The analysis of the numerical results indicates that the proposed perfectly matched layer scheme is significantly effective and more efficient in absorbing spurious reflections from the model boundaries.  相似文献   

6.
1 INTRODUCTION Numerical computation provides an easily extended and user-friendly environment with computer aided programming for the simulation of pollutants in river systems. In the most of water quality assessment and monitoring problems during water pollution control and environmental impact assessment studies of river systems, mathematical modeling has been playing a key role for the last two decades. The majority of existing water quality models are the mechanistic and are based o…  相似文献   

7.
In recent years there has been a growing interest in using Godunov-type methods for atmospheric flow problems. Godunov's unique approach to numerical modeling of fluid flow is characterized by introducing physical reasoning in the development of the numerical scheme (van Leer, 1999). The construction of the scheme itself is based upon the physical phenomenon described by the equation sets. These finite volume discretizations are conservative and have the ability to resolve regions of steep gradients accurately, thus avoiding dispersion errors in the solution. Positivity of scalars (an important factor when considering the transport of microphysical quantities) is also guaranteed by applying the total variation diminishing condition appropriately. This paper describes the implementation of a Godunov-type finite volume scheme based on unstructured adaptive grids for simulating flows on the meso-, micro- and urban-scales. The Harten-Lax-van Leer-Contact (HLLC) approximate Riemann solver used to calculate the Godunov fluxes is described in detail. The higher-order spatial accuracy is achieved via gradient reconstruction techniques after van Leer and the total variation diminishing condition is enforced with the aid of slope-limiters. A multi-stage explicit Runge-Kutta time marching scheme is used for maintaining higher-order accuracy in time. The scheme is conservative and exhibits minimal numerical dispersion and diffusion. The subgrid scale diffusion in the model is parameterized via the Smagorinsky-Lilly turbulence closure. The scheme uses a non-staggered mesh arrangement of variables (all quantities are cell-centered) and requires no explicit filtering for stability. A comparison with exact solutions shows that the scheme can resolve the different types of wave structures admitted by the atmospheric flow equation set. A qualitative evaluation for an idealized test case of convection in a neutral atmosphere is also presented. The scheme was able to simulate the onset of Kelvin-Helmholtz type instability and shows promise in simulating atmospheric flows characterized by sharp gradients without using explicit filtering for numerical stability.  相似文献   

8.
The flow patterns in confluence channel and the simulation of confluence flow are more complex than that in straight channel. Additional terms in the momentum equations, i.e. dissipation terms, denoting the impact of turbulence, and dispersion terms, denoting the vertical non‐uniformity of velocity, show great impacts on the accuracy of numerical simulations. The dissipation terms, i.e. the product of eddy viscosity coefficient and velocity gradient, are much larger than those of the flow in straight channel. In this study, the zero equation model and the depth‐averaged k‐ε model are used to analyse the impact of eddy viscosity. Meanwhile, the dispersion terms in the momentum equation, depending on the vertical non‐uniformity of velocity, are usually neglected in routine simulation. With the use of detailed experimental data for verification, this study presents the distribution of parameters of vertical non‐uniformity and the intimated connection between non‐uniformity parameters and accuracy of numerical simulations of confluence flow with depth‐averaged models. The results present that simulation accuracy of confluence flow is very sensitive to the turbulence modes, which cannot be handled by normal, simple turbulence model. On the contrary, the impact of dispersion terms is both flow‐condition‐dependent and place‐dependent, and such impact is negligible when secondary circulation is weak. The results indicate the key elements in modelling confluence flow and are helpful for selecting suitable numerical model and solving engineering problems encountered in confluence channel. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

9.
如何有效压制数值频散是有限差分正演模拟研究中的关键问题之一.近年来,许多学者对二阶声波方程的差分算子开展了大量的优化工作,在压制频散方面取得不错的效果.一阶压强-速度方程广泛用于研究地震波在地下变密度模型中传播规律,目前针对一阶方程的优化工作大多只是在空间差分算子上展开.本文在前人研究的基础上,推导出一阶声波方程中压强场与偏振速度场之间的解析关系,据此在传统交错网格基础上给出一种高精度的显式时间递推格式,该递推格式将时间差分与空间差分算子结合在一起,并采用共轭梯度法得到精确时间递推匹配系数,实现时空差分算子的同时优化.在编程实现算法的基础上,通过频散分析与三个典型模型测试表明:本文方法能够较为有效地压制时间频散与空间频散,提高数值计算精度;同时对复杂模型也有很好适用性.  相似文献   

10.
间断有限元方法(Discontinuous Galerkin method,简称DGM)在求解地震波动方程时具有低数值频散、网格剖分灵活等优点,因此,为适应数值模拟对模拟精度和复杂地质结构的要求,本文提出一种新的加权Runge-Kutta间断有限元(weighted Runge-Kutta discontinuous Galerkin,简称WRKDG)方法,用于求解三维D′Alembert介质中声波方程.本文不仅详细推导了其数值格式,特别地,根据常微分方程理论给出了满足数值稳定性条件的一般经验公式,并首次对该方法的数值频散和耗散进行了深入分析,且考虑了耗散参数对结果的影响.同时,我们也对该方法进行了精度测试,并分析了3D情形下WRKDG方法的并行加速比,结果表明3D WRKDG方法具有良好的并行性.最后,我们给出了包含均匀模型、非规则几何模型以及非均匀Marmousi模型在内的数值模拟算例.结果表明,该方法不仅计算准确,能与解析解很好地吻合,且能有效模拟包含球体在内的非规则模型及非均匀Marmousi模型中的衰减声波波场.数值模拟实验进一步验证了WRKDG方法在求解三维D′Alembert介质中声波方程时的正确性和有效性,并获得了对这种强衰减介质中波传播特征的规律性新认识.  相似文献   

11.
The success of transient storage (TS) modeling for natural streams depends, in part, on the ability to describe the dispersion process accurately. Evidence based on stream tracer data shows that solute transport processes often do not follow the classical second-order dispersion model (e.g., early breakthrough and faster than Fickian travel times were observed). While models based on space-fractional dispersion are a promising alternative, different definitions of fractional derivatives exist in the literature. Unlike integer-order derivatives, fractional derivatives represent convolutions of concentration with long-range spatial correlation and numerical approximations can produce dense matrices. Therefore issues of both accuracy and computational efficiency need to be examined to successfully identify model parameters for natural streams. In this paper, we first compare the performance of several numerical approaches for solving the space-fractional dispersion equation. We examine three different numerical approaches to approximate the space-fractional derivatives including: (a) a fully-implicit scheme based on the shifted Grünwald–Letnikov (GL) approximation (b) a three-point implicit representation based on the GL formula and (c) a three-point implicit scheme based on mass conservation and the Caputo definition of the fractional derivative. We then use an operator-splitting technique to evaluate a TS model based on space-fractional dispersion (the FSTS model) and test the model against analytical solutions and stream tracer data. A sequence acceleration method (Richardson extrapolation) significantly improves the performance of all schemes examined. Results indicate that the fully-implicit GL method with Richardson extrapolation produces the most accurate solutions while the three-point implicit GL scheme has a stringent time-step restriction to produce acceptable solutions. The three-point implicit scheme based on the Caputo derivative produces accurate solutions in a fraction of the time taken by the fully-implicit GL method and represents the best trade-off between accuracy and computational efficiency for practical applications. The scheme is suitable for parameter estimation and is used to successfully describe tracer data in a natural stream.  相似文献   

12.
Solute transport is usually modeled by the advection-dispersion-reaction equation. In the standard approach, mechanical dispersion is a tensor with principal directions parallel and perpendicular to the flow vector. Since realistic scenarios include nonuniform and unsteady flow fields, the governing equation has full tensor mechanical dispersion. When conventional grid-based numerical methods are used, approximation of the cross terms arising from the off-diagonal terms cause nonphysical solution with oscillations. As an example, for the common scenario of contaminant input into a domain with zero initial concentration, the cross-dispersion terms can result in negative concentrations that can wreak havoc in reactive transport applications. To address this issue, we use the well-known flux-corrected-transport (FCT) technique for a standard finite volume method. Although FCT has most often been used to eliminate oscillations resulting from discretization of the advection term for explicit time stepping, we show that it can be adapted for full-tensor dispersion and implicit time stepping. Unlike other approaches based on new discretization techniques (e.g., mimetic finite difference, nonlinear finite volume), FCT has the advantage of being flexible and widely applicable. Implementation of FCT requires solving an additional system of equations at each time step, using a modified “low order” matrix and a modified right-hand-side vector. To demonstrate the flexibility of FCT, we have modified the well-known and widely used groundwater solute transport simulator, MT3DMS. We apply the new simulator, MT3DMS-FCT, to several benchmark problems that suffer from negative concentrations when using MT3DMS. The new results are mass conservative and strictly nonnegative.  相似文献   

13.
1 INTRODUCTIONFor many hydraulic engineering problems, the analysis of flow and bed level variations in openchannels is a fundamental prerequisite. forcal methOds fOr alluvial rivers are well develoPednowadays as far as onediInensional descriPtions are concemed. A cOmPrhensive analysis of Ihe wellknown models is Presented by Habersack(l998). HOwever, for a number of Problems such as channelwidening, flow pattem close to sPuds and etc. a more deailed knowledge of the bed level behavio…  相似文献   

14.
Analytical solutions for the water flow and solute transport equations in the unsaturated zone are presented. We use the Broadbridge and White nonlinear model to solve the Richards’ equation for vertical flow under a constant infiltration rate. Then we extend the water flow solution and develop an exact parametric solution for the advection-dispersion equation. The method of characteristics is adopted to determine the location of a solute front in the unsaturated zone. The dispersion component is incorporated into the final solution using a singular perturbation method. The formulation of the analytical solutions is simple, and a complete solution is generated without resorting to computationally demanding numerical schemes. Indeed, the simple analytical solutions can be used as tools to verify the accuracy of numerical models of water flow and solute transport. Comparison with a finite-element numerical solution indicates that a good match for the predicted water content is achieved when the mesh grid is one-fourth the capillary length scale of the porous medium. However, when numerically solving the solute transport equation at this level of discretization, numerical dispersion and spatial oscillations were significant.  相似文献   

15.
I. Haltas 《水文研究》2012,26(22):3448-3458
Recognizing the spatial heterogeneity of hydraulic parameters, many researchers have studied the solute transport by both groundwater and channel flow in a stochastic framework. One of the methodologies used to up‐scale the stochastic solute transport equation, from a point‐location scale to a grid scale, is the cumulant expansion method combined with the calculus for the time‐ordered exponential and the calculus for the Lie operator. When the point‐location scale transport equation is scaled up to the grid scale, using the cumulant expansion method, a new dispersion coefficient emerges in the dispersive term of the solute transport equation in addition to the molecular dispersion coefficient. This velocity driven dispersion is called ‘macrodispersion’. The macrodispersion coefficient is the integral function of the time‐ordered covariance of the random velocity field. The integral is calculated over a Lagrangian trajectory of the flow. The Lagrangian trajectory depends on the following: (i) the spatial origin of the particle; (ii) the time when the macrodispersion is calculated; and (iii) the mean velocity field along the trajectory itself. The Lagrangian trajectory is a recursive function of time because the location of the particle along the trajectory at a particular time depends on the location of the particle at the previous time. This recursive functional form of the Lagrangian trajectory makes the calculation of the macrodispersion coefficient difficult. Especially for the unsteady, spatially non‐stationary, non‐uniform flow field, the macrodispersion coefficient is a highly complex expression and, so far, calculated using numerical methods in the discrete domains. Here, an analytical method was introduced to calculate the macrodispersion coefficient in the discrete domain for the unsteady and steady, spatially non‐stationary flow cases accurately and efficiently. This study can fill the gap between the theory of the ensemble averaged solute transport model and its numerical implementations. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

16.
In this paper a fuzzy dynamic wave routing model (FDWRM) for unsteady flow simulation in open channels is presented. The continuity equation of the dynamic wave routing model is preserved in its original form while the momentum equation is replaced by a fuzzy rule based model which is developed on the principle that during unsteady flow the disturbances in the form of discontinuities in the gradient of the physical parameters will propagate along the characteristics with a velocity equal to that of velocity of the shallow water wave. The model gets rid off the assumptions associated with the momentum equation by replacing it with the fuzzy rule based model. It overcomes the necessity of calculating friction slope (Sf) in flow routing and hence the associated uncertainties are eliminated. The robustness of the fuzzy rule based model enables the FDWRM to march the solution even in regions where the aforementioned assumptions are violated. Also the model can be used for flow routing in curved channels. When the model is applied to hypothetical flood routing problems in a river it is observed that the results are comparable to those of an implicit numerical model (INM) which solves the dynamic wave equations using an implicit numerical scheme. The model is also applied to a real case of flow routing in a field canal. The results match well with the measured data and the model performs better than the INM. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

17.
本文以基于改进BISQ模型的二维双相各向同性介质一阶速度-应力方程为基础,推导出了曲线坐标系下对应的方程,然后采用低频散、低耗散的同位网格MacCormack有限差分法来离散方程,并采用紧致的单边MacCormack差分格式结合牵引力镜像法来施加自由地表边界条件,实现了地震波场数值模拟.曲线网格有限差分法采用贴体网格来描述自由表面,地表的网格线紧贴地形,避免了台阶近似造成的数值散射.数值模拟结果表明,在双相介质起伏自由地表和分界面处,各类波型复杂的反射透射规律可以清晰展现,曲线网格有限差分法可以精确地解决地震波在含起伏地表的双相各向同性介质中的传播问题.  相似文献   

18.
Herrera P  Valocchi A 《Ground water》2006,44(6):803-813
The transport of contaminants in aquifers is usually represented by a convection-dispersion equation. There are several well-known problems of oscillation and artificial dispersion that affect the numerical solution of this equation. For example, several studies have shown that standard treatment of the cross-dispersion terms always leads to a negative concentration. It is also well known that the numerical solution of the convective term is affected by spurious oscillations or substantial numerical dispersion. These difficulties are especially significant for solute transport in nonuniform flow in heterogeneous aquifers. For the case of coupled reactive-transport models, even small negative concentration values can become amplified through nonlinear reaction source/sink terms and thus result in physically erroneous and unstable results. This paper includes a brief discussion about how nonpositive concentrations arise from numerical solution of the convection and cross-dispersion terms. We demonstrate the effectiveness of directional splitting with one-dimensional flux limiters for the convection term. Also, a new numerical scheme for the dispersion term that preserves positivity is presented. The results of the proposed convection scheme and the solution given by the new method to compute dispersion are compared with standard numerical methods as used in MT3DMS.  相似文献   

19.
In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.  相似文献   

20.
地震波场数值模拟在地球物理勘探和地震学中具有重要的支撑作用.本文将组合型紧致差分格式用于声波和弹性波方程的数值模拟中.根据泰勒级数展开和声波方程,建立了位移场时间四阶离散格式,并将组合型紧致差分格式用于位移场空间导数的求取,然后对该差分格式进行了精度分析、误差分析、频散分析和稳定性分析.理论研究结果表明:①该差分格式为时间四阶、空间六阶精度,与常规七点六阶中心差分和五点六阶紧致差分相比,具有更小的截断误差和更高的模拟精度;②每个波长仅需要5.6个采样点,且满足稳定性条件的库郎数为0.792,可以使用粗网格和较大时间步长进行计算.所以该方法具有占用内存少、计算效率高和低数值频散等优势.最后,本文进行了二维各向同性完全弹性介质的声波和弹性波方程的数值模拟,实验结果表明本文提出的方法具有更高的计算精度,能够大幅度的节约计算量和内存需求,对于三维大尺度模型问题具有更好的适应性.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号