A nearly analytic symplectic partitioned Runge–Kutta method based on a locally one-dimensional technique for solving two-dimensional acoustic wave equations     

Chol Sim  Chunyou Sun  Nam Yun 《Geophysical Prospecting》2020,68(4):1253-1269

In this paper, we develop a new nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique for numerically solving two-dimensional acoustic wave equations. We first split two-dimensional acoustic wave equation into the local one-dimensional equations and transform each of the split equations into a Hamiltonian system. Then, we use both a nearly analytic discrete operator and a central difference operator to approximate the high-order spatial differential operators, which implies the symmetry of the discretized spatial differential operators, and we employ the partitioned second-order symplectic Runge–Kutta method to numerically solve the resulted semi-discrete Hamiltonian ordinary differential equations, which results in fully discretized scheme is symplectic unlike conventional nearly analytic symplectic partitioned Runge–Kutta methods. Theoretical analyses show that the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique exhibits great higher stability limits and less numerical dispersion than the nearly analytic symplectic partitioned Runge–Kutta method. Numerical experiments are conducted to verify advantages of the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique, such as their computational efficiency, stability, numerical dispersion and long-term calculation capability.  相似文献   

不同非感应起电及感应起电参数化方案对青海东部一次雷暴云电荷结构影响的数值模拟研究          下载免费PDF全文

李江林  余晔  李万莉  李亚珺 《地球物理学报》2019,62(7):2366-2381

本研究利用加入起电、放电参数化方案的数值模式（Weather Research and Forecasting Model（Version 3.7.1），WRF3.7.1_ELEC），通过设计五组不同非感应起电及感应起电参数化方案敏感性试验，对发生在青藏高原东北部青海大通地区的一次雷暴过程进行模拟研究，对比分析了不同非感应起电机制及感应起电机制对雷暴云电荷结构的影响.结果表明：在雷暴云发展旺盛阶段，Saunders（S91）、Riming Rate（RR）、和Saunders和Peck（SP98）三种非感应起电方案模拟的雷暴云最低层均为负电荷区，而混合方案（Brooks and SP98，BSP）模拟的雷暴云最低层为正电荷区，主电荷区自下而上为"+-+-"排列的四层电荷结构.与甚高频辐射源定位法推算的结果对比，BSP方案模拟的本次高原雷暴云电荷结构更接近实际情况；几种不同非感应起电方案模拟的主电荷区外围与主电荷区电荷结构不同，说明在雷暴发展的不同阶段雷暴云的电荷结构是不同的；几种非感应起电方案模拟的电荷结构不尽相同，主要是由于霰、冰和雪粒子在不同高度所带电荷的极性及电量的大小不同，霰粒子的电荷密度对低层的影响较大，冰粒子和雪粒子的电荷密度对中上层的影响较大；加入感应起电机制后，雷暴云电荷结构分布几乎没有变化，但能使雷暴云发展旺盛阶段低层和中层的正负电荷区电荷密度有所加强.  相似文献   

Uncertainty in marine weather routing     

《Applied Ocean Research》2019

Weather routing methods are essential for planning routes for commercial shipping and recreational craft. This paper provides a methodology for quantifying the significance of numerical error and performance model uncertainty on the predictions returned from a weather routing algorithm. The numerical error of the routing algorithm is estimated by solving the optimum path over different discretizations of the environment. The uncertainty associated with the performance model is linearly varied in order to quantify its significance. The methodology is applied to a sailing craft routing problem: the prediction of the voyaging time for an ethnographic voyaging canoe across long distance voyages in Polynesia. We find that the average numerical error is an order of magnitude smaller than the performance model uncertainty. These results illustrate the significance of considering the influence of numerical error and performance uncertainty when performing a weather routing study.  相似文献   

Climatological Characteristics of the East Asian Winter Monsoon Simulated by CWRF Regional Climate Model     

Bingdi Wang  Qingquan Li  Xinyong Shen  Lili Dong  Fang Wang  Tao Wang  Xinzhong Liang 《地球科学进展》2020,35(3):319-330

Based on the ERA-Interim atmospheric reanalysis data from the European Medium-Term Weather Forecast Center from 1979 to 2016 and the ERSSTv4 sea surface temperature data from the US National Oceanic and Atmospheric Administration, the regional climate model CWRF was used to simulate the climate characteristics in East Asia. The results show that the CWRF model can well reproduce the average characteristics of the East Asian winter monsoon circulation, including the location and intensity of the low-level continental cold high pressure and variation characteristics of wind field in high and low levels. The occurrence area and frequency of the north wind in the simulation and the reanalysis data were further calculated and compared. It is shown that they are basically consistent. The distribution of air temperature and precipitation over China are well represented by the model. The water vapor transport is also in good agreement with the reanalysis data. The water vapor from the Bay of Bengal plays a vital role in the precipitation over South China. The simulation results of apparent heat source and apparent moisture sink show that the model can well simulate the thermal difference between the East Asian continent and the adjacent sea area. The analysis results indicate that CWRF model has the ability to simulate the main characteristics of the East Asian winter monsoon.  相似文献   

建筑结构地震易损性分析研究综述   总被引：4，自引：4，他引：0  

贾晗曦  林均岐  刘金龙 《震灾防御技术》2019,14(1):42-51

地震易损性是指在不同强度水平的地震作用下建筑结构发生各种破坏的条件概率。本文对地震易损性分析的发展历程进行回顾和总结，介绍了地震易损性的相关概念、研究目的和意义以及5个表达工具。随后，针对国内外学者对建筑结构的易损性分析方法开展了系统的梳理和分析，特别是从基于震害调查的经验分析法和基于数值模拟的理论计算法方面展开了重点综述，并简单介绍了一些不常用的方法。最后指出了国内外对建筑结构易损性分析研究的不足和未来的研究方向。  相似文献   

Numerical assessment of tsunami attack on a rubble mound breakwater using OpenFOAM®     

《Applied Ocean Research》2018

A numerical assessment study of tsunami attack on the rubble mound breakwater of Haydarpasa Port, located at the southern entrance of the Istanbul Bosphorus Strait in the Sea of Marmara, Turkey, is carried out in this study using a Volume-Averaged Reynolds-Averaged Navier-Stokes solver, IHFOAM, developed in OpenFOAM^® environment. The numerical model is calibrated with and validated against the data from solitary wave and tsunami overflow experiments representing tsunami attack. Furthermore, attack of a potential tsunami near Haydarpasa Port is simulated to investigate effects of a more realistic tsunami that cannot be generated in a wave flume with the present state of the art technology. Discussions on practical engineering applications of this type of numerical modeling studies are given focusing on pressure distributions around the crown-wall of the rubble mound breakwater, and the forces acting on the single stone located behind the crown-wall at the rear side of the breakwater. Numerical modeling of stability/failure mechanism of the overall cross-section is studied throughout the paper.The present study shows that hydrodynamics along the wave flume and over the breakwater can be simulated properly for both solitary wave and tsunami overflow experiments. Stability of the overall cross-section can only be simulated qualitatively for solitary wave cases; on the other hand, the effect of the time elapsed during tsunami overflow cannot be reflected in the simulations using the present numerical tool. However, the stability of the overall cross-section under tsunami overflow is assessed by evaluating forces acting on the rear side armor unit supporting the crown-wall of the rubble mound breakwater as a practical engineering application in the present paper. Furthermore, two non-dimensional parameters are derived to discuss the stability of this armor unit; and thus, the stability condition of the overall cross-section. Approximate threshold values for these non-dimensional parameters are presented comparing experimental and numerical results as a starting point for engineers in practice. Finally, investigations on the solitary wave and tsunami overflow experiments/simulations are extended to the potential tsunami simulation in the scope of both representation of a realistic tsunami in a wave flume and stability of the rubble mound breakwater.  相似文献   

基于FLAC3D数值模拟的清徐地裂缝与交城断裂成因关系研究     

曾金艳  李自红  陈文 《震灾防御技术》2020,15(1):114-122

基于断裂力学理论，应用FLAC3D数值模拟方法对交城断裂活动时附近土体的变形特征进行模拟计算。计算结果表明:断层活动带动上覆土层差异沉降及应力场变化，当断层错动量达到一定程度，断层附近出现拉应力区，土体在拉张应力和自身重力作用下产生垂直差异形变，从而产生地裂缝。将数值模拟结果与方山探槽揭示的地裂缝剖面特征进行对比，二者结果相互印证，说明清徐地裂缝是交城断裂活动的结果，二者具有明显的对应关系，地裂缝为构造成因裂缝。数值模拟结果为研究断裂和地裂缝成因关系提供了佐证。  相似文献   

Hydrogeophysical modelling of Hisarcik (Kütahya) geothermal field,western Turkey     

Serkan Üner  Gülçin Özürlan Ağaçgözgü  Doğa Düşünür Doğan 《Geophysical Prospecting》2019,67(8):2176-2195

Western Anatolia hosts many low-to-moderate and high-temperature geothermal sources in which active faults play a dominant role to control the recharge and the discharge of geothermal fluid. In this study, we used the two-dimensional geoelectric structure of Kütahya Hisarcık geothermal field, and created a conceptual hydrogeophysical model that includes faults, real topographical variations and geological units. The temperature distribution and fluid flow pattern are also investigated. The depth extension of Hisarcık Fault, electrical basement and low resistivity anomalies related to the presence of geothermal fluid are determined by using resistivity studies in the area. Numerical simulations suggest that Hisarcık fault functioning as a fluid conduit primarily enables hot fluid to be transported from depth to the surface. It is shown that the locations of predicted outflow vents coincide with those of hot springs in the area.  相似文献   

Travel time computations using a compact eikonal equation for vertical transverse isotropic media          下载免费PDF全文

Peyman Moghaddam  Reza Khajavi  Henk Keers 《Geophysical Prospecting》2018,66(8):1475-1484

Eikonal solvers often have stability problems if the velocity model is mildly heterogeneous. We derive a stable and compact form of the eikonal equation for P‐wave propagation in vertical transverse isotropic media. The obtained formulation is more compact than other formulations and therefore computationally attractive. We implemented ray shooting for this new equation through a Hamiltonian formalism. Ray tracing based on this new equation is tested on both simple as well as more realistic mildly heterogeneous velocity models. We show through examples that the new equation gives travel times that coincide with the travel time picks from wave equation modelling for anisotropic wave propagation.  相似文献   

Numerical study for vegetation effects on coastal wave propagation by using nonlinear Boussinesq model     

《Applied Ocean Research》2018

The vegetation has important impacts on coastal wave propagation. In the paper, the sensitivities of coastal wave attenuation due to vegetation to incident wave height, wave period and water depth, as well as vegetation configurations are numerically studied by using the fully nonlinear Boussinesq model. The model is based on the implementation of drag resistances due to vegetation in the fully nonlinear Boussinesq equation where the drag resistance is provided by the Morison’s formulation for rigid structure induced drag stresses. The model is firstly validated by comparing with the experimental results for wave propagation in vegetation zones. Subsequently, the model is used to simulate waves with different height, period propagating on vegetation zones with different water depth and vegetation configurations. The sensitivities of wave attenuation to incident wave height, wave period, water depth, as well as vegetation configurations are investigated based on the numerical results. The numerical results indicate that wave height attenuation due to vegetation is sensitive to incident wave height, wave period, water depth, as well as vegetation configurations, and attenuation ratio of wave height is increased monotonically with increases of incident wave height and decreases of water depth, while it is complex for wave period. Moreover, more vegetation segments can strengthen the interaction of vegetation and wave in a certain range.  相似文献