共查询到19条相似文献，搜索用时 140 毫秒
1.
一、引言圣维南方程组通常被用来描述河渠中的非恒定的流动,对应于动量方程的不同形式与连续方程的组合,则可得到动力波,惯性波、扩散波和运动波等不同洪水波模型。动量方程中保留所有的项,即得通常所说的动力波,本文称为全方程情形,亦指除惯性 相似文献
2.
基于改进BISQ模型的地震波场数值模拟 总被引：1，自引：0，他引：1
Biot流动和喷射流动是含流体多孔隙介质中流体流动的二种重要力学机制,Dvorkin和Nur提出了同时包含Biot流动和喷射流动力学机制的统一的BISQ(Biot-Squirt)模型。由于BISQ模型的流体压力表达式十分复杂,MamadouS.D和ErwinA又在引入Squirt机制的同时提出了不含特征喷射流动长度的改进BISQ模型。这里基于改进BISQ模型,运用伪谱法进行了波场模拟,证实了快P波、慢P波、SV波和SH波的存在。通过与Biot模型和BISQ模型的比较,改进BISQ模型计算简单,便于描述,正演模拟结果也与其它二个模型保持一致,是一种可行的新BISQ模型。 相似文献
3.
概述了模拟海岸、河口的浅水波浪数值模型研究现状、存在的问题以及用能量平衡方程预报海浪的发展历史。介绍了基于当代最新波浪理论研究成果的第三代浅水波浪数值模型SWAN模型,对模型的适用性、数值特性、功能及局限性进行了阐述。介绍动谱平衡方程数学模型、方程离散要求、边界条件的处理和源项(包括能量输入、损耗及波与波之间非线性相互作用)的处理方法,重点介绍三相波非线性相互作用。模拟海安湾有效波高、波周期场,并分析波与波之间非线性相互作用对波浪要素预报的影响,最后对SWAN模型的应用前景和研究趋势进行了展望。 相似文献
4.
二阶弹性波动方程高精度交错网格波场分离数值模拟 总被引：4，自引：0，他引：4
给出了一种等价的二阶弹性波动方程,以解决弹性波场中完全弹性波动方程不能完全分离耦合的纵、横波波场问题.应用高阶交错网格有限差分法求解该波动方程,并使用通量校正技术(FCT)进一步压制频散,采用均匀介质模型和层状介质模型进行波场分离数值试验,精确得到了混合波场、完全分离的纯纵波及纯横波波场.数值结果分析表明,本文方法在均匀介质情况下准确可靠,在分离后的纯纵、横波波场中可观察到较为丰富的能量转换信息,这对认识复杂弹性波的传播规律及弹性波理论具有重要意义. 相似文献
5.
采用复合单元法建立了模拟裂隙多孔介质变饱和流动的数值模型。该模型具有以下特点：裂隙不需要离散成特定单元，而是根据几何位置插入到孔隙基质单元中形成复合单元；在复合单元中，分别建立裂隙流和孔隙基质流的计算方程，二者通过裂隙?基质界面产生联系并整合成复合单元方程；复合单元方程具有和常规有限单元方程相同的格式，因此，可以使用常规有限单元方程的求解技术。采用欠松弛迭代、集中质量矩阵以及自适应时步调节等技术，开发了裂隙多孔介质变饱和流动计算程序。通过模拟一维干土入渗和复杂裂隙含水层内的流动问题，验证了该模型的合理性和适用性。模拟结果为进一步认识非饱和裂隙含水层地下水流动特性提供了理论依据。 相似文献
6.
7.
裂缝诱导双相HTI介质模型及其弹性波传播方程 总被引：1，自引：0，他引：1
将Biot双相介质理论与Gurevich裂缝各向异性理论相结合,建立了能够同时考虑实际裂缝性储层孔隙性和各向异性的裂缝诱导双相HTI介质模型。从本构方程、动力学方程和动力学达西定律出发,推导出了裂缝诱导双相HTI介质中弹性波传播的一阶速度-应力方程,并针对方程的刚性问题,给出了利用显式二阶时间积分法数值求解该方程时所需要满足的稳定性条件。该方程能够定量地给出双相HTI介质的波场特征与裂缝参数、背景孔隙介质参数之间的关系,描述弹性波在这种介质中的传播机理。 相似文献
8.
9.
建立了基于高阶Boussinesq水波方程的一维波浪破碎数值模型。基于一组具有二阶完全非线性特征的Boussinesq水波方程,建立了交错网格下的高精度差分格式,推导了适用于该组方程的永形波解析解,其和松弛造波技术相结合实现了数值波浪水槽中(强)非线性波浪的无反射入射。通过模拟封闭容器内水体晃动问题对数值格式进行了验证,通过模拟孤立波在斜坡海岸上的浅化过程说明了将方程从弱非线性发展到完全非线性的必要性。采用涡粘方法处理波浪破碎,利用物理模型实验数据,分析了模型中各波浪破碎参数对数值结果的影响并对参数进行了率定。应用该模型对规则波在斜坡海岸上的传播、变浅以及破碎过程进行了数值模拟研究,数值结果同实验数据吻合良好,验证了模型的有效性。 相似文献
10.
TTI介质qP波数值模拟方法因为考虑了倾角因素,可以比VTI介质qP波数值模拟方法更加准确地描述各向异性介质中地震波场的传播规律。文中用拟声波方程对TTI介质中的地震波场进行了高阶有限差分数值模拟,在改进衰减函数分布方式后,通过坐标变换,利用改进的完全匹配层(perfectly matched layer,PML)边界控制方程对波场边界进行吸收处理,取得了良好的效果;然后分析了拟声波方程数值模拟中的稳定性问题,并对波场中的伪横波进行压制。通过对不同模型的数值模拟,验证了文中使用的TTI介质拟声波波动方程的稳定性以及所采用的PML边界控制方程的可靠性和适用性。 相似文献
11.
基于常Q模型的解耦分数阶拉普拉斯算子粘滞波动方程,可以分开模拟振幅衰减和相位错动。但该方程拉普拉斯算子的阶数是随空间变化的,因此数值求解存在一定困难。这里基于截断的泰勒展开,经过一系列近似,推导出拉普拉斯算子的阶数与空间无关的解耦分数阶粘滞弹性波动方程。采用中心差分计算时间导数,使用交错网格伪谱法计算空间导数。数值算例表明,新的方程在处理非均匀介质时具有精度高,计算简便的优点。 相似文献
12.
平面波的传播问题通常可以归结为一维波动方程的定解问题。在非均匀介质中,即使简单的一维波动方程也需要借助于数值方法获得近似解。3层5点古典差分格式是计算偏微分方程一种常用算法,作为一种显式迭代格式,需要满足稳定性条件 ,其中 为波速, 为空间采样间隔, 为时间采样间隔。当 时, ,古典差分格式达到临界稳定状态。在这种情况下,平面波在 时间内的传播距离恰好等于空间采样间隔,差分格式真实地反映了平面波的传播原理,因而可以得到一维波动方程的精确解。但是,由于在非均匀介质中存在不连续的波阻抗界面,此方法不适于计算非均匀介质的波场。为了将临界稳定情况下的古典差分格式推广应用至非均匀层状介质,提出了一种能够处理波阻抗界面的有限差分格式,并应用傅里叶分析法得到其稳定性条件。模型算例验证了此算法的正确性。 相似文献
13.
In this paper, we study the properties of approximate solutions to a doubly nonlinear and degenerate diffusion equation, known
in the literature as the diffusive wave approximation of the shallow water equations (DSW), using a numerical approach based
on the Galerkin finite element method. This equation arises in shallow water flow models when special assumptions are used
to simplify the shallow water equations and contains as particular cases the porous medium equation and the p-Laplacian. Diverse
numerical schemes have been implemented to approximately solve the DSW equation and have been successfully applied as suitable
models to simulate overland flow and water flow in vegetated areas such as wetlands; yet, no formal mathematical analysis
has been carried out in order to study the properties of approximate solutions. In this study, we propose a numerical approach
as a means to understand some properties of solutions to the DSW equation and, thus, to provide conditions for which the use
of the DSW equation may be inappropriate from both the physical and the mathematical points of view, within the context of
shallow water modeling. For analysis purposes, we propose a numerical method based on the Galerkin method and we obtain a
priori error estimates between the approximate solutions and weak solutions to the DSW equation under physically consistent
assumptions. We also present some numerical experiments that provide relevant information about the accuracy of the proposed
numerical method to solve the DSW equation and the applicability of the DSW equation as a model to simulate observed quantities
in an experimental setting. 相似文献
14.
在基于波动方程的有限差分数值模拟中,会不可避免地出现数值频散（也称网格发散）问题。数值频散问题通常会给数值模拟的结果造成严重影响,因此在数值模拟中,应尽量设法消除这种现象。这里在前人的基础上,在基于各向同性介质的弹性波方程的数值模拟过程中,通过引入通量校正方法来解决数值频散问题。由数值模拟结果表明,该方法可以有效地消除数值频散现象,大大改善数值模拟的结果。 相似文献
15.
基于波作用量守恒方程建立了波流共存场中多向随机波浪传播变形数学模型,模型中考虑了波浪绕射的影响和水流引起的波浪弥散多普勒效应,应用包含水流和地形影响的激破波模式计算波浪破碎的能量耗散,采用一阶上迎风有限差分格式离散控制方程。分别计算了有无近岸流情况下单向和多向随机波浪的波高分布,考虑水流影响的数值计算结果与物理模型实验数据吻合良好,比较分析表明,所建立的数学模型能够复演由于离岸流引起的波高增大,可用于波流共存场多向随机波浪传播变形的模拟和预报。 相似文献
16.
17.
18.
19.
Acoustic imaging and sensor modeling are processes that require repeated solution of the acoustic wave equation. Solution of the wave equation can be computationally expensive and memory intensive for large simulation domains. One scheme for speeding up solution of the wave equation is the operator-based upscaling method. The algorithm proceeds in two steps. First, the wave equation is solved for fine grid unknowns internal to coarse blocks assuming the coarse blocks do not need to communicate with neighboring blocks in parallel. Second, these fine grid solutions are used to form a new problem which is solved on the coarse grid. Accurate and efficient wave propagation schemes also must avoid artificial reflections off of the computational domain edges. One popular method for preventing artificial reflections is the nearly perfectly matched layer (NPML) method. In this paper, we discuss applying NPML to operator upscaling for the wave equation. We show that although we only apply NPML to the first step of this two step algorithm (directly affecting the fine grid unknowns only), we still see a significant reduction of reflections back into the domain. We describe three numerical experiments (one homogeneous medium experiment and two heterogeneous media examples) in which we validate that the solution of the wave equation exponentially decays in the NPML regions. Numerical experiments of acoustic wave propagation in two dimensions with a reasonable absorbing layer thickness resulted in a maximum pressure reflection of 3–8%. While the coarse grid acceleration is not explicitly damped in our algorithm, the tight coupling between the two steps of the algorithm results in only 0.1–1% of acceleration reflecting back into the computational domain. 相似文献