共查询到18条相似文献，搜索用时 156 毫秒
1.
波动数值模拟中透射边界的稳定实现 总被引：21，自引：0，他引：21
从波动能量在计算区内累积增大的观点出发,通过简单的一维弹性波模型,系统地阐明了在近场波动数值模拟中透射边界两类数值失稳--"高频振荡"和"零频飘移"的机理:前者源于对波动数值模拟无意义的高频波动在人工边界上的放大和波在有限计算区内多次反射产生的反复放大;后者则源于透射边界允许零频和接近零频的分量不断进入计算区. 由此提出了稳定实现透射边界的完整方案包括两项简单措施:第一,在全部计算区内按文中建议的方法注入小阻尼,以消除高频振荡;第二, 给出一种具有明确物理意义的消除零频飘移的算子算法. 最后,提供了三维波源问题和散射问题的详细数值试验结果. 相似文献
2.
《地震工程与工程振动》2017,(1)
在波动数值模拟中,瑞利阻尼可近似描述介质耗散特性,且可用于抑制人工边界引发的高频和零频失稳,但瑞利阻尼对波动的影响尚未清晰认识。针对集中质量有限元模拟的一维波动,利用傅里叶模态分析了有阻尼离散网格中波动的性质。理论分析表明alpha阻尼必然使得波数为零及邻近的波动对应为非行进波,其使行波衰减一致,而beta阻尼不会导致波数为零及邻近的非行进波,其使行波衰减随着波数增大而增大。数值实验验证了上述结论。本文研究结果为进一步推进瑞利阻尼在波动数值模拟中的应用提供一定的理论依据。 相似文献
3.
时域高阶双渐近透射边界能够同时模拟层状介质中行波和快衰波的传播,具有很高的计算精度和计算效率.本文将高阶双渐近透射边界推广应用到多层层状地基系统弹性波传播问题的模拟,采用广义特征值分解分析该透射边界的数值稳定性,通过移谱法消除导致数值不稳定的虚假模态.将高阶双渐近透射边界以超单元的形式直接嵌入到近场有限元方程,建立了有限元-高阶双渐近透射边界时域耦合分析模型,并将其应用到重力坝-层状地基动力相互作用分析.数值算例分析结果表明,该时域耦合分析模型具有很高的精度和计算效率,适用于实际重力坝工程的地震响应分析. 相似文献
4.
高阶双渐近时域透射边界能够同时模拟行波和快衰波的传播,并且能够在全频范围内迅速逼近准确解,具有优良的收敛性能和计算效率.本文将动水压力波高阶双渐近透射边界直接嵌入到近场有限元方程中,建立了大坝-库水动力相互作用的直接耦合分析模型.该模型的整体控制方程保留了近场有限元方程系数矩阵对称稀疏的优势,可以方便地利用现有的通用有限元求解器求解.基于有限元开源软件框架体系OpenSees（Open System for Earthquake Engineering Simulation）,编程实现了直接耦合分析模型,并将其应用于二维重力坝、三维拱坝与库水动力相互作用分析.数值算例表明,该直接耦合分析模型具有很高的精度和计算效率. 相似文献
5.
本文采用广义反射系数法推导了水平层状各向异性地层中电磁场的积分解析解,并利用快速汉克尔变换技术实现了三维感应仪器测井响应的快速计算.三维感应测井响应与地层水平电导率、垂直电导率和井斜角及仪器方位角同时有关,单一分量的测井曲线不能满足资料解释的需要.通过对仪器测量分量响应特征的考察,本文提出了一种基于组合量测井曲线的资料直观解释方法.数值模拟显示,交叉分量相关组合量可准确划分地层纵向边界,并可直观识别各向异性层;与单独分量相比,主分量相关组合量提高了纵向分辨率、减弱了与地层电导率参数的非线性关系. 相似文献
6.
7.
8.
9.
含煤层地质环境下地震波场的数值模拟 总被引：3，自引：2，他引：1
本文对含低速煤层地质环境下弹性波场多波多分量地震资料进行了二维数值模拟研究,对人工边界反射进行了有效处理,频散效应得到了有效的压制,对几种不同激发与观测排列方式下的弹性波资料进行了模型计算与分析。 相似文献
10.
CSAMT三维单分量有限元正演 总被引：1，自引：0，他引：1
CSAMT的工作装置和部分探测目标体的三维特性,决定了对CSAMT方法进行三维研究的必要性.在对三维三分量CSAMT方法有限元分析初探的基础上,应用了边界场值不为零的第一类边界条件,并将三维三分量的有限元分析退化为三维单分量来研究.研究结果表明:通过应用边界场值已知的第一类边界条件,缩小了研究范围,提高了计算精度;通过减少研究的分量,使计算速度大大提高;模型的模拟结果很好地说明了电磁场在地下传播中所具有的穿透性和体积效应.以上结果表明本文所用的边界条件有效,三维单分量的研究使计算时间大大减少,从而使本文发展的三维单分量电磁场有限元正演应用于电性介质接近于均匀的三维反演成为可能. 相似文献
11.
The novel doubly asymptotic multi-directional transmitting boundary combines the advantages of the doubly asymptotic and multi-directional formulations. It is rigorous for the low-frequency limit (static case) and the high-frequency limit in the wave-propagation direction perpendicular to the artificial boundary and at all the preselected angles corresponding to various apparent velocities. It is implemented straightforwardly in the finite-element method. The doubly asymptotic multi-directional transmitting boundary is temporally local and spatially either global or local depending on the implementation of the static-stiffness matrix. It leads to much higher accuracy than other transmitting boundaries with the same finite-element mesh. 相似文献
12.
The scaled boundary finite‐element method is extended to simulate time‐harmonic responses of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. The unbounded domains and the elasticity matrices are transformed to the scaled boundary coordinates. The scaled boundary finite‐element equation in displacement amplitudes are derived directly from the governing equations of elastodynamics. To enforce the radiation condition at infinity, an asymptotic expansion of the dynamic‐stiffness matrix for high frequency is developed. The dynamic‐stiffness matrix at lower frequency is obtained by numerical integration of ordinary differential equations. Only the boundary is discretized yielding a reduction of the spatial dimension by one. No fundamental solution is required. Material anisotropy is modelled without additional efforts. Examples of two‐ and three‐dimensional non‐homogeneous isotropic and transversely isotropic unbounded domains are presented. The results demonstrate the accuracy and simplicity of the scaled boundary finite‐element method. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
13.
The assumption of spatial repetition is commonly made when producing bedform scale models of the hyporheic zone. Two popular solute transport codes, MT3DMS and PHT3D, do not currently provide the necessary boundary condition required to simulate spatial periodicity in hyporheic zone transport problems. In this study, we develop a spatially periodic boundary (SPB) for solutes that is compatible with a SPB that was previously developed for MODFLOW to simulate the flow component of spatially periodic problems. The approach is ideal for simulating groundwater flow and transport patterns under repeating surface features, such as ripples or dunes on the bottom of a lake or stream. The appropriate block‐centered finite‐difference approach to implement the boundary is presented and the necessary source code modifications are discussed. The performance of the solute SPB, operating in conjunction with the groundwater flow SPB, is explored through comparison of a multi‐bedform hyporheic‐zone model with a single bedform variant. The new boundary conditions perform well in situations where both dispersive effects and lateral seepage flux in the underflow regime beneath the hyporheic zone are minimal. 相似文献
14.
Steven D. London 《地球物理与天体物理流体动力学》2013,107(6):616-628
We consider an electrically conducting fluid confined to a thin rotating spherical shell in which the Elsasser and magnetic Reynolds numbers are assumed to be large while the Rossby number is assumed to vanish in an appropriate limit. This may be taken as a simple model for a possible stable layer at the top of the Earth's outer core. It may also be a model for the thin shells which are thought to be a source of the magnetic fields of some planets such as Mercury or Uranus. Linear hydromagnetic waves are studied using a multiple scale asymptotic scheme in which boundary layers and the associated boundary conditions determine the structure of the waves. These waves are assumed to be of the form of an asymptotic series expanded about an ambient magnetic field which vanishes on the equatorial plane and velocity and pressure fields which do not. They take the form of short wave, slowly varying wave trains. The results are compared to the author's previous work on such waves in cylindrical geometry in which the boundary conditions play no role. The approximation obtained is significantly different from that obtained in the previous work in that an essential singularity appears at the equator and nonequatorial wave regions appear. 相似文献
15.
考虑各向异性的层土-盾构隧道地震反应数值模拟 总被引：3，自引：0，他引：3
在层状各向异性土体-盾构隧道地震反应分析中,引入了横观各向同性弹塑性模型理论,建立了适合于横观各向同性介质的双渐近-多向透射边界条件。针对地铁区间盾构隧道抗震设计的特点,基于横观各向同性弹塑性模型,研制了考虑层状土体各向异性和施工开挖效应,适合于盾构隧道动力计算的各向异性弹塑性动力有限元程序。在程序中对于不同的材料采用了不同的本构关系和不同的单元形式,并采用了关联流动法则和多种屈服准则,可同时进行各向异性土体与地下结构的二维平面应力、平面应变和轴对称问题的静力和动力数值分析。最后利用所研制的程序进行了上海地铁二号线石门一路站附近区间隧道在不同超越概率地震动输入下的隧道反应计算。结果表明,在层状土体-地铁区间隧道的抗震分析中考虑土体各向异性的影响是必要的,所提出的计算模型是可行的。 相似文献
16.
17.
A general analytical solution for flow to a single horizontal well by Fourier and Laplace transforms
The objective of this paper is to present an analytical solution for describing the head distribution in an unconfined aquifer with a single pumping horizontal well parallel to a fully penetrating stream. The Laplace-domain solution is developed by applying Fourier sine, Fourier and Laplace transforms to the governing equation as well as the associated initial and boundary conditions. The time-domain solution is obtained after taking the inverse Laplace transform along with the Bromwich integral method and inverse Fourier and Fourier sine transforms. The upper boundary condition of the aquifer is represented by the free surface equation in which the second-order slope terms are neglected. Based on the solution and Darcy’s law, the equation representing the stream depletion rate is then derived. The solution can simulate head distributions in an aquifer infinitely extending in horizontal direction if the well is located far away from the stream. In addition, the solution can also simulate head distributions in confined aquifers if specific yield is set zero. It is shown that the solution can be applied practically to evaluate flow to a horizontal well. 相似文献
18.
Two-layer equatorial primitive equations for the free troposphere in the presence of a thin atmospheric boundary layer and thermal dissipation are developed here. An asymptotic theory for the resonant nonlinear interaction of long equatorial baroclinic and barotropic Rossby waves is derived in the presence of such dissipation. In this model, a self-consistent asymptotic derivation establishes that boundary layer flows are generated by meridional pressure gradients in the lower troposphere and give rise to degenerate equatorial Ekman friction. That is to say, the asymptotic model has the property that the dissipation matrix has one eigenvalue which is nearly zero: therefore the dynamics rapidly dissipates flows with pressure at the base of the troposphere and creates barotropic/baroclinic spin up/spin down. The simplified asymptotic equations for the amplitudes of the dissipative equatorial barotropic and baroclinic waves are studied by linear theory and integrated numerically. The results indicate that although the dissipation slightly weakens the tropics to midlatitude connection, strong localized wave packets are nonetheless able to exchange energy between barotropic and baroclinic waves on intraseasonal timescales in the presence of baroclinic mean shear. Interesting dissipation balanced wave-mean flow states are discovered through numerical simulations. In general, the boundary layer dissipation is very efficient for flows in which the barotropic and baroclinic components are of the same sign at the base of the free troposphere whereas the boundary layer dissipation is less efficient for flows whose barotropic and baroclinic components are of opposite sign at the base of the free troposphere. 相似文献