Properties of High-Order Finite Difference Schemes and Idealized Numerical Testing     

Daosheng XU  Dehui CHEN  Kaixin WU 《大气科学进展》2021,38(4):615-626

Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, the properties of high-order finite difference schemes are studied based on idealized numerical testing, for the purpose of their application in the Global/Regional Assimilation and Prediction System(GRAPES) model. It is found that the pros and cons due to grid staggering choices diminish with higher-order schemes based on linearized analysis of the one-dimensional gravity wave equation. The improvement of higher-order difference schemes is still obvious for the mesh with smooth varied grid distance. The results of discontinuous square wave testing also exhibits the superiority of high-order schemes. For a model grid with severe non-uniformity and non-orthogonality, the advantage of high-order difference schemes is inapparent, as shown by the results of two-dimensional idealized advection tests under a terrain-following coordinate. In addition, the increase in computational expense caused by high-order schemes can be avoided by the precondition technique used in the GRAPES model. In general, a high-order finite difference scheme is a preferable choice for the tropical regional GRAPES model with a quasi-uniform and quasi-orthogonal grid mesh.  相似文献   

基于人工地震波作用下大跨度悬索桥多点激励地震响应分析          下载免费PDF全文

易富  杜常博  王伟  孙悦  张智 《地震工程学报》2019,41(2):278-285

为研究大跨度悬索桥在多点激励作用下的地震响应规律,在已有的功率谱模型基础上提出了改进的功率谱模型,并将其应用到人工地震波合成过程中。以某水库上主跨720 m双塔单跨悬索桥为研究对象,用Midas civil建立全桥有限元模型,采用大质量法进行不同波速多点激励地震响应分析。结果表明:大跨度悬索桥在多点激励作用下的主塔内力响应、主塔位移响应及主梁位移响应均受行波效应、衰减效应、不相干效应及衰减后地震波叠加效应影响,且最终的响应值由衰减后地震波的叠加效应与衰减效应的共同作用决定;500 m/s多点激励下衰减后的地震波叠加加强效应与衰减效应共同作用后,对主塔轴力响应、主塔塔底顺桥向剪力响应、左塔上横梁处剪力和弯矩响应的加强效果最大;各波速多点激励下,主塔顶顺桥向位移响应相对变化率均大于零且几乎保持不变,主梁两端顺桥向位移响应相对变化率在1 000 m/s波速取得最大值,主梁竖向位移响应在500 m/s取得最大值并随波速的增加逐渐接近一致激励情况。  相似文献   

基于多项式拟合均匀荷载面曲率曲线的结构损伤识别方法          下载免费PDF全文

唐盛华  苏彬建  方志  张学兵 《世界地震工程》2018,34(1):121-130

为仅利用结构损伤状态的柔度矩阵对结构进行损伤程度识别，先对损伤状态的均匀荷载面曲率曲线进行最小二乘法拟合。根据曲率曲线差判断结构的损伤位置，对损伤位置的点进行剔除后，再利用未损伤位置上的点进行局部最小二乘法拟合，代替损伤前的均匀荷载面曲率曲线，用于结构的损伤定位与定量。通过一简支梁数值算例，先以理论的二次多项式进行拟合，考虑单损伤和多损伤的情况，进行损伤识别分析，再分析多项式次数、测点数目以及不同噪声水平对损伤定量精度的影响。结果表明:在一定范围内，次数越高拟合误差越小，但差别不明显，采用理论的二次多项式拟合即可满足结构损伤识别要求，无噪声的情况下，测点数目减半不影响损伤识别的精度，该方法具有一定的抗噪声能力。  相似文献   

场地地震反应分析中侧向人工边界距离选取及影响因素分析          下载免费PDF全文

白建方  董士欣 《世界地震工程》2018,34(3):179-187

结合均匀场地模型和分层均匀场地模型，采用数值试验的方法研究了场地地震反应分析中侧向人工边界距离的选取问题。分析了各主要因素对人工边界距离选取的影响，得到了如下结论:单侧人工边界距离与土层深度之比L/h是影响人工边界距离选取的主要参数；随着侧向人工边界距离的增加，有限元模型的计算误差并不一直减小，在局部范围内，会出现计算误差随着边界距离的增加而增大的现象；场地硬度、阻尼比和激振频率等因素不会影响计算误差随边界距离的变化规律，但采用自由场边界时受输入地震动频谱特性的影响较为明显；对水平激振下的一维场地模型来说，采用结合自由度边界时所需边界距离最短，自由场边界次之，采用粘性边界达到收敛时所需边界距离最长；无论哪种边界条件，当边界距离大于某一范围后，都能得到稳定的计算结果；相对于均质场地模型来说，实际的分层均匀场地建模时可适当减小侧向人工边界距离。  相似文献   

卫星遥感中风场对降雨飘移距离的影响及误差校正     

钱程程  陈戈 《遥感学报》2011,15(5):1079-1091

首先利用海表面风速与高度的对数关系、风向与高度的关系,推算出高空的风速、风向及相应风力的大小。在此基础上根据降雨的形成高度,假定在海表面10 m处风速为10 m/s的情况下,得到单位质量为1的雨水在风场影响下的飘移距离在3—100 km之间,可以判定这个偏移距离即为在风场影响下卫星观测降雨的误差。利用上述方法对10组TRMM卫星降雨产品数据进行纠正,将纠正前后的TRMM数据与GPCP数据进行对比。结果显示,经过纠正的TRMM数据与GPCP数据更为一致,整体精度提高了12%左右,验证了方法的正确性。  相似文献   

利用改进遗传算法求解整周模糊度   总被引：2，自引：0，他引：2  

邢喆  樊妙 《测绘科学》2011,36(3):110-113

本文首先通过整数高斯变换,降低了各整周模糊度分量之间的相关性,并在此基础上提出了均匀设计与遗传算法相结合的方法进行整周模糊度搜索,最后进行了基线解算实验。解算结果表明,该算法用于模糊度求解极大提高了搜索效率。  相似文献   

Earthquake source parameters of the 2009 M_W7.8 Fiordland （New Zealand） earthquake from L-band InSAR observations     

Kateline Young 《地震学报(英文版)》2011,(2):199-206

The 2009 M W 7.8 Fiordland (New Zealand) earthquake is the largest to have occurred in New Zealand since the 1931 M W 7.8 Hawke’s Bay earthquake, 1 000 km to the northwest. In this paper two tracks of ALOS PALSAR interferograms (one ascending and one descending) are used to determine fault geometry and slip distribution of this large earthquake. Modeling the event as dislocation in an elastic half-space suggests that the earthquake resulted from slip on a SSW-NNE orientated thrust fault that is associated w...  相似文献   

A universal algorithm to generate pseudo-random numbers based on uniform mapping as homeomorphism     

王福来 《海洋学报(英文版)》2010,29(9)

A specific uniform map is constructed as a homeomorphism mapping chaotic time series into [0,1] to obtain sequences of standard uniform distribution. With the uniform map, a chaotic orbit and a sequence orbit obtained are topologically equivalent to each other so the map can preserve the most dynamic properties of chaotic systems such as permutation entropy. Based on the uniform map, a universal algorithm to generate pseudo random numbers is proposed and the pseudo random series is tested to follow the standard 0-1 random distribution both theoretically and experimentally. The algorithm is not complex, which does not impose high requirement on computer hard ware and thus computation speed is fast. The method not only extends the parameter spaces but also avoids the drawback of small function space caused by constraints on chaotic maps used to generate pseudo random numbers. The algorithm can be applied to any chaotic system and can produce pseudo random sequence of high quality, thus can be a good universal pseudo random number generator.  相似文献   

潮汐河口水库中水力充填堤坝结构的渗控关键技术研究     

王生海  邓鹏 《中国勘察设计》2014,(5):66-69

随着近些年水力充填堤坝结构在潮汐河口水库中的应用和发展,渗流控制技术也随之不断完善。本文以上海青草沙水库工程中的北堤工程为例,论述了水力充填堤坝应用在潮汐河口水库中,其渗流控制技术研究分析过程,为类似工程的设计和研究提供了经验和参考。  相似文献   

Elastic solutions for stresses in a transversely isotropic half‐space subjected to three‐dimensional buried parabolic rectangular loads     

Cheng‐Der Wang  Jyh‐Jong Liao 《国际地质力学数值与分析法杂志》2002,26(14):1449-1476

In many areas of engineering practice, applied loads are not uniformly distributed but often concentrated towards the centre of a foundation. Thus, loads are more realistically depicted as distributed as linearly varying or as parabola of revolution. Solutions for stresses in a transversely isotropic half‐space caused by concave and convex parabolic loads that act on a rectangle have not been derived. This work proposes analytical solutions for stresses in a transversely isotropic half‐space, induced by three‐dimensional, buried, linearly varying/uniform/parabolic rectangular loads. Load types include an upwardly and a downwardly linearly varying load, a uniform load, a concave and a convex parabolic load, all distributed over a rectangular area. These solutions are obtained by integrating the point load solutions in a Cartesian co‐ordinate system for a transversely isotropic half‐space. The buried depth, the dimensions of the loaded area, the type and degree of material anisotropy and the loading type for transversely isotropic half‐spaces influence the proposed solutions. An illustrative example is presented to elucidate the effect of the dimensions of the loaded area, the type and degree of rock anisotropy, and the type of loading on the vertical stress in the isotropic/transversely isotropic rocks subjected to a linearly varying/uniform/parabolic rectangular load. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献