排序方式： 共有29条查询结果，搜索用时 15 毫秒
21.
利用高阶Li空间微分方案(Li, 2005),实现了时间积分为3～6阶Runge-Kutta-Li(RKL)格式的求解算法。二维线性平流方程的试验结果表明:在计算稳定的条件下,各阶算法的计算误差随时间的推移基本上是线性增加的。非转动背景场的平流算例中（高斯型的初值）,高阶RKL算法可以取得较好的计算效果。与3、4、5、6阶RK算法配合的Li空间差分方案有效阶数可以达到5、7、9、10阶。RK 算法的阶数为5（6）阶时,总误差控制在10-7（10-8）以内。随RK阶数增加Li微分的有效阶数有增加趋势,且总误差逐渐减小。定常转速的背景场算例中（偏心的高斯型初值）,当RK阶数为3时,最优空间差分阶数为10;相应的阶数为4、5、6时对应的空间最优阶为16,22,22,总计算误差可以控制在10-15～10-16。随着精度的提高,误差的绝对值减小很迅速,说明算法是非常有效的。对于圆锥型初值（定常转速的背景场）,4、5、6阶RK算法和3阶算法的效果差不多。高阶算法对此类具有导数不连续点的算例,效果不如高斯初始场好,结果不能保持正定,有些地方误差出现下冲和上翘。随着空间差分精度的提高,非正定的解数量和数值减小,误差的绝对值减小,说明了算法在一定程度上是有效的,但并不适合追求极高的算法阶数。这与谱方法中的导数不连续问题有些相似,误差的产生主要源于导数的不连续性,差分类方法仅能获得与导数连续性阶数相当的算法精度。各种算例中,采用恰当的边界条件是必要的,例如旋转背景场算例,比较适合使用无穷远边界条件,否则会出现计算不稳定或无法将计算误差控制到较小的范围内。 相似文献
22.
A high-order Rankine panel method based on Non-Uniform Rational B-Spline (NURBS) is developed for solving the three-dimensional radiation and diffraction problems with forward speed. A NURBS surface is used to precisely represent the body geometry. Velocity potential on the body surface is described by B-spline after the source density distribution on the boundary surface is determined. A collocation approach is applied to numerical computation and the integral equations are evaluated by applying Gauss–Legendre quadrature. The mj-terms are evaluated by a desingularized method which utilizes NURBS technique. In order to verify the method proposed, it is firstly applied to the unbounded flow problem of a sphere and spheroids. The numerical results are found to be in good agreement with analytical solutions. Then the method is used to solve the radiation and diffraction problems of a sphere and the diffraction problem of a spheroid moving with a forward speed beneath the free surface in frequency domain. The numerical results are satisfactory in comparison with the published analytical results and experimental results. 相似文献
23.
该文为高阶矩湍流方程组的二阶矩闭合方案提供了参数组选择。此参数组是用近地层实测资料结合湍流方程二阶矩闭合方案推算出的。然后，将其与Mellor等人在实验室所得的参数组分别代入湍流二阶矩闭合方程组拟合，得动量(ΦM(ξ))、热量(ΦH(ξ))与稳定度(ξ)的二组关系曲线，可见此二组曲线几乎重合，且与Businger野外近地层实测结果符合较好。拟合的垂直速度方差((W*2)1/2(ξ))、温度方差((θ*2)1/2(ξ))与ξ的二组关系曲线相对于UTAH资料，所得曲线偏差比Mellor曲线偏差小得多，而相对于KANSAS资料，所得曲线偏差稍大于Mellor曲线。 相似文献
24.
The aim of this experimental study is to investigate the interaction between turbulent flow and a gravelbed that mimics the actual roughness structures of a natural bed and its implication on sediment transport.In particular,the response of the Reynolds stresses and the role of intermittency to the bed roughness is the primary focus of the current study.To this end,the flow field,measured with an Acoustice Doppler Velocimeter (ADV),is thoroughly examined,considering the conditional Reynolds shea... 相似文献
25.
Kasemsan Manomaiphiboon Sun-Kyoung Park Armistead G. Russell 《Stochastic Environmental Research and Risk Assessment (SERRA)》2008,22(2):159-168
Probabilistic characterization of environmental variables or data typically involves distributional fitting. Correlations,
when present in variables or data, can considerably complicate the fitting process. In this work, effects of high-order correlations
on distributional fitting were examined, and how they are technically accounted for was described using two multi-dimensional
formulation methods: maximum entropy (ME) and Koehler–Symanowski (KS). The ME method formulates a least-biased distribution
by maximizing its entropy, and the KS method uses a formulation that conserves specified marginal distributions. Two bivariate
environmental data sets, ambient particulate matter and water quality, were chosen for illustration and discussion. Three
metrics (log-likelihood function, root-mean-square error, and bivariate Kolmogorov–Smirnov statistic) were used to evaluate
distributional fit. Bootstrap confidence intervals were also employed to help inspect the degree of agreement between distributional
and sample moments. It is shown that both methods are capable of fitting the data well and have the potential for practical
use. The KS distributions were found to be of good quality, and using the maximum likelihood method for the parameter estimation
of a KS distribution is computationally efficient. 相似文献
26.
数字岩石物理利用三维成像技术和数学方法,建立数字岩心,开展多场物理响应模拟,计算岩石的等效弹性参数,为岩石物理学研究开拓了新的领域.本文发展了一种基于高阶有限差分的岩石模量数值计算的新方法,该方法便于理解,容易实现,占用内存相对较少,计算效率高,结果合理,弥补了常规岩石物理实验周期长,成本高,误差大等不足.该方法将三维数字岩心样本嵌入一个具有与岩心样本骨架颗粒相同弹性性质的区域中扩充成一个新模型,测量由数字岩石样本非均匀结构带来的纵波或横波峰值振幅的时间差,利用该时间差(相对参考模型)估算纵横波的等效速度,进而求取等效弹性模量.理论模型和实际岩心的计算结果表明,数值模拟结果与实验结果有较高的吻合度,验证了该方法的合理性. 相似文献
27.
以先进的波动理论为基础的波动方程保幅地震偏移成像是在给出正确位置的同时也给出真实振幅的一种特殊完善.作者从保幅单程波动方程的非稳态相移公式出发,基于反问题求解中常用的摄动理论,利用单平方根算子的渐进展开,从而推导出保幅叠前深度偏移方程的高阶广义屏形式;针对散射波场计算项对于横向变速介质的不稳定性,通过数学近似提出一个有效提高稳定性的策略,应用到波场递归外推过程中,从而得到一种稳定的保幅高阶广义屏叠前深度偏移算子.理论模型试算和实际资料处理表明,该方法不但可以更精确地使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,使AVO响应更加清晰,提高了AVO资料的分析精度. 相似文献
28.
Alberto Canestrelli Michael Dumbser Annunziato Siviglia Eleuterio F. Toro 《Advances in water resources》2010
In this paper, we study the numerical approximation of the two-dimensional morphodynamic model governed by the shallow water equations and bed-load transport following a coupled solution strategy. The resulting system of governing equations contains non-conservative products and it is solved simultaneously within each time step. The numerical solution is obtained using a new high-order accurate centered scheme of the finite volume type on unstructured meshes, which is an extension of the one-dimensional PRICE-C scheme recently proposed in Canestrelli et al. (2009) [5]. The resulting first-order accurate centered method is then extended to high order of accuracy in space via a high order WENO reconstruction technique and in time via a local continuous space–time Galerkin predictor method. The scheme is applied to the shallow water equations and the well-balanced properties of the method are investigated. Finally, we apply the new scheme to different test cases with both fixed and movable bed. An attractive future of the proposed method is that it is particularly suitable for engineering applications since it allows practitioners to adopt the most suitable sediment transport formula which better fits the field data. 相似文献
29.
为解决3D AVO地震数据快速保幅重建问题,在传统3D抛物Radon变换的基础上提出一种3D快速高阶抛物Radon变换方法.该方法将传统抛物Radon变换与正交多项式相结合,通过正交多项式系数描述地震数据AVO信息,确保重建后的地震数据具有良好保幅效果.同时,该方法引入新变量λ_x=q_xf和λ_y=q_yf,通过对q_xf和q_yf的整体采样,消除了3D高阶抛物Radon变换算子对频率的依赖,使变换算子的求逆过程仅需计算一次,大大节省计算时间.理论模型和实际地震资料的处理结果表明,该方法重建效率高,保幅效果良好. 相似文献