收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

海洋声学中缓变分层波导的声波传播计算

李鹏  
【摘要】: 在声学、电磁学、地震学和一些其他的应用中,有许多在具有分层缓变波导中大尺度的波传播的问题。如果使用直接方法,譬如有限元法、有限差分法等,往往会导致求解大规模的线性方程组,这使得使得算法在内存和计算量的花费非常大,计算效率不高。在本文中,我们针对具有弯曲界面或者弯曲边界的波导,提出了适合于步进计算一般的坐标变换变换的方法,在坐标变换后的坐标系中,使用基于DtN(Dirichlet to Neumann)映射的正向基本解算子步进方法(FFOMM)和逆向基本解算子步进方法(IFOMM)分别求解Helmholtz方程正问题和反问题。 在海洋声学中,海洋环境可以看作是一种分层缓变的具有多个弯曲内部界面和多层介质的声波导。在这样的波导中如果要使用步进方法大步长的求解Helmholtz方程,首先需要用适合于步进计算坐标变换将波导中弯曲的界面或者弯曲边界拉直。本文中给出两种坐标变换的方法:第一种是局部的解析正交坐标变换法,局部解析坐标变换用一个解析的积分方程联系新旧坐标,使用牛顿法求解这一积分方程实现新旧坐标之间的转换。坐标变换必然相应的引起方程的变换,我们使用的方程变换将Helmholtz方程转化成为一个可以用步进方法求解的形式。对于具有一个平的顶部、平的底部、n个弯曲的内部界面的二维波导,我们给出了在变换方程系数在各个介质层一般形式的解析表达式,从而完成方程变换。然而,局部解析坐标变换方法需要满足在两个相邻界面之间存在一条水平直线的介质可分性条件。针对不满足可分性条件的情况,我们又提出了一种使用经典求解常微分方程的Runge-Kutta方法计算坐标变换的数值坐标变换方法,这一方法不仅具有更广的适用范围,而且不必引入辅助界面,从而更好的逼近原始问题。 对于变换之后的计算坐标系中的正问题,我们使用正向基本解步进方法来求解方程。FFOMM基于DtN算子的单向重建,将边界值问题转化为初始值问题,对于DtN算子和基本解算子满足的方程使用大步长方法离散在传播方向上的变量,用截断特征值展开方法逼近这些算子。FFOMM对于求解大尺度缓变波导中波的传播问题非常有效。 在具有分层缓变特性的波导中的波的大尺度传播问题中,反问题有非常重要的实际意义。一般说来,海洋声学中实际的反问题都很复杂。求解海洋中反问题一般使用两类方法:一种是基于正向方法(如FFOMM)的迭代法,另一种为反向方法,反向方法的运算量和计算速度都相比正向方法有很大的优势。基于FFOMM,我们提出了逆向基本解算子步进方法。在IFOMM中使用波导中存在的传播模的个数作为正则化参数,用截断的奇异值分解方法求解步进过程中出现的病态线性方程。数值计算证明,这个方法用于求解Helmholtz方程的反边界值有效,精确并且对于初始波源传播部分稳定,对于海洋声波中实际的反问题研究有着一定的意义。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 刘克英;李鹏;;基于Tikhonov正则化的逆向基本解算子步进方法[J];华北水利水电学院学报;2011年02期
2 李鹏;蔡洪涛;;波在多层介质缓变波导中的传播[J];华北水利水电学院学报;2010年05期
3 彭桓武;包含界面条件在内的变分近似法[J];科学通报;1975年09期
4 蔡煜东,姚林声;人工神经网络在白云岩类型识别中的应用[J];岩石矿物学杂志;1994年01期
5 杨光;一维线性抛物型方程边界值反演的TBG方法[J];石油化工高等学校学报;1998年03期
6 孙芳;曾周末;王晓媛;靳世久;詹湘琳;;界面条件下线型超声相控阵声场特性研究[J];物理学报;2011年09期
7 谢鸿喜;普救寺的建筑与《西厢记》传播模式[J];三晋测绘;1994年01期
8 孙雁;钟万勰;;电磁波导的通过谱计算[J];计算力学学报;2006年06期
9 R. T. Coates;断层与裂隙的有限差分模拟[J];勘探地球物理进展;1996年05期
10 岑人经;冠脉入口区域血液流动分析[J];华南理工大学学报(自然科学版);1982年01期
11 张业荣,张天仁,黄天霖,张惕远,漆兰芬;室内环境的视距传播模型[J];电波科学学报;1998年03期
12 杨永;张巧;;圆形区域内熔体凝固的交界面稳定性分析[J];科技信息;2009年16期
13 苏铭德;有间断界面的二维可压缩理想流体流场的变分原理[J];空气动力学学报;1985年03期
14 陆章基;郑广平;;牛顿-微极分层流体的狭窄流动[J];汕头大学学报(自然科学版);1990年01期
15 ;Geophysics,Vol.56,No.2,1991论文文摘[J];勘探地球物理进展;1991年05期
16 李德筠,王日新,沈国光,程虹;连续分层重力内波的数学模型及其传播模式[J];天津大学学报;1996年06期
17 黄永年,王志丹,王永利;高能太阳宇宙线行星际传播的Green函数及传播模型[J];中国科学A辑;1998年06期
18 朱少英,徐渝;基于组织学习的知识动态传播模型[J];科研管理;2003年01期
19 曾金平;叶玉其;;基于Robin界面条件的非重叠区域分解法[J];湖南大学学报(自然科学版);2007年06期
20 张玉君;环状单电极测井理论[J];物探与化探;1979年02期
中国重要会议论文全文数据库 前10条
1 唐少强;;反射系数法不能正确判别界面条件稳定性[A];中国力学学会学术大会'2009论文摘要集[C];2009年
2 唐少强;;多尺度计算数值界面条件的设计与分析[A];中国力学学会学术大会'2009论文摘要集[C];2009年
3 唐少强;;多尺度计算数值界面条件的设计与分析[A];现代数学和力学(MMM-XI):第十一届全国现代数学和力学学术会议论文集[C];2009年
4 杜舒舒;龚建华;程承祺;;Multi-Agent技术在SARS传播模拟中的应用研究[A];中国地理信息系统协会第三次代表大会暨第七届年会论文集[C];2003年
5 林伟军;王秀明;张海澜;;弹性波在非均匀各向异性介质中的传播模拟及其并行算法[A];2004年全国物理声学会议论文集[C];2004年
6 陈永强;徐鉴君;;雾化液粒凝固过程的数学模型及其渐近解[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
7 方明;唐少强;李志辉;;面心立方晶格原子模拟的人工界面条件及其应用[A];北京力学会第18届学术年会论文集[C];2012年
8 倪尔瑚;;介质谐振器材料的毫米波测量[A];中国电子学会第十三届电子元件学术年会论文集[C];2004年
9 沈星星;张吉;谭建龙;郭丽;;基于数据流管理平台的未知蠕虫发现[A];全国网络与信息安全技术研讨会’2004论文集[C];2004年
10 徐烈;徐瑞萍;;基于Cantor集分形的固体接触界面热阻研究[A];上海市制冷学会2007年学术年会论文集[C];2007年
中国博士学位论文全文数据库 前9条
1 李鹏;海洋声学中缓变分层波导的声波传播计算[D];浙江大学;2008年
2 唐敏;输运方程与波方程的渐近保持格式[D];清华大学;2008年
3 陈高洁;椭圆型变分问题的区域分解法[D];湖南大学;2006年
4 张梅;颗粒复合材料等效热导率的加权残值法和边界元法研究[D];武汉理工大学;2008年
5 刘建光;粘性介质压力成形有限元分析与实验研究[D];哈尔滨工业大学;2006年
6 石零;低温固—固接触界面热传输研究[D];华中科技大学;2006年
7 曹伟;注塑成型中粘性不可压缩非等温熔体的三维流动模拟[D];郑州大学;2004年
8 张红旗;沿涂敷介质层的导电基底传播的表面波与侧面波[D];中国科学院陕西天文台;2001年
9 董必昌;岩土工程仿真中接触单元和相关参数研究[D];华中科技大学;2005年
中国硕士学位论文全文数据库 前10条
1 解中力;带自由界面光波导中的泄漏模—渐近解分析[D];浙江大学;2007年
2 邓小毛;弯曲界面声波导中共轭算子的构造和分析[D];浙江大学;2008年
3 周巧仙;Helmholtz方程在三层平板光波导中的特征计算[D];浙江大学;2005年
4 左冬梅;本体与界面条件下有机硅氧烷水解—缩合过程簇集行为的比较研究[D];陕西师范大学;2008年
5 梁冠军;交通荷载作用下城市沥青道路路基的动力响应和塑性变形分析[D];同济大学;2007年
6 胡卫国;沥青路面基面层间界面抗剪研究[D];长沙理工大学;2007年
7 黄文星;特征模展开方法应用于声波导计算的有效性研究[D];浙江大学;2010年
8 钱微微;一类带PML声波导中的共轭特征函数构造及其应用[D];浙江大学;2006年
9 高振严;几种有效的数值方法在光波导中的应用[D];浙江大学;2007年
10 俞坚钢;Ni/Ti中子多层膜周期结构的研究[D];首都师范大学;2009年
中国重要报纸全文数据库 前6条
1 张中桥;二十一世纪医学展望[N];西安日报;2000年
2 张树;解剖蠕虫[N];中国计算机报;2004年
3 金晶 张兆杰;互联网掀起的音乐革命[N];江苏经济报;2004年
4 ;从网络层防治蠕虫[N];中国计算机报;2005年
5 王飞燕;中国科协与联合国儿童基金会项目官员在我区考察[N];内蒙古日报(汉);2000年
6 信息产业部电信研究院规划所 贺丰;3G规划软件市场大有可为[N];通信产业报;2005年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978